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

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

ID:47603468

大小:593.47 KB

页数:13页

时间:2019-09-27

应用MATLAB进行非线性回归分析报告_第1页
应用MATLAB进行非线性回归分析报告_第2页
应用MATLAB进行非线性回归分析报告_第3页
应用MATLAB进行非线性回归分析报告_第4页
应用MATLAB进行非线性回归分析报告_第5页
资源描述:

《应用MATLAB进行非线性回归分析报告》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、实用标准文档应用MATLAB进行非线性回归分析摘要早在十九世纪,英国生物学家兼统计学家高尔顿在研究父与子身高的遗传问题时,发现子代的平均高度又向中心回归大的意思,使得一段时间内人的身高相对稳定。之后回归分析的思想渗透到了数理统计的其他分支中。随着计算机的发展,各种统计软件包的出现,回归分析的应用就越来越广泛。回归分析处理的是变量与变量间的关系。有时,回归函数不是自变量的线性函数,但通过变换可以将之化为线性函数,从而利用一元线性回归对其进行分析,这样的问题是非线性回归问题。下面的第一题:炼钢厂出钢水时用的钢包,在使用过程中

2、由于钢水及炉渣对耐火材料的侵蚀,使其容积不断增大。要找出钢包的容积用盛满钢水时的质量与相应的实验次数的定量关系表达式,就要用到一元非线性回归分析方法。首先我们要对数据进行分析,描出数据的散点图,判断两个变量之间可能的函数关系,对题中的非线性函数,参数估计是最常用的“线性化方法”,即通过某种变换,将方程化为一元线性方程的形式,接着我们就要对得到的一些曲线回归方程进行选择,找出到底哪一个才是更好一点的。此时我们通常可采用两个指标进行选择,第一个是决定系数,第二个是剩余标准差。进而就得到了我们想要的定量关系表达式。第二题:给出

3、了某地区1971—2000年的人口数据,对该地区的人口变化进行曲线拟合。也用到了一元非线性回归的方法。首先我们也要对数据进行分析,描出数据的散点图,然后用MATLAB编程进行回归分析拟合计算输出利用 Logistic模型拟合曲线。关键词:参数估计, Logistic模型,MATLAB正文一、一元非线性回归分析的求解思路:文案大全实用标准文档•求解函数类型并检验。•求解未知参数。可化曲线回归为直线回归,用最小二乘法求解;可化曲线回归为多项式回归。二、回归曲线函数类型的选取和检验1、直接判断法2、作图观察法,与典型曲线比较,

4、确定其属于何种类型,然后检验。3、直接检验法(适应于待求参数不多的情况)4、表差法(适应于多想式回归,含有常数项多于两个的情况)三、化曲线回归为直线回归问题用直线检验法或表差法检验的曲线回归方程都可以通过变量代换转化为直线回归方程,利用线性回归分析方法可求得相应的参数估计值。题目:例8.5.1炼钢厂出钢水时用的钢包,在使用过程中由于钢水及炉渣对耐火材料的浸蚀,其容积不断增大。现在钢包的容积用盛满钢水时的重量y(kg)表示,相应的试验次数用x表示。数据见表8.5.1,要找出y与x的定量关系表达式。表8.5.1钢包的重量y与

5、试验次数x数据序号xy序号xy12106.42811110.5923108.20914110.6034109.581015110.9045109.501116110.7657110.001218111.0068109.931319111.20710110.49   1)1/y=a+b/x•y=a+blnx文案大全实用标准文档formatlongx=[23457810111415161819];y=[106.42108.20109.58109.5110109.93110.49110.59110.60110.9110.761

6、11111.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,‘

7、k+’,x,yc,‘r’)%数据的散点图和回归曲线图legend('散点图','回归函数')b=0.008966629680570.00082917436336R2=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.

8、/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

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。