水准网平差实例.doc

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('C:\DocumentsandSettings\陈岑\桌面\tiaojian.dat','w+');
fprintf(fid,'%s\n\n','######水准网平差算例######');
fprintf(fid,'%s','观测值改正数');fprintf(fid,'\n');
fori=1:7
fprintf(fid,'%s%1.0f','V',i);
fprintf(fid,'%15.3f\n',V(i));
end
fprintf(fid,'\n');

fprintf(fid,'%s\n','观测值平差值');
fori=1:7
fprintf(fid,'%s%1.0f','h',i);
fprintf(fid,'%15.3f',h(i));
fprintf(fid,'\n');
end
fprintf(fid,'\n');

fprintf(fid,'%s','单位权中误差');fprintf(fid,'\n');
fprintf(fid,'%s%3.2f','m0=',var0);fprintf(fid,'\n');
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
h50.657
h60.238
h7-0.595

单位权中误差
m0=2.22

--------------------------------------------------------------------------------------
间接平差（水准网）

function[B,l,P,h]=Finput

%误差方程的组成
B=[10;01;10;01;-11;-10];
l=[004372]';

%协因数阵Q的组成
Q=diag([1.11.72.32.72.44.0]);
P=inv(Q);

%观测值
h=[1.3592.0090.3631.0120.657-0.357]';

functionjianjie1

[B,l,P,h]=Finput;

%fafangchengdezucheng
Nbb=B'*P*B;

%adjxV
adjx=inv(Nbb)*B'*P*l;
V=B*adjx-l;

%计算观测值的平差值
adjh=h+V/1000;

%精度评定
VTPV=V'*P*V;
r=4;
var0=sqrt(VTPV/r);%计算单位权方差

%写输出结果文件
fid=fopen('C:\DocumentsandSettings\陈岑\桌面\jianjie1.dat','w+');
fprintf(fid,'%s\n\n','######水准网间接平差算例######');
fprintf(fid,'%s','观测值改正数');fprintf(fid,'\n');
fori=1:6
fprintf(fid,'%s%1.0f','V',i);
fprintf(fid,'%15.3f\n',V(i));
end
fprintf(fid,'\n');

fprintf(fid,'%s\n','观测值平差值');
fori=1:6
fprintf(fid,'%s%1.0f','adjh',i);
fprintf(fid,'%15.3f',adjh(i));
fprintf(fid,'\n');
end
fprintf(fid,'\n');

fprintf(fid,'%s','单