南京市明城墙监测回归分析.pptVIP

基于南京市明城墙监测数据的回归模型分析 目录 ? 回归模型函数介绍 ? 南京市明城墙监测数据回归模型分析 多项式回归函数polyconf , polyfit, (polytool) 1. 确定回归系数的点估计值 [ p, s ] = polyfit ( X, Y, m ) 次数 2. 预测和预测误差估计 [ Yhat ] = polyval ( p, xdat ) 预测值 系数估计值 指定x 值 [Yhat, DELTA] = polyconf ( p, xdat, s, alpha ) (1) Yhat=polyval ( p, xdat) 求polyfit所得的回归多项式在xdat 处的预测值Y； (2) [Yhat, DELTA]=polyconf ( p, xdat, s, alpha) 求polyfit所得的回归多项式在x 处的预测值Yhat及预测值的置信水平为1- alpha的置信区间Yhat DELTA；alpha缺省时为0.05. 3. 多项式拟合交互式工具polytool h = polytool ( X, Y, m, alpha, xname, yname ) 线对象句柄值 多项式次数 显著性水平 X轴标签 Y轴标签 南京市古城墙监测数据分析 项目简介：南京市明城墙变形监测项目，是由东南大学土木学院和测量工程系共同承接，该项目自2008年启动至今，已有近三年时间。而且该项目作为南京市的重点民生工程之一，将由南京市政府牵头，由东南大学承接，长期开展下去。 明城墙监测数据 城墙监测中采用了多种监测手段，包括江苏CORS系统，传统的全站仪观测以及自动化程度较高的裂缝仪。这里我选取的是南京市西水关到东水关的23个观测点的CORS观测数据（X，Y坐标）作为分析数据。 1—10号点前6期观测数据的X坐标值 1—10号点前6期观测数据的Y坐标值 对XD1号点的数据分析 以观测期数为横坐标，分别以XD1号点的各期观测X，Y坐标为纵坐标，利用MATLAB做回归分析。 MATLAB程序代码： function X_reg_analysis() x=[1,2,3,4,5,6]; y0=[493472.357，493472.3552，493472.3535，493472.3578，493472.3589，493472.3588]; y=(y0-493470); figure(1) plot(x,y,bo); grid on; xlabel(‘观测期数’); ylabel(‘X坐标值’)； title(‘1号点X坐标多项式回归分析); set(gca,color,w); [p,s]=polyfit(x,y,2); [yhat,delta]=polyconf(p,x,s); [y yhat y-yhat yhat-delta yhat+delta]; p figure(2) plot(x,y,bo,x,yhat,r,x,yhat-delta,b,x,yhat+delta,b); grid on; xlabel(‘观测期数’); ylabel(‘X坐标值’)； title(‘1号点X坐标多项式回归分析); set(gca,color,w); h=polytool(x,y,2,0.05, 观测期数, X坐标值); 散点分布 二次多项式拟合 直接用交互式工具polytool进行拟合 P=0.00035178571332 -0.00176535714127 2.35771000001114 拟合函数： 当alpha取0.05时，各期的置信区间为： 置信区间应该再加上493470 同理，对于Y坐标的回归分析结果如下： 散点分布： 二次多项式拟合 直接用交互式工具polytool进行拟合 P= -0.00064285712350 0.00375714270631 4.04760000021197 拟合函数： 当alpha取0.05时，各期的置信区间为： 置信区间应该再加上3545010 谢 谢 ！！