10090320凌威峰2.docVIP

衡阳师范学院 学 生 实 验 报 告 实验课程名称： 数学建模 实 验 内 容： 盐溶液变化规律的数学模型 系 别： 数学 年级： 10级 专业班：应用数学3班 学 生 姓 名 凌威峰 学 号 开 课 时 间： 2012 年 上 学 期 成 绩 实验指导教师 宫兆刚 1、设一容器内原有100L盐水,内含有盐10kg,现以3L/min的速度注入质量浓度为0.01kg/L的淡盐水,同时以2L/min的速度抽出混合均匀的盐水,求容器内盐量 变化的数学模型. 解：设时刻容器内的盐量为kg,考虑到时间内容器中盐的变化情况,在时间内 容器中盐的改变量注入的盐水中所含盐量－抽出的盐水中所含盐量 容器内盐的改变量为,注入的盐水中所含盐量为,时刻容器内溶液的质量浓度为,假设到时间内容器内溶液的质量浓度不变（事实上,容器内的溶液质量浓度时刻在变,由于时间很短,可以这样看）.抽出的盐水中所含盐量为,这样即可列出方程 , 即 . 又因为时,容器内有盐kg,于是得该问题的数学模型为 输入命令：x=dsolve(Dx+2*x/(100+t)-0.03 =0,x(0)=10,t) 得到结果： x=dsolve(Dx+2*x/(100+t)-0.03 =0,x(0)=10,t) x = 90000/(100+t)^2+1+1/100*t 即得解析式：X= 90000/(100+t)^2+1+1/100*t 下面对该问题进行一下简单的讨论,由上式不难发现:时刻容器内溶液的质量浓度为 P=X/(100+t)=0.01+9*10^4/(100+t)^3 且当t趋近于无穷大时,p趋近于0.01 ,即长时间地进行上述稀释过程,容器内盐水的质量浓度将趋 于注入溶液的质量浓度. 0 1 2 3 4 5 6 7 8 9 10 9.83 9.67 9.51 9.36 9.21 9.06 8.93 8.80 8.67 odefun=inline(0.03-2*x/(100+t),t,x); [t,x]=ode45(odefun,[0,150],10); [t,x] ans = 0 10.0000 2.9552 9.5203 5.9103 9.0826 8.8655 8.6825 11.8206 8.3160 16.8206 7.7630 21.8206 7.2828 26.8206 6.8640 31.8206 6.4976 36.8206 6.1759 41.8206 5.8929 46.8206 5.6433 51.8206 5.4228 56.8206 5.2278 61.8206 5.0552 66.8206 4.9022 71.8206 4.7667 76.8206 4.6468 81.8206 4.5406 86.8206 4.4469 91.8206 4.3642 96.8206 4.2915 101.8206 4.2278 106.8206 4.1723 111.8206 4.1241 116.8206 4.0826 121.8206 4.0473 126.8206 4.0176 131.8206 3.9929 136.8206 3.9729 141.8206 3.9573 146.8206 3.9455 151.8206 3.9375 156.8206 3.9327 161.8206 3.9311 166.8206 3.9324 171.8206 3.9363 176.8206 3.9427 181.8206 3.9514 186.8206 3.9622 191.8206 3.9751 193.8655 3.9808 195.9103 3.9869 197.9552 3.9933 200.0000