数学实验盐水混合问题..docxVIP

长 沙 理 工 大 学数学与计算科学学院实 验 报 告实验项目名称 盐水混合问题 所属课程名称 数学实验 实 验 类 型 综合 实 验 日 期 2014.10.10 班 级 学 号 姓 名 成 绩 一、实验概述：【实验目的】1、学会用 Matlab 求简单微分方程的解析解. 2、学会用 Matlab 求微分方程的数值解. 3、学会建立简单的数学模型求解实际问题.【实验原理】1、用什么方法建模图形法2、用什么原理求解计算每个微小时间变化后的盐的含量与变化前盐的含量之间的关系，建立方程，然后利用MATLAB进行求解。【实验环境】硬件:计算机软件: Matlab2010 版本二、实验内容：【实验方案】1、问题描述一个圆柱形的容器，内装350L的均匀混合的盐水溶液。如果纯水以14L/s的速度从容器顶部流入，同时，容器内的混合的盐水以10.5L/s的速度从容器底部流出。开始时，容器内盐的含量为7kg。求经过时间ts后容器内盐的含量。2、问题分析考虑到盐在不同浓度的水内的扩散情况，由于容器体积为350L而纯水进入和混合盐水流出的速度都不足15L/s所以可以忽略盐在不同浓度的水内的扩散至均匀的时间。?3、符号说明假设注水时间为 tt时刻时容器内含盐量为p(t)纯水流入容器的速度为 v1混合液流出的速度为 v2。【实验过程】（实验步骤、记录、数据、分析）1、问题假设1、不考虑结晶度问题。?2、不考虑盐水混合纯水时盐在不同浓度的水内的扩散速度，即假设任意时刻容器内混合的、流出的盐水都均匀（包含初始状态350升混合盐水）。?3、水流入及盐水流出的速度均为匀速的。?2、模型建立1、溶解度=溶质/溶剂2、纯水以v1=14L/s的速度流入容器，混合液以v2=10.5L/s的速度流出容器，则根据溶解度公式，在t时刻的任一个足够小的Δt时段内：t时刻混合液中盐含量/t时刻混合液量=Δt时段混合液中盐含量的变化量/Δt时段流出的混合液量3、依题可列出方程为：p(t)/(350+v1*t-v2*t)=(p(t)-p(t+Δt))/v2*Δt当t趋向于0时，可得微分方程dp(t)/dt=-(v2*p(t)/ (350+v1*t-v2*t))即dp(t)/dt=-(10.5*p(t)/ (350+3.5*t))3、模型求解 dsolve(Dy=-(10.5*y/ (350+3.5*t)),t) ans = C2/(t + 100)^34、模型检验又由于初始t=0，时溶液中盐含量为7kg，所以可以计算出C2=7000000。5、画图程序1）由上可得方程：y=7000000/(x+100)^32）画图过程：clear allx=0:600;y=7000000./(x+100).^3;plot(x,y);xlabel(t时间)ylabel(p(t)混合中盐含量)title(混合液中盐含量曲线)【实验结论】（结果）由上步建立的数学模型和相应的编程，得到盐含量随时间变化图形如下：由下图可以得出t越大时p(t)是接近于0的，但是不等于0,即随时间增大混合液中盐含量越来越小。【实验小结】（收获体会）1、利用MATLAB可以解决一些方程的数值解求解问题。2、对于实际问题，应该根据实际情况，建立合乎实际的有效地可解的数学模型。3、MATLAB使用有待加强。三、指导教师评语及成绩：评 语评语等级优良中及格不及格1.实验报告按时完成,字迹清楚,文字叙述流畅,逻辑性强2.实验方案设计合理3.实验过程（实验步骤详细,记录完整,数据合理,分析透彻）4实验结论正确. 成 绩： 指导教师签名： 批阅日期：附录1：源 程 序1）求微分方程dp(t)/dt=-(10.5*p(t)/ (350+3.5*t))的通解：dsolve(Dy=-(10.5*y/ (350+3.5*t)),t)ans =C2/(t + 100)^32）画图过程：clear allx=0:600;y=7000000./(x+100).^3;plot(x,y);xlabel(t时间)ylabel(p(t)混合中盐含量)title(混合液中盐含量曲线)1