流体力学实验思考题
参考答案
流体力学实验室
二○○六年
静水压强实验
1．同一静止液体内的测压管水头线是根什么线？ 测压管水头指p z +，即静水力学实验仪显示的测压管液面至基准面的垂直高度。测压管水头线指测压管液面的连线。实验直接观察可知，同一静止液面内的测压管水头线是一根水平线。
2．当0?B p 时，试根据记录数据，确定水箱内的真空区域。
0?B p ，相应容器的真空区域包括以下三个部分：
（1）过测压管2液面作一水平面，由等压面原理知，相对测压管2及水箱内的水体而言，该水平面为等压面，均为大气压强，故该平面以上由密封的水、气所占区域，均为真空区域。
（2）同理，过箱顶小不杯的液面作一水平面，测压管4中，该平面以上的水体亦为真空区域。
（3）在测压管5中，自水面向下深度某一段水柱亦为真空区域。这段高度与测压管2液面低于水箱液面的高度相等，亦与测压管4液面高于小水杯液面高度相等。
3．若再备一根直尺，试采用另外最简便的方法测定0γ。
最简单的方法，是用直尺分别测量水箱内通大气情况下，管5油水界面至水面和油水界面至油面的垂直高度h 和0h ，由式00h h w w γγ= ，从而求得0γ。
4．如测压管太细，对于测压管液面的读数将有何影响？
设被测液体为水，测压管太细，测压管液面因毛细现象而升高，造成测量误差，毛细高度由下式计算
式中，σ为表面张力系数；γ为液体容量；d 为测压管的内径；h 为毛细升高。常温的水，m N 073.0=σ，30098.0m N =γ。水与玻璃的浸润角θ很小，可以认为0.1cos =θ。于是有
d h 7.29= （h 、d 均以mm 计）
一般来说，当玻璃测压管的内径大于10mm 时，毛细影响可略而不计。另外，当水质不洁时，σ减小，毛细高度亦较净水小；当采用有机下班玻璃作测压管时，浸润角θ较大，
其h较普通玻璃管小。
如果用同一根测压管测量液体相对压差值，则毛细现象无任何影响。因为测量高、低压强时均有毛细现象，但在计算压差时，互相抵消了。
5．过C点作一水平面，相对管1、2、5及水箱中液体而言，这个水平面是不是等压面？哪一部分液体是同一等压面？
不全是等压面，它仅相对管1、2及水箱中的液体而言，这个水平面才是等压面。因为只有全部具有下列5个条件的平面才是等压面：（1）重力液体；（2）静止；（3）连通；（4）连通介质为同一均质液体；（5）同一水平面。而管5与水箱之间不符合条件（4），相对管5和水箱中的液体而言，该水平面不是水平面。
6、用该实验装置能演示变液位下的恒定水流吗？
关闭各通气阀门，开启底阀，放水片刻，可看到有空气由C进入水箱。这时阀门的出流就是变液位下的恒定水流。因为由观察可知，测压管1的液面始终与C点同高，表明作用于底阀上的总水头不变，故为恒定流动。这是由于液位的降低与空气补充使箱体表面真空度的减小处于平衡状态。医学上的点滴注射就是此原理应用的一例，医学上称这为马利奥特容器的变液位下恒定流。
文丘里实验
1．本实验中，影响文丘里管流量系数大小的因素有哪些？哪个因素最敏感？对本实验的管道而言，若因加工精度影响，误将()cm d 01.02-值取代2d 值时，本实验在最大流量下的μ值将变为多少？
由式
得
可见本实验的μ值大小与Q 、1d 、2d 、h ?有关。其中1d 、2d 影响最敏感。本实验中若文丘里管1d =1.4cm ，2d =0.71cm ，通常在切削加工中2d 比1d 测量方便，容易掌握好精度，2d 不易测量准确，从而不可避免的要引起实验误差。例如当最大流量时μ值为0.976，若2d 的误差为－0.01cm ，那么μ值将变为1.006，显然不合理。
2．为什么计算流量与实际流量不相等？
因为计算流量'Q 是在不考虑水头损失情况下，即按理想液体推导的，而实际流体存在粘性必引起阻力损失，从而减小过流能力，Q ?'Q ，即μ〈1.0。
3．文丘里流量计能否倾斜安装，为什么？
如图所示
根据流体静力学方程
得 22112
1
H h h H p p -?+?++=γγ
则 )()(222211212211γ
γγγp z H h h H p z p z p z +--?+?+++=???? ??+-???? ??
+ 由图可知 )()(4321h h h h h -+-=?
式中，1h 、2h 、3h 、4h 分别为各测压管的液面读数。
因此，无论文丘里流量计是否倾斜安装，对测压管读数都不影响。
4．文丘里管喉颈处容易产生真空，允许最大真空度为6-7mH 2O 。工程中应用文丘里管时，应检验其最大真空度是否在允许范围内。根据你的实验成果，分析本实验文丘里管喉颈处最大真空值为多少？
本实验若1d =1.4cm ，2d =0.71cm ，以管轴线高程为基准面，以水箱液面和喉颈断面分别为1-1、2-2计算断面，列能量方程得
则 ∴22.522
?-γp cmH 2O
即实验中最大流量时，文丘里管喉颈处真空度22.52?v h cmH 2O ，而由本实验实测约为
60cm。
进一步分析可知，若水箱水位高于管轴线4m左右时，实验中文丘里管喉颈处真空度可达7mH2O左右。
雷诺实验
1．为什么上临界雷诺数无实际意义，而采用下临界雷诺数作为层流和紊流的判据？实测下临界雷诺数为多少？
根据实验测定，上临界雷诺数实测值在3000-5000范围内，与操作的快慢、水箱的紊动度、外界干扰等密切相关。有关学者做了大量实验，有的得12000，有的得20000，有的甚至得到40000。实际水流中，干扰总是存在的，故上临界雷诺数为不定值，无实际意义。只有下临界雷诺数才可以作为判别流态的标准。凡水流的雷诺数小于下临界雷诺数者为层流。一般实测下临界雷诺数为2100左右。
2．雷诺实验得出的圆管流动下临界雷诺数理论值为2320，而目前一般教科书中介绍采用的下临界雷诺数是2000，原因何在？
下临界雷诺数也并非与干扰绝对无关。雷诺实验是在环境的干扰极小，实验前水箱中的水体经长时间的稳定情况下，经反复多次细心量测得出的。而后人的大量实验很难重复得出雷诺实验的准确数值，通常在2000-2300之间。因此，从工程实用出发，教科书中介绍的圆管下临界雷诺数一般是2000。
3．圆管内液流有哪几种流动状态？它们的特点如何？结合实验现象说明。
略。
5．为什么必须用一个保持溢流的恒压水箱来做雷诺实验？
略。
沿程阻力损失实验
1．本实验中，沿程阻力损失就是压差计的压差，如果管道有一定的倾角，压差计的压差是否还是沿程阻力损失？为什么？
现以倾斜等径管道上装设的水银多管压差计为例说明（图中A —A 为水平线）：
如图示O —O 为基准面，以1—1和2—2为计算断面，计算点在轴心处，设21v v =，∑=0j h ，由能量方程可得
这表明水银压差计的压差值即为沿程水头损失，且和倾角无关。
2．根据实测m 值判断本实验的流区。
f h l
g ～v lg 曲线的斜率m=1.0～1.8，即f
h 与8.10.1-v 成正比，表明流动为层流（m=1.0）、
紊流光滑区和紊流过渡区（未达阻力平方区）。
3．管道的当量粗糙度如何测得？
当量粗糙度的测量可用实验的方法测定λ及e R 的值，然后用下式求解：
（1）考尔布鲁克公式
莫迪图即是本式的图解。
（2）S ．J 公式
（3）Barr 公式
其中（3）式精度最高。在反求
d ?时，（2）式开方应取负号。也可直接由λ～
e R 关系在莫迪图上查得d
?，进而得出当量粗糙度?值。 4．实验工程中的钢管中的流动，大多为光滑紊流或紊流过渡区，而水电站泄洪洞的流动，大多数为紊流阻力平方区，其原因何在？
钢管的当量粗糙度一般为0.2mm ，常温下，s cm /01.02=ν，经济流速s cm /300，若
实用管径D=（20～100）cm ，其5106?=e R ～6103?，相应的d
?=0.0002～0.001，由莫迪图可知，流动均处在过渡区。
若需达到阻力平方区，那么相应的610=e R ～6109?，流速应达到(5～9)m/s 。这样高
速的有压管流在实际工程中非常少见。
而泄洪洞的当量粗糙度可达(1～9)mm ，洞径一般为（2～3）m ，过流速往往在(5～10)m/s 以上，其e R 大于7
10，故一般均处于阻力平方区。
局部阻力损失实验
1．结合实验结果，分析比较突扩与突缩圆管在相应条件下的局部阻力损失大小关系。
由式
及
表明影响局部阻力损失的因素是v 和21d d 。由于有 突扩：2211???? ?
?-=A A e ζ 突缩：???
? ??
-=2115.0A A s ζ 则有
当 5.021?A A
或
时，突然扩大的水头损失比相应的突然收缩的要大。在本实验最大流量Q 下，突然扩大损失较突然缩小损失约大一倍，即817.160.3/54.6==js je h h 。
21d d 接近于1时，突然扩大的水流形态接近于逐渐扩大管的流动，因而阻力损失显著减小。
2．结合流动仪演示的水力现象，分析局部阻力损失机理何在？产生突扩与突缩局部损失的主要部位在哪里？怎样减小局部阻力损损失？
流动演示仪1-7型可显示突扩、突缩、渐扩、渐缩、分流、合流、阀道、绕流等三十多种内、外流的流动图谱。据此对于局部阻力损失的机理分析如下：
从显示的图谱可见，凡流道边界突变处，形成大小不一的漩涡区。漩涡是产生损失的主要根源。由于水质点的无规则运动和激烈的紊动，相互磨擦，便消耗了部分水体的自储能量。另外，当这部分低能流体被主流的高能流体带走时，还须克服剪切流的速度梯度，经质点间的动能交换，达到流速的重新组合，这也损耗了部分能量。这样就造成了局部阻力损失。
从流动仪可见，突扩段的漩涡主要发生在突扩断面以后，而且与扩大系数有关，扩大系数越大，漩涡区也越大，损失也越大，所以产生突扩局部阻力损失的主要部位在突扩断面的后部。而突缩段的漩涡在收缩断面均有。突缩前仅在死角区有小漩涡，且强度较小，而突缩的后部产生了紊动度较大的漩涡环区。可见产生突缩水头损失的主要部位是在突缩断面后。
从以上分析可知，为了减小局部阻力损失，在设计变断面管道几何边界形状时应流线型化或昼接近流线形，以避免漩涡的形成，或使漩涡区尽可能小。如欲减小管道的局部阻力，就应减小管径比以降低突扩段的漩涡区域；或把突缩进口的直角改为圆角，以消除突缩断面后的漩涡环带，可使突缩局部阻力系数减小到原来的21～101。突然收缩实验管道使用年份长以后，实测阻力系数减小，主要原因也在这里。
3．现备有一段长度及联接方式与调节阀相同，内径与实验管道相同的直管段(见实验装置图)，如何用两点法测量阀门的局部阻力系数？
两点法是测量局部阻力系数的简便有效办法。它只需在被测流段（如阀门）前后的直管段长度大于（20～40）d 的断面处，各布置一个测压点便可。先测出整个测量段上的总水头损失21-w h ，有
式中：ji h ----分别为两测点间不干扰的各个局部阻力段的阻力损失；
jn h ----被测段的局部阻力损失；
21-w h ----两测点间的沿程阻力损失。
然后，把被测段（如阀门）换上一段长度及联接方法与被测段相同，内径与管道相同的
直管段，再测出相同流量下的总水头损失21-'w
h ，同样有 所以
孔口、管嘴出流实验
1．结合观测不同类型管嘴与孔口出流的流股特征，分析流量系数不同的原因及增大过流能力的途径。
根据实验结果可知，流股形态及流量系数如下：
圆角管嘴出流的流股呈光滑圆柱形，=μ ；
直角管嘴出流的流股呈圆柱形麻花状扭变，=μ ；
圆锥管嘴出流的流股呈光滑圆柱形，=μ ；
孔口出流的流股在出口附近有侧收缩，呈光滑圆柱形，=μ 。
影响流量系数大小的原因有：
（1）出口附近流股直径，孔口为=c d cm ，=d d c ，其余同管嘴的出口直径，1=d d c 。
（2）直角进口管嘴出流，μ大于孔口c μ，是因为前者进口段后由于分离，使流股侧收缩而引起局部真空（实际实验测得局部真空度约为16O cmH 2），产生抽吸作用从而加大过流能力。后者孔口出流流股侧面均为大气压，无抽吸力存在。
（3）直角进口管嘴的流股呈扭变，说明横向脉速大，紊动度大，这是因为在侧收缩断面附近形成漩涡这故。而圆角进口管嘴的流股为光滑圆柱形，横向脉速微弱，这是因为进口近乎流线形，不易产生漩涡之故，所以直角管嘴比圆角管嘴出流损失大，μ值小。
（4）圆锥管嘴虽亦属直角进口，但因进口直径渐小，不易产生分离，其侧收缩面面积接近出口面积（μ值以出口面积计），故侧收缩并不明显影响过流能力。另外，从流股形态看，横向脉动亦不明显，说明渐缩管对流态有稳定作用（工程或实验中，为了提高工作段水流和稳定性，往往在工作段前加一渐缩段，正是利用渐缩的这一水力特性）。能量损失小，因此其μ值与圆角管嘴相近。
从以上分析可知，为了加大管嘴的过流能力，进口形状应力求线化，只要将进口修圆，提高μ的效果就十分显著。孔口及直角管嘴的流量系数的实验值有时比经验值偏大，其主要原因亦与制作工艺上基使用上不小心将孔口、管嘴的进口棱角磨损了有关。
2．观察H d
> 0.1时，孔口出流的侧收缩率较H d <0.1时有何不同？ 当H d > 0.1时，观察知，收缩断面直径c d 增大，并接近于孔径d ，这叫不完全收缩，实验测知，μ增大，可达0.7左右。
北京交通大学海滨学院考试试题 课程名称：数学实验2010-2011第一学期出题教师：数学组适用专业: 09机械, 物流, 土木, 自动化 班级：学号：姓名： 选做题目序号： 1.一对刚出生的幼兔经过一个月可以长成成兔, 成兔再经过一个月后可以 繁殖出一对幼兔. 如果不计算兔子的死亡数, 请用Matlab程序给出在未来24个月中每个月的兔子对数。 解: 由题意每月的成兔与幼兔的数量如下表所示： 1 2 3 4 5 6 ··· 成兔0 1 1 2 3 5··· 幼兔 1 0 1 1 2 3··· 运用Matlab程序： x=zeros(1,24); x(1)=1;x(2)=1； for i=2:24 x(i+1)=x(i)+x(i-1); end x 结果为x = 1 1 2 3 5 8 13 21 3 4 5 5 89 144 233 377 610 987 1597 2584 4181 6765 1094 6 7711 2865 7 46368 2.定积分的过程可以分为分割、求和、取极限三部分, 以1 x e dx 为例, 利用
已学过的Matlab 命令, 通过作图演示计算积分的过程, 并与使用命令int() 直接积分的结果进行比较. 解：根据求积分的过程，我们先对区间[0,1]进行n 等分， 然后针对函数x e 取和，取和的形式为10 1 i n x i e e dx n ξ=≈ ∑ ? ，其中1[ ,]i i i n n ξ-?。这里取i ξ为区间的右端点，则当10n =时，1 x e dx ?可用10 101 1.805610 i i e ==∑ 来近似计算， 当10n =0时，100 100 1 01 =1.7269100 i x i e e dx =≈ ∑?，当10n =000时，10000 10000 1 1 =1.718410000 i x i e e dx =≈ ∑ ?. 示意图如下图，Matlab 命令如下： x=linspace (0,1,21); y=exp(x); y1=y(1:20); s1=sum(y1)/20 y2=y(2:21); s2=sum(y2)/20 plot(x,y); hold on for i=1:20 fill([x(i),x(i+1),x(i+1),x(i),x(i)],[0,0,y(i),y(i),0],'b') end syms k;symsum(exp(k/10)/10,k,1,10)；%n=10 symsum(exp(k/100)/100,k,1,100)；%n=100 symsum(exp(k/10000)/10000,k,1,10000)；%n=10000
工程流体力学及水力学实验报告及分析讨论
工程流体力学及水力学实验报告及分析讨论 实验一流体静力学实验 验原理 重力作用下不可压缩流体静力学基本方程 (1.1) 中： z被测点在基准面的相对位置高度； p被测点的静水压强，用相对压强表示，以下同； p0水箱中液面的表面压强； γ液体容重； h被测点的液体深度。 对装有水油（图1.2及图1.3）U型测管，应用等压面可得油的比重S0有下列关系： (1.2) 此可用仪器（不用另外尺）直接测得S0。 验分析与讨论 同一静止液体内的测管水头线是根什么线？ 测压管水头指，即静水力学实验仪显示的测管液面至基准面的垂直高度。测头线指测压管液面的连线。实验直接观察可知，同一静止液面的测压管水头线是一根。 当P B<0时，试根据记录数据，确定水箱内的真空区域。 ，相应容器的真空区域包括以下三部分：
)过测压管2液面作一水平面，由等压面原理知，相对测压管2及水箱内的水体而言，面为等压面，均为大气压强，故该平面以上由密封的水、气所占的空间区域，均为真。 )同理，过箱顶小水杯的液面作一水平面，测压管4中，该平面以上的水体亦为真空区)在测压管5中，自水面向下深度某一段水柱亦为真空区。这段高度与测压管2液面低液面的高度相等，亦与测压管4液面高于小水杯液面高度相等。 若再备一根直尺，试采用另外最简便的方法测定γ0。 最简单的方法，是用直尺分别测量水箱内通大气情况下，管5油水界面至水面和油 至油面的垂直高度h和h0，由式，从而求得γ0。 如测压管太细，对测压管液面的读数将有何影响？ 设被测液体为水，测压管太细，测压管液面因毛细现象而升高，造成测量误差，毛由下式计算 中，为表面张力系数；为液体的容量；d为测压管的内径；h为毛细升高。常温(t=20℃，=7.28dyn/mm，=0.98dyn/mm。水与玻璃的浸润角很小，可认为cosθ=1.0。于是有 单位为mm) 一般来说，当玻璃测压管的内径大于10mm时，毛细影响可略而不计。另外，当水质，减小，毛细高度亦较净水小；当采用有机玻璃作测压管时，浸润角较大，其h较普管小。 如果用同一根测压管测量液体相对压差值，则毛细现象无任何影响。因为测量高、时均有毛细现象，但在计算压差时，互相抵消了。 过C点作一水平面，相对管1、2、5及水箱中液体而言，这个水平面是不是等压面？
第一章 绪论 1-1．20℃的水2.5m 3 ，当温度升至80℃时，其体积增加多少？ [解] 温度变化前后质量守恒，即2211V V ρρ= 又20℃时，水的密度3 1/23.998m kg =ρ 80℃时，水的密度32/83.971m kg =ρ 32 1 125679.2m V V == ∴ρρ 则增加的体积为3 120679.0m V V V =-=? 1-2．当空气温度从0℃增加至20℃时，运动粘度ν增加15%，重度γ减少10%，问此时动力粘度μ增加多少（百分数）？ [解] 原原ρννρμ)1.01()15.01(-+== 原原原μρν035.1035.1== 035.0035.1=-=-原 原 原原原μμμμμμ 此时动力粘度μ增加了3.5% 1-3．有一矩形断面的宽渠道，其水流速度分布为μρ/)5.0(002.02 y hy g u -=，式中ρ、μ分别为水的密度和动力粘度，h 为水深。试求m h 5.0=时渠底（y =0）处的切应力。 [解] μρ/)(002.0y h g dy du -= )(002.0y h g dy du -==∴ρμ τ 当h =0.5m ，y =0时 )05.0(807.91000002.0-??=τ Pa 807.9= 1-4．一底面积为45×50cm 2，高为1cm 的木块，质量为5kg ，沿涂有润滑油的斜面向下作等速运动，木块运动速度u=1m/s ，油层厚1cm ，斜坡角22.620 （见图示），求油的粘度。
[解] 木块重量沿斜坡分力F 与切力T 平衡时，等速下滑 y u A T mg d d sin μθ== 001 .0145.04.062 .22sin 8.95sin ????= = δθμu A mg s Pa 1047.0?=μ 1-5．已知液体中流速沿y 方向分布如图示三种情况，试根据牛顿内摩擦定律y u d d μ τ=，定性绘出切应力沿y 方向的分布图。 [解] 1-6．为导线表面红绝缘，将导线从充满绝缘涂料的模具中拉过。已知导线直径0.9mm ，长度20mm ，涂料的粘度μ=0.02Pa ．s 。若导线以速率50m/s 拉过模具，试求所需牵拉力。（1.O1N ） [解] 2 53310024.51020108.014.3m dl A ---?=????==π y u u u u y u u y ττ= 0y ττy 0 τττ=0 y
大学数学实验 项目一 矩阵运算与方程组求解 实验1 行列式与矩阵 实验目的 掌握矩阵的输入方法. 掌握利用Mathematica (4.0以上版本) 对矩阵进行转置、加、减、数乘、相乘、乘方等运算, 并能求矩阵的逆矩阵和计算方阵的行列式. 基本命令 在Mathematica 中, 向量和矩阵是以表的形式给出的. 1. 表在形式上是用花括号括起来的若干表达式, 表达式之间用逗号隔开. 如输入 {2,4,8,16} {x,x+1,y,Sqrt[2]} 则输入了两个向量. 2. 表的生成函数 (1) 最简单的数值表生成函数Range, 其命令格式如下: Range[正整数n]—生成表{1,2,3,4,…,n }; Range[m, n]—生成表{m ,…,n }; Range[m, n, dx]—生成表{m ,…,n }, 步长为d x . (2) 通用表的生成函数Table. 例如,输入命令 Table[n^3,{n,1,20,2}] 则输出 {1,27,125,343,729,1331,2197,3375,4913,6859} 输入 Table[x*y,{x,3},{y,3}] 则输出 {{1,2,3},{2,4,6},{3,6,9}} 3. 表作为向量和矩阵 一层表在线性代数中表示向量, 二层表表示矩阵. 例如,矩阵 ??? ? ??5432 可以用数表{{2,3},{4,5}}表示. 输入 A={{2,3},{4,5}} 则输出 {{2,3},{4,5}} 命令MatrixForm[A]把矩阵A 显示成通常的矩阵形式. 例如, 输入命令: MatrixForm[A] 则输出 ??? ? ??5432 但要注意, 一般地, MatrixForm[A]代表的矩阵A 不能参与运算. 输入 B={1,3,5,7} 输出为 {1,3,5,7} 输入 MatrixForm[B] 输出为
一、实验内容 P206第六题 function f=wuyan2(c) y=[3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.41 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4 281.4] t=[0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210] f=y-c(1)/(1+c(1)/3.9-1)*exp^(-c(2)*t) c0=[1 1] c=lsqnonlin('wuyan2',c0) P206第七题 function f=wuyan1(c) q=[0.4518 0.4862 0.5295 0.5934 0.7171 0.8964 1.0202 1.1963 1.4928 1.6909 1.8548 2.1618 2.6638 3.4634 4.6759 5.8478 6.7885 7.4463 7.8345 8.2068 8.9468 9.7315 10.5172 11.7390 13.6876 ]; k=[0.0911 0.0961 0.1230 0.1430 0.1860 0.2543 0.3121 0.3792 0.4754 0.4410 0.4517 0.5595 0.8080 1.3072 1.7042 2.0019 2.2914 2.4941 2.8406 2.9855 3.2918 3.7214 4.3500 5.5567 7.0477]; l=[4.2361 4.3725 4.5295 4.6436 4.8179 4.9873 5.1282 5.2783 5.4334 5.5329 6.4749 6.5491 6.6152 6.6808 6.7455 6.8065 6.8950 6.9820 7.0637 7.1394 7.2085 7.3025 7.3470 7.4432 7.5200]; f=q-c(1)*k.^c(2).*l.^c(3) c0=[1 1 1] c=lsqnonlin('wuyan1',c0) c = 0.4091 0.6401 1.1446 a=0.4091 α=0.6401 β=1.1446 P239第五题 c=[-20 -30]; A=[1 2;5 4]; b=[20 70]; v1=[0 0]; [x,f,ef,out,lag]=linprog(c,A,b,[],[],v1) z=-f x = 10.0000 5.0000
流体力学 实验指导书与报告 静力学实验 雷诺实验 中国矿业大学能源与动力实验中心
学生实验守则 一、学生进入实验室必须遵守实验室规章制度，遵守课堂纪律，衣着整洁，保持安静，不得迟到早退，严禁喧哗、吸烟、吃零食和随地吐痰。如有违犯，指导教师有权停止基实验。 二、实验课前，要认真阅读教材，作好实验预习，根据不同科目要求写出预习报告，明确实验目的、要求和注意事项。 三、实验课上必须专心听讲，服从指导教师的安排和指导，遵守操作规程，认真操作，正确读数，不得草率敷衍，拼凑数据。 四、预习报告和实验报告必须独自完成，不得互相抄袭。 五、因故缺课的学生，可向指导教师申请一次补做机会，不补做的，该试验以零分计算，作为总成绩的一部分，累计三次者，该课实验以不及格论处，不能参加该门课程的考试。 六、在使用大型精密仪器设备前，必须接受技术培训，经考核合格后方可使用，使用中要严格遵守操作规程，并详细填写使用记录。 七、爱护仪器设备，不准动用与本实验无关的仪器设备。要节约水、电、试剂药品、元器件、材料等。如发生仪器、设备损坏要及时向指导教师报告，属责任事故的，应按有关文件规定赔偿。 八、注意实验安全，遵守安全规定，防止人身和仪器设备事故发生。一旦发生事故，要立即向指导教师报告，采取正确的应急措施，防止事故扩大，保护人身安全和财产安全。重大事故要同时保护好现场，迅速向有关部门报告，事故后尽快写出书面报告交上级有关部门，不得隐瞒事实真相。 九、试验完毕要做好整理工作，将试剂、药品、工具、材料及公用仪器等放回原处。洗刷器皿，清扫试验场地，切断电源、气源、水源，经指导教师检查合格后方可离开。 十、各类实验室可根据自身特点，制定出切实可行的实验守则，报经系(院)主管领导同意后执行，并送实验室管理科备案。 1984年5月制定 2014年4月再修订 中国矿业大学能源与动力实验中心
流体力学实验思考题 参考答案 流体力学实验室二○○六年静水压强实验1．同一静止液体内的测压管水头线是根什么线？测压管水头指z p ，即静水力学实验仪显示的测压管液面至基准面的垂直高度。测压管水头线指测压管液面的连线。实验直接观察可知，同一静止液面内的测压管水头线是一根水平线。 2．当p B 0 时，试根据记录数据，确定水箱内的真空区域。 p B 0 ，相应容器的真空区域包括以下三个部分： （1）过测压管2 液面作一水平面，由等压面原理知，相对测压管2及水箱内的水体而 言，该水平面为等压面，均为大气压强，故该平面以上由密封的水、气所占区域，均为真空区域。 （2）同理，过箱顶小不杯的液面作一水平面，测压管 4 中，该平面以上的水体亦为真 空区域。 （3）在测压管5 中，自水面向下深度某一段水柱亦为真空区域。这段高度与测压管2液面低于水箱液面的高度相等，亦与测压管4 液面高于小水杯液面高度相等。3．若再备一根直尺，试采用另外最简便的方法测定0 。 最简单的方法，是用直尺分别测量水箱内通大气情况下，管5 油水界面至水面和油水界面至油面的垂直高度h和h0 ，由式w h w 0h0 ，从而求得0 。4．如测压管太细，对于测压管液面的读数将有何影响？ 设被测液体为水，测压管太细，测压管液面因毛细现象而升高，造成测量误差，毛细高度由下式计算 式中，为表面张力系数；为液体容量；d 为测压管的内径；h 为毛细升高。常温的水， 0.073N m ，0.0098N m3。水与玻璃的浸润角很小，可以认为cos 1.0。 于是有 h 29.7 d （h 、d 均以mm 计） 一般来说，当玻璃测压管的内径大于10 mm时，毛细影响可略而不计。另外，当水质 不洁时，减小，毛细高度亦较净水小；当采用有机下班玻璃作测压管时，浸润角较大，其h 较普通玻璃管小。如果用同一根测压管测量液体相对压差值，则毛细现象无任何影响。因为测量高、低压强时均有毛细现象，但在计算压差时，互相抵消了。 5．过C 点作一水平面，相对管1、2、5 及水箱中液体而言，这个水平面是不是等压面？哪一部分液体是同一等压面？ 不全是等压面，它仅相对管1、2 及水箱中的液体而言，这个水平面才是等压面。因为只有全部具有下列5 个条件的平面才是等压面：（1）重力液体；（2）静止；（3）连通；（4）连通介质为同一均质液体；（5）同一水平面。而管5 与水箱之间不符合条件（4），相对管5 和水箱中的液体而言，该水平面不是水平面。
清华大学数学实验报告4
————————————————————————————————作者: ————————————————————————————————日期： ?
电13 苗键强2011010645
一、实验目的 1.掌握用 MATLAB 软件求解非线性方程和方程组的基本用法, 并对结果作初步分析； 2.练习用非线性方程和方程组建立实际问题的模型并进行求解。 二、实验内容 题目1 【问题描述】 (Q1）小张夫妇以按揭方式贷款买了1套价值２0万元的房子，首付了5万元，每月还款1000元，15年还清。问贷款利率是多少？ (Q2)某人欲贷款50 万元购房,他咨询了两家银行,第一家银行 开出的条件是每月还45０0元,15 年还清；第二家银行开出的条件是每年还４5０00 元,2０年还清。从利率方面看,哪家银行较优惠(简单假设:年利率=月利率×12）? 【分析与解】 假设初始贷款金额为x０,贷款利率为ｐ,每月还款金额为x,第i 个月还完当月贷款后所欠银行的金额为x i,(i=1,2，3,.．....，n）。由题意可知： x1=x0(1+p)?x x2=x0(1+p)2?x(1+p)?x x3=x0(1+p)3?x(1+p)2?x(1+p)?x ……
x n=x0(1+p)n?x(1+p)n?1???x(1+p)?x =x0(1+p)n?x (1+p)n?1 p =0 因而有： x0(1+p)n=x (1+p)n?1 p (1) 则可以根据上述方程描述的函数关系求解相应的变量。 （Ｑ1) 根据公式(１)，可以得到以下方程： 150p(1+p)180?(1+p)180+1=0 设 f(p)=150p(1+p)180?(1+p)180+1，通过计算机程序绘制ｆ（p）的图像以判断解ｐ的大致区间,在Ｍatlab中编程如下： ｆorｉ = 1:25 t = ０．0001＊i; p(i) = t； f(i) =１50*ｔ＊(１+t).^180-(１+t).^180+1; ｅｎd; ｐｌｏｔ(ｐ,ｆ),hｏld ｏn,grid on; 运行以上代码得到如下图像：
实验2 方程模型及其求解算法 一、实验目的及意义 [1] 复习求解方程及方程组的基本原理和方法； [2] 掌握迭代算法； [3] 熟悉MATLAB软件编程环境；掌握MATLAB编程语句(特别是循环、条件、控制等语句)； [4] 通过范例展现求解实际问题的初步建模过程； 通过该实验的学习，复习和归纳方程求解或方程组求解的各种数值解法（简单迭代法、二分法、牛顿法、割线法等），初步了解数学建模过程。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。 二、实验内容 1．方程求解和方程组的各种数值解法练习 2．直接使用MATLAB命令对方程和方程组进行求解练习 3．针对实际问题，试建立数学模型，并求解。 三、实验步骤 1．开启软件平台——MATLAB，开启MATLAB编辑窗口； 2．根据各种数值解法步骤编写M文件 3．保存文件并运行； 4．观察运行结果(数值或图形)； 5．根据观察到的结果写出实验报告，并浅谈学习心得体会。 四、实验要求与任务 基础实验 1．用图形放大法求解方程x sin(x) = 1. 并观察该方程有多少个根。 画出图形程序： x=-10:0.01:10; y=x.*sin(x)-1; y1=zeros(size(x)); plot(x,y,x,y1) MATLAB运行结果：
-10-8-6-4-20246810 -8-6 -4 -2 2 4 6 8 扩大区间画图程序： x=-50:0.01:50; y=x.*sin(x)-1; y1=zeros(size(x)); plot(x,y,x,y1) MATLAB 运行结果： -50-40-30-20-1001020304050 由上图可知，该方程有偶数个无数的根。
流体力学课程实验思考题解答 （一）流体静力学实验 1、 同一静止液体内的测压管水头线是根什么线？ 答：测压管水头指γ p Z + ，即静水力学实验仪显示的测压管液面至基准面的垂直高度。测 压管水头线指测压管液面的连线。从表1.1的实测数据或实验直接观察可知，同一静止液面的测压管水头线是一根水平线。 2、 当0东华大学MATLAB数学实验第二版答案(胡良剑)
东华大学M A T L A B数学实验第二版答案(胡良 剑) -CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN
数学实验答案 Chapter 1 Page20,ex1 (5) 等于[exp(1),exp(2);exp(3),exp(4)] (7) 3=1*3, 8=2*4 (8) a为各列最小值，b为最小值所在的行号 (10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture (11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10) (12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10) Page20, ex2 (1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c 相等，但他们不等于b (2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码 Page20,ex3 >> r=2;p=0.5;n=12; >> T=log(r)/n/log(1+0.01*p) Page20,ex4 >> x=-2:0.05:2;f=x.^4-2.^x; >> [fmin,min_index]=min(f) 最小值最小值点编址 >> x(min_index) ans = 0.6500 最小值点 >> [f1,x1_index]=min(abs(f)) 求近似根--绝对值最小的点 f1 = 0.0328 x1_index = 24 >> x(x1_index) ans = -0.8500 >> x(x1_index)=[];f=x.^4-2.^x; 删去绝对值最小的点以求函数绝对值次小的点>> [f2,x2_index]=min(abs(f)) 求另一近似根--函数绝对值次小的点 f2 = 0.0630 x2_index = 65 >> x(x2_index) ans =
最新大学工程流体力学实验-参考答案 参考答案 流体力学实验室 二○○六年 静水压强实验 1．同一静止液体内的测压管水头线是根什么线？ 测压管水头指γp z +，即静水力学实验仪显示的测压管液面至基准面的垂直高度。测压管水头线指测压管液面的连线。实验直接观察可知，同一静止液面内的测压管水头线是一根水平线。 2．当0?B p 时，试根据记录数据，确定水箱内的真空区域。 0?B p ，相应容器的真空区域包括以下三个部分： （1）过测压管2液面作一水平面，由等压面原理知，相对测压管2及水箱内的水体而言，该水平面为等压面，均为大气压强，故该平面以上由密封的水、气所占区域，均为真空区域。 （2）同理，过箱顶小不杯的液面作一水平面，测压管4中，该平面以上的水体亦为真空区域。 （3）在测压管5中，自水面向下深度某一段水柱亦为真空区域。这段高度与测压管2液面低于水箱液面的高度相等，亦与测压管4液面高于小水杯液面高度相等。 3．若再备一根直尺，试采用另外最简便的方法测定0γ。 最简单的方法，是用直尺分别测量水箱内通大气情况下，管5油水界面至水面和油水界面至油面的垂直高度h 和0h ，由式00h h w w γγ= ，从而求得0γ。 4．如测压管太细，对于测压管液面的读数将有何影响？ 设被测液体为水，测压管太细，测压管液面因毛细现象而升高，造成测量误差，毛细高度由下式计算 式中，σ为表面张力系数；γ为液体容量；d 为测压管的内径；h 为毛细升高。常温的水，m N 073.0=σ，30098.0m N =γ。水与玻璃的浸润角θ很小，可以认为0.1cos =θ。于是有 d h 7.29= （h 、d 均以mm 计） 一般来说，当玻璃测压管的内径大于10mm 时，毛细影响可略而不计。另外，当水质不洁时，σ减小，毛细高度亦较净水小；当采用有机下班玻璃作测压管时，浸润角θ较大，其h 较普通玻璃管小。
大学数学实验心得体会 [模版仅供参考，切勿通篇使用] 大学数学实验心得体会（一） 数学，在整个人类生命进程中至关重要，从小学到中学，再到大学，乃至更高层次的科学研究都离不开数学，随着时代的发展，人们越来越重视数学知识的应用，对数学课程提出了更高层次的要求，于是便诞生了数学实验。 学期最初，大学数学实验对于我们来说既熟悉又陌生，在我们的记忆中，我们做过物理实验、化学实验、生物实验，故然我们以为数学实验与它们一样，当我们在网上搜索有关数学实验的信息时，我们才知道，大学数学实验作为一门新兴的数学课程在近十年来取得了迅速的发展。数学实验以计算机技术和数学软件为载体，将数学建模的思想和方法融入其中，现在已经成为一种潮流。 当我们怀着好奇的心情走进屈静国老师的数学实验课堂时，我们才渐渐懂得，数学实验是一门有关计算机软件的课程，就像c语言一样，需要编辑运行程序，从而进行数学运算，它不需要自己来运算，就像计算器一样，只要我们自己记下重要程序语句，输入运行程序，便可得到运行结果，大大降低了我们的运算量，
给我们生活带来许多便捷，在大一时，我学过c语言，由于这样的基础，让我能够更快的学会并应用此软件。 时间飞逝，转眼间，我们就要结课了，这学期我们学习了mathematics的基础，微积分实验，线性代数实验，概率论与数理统计实验，数值计算方法及实验。通过这学期的学习，我也积累了些自己的学习方法和心得。首先，我们要在平时上课牢记那些mathematics语言和公式，那些东西就想单词和公式一样，只需要背诵;然后，我们要看几遍书，并多看一下例题;最后，我们要多应用mathematics软件去练习。正所谓熟能生巧，我坚信，只要我们能够做到这三步，我们就能很好的掌握这门课程。 通过学习使用数学软件，数学实验建模，使我们能够从实际问题出发，认真分析研究，建立简单数学模型，然后借助先进的计算机技术，最终找出解决实际问题的一种或多种方案，从而提高了我们的数学思维能力，为我们参加数学竞赛和数学建模打下了坚实的基础，同时也为我们进一步深造和参加工作打下一定的实践基础! 大学数学实验心得体会（二） 在此期间我充分利用研修活动时间学习，感到既有辛苦，又有收获。既有付出，又有新所得。这次远程研修让我有幸与专家和各地的数学精英们交流，面对每次探讨的主题，大家畅所欲言，
工程流体力学及水力学实验报告实验分析与讨论 1.同一静止液体内的测管水头线是根什么线？ 测压管水头指，即静水力学实验仪显示的测管液面至基准面的垂直高度。测压管水头线指测 压管液面的连线。实验直接观察可知，同一静止液面的测压管水头线是一根水平线。 2.当P B <0时，试根据记录数据，确定水箱内的真空区域。 ，相应容器的真空区域包括以下三部分： (1)过测压管2液面作一水平面，由等压面原理知，相对测压管2及水箱内的水体而言，该水平面为等压面，均为大气压强，故该平面以上由密封的水、气所占的空间区域，均为真空区域。 (2)同理，过箱顶小水杯的液面作一水平面，测压管4中，该平面以上的水体亦为真空区域。 (3)在测压管5中，自水面向下深度某一段水柱亦为真空区。这段高度与测压管2液面低于水箱液面的高度相等，亦与测压管4液面高于小水杯液面高度相等。 3.若再备一根直尺，试采用另外最简便的方法测定γ 。 最简单的方法，是用直尺分别测量水箱内通大气情况下，管5油水界面至水面和油水界面至油面的垂 直高度h和h 0，由式，从而求得γ 。 4.如测压管太细，对测压管液面的读数将有何影响？ 设被测液体为水，测压管太细，测压管液面因毛细现象而升高，造成测量误差，毛细高度由下式计算 式中，为表面张力系数；为液体的容量；d为测压管的内径；h为毛细升高。常温(t=20℃)的水，=7.28dyn/mm， =0.98dyn/mm。水与玻璃的浸润角很小，可认为cosθ=1.0。于是有(h、d单位为mm) 一般来说，当玻璃测压管的内径大于10mm时，毛细影响可略而不计。另外，当水质不洁时，减小，毛细高度亦较净水小；当采用有机玻璃作测压管时，浸润角较大，其h较普通玻璃管小。 如果用同一根测压管测量液体相对压差值，则毛细现象无任何影响。因为测量高、低压强时均有毛细现象，但在计算压差时，互相抵消了。 5.过C点作一水平面，相对管1、2、5及水箱中液体而言，这个水平面是不是等压面？哪一部分液体是同一等压面？ 不全是等压面，它仅相对管1、2及水箱中的液体而言，这个水平面才是等压面。因为只有全部具备下列5个条件的平面才是等压面：（1）重力液体；（2）静止；（3）连通；（4）连通介质为同一均质液体；（5）同一水平面。而管5与水箱之间不符合条件（4），因此，相对管5和水箱中的液体而言，该水平面不是等压面。 6.用图1.1装置能演示变液位下的恒定流实验吗？ 关闭各通气阀门，开启底阀，放水片刻，可看到有空气由c进入水箱。这时阀门的出流就是变液位下的恒定流。因为由观察可知，测压管1的液面始终与c点同高，表明作用于底阀上的总水头不变，故为恒
第1章 绪论 选择题 【1.1】 按连续介质的概念，流体质点是指：（a ）流体的分子；（b ）流体内的固体颗粒； （c ）几何的点；（d ）几何尺寸同流动空间相比是极小量，又含有大量分子的微元体。 解：流体质点是指体积小到可以看作一个几何点，但它又含有大量的分子，且具有诸如速度、密度及压强等物理量的流体微团。 （d ） 【1.2】 与牛顿内摩擦定律直接相关的因素是：（a ）切应力和压强；（b ）切应力和剪切变 形速度；（c ）切应力和剪切变形；（d ）切应力和流速。 解：牛顿内摩擦定律是 d d v y τμ =，而且速度梯度d d v y 是流体微团的剪切变形速度 d d t γ，故d d t γ τμ=。 （b ） 【1.3】 流体运动黏度υ的国际单位是：（a ）m 2 /s ；（b ）N/m 2 ；（c ）kg/m ；（d ）N·s/m 2 。 解：流体的运动黏度υ的国际单位是/s m 2 。 （a ） 【1.4】 理想流体的特征是：（a ）黏度是常数；（b ）不可压缩；（c ）无黏性；（d ）符合RT p =ρ 。 解：不考虑黏性的流体称为理想流体。 （c ） 【1.5】当水的压强增加一个大气压时，水的密度增大约为：（a ）1/20 000；（b ） 1/1 000；（c ）1/4 000；（d ）1/2 000。 解：当水的压强增加一个大气压时，其密度增大约 95d 1 d 0.51011020 000k p ρ ρ -==???= 。 （a ） 【1.6】 从力学的角度分析，一般流体和固体的区别在于流体：（a ）能承受拉力，平衡时 不能承受切应力；（b ）不能承受拉力，平衡时能承受切应力；（c ）不能承受拉力，平衡时不能承受切应力；（d ）能承受拉力，平衡时也能承受切应力。 解：流体的特性是既不能承受拉力，同时具有很大的流动性，即平衡时不能承受切应力。 （c ） 【1.7】下列流体哪个属牛顿流体：（a ）汽油；（b ）纸浆；（c ）血液；（d ）沥青。 解：满足牛顿内摩擦定律的流体称为牛顿流体。 （a ） 【1.8】 15C o 时空气和水的运动黏度6215.210m /s υ-=?空气，621.14610m /s υ-=?水，这说明：在运动中（a ）空气比水的黏性力大；（b ）空气比水的黏性力小；（c ）空气 与水的黏性力接近；（d ）不能直接比较。 解：空气的运动黏度比水大近10倍，但由于水的密度是空气的近800倍，因此水的黏度反而比空气大近50倍，而黏性力除了同流体的黏度有关，还和速度梯度有 关，因此它们不能直接比较。 （d ） 【1.9】 液体的黏性主要来自于液体：（a ）分子热运动；（b ）分子间内聚力；（c ）易变形 性；（d ）抗拒变形的能力。解：液体的黏性主要由分子内聚力决定。 （b ）
二、判断题（“对”在题号前（）中打√×）（10分） （√）1、误差是指测量值与真值之差，即误差=测量值－真值，如此定义的误差反映的是测量值偏离真值的大小和方向，既有大小又有正负符号。 （×）2、残差（偏差）是指测量值与其算术平均值之差，它与误差定义一样。（√）3、精密度是指重复测量所得结果相互接近程度，反映的是随机误差大小的程度。 （√）4、测量不确定度是评价测量质量的一个重要指标，是指测量误差可能出现的范围。 （×）7、分光计设计了两个角游标是为了消除视差。 （×）9、调节气垫导轨水平时发现在滑块运动方向上不水平，应该先调节单脚螺钉再调节双脚螺钉。 （×）10、用一级千分尺测量某一长度（Δ仪=0.004mm），单次测量结果为N=8.000mm，用不确定度评定测量结果为N=（8.000±0.004）mm。 三、简答题（共15分） 1．示波器实验中，（1）CH1（x）输入信号频率为50Hz，CH2（y）输入信号频率为100Hz；（2）CH1（x）输入信号频率为150Hz，CH2（y）输入信号频率为50Hz；画出这两种情况下，示波器上显示的李萨如图形。（8分）
差法处理数据的优点是什么？（7分） 答：自变量应满足等间距变化的要求，且满足分组要求。（4分） 优点：充分利用数据；消除部分定值系统误差 四、计算题（20分，每题10分） 1、用1/50游标卡尺，测得某金属板的长和宽数据如下表所示，求金属板的面 解：（1）金属块长度平均值：)(02.10mm L = 长度不确定度： )(01.03/02.0mm u L == 金属块长度为：mm L 01.002.10±= %10.0=B （2分） （2）金属块宽度平均值：)(05.4mm d = 宽度不确定度： )(01.03/02.0mm u d == 金属块宽度是：mm d 01.005.4±= %20.0=B （2分） （3）面积最佳估计值：258.40mm d L S =?= 不确定度：2222222 221.0mm L d d s L s d L d L S =+=??? ????+??? ????=σσσσσ 相对百分误差：B ＝%100?S s σ=0.25％ （4分） （4）结果表达：21.06.40mm S ±= B ＝0.25％ （2分） 注：注意有效数字位数，有误者酌情扣 5、测量中的千分尺的零点误差属于已定系统误差；米尺刻度不均匀的误差属于未
中国石油大学（华东）工程流体力学实验报告 实验日期：2014.12.11成绩： 班级：石工12-09学号：12021409姓名：陈相君教师：李成华 同组者：魏晓彤，刘海飞 实验二、能量方程（伯诺利方程）实验 一、实验目的 1．验证实际流体稳定流的能量方程； 2．通过对诸多动水水力现象的实验分析，理解能量转换特性； 3．掌握流速、流量、压强等水力要素的实验量测技能。 二、实验装置 本实验的装置如图2-1所示。 图2-1 自循环伯诺利方程实验装置 1.自循环供水器； 2.实验台； 3.可控硅无极调速器；4溢流板；5.稳水孔板； 6.恒压水箱； 7.测压机；8滑动测量尺；9.测压管；10.试验管道； 11.测压点；12皮托管；13.试验流量调节阀 说明 本仪器测压管有两种： （1）皮托管测压管（表2-1中标﹡的测压管），用以测读皮托管探头对准点的总水头； （2）普通测压管（表2-1未标﹡者），用以定量量测测压管水头。 实验流量用阀13调节，流量由调节阀13测量。
三、实验原理 在实验管路中沿管内水流方向取n 个过水断面。可以列出进口断面（1）至另一断面（i ）的能量方程式（i =2,3,…,n ） i w i i i i h g v p z g p z -++ + =+ + 1222 2 111 1αγυαγ 取12n 1a a a ==???==，选好基准面，从已设置的各断面的测压管中读出 z+p/r 值，测 出透过管路的流量，即可计算出断面平均流速，从而即可得到各断面测压管水头和总水头。 四、实验要求 1．记录有关常数实验装置编号 No._4____ 均匀段1d = 1.40-210m ?；缩管段2d =1.01-210m ?；扩管段3d =2.00-2 10m ?； 水箱液面高程0?= 47.6-2 10m ?；上管道轴线高程z ?=19 -2 10m ? （基准面选在标尺的零点上） 2．量测（p z γ + ）并记入表2-2。 注：i i i p h z γ =+ 为测压管水头，单位：-2 10m ，i 为测点编号。 3．计算流速水头和总水头。
高等数学数学实验报告实验人员：院（系） ___________学号_________姓名____________ 实验地点：计算机中心机房 实验一 一、实验题目： 根据上面的题目，通过作图，观察重要极限：lim（1+1/n）n=e 二、实验目的和意义 方法的理论意义和实用价值。 利用数形结合的方法观察数列的极限，可以从点图上看出数列的收敛性，以及近似地观察出数列的收敛值；通过编程可以输出数列的任意多项值，以此来得到数列的收敛性。通过此实验对数列极限概念的理解形象化、具体化。 三、计算公式（1+1/n）n 四、程序设计 五、程序运行结果 六、结果的讨论和分析 当n足够大时，所画出的点逐渐接近于直线，即点数越大，精确度越高。对于不同解题方法最后均能获得相同结果，因此需要择优，从众多方法中尽可能选择简单的一种。程序编写需要有扎实的理论基础，因此在上机调试前要仔细审查细节，对程序进行尽可能的简化、改进与完善。 实验二 一、实验题目 制作函数y=sin cx的图形动画，并观察参数c对函数图形的影响。 二、实验目的和意义 本实验的目的是让同学熟悉数学软件Mathematica所具有的良好的作图功能，并通过函数图形来认识函数，运用函数的图形来观察和分析函数的有关性态，建立数形结合的思想。 三、计算公式：y=sin cx 四、程序设计 五、程序运行结果
六、结果的讨论和分析 c 的不同导致函数的区间大小不同。 实验三 一、实验题目 观察函数f(x)=cos x 的各阶泰勒展开式的图形。 二、实验目的和意义 利用Mathematica 计算函数)(x f 的各阶泰勒多项式，并通过绘制曲线图形，来进一步掌握泰勒展开与函数逼近的思想。 三、计算公式 四、程序设计 五、程序运行结果 六、结果的讨论和分析 函数的泰勒多项式对于函数的近似程度随着阶数的提高而提高，但是对于任一确定次数的多项式，它只在展开点附近的一个局部范围内才有较好的近似精确度。 实验四 一、实验题目 计算定积分的黎曼和 二、实验目的和意义 在现实生活中许多实际问题遇到的定积分，被积函数往往不能用算是给出，而通过图像或表格给出；或虽然给出，但是要计算他的原函数却很困难，甚至原函数非初等函数。本实验目的，就是为了解决这些问题，进行定积分近似计算。 三、计算公式 四、程序设计 五、程序运行结果 六、结果的讨论和分析 本实验求的近似值由给出的n 的值的不同而不同。给出的n 值越大，得到的结果越接近准确的
第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- 程序： syms x limit((1001*x-sin(1001*x))/x^3,x,0) 结果： 程序： syms x limit((1001*x-sin(1001*x))/x^3,x,inf) 结果： 0 cos 1000 x mx y e =，求''y 程序： syms x diff(exp(x)*cos(1001*x/1000),2) 结果： -2001/1000000*exp(x)*cos(1001/1000*x)-1001/500*exp(x)*sin(1001/1000*x)
计算 2 2 11 00 x y e dxdy +?? 程序： dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) 结果： 计算4 2 2 4x dx m x +? 程序： syms x int(x^4/(1000^2+4*x^2)) 结果： (10)cos , x y e mx y =求 程序： syms x diff(exp(x)*cos(1000*x),10) 结果： 给出 0x =的泰勒展式(最高次幂为4). 程序： syms x taylor(sqrt(1001/1000+x),5) 结果： Fibonacci 数列{}n x 的定义是121,1x x ==， 12,(3,4,)n n n x x x n --=+=L 用循环语句编程给出该数列的前20项（要求将结果用向量的形式给出）。 程序： x=[1,1]; for n=3:20 x(n)=x(n-1)+x(n-2); end x 结果： Columns 1 through 10 1 1 2 3 5 8 13 21 3 4 5 5 Columns 11 through 20 89 144 233 377 610 987 1597 2584 4181 6765