數(shù)學(xué)建模實(shí)驗(yàn)答案-簡單的優(yōu)化模型Word版

整理為word格式整理為word格式整理為word格式實(shí)驗(yàn)03簡單的優(yōu)化模型（2學(xué)時(shí)）（第3章簡單的優(yōu)化模型）1.生豬的出售時(shí)機(jī)p63~65目標(biāo)函數(shù)（生豬出售純利潤，元）：Q(t)=(8–gt)(80+rt)–4t–640其中，t≥0為第幾天出售，g為每天價(jià)格降低值（常數(shù)，元/公斤），r為每天生豬體重增加值（常數(shù)，公斤）。求t使Q(t)最大。1.1（求解）模型求解p63(1)圖解法繪制目標(biāo)函數(shù)Q(t)=(8–gt)(80+rt)–4t–640的圖形（0≤t≤20）。其中，g=0.1,r=2。從圖形上可看出曲線Q(t)的最大值。(2)代數(shù)法對目標(biāo)函數(shù)Q(t)=(8–gt)(80+rt)–4t–640用MATLAB求t使Q(t)最大。其中，r,g是待定參數(shù)。（先對Q(t)進(jìn)行符號函數(shù)求導(dǎo)，對導(dǎo)函數(shù)進(jìn)行符號代數(shù)方程求解）然后將代入g=0.1,r=2，計(jì)算最大值時(shí)的t和Q(t)。要求：①編寫程序繪制題(1)圖形。②編程求解題(2).③對照教材p63相關(guān)內(nèi)容。相關(guān)的MATLAB函數(shù)見提示?！镆螈俚某绦蚝瓦\(yùn)行結(jié)果：程序：t=0:1:30;g=0.1;r=2;Q=(8-g.*t).*(80+r.*t)-4.*t-640;plot(t,Q)圖形：整理為word格式整理為word格式整理為word格式★要求②的程序和運(yùn)行結(jié)果：程序：symsgtr;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運(yùn)行結(jié)果： 整理為word格式整理為word格式整理為word格式1.2（編程）模型解的的敏感性分析p63~64對1.1中(2)所求得的符號表達(dá)式t(r,g)，分別對g和r進(jìn)行敏感性分析。(1)取g=0.1，對t(r)在r=1.5:0.1:3上求r與t的關(guān)系數(shù)據(jù)，繪制r與t的關(guān)系圖形（見教材p65）。(2)取r=2，對t(g)在g=0.06:0.01:0.15上求g與t的關(guān)系數(shù)據(jù)，繪制g與t的關(guān)系圖形（見教材p65）。要求：分別編寫(1)和(2)的程序，調(diào)試運(yùn)行?！锝o出(1)的程序及運(yùn)行結(jié)果：程序：symsgtr;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]數(shù)值結(jié)果：圖形結(jié)果：整理為word格式整理為word格式整理為word格式★給出(2)的程序及運(yùn)行結(jié)果：程序：symsgtr;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(q);plot(g,t)[g;t]數(shù)值結(jié)果：整理為word格式整理為word格式整理為word格式圖形結(jié)果：2.（編程）冰山運(yùn)輸模型求解p77~81按函數(shù)調(diào)用順序。(1)每立方米水所需費(fèi)用整理為word格式整理為word格式整理為word格式u為船速，V0為冰山的初始體積。(2)冰山運(yùn)抵目的地后可獲得水的體積為冰山抵達(dá)目的地所需天數(shù)。(3)第t天冰山球面半徑融化速率：(4)運(yùn)送冰山費(fèi)用為冰山抵達(dá)目的地所需天數(shù)。(5)船的日租金參照教材p81的表4，求不同V0，u下每立方米水的費(fèi)用。下面是不完整的MATLAB程序：functiony=mainfun()clc;VV0=[10^75*10^610^6];%冰山的初始體積，3種uu=[33.544.55];%船速，5種y=zeros(length(VV0),length(uu));%初始化fori=1:length(VV0)forj=1:length(uu)y(i,j)=Y(uu(j),VV0(i));endendy=round(10000*y)/10000;%四舍五入取整。取小數(shù)點(diǎn)后4位數(shù)字整理為word格式整理為word格式整理為word格式%以下函數(shù)的輸入輸出均為標(biāo)量functiony=Y(u,V0)%(1)y=S(u,V0)/W(u,V0);functiony=W(u,V0)%(2)編寫該程序functiony=r(t,u)%(3)ift>=0&&t<=1000/6/uy=1.56*10^(-3)*u*(1+0.4*u)*t;elseift>1000/6/uy=0.2*(1+0.4*u);elseerror('k不能小于0！');%顯示出錯(cuò)信息并退出運(yùn)行endfunctiony=S(u,V0)%(4)T=400/u;y=0;fort=1:Trr=(3*V0/4/pi)^(1/3);fork=1:trr=rr-r(k,u);endy=y+log10(rr);endy=400*f(V0)/u+7.2*u*(u+6)*(3*y-151/u);functiony=f(V0)%(5)編寫該程序要求：①編寫所要求的程序。②運(yùn)行。注：第一個(gè)函數(shù)為主函數(shù)，沒有輸入?yún)?shù)，可直接執(zhí)行③結(jié)果與教材p81表4比較?！锿暾某绦颍篺unctiony=mainfun()VV0=[10^75*10^610^6];uu=[33.544.55];y=zeros(length(VV0),length(uu));fori=1:length(VV0)forj=1:length(uu)y(i,j)=Y(uu(j),VV0(i));end整理為word格式整理為word格式整理為word格式endy=round(10000*y)/10000;functiony=Y(u,V0)%(1)y=S(u,V0)/W(u,V0);functiony=W(u,V0)%(2)編寫該程序T=400/u;rr=0;fort=1:Trr=rr+r(t,u);endy=3.4*pi/3*(((3*V0)/(4*pi))^(1/3)-rr)^3;functiony=r(t,u)%(3)ift>=0&&t<=1000/6/uy=1.56*10^(-3)*u*(1+0.4*u)*t;elseift>1000/6/uy=0.2*(1+0.4*u);elseerror('k不能小于0！')endfunctiony=S(u,V0)%(4)T=400/u;y=0;fort=1:Trr=(3*V0/4/pi)^(1/3);fork=1:trr=rr-r(k,u);endy=y+log10(rr);endy=400*f(V0)/u+7.2*u*(u+6)*(3*y-151/u);functiony=f(V0)%(5)編寫該程序ifV0<=5*10^5y=4.0;elseifV0>5*10^5&&V0<=10^6y=6.2;elseifV0>10^6&&V0<=10^7y=8.0;elseerror('k超出取值范圍！');end整理為word格式整理為word格式整理為word格式★程序運(yùn)行結(jié)果：整理為word格式整理為word格式整理為word格式附1：實(shí)驗(yàn)提示第1.1題MATLAB函數(shù)：@,fplot,syms,sym,diff,solve,eval創(chuàng)建函數(shù)句柄符號 @繪制函數(shù)圖函數(shù) fplot定義多個(gè)符號對象命令 syms生成符號對象函數(shù) sym微分函數(shù) diff代數(shù)方程的符號求解函數(shù) solve符號表達(dá)式賦值函數(shù) eval把表達(dá)式或語句表示成一個(gè)字符串s，eval(s)先把s轉(zhuǎn)換回表達(dá)式或語句，再執(zhí)行。整理為word格式整理為word格式整理為word格式附2：第3章簡單的優(yōu)化模型3.2生豬的出售時(shí)機(jī)整理為word格式整理為word格式整理為word格式整理為w