实验2差分方程和数值微分.docxVIP

实验 2 差分方程和数值微分 题目 3 某湖泊每天有104 m3 的河水流入，河水中污物浓度为0.02 g / m3 ，经渠道排水后湖泊容 积保持200 积保持200 ?104 m3 g m 不变，现测定湖泊中污物浓度为 ，建立差分方程计算湖泊中1 年内逐月（按 30 天计算）下降的污物浓度，问要多长时间才能达到环保要求的浓度 0.04 g / m3 为了把这个时间缩短为一年，应将河水中污物浓度降低到多少 【模型及其求解】 0记湖泊总水量V 0 ，每天流入湖泊的河水量V ，湖泊最初污物浓度x ，目前河流的污物 in0k k ?1 rV浓度 pr0 ，假定外部环境没有任何其他的污染物进入湖泊，并且湖泊本身不会制造污染物， 流入湖泊的河水的污物浓度不变，不考虑每天任意时刻河水与湖水的混合过程，另外，假定湖很大，河流入水口与湖泊出水口相距很远（认为当天从湖泊流出的水的浓度是不变的），同时，湖水的总量任意时刻都保持不变（不会出现先蓄积再排出的现象）。基于以上假定， 则第 k 天湖水的污物浓度 x 只和第 k-1 天湖水的污物浓度 x in 0 k k ?1 r V VinVa ? 0 ? V V in V x k ?V p ? b in r 0 ?V ? 0 in V 0 x k ?1 V p ?in r 0 ? V 0 记 0 ， V ，则： x ? ax b k ?1 k0代入题目中给定的数据，在 MATLAB 中计算。源程序为： k 0 %第二章第 3 题源程序clear all V0 = 200*10^4; %参数初始值设定Vin = 10^4; K = 0; pr0 = 0.02; x0 = 0.2; a = (V0-Vin)/V0; P = [0 0.01 0.02 0.04 0.06 0.08 0.1 0.14 0.16 0.2 0.24] %河水不同的污物浓度 for j = 1:11 pr = P(j); b = Vin*pr/V0; x(1) = x0; x(2) = a*x(1) + b; for k = 2 : 720 x(k) = a*x(k-1) + b; %差分关系式 C(j,k) = x(k); %湖水污物浓度变化矩阵 end C(j,1) = x0; if(j == 3) %题目给出的河水污染物浓度下的情况 m = 1:12; end  end B(m+1) = x(30*m); B(1) = x0; figure, plot(0:12,B,b*,0:12,B,r-);%一年内按月份污物浓度变化(pr=0.02) title(图 1.一年内湖泊污染物浓度(按月)) xlabel(月份) ylabel(污物浓度(g/m^{3})) gtext(pr=0.02); figure, plot(1 : 720,C(1,:),m,1 : 720,C(2,:),b,...%pr 取不同值污染浓度曲线比较 1 : 720,C(3,:),y,1 : 720,C(4,:),r,1 : 720,C(5,:),g,... 1 : 720,C(6,:),c,1 : 720,C(7,:),k,1 : 720,C(8,:),m,... 1 : 720,C(9,:),g,1 : 720,C(10,:),r,1 : 720,C(11,:),k),... axis([0 720 0 0.25]); title(图 2.河水污物浓度不同时湖泊污染物浓度变化情况(按天)); xlabel(天) ylabel(污物浓度(g/m^{3})),gtext(pr=0); gtext(pr=0.01);gtext(pr=0.02);gtext(pr=0.04); gtext(pr=0.06); gtext(pr=0.08); gtext(pr=0.1); gtext(pr=0.14); gtext(pr=0.16); gtext(pr=0.2); gtext(pr=0.21); a = 0.9950; k = 360; xk = 0.04; x1 = 0.2; %如果一年内污染物浓度达要求求解 prX b = (1-a)*(xk-x1*a^(k-1))/(1-a^(k-1)); prX = b*V0/Vin; 问题解决 Step1： 首先，得到一年内 1～12 月份的湖水污物浓度数据，如表 1。 月份 1 2 3 4 5 6 污物浓度 0.1756 0.1539 0.1352 0.1191 0.1053 0.0934 月份 7 8 9 10 11 12 污物浓度 0.0831 0.0743 0.0667 0.0602 0.0546 0.0498 表 1. 一年内湖水污物浓度（月份） （浓度单位：g/m3） MATL