From 1b3b07a7f62e351da7467a60edfae69dbcd1afb8 Mon Sep 17 00:00:00 2001
From: bisnsh <bishsh2006@outlook.com>
Date: Mon, 17 Nov 2025 11:16:02 +0800
Subject: [PATCH] 20251117

---
 difun.m            |   8 ++++++
 dms2mat.m          |   7 ++++++
 dms2rad.m          |   2 ++
 jianjiepingcha.m   |  23 ++++++++++++++++++
 jjpc.m             |  11 +++++++++
 jjpcs.m            |   3 +++
 lab1.m             |  13 ++++++++++
 lab1/exercise1.m   |  59 +++++++++++++++++++++++++++++++++++++++++++++
 lab2/dtor.m        |   7 ++++++
 lab2/exercise2.m   |  25 +++++++++++++++++++
 lab2/exercise2.mlx | Bin 0 -> 4268 bytes
 lab2/exercise2_5.m |  14 +++++++++++
 lab2/exercise2_6.m |  16 ++++++++++++
 lab2/exercise2_7.m |  16 ++++++++++++
 lab2/exercise2_9.m |   3 +++
 lab2/fn.m          |  20 +++++++++++++++
 lab2/pmzbfs.m      |  10 ++++++++
 lab2/rad2deg.m     |   7 ++++++
 lab3/exercise3_1.m |   8 ++++++
 lab3/exercise3_2.m |   6 +++++
 lab3/exercise3_3.m |   7 ++++++
 lab3/exercise3_4.m |  10 ++++++++
 lab3/lab3.m        |  42 ++++++++++++++++++++++++++++++++
 lab3/nikon.txt     |   7 ++++++
 lab3/zhuanhuan.txt |   9 +++++++
 mat2dms.m          |   5 ++++
 num2p.m            |   8 ++++++
 qwe.mat            | Bin 0 -> 375 bytes
 rad2dms.m          |   2 ++
 tiaojianpingcha.m  |  31 ++++++++++++++++++++++++
 tjpc.m             |  12 +++++++++
 ttry.m             |   7 ++++++
 untitled3.m        |  10 ++++++++
 untitled5.m        |  30 +++++++++++++++++++++++
 veryfun.m          |  16 ++++++++++++
 35 files changed, 454 insertions(+)
 create mode 100644 difun.m
 create mode 100644 dms2mat.m
 create mode 100644 dms2rad.m
 create mode 100644 jianjiepingcha.m
 create mode 100644 jjpc.m
 create mode 100644 jjpcs.m
 create mode 100644 lab1.m
 create mode 100644 lab1/exercise1.m
 create mode 100644 lab2/dtor.m
 create mode 100644 lab2/exercise2.m
 create mode 100644 lab2/exercise2.mlx
 create mode 100644 lab2/exercise2_5.m
 create mode 100644 lab2/exercise2_6.m
 create mode 100644 lab2/exercise2_7.m
 create mode 100644 lab2/exercise2_9.m
 create mode 100644 lab2/fn.m
 create mode 100644 lab2/pmzbfs.m
 create mode 100644 lab2/rad2deg.m
 create mode 100644 lab3/exercise3_1.m
 create mode 100644 lab3/exercise3_2.m
 create mode 100644 lab3/exercise3_3.m
 create mode 100644 lab3/exercise3_4.m
 create mode 100644 lab3/lab3.m
 create mode 100644 lab3/nikon.txt
 create mode 100644 lab3/zhuanhuan.txt
 create mode 100644 mat2dms.m
 create mode 100644 num2p.m
 create mode 100644 qwe.mat
 create mode 100644 rad2dms.m
 create mode 100644 tiaojianpingcha.m
 create mode 100644 tjpc.m
 create mode 100644 ttry.m
 create mode 100644 untitled3.m
 create mode 100644 untitled5.m
 create mode 100644 veryfun.m

diff --git a/difun.m b/difun.m
new file mode 100644
index 0000000..8a1b267
--- /dev/null
+++ b/difun.m
@@ -0,0 +1,8 @@
+function y=difun(x)
+if x<-pi
+    y=-x-pi
+elseif x>-pi && x<pi
+    y=sin(x)
+else
+    y=(x-pi)/2
+end
\ No newline at end of file
diff --git a/dms2mat.m b/dms2mat.m
new file mode 100644
index 0000000..b67e652
--- /dev/null
+++ b/dms2mat.m
@@ -0,0 +1,7 @@
+function dmsv=dms2mat(dv)
+[l,c]=size(dv);
+for i=1:l
+    dmsv(i,1)=fix(dv);
+    dmsv(i,2)=fix((dv-fix(dv))*60);
+    dmsv(i,3)=(dv-dmsv(i,1)-dmsv(i,2)*60);
+end
\ No newline at end of file
diff --git a/dms2rad.m b/dms2rad.m
new file mode 100644
index 0000000..18157ec
--- /dev/null
+++ b/dms2rad.m
@@ -0,0 +1,2 @@
+function radv=dms2rad(dms)
+radv=dms*pi/180;
\ No newline at end of file
diff --git a/jianjiepingcha.m b/jianjiepingcha.m
new file mode 100644
index 0000000..db0e03c
--- /dev/null
+++ b/jianjiepingcha.m
@@ -0,0 +1,23 @@
+disp('水准网间接平差')
+disp('已知高程')
+Ha=5.015
+Hb=6.016
+disp('观测高差（m）')
+L=[1.359;2.009;0.363;1.012;0.657;-0.357]
+disp('系数矩阵')
+B=[1,0;0,1;1,0;0,1;-1,1;-1,0]
+l=[0;0;4;3;7;2]
+disp('C:')
+C=1*ones(1,6);
+S=[1.1,1.7,2.3,2.7,2.4,4.0]
+P=C./S
+P=diag(P)
+disp('参数的解')
+x=inv(B'*P*B)*B'*P*L
+disp('wucha')
+V=B*x-l
+L1=L+V/1000
+disp('jingdu')
+n=6;
+t=2;
+delta=sqrt(V'*P*V/(n-t))
\ No newline at end of file
diff --git a/jjpc.m b/jjpc.m
new file mode 100644
index 0000000..f397d6d
--- /dev/null
+++ b/jjpc.m
@@ -0,0 +1,11 @@
+function [csgz,qxx,dert0,dlj]=jjpc(bb,ll,pp)
+pn=diag(pp);
+nbb=bb'*pn*bb;
+ww=bb'*pn*ll;
+qxx=inv(nbb);
+csgz=qxx*ww;
+vv=bb*csgz-ll;
+vtpv=ll'*pn*ll-ww'*csgz;
+[rr,nn]=size(bb);
+dert0=sqrt(vtpv/rr);
+dlj=der0*dert0*bb*qxx*bb';
\ No newline at end of file
diff --git a/jjpcs.m b/jjpcs.m
new file mode 100644
index 0000000..2dc420c
--- /dev/null
+++ b/jjpcs.m
@@ -0,0 +1,3 @@
+bbb=[]
+lll=[]
+quan=[]
diff --git a/lab1.m b/lab1.m
new file mode 100644
index 0000000..e086bbd
--- /dev/null
+++ b/lab1.m
@@ -0,0 +1,13 @@
+a=[2,3,6,7;4,2,8,6;3,6,9,2;4,7,9,3]
+a1=a([1:3],[2:3])
+a2=a(2,3)
+zeros(3,4)
+ones(5)
+eye(3)
+rand(4)
+randn(3)
+% 删除行列
+a(:,[3])=[]
+% 增加元素和列
+a(5,6)=45
+a(:,4)=[90,89,34,56,98]
\ No newline at end of file
diff --git a/lab1/exercise1.m b/lab1/exercise1.m
new file mode 100644
index 0000000..869105a
--- /dev/null
+++ b/lab1/exercise1.m
@@ -0,0 +1,59 @@
+% 练习1
+x=linspace(-pi,pi,21)
+y=sin(x)
+% 练习2
+clear
+k=1:1000;
+kk=1./k;
+kkk=kk./k;
+res=sum(kkk)-(pi^2)/6
+% 练习3
+sum(kk-log(1000))
+% 练习4
+clear
+a=rand(2,3)
+x=a(1,:)
+y=a(2,:)
+norm(x)
+norm(y)
+acos(dot(x,y)/norm(x)/norm(y))
+% 练习5
+clear
+rand(3)
+det(ans)  % 线性无关
+% 练习6
+clear
+a=rand(3,2)
+x=a(1,:)
+y=a(2,:)
+z=a(3,:)
+alpha=x-z
+beta=y-z
+alpha=[alpha 0]
+beta=[beta 0]
+cross(alpha,beta)
+% 练习7
+clear
+a=11:19
+b=a
+for k=1:8
+    b=[b;a+10*k];
+end
+rank(b)
+% 练习8
+y=[]
+for x=-6:0.05:6
+    y=[y difun(x)];
+end
+plot(x,y)
+plot([-6:0.05:6],y)
+% 练习9
+A=[6,2,1,-1;2,4,1,0;1,1,4,-1;-1,0,-1,3];
+b=[6;1;5;-5]
+x=A\b
+% 练习10
+[diag(1:4) eye(4)]
+% 练习11
+num2p(123453656345)
+% 练习12
+ttry(100)
\ No newline at end of file
diff --git a/lab2/dtor.m b/lab2/dtor.m
new file mode 100644
index 0000000..35a6c3c
--- /dev/null
+++ b/lab2/dtor.m
@@ -0,0 +1,7 @@
+% 度分秒转弧度
+function jrad=dtor(deg)
+degd=fix(deg);
+degf=(deg-degd)*100;
+degm=fix(degf);
+degs=(degf-degm)/36.0;
+jrad=(degd+degm*60.00+degs)*pi/180;
\ No newline at end of file
diff --git a/lab2/exercise2.m b/lab2/exercise2.m
new file mode 100644
index 0000000..c02da2e
--- /dev/null
+++ b/lab2/exercise2.m
@@ -0,0 +1,25 @@
+% 练习2第1-4题
+% 1
+A=[6,9,3;2,7,5];
+B=[2,4,1;4,6,8];
+A+B
+A*B
+A*B'
+
+% 2
+clear
+A=[3,2,0,1,5];
+roots(A)
+
+% 3
+clear
+A = [4,9,2;7,6,4;3,5,7];
+B = [37;26;28];
+A\B
+
+% 4
+syms x
+f = 3*x^4-12*x^3-5*x+9;
+int (f,x)
+diff(f,x)
+int(f,x,0,2)
diff --git a/lab2/exercise2.mlx b/lab2/exercise2.mlx
new file mode 100644
index 0000000000000000000000000000000000000000..ffb6f40c8bf4e3131452f2adfd5e04e800411a2d
GIT binary patch
literal 4268
zcmaKv2UHXLw#GvXz4s;{E%eZPFAAZAPNW71HH6+oP<lBCQWOQG1p){h=?E%C5DXwi
zs))d$3qrtyND+MDy}O=Q?zzsa$*h@`{jJ%*y=Tw=gF=XinE(I)8UEk{yzf`9k|zKF
z>WBaUM*JOnU4K*v0u|zH6&ZjCc9e|pL$>6+HSGn{CE%hdhHP`61!$28?!@w>QtDNy
zbuS|esxtiE7brwWDVT_m&=8ZO0>d^6_Vuvqj&6NG-P*@**CJK!6V`<x{A}BA!RCc+
zE%XTO;{sWN+j~K;Zz*6|R5Q)%=lUIN^v-WoSBW8cF(c?T?7Ku@N+!HQZBxD$>S$b<
zXQT@85-Ba>^JL^}w0F?)y*^KcE-HfTElhHxG*#b~=nBXolBM=KRBd<?t+qlUva5oj
z3Sail^?5;#uBoHkp4DWJFlnbYhE#q?l3(0Ai!9L$yjv={FI!a4-@?jFw@J?tQweMg
zk~n_wu9%MBcTCb{naGT2Dgv}T$wY3%3q5lUoIh<jyQW<^@|tYL<8*X-04Jl-YyM>e
zUNp|DcdH(vs4<edsSd>){+L640YZPF+<u_41j19kLJR;<;wd`^A&|jRk{7>qIX!r~
zR9d*RG|)GBUqvA;zi`58#e{y+hkd>^Erz+}$kQm9ScBdYQO1Ec`+o;>d2Z{@WGl+G
zn#)ugs#OvT+LmQmn>s6GHrX515{Z;pGX&q%gJ3&S8a<Ic#{%|>a|~u#3Y9IvCq>dj
zzMs4MjnD^Kbu!#ax8;~8y2GJ-kd42YWp#lPS9J)*Mf6QgL>P<|Ttq{ZzZSAMJ%YZ#
zK|5wAqW7xy%UZi&xcl$lKdi}!UO(<Pt;o^0I~hiQwoN_I(fMln?FVL0o0Hy-ig+v*
zKwA4D7_o;vN>=z&TNSqMOJ#UsA|Q}KCLv6w&c!u2vDnd%RQ3ro$OJ7ILk}X!SRn1z
zGuIdH8qqBA=hyJa3DX{$LHY^LjY)$??CK}=4`_Y`%B8jW=r$fG8V~fMZT;XONVuz%
zyT4ngAHI7p+V-v0xalh}Q}m+kaAyh>v<-~EtF3T%RA70c&mF}*j1==0;~YIt>4Jaf
zX1v{tBL1o^DE)9R$Uh1d<&XP><u(}<x1LGZPgtiJQ*kWGRx3Bpmf{Lfe~zT>c2^Uv
zO-Z6<$0mx7Zp@2v)ib#w`1*PUv4jhyV%bA>#xNGLb_Asrl%<Z1ntVgX@xj&kokMdL
z0p%)-g>f&-V{1Q5ZJ4cV4ns@*wDuN2+lT;~xR|UeE{Jc+Ow|_pmA1TrnGIp<nyY0Q
zO*_?#zHgkcZN~<}=f{bv!O=7OrvYT4oVtvM>#fVm8cD!OB+=ntT7;+}%Xw4Fd`!4X
zappIT*c_<w`@AQoi^r^zdp+%BmJDEhnKE6f=g^vCmt2a4fYFSq;Il3jTka{va*lex
z!BIsd_hAXlyNqHJbN5u?Pr<+yZhG$io=-~!=`|G$L#<yuB7K|X#UC&1$`4JC*)S-J
za7c@f=7LgJc8px7bOwPxc$u3yDK>+b(iaJSt(T4BR!9;F06-)M0MP%emw#wTK<NKg
zYsda=5L$_@b3+8gv)pS%MfhF!s%(p9i?yYEQ}FdL99wPz#ch(St>M1KZdbgVkb<#(
z^(1jh3RbU`U8>C?XY$T%ll`&9ww=N+%aj6pT{Hmv2;S0e8g9jieA@GyM=Kr*M@XGj
zse3ED8&PR!CSjoi`qr4XdVSQ(8%{p(vY1SXbmK(|i=p^0om@AkRg43t`6ziotb$;7
z%L&u`htiFO^0^S1hk8YSNT&rSf@siS4CTYY>&Zc-v}-RLzKkp&N1T^uZ$KirOOLzD
zb_leW`FKcX6ReB_mS!EbAB6_hA(K2^xze!Z1qLq1Z<^zX2Wv^3>y8%ANFQyD@^IuT
zs<R~z^;wk;Ym{=|tmY9AVsG%acum?-$cUk&mA0Nvp3HlbySbR|bf<srg!%O0XFaXw
zWj-MxKsQ-yAip82;3te(VQLxJNPSr%xoq`0<`w)0&1_oyhaTh}PvjR)%_|F}`7U>y
zI&On4yt-HIhz@~FH%=vG6u%Fv6mgblS~3U+TEHu83hTgy@lA_d#G)^K(!;XDTF@0E
zs7<3vF5)Lk^uAZKt<0BRmT9rVs>Fo(OBd+Lt@+5O!UVX82KDq22NReBpAg6kX}UEQ
zef4bEoK<~7pvrpsbC-<sze#63Gdf4Py5SkzCt1lhTwC)X{$8lur=0S1p(Ksar(ul+
zap?*JmWg)@J^`2ziuSNa0UxY~(m$lN<@$a)g|T+dD0R-{A9bquq!teVhV6#8`1Nfo
zYsn4r4a$Y77Iacp6;dp~O{Pp^!@EI!+}xg8erNSLEnQ^lGwbTR!xb}EI{9-W9=r+E
zyuRM+nU1!HIZ5WBD-<YJgxT5pjn55|>WA8%>S6=2B6BxJlPRT|os%AI)pZRw=Qunb
z>7ZK4bJw#g+8*>uG}nTpxksEO;NIwBfC4E@q_ZtK4i5~YbwF6>1y{Yega|*4^cQmM
z7FoJ<-wpK4;Ksw(0}AJq?v4l4I2V<w6&ib>H^%t6q(oOmGaC?jwxRR3GukFI(JOOP
zB@Nr?{;JGoml}>Sdk()rSb?T^?v>uzV%m>O<BoxT0*OU$r`ub<Z|*|&wuB^8CrBx_
zBnvcai!=?<>c9kuYdFM{>(jlu4v7C#8T7~eQ1^~8@*NjicZY8QdSUM2so`Q(y)(5P
zmcqAv*(TwuPxB1DkFEga(X-T<C>&GpHQl=O7R&zks{OMK<ky(&LvO|8xmVf_<bbZy
ztW;%`s#T|y<e5~%7M>%GKBA1aNN00J!JUDYCzBf{*iEG66SvU6S8sw=Aho<?(oup9
zQ3KM$PTt=iza$X3NuHkj(Lr@;C*?&|){wPHt4wvAoxlC~?#!&8@2f-SP6wY?)ADdy
z9eOLw^bkRVkl2gSxiX(;r$nx_Y)?!isdgTe%2z(!Wxo*`@dFTbu-bU7I3{y9lH{@9
zveD}{Dq%63j32-qqP~r0?OTJ1#iWnG24;P8bBh$?&N7d{+6u3~`li>i=Io_)mWZbD
zQWWAczdd?$tNBGjItcrvMJC6XWUxXb8rwhJW0T&+F!3p>%`w@uD+`q$DRa!A<drz6
zk=4$g8ZV%D00L=|>96USpOWe`>x@5xHXA_04#koahM7lFZ5Yd@+P%I$@Y0qV+tO)v
zlbj2LH2X-6Wh8oeIy5((ldS2Hy<XChJ|RVjj+{_v?ubdis=mf0%1S%M=#adX*@E#D
z@d0~B^GaZ^)Nq)OOMq9S$t5ryH7SmL_5>DxHMyJzu@&F0*55fQCp?1DAAltY_YL*d
zHuPIyP8WYN)LrA*o${zWO>;(uT{GB+>(+q?@}~f*WI3{cktfepUMw@F4lYhlwkE!6
zrgv(G_FEh;{P^;n_Sd0pq5<D#!GC1g2>}4ki=mAOfxE*);8Je>K?rD&e*hvV#2XQO
z@v*(1rEL}uW=q1IGMnI_6mHVg_3<9UCdPeUFXY!aeakDnYmh${Xh=~Ve+;ZgcM3aq
zwb|jUlQnu&-fuHZDD{vvZ*o?DV9=Z?*aY%#w4)+sd_nsJTn5`d=48Ek{(p@h0&DY@
zpkVKT1d;ubs9v%_SJJfSN^~j-9~MM+dna1^0`^DjdA+l1QwT>nca@bX_+>fw4d180
zyfu)D6dc+mL?m;nI&MZWc7kU;-BYO|YG`%4wLf3E$>c-{d-|{?MZ7VP{A`BZ?me40
zn!3KVQb)Yx6AKM&l!eDec|il}aUU_Og3n*pC_XdFk;??#-E1_qcu;4T-Nn2*^F7fr
zF|G9r&%0aL#?Qto<U03P&GM#}DSySefRJDx!DBSRBjx^goPOc|3Dj|{#|OVG1TIAt
zy*lbtqE=7;>sfPrQHQkAA>fGwho!wjt9M>o@qnZWYWdoa<a_gH5d*sq8>!{?Og;oA
zmB%-Pa36h(xD~l3&dN(9Nk?<5fB3zr{5#E%>qqx>`mDq~Z1@I=ij7kJb==iUljm=O
zCiu(kQK5TR<2rko^%wnaHhd1A0t<zf#98$#m&LI?o`P;6mQB}FpX9_tn5m8S9bohW
z?pf&tnmy=~+04qb*q2Wse^h8LIcuVI-&Bo6&-hDOh)-aEHKbMRtiOWTUkTQ##sgTu
zA40zap&x;_BEjDNsEZbJ%ur>B12X}lOtTDm;@Id4@7$kexp`+`g&LflV7QDsj8WMM
z54XWx@$gXm0uO5at43}mUA!ZvQMY%rGQLWQk{(Ih)HJy#OfnMYSN4|saV0vO@>t2-
z1$U2{n#awd)dI+c)w)jIx4-*cuzWM!0U!mwH{5hTyX*X%9vdSWQ2@p~ovHRGsJX(2
z)3cf6EVW85lVna|Ah3>YvDh{qlFIiwzC9@<^=q=yY}1Tw_#`FpRR{epSqlUb0S`tP
zqCEUB(hW6Lrp2ebaLjCWXc3_^R-XIpSvy^37Oxer1BPT}G;ffewG9`OnTW$A<=uPo
zd9FJ|`~&eg#?huV(#2`yTBEEks-BuJ;IRd?HPckV@acZaQtrX)x}5uQWPTwyy3`v2
zMgTsy;Lw$0quSGXfJ;;Nc*|I$!~u8nV|DZ&f64~UEl<dyDQBF7tVD+@V-jRfhB{I2
zC}Y@yw1^$AFy0uqxslRrq1zix+onnzZD5OaIanOxfkFuI)8fDTOFRJpp2WYle{h?Z
z1uk2si^%^I@W7w`%}8AiylfGE2Dag!@E;7rWr53U)}I2Nc-QeSfuAeb%i))&=g)91
zir>!C%R-mu#!n$%yqEm{q5Mx%<g(P|H}9uZ2-R<-et8csYhG^fpPF#$-)LTHcPNAe
Rco9j4-<9z7aiY1n`X7FA77PFY

literal 0
HcmV?d00001

diff --git a/lab2/exercise2_5.m b/lab2/exercise2_5.m
new file mode 100644
index 0000000..0fe85fe
--- /dev/null
+++ b/lab2/exercise2_5.m
@@ -0,0 +1,14 @@
+% 练习2第五题
+subplot(1,2,1)
+t=0:0.1:2*pi;
+x=sin(3*t).*cos(t);
+y=sin(3*t).*sin(t);
+plot(x,y,'-r','Linewidth',5);
+subplot(1,2,2)
+y1=sin(x);
+y2=sin(x).*exp(-0.1*x);
+plot(x,y1,'+b',x,y2,'-g')
+title('曲线y1 y2')
+xlabel('变量x')
+ylabel('变量y')
+legend('y1','y2');
\ No newline at end of file
diff --git a/lab2/exercise2_6.m b/lab2/exercise2_6.m
new file mode 100644
index 0000000..b86fffe
--- /dev/null
+++ b/lab2/exercise2_6.m
@@ -0,0 +1,16 @@
+x=1:10;
+y=[4.9,3.3,4.2,4.2,5.7,6.5,7.7,9.4,14.7,19.0];
+a=polyfit(x,y,2);
+b=polyfit(x,y,3);
+c=polyfit(x,y,4);
+xx=1:0.1:10;
+yy1=polyval(a,xx);
+yy2=polyval(b,xx);
+yy3=polyval(c,xx);
+plot(x,y,'*b')
+hold on
+plot(xx,yy1,'-r',xx,yy2,'-g',xx,yy3,'-y')
+title('2,3,4次曲线拟合')
+xlabel('变量x')
+ylabel('变量y')
+legend('y','yy1','yy2','yy3')
\ No newline at end of file
diff --git a/lab2/exercise2_7.m b/lab2/exercise2_7.m
new file mode 100644
index 0000000..9d2441a
--- /dev/null
+++ b/lab2/exercise2_7.m
@@ -0,0 +1,16 @@
+subplot(1,2,1);
+x=-3:0.1:3;
+y=-3:0.1:3;
+z=x.*exp(-x.^2-y.^2);
+plot3(x,y,z);
+subplot(1,2,2);
+[x,y]=meshgrid(-3:0.1:3);
+z=x.*exp(-x.^2-y.^2);
+surf(x,y,z);
+shading interp;
+colormap("spring");
+k=find(x<=0 & y<=0);
+z1=z;
+z1(k)=NaN;
+surf(x,y,z1);
+shading interp;
\ No newline at end of file
diff --git a/lab2/exercise2_9.m b/lab2/exercise2_9.m
new file mode 100644
index 0000000..ec1c40e
--- /dev/null
+++ b/lab2/exercise2_9.m
@@ -0,0 +1,3 @@
+axes(handles.axes1);
+[X,Y,Z]=sphere(25);
+mesh(X,Y,2*(Z+1));
diff --git a/lab2/fn.m b/lab2/fn.m
new file mode 100644
index 0000000..818cff5
--- /dev/null
+++ b/lab2/fn.m
@@ -0,0 +1,20 @@
+function t=fn(n);
+m=3;
+t=1:n;
+a=1;
+j=1;
+ind=0;
+while a<n
+    while j<m
+        ind =ind + 1;
+        if ind > length(t)
+            ind=1;
+        end
+        j=j+1;
+    end
+    t(ind)=[]
+    ind=ind-1;
+    j=1;
+    a=a+1;
+end
+t
diff --git a/lab2/pmzbfs.m b/lab2/pmzbfs.m
new file mode 100644
index 0000000..0bef890
--- /dev/null
+++ b/lab2/pmzbfs.m
@@ -0,0 +1,10 @@
+function [raf,jl]=pmzbfs(x1,y1,x2,y2)
+detx=x2-x1;dety=y2-y1;
+faw1=atan(dety/detx);
+j1=sqrt(detx*detx+dety*dety);
+if detx<0
+    raf=faw1+pi;
+end
+if (detx>0)*(det<0)
+    raf=faw1+2*pi;
+end
\ No newline at end of file
diff --git a/lab2/rad2deg.m b/lab2/rad2deg.m
new file mode 100644
index 0000000..b7a318f
--- /dev/null
+++ b/lab2/rad2deg.m
@@ -0,0 +1,7 @@
+function jdeg=rad2deg(rad)
+rad=rad*180/pi;
+deg=fix(rad);
+xsbf=(rad-deg)*60;
+degf=fix(xsbf);
+degm=(xsbf-degf)*60;
+jdeg=deg+degf/100.0+degm/10000.0;
\ No newline at end of file
diff --git a/lab3/exercise3_1.m b/lab3/exercise3_1.m
new file mode 100644
index 0000000..6d4344a
--- /dev/null
+++ b/lab3/exercise3_1.m
@@ -0,0 +1,8 @@
+% 第一题
+
+A=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];
+B=[4;-3;9;-8];
+% 系数矩阵的秩
+rank(A)
+% 方程的解
+C = A\B
\ No newline at end of file
diff --git a/lab3/exercise3_2.m b/lab3/exercise3_2.m
new file mode 100644
index 0000000..6f4b642
--- /dev/null
+++ b/lab3/exercise3_2.m
@@ -0,0 +1,6 @@
+% 练习3 第二题
+x=-10:1:10;
+y=-10:1:10;
+[X,Y]=meshgrid(x,y);
+Z=sin(sqrt(X.^2+Y.^2))./(sqrt(X.^2+Y.^2));
+surf(X,Y,Z)
\ No newline at end of file
diff --git a/lab3/exercise3_3.m b/lab3/exercise3_3.m
new file mode 100644
index 0000000..c150458
--- /dev/null
+++ b/lab3/exercise3_3.m
@@ -0,0 +1,7 @@
+num=[1];
+den=[1 2 2 1];
+sys=tf(num,den);
+figure(1); pzmap(sys); title('零极点图');
+figure(2); impulse(sys);
+w=0:0.1:10*pi;
+figure(3); freqs(num,den,w)
\ No newline at end of file
diff --git a/lab3/exercise3_4.m b/lab3/exercise3_4.m
new file mode 100644
index 0000000..8768c19
--- /dev/null
+++ b/lab3/exercise3_4.m
@@ -0,0 +1,10 @@
+subplot(2,2,2);
+x=-pi/2:pi/10:pi/2;
+y=sqrt(cos(x));
+plot(x,y);
+subplot(2,2,3);
+x=-2:0.1:2;
+y=-4:0.1:4;
+[x,y]=meshgrid(x,y);
+z=x.^2/4+y.^2/16;
+surf(x,y,z)
\ No newline at end of file
diff --git a/lab3/lab3.m b/lab3/lab3.m
new file mode 100644
index 0000000..1b5841b
--- /dev/null
+++ b/lab3/lab3.m
@@ -0,0 +1,42 @@
+% 实验3
+% 读取数据
+[filename1,p1]=uigetfile('*.txt', '选择尼康全站仪坐标文件名');
+fp=fopen(strcat(p1,filename1),'r');
+i=1;bss=[];
+while ~feof(fp)
+    line=fgetl(fp);
+    data=findstr(line,',');
+    bss(i).dh=line(1:data(1)-1); 
+    bss(i).x=str2num(line(data(1)+1:data(2)-1));
+    bss(i).y=str2num(line(data(2)+1:data(3)-1));
+    bss(i).h=str2num(line(data(3)+1:data(4)-1));
+    bss(i).dm=line(data(4)+1:end);
+    i=i+1;
+end
+fclose(fp);
+
+% 计算国家坐标系的坐标
+a=3976223.453;b=39512553.524;rad=deg2rad(0.433267);
+c=cos(rad);d=-sin(rad);lmd=1.000034336;
+m=i-1;gjzb=[]; 
+for i=1:m
+    gjzb(i).dh=bss(i).dh;
+    gjzb(i).dm=bss(i).dm;
+    gjzb(i).h=bss(i).h;
+    gjzb(i).x=a+lmd*(c*bss(i).x-d*bss(i).y);
+    gjzb(i).y=b+lmd*(d*bss(i).x+c*bss(i).y);
+end
+
+% 输出数据
+[filename2,p2]=uigetfile('*.txt', '输出文件');
+fn=fopen(strcat(p2,filename2),'w');
+fprintf(fn,'%s\n\n','   点名       X(m)        Y(m)     H(mm)   code');
+fclose(fp);fn=fopen(strcat(p2,filename2),'a');
+for i=1:m
+    fprintf(fn,'%8s', gjzb(i).dh);
+    fprintf(fn,'%15.3f', gjzb(i).x);
+    fprintf(fn,'%16.3f', gjzb(i).y);
+    fprintf(fn,'%10.3f', gjzb(i).h);
+    fprintf(fn,'%8s\n', strcat('   ',gjzb(i).dm));
+end
+fclose(fn);
\ No newline at end of file
diff --git a/lab3/nikon.txt b/lab3/nikon.txt
new file mode 100644
index 0000000..fa7a74a
--- /dev/null
+++ b/lab3/nikon.txt
@@ -0,0 +1,7 @@
+3,49510.3110,80362.3190,405.8240,s50
+63,49488.9080,80426.0080,408.8490,S77
+164,49461.2960,80462.3610,408.9950,S77
+165,49457.2500,80485.8770,409.4490,S81
+166,49476.6650,80395.8870,408.1090,S66
+167,49475.1380,80398.7300,408.1570,S77
+298,49465.2360,80593.9700,415.1290,S99
\ No newline at end of file
diff --git a/lab3/zhuanhuan.txt b/lab3/zhuanhuan.txt
new file mode 100644
index 0000000..f9ff4c0
--- /dev/null
+++ b/lab3/zhuanhuan.txt
@@ -0,0 +1,9 @@
+   点名       X(m)        Y(m)     H(mm)   code
+
+       3    4026749.416    39592285.024   405.824     s50
+      63    4026728.820    39592348.981   408.849     S77
+     164    4026701.670    39592385.683   408.995     S77
+     165    4026697.922    39592409.249   409.449     S81
+     166    4026716.196    39592319.017   408.109     S66
+     167    4026714.706    39592321.879   408.157     S77
+     298    4026707.277    39592517.236   415.129     S99
diff --git a/mat2dms.m b/mat2dms.m
new file mode 100644
index 0000000..8172cb2
--- /dev/null
+++ b/mat2dms.m
@@ -0,0 +1,5 @@
+function dmsv=mat2dms(D)
+[l,c]=size(D);
+for i=1:l
+    dmsv(i)=D(i,1)+D(1,2)/60+D(i,3)/3600;
+end
\ No newline at end of file
diff --git a/num2p.m b/num2p.m
new file mode 100644
index 0000000..800ed33
--- /dev/null
+++ b/num2p.m
@@ -0,0 +1,8 @@
+function a=num2p(n)
+strr=num2str(n);
+index=size(strr);
+indexx=index(2);
+a=[];
+for k=indexx-1:-1:0
+    a=[a rem(fix(n/(10^k)),10)];
+end
diff --git a/qwe.mat b/qwe.mat
new file mode 100644
index 0000000000000000000000000000000000000000..ecb75a862600b78ba8efd830613c0e4d1f88bfcc
GIT binary patch
literal 375
zcmeZu4DoSvQZUssQ1EpO(M`+DN!3vZ$Vn_o%P-2cQV4Jk_w+L}(NS<NN=+<DO;O0t
zvr=#?%2e>nFH=x3S1>fRGBL6;HdHV&Ffvsj5-`93qo*%Fko_KrE9N{-PDtQLNJ&Uw
zz9hpVAjJ}~u*OdzjWNlMEv;?I5*34Z#uass4^Iy^-y2vkW1{+%ng>^w|GjMXQDj}^
z+R_6~B})%(;}zm--4=6QJb!}4hg}t`U%bD3if_3S-@1sRzV8QihCW~4c2vUB;VA!;
z<}KQ<td3;Q6SK%mEmEkMl@y=(@cZ97arI+=vZoa`Tw<TaZ2R17#YedZ8|KPd?T^0N
zm_PrZ&synKd;x{*Y)lJQzdxz|)BFB+Z@Vsm`*neg@API^e^Jx<k@(}ybX&#Ze;y+7
j;zmIw2JcIJDykeMp6Ttdma_SF{nR@FI|c^5t$FPLCg6+)

literal 0
HcmV?d00001

diff --git a/rad2dms.m b/rad2dms.m
new file mode 100644
index 0000000..6deb25b
--- /dev/null
+++ b/rad2dms.m
@@ -0,0 +1,2 @@
+function radv=rad2dms(dv)
+radv=rem(dv*pi/180,2*pi);
\ No newline at end of file
diff --git a/tiaojianpingcha.m b/tiaojianpingcha.m
new file mode 100644
index 0000000..d27fec5
--- /dev/null
+++ b/tiaojianpingcha.m
@@ -0,0 +1,31 @@
+clc
+disp('条件平差示例2')
+disp('三角形内角观测值')
+L1=[42,12,20]
+L2=[78,9,9]
+L3=[59,38,40]
+L=[L1;L2;L3]
+disp('将角度转换为弧度')
+% LL=dms2rad(mat2dms(L));
+% A=[1,1,1];
+% w=sum(LL(:))-pi;
+% w=dms2mat(rad2dms(w));
+% P=eye(3);
+% Naa=A*inv(P)*A';
+% Ka=-inv(Naa)*w;
+% V=A'*Ka;
+% L1=L+V;
+% LL=dms2rad(mat2dms(L1));
+% sumLL=sum(LL);
+% if (abs(sum(LL)-pi)<1e-4)
+%     disp('正确')
+% else
+%     disp('错误')
+% end
+
+SumLL=tjpc(L);
+if (abs(SumLL-pi)<1e-4)
+    disp('正确')
+else
+    disp('错误')
+end
\ No newline at end of file
diff --git a/tjpc.m b/tjpc.m
new file mode 100644
index 0000000..7ff618c
--- /dev/null
+++ b/tjpc.m
@@ -0,0 +1,12 @@
+function sumLL=tjpc(L)
+LL=dms2rad(mat2dms(L))
+A=[1,1,1];
+w=sum(LL(:))-pi;
+w=dms2mat(rad2dms(w));
+P=eye(3);
+Naa=A*inv(P)*A';
+Ka=-inv(Naa)*w;
+V=A'*Ka;
+L1=L+V;
+LL=dms2rad(mat2dms(L1));
+sumLL=sum(LL);
\ No newline at end of file
diff --git a/ttry.m b/ttry.m
new file mode 100644
index 0000000..6f2b3d6
--- /dev/null
+++ b/ttry.m
@@ -0,0 +1,7 @@
+function ppi=ttry(n)
+pointx=-1+2*rand(1,n);
+pointy=-1+2*rand(1,n);
+index=pointx.*pointx+pointy.*pointy
+indexx=index<=1;
+ppi=4*sum(indexx)/n
+end
diff --git a/untitled3.m b/untitled3.m
new file mode 100644
index 0000000..6d62bec
--- /dev/null
+++ b/untitled3.m
@@ -0,0 +1,10 @@
+A=[1,1,1];
+L=[1,2,2];
+H=[1.003,0.501,-0.503];
+w=7;
+P=diag(2./L);
+P1=inv(P);
+N=A*P1*A';
+Ka=7/N*8;
+V=P1*A'*Ka;
+h=V'+H
\ No newline at end of file
diff --git a/untitled5.m b/untitled5.m
new file mode 100644
index 0000000..7658ec3
--- /dev/null
+++ b/untitled5.m
@@ -0,0 +1,30 @@
+clc
+clear
+h1=-1.004;
+h2=1.516;
+h3=2.512;
+h4=1.520;
+HA=12.013;
+HB=10.013;
+h=[h1,h2,h3,h4]'
+s1=2;
+s2=1;
+s3=2;
+s4=1.5;
+s=[s1,s2,s3,s4]'
+A=[1,1,-1,0;0,1,1,-1];
+w1=h1+h2-h3+Ha-HB;
+w2=h2-h4;
+w=[w1;w2];
+P=diag(1./S)
+Naa=A*inv(P)*A';
+Ka=-inv(P)*A'*Ka;
+H=h+V;
+if H(1,1)+H(2,1)-H(3,1)+HA-HB == 0 && H(2,1)-H(4,1)==0
+    disp('正确')
+else
+    disp('错误')
+end
+disp('平差结果');
+HC=HA+H(1,1)
+HD=HA+H(1,1)+H(4,1)
\ No newline at end of file
diff --git a/veryfun.m b/veryfun.m
new file mode 100644
index 0000000..9b2055e
--- /dev/null
+++ b/veryfun.m
@@ -0,0 +1,16 @@
+function f=veryfun(n)
+if n==6174
+    return
+else
+    index(1)=(n-rem(n,1000))/1000;
+    index(2)=fix(n/100)-index(1)*10;
+    index(3)=rem(fix(n/10),10);
+    index(4)=rem(n,10);
+    index2=sort(index,'ascend');
+    index3=sort(index,'descend');
+    a=[1000 100 10 1];
+    max1=index3*(a');
+    min1=index2*(a');
+    max1-min1
+    veryfun(max1-min1)
+end
\ No newline at end of file