程序及运算结果.docxVIP

MATLAB程序展示——水准网间接平差function [v,ch,cx]=szw(s,h,B,x0,d,n,t,f)%此函数用于计算改正数，高差中误差，高程中误差p=diag(1./s);%定义权阵disp(P=)disp(p)l=h-B*x0-d;W=B*p*l;Nbb=B*p*B;disp(NBB=)disp(Nbb)x=(inv(Nbb))*W;disp(改正数)v=B*x-l;%改正数disp(v);c0=sqrt((v*p*v)/(n-t));%单位中误差Nbb=B*p*B;Qh=f*(inv(Nbb))*f;%h5的协因数阵disp(高差平差值中误差)ch=c0*sqrt(Qh);%高差平差值中误差disp(ch);Nbbn=inv(Nbb);%求逆矩阵disp(高程平差值中误差)cx=c0*sqrt(diag(Nbbn));%高程平差值中误差disp(cx);returnload B.txt;load d.txt;load s.txt;load x0.txt;load h.txt;load n.txt;load t.txtload f.txt[v,ch,cx]=szw(s,h,B,x0,d,n,t,f);计算结果展示MATLAB程序展示——导线网间接平差load n.txt%观测总数load t.txt%必要观测数load dws.txt%导入点位数load csgs.txt%导入参数个数load bs.txt%导入待测边数load dcdh1.txt%%第一个待测点号load x1y1.txt%导入已知点位一的坐标load xydws.txt%导入末点位的坐标load deg.txt%导入由观测角计算而得的方位角load s.txt%导入边长观测值load gcj.txt%各点位观测角digits(12);digits(15);x1=x1y1(1,1);y1=x1y1(1,2);x0=zeros(dws,1);y0=zeros(dws,1);a=zeros(1,3);x0(1)=x1;y0(1)=y1;for i=2:(dws-1) a=deg(i-1,:); [x2,y2]=zbzs(x1,y1,s(i-1),a); x1=x2; y1=y2; x0(i)=x1; y0(i)=y1;endx1=xydws(1,1);y1=xydws(1,2);x0(dws)=x1;%各点位x的近似值y0(dws)=y1;%各点位y的近似值s0=zeros(dws-1,1);%各各点位之间距离的近似值for i=1:(dws-1) s0(i)=sqrt((x0(i)-x0(i+1)).^2+(y0(i)-y0(i+1)).^2);end[B1,L1]=jsB1L1(x0,y0,s0,s,bs,dcdh1,csgs);[B2,L2,ajs]=jsB2L2(x0,y0,s0,dws,dcdh1,csgs,gcj,bs);B=[B1;B2];L=[L1;L2];juw=5;%测角中误差buw=0.5*sqrt(s);%测边中误差P=zeros(1,length(gcj));for i=1:length(gcj) P(i)=1./(juw.^2);endP=[1./(buw.^2),P];P=diag(P)*(juw.^2);NBB=B*P*B;W=B*P*L;x=pinv(NBB)*W;disp([xE yE xF yF]=)disp(x)X=zeros(1,length(x));X=[x(1)+x0(2) x(2)+y0(2) x(3)+x0(3) x(4)+y0(3)];disp(E点坐标平差值=)disp(X(1,1:2))disp(F点坐标平差值=)disp(X(1,3:4))V=B*x-L;uwef=zwc(V,P,n,t,B);%各点位精度disp(点位精度uwe uwf=)disp(uwef)disp(各观测量改正数[Vb1 Vb2 Vb3 Vb4 Vj1 Vj2 Vj3]=)disp(V)spcz=s+V(1:3,1);%边长观测值平差值disp(边长观测平差值S1 S2 S3=)disp(spcz)gcjpcz=gcj+V(4:7,1);disp(角度观测平差值J1 J2 J3 J4=)disp(gcjpcz)function A=deg2rad(a)%角度转弧度A=a*[pi/180;pi/10800;pi/648000];Returnfunction [x2,y2]=zbzs(x1,y1,s,a)%坐标正算A=deg2rad(a);x2=x1+s*cos(A);y2=y1+s*sin(A);Returnfunction [B1,L1]=jsB1L1(x0,y0,s0,s,bs,dcdh