数学建模实验1,2章

合集下载

数学建模实验报告(一)MATLAB中矩阵的基本操作

数学建模实验报告(一)MATLAB中矩阵的基本操作
Column 6
1.5270
j =
2 3 2 5 1 4
>> min(a,[],1)
ans =
Columns 1 through 5
-2.3299 -0.1303 -1.3617 -1.1176 -0.3031
Column 6
0.0230
>> min(a,[],2)
ans =
-0.4762
-0.0679
-2.3299 -0.1303 0.4550 -1.1176 -0.2176
-1.4491 0.1837 -0.8487 1.2607 -0.3031
Column 6
0.0230
0.0513
0.8261
1.5270
0.4669
>> size(a)
ans =
5 6
>> [i,j]=find(a==max(max(a)))
-1.0000 2.5000 1.0000
3.5000 5.5000 2.5000
>> X=D
X =
4.0000 1.5000 -1.0000
-1.0000 2.5000 1.0000
3.5000 5.5000 2.5000
5、利用randn(5,6)命令生成一个随机矩阵T,求T的矩阵大小,每一行、每一列的最大值和最小值,整个矩阵的最大值与最小值;然后将整个矩阵的最大值所在位置的元素换为100,将最小值所在位置的元素取为-100。
(2):>> a=[2 5 8;7 1 9]
a =
2 5 8
7 1 9
>> b=[4 2 1 3;0 7 6 2;-3 5 9 -1]

数学建模 第二篇1 MATLAB作图讲解

数学建模 第二篇1 MATLAB作图讲解

MATLAB作图
(2) mesh(x,y,z) 画网格曲面
数据矩阵。分别表示数据点 的横坐标、纵坐标、函数值
例 画出曲面Z=(X+Y).^2在不同视角的网格图. 解 x=-3:0.1:3;y=1:0.1:5; [X,Y]=meshgrid(x,y); Z=(X+Y).^2; mesh(X,Y,Z)
MATLAB作图
(2) figure(h) 新建h窗口,激活图形使其可见,并置于其它图形之上


区间[0,2*pi]新建两个窗口分别画出 y=sin(x);z=cos(x)。
x=linspace(0,2*pi,100); y=sin(x);z=cos(x); plot(x,y); title('sin(x)'); pause figure(2); plot(x,z); title('cos(x)'); 返回
hh = zlabel(string) hh = title(string)
MATLAB作图
例 在区间[0,2*pi]画sin(x)的图形,并加注图例 “自变量X”、“函数Y”、“示意图”, 并加格栅.
解 x=linspace(0,2*pi,30); y=sin(x); plot(x,y) xlabel('自变量X') ylabel('函数Y') title('示意图') grid on
3.图形保持 hold off 释放当前图形窗口
MATLAB作图
(1) hold on 保持当前图形, 以便继续画图 例 将y=sin(x),y=cos(x)分别用点和线画在一图上
解 x=linspace(0,2*pi,30); y=sin(x); z=cos(x) plot(x,z,:) hold on Plot(x,y) Matlab liti 5

数学建模 -实验报告1

数学建模 -实验报告1
推导出了动力学方程
������������⁄������������ = ������������(1 − (������ + ������)) − ������1������∗������,
(4 − 3)
������������∗⁄������������ = −������1������∗������ + ������2������
二、 问题分析
建立肿瘤细胞增长模型时,我们可以从自由增长模型开始分析,引进 Logistic 阻滞增长模型,构成肿瘤细胞增长初步框架。再者肿瘤细胞不同于普 通细胞,其生长受到人体自身免疫系统的制约。于是综合考虑正常细胞转化,癌 细胞增殖,癌细胞死亡,癌细胞被效应细胞消除等情况,建立动力学方程。并对 模型进行适当简化求解。在放射治疗方案的设计中,我们可以引入放射生物学中 广泛接受的 LQ 模型对问题进行分析,由于放疗对人体伤害相当大,因此我们采 取分次逐次放疗的方式进行治疗。我们具体分两种情形进行讨论,一是在总剂量 一定的条件下,不同的分次剂量组合对生物效应的影响;二是在产生相同生物效 应的情况下,分析最优的分次剂量组合。
易算出癌细胞转入活动期已有 300 多天,故如何在早期发现癌症是攻克癌症的关键之一 (2)手术治疗常不能割去所有癌细胞,故有时需进行放射疗法。射线强度太小无法杀
死癌细胞,太强病人身体又吃不消且会使病人免疫功能下降。一次照射不可能杀死全部癌细 胞,请设计一个可行的治疗方案(医生认为当体内癌细胞数小于 100000 个时即可凭借体内 免疫系统杀灭)。
进一步简化,根据(4-4),(4-5)式可知,效应细胞������∗和复合物������有出有进.假 设出入保持平衡,则有
������ + ������∗ = C (C 为常数)

数学建模精讲第2章 ppt课件

数学建模精讲第2章 ppt课件
qi=Npi /P不全为整数时,ni 应满足的准则: 记 [qi]– =floor(qi) ~ 向 qi方向取整; [qi]+ =ceil(qi) ~ 向 qi方向取整.
1) [qi]– ni [qi]+ (i=1,2, … , m), 即ni 必取[qi]– , [qi]+ 之一
2) ni (N, p1, … , pm ) ni (N+1, p1, … , pm) (i=1,2, … , m) 即当总席位增加时, ni不应减少
《数精学品课建程模》
“公平”分配方 将绝对度量改为相对度量 法若 p1/n1> p2/n2 ,定义
p1/n1p2/n2 p2/n2
rA(n1,n2)~
对A的相对不公平度 公平分配方案应
类似地定义 rB(n1,n2)
使 rA , rB 尽量小
将一次性的席位分配转化为动态的席位分配, 即
设A, B已分别有n1, n2 席,若增加1席,问应分给A, 还是B
模型建立
《数精学品课建程模》
建立t与n的函数关系有多种方法 1. 右轮盘转第 i 圈的半径为r+wi, m圈的总长度
等于录像带在时间t内移动的长度vt, 所以
m
2(r wi) vt
mk n
i1
twk 2 n2 2rkn
v
v
模型建立
《数精学品课建程模》
2. 考察右轮盘面积的 变化,等于录像带厚度
乘以转过的长度,即
Q1最大,第20席给甲系
第21席 Q11110123280.4, Q2, Q3同上
Q3最大,第 21席给丙系
Q值方法 分配结果
甲系11席,乙系6席,丙系4 席
公平吗?

1数学建模实验-圆周率的计算

1数学建模实验-圆周率的计算

11-23π
1 1 = 4 arctan − arctan (Machin公式) 公式) 公式 4 5 239
此式求得了π的第 位小数且全部正确 此式求得了 的第100位小数且全部正确 的第
12-23
方法1 利用幂级数表达式
1 2 4 n−1 2n−2 =1− x + x −L+ (−1) x +L 2 1+ x
4-23
古 典 方 法 分 析 方 法 其 它 方 法 Ä 概率方法 Ä 数值积分方法
5-23
古典方法
的近似值呢?显然, 用什么方法来计 算π的近似值呢?显然,不可能仅根 据圆周率的定义,用圆的周长去除以直径。起先, 据圆周率的定义,用圆的周长去除以直径。起先,人们 采用的都是用圆内接正多边形和圆外切正多边形来逼近 的古典方法。 的古典方法。
19-23
设计方案
在正方形 0< x <1, 0< y<1 上随机的投大量的点,那么 落在四分之一园内的点数 数m与在正方形内的点数n 之比m/n应为这两部分图形 面积之比=π/4,故 π=4 m/n 计算机模拟:产生区间[0,1]上数目为n的一组 随机数(x,y),计算满足x2+y2<1的点数m
3.1415926535897932384626433832795028841971 6939937510582097494459230781640628620899 8628034825342117068 但是你会计算π的值吗?你又能用几种 方法计算π的值?
3-23
π的计算 的计算
Ä 圆周率是人类获得的最古老的数学概念 之一,早在大约3700年前(即公元前1700 之一,早在大约3700年前(即公元前1700 3700年前 年左右) 256/81( 年左右)的古埃及人就已经在 用256/81( 约3.1605)作为π的近似值了。几千年来 3.1605)作为π的近似值了。 ,人们一直没有停止过求π的努力。 人们一直没有停止过求π的努力。

第二章-1-建模基本概念-电路-传递函数-方块图

第二章-1-建模基本概念-电路-传递函数-方块图

2
1 RCs RC 1
电路及组成
例2:电阻电感电容(RLC)串联电路
1 LDi Ri ie CD
uR
L 1 DuR uR e R RCD
d 2uc (t ) duc (t ) T1T2 T2 uc (t ) e 2 dt dt
• 上述方程是线性定常微分方程。由这种方程描述的系统又称为 线性时不变( linear time-invariant, LTI )系统。由二阶微 分方程描述的系统称为二阶系统。
的方块图。
U
ei
i
R
U
o
I
1
Cs
e0
1 U0 I , Cs
U
i
Ui Uo I R
I
1 Cs
1 R
U
o
传递函数
U o (s) 1 U i(s) RCs 1
电路及组成
一阶系统的阶跃响应
考察标号为***的方程( 称为一阶微分方程 )
de0 T e0 ei dt
控制轨迹
***
19
电路及组成
一阶系统的阶跃响应
y x
A KA 0.632KA
de0 T e0 ei d dt
***
y (t ) KA(1 e
T1 T2
t
T
)
t 时域响应分析: 当 t=0, y(0)=0, 当 t=T, , 当
t
dy dt
t
t 0

KA T dy dt 0
y (T ) KA(1 e 1 ) 0.632 KA
图 2.1
va
LD R LD
vb
16

数学建模实验教学大纲

数学建模实验教学大纲

数学建模实验教学大纲一、引言数学建模是一门涉及数学、计算机科学和实际问题解决的跨学科课程。

通过数学建模实验教学,学生将学习如何将实际问题抽象化、建立模型,并运用数学方法进行问题求解。

本教学大纲旨在为数学建模实验课程提供指导,帮助教师和学生达到教育目标。

二、课程目标1. 培养学生的科学思维和实际问题解决能力。

2. 掌握各种数学模型的建立与求解方法。

3. 学习数据分析技术和模型验证方法。

4. 提高学生的团队合作和沟通能力。

三、教学内容1. 数学建模的基础知识(1) 数学建模的定义和基本步骤。

(2) 常见数学模型的分类和特点。

2. 实际问题抽象化和模型建立(1) 学习如何从实际问题中提取关键信息。

(2) 学习如何建立数学模型,选择合适的数学方法和假设。

3. 数学模型求解(1) 学习常见数学方法的应用,如线性规划、微分方程等。

(2) 掌握数学软件工具的使用,如Matlab、Python等。

4. 数据分析和模型验证(1) 学习数据收集和处理的基本技巧。

(2) 学习如何验证数学模型的准确性和可靠性。

5. 团队合作和沟通(1) 学习如何分工合作,形成高效的团队。

(2) 提高表达和演示能力,培养良好的沟通能力。

四、教学方法1. 理论授课:通过讲授基础知识,引导学生了解数学建模的概念和步骤。

2. 实践操作:组织学生动手实践,参与实际问题的建模和求解过程。

3. 小组讨论:鼓励学生在小组内讨论并解决问题,加强团队合作和沟通能力。

4. 作业练习:布置作业练习,提供问题求解的机会,巩固学生的知识和技能。

五、教学评估1. 课堂表现:考察学生的参与度、思维逻辑和问题解决能力。

2. 作业考核:通过作业的完成情况,评估学生对知识的掌握程度。

3. 实践项目:组织学生实施实际项目,并对项目结果进行评估。

4. 小组评价:学生之间进行互评,评估团队合作和沟通效果。

六、教学资源1. 教材:提供适合教学内容的教材,包括数学建模原理和实例分析。

数学建模作业答案

数学建模作业答案

一、教材76页第1章习题1第7题(来自高中数学课本的数学探究问题,满分10分) 表1.17是某地一年中10天的白昼时间(单位:小时),请选择合适的函数模型,并进行数据拟合.一、解:根据地理常识,某地的白昼时间是以一年为周期而变化的,以日期在一年中的序号为自变量x ,以白昼时间为因变量y ,则根据表1.17的数据可知在一年(一个周期)内,随着x 的增加,y 先增后减,y 大约在6月21日(夏至)达到最大值,在12月21日(冬至)达到最小值,在3月21日(春分)或9月21日(秋分)达到中间值。

选择正弦函数sin(2/365)y A x b πϕ=++作为函数模型。

根据表1.17的数据,推测,A b ϕ和的值,作非线性拟合得26.9022sin( 1.3712)12.385365y x π=-+,预测该地12月21日的白昼时间为5.49小时。

二、教材100页第2章习题2第1题(满分10分)继续考虑第2.2节“汽车刹车距离”案例,请问“两秒准则”和“一车长度准则”一样吗?“两秒准则”是否足够安全?对于安全车距,你有没有更好的建议?二、解:“两秒准则”表明亲厚车距D 与车速v 成正比例关系2D K v =,其中22K s =。

对于小型汽车,“一车长度准则”与“两秒准则”不一致。

由221[()]d D v k v K k -=--可以计算得到当212()/54.428/v K k k km h <-=时有d D <,“两秒准则”足够安全,或者把刹车距离实测数据和“两秒准则”画在同一幅图中,根据图形指出“两秒准则”足够安全的车速范围。

用最大刹车距离除以车速,得到最大刹车距离所需要的尾随时间,并以尾随时间为依据,提出更安全的准则,如“3秒准则”,“4秒准则”或“t 秒准则”(见下图)三、教材100页第2章习题2第3题(满分10分)继续考虑第2.3节“生猪出售时机”案例,做灵敏度分析,分别考虑农场每天投入的资金对最佳出售时机和多赚的纯利润的影响.三、解:(,)2,(,)4dt c dQ x S t c S Q c dc tdcQ =∙=-=∙=-四、教材143页第3章习题3第2题(满分10分)某种山猫在较好、中等及较差的自然环境下,年平均增长率分别为 1.68%、0.55%和-4.5%. 假设开始时有100只山猫,按以下情况分别讨论山猫数量逐年变化的过程及趋势:(1) 三种自然环境下25年的变化过程,结果要列表并图示;(2) 如果每年捕获3只,山猫数量将如何变化?会灭绝吗?如果每年只捕获1只呢? (3) 在较差的自然环境下,如果要使山猫数量稳定在60只左右,每年要人工繁殖多少只?解:(1)设第k 年山猫的数量为k x ,列式得1(1)(0,1,2,k k x r x k +=+=…),用循环语句计算,并列表和作图。

数学建模实验报告1

数学建模实验报告1

数学建模实验报告1桂林电⼦科技⼤学2017-2018学年第1学期数学建模⼀、实验⽬的1. 熟悉MATLAB 软件的⽤户环境;2. 了解MATLAB 软件的⼀般命令;3. 掌握MATLAB 向量、数组、矩阵操作与运算函数;4. 掌握MATLAB 软件的基本绘图命令;5. 掌握MATLAB 语⾔的⼏种循环、条件和开关选择结构及其编程规范。

⼆、实验内容1. MATLAB 软件的矩阵输⼊和操作2. ⽤MA TLAB 语⾔编写命令M ⽂件和函数M ⽂件3. 直接使⽤MATLAB 软件进⾏作图练习;三、实验任务1. 有⼀个4×5的矩阵,编程求出其元素最⼤值及其所在的位置。

Jm.m ⽂件代码: clear;a=input('请输⼊⼀个4*5矩阵'); max=a(1,1); maxi=0; maxj=0; for i=1:4 for j=1:5if a(i,j)>max max=a(i,j); maxi=i; maxj=j;end end endfprintf('最⼤值为:%d 位置:o%d %d \n',max,maxi,maxj); 实验结果:2. 有⼀函数f(x,y)=x 2+sin xy+2y,写⼀程序,输⼊⾃变量的值,输出函数值。

Jm_5.m ⽂件代码: function f=Jm_5(x,y) f=x.^2+sin(x*y)+2*y;实验结果:3.⽤surf,mesh绘制曲⾯z=2x2+y2。

Jm5.m代码:x=-3:0.1:3;y=1:0.1:5;[X,Y]=meshgrid(x,y);Z=2*X.^2+Y.^2;subplot(1,2,1);surf(X,Y,Z);title('surf(x,y)');subplot(1,2,2);mesh(X,Y,Z);title('mesh(x,y)');实验结果:4.在同⼀平⾯的两个窗⼝中分别画出⼼形线和马鞍⾯。

数学建模技巧与实例

数学建模技巧与实例

● 02
第2章 数学建模的数学基础
数学分析
数学分析是数学的基 础学科之一,主要研 究函数的极值、曲线 拟合以及导数求解等 内容。通过数学分析, 可以深入理解数学模 型中的各种数学规律 和运算方法。
概率与统计
概率分布
了解随机变量的 分布规律
参数估计
估计总体参数的 值
统计分析
对数据进行统计 性分析
数学建模的步骤
确定问题
明确问题目标和 限制条件
模型求解
利用数学技巧解 决问题
验证和分析
检验模型的有效 性和可靠性
建立模型
选择合适的数学 工具和方法
数学建模的应用范围
经济学
市场预测、资源 配置等
工程学
结构设计优化、 流体力学模拟等
管理学
风险评估、决策 支持等
生物学
生态环境模拟、 基因分析等
数学建模的意义
季节性分析
研究数据在不同 季节中的波动情

数学建模的常用技术
数学建模技术是将数学方法应用于实际问题的过 程,包括最优化方法、随机过程、图论和时间序 列分析等内容,通过数学建模可以解决实际问题、 优化决策和预测趋势。掌握这些技术可以帮助分 析复杂问题、优化资源配置和提升决策效率。
● 04
第4章 数学建模实例分析
数学建模技巧与实例
汇报人:XX
2024年X月
目录
第1章 数学建模概述 第2章 数学建模的数学基础 第3章 数学建模的常用技术 第4章 数学建模实例分析 第5章 数学建模的工具与软件 第6章 数学建模实践与展望
● 01
第1章 数学建模概述
什么是数学建模
数学建模是一种将现 实生活中的复杂问题 转化为数学模型的方 法。通过建立数学模 型,我们可以更好地 理解问题的本质,从 而提出解决方案。数 学建模在科学研究和 工程应用中起着至关 重要的作用。

电大数学建模复习资料

电大数学建模复习资料

第一次作业(第1-2章)一、填空题1.设年利率为0.05,则10年后20万元的现值按照复利计算应为 . 2.所谓数学建模的五步建模法是指下列五个基本步骤,按一般顺序可以写出为 .3.设某种商品的需求量函数是,1200)(25)(+-=t p t Q 而供给量函数是3600)1(35)(--=t p t G ,其中)(t p 为该商品的价格函数,那麽该商品的均衡价格是 .4.在夏季博览会上,商人预测每天冰淇淋销量N 将和下列因素有关:(1) 参加展览会的人数n ;(2)气温T 超过C10;(3)冰淇淋的售价p .由此建立的冰淇淋销量的比例模型应为 .5.有人观察到鱼尾每摆动一次,鱼所移动的距离几乎与鱼身的长度相等,则鱼尾摆动的次数T (次/秒)、鱼身的长度L 和它的速度V 的关系式为 .二、分析判断题1. 从下面不太明确的叙述中确定要研究的问题,需要哪些数据资料(至少列举3个),要做些甚麽建模的具体的前期工作(至少列举3个) ,建立何种数学模型:一座高层办公楼有四部电梯,早晨上班时间非常拥挤,该如何解决.2. 一条公路交通不太拥挤,以至人们养成“冲过”马路的习惯,不愿意走临近的“斑马线”.交管部门不允许任意横穿马路,为方便行人,准备在一些特殊地点增设“斑马线”,以便让行人可以穿越马路.那末“选择设置斑马线的地点”这一问题应该考虑哪些因素?试至少列出3种.3.怎样解决下面的实际问题.包括需要哪些数据资料,要作些什么观察、试验以及建立什么样的数学模型等. (1)估计一个人体内血液的总量.(2)为保险公司制定人寿保险计划(不同年龄的人应缴纳的金额和公司赔偿的金额). (3)估计一批日光灯管的寿命.(4)确定火箭发射至最高点所需的时间. (5)决定十字路口黄灯亮的时间长度.(6)为汽车租赁公司制订车辆维修、更新和出租计划.(7)一高层办公楼有4部电梯,早晨上班时间非常拥挤,试制订合理的运行计划4.在2.5节森林救火模型中,如果考虑消防队员的灭火速度λ与开始救火时的火势b 有关,试假设一个合理的函数关系,重新求解模型.5. 假设某个数学模型建成为如下形式:.e ])1(1[)(22122x ax x M x P --= 试在适当的假设下将这个模型进行简化.三、计算题 1.在超市购物时你注意到大包装商品比小包装商品便宜这种现象了吗?比如洁银牙膏50g 装的每支1.50元,120g 装的每支3.00元,二者单位重量的价格比是1.2:1.试用比例方法构造模型解释这个现象.(1)分析商品价格c 与商品重量w 的关系.价格由生产成本、包装成本和其它成本决定,这些成本中有的与重量w 成正比,有的与表面积成正比,还有与w 无关的因素.(2)给出单位重量价格c 与w 的关系,画出它的简图,说明w 越大c 越小,但是随着w 的增加c 减小的程度变小.解释实际意义是什么?2.设某产品的售价为p ,成本为q ,售量为x (与产量相等),则总收入与总支出分别为px I =,qx C =.试在产销平衡的情况下建立最优价格模型.四、综合题某人身高1.70 m , 以适当的初速度在地球表面上可跳过与其身高相同的高度.试利用类比建模法说明:若该人以相同的初速度在月球上跳,试问他能跳多高?(地球与月球的重力加速度之比为6:1)第一次作业(第1-2章)讲评一、填空题1.解:根据现值计算公式:10)05.01(20)1(+=+=n R S Q 2783.1221201011≈=(万元) 应该填写:12.2783万元.2.解:应该填写:问题分析,模型假设,模型建立,模型求解,模型分析. 3.解: 由商品的均衡价格公式:80352536001200)(=++=++=c ad b t p 应该填写:80.4.解:因为冰淇淋的销量与人数n 、气温T 成正比,与售价p 成反比,因此应该填写: ),10(,/)10(0C T p T Kn N ≥-= K 是比例常数;5.解:因为鱼尾摆动的次数T 、鱼身的长度L 与它的速度V 成正比,因此应该填写:kTL V = (k 是常数);二、分析判断题 1.解:(1)要研究的问题:如何设置四部电梯的停靠方式,使之发挥最大效益.(2)所需资料为:每天早晨乘电梯的总人数、各层上、下电梯的人数、电梯的速度、楼层的高度、层数等. (3)要做的具体建模前期工作:观察和统计所需资料,一般讲,需要统计一周内每天的相关资料. (4)可以建立概率统计模型,亦可在适当的假设下建立确定性模型. 2.解:(1)车流的密度 (2)车的行驶速度 (3)道路的宽度(4)行人穿越马路的速度 (5)设置斑马线地点的两侧视野等. 3.解:(1)注射一定量的葡萄糖,采集一定容量的血样,测量注射前后葡萄糖含量的变化,即可估计人体的血液总量.注意采集和测量的时间要选择恰当,使血液中的葡萄糖含量充分均匀,又基本上未被人体吸收.(2)调查不同年龄的人的死亡率,并估计其在未来一定时期的变化,还应考虑银行存款利率和物价指数,保险金与赔偿金之比大体上应略高于死亡率.(3)从一批灯管中取一定容量的样本,测得其平均寿命,可作为该批灯管寿命的估计值.为衡量估计的精度,需要从样本寿命确定该批灯管寿命的概率分布,即可得到估计值的置信区间.还可试验用提高电压的办法加速寿命测试,以缩短测量时间.(4)根据牛顿第二定律建立火箭向上发射后的运动方程,初速已知,若不考虑空气阻力,很容易算出到达最高点(即速度为零)时间;若考虑空气阻力,不妨设其与火箭速度(或速度的平方)成正比,并有试验及拟合方法确定阻力系数,再解方程得到结果.(5)司机看到黄灯后停车要有一定的刹车距离S 1,设通过十字路口的距离为S 2,汽车行驶速度为v ,则黄灯的时间长度t 应使距停车线S 1之内的汽车能通过路口,即t ≈(S 1+S 2)/v .S 1可由试验得到,或按照牛顿第二定律解运动方程,进一步可考察不同车重、不同路面及司机反应灵敏程度等因素的影响.(6)根据资料和经验确定维修费用随着车龄和行驶里程的增加而增加的关系,再考虑维修和更新费用,可以以一年为一个时段,结合租金决定应该维修或更新.(7)统计在各层上班的人数,通过数据或计算确定电梯运行时间,以等待的人数与时间乘积为目标,建立优化模型,确定每部电梯运行的楼层(有的从大厅直接运行到高层). 4.解:不妨设1)(+'=b b λλ,表示火势b 越大,灭火速度λ越小,分母b +1中的1是防止b →0x 时λ→∞而加的.最优解为λβλβλ'++'+++'=)1()(21]()1(2[23221b c b b b c b c x . 5.解:当ax较小的时候,可以利用二项展开式将小括号部分简化为 ,21)1(222122ax a x -≈- 从而有2e 2)(2x x a M x P =.若x 也很小,则可以利用x x+≈1e 将其进一步化简为 ).1(2)(22x x aMx P +=三、计算题 1.在超市购物时你注意到大包装商品比小包装商品便宜这种现象了吗?比如洁银牙膏50g 装的每支1.50元,120g 装的每支3.00元,二者单位重量的价格比是1.2:1.试用比例方法构造模型解释这个现象.(1)分析商品价格c 与商品重量w 的关系.价格由生产成本、包装成本和其它成本决定,这些成本中有的与重量w 成正比,有的与表面积成正比,还有与w 无关的因素.(2)给出单位重量价格c 与w 的关系,画出它的简图,说明ww 的增加c 减小的程度变小.解释实际意义是什么? 解:(1)生产成本主要与重量w 成正比,包装成本 主要与表面积s 成正比,其它成本也包含与w 和s 成正 比的部分,上述三种成本中都含有与w 和s 无关的成分. 又因为形状一定时一般有s ∝w 2/3,故商品的价格可表为 C = αw +β w 2/3+γ(α,β,γ为大于0的常数).(2)单位重量价格131--++==w w wCc γβα,其 图2 简图如图2所示.显然c 是w 的减函数,说明大包装商品比小包装商品便宜;曲线是下凸的,说明单价的减少值随包装的变大是逐渐降低的,不要追求太大包装的商品.2.解:因为售量x 依赖于价格p ,记作)(p f x =,称为需求函数,它是p 的减函数.由此可知收入I 和支出C 都是价格p 的函数,所以利润U 可以表示为)()()(p C p I p U -= (2.8)使利润U (p )达到最大的最优价格p *可以由0d d *==p p pU 得到,即**d d d d p p p p pC pI ===(2.9)其中p I d d 称为边际收入,pC d d 称为边际支出.(2.9)式表明最大利润在边际收入等于边际支出时达到.假设需求函数是线性函数,即bp a p f -=)(,0,>b a (2.10)并且每件产品的成本q 与产量x 无关,将总收入函数、总支出函数、需求函数和(2.10)式代入(2.8)式可得))(()(bp a q p p U --=用微分法求出使U (p )达到最大的最优价格p *为baq p 22*+=(2.11) 在(2.10)式中a 可以理解为这种产品免费供应时(p = 0)社会的需求量,称为“绝对需求量”.pxb d d -=表示价格上涨一个单位时销售量下降的幅度.在实际工作中a ,b 可以由价格p 和售量x 的统计数据用最小二乘法拟合来确定.(2.11)式表明最优价格是两部分之和,一部分是成本q 的一半,另一部分与“绝对需求量”成正比,与市场需求对价格的敏感系数成反比.四、综合题解:问题分析 由于月球上的情况不了解,可先建立我们所熟悉的在地球上的有关结论,然后通过类比来加以解决.模型假设(1) 人在地球上跳高与空气阻力关系微弱,故可忽略空气阻力不计; (2) 在地面上跳高,实际上就是克服地球引力把身体“抛”到高处.其实质是把人体的重心提高到了1.70米,故可视人体为一质点.一般地,人体的重心约在身高的一半处. 模型构建与求解 依假设,可视跳高为以初速 v 0 把位于身高一半处的一质点铅直上抛.为了求出所跳高度x 与时间 t 的函数关系,可建立起跳处为原点,水平方向为 x 轴,铅直向上为 y 轴正向的平面直角坐标系.则由g tv-=d d ,0)0(v v = 知 v (t ) = -gt + v 0 (2.12) 又由)(d d t v t x =,85.0270.1)0(==x 得 85.02)(02++=t v gt t x (2.13) 类比建模: 在月球上跳高与在地球上跳高相比是完全类似的,所差的仅是重力加速度.设月球上的重力加速度为g m ,若记月球上的速度及位置函数分别为v 0,x m (因题设初始速度相同,故仍记月球上的初速度为v 0), 则应有 v m (t ) = - g m t + v 0 (2.14)85.02)(02++=t v t g t x m m (2.15) 由(2.15)知,为求出此人在月球上能跳多高,只需求出初速v 0及跳到最高处所需时间.注意到初速与地球的相同,故可由式(2.12),(2.13)求之:因跳到最高处时v 0= 0,故v 0 = gt ,于是 t = v 0/g .又此人在地球上跳了1.70m 高,故有85.0)()(2170.10020++-=gv v g vg由此得v 0=g 7.1= 4.082 m/s (2.16)于是该人在地球上跳到1.70m 高处时所用的时间为t = v 0/g = 0.42s . 以下再求在月球上以相同的初速跳到最高处所用的时间,即由 (2.14) 式及v m (t)=0,得v 0= g m t m ,即g 7.1=g m t m ,由此可得t m =g 7.1/g m . (2.17)将(2.16),(2.17)两式代入(2.15)式,便有x m ≈-21g m (m g g 7.1)2+g 7.1(mg g 7.1)+0.85 =27.1mg g+0.85=5.9 (m)即在月球上能跳过的高度约为5.9米.(m g g 6=)模型分析为求出人在月球上的活动结论,与同类活动在地球上的相应结论通过类比方法加以解决,这是类比法的又一个成功范例.同样,利用地球上的初速及相应公式求得月球上所需数据也是很关键的一步,亦是巧妙之举.第二次作业(第3-5章)一、填空题1.设开始时的人口数为0x ,时刻t 的人口数为)(t x ,若人口增长率是常数r ,那麽人口增长问题的马尔萨斯模型应为 ,其解为 .2.设某种物资有两个产地21,A A ,其产量分别为10、20,两个销地21,B B 的销量相等均为15.如果从任意产地到任意销地的单位运价都相等为a ,则最优运输方案与运价具有 两个特点.二、分析判断题1.对于技术革新的推广,在下列几种情况下分别建立模型.(1)推广工作通过已经采用新技术的人进行,推广速度与采用新技术的人数成正比,推广是无限的. (2)总人数有限,因而推广速度还会随着尚未采用新技术人数的减少而降低. (3)在(2)的前提下考虑广告等媒介的传播作用.三、计算题1.与Logistic 模型不同的另一种描述种群增长规律的是Gompertz 模型:xNrx t xln )(= ,其中r 和N 的意义与Logistic 模型相同.设渔场鱼量的自然增长服从这个模型,且单位时间捕捞量为h =Ex .讨论渔场鱼量的平衡点及其稳定性,求最大持续产量h m 及获得最大产量的捕捞强度 E m 和渔场鱼量水平x *0.2. 试求如表4所示运输问题的最优运输方案和最小运输费用:B 1、B 2、B 3供选择,运费依次为20,40和30;而进口港也有三个可供选择,代号为C 1,C 2和C 3,运费为:B 1到C 1、C 2、C 3依次为70、40、60,B 2到C 1、C 2、C 3依次为30、20、40,B 3到321,,C C C 依次为40、10、50;进口后可经由两个城市D 1, D 2运抵目的地E ,从C 1、C 2、C 3到D 1、D 2的运费为10和40,60和30,30和30;从D 1、D 2到E 的运费则为30和40. 试利用图模型协助策划一个运输路线,使总运费最低.四、综合题一个星级旅馆有150个客房,经过一段时间的经营实践,旅馆经理得到一些数据:若每间客房定价为160元,住房率为55%;每间客房定价为140元,住房率为65%;每间客房定价为120元,住房率为75%;每间客房定价为100元,住房率为85%.欲使每天收入最高,每间客房定价应为多少?注:本题要求按照五步建模法给出建模全过程.第二次作业(第3-5章)讲评一、填空题1.解 应该填写:⎪⎩⎪⎨⎧==0)0(d d x x rx t x ,.e )(0rtx t x =2.解:因为该问题从任意产地到任意销地的单位运价都相等故其具有最优运输方案不惟一;总运费均相等特点.应该填写: 最优运输方案不惟一;总运费均相等.二、分析判断题1.解:设t 时刻采用新技术的人数为x(t ). (1)指数模型xtxλ=d d . (2)Logistic 模型)(d d x N ax tx-=,N 为总人数. (3)广告等媒介在早期作用较大,它对传播速度的影响与尚未采用新技术的人数成正比,在模型(2)的基础上,有))((d d x N b ax tx-+= (2)和(3)区别见图1.图1三、计算题1.解: 模型为Ex xNrx x F x-==ln )( , 如图3所示,有2个平衡点:x = 0和x 0 =rE N -e.可证x = 0不稳定,x 0稳定(与E ,r 的大小无关).最大持续 产量为h m = rN/e ,获得h m 的E m = r ,x *0 =e /N .图32. 解:易见,这是一个产销平衡且为最小值类型的运输问题.我们有 (1) 利用最小元素法可得初始方案如表5,(2)使用闭回路法可得负检验数为12λ= -1,故令12x 进基(3)使用闭回路法进行调整知11x 出基,便得新的运输方案如表6表6(4)再进行检验知,所有检验数0≥ij λ,故得最优运销图如图2:图2 最小费用为385(百元).3.某公司自国外A 厂家进口一部分精密机器.由厂家到出口港有三个港口B 1、B 2、B 3供选择,运费依次为20,40rN/A 1 B 3 B 2 5 15 A 2 B 2 B 1 10 5 A 3 B 4 B 2 10 15和30;而进口港也有三个可供选择,代号为C 1,C 2和C 3,运费为:B 1到C 1、C 2、C 3依次为70、40、60,B 2到C 1、C 2、C 3依次为30、20、40,B 3到321,,C C C 依次为40、10、50;进口后可经 由两个城市D 1, D 2运抵目的地E , 从C 1、C 2、C 3到D 1、D 2的运费 为10和40,60和30,30和30; 从D 1、D 2到E 的运费则为30和 40. 试利用图模型协助策划一个 运输路线,使总运费最低.解:首先建立图模型如图7. 图7 利用双标号法求最短路线过程如图8.图8利用逆向搜索法可得最优运输方案为方案1 ,223E D C B A ⇒⇒⇒⇒ 方案2 ,113E D C B A ⇒⇒⇒⇒方案3 .112E D C B A ⇒⇒⇒⇒ .110min =l四、综合题 解:(一)问题分析1. 易于看出,定价每降低20元,住房率便增加10%,呈线性增长趋势;2. 160元的定价是否为最高价应给予确定;3. 是否所有客房定价相同需要确定. (二) 模型假设1. 在无其他信息时,每间客房的最高定价均为160元;2. 所有客房定价相同. (三)模型建立根据假设1.,如果设y 代表旅馆一天的总收入,而x 表示与160元相比降低的房价,则可得每降低1钱元的房价,住房率增加为10%/20=0.005.由此便可以得到)005.055.0)(160(150x x y +-= (1) 注意到,1005.055.0≤+x 又得到,900≤≤x 于是得到所求的数学模型为: max )005.055.0)(160(150x x y +-=,.900≤≤x (四)模型求解这是一个二次函数的极值问题,利用导数方法易于得到]90,0[25∈=x 为唯一驻点,问题又确实存在最大值,故25=x (元)即为价格降低幅度,也即160-25=135(元)应为最大收入所对应的房价.(五)模型分析 1. 将房价定在135元时,相应的住房率为%,5.6725005.055.0=⨯+最大收入为75.13668%5.67135150max =⨯⨯=y (元).表面上住房率没有达到最高,但是总收入达到最大,这自然是住房率与价格相互制约造成.2. 可以将五种定价的总收入求出以做比较(从略)和检验,知我们的结果是正确的.3. 为了便于管理,将价格定在140元/(天.间)也无妨,因为此时的总收入与最高收入仅差18.75元.4. 假如定价是180元,住房率应为45%,其相应的收入只有12150元,由此可知,我们的假设1.是正确的.13春综合练习题1一、填空题 1.设开始时的人口数为0x ,时刻t 的人口数为)(t x ,若人口增长率是常数r ,那麽人口增长问题的马尔萨斯模型应为 ,其解为 .2.设某种物资有两个产地21,A A ,其产量分别为10、20,两个销地21,B B 的销量相等均为15.如果从任意产地到任意销地的单位运价都相等为a ,则最优运输方案与运价具有 两个特点.二、分析判断题1.对于技术革新的推广,在下列几种情况下分别建立模型.(1)推广工作通过已经采用新技术的人进行,推广速度与采用新技术的人数成正比,推广是无限的. (2)总人数有限,因而推广速度还会随着尚未采用新技术人数的减少而降低. (3)在(2)的前提下考虑广告等媒介的传播作用.三、计算题1. 试求如表4所示运输问题的最优运输方案和最小运输费用:B 1、B 2、B 3供选择,运费依次为20,40和30;而进口港也有三个可供选择,代号为C 1,C 2和C 3,运费为:B 1到C 1、C 2、C 3依次为70、40、60,B 2到C 1、C 2、C 3依次为30、20、40,B 3到321,,C C C 依次为40、10、50;进口后可经由两个城市D 1, D 2运抵目的地E ,从C 1、C 2、C 3到D 1、D 2的运费为10和40,60和30,30和30;从D 1、D 2到E 的运费则为30和40. 试利用图模型协助策划一个运输路线,使总运费最低.一、填空题1.解 应该填写:⎪⎩⎪⎨⎧==0)0(d d x x rx t x ,.e )(0rtx t x =2.解:因为该问题从任意产地到任意销地的单位运价都相等故其具有最优运输方案不惟一;总运费均相等特点. 应该填写: 最优运输方案不惟一;总运费均相等.二、分析判断题1.解:设t 时刻采用新技术的人数为x (t ). (1)指数模型x txλ=d d . (2)Logistic 模型)(d d x N ax tx-=,N 为总人数. (3)广告等媒介在早期作用较大,它对传播速度的影响与尚未采用新技术的人数成正比,在模型(2)的基础上,有))((d d x N b ax tx-+= (2)和(3)区别见图1.图1三、计算题1. 解:易见,这是一个产销平衡且为最小值类型的运输问题.我们有 (1) 利用最小元素法可得初始方案如表5,表5(2)使用闭回路法可得负检验数为12λ= -1,故令12x 进基(3)使用闭回路法进行调整知11x 出基,便得新的运输方案如表6表6(4)再进行检验知,所有检验数0≥ij λ,故得最优运销图如图2:图2 最小费用为385(百元).2.某公司自国外A 厂家进口一部分精密机器.由厂家到出口港有三个港口B 1、B 2、B 3供选择,运费依次为20,40和30;而进口港也有三个可供选择,代号为C 1,C 2和C 3,运费为:B 1到C 1、C 2、C 3依次为70、40、60,B 2到C 1、C 2、C 3依次为30、20、40,B 3到321,,C C C 依次为40、10、50;进口后可经由两个城市D 1, D 2运抵目的地E ,从C 1、C 2、C 3到D 1、D 2的运费 为10和40,60和30,30和30; 从D 1、D 2到E 的运费则为30和40. 试利用图模型协助策划一个运输路线,使总运费最低. 解:首先建立图模型如图7. 利用双标号法求最短路线过程如图8.图8利用逆向搜索法可得最优运输方案为方案1 ,223E D C B A ⇒⇒⇒⇒ 方案2 ,113E D C B A ⇒⇒⇒⇒方案3 .112E D C B A ⇒⇒⇒⇒ .110min =l13春综合练习题2一、填空题1.设年利率为0.05,则10年后20万元的现值按照复利计算应为 .A 1B 3 B 2 5 15 A 2 B 2 B 1 10 5 A 3 B 4 B 2 10 152.所谓数学建模的五步建模法是指下列五个基本步骤,按一般顺序可以写出为 .3.设某种商品的需求量函数是,1200)(25)(+-=t p t Q 而供给量函数是3600)1(35)(--=t p t G ,其中)(t p 为该商品的价格函数,那麽该商品的均衡价格是 .4.在夏季博览会上,商人预测每天冰淇淋销量N 将和下列因素有关:(2) 参加展览会的人数n ;(2)气温T 超过C10;(3)冰淇淋的售价p .由此建立的冰淇淋销量的比例模型应为 .5.有人观察到鱼尾每摆动一次,鱼所移动的距离几乎与鱼身的长度相等,则鱼尾摆动的次数T (次/秒)、鱼身的长度L 和它的速度V 的关系式为 .二、分析判断题1. 从下面不太明确的叙述中确定要研究的问题,需要哪些数据资料(至少列举3个),要做些甚麽建模的具体的前期工作(至少列举3个) ,建立何种数学模型:一座高层办公楼有四部电梯,早晨上班时间非常拥挤,该如何解决.2. 一条公路交通不太拥挤,以至人们养成“冲过”马路的习惯,不愿意走临近的“斑马线”.交管部门不允许任意横穿马路,为方便行人,准备在一些特殊地点增设“斑马线”,以便让行人可以穿越马路.那末“选择设置斑马线的地点”这一问题应该考虑哪些因素?试至少列出3种.3.怎样解决下面的实际问题.包括需要哪些数据资料,要作些什么观察、试验以及建立什么样的数学模型等. (1)估计一个人体内血液的总量.(2)为保险公司制定人寿保险计划(不同年龄的人应缴纳的金额和公司赔偿的金额). (3)估计一批日光灯管的寿命.(4)确定火箭发射至最高点所需的时间. (5)决定十字路口黄灯亮的时间长度.(6)为汽车租赁公司制订车辆维修、更新和出租计划.(7)一高层办公楼有4部电梯,早晨上班时间非常拥挤,试制订合理的运行计划三、计算题 1.在超市购物时你注意到大包装商品比小包装商品便宜这种现象了吗?比如洁银牙膏50g 装的每支1.50元,120g 装的每支3.00元,二者单位重量的价格比是1.2:1.试用比例方法构造模型解释这个现象.(1)分析商品价格c 与商品重量w 的关系.价格由生产成本、包装成本和其它成本决定,这些成本中有的与重量w 成正比,有的与表面积成正比,还有与w 无关的因素.(2)给出单位重量价格c 与w 的关系,画出它的简图,说明w 越大c 越小,但是随着w 的增加c 减小的程度变小.解释实际意义是什么?第一次作业(第1-2章)讲评一、填空题1.解:根据现值计算公式:10)05.01(20)1(+=+=n R S Q 2783.1221201011≈=(万元)应该填写:12.2783万元.2.解:应该填写:问题分析,模型假设,模型建立,模型求解,模型分析. 3.解: 由商品的均衡价格公式:80352536001200)(=++=++=c ad b t p 应该填写:80.4.解:因为冰淇淋的销量与人数n 、气温T 成正比,与售价p 成反比,因此应该填写:),10(,/)10(0C T p T Kn N ≥-= K 是比例常数;5.解:因为鱼尾摆动的次数T 、鱼身的长度L 与它的速度V 成正比,因此应该填写:kTL V = (k 是常数);二、分析判断题 1.解:(1)要研究的问题:如何设置四部电梯的停靠方式,使之发挥最大效益.(2)所需资料为:每天早晨乘电梯的总人数、各层上、下电梯的人数、电梯的速度、楼层的高度、层数等. (3)要做的具体建模前期工作:观察和统计所需资料,一般讲,需要统计一周内每天的相关资料. (4)可以建立概率统计模型,亦可在适当的假设下建立确定性模型. 2.解:(1)车流的密度 (2)车的行驶速度 (3)道路的宽度(4)行人穿越马路的速度 (5)设置斑马线地点的两侧视野等. 3.解:(1)注射一定量的葡萄糖,采集一定容量的血样,测量注射前后葡萄糖含量的变化,即可估计人体的血液总量.注意采集和测量的时间要选择恰当,使血液中的葡萄糖含量充分均匀,又基本上未被人体吸收.(2)调查不同年龄的人的死亡率,并估计其在未来一定时期的变化,还应考虑银行存款利率和物价指数,保险金与赔偿金之比大体上应略高于死亡率.(3)从一批灯管中取一定容量的样本,测得其平均寿命,可作为该批灯管寿命的估计值.为衡量估计的精度,需要从样本寿命确定该批灯管寿命的概率分布,即可得到估计值的置信区间.还可试验用提高电压的办法加速寿命测试,以缩短测量时间.(4)根据牛顿第二定律建立火箭向上发射后的运动方程,初速已知,若不考虑空气阻力,很容易算出到达最高点(即速度为零)时间;若考虑空气阻力,不妨设其与火箭速度(或速度的平方)成正比,并有试验及拟合方法确定阻力系数,再解方程得到结果.(5)司机看到黄灯后停车要有一定的刹车距离S 1,设通过十字路口的距离为S 2,汽车行驶速度为v ,则黄灯的时间长度t 应使距停车线S 1之内的汽车能通过路口,即t ≈(S 1+S 2)/v .S 1可由试验得到,或按照牛顿第二定律解运动方程,进一步可考察不同车重、不同路面及司机反应灵敏程度等因素的影响.(6)根据资料和经验确定维修费用随着车龄和行驶里程的增加而增加的关系,再考虑维修和更新费用,可以以一年为一个时段,结合租金决定应该维修或更新.(7)统计在各层上班的人数,通过数据或计算确定电梯运行时间,以等待的人数与时间乘积为目标,建立优化模型,确定每部电梯运行的楼层(有的从大厅直接运行到高层).三、计算题 1.在超市购物时你注意到大包装商品比小包装商品便宜这种现象了吗?比如洁银牙膏50g 装的每支1.50元,120g 装的每支3.00元,二者单位重量的价格比是1.2:1.试用比例方法构造模型解释这个现象.(1)分析商品价格c 与商品重量w 的关系.价格由生产成本、包装成本和其它成本决定,这些成本中有的与重量w 成正比,有的与表面积成正比,还有与w 无关的因素.(2)给出单位重量价格c 与w 的关系,画出它的简图,说明w 越大c 越小,但是随着w 的增加c 减小的程度变小.解释实际意义是什么? 解:(1)生产成本主要与重量w 成正比,包装成本 主要与表面积s 成正比,其它成本也包含与w 和s 成正 比的部分,上述三种成本中都含有与w 和s 无关的成分. 又因为形状一定时一般有s ∝w 2/3,故商品的价格可表为 C = αw +β w 2/3+γ(α,β,γ为大于0的常数).(2)单位重量价格131--++==w w wCc γβα,其 图2 简图如图2所示.显然c 是w 的减函数,说明大包装商品比小包装商品便宜;曲线是下凸的,说明单价的减少值随包装的变大是逐渐降低的,不要追求太大包装的商品.数学建模13春综合练习题3一、填空题。

数学建模实验(传染病模型)

数学建模实验(传染病模型)

实验二:传染病模型1、SI 模型的建立基于以下三个假设,求出平衡点,给出参数,图示模型曲线。

(1)不考虑人口的出生、死亡、流动等种群动力因素。

人口始终保持一个常数,即()K t N ≡。

(2)一个病人一旦与易感者接触就必然具有一定的传染力。

假设t 时刻单位时间内,一个病人能传染的易感者数目与此环境内易感者总数()t S 成正比,比例系数为β,从而在t 时刻单位时间内被所有病人传染的人数为()()t I t S β。

2、SIS 模型的建立基于以下三个假设,求出平衡点,给出参数,图示模型曲线。

(1)不考虑人口的出生、死亡、流动等种群动力因素。

人口始终保持一个常数。

即()K t N ≡。

(2)一个病人一旦与易感者接触就必然具有一定的传染力。

假设t 时刻单位时间内,一个病人能传染的易感者数目与此环境内易感者总数()t S 成正比,比例系数为β,从而在t 时刻单位时间内被所有病人传染的人数为()()t I t S β。

(3)t 时刻,单位时间内从染病者中治愈的人与病人数量成正比,比例系数为γ,单位时间内治愈的人不具有免疫,将再成为易感者。

3、SIR 模型的建立基于以下三个假设,求出平衡点、给出参数、图示模型曲线。

(1)不考虑人口的出生、死亡、流动等种群动力因素。

人口始终保持一个常数,即()K t N ≡。

(2)一个病人一旦与易感者接触就必然具有一定的传染力。

假设t 时刻单位时间内,一个病人能传染的易感者数目与此环境内易感者总数()t S 成正比,比例系数为β,从而在t 时刻单位时间内被所有病人传染的人数为()()t I t S β。

(3)t 时刻,单位时间内从传染者中移出的人数与病人数量成正比,比例系数为γ,单位时间内移出者的数量为γ)(t I 。

求解过程1、SI 模型:由题目条件假设可以得到微分方程:K()()dIK S t I t dtβ=,又因为()()1S t I t +=, 令初始时刻病人的比例为0I ,则有:0()(1()),(0)dII t I t I I dtβ=-= %求平衡点,r 为有效传染率,x 病人比例 syms r xsolve('r*x*(1-x)','x') ans = 0 1 %方程求解syms i r t dsolve('Di=r*i*(1-i)','i(0)=i0','t')ans =1/(1-exp(-r*t)*(-1+i0)/i0) %绘制图形r=0.5,i0=0.01 fplot('1/(1-exp(-r*t)*(-1+i0)/i0)',[0,40]) fplot('1/(1-exp(-0.5*t)*(-1+0.01)/0.01)',[0,40]) function di=isf(t,i)di=0.5*i*(1-i); [t,i]=ode45(@isf,[0 40],[0.01]);plot(t,i)t ♓i♎♓ ♎♦图示4 SI 模型的i~t 曲线 图示5 SI 模型的di/dt~i 曲线2、SIS 模型 根据SI 模型及增加的假设条件,可得:)()()(t KI t I t KS dtdiKγβ-=,即: 0)0(),())(1)((I I t I t I t I dtdi=--=γβ 记 γβσ=, 则方程改写为 )]1([σβ---=i i i dt di%求解方程syms r b i t % b 为有效传染率,r 为治愈率dsolve('Di=b*i*(1-i)-r*i','i(0)=i0','t')ans =(b-r)/(b-exp(-(b-r)*t)*(-b+r+i0*b)/i0/(b-r)*b+exp(-(b-r)*t)*(-b+r +i0*b)/i0/(b-r)*r)%求平衡点syms x %(b=0.5,r=0.2)solve('0.5*x*(1-x)-0.2*x; ')ans =0..60000000000000000000000000000000%绘制图形function di=sisf(t,i)di=0.5*i*(1-i)-0.2*i;[t,i]=ode45(@sisf,[0 40],[0.01]);plot(t,i)t♓t ♓图示6 SIS 模型的i~t 曲线(σ>1) 图示7 SIS 模型的i~t 曲线(σ≤1)fplot('-0.5*x*[x-(1-1/20)]',[0,1]) fplot('-0.5*x*[x-(1-2)]',[ 0,1])i♎♓ ♎♦i♎♓ ♎♦图示8SIS 模型的di/dt~i 曲线(σ>1) 图示9SIS 模型的di/dt~i 曲线(σ≤1) 3、 SIR 模型模型的方程为{00()()(),(0)()(),(0)dIS t I t I t I I dtdSS t I t S S dtβγβ=-==-=function dx=sirf(t,x)dx=zeros(2,1);dx(1)=0.5*x(1)*x(2)-0.2*x(1); %x(1)表示i,x(2)表示s dx(2)=-0.5*x(1)*x(2);[t,x]=ode45(@sirf,[0 50],[0.01 0.99]);plot(t,x(:,1),t,x(:,2)),grid,pauseplot(x(:,2),x(:,1)),grid00.20.40.60.81s图示10 SIR模型的图形)(),(tStI图示11 SIR模型的相轨线备注:由于Matlab与Word连接不好,所绘制的图形上标的字符在Word中看不清楚。

2022年Python数学实验与建模第1章 线性规划

2022年Python数学实验与建模第1章  线性规划

数学建模算法与应用
第1章 线性规划
5.灵敏度分析
灵敏度分析是指对系统因周围条件变化显示出 来的敏感程度的分析。
实际问题中aij ,bi ,c j是怎么来的?是确定的常数吗?
航空基础学院数学第教12研页室
数学建模算法与应用
第1章 线性规划
两个问题 (1)如果参数aij ,bi ,c j 中的一个或者几个发生了 变化,现行最优方案会有什么变化?
第三步:根据问题的目标,构造关于决策变量的 一个线性函数,即为目标函数。
航空基础学院数学第教6研页室
数学建模算法与应用
第1章 线性规划
3.线性规划模型的形式
线性规划模型的一般形式(代数形式)为
max(或min)z c1 x1 c2 x2 cn xn ,
a11 x1 a12 x2
s.t.
第1章 线性规划
约束条件——①生产甲、乙两种机床所花费的加工 时间不能超过 A、B、C 机器每天的最大可用加工时间, 故
2 x1 x2 10, x1 x2 8, x2 7.
约束条件——②甲乙两种机床的产量还应该满足非 负约束,即
xi 0, i 1, 2。
航空基础学院数学第教3研页室
数学建模算法与应用
求解数学规划模型的常用软件有 MATLAB、Python、LINGO 等
MATLAB 求解数学规划问题采用两种模式: (1)基于求解器的求解方法 (2)基于问题的求解方法
航空基础学院数学第教14研页室
数学建模算法与应用
第1章 线性规划
Naval Aeronautical University
02 线性规划模型求解及 应用
航空基础学院数学第教15研页室
数学建模算法与应用

数学建模线性方程组的数值解法

数学建模线性方程组的数值解法
或 AX=b
直接法 经过有限次算术运算求出精确解(实际上 由 于 有 舍 入 误 差 只 能 得 到 近 似 解 ) ---- 高 斯 (Gauss)消元法及与它密切相关的矩阵LU分解 迭代法 从初始解出发,根据设计好的步骤用逐次 求出的近似解逼近精确解 ---- 雅可比(Jacobi) 迭代法和高斯—塞德尔(Gauss—Seidel)迭代法
(k )
0.1x1
( k 1)
0.3x2
( k 1)
1.4
Gauss-Seideil迭代公式 Dx ( k 1) Lx ( k 1) Ux ( k ) b
用它作除数会导致舍入误 差的很大增加 解决 办法 选
(k ) aik
(i k , n) 最大的一个(列主元)
将列主元所在行与第k行交换后, 再按上面的 高斯消元法进行下去,称为列主元消元法。
直接法 - 高斯消元法的矩阵表示
高斯消元法的第一次消元
a11 x1 a12 x2 a1n xn b1 a21 x1 a22 x2 a2 n xn b2 an1 x1 an 2 x2 ann xn bn
数值解法(迭代解法)的收敛性
实验5的主要内容
1. 两类数值解法: 直接方法;迭代方法
2. 超定线性方程组的最小二乘解 3. 线性方程组数值解法的MATLAB实现 4. 实际问题中方程组的数值解
线性方程组的一般形式、两类解法
a11 x1 a12 x 2 a1n x n b1 a 21 x1 a 22 x 2 a 2 n x n b2 a n1 x1 a n 2 x 2 a nn x n bn
大学数学实验
Mathematical Experiments 实验5 线性代数方程组的数值解法

数学建模第1章线性规划

数学建模第1章线性规划

数学
建模
例 1.6
min{max
xi
yi
|
ei
|},其中e i
=
xi -
yi 。
取v
=
max yi
|
e
i
|,这样,上面的问题就变换成
min v,
s.t.
ìïïíïïî
x1 y1
-
y1 ? x1 ?
v,L , xn v,L , yn
yn ? v, n ? v.
25/39
基础部数学教研室
数学 建模
2x1 - 5x2 + x3 ? 10, x1 + 3x2 + x3 ? 12, x1, x2 , x3 ³ 0.
11/39
基础部数学教研室
数学 建模
解 (1)化成 Matlab 标准型
min w = - 2x1 - 3x2 + 5x3,
s.t.
轾 犏- 2 犏 臌1
5 3
-1 1
轾 犏x1 犏 犏x2 犏 臌x3
a=1 -1 -1 1 1 -1 1 -3 1 -1 -2 3;
enddata
min=@sum(col:c*@abs(x));
@for(row(i):@sum(col(j):a(i,j)*x(j))<b(i));
@for(col:@free(x)); !x的分量可正可负;
end
24/39
基础部数学教研室
@for(row(i):@sum(col(j):a(i,j)*x(j))<b(i));
@sum(col:x)=7;
14/39
end
基础部数学教研室
数学 建模
例 1.2 求解下列线性规划问题 max z = 2x1 + 3x2 - 5x3, s.t. x1 + x2 + x3 = 7, 2x1 - 5x2 + x3 ? 10, x1 + 3x2 + x3 ? 12, x1, x2 , x3 ³ 0.

数学建模-系统可靠性分析

数学建模-系统可靠性分析

2021/3/18
6
1.3 可靠性内涵
(1)可靠性按学科分类: 一般可分为:可靠性数学;可靠性工程;可靠性管理;可 靠性物理等。
(2)可靠性的技术基础: 概率论和数理统计;材料、结构、物理学;故障物理学; 基础试验技术;环境技术等。
(3)可靠性学科特点: 可靠性学科特点是:管理与技术高度结合;众多学科的综 合;反馈和循环(通过反馈与循环不断提高产品的可靠性)。
f(t)d F (t)F '(t); 或 F (t)tf(x)d x
d t
0
2021/3/18
14
假设n(t)表示t时刻失效的产品数,△n(t)表示在(t, t+△t)时间内失效的产品数。
累 积 失 效 概 率 为 : F ˆ(t)= 到 t时 试 刻 验 失 产 效 品 的 总 产 数 品 数 = n N (t)
R (t) 1 F (t) R (t) R (t)
t 0
t
lim P (t T t t) t 0 P (T t) t
lim F (t t) F (t)
t 0
R (t) t
F '(t) R (t)
f (t) R (t)
2021/3/18
17
系列关系式:
R(t)1F(t)
失 效 率 : (t)F '(t)F '(t)f(t)R '(t)
100
99.90% 漏导致爆炸,直到2000
1000
99.01% 年12月完全关闭,14年
1万
90.48%
里乌克兰共有336万人遭
10万
36.79%
到核辐射侵害。
201210/30/1万8
<0.1%

二阶环节的阶跃响应及时间参数的影响

二阶环节的阶跃响应及时间参数的影响

实验1-2二阶环节的阶跃响应及时间参数的影响一、实验目的:1. 学会建立典型的二阶系统数学模型与传递函数。

2、加深对系统瞬态误差与稳态误差等概念的理解。

3、研究二阶系统的特征参数,阻尼比ζ和无阻尼自然频率ωn对系统动态性能的影响。

定量分析ζ和ωn 与最大超调量Mp和调节时间tS之间的关系。

4. 掌握二阶系统时域性能指标的测量方法。

二、实验原理:1、二阶系统的数学建模二阶系统是由一个比例环节和两个惯性所构成,参考电路图如下:2、二阶系统单位反馈方块图为:3、二阶系统单位反馈传递函数为:ω为系统固有频率。

其中:K为开环增益、k为闭环增益、ξ为系统阻尼比、n4、不同阻尼比下二阶系统的单位阶跃响应曲线三、实验方法与步骤:实验方法同上,在参数设置对话框中设置目的电压U1=1000mV。

1.先做二阶系统的开环时域响应,观察其曲线的变化。

2.将二阶开环系统进行单位反馈,组成二阶闭环系统,观察闭环响应曲线(注意:单位反馈的接连)。

3.在模拟机上,自行设计二阶系统、观察开环系统时域响应曲线;4.采用单位反馈构成闭环系统。

5. 计算该二阶系统模拟电路的阻尼比ζ与固有频率ωn 。

6. 改变ζ记录两种不同ζ下,二阶系统的单位阶跃响应曲线 ,将理论值与实测值进行比较。

7. 讨论典型二阶系统性能指标p M 、s t 与ζ,ωn 的关系。

四、实验内容1. 计算出该二阶系统模拟电路的阻尼比ζ与固有频率ωn 。

(见下表)2. 图示二组不同ζ下,二阶系统的单位阶跃响应曲线 ,将理论值与实测值进行比较。

实验曲线:(见下页)模型1:比例环节反馈电阻900K Ω,可变电容电容值1μF ,系统反馈电阻为100K Ω模型2:比例环节反馈电阻800KΩ,可变电容电容值1μF,系统反馈电阻为100KΩ3. 分析典型二阶系统性能指标p M 、s t 与ζ,的关系(比较1组数据即可)。

(1)p M 仅仅与阻尼比ζ有关。

ζ越大,p M 越小。

(2)s t 的大小与ζ和ωn 均有关。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

1.1 解:(1)注射一定量的葡萄糖,采集一定容积的血样,测量注射前后葡萄糖含量的变化,即估计出人体的血液总量。

注意采集和测量的时间要选择适当,使血液中的葡萄糖含量充分均匀,又基本上未被人体吸收。

(2)调查不同年龄的人的死亡率,并估计出其在未来一定时期的变化,还应考虑银行存款利率和物价指数,保险金与赔偿金之比大体上应略高于死亡率。

(3)从一批灯管中取一定容量的样本,测得其平均寿命,可作为该批灯管寿命的估计值。

为衡量估计的精度,需要从样本寿命确定该批灯管寿命的概率分布,即可得到估计值的置信区间。

还可试验用提高电压的办法加速寿命测试,以缩短测量时间。

(4)根据牛顿第二定律建立火箭向上发射后的运动方程,初速度已知,若不考虑空气阻力,很容易算出到达最高点的时间;若考虑空气阻力,不妨设其与火箭速度成正比,并由试验及拟合方法确定阻力系数,再解方程得到结果(5)司机看到黄灯后停车要有一定的刹车距离1S ,设通过十字路口的距离为2S ,汽车行驶的速度为v ,则黄灯的时间长度t 应使距停车线1S 之内的车能通过路口,即12()/t S S v ≈+。

1S 可由试验得到,或按照牛顿第二定律解运动方程,进一步可考察不同车重、不同路面及司机反应灵敏程度等因素的影响。

(6)根据资料和经验确定维修费用随着车龄和行驶里程的增加而增加的关系,再考虑维修和更新费用,可以以一年为一个时段,结合租金觉得应该维修或更新(7)统计在各层上班的人数,通过数据或计算确定电梯运行时间,以等待的人数与时间乘积为目标,建立优化模型,确定每部电梯运行的楼层。

1.2(1)分析商品的价格C 与商品重量W 的关系.价格由生产成本、包装成本和其它成本等决定,这些成本中有的与重量W 成正比,有的与表面积成正比,还有与W 无关的因素。

(2)给出单位重量价格C 与W 的关系。

画出它的简图,说明W 越大C 越小,但是随着W 的增加C 减小的程度变小。

解释实际意义是什么。

解:(1)生产成本主要与重量ω成正比,包装成本主要与表面积s 成正比,其他成本也包含于ω和s 成正比的部分,上述三种成本中都含有与ω,s 均无关的成分。

又因为形状一定时时一般有32ω∝s ,故商品的价格可表示为γβωαω++=32C (的常数)为大于0,,γβα。

(2)单位重量价格131--++=γωβωαc ,实际意义是:不要追求太大包装的商品。

图形如下: c0 w 1.3 解:(1)分段指数增长模型。

可以利用数值微分近似计算出人口的增长率。

可以大致分为3段:1790到1880年的平均年增长率为2.83%;1890到1960年的平均年增长率为1.53%;1970到2000年的平均年增长率为1.12%.3段模型为(1790年为t=0,1800年为t=1…)0.28310.153(10)210.112(18)32() 3.9,0,1,,10()(10),11,,18()(18),19,,21t t t x t e t x t x e t x t x e t --====== (2)阻滞增长模型。

可以适当实际增长率中的值做为固有增长率r ,取某些专家的估计400百万为最大容量m x ,以1790年的实际人口为0x 。

由公式0()1(1)mrtm x x t e x -=+-可以进行计算。

1.4 解:由于0t t =时,/2m x x =,可以有0()()1m r t t x x t e--=+,且0001ln m x x t r x -= 1.5 解:由题意知 ))(()()(t x x r t t x t t x m -=∆-∆+,其中r 为比例系数,令0)0(x x = 令0→∆t ,则有 ))(()(t x x r dt t dx m -= ,解得 rt m m e x x x t x ---=)()(0,如下图粗线所示:由上图知,当t 充分大时,它与Logistic 模型相近。

1.6 解:模型我们假定池中只有一种鱼。

对于这一种鱼其体型和形状是相似的,密度也大体上是相同的。

一、模型建立主要符号说明如下:W——鱼的重量、l——鱼的身长、d----鱼的胸围即鱼的最大周长、K1---第一种数学估计模型中的系数K2---第二种数学估计模型中的系数1,建立的第一种数据估计模型为:重量w与身长l的立方成正比,即W=K13l2,建立的第二种数据估计模型为:d l横截面积与鱼身最大周长的平方成正比,即W=K22(一)第一种数据估计模型对于同一种鱼,不访认为其整体形状是相似的,密度也大体上相同,所以重量w与身长l的立方成正比,即W=K13l,K1为比例系数。

把实际测得的数据代入W=K13l计算比例系数K1。

计算出实际测得的身长的平均值为: 36.8计算出实际测得的重量的平均值为:765.375把W=765.375,l=36.8代入W=K13l计算得:K1≈0.0153(二)第二种数据估计模型常调得较肥的鱼的垂钓者不一定认可上述模型,因为它对肥鱼和瘦鱼同等看待,如果d l,只假设鱼的模截面是相似的,则横截面积与鱼身最大周长的平方成正比,于是W=K22K2为比例系数。

d l计算比例系数K2。

把实际测得的数据代入W=K22计算出实际测得的胸围的平均值为:24.5875d l计算得:K2≈0.0344把W=765.375,d=24.5875, l=36.8,代入W=K22(三)第一种数据估计模型和第二种数据估计模型与实际情况的比较比较第一种数据估计模型和第二种数据估计模型与实际情况的差别,并计算误差。

计算结果如表1所示:表1计算出第一种数据估计模型与实际情况的误差为3.47% 计算出第二种数据估计模型与实际情况的误差为2.50%由上述计算结果知,第二种数据估计模型与实际情况的误差小于第一种数据估计模型与实际情况的误差。

所以,应选用W=K22d l 来建立数学模型。

1.7图1 图2 模型假设:1、降雨的速度和降水的强度保持不变; 2、人在雨中行走的速度是定量; 3、风速保持恒定; 4、人体视为一个长方体;5、假设产生的影响各个因素相互独立;问题一:不考虑降雨的角度影响:模型一:当不考虑降雨角度时,假设淋雨的部位时全身所有部位,因此淋雨的面积为)(2ac ab bc S ++=。

淋雨时间为:v D t /=,淋雨量为:[]v wac ab bc D v DSw twS C )(2++===求解:s m v h cm w m mc b m a m D /5,/2.0,2.05.0,5.1,1000======将上述数据代入模型一进行求解,有[][]Lvwac ab bc D v DSw twS C 444.23600001000*52.0*)2.0*5.15.0*5.1(22.0*5.0*1000)(2=++=++===问题二:考虑降雨的角度影响(迎面): 模型二:20πθ≤<,淋雨的部位为顶部和前方。

头顶部淋雨量为:)sin ()/(1θpu bc v D C =,前方的淋雨量为:))sin (()/(2v u p ab v D C +=θ,总的淋雨量为:[])sin (sin 21v u a cu v DbpC C C ++=+=θθ,从表达式可以看出,雨中行走的速度越快,淋雨量越小。

求解:雨中行走的速度越快,淋雨量越小。

所以取sm v m /5=时淋雨量最少。

当︒=0θ时:[]L v u a cu v DbpC C C 069.0)sin (sin 21=++=+=θθ 当︒=30θ时:[]L v u a cu v DbpC C C 606.0)sin (sin 21=++=+=θθ。

问题三:考虑降雨的角度影响(背面): 模型三:πθπ<<2/,雨滴将从身后落下(设απθ+=2/)(1):αsin u v ≤ 行走的速度慢于雨滴的水平运动速度。

淋在后背上的雨水量为vv u pabD C /)sin (3-=α,淋雨总量为:[]vv u a uc pbD C /)sin (cos -+=αα。

由于v 的限制,最大的行走速度为αsin u v m =,雨水量mv uc pbD C /)cos (α=求解:由模型三可知当速度越快时,淋雨量越少。

所以取sm v m /5=时淋雨量最少。

当︒=30α:(1) αsin u v ≤:;096.0/)cos (L v uc pbD C m ==α(2) αsin u v >时ααsin cos a c -小于零,是关于v 的增函数,所以速度应当最小,即速度为s m v /2=,此时[]L v u v a uc pbD C 024.0/)sin (cos =-+=αα问题四、问题五的求解:人行走的路线为直线,行走距离为d, 跑步最大速度sm u v 5m ax ==选择适当的直角坐标系,使人行走速度为:)0,0,(u v =,则行走的时间为 u l 。

2. 雨的速度不变,记为:()z y x v v v w ,,=相对速度:()z y x v v u v u w v ,,-=-=3. 人体为长方体,其前、侧、顶的面积之比为 c b a :: 单位时间内的淋雨量正比于cv b v a u v z y x ++-,从而总淋雨量正比于()()Tc v b v a u v u R z y x ++-= (行走的时间为u l )()ul k a u vx+-=()0>+=c v b vk z y已知 av L x ,,求u 为何值是()u R 最小?1.>x v()()()⎪⎪⎩⎪⎪⎨⎧>+-<-+=xx x x v u L u v a L v u L u v a L u R ;a v x >av x <;a v x >的情形(有最小值) av x <的情形(无最小值)当 av x >时,()xx v La v u ==m in R u R 取最小值才使。

当av x <时,u 尽可能大时,()u R 才会尽可能小。

<x v()()()L v a L a v u u l u R xx ++=++=其图像为下图易知无最小值. 同样有对a v x ==x v 0及情形的讨论。

结论:仅当0>>a v x 时,取xv u =可使前后不淋雨,其淋雨总量最小,其它情况下,应使u 尽可能的大,才能使淋雨量尽可能小,这比较符合人们生活的常识。

1.8 解:(1)令y x xt +=,则1)1()(=-+t f t f 。

画出)(t f 的示意图。

(2)写出甲公司利润的表达式)(x p 。

对于一定的y ,使)(x p 最大的x 的最优值应满足什么关系。

用图解法确定这个最优值。

相关文档
最新文档