常微分方程研究混合溶液问题.PDF

 实 验 报 告 书 实验项目：用常微分方程研究混合溶液问题 姓名:周弘斌居文瑜 学号:50907190315090719005  班级:F0907102F0907101  学院(系):理学院数学系  所属课程:常微分方程        一、 实验目的 学会对实际问题建立常微分方程的数学模型，并用MATLAB方法进行讨论。 二、 实验内容 实际问题： 设一容器内原有100L盐，内含有盐10kg，现以3L/min的速度注入质量浓度为 0.01kg/L的淡盐水，同时以2L/min的速度抽出混合均匀的盐水，求容器内炎凉变化的 数学模型。 三、 实验步骤 设t时刻容器内的盐量为x(t)kg，考虑t到t+dt时间内容器中盐的变化情况，在dt 时间内：容器中演的改变量=注入的盐水中所含盐量-抽出的盐水中所含盐量 容器内盐的改变量为dx，注入的盐水中所含盐量为0.01*3dt，t时刻容器内溶液 的质量浓度为 ，假设t到t+dt时间内容器内溶液的质量浓度不变（此假设 参考注意事项2），于是抽出的盐水中所含盐量为 ，这样既可列出方 程 ， 即 。 又因为t=0时，容器内有盐10kg，于是得到该问题的数学模型为 这是一阶非齐次线性方程的初值问题，由MATLAB作图，可以得到其积分曲线并进 行进一步讨论。 四、 实验作图 根据以上一阶非齐次线性方程的初值问题，利用MATLAB如下： - 2 - 因此得出x(t)的图像如下： 2 y=1+1/100 t+90000/(100+t) 14 12 10 y 8 6 4 0 500 1000 1500 t 因此得出p(t)的图像如下： - 3 -