数学建模报告-水的流出时间

《数学建模报告-水的流出时间》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、数学建模实验报告水的流出时间软件84周天牧08161102同组完成者俞乐晨参与讨论课讨论【问题描述】三个横截面积为常数A，高分别为H1，H2，H3的水池内都盛满了水，都由池底一横截面积为B的小孔放水。求在任意时刻的水面高度和将水放空所需的时间。【模型假设】由于我们建立的是理想化的模型，做如下假设：1.题目中的水流从前一个水池流出后立即流入后一个水池化为水位的上升，在中间飞行的时间忽略不计。2.题目中水池底部的出水孔很小，无论水位多么低都可以充满以保证流出速度，除非水位为0有了这两个假设，我们可以开始建立模型。【模型建立】

2、将此问题抽象成数学问题：求初值分别为H1，H2，H3的函数h1，h2，h3随时间变化的函数，以及他们变为0所需的时间（设水池1，2，3的流出速度为v1，v2，v3）。由机械能守恒定律可知，水从第i个水池小孔流出的速度为v（i）=sqrt（2gh（i））根据任何一个水池内水量在一个时间微元内的减少量等于流出量减去流入量可以得到如下关系：水池1：-Dh1A=Ds1B将D写成d，两边同除dt可得两边积分得把水放空所需时间为水池2：-Dh2A=Ds2B-Dh1A注意到Dh1A=Ds1B,故两边写D为d并同除以dt可得化简并带入v

3、1的函数表达式可得这是一个非线性常微分方程，难以得到解析解（并非不可求，但是在此处求出解析解并不是建模的重点，因为即使h2存在解析解，h3也不一定存在。后面给出了求此方程解析解的方法。），在这里我们采用数值方法进行分析。由于第二步的数值会对后续的水池产生影响，我们选用了精度较高的Runge-Kutta法（龙格-库塔法，以下简称RK法，具体方法参见数值计算方法第六章），matlab中给出了实现的函数：[t,y]=ODE45(odefun,tspan,y0)其中t，y用于输出广义时间与相空间向量，odefun为要求的微分方程

4、，tspan为广义时间区间，y0为初值向量。定义函数：functiony=funt(t,x)A=2;B=1;g=10;H=10;y=-1B/A(sqrt(2gx)-sqrt(2gH)+B/Agt)end再调用ode45函数，即可得到h2的数值解，即一组h2关于t的向量。需要注意的是，由于非线性微分方程解的稳定性问题，此迭代法是局部收敛而非全局收敛的，故当第二个水池流空了而第一个水池还未流空时会产生复数解，需舍去。水池3：-Dh3A=Ds3B-Dh2A同理可得由于只有三步，此处我们采用简单一些的欧拉法(Euler法，具体方

5、法见数值计算方法第六章)进行数值求解。采用此方法可以得到h3的数值解，即h3关于时间的离散函数。【代码实现】在这里给出欧拉法的代码实现j=1;h3(1)=H;arrayt3(1)=0;temp3=0while(j

6、表格（前提是给定初值H1，H2，H3。其实最理想的情况其实是得到h1，h2，h3关于时间t的带参数H1，H2，H3的函数，但由于水平所限，无法用mathmatica软件实现RK方法和Euler方法，而matlab又不支持字母运算，故用数值代替之）。取H0=10，讨论以下初始情况下三个水池水位的关系：分别是三个水池初始量相等的情况，一多两少的情况和一少量多的情况 H1H2H31H0H0H02H0H0/2H0/23H0/2H0H0/24H0/2H0/2H05H0/2H0H06H0H0/2H07H0H0H0/2情况一：H1=1

7、0,H2=10,H3=10•t1=2.83•t2=3.6492•t3=4.6333情况二：H1=10,H2=5,H3=5•t1=2.83•t2=3.4458•t3=4.3738情况三：H1=5,H2=10,H3=5•t1=2•t2=3.6•t3=4.6633情况四：H1=5,H2=5,H3=10•t1=2•t2=3.4599•t3=4.4581情况五：H1=5,H2=10,H3=10•t1=2•t2=3.6•t3=4.6751情况六：H1=10,H2=5,H3=10•t1=2.83•t2=3.4458•t3=4.4707

8、情况七：H1=10,H2=10,H3=5•t1=2.83•t2=3.6492•t3=4.6233【模型讨论】由以上的建模过程可以得到这样一些结论h1是单调减函数，并且其变化规律只与其初值有关系，给定初值的条件下流空的时间是一常数。由通式-dhi/dt=(v(i)-v(i-1))B/A可知，任一时刻当第i个水池的水面高