22一些数值方法的建模问题

《22一些数值方法的建模问题》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、2.2一些数值方法的建模问题2.2.1.排水渠道的设计问题1)问题的提出土木工程师和环境工程师在设计一条排水渠道时必须考虑渠道的宽度、深度、内壁光滑度等参数及水流的速度、流量、水深等物理量之间的关系。现在要设计一条横断面为矩形的水渠，其宽度为20米，水流量为5立方米/秒，水渠的斜度系数为0.0002,Manning粗糙系数为0.03,试确定渠道中水的深度。2)符号约定H：渠道中水的深度，单位为米B：渠道宽度为，单位为米，本题中B=20米Q：水流量，单位为立方米/秒，木题中Q=5立方米/秒S:水渠的斜度系数，无量刚量，木题'l>S=0.0002N：Manning粗糙系数

2、,是一个与水渠内壁材料的光滑度有关的无量刚量，木题中N=0.03U：水的流速，单位为米/秒3)问题分析与建立模型查阅水工学应用的有关资料，获得一个有关流速的关系式：小』_

3、NB+2H(1)而质量守衡处律告诉我们Q=UBH(2)将式(1)代入式(2),同时将S,N,B,Q用给定的数值替代,有J0.0002().035_2・(2(W)3(20+2H)3当晋(2oh)§(2o+2h)刁一5令003,于是问题归结为求方程丿(已)一°的正根问题。4)模型求解由于方=。没有求根公式，故采用数值方法求根，这里采用Newton迭代法來求根。为找到迭代初始值心，先画出的图形观察：Tn[

4、l]:=f[h_]=Sqrt[0.0002]/0.03*(20*h厂(5/3)*(20+2h厂(-2/3)-5;In[2]:=Plot[f[h],{h,0,5}]输出图形为从图中可以看到函数在1附近冇根，选取迭代初始值从=1，根的粘:度为1(T做Newton迭代，对应的程序为Clear[x,f,g];f[x_]=Sqrt[0.0002]/0.03*(20*x)"(5/3)*(20+2x厂(-2/3)-5;g[x_]=D[f[x],x];x0=l；epsl=O.000000000001;nmax=500;eps=O.000001;Do[ul=g[x0];If[Abs[u

5、l//N]eps,Print["迭代失败"]]执行程序后得输出结果为11=0.7292262368n=leps=0.270774H=0.7025912396n=2eps二0.026635H-0.702293294n-3cps=0.000297946H=0.7022932563n=4eps=3.7772x10"从

6、输岀结果说明所求的水渠深度约为0.7022932563,其误差为3.7772x102.2.2.男大学生的身高问题1)问题的提出有关统计资料表明，我国大学生男性群体的平均身高约为170cm,H.该群体中约有99.7%的人身高在150cm至190cm之间。如果将［150cm,190cm］等分成20个高度区间，试问该群体身高在每_高度区间的分布情况怎样？特别地，身高中等(165cm至175cm之间)的人占该擀体的百分比会超过60貉吗？解：2)问题分析与建立模型因为-个人的身高涉及很多因素，通常它是一个服从正态分布N@,a)的随机变量。正态分布的概率密度函数

7、x-A)2(x)H«的两个参数卩、o分別为正态分布的均值与标准差。根据题目中我国大学生男性群体的平均身高约为170cm,可选取正态分布的均值

8、i=170cm,而由“该群体中约冇99.7%的人身高在150cm至190cm之间”和正态分布N(p,B的“3c规则”，有卩-3

9、=——^-==e妣20厉将［150cm,190cm］等分成20个高度区间后，得到高度区间为[150,152],[152,154],…，[188,190J对应的分布为9(—170)2800必北=150,152,154,…,188身高在165cni至175cmZ间的人占该群体的帀分比为9(—170)2［，75——e80()dxJl6520厉(2)如上式(1)和(2)的定积分是不能用定积分基木公式方法求出的，但用计算方法屮的数值积分可以算出。1)模型求解选用数值积分中的复化梯形公式求积方法编程可以计算出求误差小于10'的式(1)的定积分值从而得出相应分布依