数学实验 课程设计

1、.安徽工业大学大学数学实验课程设计姓名：班级：任课老师：数学实验 课程设计问题提出：某容器盛满水后，低端直径为的小孔开启（图）。根据水力学知识，当水面高度h时，水冲小孔中流出的速度 (g为重力加速度，0.6为孔口的收缩系数)。若容器为倒圆锥形（如图1），现测得容器高和上底面直径均为1.2m，小孔直径为3cm，问水从小孔中流完需要多长时间；2min水面高度是多少。若容器为倒葫芦形（如图2），现测得容器高为1.2m，小孔直径为3cm，有低端（记作x=0）向上每隔0.1m测出容器的直径D（m）如表所示，问水从小孔中流完需要多少时间；2min时水面的高度是多少。图1 ： 图2：表：x/m00.10.2

2、0.30.40.50.60.70.80.91.01.11.2D/m0.030.050.080.140.190.330.450.680.980.101.201.131.00问题分析：（1） 倒圆锥形容器流水问题中随时间t液面高度h也在变化，同时水的流速也在变化，再写变化难以用普通的方程进行模拟求解，考虑建立常微分方程竟而代入数值求解。水面的直径等于液面的高度。可以建立容器中水流失的液面高度对时间t的变化率。假设t时，液面的高度h，此时水的流速流量Q为： ； 则在时间内液面下降高度为，可得到关系式：；由此可知水下降时需要的时间：；根据此关系式知道。（2） 在第二问中，考虑倒葫芦形容器时因为他的高度

3、h不同容器直径D变化没有规律可循，同第一题相比我们只知道他的一些数值，这就需要我们建立高度h和容器直径D之间的关系矩阵，然后再欧拉方程和龙格库塔方法找出时间t和液面高度之间的分量关系。由（1）可同理推知：假设在时间t时，液面高度为h，此时流量为；经过时，液面下降，若我们取的t是在t(n)和t(n+1)之间的某一时刻，于是就可在误差范围内得到；可以得到 ；建立模型：（1） 在试验中我们不考虑圆锥的缺省对流水的影响，以及其他外界因素和玻璃的毛细作用，试验中水可以顺利流完。实验中重力加速度g=9.8；倒圆锥的液面最初高度为H=1.2m，液面直径D=1.2m=0.03，小孔的直径为=0.03m；接上文

4、中分析结论代入数据：即在T时间内将1.2m的液面高度放完，（matlab不支持一些运算符号，故用matlab运算格式）dt=-(pi/4)h2*dh)/(0.6*(pi/4)*d2*sqrt(gh)=-(h1.5*dh)/(0.6*d2*sqrt(g)h是由01.2m对t积分用matlab计算上式编辑文件：a1.m，d0=0.03;g=9.8;syms ht=(h1.5)/(0.6*d02*sqrt(g);T=int(t,0,1.2);eval(T)运行结果：>> a1ans = 373.2556结果：水从倒圆锥中流完需要373.26s；2mine之后液面的高度为h1;373.26

5、-120= h12.5/(1.5*d2*(g)0.5)=153.26可知h1=(1.5*d2*sqrt(g)*153.26)(2/5);Malab计算>> g=9.8g = 9.8000>> d0=0.03d0 = 0.0300>> h1=(1.5*d02*sqrt(g)*153.26)(2/5)h1 =0.8405即2mine之后液面的高度为0.84m；上述运行结果可知：谁需要373.26s流完，2mine之后液面高度为0.84m；（2） 在与（1）同样的条件下，倒葫芦形容器的液面最初高度H=1.2m，小孔的直径为d0=0.03m，液面直径和液面高度关系如

6、表。在分析中已经讨论出t和h的关系。dt=t(k+1)-t(k)=-(pi/4)h2*dh)/( 0.6(pi/4)d2*sqrt(gh)=-(h1.5*dh)/(0.6d2*sqrt(g) k1=0.15*sqrt(g*(x(n)*d2/(-43.6359*x(n)8+213.0457*x(n)7-414.873*x(n)6+410.2075*x(n)5-218.8936*x(n)4+62.553*x(n)3-8.3215*x(n)2+0.49619*x(n)+.014892)2; k2=0.15*sqrt(g*(x(n)-h*k1)*d2/(-43.6359*(x(n)-h*k1)8+21

7、3.0457*(x(n)-h*k1)7-414.873*(x(n)-h*k1)6+410.2075*(x(n)-h*k1)5-218.8936*(x(n)-h*k1)4+62.553*(x(n)-h*k1)3-8.3215*(x(n)-h*k1)2+0.49619*(x(n)-h*k1)+.014892)2;x(n+1)=x(n)-h*(k1+k2)/2;matlab程序编写及计算：g=9.8;d=0.03;k1=0;k2=0;h=0.4;x(1)=1.2;for n=1:1000 k1=0.15*sqrt(g*(x(n)*d2/(-43.6359*x(n)8+213.0457*x(n)7-4

8、14.873*x(n)6+410.2075*x(n)5-218.8936*x(n)4+62.553*x(n)3-8.3215*x(n)2+0.49619*x(n)+.014892)2; k2=0.15*sqrt(g*(x(n)-h*k1)*d2/(-43.6359*(x(n)-h*k1)8+213.0457*(x(n)-h*k1)7-414.873*(x(n)-h*k1)6+410.2075*(x(n)-h*k1)5-218.8936*(x(n)-h*k1)4+62.553*(x(n)-h*k1)3-8.3215*(x(n)-h*k1)2+0.49619*(x(n)-h*k1)+.014892

9、)2; x(n+1)=x(n)-h*(k1+k2)/2;end x(300) t=0:h:1000*h; plot(t,x); axis(0,400,0,1.21); p=polyfit(x,t,5);T=real(polyval(p,0)运行结果：>> a3ans = 1.0278Warning: Imaginary parts of complex X and/or Y arguments ignored > In a3 at 14T = 396.5254运行结果得知：倒葫芦形容器的水大约396.53s流完，2mine后液面高度为1.03m。实验总结：在研究实际问题时，建

10、立模型找出微分方程是分析客观对象变化规律的有力工具，但不是所有问题多可以建立解析表达式（如（2）问），即使获得了解析表达式也十分复杂，因此就必须通过数值求解的方法算出微分方程的某些离散点处得近似解析，进而分析微分方程反映客观问题。指点迷津：在上述问题的（1）中仍有疑问，之前思考时建立模型是倒圆锥内剩余的水的体积变化量dV对时间dt的关系，即dV/dt=Q=vt();假设t时间内水流了VV=；dV/dt=;求次微分方程：>> h=dsolve('Dh=(0.6*sqrt(9.8*h)/(pi/4)*(1.2-h)2)','h(0)=1.2','t') h = RootOf(X2245 - 8*X2244 + (152*X2243)/5 - (288*X2242)/5 + (1296*X224)/25 - (55296*pi2 + 1102500*t2 - 201600*pi*6(1/2)*t)/(3125*pi2), X224) RootOf(X23