原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验 03 简单的优化模型（ 2 学时）
（第3 章 简单的优化模型）
1. 生猪的出售时机 p63~65
目标函数（生猪出售纯利润，元） ：
Q(t) = ( 8 – gt )( 80 + rt ) – 4t –640
其中， t ≥0 为第几天出售， g 为每天价格降低值（常数，元 / 公斤），r 为每
天生猪体重增加值（常数，公斤） 。
求 t 使 Q(t)最大。
1.1 （求解）模型求解 p63
(1) 图解法
绘制目标函数
Q(t) = ( 8 – gt )( 80 + rt ) – 4t –640
的图形（ 0 ≤t ≤ 20）。其中， g=0.1, r=2 。
从图形上可看出曲线 Q(t)的最大值。
(2) 代数法
对目标函数
Q(t) = ( 8 – gt )( 80 + rt ) – 4t –640
用 MATLAB 求 t 使 Q(t)最大。其中， r, g 是待定参数。 （先对Q(t)进行符号函数
求导，对导函数进行符号代数方程求解）
然后将代入 g=0.1, r=2，计算最大值时的 t 和 Q(t) 。
要求：
① 编写程序绘制题 (1) 图形。
② 编程求解题 (2).
③对照教材 p63 相关内容。
相关的 MATLAB 函数见提示。
★ 要求①的程序和运行结果：
程序：
t=0:1:30;
g=0.1;r=2;
Q=(8-g.*t).*(80+r.*t)-4.*t-640;
plot(t,Q)
图形：
1
★ 要求②的程序和运行结果：
程序：
syms g t r ;
Q=(8-g.*t).*(80+r.*t)-4.*t-640;
q=diff(Q,t);
q=solve(q);
g=0.1;r=2;
tm=eval(q)
Q=(8-g.*tm).*(80+r.*tm)-4.*tm-640
运行结果：
2
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) 的程序及运行结果：
程序：
syms g t r ;
Q=(8-g.*t).*(80+r.*t)-4.*t-640;
q=diff(Q,t);
q=solve(q);
g=0.1;r=1.5:0.1:3;
t=eval(q);
plot(r,t)
[r;t]
数值结果：
图形结果：
3
★ 给出 (2) 的程序及运行结果：
程序：
syms g t r;
Q=(8-g.*t).*(80+r.*t)-4.*t-640;
q=diff(Q,t);
q=solve(q);
r=2;g=0.06:0.01:0.15;
t=eval
文档评论（0）