《测量数据处理实验报告》由会员分享，可在线阅读，更多相关《测量数据处理实验报告（27页珍藏版）》请在金锄头文库上搜索。
1、ftl院系目 录（页码根据自己需要填写）一、实验项目须知二、实验一 秩亏自由网上机实验三、实验二 极大验后滤波、推估上机实验四、实验三 最小二乘配置上机实验五、实验四 平差系统的统计假设检验上机实验六、编写实验报告要求实验项目须知一、课程介绍本课程是测绘工程类专业的专业选修课程之一，是在学生学习了误差理论与测量平差知 识之后又开设的一门专业的数据处理课程，该课程主要讲授近现代测量数据处理理论、模型 和方法，同时紧密结合计算机编程，从而提高学生的计算机应用能力、测量数据处理能力， 培养学生养成严谨的科学态度，提高分析和解决问题的能力，更好的适应今后工作和学习。二、测量数据处理实验教学要求1. 掌
2、握稳健估计方法处理测量数据；2. 掌握极大验后滤波、推估理论和方法处理测量数据；3. 初步掌握最小二乘配置理论和方法处理测量数据；4. 初步掌握工程或科研项目中测量数据处理的理论、步骤、方法和技巧；5. 掌握测量数据处理中常用的计算机语言和软件；6. 服从实验指导教师的规定，认真、按时、独立完成任务每次实验结束，提交书写工整或按 照规范打印实验报告。7. 在实验过程中，还应遵守纪律，爱护实验室设备、离开后将自己的随身物品带走，并记录 仪器使用状况、将凳子整理好并放到电脑桌的下方。8. 上机实验后，提交实验报告，成绩评定重点考核算法、处理结果与结果分析。三、实验选择要求可根据课堂讲课具体内容，完
3、成4 个实验，并要求写出实验报告。实验一 秩亏自由网上机实验实验目的1. 掌握各种秩亏自由网平差的函数模型和随机模型2. 熟悉各种秩亏自由网平差解算的若干公式；3. 结合算例能利用 matlab 进行编程计算。二、 实验内容在下图水准网中，观测高差、距离和各待定点高程近似值见教材P43表2-1,分别进行下列自由网平差：(1) 以 6 号点为固定点的经典自由网平差；(2) 以重心基准的自由网平差(Px=E)；(3) 以 1， 2， 5， 6 四个点为拟稳基准的拟稳平差.参考程序：(1)首先建立误差方程V二Bx-1B=-1 1 0 0 0 ;-1 0 0 1 0;-1 0 0 0 1 ;0 -1
4、1 0 0 ;0 -1 0 1 0 .;0 -1 0 0 0 ;0 0 -1 1 0 ;0 0 -1 0 1 ;0 0 0 0 -1 ; l=-1;4;5;8;16;2;-13;-20;-1;S=204.1 188.7 344.8 149.2 142.9 250 128.2 98 196.1; P=diag(10./S);N=B*P*B; x=inv(N)*B*P*l V=B*x-l Qx=inv(N) 首先建立误差方程V二B X-1B=-1 1 0 0 0 0;-1 0 0 1 0 0;-1 0 0 0 1 0;0 -1 1 0 0 0;0 -1 0 1 0 0 . ;0 -1 0 0 0
5、1;0 0 -1 1 0 0;0 0 -1 0 1 0;0 0 0 0 -1 1;l=-1;4;5;8;16;2;-13;-20;-1S=204.1 188.7 344.8 149.2 142.9 250 128.2 98 196.1; P=diag(10./S);S1=1;1;1;1;1;1;W=B*P*l;Xr=inv(B*P*B+S1*S1)*WV=B*Xz-l Qxr=inv(B*P*B+S1*S1)-inv(B*P*B+S1*S1)*S1*S1*inv(B*P*B+S1*S1)(3)首先建立误差方程V二Bx-1B=-1 1 0 0 0 0;-1 0 0 1 0 0;-1 0 0 0
6、1 0;0 -1 1 0 0 0;0 -1 0 1 0 0 .;0 -1 0 0 0 1;0 0 -1 1 0 0;0 0 -1 0 1 0;0 0 0 0 -1 1; l=-1;4;5;8;16;2;-13;-20;-1S=1;1;1;1;1;1;S1=204.1 188.7 344.8 149.2 142.9 250 128.2 98 196.1;P=diag(10./S1);N=B*P*BW=B*P*l;px=1 1 0 0 1 1Px=diag(px)Ss=Px*SXs=inv(B*P*B+Ss*Ss)*WV=B*Xs-lQxs=inv(B*P*B+Ss*Ss)*N*inv(B*P*B
7、+Ss*Ss)三、 实验心得和体会实验二 极大验后滤波、推估上机实验一、实验目的和要求1. 掌握极大验后滤波与推估的基本原理和计算公式；2. 能根据观测值和给定的数学模型计算滤波信号的协方差矩阵的估计矩阵；3. 能使用 matlab 进行有关的计算；二、实验内容:1.（基本要求题目）沿A、B连线在A、l、2、3、4等五个点测定了大气温度，获得18组观测值 x ,其结果见下表，已知各点的距离 S （j 1、2、3、4）分别为 ij Aj4.511,10.747,16.753,22.220（km）.假定AB连线上的大气温度，是一个以各点至A的距离S为因素的平稳随机过程，试估计此随机过程的协方差函数
8、。观测数据序号A1234116.415.916.315.815.4216.316.116.215.915.5316.816.415.916.316.3416.616.216.316.216.3516.816.716.116.216.3616.616.616.416.216.3716.716.617.016.316.2816.816.516.916.516.9916.916.416.116.616.21016.116.616.217.117.01116.516.716.716.816.51217.116.916.616.816.71316.816.916.817.117.01417.016.91
9、6.716.917.01517.116.916.716.917.01616.216.016.316.516.41716.115.916.216.416.41816.115.916.016.116.51916.316.016.616.517.02016.116.015.916.516.3参考程序：X=16.4 15.9 16.3 15.8 15.4;16.3 16.1 16.2 15.9 15.5;16.8 16.4 15.9 16.3 16.3; .16.6 16.2 16.3 16.2 16.3;16.8 16.7 16.1 16.2 16.3; .16.6 16.6 16.4 16.2 1
10、6.3;16.7 16.6 17.0 16.3 16.2; .16.8 16.5 16.9 16.5 15.9;16.9 16.4 16.1 16.6 16.2; .16.1 16.6 16.2 17.1 17.0;16.5 16.7 16.7 16.8 16.5; .17.1 16.9 16.6 16.8 16.7;16.8 16.9 16.8 17.1 17.0; .17.0 16.9 16.7 16.9 17.0;17.1 16.9 16.7 16.9 17.0; .16.2 16.0 16.3 16.5 16.4;16.1 15.9 16.2 16.4 16.4; .16.1 15.9
11、 16.0 16.1 16.5;16.3 16.0 16.5 16.6 17.0;16.1 16.0 15.9 16.5 16.3 Dx1=cov(X)s1=0;4.510;10.746;16.754;22.230A=ones(5,1),s1,s1.A2,s1.A3,s1.A4s2=0;6.236;12.242;17.709B=ones(4,1),s2,s2.A2,s2.A3,s2.A4 s3=0;6.006;11.473C=ones(3,1),s3,s3.A2,s3.A3,s3.A4 s4=0;5.467D=ones(2,1),s4,s4.A2,s4.A3,s4.A4E=1 0 0 0 0Y
12、=A;B;C;D;E L=0.1168;0.0958;0.0506;0.0466;0.0394;0.1362;0.0705;0.1023;0.1047;0.1073;0.0567;0.0346; .0.1529;0.1561;0.2203 b=inv(Y*Y)*Y*Lb=b(l),b(2)*10,b(3)*100,b(4)*1000,b(5)*10000%将 km 单位化为 10km 单位S=0 4.511 10.747 16.753 22.220;4.511 0 6.236 12.242 17.709;10.747 6.236 0 6.006 11.473; . 16.753 12.242
13、6.006 0 5.467;22.220 17.709 11.473 5.467 0b=fliplr(b)for i=1:25 S(i)=polyval(b,S(i)/10);endS %即为滤波信号协方差阵的估计矩阵Warning: Unable to create personal MATLAB work folder:C:yonghu?MATLAB Warning: Userpath must be an absolute path and must exist on disk. X=16.4 15.9 16.3 15.8 15.4;16.3 16.1 16.2 15.9 15.5;16
14、.8 16.4 15.9 16.3 16.3; .16.6 16.216.316.216.3;16.816.716.116.216.3;.16.6 16.616.416.216.3;16.716.617.016.316.2;.16.8 16.516.916.515.9;16.916.416.116.616.2;.16.1 16.616.217.117.0;16.516.716.716.816.5;.17.116.916.616.816.7;16.816.916.817.117.0; .17.016.916.716.917.0;17.116.916.716.917.0; .16.216.016.316.516.4;16.115.916.216.416.4; .16.115.916.016.116.5;16.316.016.516.617.0;16.1 1
