数学建模实验报告

本科学生综合性实验报告项目组长杨维学号0113239成员杨维专业电子信息工程班级电信111班实验项目名称统计方法回归分析建模实验指导教师及职称党建武（教授）开课学期2013至2014学年第一学期上课时间2013年11月26日一、实验设计方案实验名称统计方法回归分析建模实验实验时间20131126小组合作是否小组成员无实验目的及要求学习统计方法回归分析的思想和基本原理；掌握建立回归模型的基本步骤，明确回归分析的主要任务；熟悉MATLAB软件进行回归模型的各种统计分析；通过范例学习，熟悉统计分析思想和建立回归模型的基本要素。通过该实验的学习，使学生掌握回归分析的统计思想，认识面对什么样的实际问题可以建立回归模型，并且对回归模型作统计分析，同时使学生学会使用MATLAB软件进行回归分析和计算的基本命令，了解统计软件的功能和作用；熟悉处理大量数据的要领和方法。2实验场地及仪器、设备和材料实验场地H123仪器、设备和材料计算机，MATLAB软件3实验学时数与实验类型3个学时，综合性实验4实验思路（实验内容、数据处理方法及实验步骤等）A实验内容1线性回归模型的建立与分析步骤（问题假设模型参数估计模型检验确定最优回归方程预测）；2非线性回归模型的建立与分析步骤；3使用MATLAB命令对回归模型进行计算与分析（包括模型检验与预测）；4利用某些数值与图形对统计特征作定性分析。B数据处理方法通过MATLAB软件的使用，实现一些线性回归方程问题求解的功能，得出数据及分析结果C实验步骤1开启软件平台MATLAB，开启MATLAB编辑窗口；2根据微分方程求解步骤编写M文件3保存文件并运行；4观察运行结果数值或图形；5根据观察到的结果和体会写出实验报告。D实验要求与任务1某校60名学生的一次考试成绩如下937583939185848277767795948991888683968179977875676968848381756685709484838280787473767086769089716686738094797877635355（1）计算均值、标准差、极差、偏度、峰度，画出直方图；（2）检验分布的正态性；（3）若检验符合正态分布，估计正态分布的参数并检验参数。2混凝土的抗压强度随养护时间的延长而增加，现将一批混凝土作成12个试块，记录了养护日期X（日）及抗压强度Y（KG/CM2）的数据养护时间X234579121417212856抗压强度Y354247535965687376828699试求型回归方程。XBALN3电影院调查电视广告费用和报纸广告费用对每周收入的影响，得到下面的数据，试建立回归模型，并进行检验（写出模型检验的依据），并预测电视广告费用为1，报纸广告费用为6时的周收入（写出预测的程序指令）。每周收入9690959295959494电视广告费1520152533234225报纸广告费5020402530352530二、实验结果与分析1、实验目的、场地及仪器、设备和材料、实验思路等见实验设计方案2、实验现象、数据及结果第一题某校60名学生的一次考试成绩如下9375839391858482777677959489918886839681799778756769688483817566857094848382807874737670867690897166867380947978776353551计算均值、标准差、极差、偏度、峰度，画出直方图；2检验分布的正态性；3若检验符合正态分布，估计正态分布的参数并检验参数。解输入数据S193758393918584827776S277959489918886839681S379977875676968848381S475668570948483828078S574737670867690897166S686738094797877635355SS1S2S3S4S5S6SAVEDATAS保存在当前工作目录下调用LOADDATA1）程序MEANS均值STDS标准差DMAXSMINS极差SKEWNESSS偏度KURTOSISS峰度HISTS,10直方图结果ANS801000ANS97106D44ANS04682ANS31529505560657075808590951000246810122）NORMPLOTS检验分布的正态性556065707580859095001002005010025050075090095098099DATAPROBABILITYNORMALPROBABILITYPLOT初步可以断定为正态分布3）MUHAT,SIGMAHAT,MUCI,SIGMACINORMFITS估计正态分布的参数MUHAT801000SIGMAHAT97106MUCI775915826085SIGMACI82310118436H,SIG,CITTESTS,801检验参数H0SIG1CI775915826085可以认为学生的考试成绩服从正态分布，平均成绩为8012混凝土的抗压强度随养护时间的延长而增加，现将一批混凝土作成12个试块，记录了养护日期X（日）及抗压强度Y（KG/CM2）的数据养护时间X234579121417212856抗压强度Y354247535965687376828699试求型回归方程。XBALN解一对将要拟合的非线性模型，建立M文件VOLUMM如下XBAYLNFUNCTIONYHATVOLUMBETA,XYHATBETA1BETA2LOGX输入数据X234579121417212856Y354247535965687376828699BETA051求回归系数BETA,R,JNLINFITX,Y,VOLUM,BETA0BETA结果BETA210053195287即得回规模型为XYLN528719032预测和作图010203040506030405060708090100解二令LNXU，则XE原回归方程为UBUAY输入数据X234579121417212856ULOGXUONES12,1UY354247535965687376828699回归分析及检测B,BINT,R,RINT,STATSREGRESSY,UB,BINT,STATS结果B210058195285BINT194463225653188943201627STATS10E003000104706900000作残差图RCOPLOTR,RINT作残差图246810123210123RESIDUALCASEORDERPLOTRESIDUALSCASENUMBER预测及作图ZB1B2LOGX预测及作图PLOTX,Y,K,X,Z,RZ3454194246004807805243575900656391436953237254267633428046088607889961490102030405060304050607080901003电影院调查电视广告费用和报纸广告费用对每周收入的影响，得到下面的数据，试建立回归模型，并进行检验（写出模型检验的依据），并预测电视广告费用为1，报纸广告费用为6时的周收入（写出预测的程序指令）。每周收入9690959295959494电视广告费1520152533234225报纸广告费5020402530352530解选择纯二次模型，即21210XXY方法一用多元二项式回归命令数据输入X11520152533234225X25020402530352530Y9690959295959494XX1X2回归、检测及预测RSTOOLX,Y,PUREQUADRATIC225335909294969810025335445BETA,RMSE结果BETA76301915280664540077906252RMSE02496故回归模型为221216507964552803976XXY剩余标准差为02496，显著性很好。预测电视广告费用为1，报纸广告费用为6时的周收入为951186修改RSTOOL函数MAXXMAXX1MINXMINX1在左边图形下方的方框中输入1，在右边图形下方的方框中输入6，则画面左边的“PREDICTEDY1”下方的数据即是周收入。1152253354458085909510010515225335445555方法二化为多元线性回归程序X11520152533234225X25020402530352530Y9690959295959494XONES8,1X1X2X12X22B,BINT,R,RINT,STATSREGRESSY,XB,STATS结果B76301915280664540077906252STATS09930107054000014006232混凝土的抗压强度随养护时间的延长而增加，现将一批混凝土作成12个试块，记录了养护日期X（日）及抗压强度Y（KG/CM2）的数据养护时间X234579121417212856抗压强度Y354247535965687376828699试求型回归方程。XBALN解一对将要拟合的非线性模型，建立M文件VOLUMM如下XBAYLNFUNCTIONYHATVOLUMBETA,XYHATBETA1BETA2LOGX输入数据X234579121417212856Y354247535965687376828699BETA051求回归系数BETA,R,JNLINFITX,Y,VOLUM,BETA0BETA结果BETA210053195287即得回规模型为XYLN528719032预测和作图010203040506030405060708090100解二令LNXU，则XE原回归方程为UBUAY输入数据X234579121417212856ULOGXUONES12,1UY354247535965687376828699回归分析及检测B,BINT,R,RINT,STATSREGRESSY,UB,BINT,STATS结果B210058195285BINT194463225653188943201627STATS10E003000104706900000作残差图RCOPLOTR,RINT作残差图246810123210123RESIDUALCASEORDERPLOTRESIDUALSCASENUMBER预测及作图ZB1B2LOGX预测及作图PLOTX,Y,K,X,Z,RZ3454194246004807805243575900656391436953237254267633428046088607889961490102030405060304050607080901003电影院调查电视广告费用和报纸广告费用对每周收入的影响，得到下面的数据，试建立回归模型，并进行检验（写出模型检验的依据），并预测电视广告费用为1，报纸广告费用为6时的周收入（写出预测的程序指令）。每周收入9690959295959494电视广告费1520152533234225报纸广告费5020402530352530解选择纯二次模型，即21210XXY方法一用多元二项式回归命令数据输入X11520152533234225X25020402530352530Y9690959295959494XX1X2回归、检测及预测RSTOOLX,Y,PUREQUADRATIC225335909294969810025335445BETA,RMSE结果BETA76301915280664540077906252RMSE02496故回归模型为221216507964552803976XXY剩余标准差为02496，显著性很好。预测电视广告费用为1，报纸广告费用为6时的周收入为951186修改RSTOOL函数MAXXMAXX1MINXMINX1在左边图形下方的方框中输入1，在右边图形下方的方框中输入6，则画面左边的“PREDICTEDY1”下方的数据即是周收入。1152253354458085909510010515225335445555方法二化为多元线性回归程序X11520152533234225X25020402530352530Y9690959295959494XONES8,1X1X2X12X22B,BINT,R,RINT,STATSREGRESSY,XB,STATS结果B76301915280664540077906252STATS09930107054000014006232混凝土的抗压强度随养护时间的延长而增加，现将一批混凝土作成12个试块，记录了养护日期X（日）及抗压强度Y（KG/CM2）的数据养护时间X234579121417212856抗压强度Y354247535965687376828699试求型回归方程。XBALN解一对将要拟合的非线性模型，建立M文件VOLUMM如下XBAYLNFUNCTIONYHATVOLUMBETA,XYHATBETA1BETA2LOGX输入数据X234579121417212856Y354247535965687376828699BETA051求回归系数BETA,R,JNLINFITX,Y,VOLUM,BETA0BETA结果BETA210053195287即得回规模型为XYLN528719032预测和作图010203040506030405060708090100解二令LNXU，则XE原回归方程为UBUAY输入数据X234579121417212856ULOGXUONES12,1UY354247535965687376828699回归分析及检测B,BINT,R,RINT,STATSREGRESSY,UB,BINT,STATS结果B210058195285BINT194463225653188943201627STATS10E003000104706900000作残差图RCOPLOTR,RINT作残差图246810123210123RESIDUALCASEORDERPLOTRESIDUALSCASENUMBER预测及作图ZB1B2LOGX预测及作图PLOTX,Y,K,X,Z,RZ3454194246004807805243575900656391436953237254267633428046088607889961490102030405060304050607080901003电影院调查电视广告费用和报纸广告费用对每周收入的影响，得到下面的数据，试建立回归模型，并进行检验（写出模型检验的依据），并预测电视广告费用为1，报纸广告费用为6时的周收入（写出预测的程序指令）。每周收入9690959295959494电视广告费1520152533234225报纸广告费5020402530352530解选择纯二次模型，即21210XXY方法一用多元二项式回归命令数据输入X11520152533234225X25020402530352530Y9690959295959494XX1X2回归、检测及预测RSTOOLX,Y,PUREQUADRATIC225335909294969810025335445BETA,RMSE结果BETA76301915280664540077906252RMSE02496故回归模型为221216507964552803976XXY剩余标准差为02496，显著性很好。预测电视广告费用为1，报纸广告费用为6时的周收入为951186修改RSTOOL函数MAXXMAXX1MINXMINX1在左边图形下方的方框中输入1，在右边图形下方的方框中输入6，则画面左边的“PREDICTEDY1”下方的数据即是周收入。1152253354458085909510010515225335445555方法二化为多元线性回归程序X11520152533234225X25020402530352530Y9690959295959494XONES8,1X1X2X12X22B,BINT,R,RINT,STATSREGRESSY,XB,STATS结果B76301915280664540077906252STATS09930107054000014006233、对实验现象、数据及观察结果的分析与讨论对于第1题的实验结果的分析，第1题是这4道题中最简单的，我是根据书上的例题及自己的理解将图形及运行代码写出来，直接输入命令语句以及画图语句就可将图形输出，操作较简单。对于第2道题的实验结果及数据的分析，虽然第2题的程序我自己按照书上的要求写的，但是系统总是显示有错误，原来我将主函数及输入命令放在了同一个窗口，所以系统总是显示错误，我只好将程序重新书写并做了修改，将他们分别放入不同的窗口，最后得出了实验结果，并正确得出了图形。对于第3题是这几道题目中算较难的，因为它所给的微分方程中含有3个参数，而它只给了我们其中的两个，第三个参数需要我们自己根据它的范围来了解方程解得变化，我在书本上都没有找到这样的例子，不知道该如何做，后来问了旁边的同学，才了解了这道题的意义，写出了代码，画出了图形，由于是需要需找规律，所以需要我们输入多个A的取值，得到多个图形，最后也得出了实验结果当A增大时，X，Y，Z慢慢不收敛而发散，当到一定值后而趋于稳定。对于第4道题4、结论通过本次实验的学习，让我掌握微分方程组求解方法（解析法、欧拉法、梯度法、改进欧拉法等），对常微分方程的数值解法有一个初步了解，同时学会使用MATLAB软件求解微分方程的基本命令，学会建立微分方程方面的数学模型。这对于我想深入理解微分、积分的数学概念，掌握数学的分析思维方法，熟悉处理大量的工程计算问题的方法是很有帮助的。5、实验总结（1）本次实验成败之处及其原因分析本实验的成功之处在于对微分方程组模型及求解微分方程的基本命令等的知识的运用；失败之处在于对MATLAB软件掌握的还不是特别的熟练以及对函数调用的方法还不是特别熟悉，尤其是函数调用这方面的知识还不是特别的熟悉，因此要加强这方面的练习。（2）本实验的关键环节及改进措