数学建模malab第四章作业(1).pdfVIP

数学建模习题4 作业 第2 题 (1) 程序如下： a=[2.0079,0.1855,2.0079*155.79/(2.0079-0.1855)]; f=@(t)a(3).*(exp(-a(2).*t)-exp(-a(1).*t)); g=@(t)-f(t); h_1=@(t)f(t)-20; h_2=@(t)f(t)-80; [smax,fval]=fminbnd(g,0,20) x1=fzero(h_1,[0,smax]) x2=fzero(h_2,[0,smax]) x3=fzero(h_1,[smax,20]) x4=fzero(h_2,[smax,20]) 运行结果： smax = 1.3069 fval = -122.2501 x1 = 0.0689 x2 = 0.3805 x3 = 11.5887 x4 = 4.1125 (2) 程序如下： a=[2.0079,0.1855,3/4*103.86,0,0,0,0,0,0,0]; a(4)=a(3)/(a(1)-a(2)); a(5)=a(1)*a(4)/a(2); a(6)=a(3)/a(2); a(7)=a(3)*(exp(2*a(1))-1); a(8)=a(4)*exp(-2*a(1))-a(5)*exp(-2*a(2))+a(6); a(9)=a(7)/(a(1)-a(2)); a(10)=a(8)*exp(2*a(2))+a(9)*exp(2*(a(2)-a(1))); f=@(t)a(4).*exp(-a(1).*t)-a(5).*exp(-a(2).*t)+a(6); g=@(t)a(10).*exp(-a(2).*t)-a(9).*exp(-a(1).*t); f1=@(t)-f(t); g1=@(t)-g(t); [smax,fval]=fminbnd(f1,0,2) [smax,fval]=fminbnd(g1,0,20) h_1=@(t)f(t)-20; h_2=@(t)f(t)-80; h_3=@(t)g(t)-20; h_4=@(t)g(t)-80; x1=fzero(h_1,[0,2]) x2=fzero(h_2,[0,2]) x3=fzero(h_3,[2,20]) x4=fzero(h_4,[2,20]) 运行结果： smax = 2.0000 fval = -101.4297 smax = 2.6327 fval = -115.7418 x1 = 0.6233 x2 = 1.6366 x3 = 12.6196 x4 = 5.1412 程序如下： a=[2.0079,0.1855,3/4*103.86,0,0,0,0,0,0,0]; a(4)=a(3)/(a(1)-a(2)); a(5)=a(1)*a(4)/a(2); a(6)=a(3)/a(2); a(7)=a(3)*(exp(2*a(1))-1); a(8)=a(4)*exp(-2*a(1))-a(5)*exp(-2*a(2))+a(6); a(9)=a(7)/(a(1)-a(2)); a(10)=a(8)*exp(2*a(2))+a(9)*exp(2*(a(2)-a(1))); h=@(t)(2.0079*155.79/(2.0079-0.1855)).*(exp(-a(2).*t)-exp(-a(1).*t)); f=@(t)a(4).*exp(-a(1).*t)-a(5).*exp(-a(2).*t)+a(6); g=@(t)a(10).*exp(-a(2).*t)-a(9).*exp(-a(1).*t); x1=linspace(0,2,300); x2=linspace(2,20,3000);x3=linspace(0,20,3300); plot([2],[f(2)],k.,x3,h(x3),k:,x2,g(x2),k,x1,f(x1),k) legend(函数的分段点,2 小时内匀速喝三瓶啤酒,很短时间内喝三瓶啤酒) 运行结果： 140 函数的分段点 120