大地测量学实验报告.docVIP

大地测量学实验报告 专业班级：测绘1103 学号：0145110608 姓名：高晗 BJ54: X Y Z A10 - 2570443.77096 4766022.43416 3359186.13435 A09 - 2576999.80818 4762622.28087 3359019.58617 A08 - 2563333.15300 4770217.78477 3358671.68519 XA80: X Y Z A10 - 2600875.54000 4961731.34539 3058234.18399 A09 - 260 7588.26103 4958655.66300 3057611.15590 A08 - 2593513.99858 4965560.87612 3058204.67725 （2）未知点数据 BJ54: X Y Z A -2570448.7710 4766033.4342 3359199.1343 二、计算过程 1. 已知数据导入； 2. 计算矩阵系数B； 3. 利用最小二乘求七参数值； 4. 求残差值； 5. 精度评定：求单位权中误差σ，七参数的中误差； 6. 带入未知点数据，计算未知点的转换坐标。 三、计算结果 1. 4.未知点坐标转换结果： 四、实验心得： 我通过编程，既了解了坐标转换的内容，也加强了编程的本领。在编程的过程中，我开始用reshape函数变换矩阵，没有得到想要的结果，后来在同学的帮助下，用循环函数代替，得到了相应的结果。这次实验我受益匪浅。 附录：程序列表： clc; format long g;%显示小数点后12的数据，并用小数位表示 load shuju1.txt; load shuju2.txt; A=shuju1; for i=1:3 x=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)]; end B=-[A1(1,:);A2(1,:);A3(1,:); A1(2,:);A2(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^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旋转参数为：,num2str(vX)]); disp([Y旋转参数为：,num2str(vY)]); disp([Z旋转参数