《南京市明城墙监测回归分析》由会员分享，可在线阅读，更多相关《南京市明城墙监测回归分析（21页珍藏版）》请在金锄头文库上搜索。
1、基于南京市明城墙监测数据的回归模型分析目录 回归模型函数介绍回归模型函数介绍 南京市明城墙监测数据回归模型分析南京市明城墙监测数据回归模型分析多项式回归函数多项式回归函数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=poly
2、val ( 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 )线线对对象象句句柄
3、柄值值多多项项式式次次数数显显著著性性水水平平X轴标签轴标签Y轴标签轴标签南京市古城墙监测数据分析南京市古城墙监测数据分析项目简介：项目简介：南京市明城墙变形监测项目，是由东南大学土木学院和测量工程系共同承接，该项目自2008年启动至今，已有近三年时间。而且该项目作为南京市的重点民生工程之一，将由南京市政府牵头，由东南大学承接，长期开展下去。明城墙监测数据明城墙监测数据城墙监测中采用了多种监测手段，包括江苏城墙监测中采用了多种监测手段，包括江苏CORS系统，传统的全站仪观测系统，传统的全站仪观测以及自动化程度较高的裂缝仪。这里我选取的是南京市西水关到东水关的以及自动化程度较高的裂缝仪。这里我选
4、取的是南京市西水关到东水关的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，
5、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;pfigure(2)plot(x,y,bo,x,yhat,r,x,yhat-delta,b,x,yhat+delta,b);grid on;xlabel(观测期数); ylabel(X坐标值)
6、；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谢谢 谢谢 ！
