测量平差课程设计最终版.doc

测 量 平 差 课 程 设 计 报告书 姓名:梁永健 学号:201131080116 华南农业大学信息（软件）学院 《误差理论与测量平差基础课程设计》成绩单 班级 11测绘一班 学号 201131080116 姓名 梁永健 课 程 设 计 内 容 1、用条件平差计算： 5.1.06 有水准网如图所示，其中A B C 三点高程未知，现在其间进行了水准测量，测得高差及水准路线长度为： h1=1.335m S1=2km ;h2 = 1.057m ; S2 = 2km h3=-2.396m S3=2km 试求各高差的平差值。 解：由题目可得：n=3，t=2, 故r=1，可列出如下一个条件方程： 以?i=hi+vi 代入上式， 经计算可得条件方程的最后形式为： 式中闭合差的单位是mm。 令1km的观测高差为单位权观测值，即，于是有 ，，， 法方程系数阵为： 由此得法方程为： 解得k=1，代入改正数方程计算V 得： 观测量的平差值为： ，， 用平差值重新列出平差值条件方程，得： 经检核计算无误，最后计算出各高差平差值为： ，， 2、代码如下 #includestdio.h main() { double h1=1.335,h2=1.055,h3=-2.396,c=1.0,s1=2.0,s2=2.0,s3=2.0; double A[1][3]={{1,1,1}},Q[3][3]={{s1/c,0,0},{0,s2/c,0},{0,0,s3/c}}; double w[1][1],v[3][1],k[1][1],Naa[1][1],INaa[1][1],AT[3][1],m[1][3],n[3][1]; double L1,L2,L3; int i,j; printf(本题n=3,t=2,r=1\n); printf(条件方程：\nh1+h2+h3=0\n即：v1+v2+v3-0.006=0\n); printf(系数矩阵A[%f %f %f],A[0][0],A[0][1],A[0][2]); printf(\n); printf(各高差的权为：p1=c/s1=%f,p2=c/s2=%f,p3=c/s3=%f\n,c/s1,c/s2,c/s3); printf(协因数阵Q为：\n); for(i=0;i3;i++) { if(i==1||i==2) printf(\n); for(j=0;j3;j++) printf(%f ,Q[i][j]); } printf(\n); w[0][0]=h1+h2+h3; for(i=0;i1;i++) { for(j=0;j3;j++) AT[j][i]=A[i][j]; } for(i=0;i1;i++) { for(j=0;j3;j++) m[i][j]=A[0][0]*Q[i][j]+A[0][1]*Q[i+1][j]+A[0][2]*Q[i+2][j]; } Naa[0][0]=m[0][0]*AT[0][0]+m[0][1]*AT[1][0]+m[0][2]*AT[2][0]; printf(Naa=%f\n,Naa[0][0]); INaa[0][0]=1/Naa[0][0]; printf(Naa的逆：\nINaa=%f\n,INaa[0][0]); k[0][0]=-INaa[0][0]*w[0][0]; printf(联系数k=%f\n,k[0][0]); for(i=0;i3;i++) { for(j=0;j1;j++) n[i][j]=Q[i][j]*AT[0][0]+Q[i][j+1]*AT[0][1]+Q[i][j+2]*AT[0][2]; } for(i=0;i3;i++) { v[i][0]=n[i][0]*k[0][0]; } printf(改正数v：\nv1=%f m\nv2=%f m\nv3=%f m\n,v[0][0],v[1][0],v[2][0]); L1=h1+v[0][0]; L2=h2+v[1][0]; L3=h3+v[2][0]; if((L1+L2+L3)-0-1e-6(L1+L2