数学模型实验4个讲解Word文件下载.docx
《数学模型实验4个讲解Word文件下载.docx》由会员分享,可在线阅读,更多相关《数学模型实验4个讲解Word文件下载.docx(16页珍藏版)》请在冰豆网上搜索。
安装了VC++、mathematica、matlab的计算机
实验日期:
2015年3月13日实验地点:
第五教学楼北802
实验目的:
掌握插值与拟合的原理,熟悉插值与拟合的软件实现。
实验准备:
1.在开始本实验之前,请回顾教科书的相关内容;
2.需要一台准备安装WindowsXPProfessional操作系统和装有VC++6.0的计算机。
实验内容及要求
下表给出了某工厂产品的生产批量与单位成本(元)的数据。
从散点图可以明显地发现,生产批量在500以内时,单位成本对生产批量服从一种线性关系,生产批量超过500时服从另一种线性关系,此时单位成本明显下降。
生产批量
650
340
400
800
300
600
单位成本
2.48
4.45
4.52
1.38
4.65
2.96
720
480
440
540
750
2.18
4.04
4.20
3.10
1.50
要求:
1、构造合适的模型全面地描述生产批量与单位成本的关系;
2、对于这种关系,试采用分段函数进行详细分析。
另外,从误差的角度出发,定量与定性相结合的方式来说明采用分段函数来描述这种关系的优点。
实验过程:
运用matlab进行以下程序运行:
n=[300340400440480540600650720750800];
t=[4.654.454.524.204.043.102.962.482.181.501.38];
plot(n,t,'
bo:
'
)
holdon
x=log(n);
y=log(t);
c=polyfit(x,y,1);
a=exp(c
(2));
b=c
(1);
a,b
n1=300:
20:
480;
t1=a*n1.^b;
plot(n1,t1,'
r-'
n1=540:
800;
b-'
实验总结(由学生填写):
对数学模型有了初步的了解,大致简单的掌握了MATLAB的用法
以及和数学模型的关系。
对学习数学模型有了很大的帮助。
实验等级评定:
数学规划模型(实验二)指导教师:
4实验设备:
2015年3月24日实验地点:
第五教学楼北802
掌握优化问题的建模思想和方法,熟悉优化问题的软件实现。
1.在开始本实验之前,请回顾教科书的相关内容;
2.需要一台准备安装WindowsXPProfessional操作系统和装有数学软件的计算机。
一家保姆服务公司专门向雇主提供保姆服务。
根据估计,下一年的需求是:
春季6000人日,夏季7500人日,秋季5500人日,冬季9000人日。
公司新招聘的保姆必须经过5天的培训才能上岗,每个保姆每季度公司(新保姆包括培训)65天。
保姆从该公司而不是从雇主那里得到报酬,每人每月工资800元。
春季开始时公司拥有120名保姆,在每个季度结束后,将有15%的保姆自动离职。
(1)如果公司不允许解雇保姆,请你为公司制定下一年的招聘计划;
哪些季度需求的增加不影响招聘计划?
可以增加多少?
(2)如果公司在每个季度结束后允许解雇保姆,请为公司制定下一年的招聘计划。
1,假设每个季度新招的人数分别为:
x1,x2,x3,x4
每个季度拥有的人数分别为:
s1,s2,s3,s4
根据题意在lingo中输入:
min=s1+s2+s3+s4;
65*s1>
=6000+5*x1;
65*s2>
=7500+5*x2;
65*s3>
=55005*x3+;
65*s4>
=9000+5*x4;
s1-x1=120;
s1-0.15*s1+x2=s2;
s2-0.15*s2+x3=s3;
s3-0.15*s3+x4=s4;
end
结果:
VariableValueReducedCost
S1120.00000.000000
S2116.50000.000000
S399.025000.000000
S4142.98570.000000
X10.0000000.8732231
X214.500000.000000
X30.0000000.9291667
X458.814480.000000
根据结果可以看出公司在春季和秋季不招人,在夏季招15人,在冬季招59人。
2,假设公司每个季度解雇的人数分别为:
y1,y2,y3
在lingo中输入:
model:
=5500+5*x3;
s1-0.15*s1-y1+x2=s2;
s2-0.15*s2-y2+x3=s3;
s3-0.15*s3-y3+x4=s4;
结果;
VariableValueReducedCost
S384.615380.000000
S4144.00640.000000
X10.0000000.9291667
X30.0000000.7147436E-01
X472.083330.000000
Y10.0000000.8333333E-01
Y214.409620.000000
Y30.0000000.8333333E-01
可以看出公司在春季不招人也不解雇人,在夏季招15人,解雇15人,在秋季不招人也不解雇人,在冬季只招73人。
通过该实验,有掌握了数学模型中的一类题型,并且学会了使用lingo,该实验的实际应用性很强,对线性规划问题也有了一定的了解。
微分方程模型(实验三)
指导教师:
4
实验设备:
年
月
日
实验地点:
第五教学楼北802
通过本次实验,让学生巩固微分方程的理论知识,掌握用微分方程来建立数学模型,会运用微分方程解决实际问题,并学会如何用Matlab软件来求解微分方程模型.
2.需要一台准备安装WindowsXP
Professional操作系统和装有Matlab软件的计算机。
意大利生物学家Ancona曾致力于鱼类种群相互制约关系的研究,他从第一次世界大战期间,地中海各港口捕获的几种鱼类捕获量百分比的资料中,发现鲨鱼等的比例有明显增加(见下表),而供其捕食的食用鱼的百分比却明显下降。
显然战争使捕鱼量下降,食用鱼增加,鲨鱼等也随之增加,但为何鲨鱼的比例大幅增加呢?
他无法解释这个现象。
希望你建立一个食饵—捕食系统的数学模型,定量地回答这个问题。
1.1符号说明:
x1(t)-------食饵在t时刻的数量;
x2(t)--------捕食者在t时刻的数量;
r1------食饵独立生存时的增长率;
r2-----捕食者独立存在时的增长率;
m1-----捕食者掠取食物的能力;
m2-----食饵对捕食者的供养能力;
e------捕获能力系数;
1.2
模型假设
1:
食饵由于捕食者的存在始增长率降低,假设降低的程度与捕食者的数量成正比。
2:
捕食者由于食饵为他提供食物的作用使其死亡率降低或使之增长,假定增长的程度与食饵数量成
正比
1.3模型建立
不考虑人工捕获
dx1/dt=x1(r1-m1x2)
dx2/dt=x2(-r2+m2x1)
设食饵和捕食者的初始数量为x1(0)=x10,x2(0)=x20
对于数据r1=1,m1=0.1,r2=0.5,m2=0.02,x10=25,x20=2,t的终值经试验后确定为15,即:
x1=x1(1-0.1x2)
x2=x2(-0.5+0.02x1)
x1(0)=25,x2(0)=2
1.4
模型求解
function
dx=shier(t,x)
dx=zeros(2,1);
dx
(1)=x
(1)*(1-0.1*x
(2));
dx
(2)=x
(2)*(-0.5+0.02*x
(1));
[t,x]=ode45('
shier'
[015],[252]);
>
plot(t,x(:
1),'
-'
t,x(:
2),'
*'
plot(x(:
1),x(:
2))
结果如下:
上图反应了x1(t)与x2(t)的关系。
可以猜测x1(t)与x2(t)都是周期函数。
考虑人工捕获:
设表示捕获能力的系数为e,相当于食饵的自然增长率由r1降为r1-e,捕食者的死亡率由r2
增为r2+e
dx1/dt=x1[(r1-e)-m1x2]
dx2/dt=x2[-(r2+e)+m2x1]
初始取值仍一样
设战前捕获能力系数e=0.3,战争中降为0.1,则模型为:
dx1/dt=x1(0.7-0.1x2)
dx2/dt=x2(-0.8+0.02x1)
dx1/dt=x1(0.9-0.1x2)
dx2/dt=x2(-0.6+0.02x1)
用以上方法可以求解这两个方程,即可以得出以下结论:
战争中的鲨鱼比例比战前高
通过该实验,我们掌握了微分方程的求法以及新的一个数学模型的题型,使我们在学习数学模型的路上又前进了一大步。
层次分析法(实验四)指导教师:
安装了VC++、mathematica、matlab的计算机
年月日实验地点:
第五教学楼北902
熟悉有关层次分析法模型的建立与计算,熟悉Matlab的相关命令。
2.需要一台准备安装WindowsXPProfessional操作系统和装有Matlab的计算机。
试用层次分析法解决一个实际问题。
问题可参考教材P296第4大题。
1.1问题分析
用层次分析法解决一个实际问题:
你要购置一台个人电脑,考虑功能、价格等的因素,如何做出决策。
1.2模型假设
假设有联想、华硕、戴尔、三星四款电脑供你选择,你会根据诸如功能、价格、性能、外观、口碑等一些准则来反复比较这四款电脑
1.3模型建立
将决策问题分解为3个层次:
目标层A,表示“购买电脑”;
准则层B,表示购买电脑所考虑的因素,有价格、售后、配置、外观、功能等;
方案层C,表示所购买的电脑的品牌,联想、华硕、戴尔、三星。
可得其层次结构为:
1.4模型求解及结果解释
构造成对比较矩阵
B对A:
C对B:
计算权向量并做一致性检验:
对于Mba,有
,可得
,归一化的特征向量
。
则一致性指标
,所以该矩阵为一致阵。
同理,对于Mcb1、Mcb2、Mcb3、Mcb4和Mcb5,可得第三层计算结果,如下表:
k
1
2
3
4
5
wk(3)
λk
CIk
当n=4时随机一致性指标RI=0.90,所以上面的CIk均可通过一致性检验。
计算组合权向量并做组合一致性检验
对于“购买联想电脑”方案,其对于目标的组合权重为
同样可以算出“购买华硕电脑”、“购买戴尔电脑”、“购买三星电脑”在目标中的组合权重为0.213、0.132、0.076所以组合权向量
结果表明“购买联想电脑”方案在电脑品牌选择中超过一半,远大于其他方案,所以应作为第一选择方案。
由以上数据可得第3层的一致性指标为CI(3)=0,随机一致性指标为RI(3)=0.90,组合一致性比率为
,所以第3层对第1层的组合一致性比率为CR*=0,组合一致性检验通过。
明确用层次分析法来做本题,确定目标层、准则层和方案层,
构造成对比较矩阵,分别计算权向量和组合权向量并做一致性检验