matlab生猪的出售时机,数学模型程序代码-Matlab-姜启源-第三章-简单的优化模型.doc...

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

第3章 简单的优化模型1. 生猪的出售时机p6365目标函数(生猪出售纯利润，元)Qt 8 gt 80 rt 4t 640其中，t 0为第几天出售，g为每天价格降低值(常数，元/公斤)，r为每天生猪体重增加值(常数，公斤)。求t使Qt最大。1.1(求解)模型求解p631 图解法绘制目标函数Qt 8 gt 80 rt 4t 640的图形(0 t 20)。其中， g0.1, r2。从图形上可看出曲线Qt的最大值。2 代数法对目标函数Qt 8 gt 80 rt 4t 640用MATLAB求t使Qt最大。其中，r, g是待定参数。(先对Qt进行符号函数求导，对导函数进行符号代数方程求解)然后将代入g0.1, r2，计算最大值时的t和Qt。要求 编写程序绘制题1图形。 编程求解题2.对照教材p63相关内容。相关的MATLAB函数见提示。 要求的程序和运行结果程序t0130;g0.1;r2;Q8-g.*t.*80r.*t-4.*t-640;plott,Q图形 要求的程序和运行结果程序syms g t r ;Q8-g.*t.*80r.*t-4.*t-640;qdiffQ,t;qsolveq;g0.1;r2;tmqQ8-g.*tm.*80r.*tm-4.*tm-640运行结果1.2(编程)模型解的的敏感性分析p6364对1.1中2所求得的符号表达式tr,g，分别对g和r进行敏感性分析。1 取g0.1，对tr在r1.50.13上求r与t的关系数据，绘制r与t的关系图形(见教材p65)。2 取r2，对tg在g0.060.010.15上求g与t的关系数据，绘制g与t的关系图形(见教材p65)。要求分别编写1和2的程序，调试运行。 给出1的程序及运行结果程序syms g t r ;Q8-g.*t.*80r.*t-4.*t-640;qdiffQ,t;qsolveq;g0.1;r1.50.13;tq;plotr,tr;t数值结果图形结果 给出2的程序及运行结果程序syms g t r;Q8-g.*t.*80r.*t-4.*t-640;qdiffQ,t;qsolveq;r2;g0.060.010.15;tq;plotg,tg;t数值结果图形结果2.(编程)冰山运输模型求解p7781按函数调用顺序。1 每立方米水所需费用u为船速，V0为冰山的初始体积。2 冰山运抵目的地后可获得水的体积为冰山抵达目的地所需天数。3 第t天冰山球面半径融化速率4 运送冰山费用为冰山抵达目的地所需天数。5 船的日租金参照教材p81的表4，求不同V0，u下每立方米水的费用。下面是不完整的MATLAB程序function ymainfunclc;VV0107 5*106 106;冰山的初始体积，3种uu3 3.5 4 4.5 5;船速，5种yzeroslengthVV0, lengthuu;初始化for i1lengthVV0 for j1lengthuu yi,jYuuj,VV0i; endendyround10000*y/10000;四舍五入取整。取小数点后4位数字以下函数的输入输出均为标量function yYu,V0 1ySu,V0/Wu,V0;function yWu,V02 编写该程序function yrt,u3if t0t1000/6/u y1.56*10-3*u*10.4*u*t;elseif t1000/6/u y0.2*10.4*u;else errork不能小于0;显示出错信息并退出运行endfunction ySu,V04T400/u;y0;for t1T rr3*V0/4/pi1/3; for k1t rrrr-rk,u; end yylog10rr;endy400*fV0/u7.2*u*u6*3*y-151/u;function yfV05 编写该程序要求 编写所要求的程序。 运行。注第一个函数为主函数，没有输入参数，可直接执行 结果与教材p81表4比较。 完整的程序function ymainfunVV0107 5*106 106;uu3 3.5 4 4.5 5;yzeroslengthVV0, lengthuu;for i1lengthVV0 for j1lengthuu yi,jYuuj,VV0i; endendyround10000*y/10000;function yYu,V0 1ySu,V0/Wu,V0;function yWu,V02 编写该程序T400/u;rr0;for t1T rrrrrt,u;endy3.4*pi/3*3*V0/4*pi1/3-rr3;function yrt,u3if t0t1000/6/u y1.56*10-3*u*10.4*u*t;elseif t1000/6/u y0.2*10.4*u;else errork不能小于0endfunction ySu,V04T400/u;y0;for t1T rr3*V0/4/pi1/3; for k1t rrrr-rk,u; end yylog10rr;endy400*fV0/u7.2*u*u6*3*y-151/u;function yfV05 编写该程序if V05*105 y4.0;elseif V05*105V0106 y6.2;elseif V0106V0107 y8.0;else error k超出取值范围;end 程序运行结果附1实验提示第1.1题MATLAB函数,fplot,syms,sym,diff,solve,创建函数句柄符号绘制函数图函数fplot定义多个符号对象命令syms生成符号对象函数sym微分函数diff代数方程的符号求解函数solve符号表达式赋值函数把表达式或语句表示成一个字符串s，s先把s转换回表达式或语句，再执行。附2第3章 简单的优化模型3.2 生猪的出售时机3.7 冰山运输