《数学模型程序代码-Matlab-姜启源-第三章-简单的优化模型》由会员分享，可在线阅读，更多相关《数学模型程序代码-Matlab-姜启源-第三章-简单的优化模型（18页珍藏版）》请在金锄头文库上搜索。
1、第3章 简朴的优化模型 生猪的发售时机p65目的函数（生猪发售纯利润,元):（） = ( 8 gt )( 8 r ) t 640其中，t0为第几天发售，g为每天价格减少值（常数,元/公斤）,为每天生猪体重增长值（常数，公斤)。求t使(t)最大。.1(求解）模型求解p63(1)图解法绘制目的函数Q（t) ( 8 g ）( 80 r) 4t 60的图形(0 t 20）。其中， g=.1,r=2。从图形上可看出曲线Q(t）的最大值。(） 代数法对目的函数(t) = ( gt)（80 + t ) 4t 640用MATLAB求t使（t)最大。其中,r, g是待定参数。(先对Q(t)进行符号函数求导,对导
2、函数进行符号代数方程求解)然后将代入g=.1， r2，计算最大值时的t和（t)。规定: 编写程序绘制题(1)图形。 编程求解题(2)对照教材p3有关内容。有关的MATLAB函数见提示。 规定的程序和运营成果:程序：=0：1:30；g=.；r2;=(8-g.*t).*（8+.*t)-4.-64;plot(,)图形: 规定的程序和运营成果：程序:syms t r ;(8-g.*t).*(0+.*）-.*t-60;=diff(,);=slve(q);01；r=2；meal（q)Q=(8-g.*tm).*(80.m)4.tm-640运营成果:1.2（编程)模型解的的敏感性分析p634对1.1中(）所求
3、得的符号体现式t(，g），分别对g和r进行敏感性分析。(1) 取=0.,对t(r）在r=1.5:1:3上求r与的关系数据,绘制r与t的关系图形(见教材p65）。(2）取=，对t()在g=0.06：01：.5上求g与t的关系数据,绘制g与t的关系图形(见教材p65）。规定:分别编写(1)和(2)的程序,调试运营。 给出(1)的程序及运营成果:程序：syms gt r ;Q(8g.*t).*（0+r.*t)-4.t64;=diff(Q,t);=solve();=0.1;r=15:.1：；teval(q)；lo(,t);数值成果：图形成果: 给出（）的程序及运营成果：程序：yms g r;=(-g.
4、*).*（+r.*)4.*-64;q=diff（Q，);q=sol(q);r=2;g=006:0.01:0.5;t=eval(q);plot(g,)g;t数值成果:图形成果:.(编程）冰山运送模型求解p771按函数调用顺序。（1)每立方米水所需费用u为船速,0为冰山的初始体积。(2) 冰山运抵目的地后可获得水的体积为冰山达到目的地所需天数。(3) 第t天冰山球面半径融化速率：() 运送冰山费用为冰山达到目的地所需天数。() 船的日租金参照教材p1的表4,求不同V，下每立方米水的费用。下面是不完整的MATAB程序:fucin y=mainfn(）clc；=107 5*6 1;%冰山的初始体积，3
5、种uu=3 35 4.5 5;%船速,5种y=zero（legth(V0）, egth（uu);%初始化fori=1：legth（VV0) for j=1:lengt(u) y(,j)=Y（uu(),VV0(i）； edndy=run(1000*y)10000;四舍五入取整。取小数点后位数字%如下函数的输入输出均为标量function =(u,V0)%(1)y=(u，V0)/W(u,V0)；fnon y=W(,V)%(2) 编写该程序funcion y=r（t,u）%（)if =0&000/u =1.50()u*（1+0.4*u)*;elsif10/6 y=.2*(+0.4);els rr(k
6、不能不不小于！);%显示出错信息并退出运营ndunction y=S（u,V0)%（)T=400/u;y0；for =1:T rr=(3*0/4/）(1/); or k=1:t r=rr-r(k，u）; nd y=log0(rr);ndy=400f()/u+7.2*u*（u+)(3*y-151);untn =f（V0)(5) 编写该程序规定：编写所规定的程序。 运营。注:第一种函数为主函数,没有输入参数,可直接执行 成果与教材p表4比较。 完整的程序：funct =mifun（)VV0107 5*10610;uu= 3. 4. 5;yeros(lngth（V)， lnt（uu);for i=1
7、：en(V0) fr j=1：lenth（uu) y（i，)=(uu(），VV0(i)）； endnud（10y)10000;function y=(,V0) %(1）y(u,V0)/W(u,V0);funciony=W(u,V)%() 编写该程序T=4/u;r=0;for =1: rrrr+r(t，u);eny=.*pi/（(3*0)(p)(/)-rr);functin y=（,)%(3）f t&t1000/6/u y=156*1(-3)*u*（1.4u）*t;eleif 106u y=0.2*(1+0.4*u);else e(不能不不小于0！)edfuon y=(,V0)%(4)=40u;
8、y;for t=:T rr(3*0/4/)(13）; or1： rr=rr-r（k,u); en yy+og10(r);nd=400f(V）/u+2u*(u+)*(3y-151/u）;unctiy(V）%(5） 编写该程序iV0=5 y=4.0;ele V5*105&106&V=1 y=8.0;lse rrr( k超过取值范畴!);n 程序运营成果:附1:实验提示第.题MTLAB函数：,fplt，sym,sym,dif,sove,eval创立函数句柄符号绘制函数图函数fpot定义多种符号对象命令syms生成符号对象函数sy微分函数diff代数方程的符号求解函数solve符号体现式赋值函数eal把体现式或语句表达到一种字符串s，val(s)先把s转换回体现式或语句,再执行。附：第章 简朴的优化模型3.2 生猪的发售时机3.7 冰山运送
