MATLAB水准网间接平差课程作业程序.docx

条件平差A-(1.-1OO1O0;0O1-11O0;0O1OO11;01O-1OO0);Q=diag(三);W=zeros(4,1);W(1.)-h(1.)-h(2>fh(5);W<2)三h(3)-h(4>÷h(5);W(3)=h(3)+h(6>+h(7);W(4)=h(2)-h(4)+H(1.)-H(2);W-W.*1000;Naa-A*Q(A,>K=-1.*inv(Naa)*W;V=Q*(A,)*K;1.-h,+V.1000;x-zeros(1/3);X<1)=H(1)÷1.(1);X(2)=H(1)÷1.(2);X(3)-H(2)-1.(7);ZWC-sqrt(V,Inv(Q)*V4);Q1.1.=Q-QA'*inv(Naa)*A*Q;ZWC_h5=ZWC*sqrt(Q1.1.(5,5);间由差h1.-1.3S9;h2=2.009;h3=0.363;h4-1.012;h5-0.6S7;h6=O.230;h7=-0.595;H1-5.016H26.016h=(h1.h2h3Mh5h6h7)s=(1.1.1.72.32.72.41.426,B-(1.OO;010;1O0;01O;-11O;-1O1;0O-1p=diag(1.s)1-(0;0;4;3;7;2;0)W-B,*p*1.Nbb=B,*p*Bx=inv(Nbb)*WV-(B*x-1.)H-h+V/1000Q=inv(Nbb)n=7；t-3;j-V,p*Vd=Sqrt(j4)f=-1.10,q-f*Q*fD-d*sqrt(q)D1.=d*sqrt(Q)disp1观测值单位m)h1.-1.359;h2-2.009;h3=0.363;h4=1.012;h5-0.657;h6-0238;h7=-0.595;H1.=5.016;%1.2如点高程,单位mH2-6.016巳知点而租.单(h-(h1.h2h3h4h5h6h7)S=1,11.72.32.72.41.42.6S是戏跖长度d工SPr系数矩阵B、1,)B-(1OO;010;100;010;-110;-101;00-1p=diag(1,s定义权阵1-(0;0;4;3;7;2;0)R-Bp1.Nbb=B,*p*Bdispr参数的解，)×-inv(Nbb)*W;v-(B*x-1.)%误差方程(mm)H=h+V1000%观测值的平茶俏Q=inv(Nbb);3观测值协因数阵j-V*p*V;)计算单位权中误差的聚数d-sqrt(j4);、曲位权中误差f=1.-1.10)'权乐数系数住q=f*Q*f;权由教协因数阵dispr输出结果,)H-h+v/1000D=d*sqrt(q)Dc=d*sqrt(Q(1.r1)Dd-d*sqrt(Q(2,2)