测绘工程matlab论文

课程编号：                      课程性质： MATLAB及其应用 课 程 论 文 学院：  测绘学院  专业：  测绘工程  姓名：               学号：                目录 一、题目内容    3 二、程序编写思路    4 1、数学模型    4 2、程序解析    4 三、程序运行结果    5 四、程序源代码    6 注：这是《误差理论与测量平差基础》这一课程在学习其基本原理之后为解决实际问题而编写的一个程序。 一、题目内容 要求：对于给出的导线控制网图2 ，根据已知条件和观测数据，设计平差方案，编制平差程序，给出平差结果，评定精度。 题目：A、B为已知点，已知点坐标已给出如下。      为待定点。同精度观测了24个角度      ，测角中误差        。测量了17条边长      ，其观测结果及中误差列于表4中，试按间接平差法求： （1）平差后单位权中误差； （2） 14个待定点的的坐标平差值及中误差； （3）平差后  边的相对中误差。 已知点坐标： =730.024 ， =126.040 ， ， 表一：导线网角度观测值 编号 角度观测值L 编号 角度观测值L 编号 角度观测值L 1 87 33 15.3 9 220 11 01.2 17 202 00 07.5 2 51 37 34.2 10 197 31 07.6 18 35 55 52.0 3 220 49 01.1 11 223 39 58.3 19 211 14 24.1 4 197 35 58.9 12 98 07 30.0 20 112 49 30.8 5 154 10 39.2 13 218 12 13.1 21 148 24 03.7 6 89 18 20.3 14 43 40 02.3 22 163 09 32.5 7 41 30 40.1 15 199 11 38.9 23 170 33 56.2 8 229 10 48.3 16 158 23 46.5 24 138 40 51.0             表二：导线网边长观测值及精度 编号 边观测值S（m） 中误差（mm） 编号 边长观测值S（m） 中误差（mm） 编号 边长观测值S（m） 中误差（mm） 1 169.435 13 7 144.605 12 13 110.02 10.5 2 195.782 14 8 136.578 11.7 14 119.1 10.9 3 184.984 13.6 9 112.502 10.6 15 123.94 11.1 4 110.036 10.5 10 131.074 11.5 16 122.84 11.1 5 163.147 12.8 11 140.88 11.9 17 139.46 11.8 6 219.482 14.9 12 132.608 11.5                   图一： 二、程序编写思路 1、数学模型 ◆函数模型  观测方程    误差方程  ◆随机模型  ◆方程的解与精度评定 令       法方程即为 可得解为：     故可得平差结果： 从而解得单位权中误差为    2、程序解析 ◆计算框图  三、程序运行结果 ◆平差后单位权中误差： =  5.30 (由于程序取的先验单位权中误差为10，两者相比较，故可以认为这个结果是可以接受的) ◆14个待定点的坐标平差值及中误差见下表： 表三： 点号 坐标平差值（m） 中误差（mm） 点号 坐标平差值（m） 中误差（mm） P1 X 678.1642 3.3638 P8 X 767.8779 4.8043 Y 287.3411 5.8339 Y 401.851 6.3376 P2 X 564.6921 5.9695 P9 X 745.7234 6.8522 Y 446.8827 7.4627 Y 531.041 6.5826 P3 X 533.8472 8.7390 P10 X 671.6262 9.0076 Y 629.2749 7.3425 Y 650.8636 6.0122 P4 X 475.6779 7.3449 P11 X 794.2843 8.4415 Y 535.8728 8.4410 Y 600.4687 6.1380 P5 X 499.141 5.7683 P12 X 859.0541 6.6941 Y 374.2314 7.6319 Y 511.5346 6.6419 P6 X 594.6191 5.1952 P13 X 898.2695 4.7257 Y 176.7976 3.3838 Y 399.0808 6.3770 P7 X 826.3848 2.7253 P14 X 919.3467 3.1685 Y 305.7569 5.4103 Y 276.9424 4.7680                 表七： 点号 中误差（mm） 点号 中误差（mm） 点号 中误差（mm） 1 6.7342 6 6.2000 11 10.4371 2 9.5565 7 6.0579 12 9.4300 3 11.4141 8 7.9527 13 7.9371 4 11.1892 9 9.5017 14 5.7248 5 9.5665 10 10.8297                 ◆平差后 边的相对中误差为： 四、程序源代码 这是要输入的文件 ◆%%从文件中输入数据 fid1=fopen('E:\学习\平差\S3.txt','rt'); [S1,count]=fscanf(fid1,'%f ',[7 24]); %#ok fclose(fid1); S1=S1'; fid2=fopen('E:\ \学习\平差\S4.txt','rt'); [S2,count]=fscanf(fid2,'%f ',[5 17]); fclose(fid2); S2=S2'; [S3]=convert(S1); [X0]=daoxian(S2,S3); [B,l,P]=BBll(X0,S2,S1,S3); x=inv(B'*P*B)*B'*P*l; X=X0'+x/1000; disp(X); ◆%%计算B，l，P 矩阵 function [B,l,P]=BBll(X0,S2,S1,S3) %#ok X10=[ 730.024 126.040 855.111  172.232 ]; for i=1:1:7 if i<2 S=sqrt((X0(1)-X10(1))^2+(X0(2)-X10(2))^2); B(1,1)=(X0(1)-X10(1))/S; B(1,2)=(X0(2)-X10(2))/S; l(1,1)=1000*(S2(1,2)-S); elseif i>6 S=sqrt((X10(1)-X0(11))^2+(X10(2)-X0(12))^2); B(i,11)=-(X10(1)-X0(11))/S; %#ok B(i,12)=-(X10(2)-X0(12))/S; %#ok l(i,1)=1000*(S2(7,2)-S); %#ok elseif i>1&&i<7 S=sqrt((X0(2*i-1)-X0(2*(i-1)-1))^2+(X0(2*i)-X0(2*(i-1)))^2); B(i,2*(i-1)-1)=-(X0(2*i-1)-X0(2*(i-1)-1))/S; %#ok B(i,2*(i-1))=-(X0(2*i)-X0(2*(i-1)))/S; %#ok B(i,2*i-1)=(X0(2*i-1)-X0(2*(i-1)-1))/S; %#ok B(i,2*i)=(X0(2*i)-X0(2*(i-1)))/S; %#ok l(i,1)=1000*(S2(i,2)-S); %#ok end end for i=8:1:16 if i<9 S=sqrt((X0(13)-X10(3))^2+(X0(14)-X10(4))^2); B(8,13)=(X0(13)-X10(3))/S; B(8,14)=(X0(14)-X10(4))/S; l(8,1)=1000*(S2(8,2)-S); elseif i>15 S=sqrt((X10(3)-X0(27))^2+(X10(4)-X0(28))^2); B(16,27)=-(X10(3)-X0(27))/S; B(16,28)=-(X10(4)-X0(28))/S; l(16,1)=1000*(S2(16,2)-S); elseif i<16&&i>8 S=sqrt((X0(2*(i-1)-1)-X0(2*(i-1)-3))^2+(X0(2*(i-1))-X0(2*(i-1)-2))^2); B(i,2*(i-1)-3)=-(X0(2*(i-1)-1)-X0(2*(i-1)-3))/S; B(i,2*(i-1)-1)=(X0(2*(i-1)-1)-X0(2*(i-1)-3))/S; B(i,2*(i-1)-2)=-(X0(2*(i-1))-X0(2*(i-1)-2))/S; B(i,2*(i-1))=(X0(2*(i-1))-X0(2*(i-1)-2))/S; l(i,1)=1000*((S2(i,2)-S)); end end S=sqrt((X0(19)-X0(5))^2+(X0(20)-X0(6))^2); B(17,5)=-(X0(19)-X0(5))/S; B(17,19)=(X0(19)-X0(5))/S; B(17,6)=-(X0(20)-X0(6))/S; B(17,20)=(X0(20)-X0(6))/S; l(17,1)=1000*(S2(17,2)-S); row=206.265; for i=1:1:24 [M1]=fangweijiao(S1(i,6),S1(i,5),X0); [M2]=fangweijiao(S1(i,6),S1(i,7),X0); if S1(i,7)>0 B(17+i,2*S1(i,7)-1)=-row*M2(1); B(17+i,2*S1(i,7))=row*M2(2);