数学模型报告新作_第1页
数学模型报告新作_第2页
数学模型报告新作_第3页
数学模型报告新作_第4页
数学模型报告新作_第5页
已阅读5页,还剩4页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

9 实 验 报 告课 程 数学模型实验 系 别 数统学院 班 级 2012级 学 号 2012325433 姓 名 王静怡 实验报告一:Matlab解决预测问题1实验课程名称 数学实验2实验项目名称 Matlab解决预测问题3 实验目的和要求 了解用Matlab软件求解图论模型及层次分析模型的方法。4 实验内容和原理内容:某公司想用全行业的销售额作为自变量来预测公司的销售额,下表给出了19771981年的公司销售额和行业销售额的分季度数据(单位:百万元).(1) 画出数据的散点图,观察用线性回归模型拟合是否合适.(2) 建立公司销售额对全行业销售额的回归模型,并用DW检验诊断随机误差项的自相关性.(3) 建立消除了随机误差项自相关性后的回归模型.年季t公司销售额y行业销售额x年季t公司销售额y行业销售额x19771120.96127.3197931124.54148.32221.40130.041224.30146.43321.96132.7198011325.00150.24421.52129.421425.64153.119781522.39135.031526.36157.32622.76137.141626.98160.73723.48141.2198111727.52164.24823.66142.821827.78165.619791924.10145.531928.24168.721024.01145.342028.78171.75主要仪器设备 计算机与Windows 2000/XP系统;Matlab等软件。6操作方法与实验步骤步骤:1.(1)打开Matlab,新建file-M文件。 (2)在M文件编辑窗口输入以下程序。命令如下:y=20.96;21.40;21.96;21.52;22.39;22.76;23.48;23.66;24.10;24.01;24.54;24.30;25.00;25.64;26.36;26.98;27.52;27.78;28.24;28.78;x=127.3;130.0;132.7;129.4;135.0;137.1;141.2;142.8;145.5;145.3;148.3;146.4;150.2;153.1;157.3;160.7;164.2;165.6;168.7;171.7;n=length(y);plot(x,y,o);b,bint,r,rint,stats=regress(y,ones(n,1) x);b bintstatspauserstool(x,y);pausee=y-b(1)-b(2)*x;et=zeros(n-1,1);et_1=zeros(n-1,1);for i=1:n-1 et(i)=e(i+1); et_1(i)=e(i);endplot(et_1,et,o)pauseDW=0;sum2=0;for i=2:n DW=DW+(e(i)-e(i-1)2; sum2=sum2+e(i)*e(i);endDW=DW/sum2pauserho=1-DW/2;yh=zeros(n-1,1);xh=zeros(n-1,1);for i=2:n yh(i-1)=y(i)-rho*y(i-1); xh(i-1)=x(i)-rho*x(i-1);endbh,bhint,rh,rhint,statsh=regress(yh,ones(n-1,1) xh);bh,bhintstatshpauserstool(xh,yh);pausee=yh-bh(1)-bh(2)*xh;et=zeros(n-1,1);et_1=zeros(n-1,1);for i=1:n-2 et(i)=e(i+1); et_1(i)=e(i);endpauseDW=0;sum2=0;for i=2:n-1 DW=DW+(e(i)-e(i-1)2; sum2=sum2+e(i)*e(i);endDW=DW/sum2图像及结果,执行代码可以得出如下结果:ans = -1.4548 -1.9047 -1.0048 0.1763 0.1732 0.1793stats = 1.0e+004 * 0.0001 1.4888 0 0.0000(3)点击执行按钮,运行其代码。y对x的散点图如上拟合图如下所示:9实验二 Lingo的简单用法1. 实验课程名称 数学实验2. 实验项目名称 利用Lingo解决简单的线性问题3. 实验目的和要求 了解Lingo的基本知识,熟悉上机环境,掌握Lingo的基本用法解决简单的线性问题。4. 实验内容和原理内容:某储蓄所每天的营业时间是上午9:00到下午5:00。根据经验,每天不同时间段所需要的服务员数量如下:时间段9101011111212112233445服务员数量43465688储蓄所可以雇佣全时和半时两类服务员,全时服务员每天报酬100元,从上午9:00到下午5:00工作,但中午12:00到下午2:00之间必须安排1h的午餐时间。储蓄所每天可以雇佣不超过三名的半时服务员,每个半时服务员必须连续工作4h,报酬40元。问该储蓄所应如何雇佣全时和半时两类服务员?如果不能雇佣半时服务员,每天至少增加多少费用?如果雇佣半时服务员的数量没有限制,每天可以减少多少费用?原理:利用Lingo编写的程序得到最优解,得出问题的答案。5. 主要仪器设备计算机与Windows 2000/XP系统,Lingo软件。6.操作方法与实验步骤步骤:1:设储蓄所每天雇佣的全时服务员中以12:001:00为午餐时间的有x1名,以1:002:00为午餐时间的有x2名;半时服务员中从9:00,10:00,11:00,12:00,1:00开始工作的分别为y1,y2,y3,y4,y5名。 2:打开Lingo,输入以下程序Min=100*x1+100*x2+40*y1+40*y2+40*y3+40*y4+40*y5 ; ; ; ; ; ; ; ;且为整数。7. 实验结果和分析实验结果:1:经过Lingo最

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论