欢迎来到天天文库
浏览记录
ID:56778001
大小:300.50 KB
页数:9页
时间:2020-07-09
《数学建模练习题.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、一.摘要在分析和研究了这种水泥凝固时放出的热量与这种水泥的四种化学成分有关,通过对所给的数据研究之后,提出了简单的多元线性回归模型,且在通过多种方法建立了回归模型,综合这几种方法建立的多元线性回归模型解决了我们面临的实际问题。模型求解和模型检验的结果表明,我们建立的模型是非常符合所求解的问题的,而且简单易懂,可操作性较高。以下这个方程为上述模型的结果:多元线性回归模型:y=62.4054+1.5511x1+0.5102x2+0.1019x3-0.1441x4关键词:多元线性回归模型模型求解模型检验二.问题重述题目1:某种水泥在凝固时放出的热量y(卡/克)与水泥种的下列
2、四种化学成分有关:X1:3CAO.AL2O3的成分(%)X2:3CAO.SiO2的成分(%)X1:4CAO.AL2O3.Fe2O3的成分(%)X1:2CAO.SiO2的成分(%)考虑y对这四个变量的线性回归,其实验数据如下表:编号(NO.)X1热量百分比X2热量百分比X3热量百分比X4Y卡克172666078.52129155274.331156820104.34113184787.6575263395.961155922109.27371176102.78131224472.59254182293.1102147426115.911140233483.8121166
3、912113.3131068812109.4请用四种方法为发热量建立回归方程。三.模型假设1.水泥在凝固时放出热量为固定值,收集的数据准确无误2.假设x1,x2,x3,x4为自变量,y为因变量3.假设y与诸x之间的线性关系可实际表示为4.是实际回归常数,是实际回归系数(j=1.2.3.4……)e是回归余四.问题分析与模型准备1.问题分析回归分析法是一种处理变量间相关关系的数理统计方法,不仅可以提供变量间相关关系的数学表达式,而且可以利用概率统计知识对此关系进行分析,以判别其有效性;还可以利用关系式,由一个或多个变量值,预测和控制另一个因变量的取值,进一步可以知道这种预
4、测和控制达到了何种程度,并进行因素分析。回归分析法就是以统计回归概念为基础,采用多种类型的回归法建立预测方程,包括一元线性、多元线性、非线性等。多元线性回归时要确定因变量与多个自变量之间的定量关系,它的数学模型为:其中,,,为待定参数;为随机变量,是除x以外其他随机因素对y影响的总和。其中,称E(y)=b0+b1x1+…+bmxm为理论回归方程。在实际问题的研究中,事先并不能断定随机变量y与变量x1,x2,…,xm之间是否有线性关系,在进行回归参数的估计前,用多元线性回归方程去拟合随机变量y与变量x1,x2,…,xm之间的关系,只是根据一些定性分析所作的一种假设。因此
5、,当求出线性回归方程后,还需对回归方程进行显著性检验,一般采用两种统计方法对回归方程进行检验,一种是回归方程显著性的F检验;另一种是回归系数显著性的t检验。2.图一给出所要求的数据考虑y对这四个变量的线性回归,其实验数据如下表:编号(NO.)X2热量百X3热量百X4Y卡克X1热量百分比分比分比172666078.52129155274.331156820104.34113184787.6575263395.961155922109.27371176102.78131224472.59254182293.1102147426115.911140233483.812116
6、6912113.3131068812109.4五.模型建立5.1最小二乘法参数估计多元线性回归模型为对上式求偏导并令其为0,然后求含n个未知参数的线性方程组,可以得到所需的参数估计值。得到:整理上式得:求解上面的方程组可求出,,,,得值。5.2方程组法进行参数估计可将问题写成y=xb+e其中式中:b是特定参数向量,e是残差向量遵从正态分布,即e~N(0,);I为n×n单位阵,满足E(e)=0,E(e)=X为已知的n×k+1常数矩阵,有自变量的各个观察值构成。Y是已知的n×1常数矩阵,有因变量的各个观察值构成。用最小二乘法可得到参数向量b的估计值为b=式中:是矩阵x的转
7、置矩阵,是矩阵的逆矩阵。回归系数,,,,通过偏导得到他们最小二乘估计值!5.3MATLAB软件参数估计由已知框图可编写以下MATLAB程序:结果说=明:b为回归模型中的常数系数和回归线系数;bint为个系数的95%置信区间,R与rint为对应实际值的残差和残差置信区间。Stats向量值为拟合优度·f值与显著性概率p。所以此函数方程为y=62.4054+1.5511x1+0.5102x2+0.1019x3-0.1441x4=0.9824说明模型拟合程度非常高。t检验:在回归模型区域中,给出了回归系数,,,,的估计值及其标准误差、检验值和回归系数估计区间
此文档下载收益归作者所有