水准网平差实例.pdfVIP

水准网平差实例--第1页

functiontiaojian

%条件方程的组成

A=[1-100100;001-1100;0010011;010-1000];

W=[786-3];

%协因数阵的组成

Q=diag([1.11.72.32.72.41.42.6]);

P=inv(Q);

%法方程组成

Naa=A*Q*A;

%求解联系数K和改正数

K=-inv(Naa)*W;

V=Q*A*K;

%计算观测值的平差值

h=[1.3592.0090.3631.0120.6570.238-0.595];

adjh=h+V/1000;

%精度评定

VTPV=V*P*V;

r=4;

var0=sqrt(VTPV/r);%计算单位权方差

%评定观测值平差值函数的精度

fs=[0000100];%函数的系数

Qss=fs*Q*fs-(A*Q*fs)*inv(Naa)*A*Q*fs;%观测值平差值函数的协因数

varS=Qss*var0;%观测值平差值函数的方差

%写输出结果文件

fid=fopen(陈岑桌面 iaojian.dat,w+);

fprintf(fid,%s

,######水准网平差算例######);

fprintf(fid,%s,观测值改正数);fprintf(fid,

);

fori=1:7

fprintf(fid,%s%1.0f,V,i);

fprintf(fid,%15.3f

,V(i));

end

fprintf(fid,

);

fprintf(fid,%s

,观测值平差值);

fori=1:7

fprintf(fid,%s%1.0f,h,i);

fprintf(fid,%15.3f,h(i));

fprintf(fid,

);

end

fprintf(fid,

);

水准网平差实例--第1页

水准网平差实例--第2页

fprintf(fid,%s,单位权中误差);fprintf(fid,

);

fprintf(fid,%s%3.2f,m0=,var0);fprintf(fid,

);

fclose(fid);

######水准网平差算例######

观测值改正数

V1-0.243

V22.855

V3-4.243

V4-0.145

V5-3.902

V6-0.615

V7-1.142

观测值平差值

h11.359

h22.009

h30.363

h41.012