《完整word版matlab回归分析方法》由会员分享，可在线阅读，更多相关《完整word版matlab回归分析方法（17页珍藏版）》请在金锄头文库上搜索。
1、第八章回归分析方法当人们对研究对象的内在特性和各因素间的关系有比较充分的认识时，一般用机理分析方法建立数学模型。如果由于客观事物内部规律的复杂性及人们认识程度的限制，无法分析本章讨论其中用途非常广泛的一类模型一一统计回归模 生产工艺优化等问题。一类叫确定性关系， 也叫函数关系，其特征是：一个变 -类关系叫相关关系， 变量之间的关系很难用一种精确的 但人的年龄和血压之间没有确定的数量 回归分析就是处理变量之间的相关关系的实际对象内在的因果关系，建立合乎机理规律的数学模型，那么通常的办法是搜集大量数据， 基于对数据的统计分析去建立模型。 型。回归模型常用来解决预测、控制、变量之间的关系可以分为两类
2、： 量随着其它变量的确定而确定。另 方法表示出来。例如，通常人的年龄越大血压越高， 关系，人的年龄和血压之间的关系就是相关关系。 一种数学方法。其解决问题的大致方法、步骤如下：(1) 收集一组包含因变量和自变量的数据；(2) 选定因变量和自变量之间的模型，即一个数学式子，禾U用数据按照最小二乘准则 计算模型中的系数；(3) 利用统计分析方法对不同的模型进行比较，找出与数据拟合得最好的模型；(4) 判断得到的模型是否适合于这组数据；(5) 利用模型对因变量作出预测或解释。应用统计分析特别是多元统计分析方法一般都要处理大量数据，工作量非常大，所以在计算机普及以前，这些方法大都是停留在理论研究上。运
3、用一般计算语言编程也要占用大量时间，而对于经济管理及社会学等对高级编程语言了解不深的人来说要应用这些统计方法更 是不可能。MATLAB等软件的开发和普及大大减少了对计算机编程的要求，使数据分析方 法的广泛应用成为可能。MATLAB统计工具箱几乎包括了数理统计方面主要的概念、理论、方法和算法。运用 MATLAB统计工具箱，我们可以十分方便地在计算机上进行计算，从而 进一步加深理解，同时，其强大的图形功能使得概念、过程和结果可以直观地展现在我们面 前。本章内容通常先介绍有关回归分析的数学原理，主要说明建模过程中要做的工作及理由，如模型的假设检验、参数估计等，为了把主要精力集中在应用上，我们略去详细
4、而繁杂的理论。在此基础上再介绍在建模过程中如何有效地使用MATLAB软件。没有学过这部分数学知识的读者可以不深究其数学原理，只要知道回归分析的目的，按照相应方法通过软件显示的图形或计算所得结果表示什么意思，那么，仍然可以学到用回归模型解决实际问题的基本方法。包括：一元线性回归、多元线性回归、非线性回归、逐步回归等方法以及如何利用 MATLAB软件建立初步的数学模型，如何透过输出结果对模型进行分析和改进，回归模型 的应用等。&1 一元线性回归分析回归模型可分为线性回归模型和非线性回归模型。非线性回归模型是回归函数关于未知参数具有非线性结构的回归模型。某些非线性回归模型可以化为线性回归模型处理；如
5、果知道函数形式只是要确定其中的参数则是拟合问题，可以使用MATLAB软件的curvefit命令或nlinfit命令拟合得到参数的估计并进行统计分析。本节主要考察线性回归模型。8.1.1 一元线性回归模型的建立及其MATLAB实现y = Po + Pix + S N(ob)其中00,卩1是待定系数，对于不同的 X, y是相互独立的随机变量。假设对于x的n个值Xi，得到y的n个相应的值yi，确定P1的方法是根据最小二乘准则，要使nQ(p0,p1)=2 Ei2i丄n=2 yi -(p0 + 陥)2i 4取最小值。利用极值必要条件令 孚=0,孚=0，求p0, P1的估计值（?0，f?，从而得到cP。C
6、r1回归直线y =氏+ ?x。只不过这个过程可以由软件通过直线拟合完成,而无须进行繁杂的运算。（1）参数的区间估计由于我们所计算出的（?0, ?仍然是随机变量，因此要对（?0, （?取值的区间进行估计,如果区间估计值是一个较短的区间表示模型精度较高。（2）对误差方差的估计 设？为回归函数的值，yi为测量值，残差平方和剩余方差s2 =Qn 2（3）线性相关性的检验由于我们采用的是一元线性回归，因此，如果模型可用的话，应该具有较好的线性关系。反映模型是否具有良好线性关系可通过相关系数R的值及F值观察（后面的例子说明）。（4）一元线性回归的 MATLAB 实现MATLAB工具箱中用命令regress
7、实现，其用法是：b=regress（y,x）b ,bint , r ,nnt , s=regress（y , x , alp ha）输入y （因变量，列向量）、x（ 1与自变量组成的矩阵，见下例），alpha是显著性水平（缺 省时默认0.05）。输出b =（?）, f?），注意：b中元素顺序与拟合命令polyfit的输出不同，bi nt是卩i的置信区间，r是残差（列向量），rint是残差的置信区间，s包含4个统计量：决定系数 R2 （相2关系数为R）; F值；F（1,n-2）分布大于F值的概率p；剩余方差s的值（MATLAB7.0以后版本）。s2也可由程序sum（rA2）/（n-2）计算。其意
8、义和用法如下： R2的值越接近1,变量的线性相关性越强，说明模型有效；如果满足F1t（1, n- 2KF，则认为变量y与x显著地有线性关系，其中F（1, n-2）的值可查F分布表，或直接用 MATLAB命令finv(1-a ,1, n-2)计算得到；如果 p吒a表示线性模型可用。这三个值可以相互印证。s2的值主要用来比较模型是否有改进，其值越小说明模型精度越 咼。8.1.2身高与腿长例1测得16名成年女子身高 y与腿长X所得数据如下：表8-116名女子身高(cm)腿长(cm)数据888588 9192939395969897969899100102143 145146147 149150153
9、154155156157158159160162164首先利用命令 plot(x,y,r*)画出散点图，从图形可以看出，这些点大致分布在一条直线 的左右，因此，可以考虑一元线性回归。可编制程序如下：y=143 145 146147149150x=888588919293n=16;X=o nes( n,1),x;b,b in t,r,ri nt,s=regress(y,X,0.05);b,b in t,s,rcopi ot(r,ri nt)运行后得到1531541551561571581591609395969897969899162164;100 102;b = 31.77131.2903bi
10、nt = 12.31961.084651.22291.4960s = 0.9282180.95310.00003.12772R =0.9282，由 finv(0.95,1,14)= 4.6001，即 F(1,n 2) = 4.6001F=180.9531 ,p0.0001，可以通过残差图发现，第二个数据为奇异数据，去掉该数据后运行后得到b = 17.65491.4363bint = -0.598635.90831.24451.6281s = 0.9527261.63890.00001.9313p0.0001，y =17.6549 +1.4363X。R2 =0.9527，由 fin v(0.95
11、,1,13)= 4.6672，即 F(1, n 2) = 4.6672m）得n组观察值，采用最小 二乘估计求得回归方程y=?o+f?xi+lli+f?kxm建立回归模型是一个相当复杂的过程，概括起来主要有以下几个方面工作（1）根据研究目的收集数据和预分析；（2）根据散点图是否具有线性关系建立基本回归模型；（3）模型的精细分析；（4 ）模型的确认与应用等。收集数据的一个经验准则是收集的数据量（样本容量）至少应为可能的自变量数目的610倍。在建模过程中首先要根据所研究问题的目的设置因变量，然后再选取与该因变量有统计关系的一些变量作为自变量。我们当然希望选择与问题关系密切的变量，mu 量相不可以在得
12、到初步的模型后利用MATLAB软件进行相关性检验。下面通过一个案例探讨 MA TLAB软件在回归分析建模各个环节中如何应用。多元线性回归的 MATLAB实现仍然用命令regress（y , X），只是要注意矩阵 X的形式，将通过如下例子说明其用法。8.2.2某类研究学者的年薪1.问题例2工薪阶层关心年薪与哪些因素有关，以此可制定出它们自己的奋斗目标。某科学基金会希望估计从事某研究的学者的年薪丫与他们的研究成果（论文、著作等）的质量指标Xi、从事研究工作的时间X2、能成功获得资助的指标X3之间的关系，为此按一定的实验设计方法调查了24位研究学者，得到如下数据（i为学者序号）：表8-2从事某种研究的学者的相关指标数据i123456789101112Xi13.55.35.15.84.26.06.85.53.17.24.54.9X29201833311325305472511x36.16.47.46.77.55.96.04.05.88.35.06.4yi33.240.338.746.841.437.539.040.730.152.938.231.8i131415161718192021222324Xi18.06.56.63.76.27.04.04.55.95.64.83.9xi2233539217403523
