多项式回归、非线性回归模型

1、多项式回归、非线性回归模型关键词：回归方程的统计检验、拟合优度检验、回归方程的显著性检验、F检验、回归系数的显著性检验、残差分析、一元多项式回归模型、一元非线性回归模型一、回归方程的统计检验1 .拟合优度检验1 .概念介绍SST总离差平方和totalSSR回归平方和regressionSSEtiJ余平方和errornn.二（yi-?i）.二（？i-yi）2R2=1_i=1_i=1n-n.二(yi-yi).二(yi-yiVidi12.例题1存在四点（-2,-3）、2.回归方程的显著性检验F例6（F检验）在合金钢强度的例1中，我们已求出了回归方程，这里考虑关于回归方程的显著性检验,经计算有：(-1

2、,-1)、nzimn（1,2）、（4,3）求拟合直线与决定系数。(4-Yi)2二(yi-yi)2/(n-2)i1SSASSE/(n-2)表5X射线照射次数与残留细菌数的方差分析表来源平方和自由度均方F比p值回归SR=327.34fR=1MSR=327.34184.940.0000残差Se=17.72fe=10MSe=1.77总计ST=345.06fT=11这里值很小，因此，在显著性水平0.01下回归方程是显著的。3 .回归系数的显著性检验4 .残差分析二、一元多项式回归模型模型如以下形式的称为一元多项式回归模型:nn1y/xanx-aiXao例1(多项式回归模型)为了分析X射线的杀菌作用，用2

3、00千伏的X射线来照射细菌，每次照射6分钟，用平板计数法估计尚存活的细菌数。照射次数记为t,照射后的细菌数为y见表1。试求：(1)给出y与t的二次回归模型。(2)在同一坐标系内作出原始数据与拟合结果的散点图。(3)预测t=16时残留的细菌数。(4)根据问题的实际意义，你认为选择多项式函数是否合适？表1X射线照射次数与残留细菌数t123456789101112131415y3522111971601421061046056383632211915程序1t=1:15;y=3522111971601421061046056383632211915;p=polyfit(t,y,2)%作二次多项式回归y

4、1=polyval(p,t);%模型估计与作图plot(t,y,'-*',t,y1,'-o');%在同一坐标系中做出两个图形legend('原始数据',二次函数')xlabel('t(照射次数)')%横坐标名ylabel('y(残留细菌数)')纵坐标名t0=16;yc1=polyconf(p,t0)%预测t0=16时残留的细菌数，方法1yc2=polyval(p,t0)%预测t0=16时残留的细菌数，方法2即二次回归模型为：y1=1.9897t2-51.1394t347.8967O40oooooO50505

5、0332211数菌细留声y原始数据二次函数+4-51015t（照射次数）图1原始数据与拟合效果的散点图原始数据与拟合结果的散点图如图所示，从图形可知拟合效果较好。照射16次后，用二次函数计算出细菌残留数为39.0396,显然与实际不符。由实际问题的意义可知，尽管二次多项式拟合效果较好，但是用于预测并不理想。因此如何根据原始数据散点图的规律，选择适当的回归曲线是非常重要的，这样就有必要给出非线性回归模型。三、一元非线性回归模型为了便于正确选择合适的函数进行回归分析建模，我们给出通常选择的6类曲线：1b（1）双曲线1-a-（如图所不）yx（2）哥函数曲线y=axb,其中x>0,a>0（

6、如图所示）（3）指数曲线y=aebx,其中参数a>0（如图所示）（4）倒指数曲线y=aeb/x,其中a>0（如图所示）（5）对数曲线y=a+blnx（如图所示）一，,1,._,（6）S型曲线y=一，其中ab0（如图所不）abe非线性回归建模通常有两种方法：一是通过适当的变换转化为线性回归模型，例如双曲1b11.线模型一=a十一（如图1所不）,如果作变换y=,x=一则有y=a+bx,此时就yxyx是线性回归模型。如果无法实现线性化，可以利用最小二乘法直接建立非线性回归模型，求解最佳参数。例2（非线性回归模型、置信区间）炼钢厂出钢时所用盛钢水的钢包，由于钢水对耐火材料的侵蚀，容积不断增

7、大，我们希望找出使用次数与增大容积之间的函数关系。实验数据见表2。1b(1)建立非线性回归模型1=ab;yx(2)预测钢包使用xo=17次后增大的容积y0；(3)计算回归模型参数的置信度为95%的置信区间。表2钢包使用次数与增大容积使用次数(x)2345678910111213141516增大容枳6.428.29.589.59.7109.939.9910.4910.5910.610.810.610.910.76(y)解：(1)建立非线性回归模型：程序2x=2:16;y=6.428.29.589.59.7109.939.9910.4910.5910.610.810.610.910.76;%建立非

8、线性双曲线回归模型b0=0.084,0.1436;%回归系数初值fun=inline('x./(b(1)*x+b(2)','b','x');%建立函数beta,r,J=nlinfit(x,y,fun,b0);%非线性拟合命令；其中，beta表示最佳回归系数的估计值，r是残差，J是雅可比矩阵beta%输出最佳参数y1=x./(0.0845*x+0.1152);%拟合曲线plot(x,y,'*',x,y1,'-or')legend('原始数据',拟合曲线')legend为图例命令初始值要先计算后

9、才能得到上面程序中的b0,选择已知程序中的点(2,6.42)和点(16,10.76),可选择手工方法解方程，也可利用以下MATLAB程序求解。程序3a,b=solve('1/6.42=a+b/2','1/10.76=a+b/16')%解方程注：当所求解的方程过于复杂时，MATLAB运行会出现错误，此时需将方程尽量化简后再进行求解，如以下形式：a,b=solve('6.42*(2*a+b)=2','10.76*(16*a+b)=16')运行程序3可得到最佳参数为a=0.0845、b=0.1152,求解得到钢包使用次数与增大容积的非线

10、性拟合图，如图2所示。11,r,."原始数据I一105T*拟合曲线-10士*r9.5T*.9.8.5.8.-7.5,一7,-6.5,一6246810121416图2钢包使用次数与增大容积的非线性拟合图(2)预测钢包使用17次后增大的容积：程序4ypred=nlpredci(fun,17,beta,r,J)%预测钢包使用17次后增大的容积(3)置信区间：程序5ci=nlparci(beta,r,J)%置信区间运行后得到ci=0.08140.08760.09340.1370即回归模型中参数的置信度为的置信区间分别为0.0814,0.0876与0.0934,0.1370。我们求出的最佳参数分别为a=0.0845和b=0.1152,均属于上述置信区间