当前位置:文档之家› 数学模型程序代码-Matlab-姜启源-第三章-简单的优化模型

数学模型程序代码-Matlab-姜启源-第三章-简单的优化模型

数学模型程序代码-Matlab-姜启源-第三章-简单的优化模型
数学模型程序代码-Matlab-姜启源-第三章-简单的优化模型

第3章简单的优化模型

1. 生猪的出售时机p63~65

目标函数(生猪出售纯利润,元):

Q(t) = ( 8 – g t )( 80 + rt ) – 4t–640

其中,t≥0为第几天出售,g为每天价格降低值(常数,元/公斤),r为每天生猪体重增加值(常数,公斤)。

求t使Q(t)最大。

1.1(求解)模型求解p63

(1) 图解法

绘制目标函数

Q(t) = ( 8 – g t )( 80 + rt ) – 4t–640

的图形(0 ≤t≤ 20)。其中,g=0.1, r=2。

从图形上可看出曲线Q(t)的最大值。

(2) 代数法

对目标函数

Q(t) = ( 8 – g t )( 80 + rt ) – 4t–640

用MATLAB求t使Q(t)最大。其中,r, g是待定参数。(先对Q(t)进行符号函数求导,对导函数进行符号代数方程求解)

然后将代入g=0.1, r=2,计算最大值时的t和Q(t)。

要求:

①编写程序绘制题(1)图形。

②编程求解题(2).

③对照教材p63相关内容。

相关的MATLAB函数见提示。

★要求①的程序和运行结果:程序:

图形:

★要求②的程序和运行结果:程序:

运行结果:

1.2(编程)模型解的的敏感性分析p63~64

对1.1中(2)所求得的符号表达式t(r,g),分别对g和r进行敏感性分析。

(1) 取g=0.1,对t(r)在r=1.5:0.1:3上求r与t的关系数据,绘制r与t的关系图形(见教材p65)。

(2) 取r=2,对t(g)在g=0.06:0.01:0.15上求g与t的关系数据,绘制g与t 的关系图形(见教材p65)。

要求:分别编写(1)和(2)的程序,调试运行。

★给出(1)的程序及运行结果:

程序:

数值结果:

图形结果:

★给出(2)的程序及运行结果:程序:

数值结果:

图形结果:

2.(编程)冰山运输模型求解p77~81

按函数调用顺序。 (1) 每立方米水所需费用

)

,()

,(),(000V u W V u S V u Y =

u 为船速,V 0为冰山的初始体积。

(2) 冰山运抵目的地后可获得水的体积

3

013.4(,)(,)3T

t W u V r t u π=?=???

∑ 400

T u

=

为冰山抵达目的地所需天数。 (3) 第t 天冰山球面半径融化速率:

3

100015610(104)06(,)10000.2(10.4),6.u .u t,t u

r t u u t u -??+≤≤??=??+>

??

(4) 运送冰山费用

0011400()151(,)7.2(6)3lg (,)T t t k f V S u V u u r k u u u ==???=++- ??? ????

∑∑ 400

T u

=

为冰山抵达目的地所需天数。 (5) 船的日租金

??

?

??≤<≤

06655001010,0.8100105,2.6105,0.4)(V V V V f

参照教材p81的表4,求不同V 0,u 下每立方米水的费用。 下面是不完整的MATLAB 程序:

要求:

①编写所要求的程序。

②运行。注:第一个函数为主函数,没有输入参数,可直接执行

③结果与教材p81表4比较。

★完整的程序:

★程序运行结果:

附1:实验提示

第1.1题

MATLAB函数:@,fplot,syms,sym,diff,solve,eval

附2:第3章简单的优化模型3.2 生猪的出售时机

3.7 冰山运输

数学模型第四版(姜启源)作业对于6.4节蛛网模型讨论下列问题:

对于6.4节蛛网模型讨论下列问题: (1)因为一个时段上市的商品不能立即售完,其数量也会影响到下一时段的价格,所以第k+1时段的价格1+k y 由第k+1和第k 时段的数 量1+k x 和k x 决定。如果设1+k x 仍只取决于k y ,给出稳定平衡的条件,并 与6.4的结果进行比较。 (2)若除了1+k y 由1+k x 和k x 决定之外,1+k x 也由前两个时段的价格k y 和 1-k y 决定,试分析稳定平衡的条件是否还会放宽。 解:(1) 设1+k y 由1+k x 和k x 的平均值决定,即价格函数表示为: )2 (11k k k x x f y +=++ 则 0),2 (0101>-+-=-++ααx x x y y k k k 0),(001>-=-+ββy y x x k k 消去y, 得到 012)1(22x x x x k k k +=++++αβαβαβ ,k=1,2,…. 该方程的特征方程为 022=++αβαβλλ 与6.4节中 )2 (11-++=k k k y y g x 时的特征方程一样, 所以0<αβ<2, 即为0p 点的稳定条件。

(2)设 )2 (11k k k x x f y +=++ )2 (11-++=k k k y y g x , 则有 0),2 (0101>-+-=-++ααx x x y y k k k 0),2 (0101>-+=--+ββy y y x x k k k 消去y,得到 0123)1(424x x x x x k k k k +=++++++αβαβαβαβ 该方程的特征方程为 02423=+++αβαβλαβλλ 令λ=x ,αβ=a , 即求解三次方程 0a 2ax ax 4x 23=+++ 的根 在matlab 中输入以下代码求解方程的根x : syms x a solve(4*x^3+a*x^2+2*a*x+a==0,x) 解得 1x = (36*a^2 - 216*a - a^3 + 24*3^(1/2)*(-a^2*(a - 27))^(1/2))^(1/3)/12 - a/12 + (a*(a - 24))/(12*(36*a^2 - 216*a - a^3 + 24*3^(1/2)*(-a^2*(a - 27))^(1/2))^(1/3)); 2x = -(2*a*(36*a^2 - 216*a - a^3 + 24*3^(1/2)*(-a^2*(a - 27))^(1/2))^(1/3) - 3^(1/2)*a*24*i - 3^(1/2)*(36*a^2 - 216*a - a^3 + 24*3^(1/2)*(-a^2*(a - 27))^(1/2))^(2/3)*i - 24*a + 3^(1/2)*a^2*i +

数学模型第四版课后答案姜启源版

《数学模型》作业答案 第二章(1)(2012年12月21日) 1. 学校共1000名学生,235人住在A 宿舍,333人住在B 宿舍,432人住在C 宿舍.学生们 要组织一个10人的委员会,试用下列办法分配各宿舍的委员数: (1). 按比例分配取整数的名额后,剩下的名额按惯例分给小数部分较大者; (2). §1中的Q 值方法; (3).d ’Hondt 方法:将A 、B 、C 各宿舍的人数用正整数n=1,2,3,……相除,其商数如下表: 将所得商数从大到小取前10个(10为席位数),在数字下标以横线,表中A 、B 、C 行有横线的数分别为2,3,5,这就是3个宿舍分配的席位.你能解释这种方法的道理吗? 如果委员会从10个人增至15人,用以上3种方法再分配名额,将3种方法两次分配的结果列表比较. 解:先考虑N=10的分配方案, ,432 ,333 ,235321===p p p ∑==3 1 .1000i i p 方法一(按比例分配) ,35.23 1 11== ∑=i i p N p q ,33.33 1 22== ∑=i i p N p q 32.43 1 33== ∑=i i p N p q 分配结果为: 4 ,3 ,3321===n n n 方法二(Q 值方法) 9个席位的分配结果(可用按比例分配)为: 4 ,3 ,2321===n n n

第10个席位:计算Q 值为 ,17.92043223521=?=Q ,75.92404333322=?=Q 2.9331544322 3=?=Q 3Q 最大,第10个席位应给C.分配结果为 5 ,3 ,2321===n n n 方法三(d ’Hondt 方法) 此方法的分配结果为:5 ,3 ,2321===n n n 此方法的道理是:记i p 和i n 为各宿舍的人数和席位(i=1,2,3代表A 、B 、C 宿舍). i i n p 是每席位代表的人数,取,,2,1Λ=i n 从而得到的i i n p 中选较大者,可使对所有的,i i i n p 尽量接近. 再考虑15=N 的分配方案,类似地可得名额分配结果.现将3种方法两次分配的结果列表如下: 2. 试用微积分方法,建立录像带记数器读数n 与转过时间的数学模型. 解: 设录像带记数器读数为n 时,录像带转过时间为t.其模型的假设见课本. 考虑t 到t t ?+时间内录像带缠绕在右轮盘上的长度,可得,2)(kdn wkn r vdt π+=两边积分,得 ?? +=n t dn wkn r k vdt 0 )(2π )22 2 n wk k(r n πvt +=∴ .2 22n v k w n v rk t ππ+=∴ 《数学模型》作业解答 第三章1(2008年10月14日)

姜启源《数学模型》第三版课件

第一章建立数学模型1.1 从现实对象到数学模型1.2 数学建模的重要意义1.3 数学建模示例 1.4 数学建模的方法和步骤1.5 数学模型的特点和分类1.6 怎样学习数学建模

1.1从现实对象到数学模型 我们常见的模型 玩具、照片、飞机、火箭模型… …~ 实物模型水箱中的舰艇、风洞中的飞机… …~ 物理模型地图、电路图、分子结构图… …~ 符号模型 模型是为了一定目的,对客观事物的一部分 进行简缩、抽象、提炼出来的原型的替代物 模型集中反映了原型中人们需要的那一部分特征

你碰到过的数学模型——“航行问题” 用x 表示船速,y 表示水速,列出方程: 75050)(750 30)(=?-=?+y x y x 答:船速每小时20千米/小时. 甲乙两地相距750千米,船从甲到乙顺水航行需30小时,从乙到甲逆水航行需50小时,问船的速度是多少? x =20y =5求解

航行问题建立数学模型的基本步骤?作出简化假设(船速、水速为常数); ?用符号表示有关量(x, y表示船速和水速); ?用物理定律(匀速运动的距离等于速度乘以时间)列出数学式子(二元一次方程); ?求解得到数学解答(x=20, y=5); ?回答原问题(船速每小时20千米/小时)。

数学模型(Mathematical Model) 和 数学建模(Mathematical Modeling) 对于一个现实对象,为了一个特定目的, 根据其内在规律,作出必要的简化假设, 运用适当的数学工具,得到的一个数学结构。建立数学模型的全过程(包括表述、求解、解释、检验等)数学模型 数学 建模

第四版姜启源数学模型复习总结(2015年春)

第四版姜启源数学模型复习总结(2015年春) 【内容总结与思考】 第1章:了解模型的概念与分类,熟练掌握数学模型的定义,数学模型的重要应用,建模的重要例子-指数模型,Logist模型。建模的一般方法及其在建模中的应用。建模的一般步骤(每步的主要内容与问题)。建模的全过程(框图)4个环节的含义。模型的特点(技艺性)。模型分类(表现特征),建模中的能力培养。 数学建模实例的建模思想及其步骤 §1 数学模型的概念: 模型:模型是为了一定目的,对客观事物的一部分信息进行简缩、抽象、提炼出来的原型的替代物。 模型的分类:具体模型(或物质模型,实的),包括直观模型,物理模型。抽象模型(或理想模型,虚的),包括思维模型,符号模型,数学模型。 数学模型:对于一个现实对象,为了一个特定目的,根据其内在规律,作出必要的简化假设,运用适当的数学工具,得到的一个数学结构。 1-1-1 模型是为了特定的目的,将原型的()而得到的原型替代物。 1-1-2数学模型可以描述为:对于一个现实对象,( )。

1-1-3 关于数学模型的如下论述中正确的是() A。数学模型是以现实世界的特定问题为研究对象。 B。数学模型只是对实际问题的近似表示,其中包含一些简化假设。C。数学模型表示是某一特定问题的内在规律的数学表示,是以方程和函数关系表示的数学结构。 D。数学模型是现实问题的真实的描述,不能做任何假设和简化。 1-1-4 关于数学建模的如下论述中正确的是() A。数学模型和数学建模是完全相同的概念。 B。数学建模是一个全过程,包括表述、求解、解释和验证四个环节。C。数学建模全过程涉及两个世界是现实世界和虚拟世界,涉及的“双向翻译”是同声翻译和文献翻译。 D.数学建模过程是一个从理论-实践-再理论-再实践不断改进的过程。 §2 建模的重要意义 (1)数学以空前的广度和深度向一切领域渗透 在一般工程技术领域数学建模仍然大有用武之地;在高新技术领域数学建模几乎是必不可少的工具了;数学进入一些新领域,为数学建模开辟了许多处女地. 数学建模的具体应用:分析与设计,预测与决策,优化与控制,规划与管理。 例1-2-1 数学建模的具体应用为()。§3实例1:椅子问题:实际问题转换为数学问题的方法:位

数学模型姜启源第四版答案

数学模型姜启源第四版答案 【篇一:姜启源数学模型课后答案(3版)】 t>第二章(1)(2008年9月16日) 1.学校共1000名学生,235人住在a宿舍,333人住在b宿舍,432人住在c宿舍.学生们 要组织一个10人的委员会,试用下列办法分配各宿舍的委员数:(1). 按比例分配取整数的名额后,剩下的名额按惯例分给小数部分 较大者; (2). 1中的q值方法; (3).d’hondt方法:将a、b、c各宿舍的人数用正整数n=1,2,3,??相除,其商数如下表: 将所得商数从大到小取前10个(10为席位数),在数字下标以横线,表中a、b、c行有横线的数分别为2,3,5,这就是3个宿舍 分配的席位.你能解释这种方法的道理吗? 如果委员会从10个人增至15人,用以上3种方法再分配名额,将 3种方法两次分配的结果列表比较. 解:先考虑n=10的分配方案, 3 p1?235,p2?333,p3?432, ?pi?1000. i?1 方法一(按比例分配) q1? p1n 3 ?2.35,q2? p2n 3 ?3.33, q3? p3n 3 ?4.32 ? i?1 pi ? i?1 pi

i?1 pi 分配结果为: n1?3, n2?3, n3?4 方法二(q值方法) 9个席位的分配结果(可用按比例分配)为: n1?2,n2?3, n3?4 第10个席位:计算q值为 q1? 235 2 2?3 ?9204.17, q2? 333 2 3?4 ?9240.75, q3? 432 2 4?5 ?9331.2 q3最大,第10个席位应给c.分配结果为 n1?2,n2?3,n3?5 方法三(d’hondt方法) 此方法的分配结果为:n1?2,n2?3,n3?5 此方法的道理是:记pi和ni为各宿舍的人数和席位(i=1,2,3代表a、b、c宿舍). pini pini pini 是 每席位代表的人数,取ni?1,2,?,从而得到的近. 中选较大者,可使对所有的i,尽量接 再考虑n?15的分配方案,类似地可得名额分配结果.现将3种方法两次分配的结果列表如下: 2.试用微积分方法,建立录像带记数器读数n与转过时间的数学模型. 解:设录像带记数器读数为n时,录像带转过时间为t.其模型的假设见课本. 考虑t到t??t时间内录像带缠绕在右轮盘上的长度,可得 vdt?(r?wkn)2?kdn,两边积分,得 ?vdt?2?k?(r?wkn)dn t

相关主题
文本预览
相关文档 最新文档