资源描述：
《南京市明城墙监测回归分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
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=polyval(p,xdat)求polyfit所得的回归多项式在xdat处的预测值Y；(2)[Y
2、hat,DELTA]=polyconf(p,xdat,s,alpha)求polyfit所得的回归多项式在x处的预测值Yhat及预测值的置信水平为1-alpha的置信区间YhatDELTA；alpha缺省时为0.05.3.多项式拟合交互式工具polytoolh=polytool(X,Y,m,alpha,xname,yname)线对象句柄值多项式次数显著性水平X轴标签Y轴标签南京市古城墙监测数据分析项目简介：南京市明城墙变形监测项目，是由东南大学土木学院和测量工程系共同承接，该项目自2008年启动至今，已有
3、近三年时间。而且该项目作为南京市的重点民生工程之一，将由南京市政府牵头，由东南大学承接，长期开展下去。明城墙监测数据城墙监测中采用了多种监测手段，包括江苏CORS系统，传统的全站仪观测以及自动化程度较高的裂缝仪。这里我选取的是南京市西水关到东水关的23个观测点的CORS观测数据（X，Y坐标）作为分析数据。1—10号点前6期观测数据的X坐标值1—10号点前6期观测数据的Y坐标值对XD1号点的数据分析以观测期数为横坐标，分别以XD1号点的各期观测X，Y坐标为纵坐标，利用MATLAB做回归分析。MATLAB程
4、序代码：functionX_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');gridon;xlabel(‘观测期数’);ylabel(‘X坐标值’)；title(‘1号点X坐标多项式回归分析');set(gca,'color','w');[p,s]=polyfit(x
5、,y,2);[yhat,delta]=polyconf(p,x,s);[yyhaty-yhatyhat-deltayhat+delta];pfigure(2)plot(x,y,'bo',x,yhat,'r',x,yhat-delta,'b',x,yhat+delta,'b');gridon;xlabel(‘观测期数’);ylabel(‘X坐标值’)；title(‘1号点X坐标多项式回归分析');set(gca,'color','w');h=polytool(x,y,2,0.05,'观测期数','X坐标值
6、');散点分布二次多项式拟合直接用交互式工具polytool进行拟合P=0.00035178571332-0.001765357141272.35771000001114拟合函数：当alpha取0.05时，各期的置信区间为：置信区间应该再加上493470同理，对于Y坐标的回归分析结果如下：散点分布：二次多项式拟合直接用交互式工具polytool进行拟合P=-0.000642857123500.003757142706314.04760000021197拟合函数：当alpha取0.05时，各期的置信区间为
7、：置信区间应该再加上3545010谢谢！！
