应用MATLAB进行非线性回归分析

应用MATLAB进行非线性回归分析摘要早在十九世纪，英国生物学家兼统计学家高尔顿在研究父与子身高的遗传问题时，发现子代的平均高度又向中心回归大的意思，使得一段时间内人的身高相对稳定。之后回归分析的思想渗透到了数理统计的其他分支中。随着计算机的开展，各种统计软件包的出现，回归分析的应用就越来越广泛。回归分析处理的是变量与变量间的关系。有时，回归函数不是自变量的线性函数，但通过变换可以将之化为线性函数，从而利用一元线性回归对其进行分析，这样的问题是非线性回归问题。下面的第一题：炼钢厂出钢水时用的钢包，在使用过程中由于钢水及炉渣对耐火材料的侵蚀，使其容积不断增大。要找出钢包的容积用盛满钢水时的质量与相应的实验次数的定量关系表达式，就要用到一元非线性回归分析方法。首先我们要对数据进行分析，描出数据的散点图，判断两个变量之间可能的函数关系，对题中的非线性函数，参数估计是最常用的“线性化方法〞，即通过某种变换，将方程化为一元线性方程的形式，接着我们就要对得到的一些曲线回归方程进行选择，找出到底哪一个才是更好一点的。此时我们通常可采用两个指标进行选择，第一个是决定系数，第二个是剩余标准差。进而就得到了我们想要的定量关系表达式。第二题：给出了某地区1971—2000年的人口数据，对该地区的人口变化进行曲线拟合。也用到了一元非线性回归的方法。首先我们也要对数据进行分析，描出数据的散点图，然后用MATLAB编程进行回归分析拟合计算输出利用

Logistic模型拟合曲线。关键词：参数估计，

Logistic模型，MATLAB正文一、一元非线性回归分析的求解思路：求解函数类型并检验。求解未知参数。可化曲线回归为直线回归，用最小二乘法求解；可化曲线回归为多项式回归。二、回归曲线函数类型的选取和检验1、直接判断法2、作图观察法，与典型曲线比拟，确定其属于何种类型，然后检验。3、直接检验法〔适应于待求参数不多的情况〕4、表差法〔适应于多想式回归，含有常数项多于两个的情况〕三、化曲线回归为直线回归问题用直线检验法或表差法检验的曲线回归方程都可以通过变量代换转化为直线回归方程，利用线性回归分析方法可求得相应的参数估计值。题目：例炼钢厂出钢水时用的钢包，在使用过程中由于钢水及炉渣对耐火材料的浸蚀，其容积不断增大。现在钢包的容积用盛满钢水时的重量y(kg)表示，相应的试验次数用x表示。数据见表，要找出y与x的定量关系表达式。表8.5.1钢包的重量y与试验次数x数据序号xy序号xy12106.42811110.5923108.20914110.6034109.581015110.9045109.501116110.7657110.001218111.0068109.931319111.20710110.49

1)1/y=a+b/xy=a+blnxformatlongx=[23457810111415161819];y=[106.42108.20109.58109.5110109.93110.49110.59110.60110.9110.76111111.20];plot(x,y,‘k+’);%数据的散点图x1=1./x;y1=1./y;plot(x1,y1,‘k+’);%变换后数据的散点图x2=[ones(13,1)x1'];[b,bint,rint,stats]=regress(y1',x2);z=b(1)+b(2)*x1;yc=1./z;plot(x1,y1,‘k+’,x1,z,‘r’)%变换后数据的散点图和回归直线图变换后数据的散点图及回归直线图R2=1-sum((y-yc).^2)/lyy;%模型的拟合优度系数

plot(x,y,‘k+’,x,yc,‘r’)%数据的散点图和回归曲线图

legend('散点图','回归函数')b=0.00896662968057

0.00082917436336

R2=0.97292374957556第一种方法的程序：formatlongx=[23457810111415161819];y=[106.42108.20109.58109.5110109.93110.49110.59110.60110.9110.76111111.20];plot(x,y,‘k+’);%数据的散点图x1=1./x;y1=1./y;plot(x1,y1,‘k+’);%变换后数据的散点图x2=[ones(13,1)x1'];[b,bint,rint,stats]=regress(y1',x2);z=b(1)+b(2)*x1;yc=1./z;plot(x1,y1,‘k+’,x1,z,‘r’)%变换后数据的散点图和回归直线图n=length(x);lyy=sum(y.^2)-n*(mean(y))^2;R2=1-sum((y-yc).^2)/lyy;%模型的拟合优度系数b=0.00896662968057

0.00082917436336

R2=0.97292374957556用类似的方法可以得出其它三个曲线回归方程，它们分别是：第二种方法的程序：formatlongx=[23457810111415161819];y=[106.42108.20109.58109.5110109.93110.49110.59110.60110.9110.76111111.20];x1=log(x);y1=y;x2=[ones(13,1)x1'];[b,bint,rint,stats]=regress(y1',x2);bz=b(1)+b(2)*x1;yc=z;n=length(x);lyy=sum(y.^2)-n*(mean(y))^2;R2=1-sum((y-yc).^2)/lyy;plot(x,y,'k+',x,yc,‘c');legend('散点图','回归函数')b=1.0e+002*R2=0.87731500489620第三种方法的程序：formatlongx=[23457810111415161819];y=[106.42108.20109.58109.5110109.93110.49110.59110.60110.9110.76111111.20];x1=sqrt(x);y1=y;x2=[ones(13,1)x1'];[b,bint,rint,stats]=regress(y1',x2);bz=b(1)+b(2)*x1;yc=z;n=length(x);lyy=sum(y.^2)-n*(mean(y))^2;R2=1-sum((y-yc).^2)/lyy;plot(x,y,'k+',x,yc,'k');legend('散点图','回归函数')b=1.0e+002*1.06301275014382三种方法的拟合效果比拟：R2=0.97292374957556R2=0.877315004896201.原始数据下表给出了某地区1971—2000年的人口数据〔表1〕。试分别用Matlab和SPSS软件，对该地区的人口变化进行曲线拟合。表1某地区人口变化数据年份时间变量t=年份-1970人口y/人197113381519722339811973334004197443416519755342121976634327197773434419788344581979934498198010344761981113448319821234488198313345131984143449719851534511198616345201987173450719881834509198919345211990203451319912134515199222345171993233451919942434519199525345211996263452119972734523199828345251999293452520003034527根据上表中的数据，做出散点图，见图1。图1某地区人口随时间变化的散点图从图1可以看出，人口随时间的变化呈非线性过程，而且存在一个与横坐标轴平行的渐近线，故可以用Logistic曲线模型进行拟合。因为Logistic曲线模型的根本形式为：所以，只要令：，就可以将其转化为直线模型：下面，我们分别用Matlab和SPSS软件进行回归分析拟合计算。2．用Matlab编程进行回归分析拟合计算源程序〔Nonlinear-Regression-Model.m〕，如下：clearclc%读入人口数据〔1971－2000年〕y=[338153398134004341653421234327343443445834498344763448334488345133449734511345203450734509345213451334515345173451934519345213452134523345253452534527];%读入时间变量数据〔t＝年份－1970〕t=[123456789101112131415161718192021222324252627282930];%线性化处理fort=1:30,x(t)=exp(-t);y(t)=1/y(t,1);end%计算，并输出回归系数Bc=zeros(30,1)+1;X=[c,x'];B=inv(X'*X)*X'*yfori=1:30,%计算回归拟合值z(i)=B(1,1)+B(2,1)*x(i);%计算离差s(i)=y(i)-sum(y)/30;%计算误差w(i)=z(i)-y(i);end%计算离差平方和SS=s*s';%回归误差平方和QQ=w*w';%计算回归平方和UU=S-Q;%计算，并输出F检验值F=28*U/Q%计算非线性回归模型的拟合值forj=1:30,Y(j)=1/(B(1,1)+B(2,1)*exp(-j));end%输出非线性回归模型的拟合曲线〔Logisic曲线〕plot(T,Y)上述程序运行后，输出〔1〕输出回归系数B及F检验值如下：B=1.0e-0