气体的流量与压力之间的回归分析

1、学院理学院专 业信息与计算科学学生姓名班级学号课程名称概率论与数理统计课程设计课程设计（论文）题目气体的流量与压力之间的回归分析设计要求（技术参数）：设计通过该课程，使学生进一步理解概率论与数理统计的基本概念、理论和方法；初步掌握Excel统计工作表在随机模拟中是应用，MATLAB统计软件包作常见的统计检验和统计分析；具备初步的运用计算机完成数据处理的技能，使课堂中学习到理论得到应用。设计任务（至少一个）：1数据整理：收集数据，录入数据，画出相应图形；建立数学模型，数据的输入与整理，各种数据的图形显示。2一元、多元线性回归模型：回归系数的估计与检验，数据散点与回归直线的图示，残差图。运用MAT

2、LAB统计软件，对给定的数据拟合回归方程。3单因素、多因素方差分析：正态总体的方差分析问题； MATLAB统计软件中关于方差分析的相关命令，做出方差分析表，box图，能对结果进行简单分析。4假设检验： MATLAB绘制出直方图，做数据分布的推测；参数估计，假设检验，绘制概率密度图。计划与进度安排：周三12节：选题，设计解决问题方法 周三58节：调试程序 周五38节：完成论文，答辩成绩：指导教师（签字）： 2010年7月8日 专业负责人（签字）：2010年7月18日主管院长（签字）： 2010年7月19日摘要数理统计是具有广泛应用的数学分支，在生产过程和科学实验中，总会遇到多个变量，同一过程中的

3、这些变量往往是相互依赖，相互制约的，也就是说他们之间存在相互关系，这种相互关系可以分为确定性关系和相关关系。变量之间的确定性关系和相关关系在一定条件下是可以相互转换的。本来具有函数关系的变量，当存在试验误差时，其函数关系往往以相关的形式表现出来相关关系虽然是不确定的，却是一种统计关系，在大量的观察下，往往会呈现出一定的规律性，这种函数称为回归函数或回归方程。回归分析是一种处理变量之间相关关系最常用的统计方法，用它可以寻找隐藏在随机后面的统计规律。确定回归方程，检验回归方程的可信度等是回归分析的主要内容。按回归模型类型可划分为线性回归分析和非线性回归分析。本文主要探究一元非线性回归分析，对于一元

4、非线性问题，可用回归曲线y = f(x)来描述。在许多情形下，通过适当的线性变换，可将其转换为一元线性回归问题。气体的流量与压力之间的关系一般由经验公式M= cpb，此公式不是线性方程，如果对其进行如下变换：对其两边同时取常用对数lgM = lg c + b*lg p令 y = lg M ,x = lg p, a = lg c,则上述经过公式可以变换成一元线性方程： y = a+bx用线性回归方法求出x与y之间的线性回归方程y = 0.5426x-0.2827，返回原来的M和p之间的经验公式可表示为M= 0.5216p0.5426，回归方程的显著性检测：a.相关系数检验，b.F检验。所求得的经

5、验公式非常显著。关键词：相互关系；回归分析；一元非线性回归分析；线性回归方程；显著性检测目录1 设计目的12 设计原理12.1 线性回归方程的建立12.2 相关系数检验2 3 设计题目54 实现过程54.1 回归方程建立54.2 回归方程显著性检验75 设计总结8致谢10参考文献101、设计目的在许多实际问题中，变量之间的关系并不是线性的，这时就应该考虑采用非线性回归分析模型。在进行非线性回归分析时，必须着重解决两方面的问题：一是如何确定非线性函数的具体形式，与线性回归不同，非线性回归函数有多种多样的具体形式，需要根所研究的实际问题的性质和实验数据特点作出适当的选择；二是如何估计函数中的参数，

6、非线性回归分析最常用的方法仍然是最小二乘法，但需要根据函数的不同情形，做适当的处理。使其转化为线性回归。是得问题简单化。熟悉Excel在试验数据处理中的应用，2、设计原理对于一元非线性问题，可用回归曲线 y = f(x)来描述。在许多情况下，通过适当的线性变换，可将其转化为一元线性回归问题。具体做法如下： 根据实验数据，在直角坐标系中画出散点图； 根据散点图，推测y 和x之间的函数关系； 选择适当的变换，使之变成线性关系； 用线性回归方法求出线回归方程； 返回原来的函数关系，得到要求的回归方程。如果凭借以往的经验和专业知识，预先知道变量之间存在一定的形式的非线性关系，上述的前两步可以省略；如果

7、预先不清楚变量之间的函数类型，则可以依据实验数据的特点活散点图来选择对应的函数表达式。需要指出的是，在一定的试验范围内，可能用不同的函数拟合实验数据，都可以得到显著性较好的回归方程，这时应该选择其中数学形式较简单的一种。在回归方程显著性检测时也是根据线性回归方程而做出的显著性检测。2.1线性回归方程的建立：设有一组实验数据，实验值为xi，yi (i = 1,2,，n),其中x是自变量，y是因变量。若x，y符合线性关系，则可以你何为直线方程，即： = a + bxi 其中a，b称为回归系数。如果将与之间的偏差称为残差，用ei表示，则有：= -显然，只有各残差平方值之和最小时，回归方程预实验值的拟

8、合程度最好。令SSe=Q=，其中xi,yi是已知实验值，故残差平方值SSe为a，b的函数，为使SSe值到达极小， 根据极值原理，只要将上式分别对a,b求偏导数,，并令其等于零，即可求得a，b之值，根据最小二乘法，可以得到对方程组求解可得回归系数a，b的计算式：a =-b式中，分别为试验值xi，yi(i= 1，2，3，n)的算术平均值。从式中可以看出，回归直线通过点(，)，为了方便计算。令： = -n()= - n，则b = 。2.2相关系数检验相关系数是用于描述变量x和y的线性相关程度的，常用r来表示。设有n(n>2)对试验值xi，yi(i = 1，2，n)，则相关系数的计算式为：r =

9、 比较回归系数b与相关系数r的计算式，可得：r = =b，所以r与b 有相同的符号。相关系数r的平方未决系数r2具有以下特点：(1);(2)如果，则表明x和y完全线性相关，这时x和y有精确的线性关系；(3) 大多数情况下，即x和y之间存在一定的线性关系，当r>0时称x与y正相关，这时直线的斜率为正值，y随x的增加而增加，当r <0时称x与y负线性相关，这时的斜率为负值，y随x的增加而减小；(4)r = 0时，则表明x和y没有线形关系，但并不意味着x与y之间不存在其它类型的关系，所以先关系数更加精确的说法以偶那个该市线性相关系数。上面的分析可知，相关系数r越接近1，x与y的线性先关程

10、度越高，然而r的大小为你呢个回答其值达到多大时，x与y之间才存在线性相关，采用线性关系才属合理，所以须对相关系数r进行显著性检验。对于给定的显著性水平，显著性检验要求时，才说明x与y之间存在密切的线性关系，其中rmin称为相关系数临界值，它与给定的显著性水平和实验数据组n(n>2)有关，可查表得出。 F检验F检验实际就是方差分析，包括一下的主要内容。(1) 离差平方和试验值yi(i = 1，2，n)之间存在差异，这种差异可用试验值yi与其算术平均数的偏差平方和来表示，称为总离差平方和，即：SST = - nLyy试验值yi的这种波动是由于两个因素造成的：一个是由的变化而引起的变化，它可以

11、用回归平方和来表示，即：SSR=它表示的是回归值与yi的算术平均值之间的偏差平方和；另一个因素是随机误差，它可以用残差平方和来表示，即 SSe= ,它表示的是试验值yi与对应的回归值之间的偏差平方和。显然，这三种平方和之间有下述关系：SSY=SSR+SSe 。(2)自由度总离差平方和SSy的自由度为： = n-1回归平方和SSR的自由度为： = 1残差平方和SSe的自由度为： = n -2显然三种自由度之间的关系为：(3)均方MSR= MSe=(4)F检验F = F服从自由度为(1,n-2)F分布。在给定的显著性水平下，从F分布表种查得。一般取0.05和0.01，1-表示检验的可靠程度。若F&

12、lt;,则称x与y没有明显的线性关系，回归方程不可信；若<F<,则称x与y有显著的线性关系，用“*”表示；若F>,则称x与y有十分显著的线形关系，用“*”表示。后两种情况说明y的变化主要是由于x的变化造成的。最后将计算结果列成方差分析表差异源SSdfMSF显著性回归误差1n -2MSR=SSRMSe = SSE/(n-2)F = 总和n-1残差分析与之间的偏差称为残差，表示为= -，它能提供许多有用途的信息。根据试验数值可以计算出参差的标准误差=如果实验的随机误差服从正态分布，则试验值落在之内的概率为95%，对于一次试验数值，残差标准差越小，说明曲线拟合的越好。3、设计题目气

13、体的流量与压力之间的关系一般由经验公式表示，式中M是压强为P时每分钟流过流量计的空气物质的量，c,b为常数，现进行一批试验，得到如表 课设-1所示的一组数据。试由这组数据定出常数c,b，建立M和p之间的经验关系式，并检测其显著性。（=0.05） P/atm2.011.781.751.731.681.621.401.360.930.53M/(mol/min)0.7630.7510.7100.6950.6980.6730.6300.6120.4980.371课设-1 实验数据4、实现过程y = a+bx已知试验次数 n = 10，根据上述变换，使用Excel对试验数据进行整理计算，在Excel中建

14、立如图课设-1的工作4.1回归方程的建立经验公式不是线性方程，如果对其两边同时取常用对数，可得lgM = lg c + b*lg p如果令 y = lg M ,x = lg p, a = lg c,则上述经过公式可以变换成一元线性方程： 表格。 从Excel的常用工具栏中的“图表向导”按钮进入图表向导。做出XY散点图如课设-2所示， 在Excel中建立如图课设-3所示的工作表格。并通过Excel中自带的数据处理工具进行相关数据的计算课设-3序号lglg12.010.7630.303196-0.117480.0919280.0138-0.0356221.780.7150.25042-0.1456

15、90.062710.021227-0.0364831.750.710.243038-0.148740.0590670.022124-0.0361541.730.6950.238046-0.158020.0566660.024969-0.0376151.680.6980.225309-0.156140.0507640.024381-0.0351861.620.6730.209515-0.171980.0438970.029579-0.0360371.40.630.146128-0.200660.0213530.040264-0.0293281.360.6120.133539-0.213250.0

16、178330.045475-0.0284890.930.498-0.03152-0.302770.0009930.091670.009542100.530.371-0.27572-0.430630.0760240.1854390.11873414.796.3651.44195-2.045360.4812350.498928-0.14661.4790.63650.144195-0.20454实验数据散点图00.10.20.30.40.50.60.70.80.90.511.522.5p/atmM/(mol/min)课设-2从表课设-3，可知： = 0.1442, = -0.2045,= 0.481

17、2, = 0.4989， = -0.1466 = -n()-0.4812-100.1442=0.2733L = - n = -0.1466 - 10(-0.2045)=0.1483L= -n()=0.4989-10(-0.2045) = 0.0807 b = = 0.5426a = - b0.1442= -0.2827x与y之间的线性方程为：y = 0.5426x 0.2827又a = lgcc= 10= 10-0.2827 = 0.5216气体的流量M与压强p之间的经验公式可表示为：M = 0.5216 p0.54264.2回归方程显著性检验a.相关系数检验r = = = 0.9986根据=

18、 0.05，n = 10，m = 1查相关系数临界值表你，得rmin= 0.6319<r,所以求得的经验公式有意义。b.F检验已知： =0.2733L=0.1483SST = Lyy = 0.0807SSR = bLxy = 0.54260.1483 = 0.0805SSe = SST SSR = 0.0807 0.0805 = 0.0002方差分析结果见方差分析表课设-4。课设-4 方差分析表差异源SSdfMSFF0.01（1,8）显著性回归残差0.08050.0002180.08050.000025322011.3*总和0.08079所求得的经验公式非常显著。5、 设计总结结果分析：在本例中，非线性回归分析问题经过变换，转化成线形回归分析，求解的回归方程很好的拟合了试验数据值，在一定的试验的范围内，可能用不同的函数拟合试验数据，都可以得到显著性较好的回归方程，这时就应该选择其