2013建模作业_6.docVIP

数学建模短学期作业6 1．设一容积为（单位：）的大湖受到某种化学废料的污染，污染物均匀地分布在湖中。若某时刻起污染源被切断，设湖水更新的速率是（单位是：／天）。试建立求污染物的浓度下降至原来的5%所需时间的数学模型。美国密西根湖的容积为4871×（），湖水的流量为3.663959132×（/天），求污染中止后，污染物浓度下降到原来湖水污染浓度的3%所需要的时间。 解： 由题意可知，污染源被切断后，湖内污染物的含量是一定的，设切断的时刻为0时刻，设污染物随时间变化的密度为p(t)，则污染物的总量为V*p(0)。湖水更新速率为r（单位是：／天），湖水的容积是不变的，，由于减少体积浓度与容积比与原浓度与容积比相同，则有： V*[p(t +∆t)-p(t)]=r* ∆t-tt +∆tr*psds p(t)=-p(0)*e-r*tV 则原问题可改为求p(t)=0.05p(0)时的t值。 当美国密西根湖的容积为4871×（），湖水的流量为3.663959132×（/天）时，把r=3.663959132×，V=4871×，p(t)=0.03p(0)代入，则有t= 465.3120。 2．肿瘤大小生长的速率与的次方成正比，其中为形状参数，0≤≤1；而其比例系数随时间减小，减小速率又与当时的值成正比，比例系数为环境参数。设某肿瘤参数＝1，＝0.1，的初始值为2，的初始值为1，问 （1）此肿瘤生长不会超过多大？ （2）过多长时间肿瘤大小翻一倍？ （3）何时肿瘤生长速率由递增转为递减？ （4）若参数＝2/3呢？ 解： V(t+∆t)-V(t)=K(t)*V(t)a*∆t K(t)-K(t+∆t)= K(t)*b*∆t dVdt=K*V^a dKdt=-1*b*K 先保存一个myfun.m文件： function dy = myfun(t,y) dy = zeros(2,1); a=1; b=0.1; dy(1)=y(2)*(y(1))^a; dy(2)=-1*b*y(2); end 然后在command window中输入： [t,y] = ode45(@myfun,[0:0.1:12],[1 2]); （1）[t,Y] = ode45(@myfun,[0:1:1000],[1 2]) 由上图中数据可以看出，此肿瘤生长不会超过4.8567×108。 （2）[t,y] = ode45(@myfun,[0:0.01:2],[1 2]) 由上图中数据可以看出，37天时肿瘤大小翻一倍。 （3）一阶导数图： [t,Y] = ode45(@myfun,[0:1:100],[1 2]); a=diff(Y(:,1))./diff(t); t(100)=[]; plot(t,a); 由上图中数据可以看出，在29天肿瘤生长速率由递增转为递减。 （4）[t,Y] = ode45(@myfun,[0:1:1000],[1 2]) [t,y] = ode45(@myfun,[0:0.01:2],[1 2]) [t,Y] = ode45(@myfun,[0:1:100],[1 2]); a=diff(Y(:,1))./diff(t); t(100)=[]; plot(t,a); 由上图中数据可以看出，a=2/3时，此肿瘤生长不会超过450.7963，41天时肿瘤大小翻一倍，在9天肿瘤生长速率由递增转为递减。 3．第一次世界大战中，因为战争很少捕杀鲨鱼，按理战后应能捕到很多的鲨鱼才是。可是世界大战后，在地中海却捕不到鲨鱼，因而渔民们大惑不解。 令为鱼饵的数量，是鲨鱼的数量，为时间。微分方程为 　　　　　　　　＝（－） 　　　　　　　　＝－（－） 式中、、、都是正常数。第一式中鱼饵的增长速度大体上与成正比，即按速率增加，而被鲨鱼吃掉的部分按的速率减少；第二式中鲨鱼的增长速度由于生存竞争的自然死亡互相咬食按的速率减少，但又根据鱼饵的量的变化按的速率增加。对＝3，＝2，＝2.5，＝1，＝＝1求解。画出解的曲线图观察鲨鱼和鱼饵数量的变化。 解： 先保存一个myfun.m文件： function dx = myfun(t,x) dx = zeros(2,1); a1=3; a2=2.5; b1=2; b2=1; dx(1)=x(1)*(a1-b1*x(2)); dx(2)=-1*x(2)*(a2-b2*x(1)); end 然后在command window中输入： [t,x] = ode45(@myfun,[0 10],[1 1]); plot(t,x(:,1),b); hold on; plot(t,x(:,2),g); 图中，蓝线代表鱼饵的数量，绿线代表鲨鱼的数量。由图像可知，当鲨鱼数量猛增时，鱼饵的数量下降，鱼饵供应量不足，又导致鲨鱼数量急剧下降，鲨