[文学]实验2差分方程和数值微分.docxVIP

实验 2 差分方程和数值微分 题目 3 某湖泊每天有 104 m3 的河水流入， 河水中污物浓度为 0.02g / m3 ，经渠道排水后湖泊容 积保持 200 104 m3 不变，现测定湖泊中污物浓度为 0.2g / m3 ，建立差分方程计算湖泊中 1 年内逐月（按 30 天计算）下降的污物浓度，问要多长时间才能达到环保要求的浓度 3 0.04g / m ？为了把这个时间缩短为一年，应将河水中污物浓度降低到多少？ 记湖泊总水量 V0 ，每天流入湖泊的河水量 Vin ，湖泊最初污物浓度x0 ，目前河流的污物 浓度 pr0 ，假定外部环境没有任何其他的污染物进入湖泊，并且湖泊本身不会制造污染物， 流入湖泊的河水的污物浓度不变， 不考虑每天任意时刻河水与湖水的混合过程， 另外，假定 湖很大，河流入水口与湖泊出水口相距很远（认为当天从湖泊流出的水的浓度是不变的） ， 同时，湖水的总量任意时刻都保持不变（不会出现先蓄积再排出的现象） 。基于以上假定， 则第 k 天湖水的污物浓度 xk 只和第 k-1 天湖水的污物浓度xk 1 及流过的河水污物浓度 pr 有 关，可以得到如下关系式： xk V0 Vin xk 1 Vin pr 0 V0 V0 a V0 Vin b Vin pr 0 V0， V0 ，则： xk axk 1 b 记 代入题目中给定的数据，在 MATLAB 中计算。源程序为： %第二章第 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; B(m+1) = x(30*m); B(1) = x0; end end 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