基于南京市明城墙监测 数据的回归模型分析目录 回归模型函数介绍 南京市明城墙监测数据回归模型分析多项式回归函数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. 多项式拟合交互式工具polytoolh = polytool ( X, Y, m, alpha, xname, yname )线对象句柄值多项式次数显著性水平X轴标签Y轴标签南京市古城墙监测数据分析项目简介：南京市明城 墙变形监测项目，是由 东南大学土木学院和测 量工程系共同承接，该 项目自2008年启动至今 ，已有近三年时间。而 且该项目作为南京市的 重点民生工程之一，将 由南京市政府牵头，由 东南大学承接，长期开 展下去。明城墙监测数据城墙监测中采用了多种监测手段，包括江苏CORS系统，传统的全站仪观测 以及自动化程度较高的裂缝仪。这里我选取的是南京市西水关到东水关的23 个观测点的CORS观测数据（X，Y坐标）作为分析数据。110号点前6期观测数据的X坐标值110号点前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.000642857123500.003757142706314.04760000021197拟合函数：当alpha取0.05时，各期的置信区间为：置信区间应该再加上3545010谢 谢 ！