大地测量学实验报告

1、大地测量学实验报告专业班级：测绘1103学号：0145110608 姓名：高啥一、采用的数据(1)已知点数据:BJ54:XYZA10-2570443.770964766022.434163359186.13435A09-2576999.808184762622.280873359019.58617A08-2563333.153004770217.784773358671.68519XA80:XYZA10-2600875.540004961731.345393058234.18399A09-260 7588.261034958655.663003057611.15590A08-2593513.9

2、98584965560.876123058204.67725(2)未知点数据BJ54:XYZA-2570448.77104766033.43423359199.1343二、计算过程已知数据导入；计算矩阵系数B;利用最小二乘求七参数值；求残差值；精度评定：求单位权中误差。，七参数的中误差；带入未知点数据，计算未知点的转换坐标。三、计算结果转换参数为：X参数为：-6. 4043e-08Y参数为：1. 2219e-07参数为：7. 5318e-08X旋转参数为：QgY旋转参数为:(L08E旋转参数为：。.鹏参数为:5. 5642e-07参数的精度为：X参数精度为：1. 1849e-30Y参数精度为：

3、1.4152e-30E参数精度为：1.8618e-30X旋转参数精度为：7.6999e-12Y流转参数精度为：2.6478e-U部转参数精度为：1.3156-11m精度为：4. 7604e-17.未知点坐标转换结果：计算未知点G明坐标：未知点GS4坐标：-对也44&苜 14766033. 43423359199. 1343未知点巽口坐标：-26口口881.口3爵4961742. 85522915262. 4823四、实验心得：我通过编程，既了解了坐标转换的内容，也加强了编程的本领。在编程的过 程中，我开始用reshape函数变换矩阵，没有得到想要的结果，后来在同学的 帮助下，用循环函数代替，得

4、到了相应的结果。这次实验我受益匪浅。附录：程序列表：clc;format long g;%显示小数点后12的数据，并用小数位表示load shuju1.txt;load shuju2.txt;A=shuju1;for i=1:3x=A(i,1);y=A(i,2);z=A(i,3);A1(i,:)=1 0 0 x 0 -z y;A2(i,:)=0 1 0 y z 0 -x;A3(i,:)=0 0 1 z -y x 0;L(3*i-2:3*i,:)=shuju2(i,1);shuju2(i,2);shuju2(i,3);endB=-A1(1,:);A2(1,:);A3(1,:);A1(2,:);A

5、2(2,:);A3(2,:);A1(3,:);A2(3,:);A3(3,:);K=-(pinv(B*B)*B*L;dX=K(1);dY=K (2);dZ=K(3);a1=K 4);a2=K(5);a3=K(6);a4=K(7);m=a1-1;vX=a2/a1;vY=a3/a1;vZ=a4/a1;%计算各个已知像点的残差V=B*K+L;m0=sqrt(V*V/2);%计算单位权中误差Q=pinv(B*B);qX=Q(1,1);qY=Q(2,2);qZ=Q(3,3);qm=Q(4,4);qvX=-a2/(a1八2) 1/a1*Q(4,4) Q(4,5);Q(5,4) Q(5,5)*-a2/(a1八

6、2);1/a1;qvY=-a3/(a1八2) 1/a1*Q(4,4) Q(4,6);Q(6,4) Q(6,6)*-a3/(a1八2);1/a1;qvZ=-a4/(a1八2) 1/a1*Q(4,4) Q(4,7);Q(7,4) Q(7,7)*-a4/(a1八2);1/a1;mX=m0*qX;mY=m0*qY;mZ=m0*qZ;mm=m0*qm;mvX=m0*qvX;mvY=m0*qvY;mvZ=m0*qvZ;disp(转换参数为：)；disp(X 参数为：,num2str(dX);disp(Y 参数为：,num2str(dY);disp(Z 参数为：,num2str(dZ);disp(X 旋转参

7、数为：,num2str(vX);disp(Y 旋转参数为：,num2str(vY);disp(Z 旋转参数为：,num2str(vZ);disp(m 参数为：,num2str(m);disp(参数的精度为：)；disp(X 参数精度为：,num2str(mX);disp(Y 参数精度为：,num2str(mY);disp(Z 参数精度为：,num2str(mZ);disp(X 旋转参数精度为：,num2str(mvX);disp(Y 旋转参数精度为：,num2str(mvY);disp(Z 旋转参数精度为：,num2str(mvZ);disp(m 精度为：,num2str(mm);disp(计算未知点G80坐标：)；disp(未知点G54坐标：,num2str(-2570448.77096,4766033.43416,3359199.13435);E=