数学建模malab第四章作业.docVIP

数学建模习题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小时内匀速喝三瓶啤酒,很短时间内喝三瓶啤酒) 运行结果： 3. 继续考虑3.4.2小节的“酵母培养物的增长”案例，建立微分方程模型，模拟酵母培养物的增长。 解：为了刻画酵母培养物的阻滞增长，在微分方程 中，假设在时刻t 酵母培养物增长率为相应的酵母培养物的数量的线性递减函数，即，其中参数，称为酵母培养物固有增长率，参数，称为酵母培养物最大容量，另外，称为酵母培养物数量尚未实现部分，是酵母培养物尚未实现部分占最大容量的比例，表明，酵母培养物与人口尚未实现部分成正比。将模型假设代入中，得到阻滞增长方程，，用分离变量法容易得到满足初始条件的解为 % b(1) =，b(2)=， b(3) =, 程序如下： t=0:18; x=[9.6,