matlab第二次实验报告.docVIP

《数学建模实验》 上机实验报告 班级：计算机35班 学号： 2130505099 姓名： 田博文 【实验一】 　　一上底面半径2米、下底面半径4米、高4米的圆台形水池内盛满了水，由池底一横截面积为0.001平方米的小孔放水。求在任意时刻的水面高度和将水放空所需的时间。 一、问题分析： 在dt很小一段时间内，水的流速v可以看成是不变的与水面高度h有关的函数， 即v=； 在dt时间内流出的水量，可以近似的为一圆柱， 对应的圆柱高为dh； 由于容器下降的水量与流出的水量相等，令小孔面积为b， 即dV=πdh=bvdt； 可得dh与dt的关系为：dh=bvdt/(π)； 根据相似三角形原理可得， r=4-h/2； 初始条件有； h(0)=4； 为表示方便，引入变量L=4-h 故r=2+L/2 代码： T=0; for L=0:0.001:4-0.001 %选取积分微元为0.001，即dh V=0.001*pi*(2+(1/2*L))^2; %水面处的dV，b=0.001 t=V/(sqrt(2*(4-L)*9.8)*0.001) %dV与dt关系 T=T+t; plot(T,(4-L),red); hold on end 运行结果： 结论： 　　在任意时刻的水面高度和将水放空所需的时间如上图红线所示，横轴为时间，纵轴为水面高度。