测量数据处理实验一回归和拟合.docVIP

《测量数据处理》实验报告 （测绘工程专业使用） 班级： 测绘112班 学号： 2011122712 姓名： 许艳超 淮海工学院测绘工程学院 实验一 回归和拟合 一、实验目的和要求 掌握线性回归模型的建立、解算和回归假设检验； 初步提高利用Matlab进行编制程序解决测量计算问题的能力； 能使用回归方法处理测量数据。 二、实验时间、地点 测量数据处理实验室 三、实验内容: 1.在某地区，一个时期内，对20个测点进行了重复重力测量，其结果列于下表，其中（）为第i测点的坐标值，为第i测点上的重力变化值。 测点的X，Y坐标及重力变化值 序号 1 2 3 4 5 6 7 8 9 10 x 9.92 10.46 12.45 13.29 15.00 14.98 14.32 12.13 15.74 16.10 y 7.17 5.13 3.91 2.84 2.87 5.52 7.64 7.31 8.02 9.72 2.20 1.24 1.50 0.44 0 0.90 1.60 1.64 1.24 1.18 序号 11 12 13 14 15 16 17 18 19 20 x 15.25 14.76 13.77 11.89 12.09 9.84 7.94 7.55 8.90 8.57 y 14.54 11.77 12.29 11.62 9.57 10.14 10.80 8.90 8.08 7.15 1.56 1.35 1.59 1.29 1.87 1.13 2.41 2.09 1.08 2.70 1）试求该区域重力变化一次趋势面： 和二次趋势面： 2）讨论本地区取哪种趋势面效果好？说出你选择的原因。 答： （1）求重力变化一次趋势面： 根据题目要求，我利用matlab的解算步骤如下： 输入已知数据X、Y的值 X=[1,9.92,7.17;1,10.46,5.13;1,12.45,3.91;1,13.29,2.84;1,15,2.87;1,14.98,5.52;1,14.32,7.64;1,12.13,7.31;1,15.74,8.02;1,16.1,9.72;1,15.25,14.54;1,14.76,11.77;1,13.77,12.29;1,11.89,11.62;1,12.09,9.57;1,9.84,10.14;1,7.94,10.8;1,7.55,8.9;1,8.9,8.08;1,8.57,7.15]; Y=[2.2;1.24;1.5;0.44;0;0.9;1.6;1.64;1.24;1.18;1.56;1.35;1.59;1.29;1.87;1.13;2.41;2.09;1.08;2.7]; 然后进行解算 软件中输入代码如下 N=inv(X*X)； A=N*X*Y；YY=X*A； Q1=Y-YY； q1=0； for i=1:20 q1=q1+Q1(i,1)*Q1(i,1) ; end q1 ym=0； for i=1:20 ym=ym+Y(i,1) End ym=ym/20 q2=0 for i=1:20 q2=q2+(Y(i,1)-ym)*(Y(i,1)-ym)； end q2 q3=0； for i=1:20 q3=q3+(YY(i,1)-ym)*(YY(i,1)-ym)； end q3 m=sqrt(q1/17) F=(q3/2)/(q1/17) T=[1;1;1]； for i=2:3 T(i,1)=A(i,1)/(m*sqrt(N(i,i))) end 求二次趋势面： 根据题目要求，我利用matlab的解算步骤如下： 先输入X的值： X=[1,9.92,7.17;1,10.46,5.13;1,12.45,3.91;1,13.29,2.84;1,15,2.87;1,14.98,5.52;1,14.32,7.64;1,12.13,7.31;1,15.74,8.02;1,16.1,9.72;1,15.25,14.54;1,14.76,11.77;1,13.77,12.29;1,11.89,11.62;1,12.09,9.57;1,9.84,10.14;1,7.94,10.8;1,7.55,8.9;1,8.9,8.08;1,8.57,7.15]; 再重新组成矩阵X，并输入输入Y的值 最后利用如下代码进行解算 N=inv(X*X) ； A=N*X*Y；YY=X*A； Q1=Y-YY ； q1=0； for