1、 论文题目：雨中行走淋雨量分析1雨中行走淋雨量分析摘要本文在给定的降雨条件下，分别建立相应的数学模型，分析人体在雨中奔跑时淋雨多少与奔跑速度、降雨方向等因素的关系。其中文中所涉及到的降雨量是指从天空降落到地面上的雨水，未经蒸发、渗透、流失而在水面上积聚的水层深度，它可以直观地表示降雨的多少。淋雨量，是指人在雨中行走时全身所接收到得雨的体积，可表示为单位时间单位面积上淋雨的多少与接收雨的面积和淋雨时间的乘积。利用 MATLAB 软件对各个问题进行了求解。针对问题一，设降雨淋遍全身不考虑雨的方向，经简化假设得人淋雨面积为前后左右及头顶面积之和。人以最大速度奔跑 1000m，用 MATLAB 求解可
2、得淋雨量近似为 0.0024 。3m针对问题二，雨迎面吹来，雨线方向与跑步方向在同一平面，人淋雨面积为前方和头顶面积之和。因各个方向上降雨速度分量不同，故分别计算头顶和前方的淋雨量后相加即为总的淋雨量。据此可列出总淋雨量 W 与跑步速度 v 之间的函数关系。分析表明当跑步速度为 maxv时，淋雨量最少。并计算出当雨与人体的夹角 =0 时，淋雨量近似为 0.0012 ；当 =30时，淋雨量近似为30.0016 。3m针对问题三，雨从背面吹来，雨线与跑步方向在同一平面内，人淋雨量与人和雨相对速度有关。列出函数关系式分析并求解，可知当人速度 v=2 时ms淋雨量最少， =30时的总淋雨量近似为 0.
3、2405556E-03 。3针对问题四，列出淋雨量 W 和跑步速度 v 之间的函数关系式，利用MATLAB 画出 分别为 0，10，.90的曲线图。针对问题五，雨线与人跑步方向不在同一平面内，则考虑人的淋雨面积为前后左右以及头顶。分别列式表示，总的淋雨量即为三者之和。关键词淋雨量；降雨的大小；降雨的方向（风） ；路程的远近；行走的速度；2一、 问题重述生活中我们常常会遇到下雨却没有遮雨工具的时刻，我们在那时会有很多选择，其中之一就是淋雨，往往很多人会在雨中快走或奔跑以使自己身体淋雨量最小化，但往往很多人会感觉到淋雨量并不会因为快走或奔跑而减少多少，反而有时候淋雨量倒有所增加，淋雨量和速度等有关
4、参数的关系如何，让我们假设一数学模型模拟计算真实情况。当我们在雨中从一处沿直线跑到另一处时，如果雨速为常数，走的时候身体的动作的大小和暴露在雨中的面积大小影响着淋雨的多少，并且行走速度也同样影响着淋雨量，将人体简化成一个长方体，高 a=1.5 米，宽 b=0.5 米，厚 c=0.2m，跑步距离 d=1000m，跑步最大速度 maxv =5m/s，雨速 u=4m/s，降雨量 w=2 /h= ，记跑步速度为 。65.10/s1.当我们不考虑雨的方向时，假设降雨会淋遍全身，这时如果我们最大速度奔跑会淋多少雨？2.雨从迎面吹来，雨线方向与跑步方向在同一平面内，且与人体的夹角为，建立总淋雨量与速度 及参
5、数 a b c d u 之间的关系。问速度 多大，总淋雨量最少。 计算 =0， =30时的总淋雨量。3.雨从背面吹来，设雨线方向与跑步方向在同一平面内，且与人体的夹角为 ，建立总淋雨量与速度 及参数 a b c d u w 之间的关系。问速度 多大，总淋雨量最少。计算 =30时的总淋雨量。4.以总淋雨量为纵轴，速度 为横轴对第 3 问作图（考虑 的影响），并解释结果的实际意义。5.若雨线方向与跑步方向不在同一平面内，模型会有什么变化？ 二、问题分析2.1 问题一分析3若不考虑雨的方向，雨以降雨量 w 均匀地淋遍全身。将人体简化成长方体，求出人接受雨的总面积，人以最大速度跑步，并计算淋雨时间、单
6、位时间、单位面积上的降雨量，求出人跑完全程的总淋雨量 W。2.2 问题二分析雨迎面吹来，雨线方向与跑步方向在同一平面内且与人体夹角为 ，如图1 所示。根据实际情况估计人体淋雨可分为头顶和前后左右几个方向上。雨迎面吹来时，由于雨相对于人的速度有变化，因此人单位时间内接收雨量变化，且与相对速度成正比。据此，推算出前后侧上单位时间接受雨量。同理，头顶部位接雨量与雨速垂直于头顶平面的分速度成正比。分别计算出头顶侧与前后侧单位时间接雨量，并分别乘以各自面积以及时间 d/t,即得到头顶及两侧淋雨的总量。在人体总的淋雨量.据此可得 W 与 v 之间关系，并能求出 =0 和 =30时的总淋雨量。图 12.3
7、问题三分析雨从背面吹来，雨线与跑步方向在同一平面内且与人体夹角为 ，如图 2 所示。左右方向上淋雨量为 0。头顶上单位时间内接收雨的量 与雨速垂直方向上的1w4分量成正比， 为头顶面积 bc 与时间的 d/v 以及 之积。当 时，前1W1wsinvu方不受雨，前后方向上单位时间内淋雨量 与人前进方向上人相对于雨的速度2（usin -v）成正比，据此推算出 ；而当 时，后方不受雨，由于人2sinvu速已经高于雨速，这时前面会向前撞上雨滴，即 与 成正比。 为2wsiv2W人体前面积 ab 和跑步时间 d/v 顶淋雨量以及 之积。由此可计算出总的淋雨量。 12W据此可得 W 与 v 之间关系，并能
8、求出 =30时的总淋雨量。图 22.4 问题四分析以总淋雨量 W 为纵轴、速度 为横，针对问题三的求解，利用 MATLAB作出当 分别为 0，10 ， 20，30，40，50，60 ，70 ，80，90时的曲线图并加以分析。2.5 问题五分析5csinbcos图 3 俯视图如图三，为人体模型的俯视图。需要分三部分计算，在前后面上，雨垂直方向分速度为 ，相对速度为 ，乘上垂直受雨的面积 ab 以cosusincovu及时间 即为前后侧受雨量 。因为垂直于左右面人的分速度为 0，左右两面dv2W上相对速度为 乘上面积 ac 以及时间 极为左右受雨量 .而头顶受sindv3W雨与雨速和人速的夹角大小
9、无关，因此 仍按（2） 、 （3）问的算法做。由1可得雨量求法公式。123W二、 模型假设1. 人在奔跑过程中， 大小与方向恒定，即沿直线匀速前进。2. 对问题 1 人体各个方向均匀接受雨量，即单位时间、单位面积上接受雨量恒定。3. 对问题 2、3 雨线与跑步方向在同一平面内，并且雨线与人体夹角不变。在此过程中左右两次因与雨速平行而不沾雨。4. 假设雨的密度相同，雨滴大小、形状相同，雨速均匀不变5. 假设单位时间内接收雨的量与雨速成正比。6. 将人体理想化为一个长、宽、高、已知的长方体模型，且人体行走过程中的震荡引起的误差可忽略不计。6三、 符号说明a 人体高度b 人体宽度c 人体厚度d 跑步
10、距离u 雨速w 降雨量 雨迎面吹来时与人体的夹角 俯视图中雨速与人速的夹角maxv跑步最大速度W 总淋雨量1s头顶面积2人前或后表面积1u雨点相对人头顶速度的垂直分量2雨点相对人前后面速度的垂直分量1w头顶单位时间接收雨量2前后面单位时间接收雨量1W头顶接收雨量2人体前后面接收雨量3人体左右面接收雨量maxv7五、 模型建立与求解5.1问题一不考虑雨的方向，因为降雨量 w 均匀地淋遍全身，所以在将人体简化成长方体的情况下，忽略次要因素，人以最大速度跑步，根据淋雨时间、单位时间、单位面积上的降雨量等有关条件，列出总淋雨量 W 的求解公式如下：max2dWabcv利用 MATLAB 编程求解（见附
11、录一） ，可得：0.2435.2问题二根据题意，将降落在人体上的雨滴分成两部分， （顶部） （前面） ，1s2s人体接收的雨量和头顶面积、头顶部分与雨滴垂直下落方向分量 、行走时间1u有关。列式求解如下：头顶：1cosub8假设降雨量 w 与与点密度（均匀不计）淋雨量与人相对速度有关，所以：11coscosudbwWtwv 正面： 2sinvu而 22212sinsinisinvuwvwuvabdWtvbdwcu +1利用 MATLAB 编程求解（见附录二） ，可得：当 v=5m/s 时，淋雨量 W 最小；当 =0时，W=0.0012 3m当 =30时， W=0.00165.3问题三根据题意，
12、根据题意，将降落在人体上的雨滴分成两部分， （顶部）1s（前后两面） ， 面积为 假设： 与雨点密度，雨点与人的相对速度2s1s1bc1w成正比而雨点均匀分布。头顶：9111coscoswvuwdtbW正面：当 时，人速大于垂直于人前后面的雨速，雨会沾到人的前面sinuv用 lingo 编程（见附录三）求解可得：当 v=2m/s 时，总淋雨量最少;雨线方向与人体夹角为 30时，淋雨量为 0.2405556E-03 。3m10图 45.4问题四根据问题三的结论，列出总的淋雨量 W 和人速度 v 之间的关系式，利用MATLAB 画出 取不同值时的函数图像如下：0 0.5 1 1.5 2 2.5 3
13、 3.5 4 4.5-0.01-0.00500.0050.010.0150.020.0250.03v轴W轴分析图像可知，当 v=2 时，总淋雨量最少。5.5问题五应用（3）中的结论 cosdbwv1W前后侧，当 时，相对速度 ，sincovuinu112sincovudWwabv可总结为 2sincovudabv同理，可得左右侧接收雨量 3siWwcuv三者相加得 sincosincovudddbwabwacvv +六、 模型评价 通过对本题的分析求解，可知道人在雨中奔跑的淋雨量不仅与跑步速度有关，还与雨线与人跑步方向的夹角，雨速以及人跑步速度等因素有关。本文忽略了降雨密度不均匀，风向不稳定等
14、次要因素，以便更好的对问题进行分析和研究。但在实际问题中的限制性因素远远超过这些，因此此文的分析方法仍存在一定的局限性，有待改进和提高。参考文献1 薛定宇，陈阳泉.高等应用数学问题的 MATLAB 求解.北京：清华大学出版社，2008 年 10 月。2 陈杰.MATLAB 宝典.北京：电子工业出版社，2007 年。12附录附录 1：问题一求解程序clear;a=1.5;b=0.5;c=0.2;w=0.02/3600;d=1000;Vm=5;W=(2*a*b+2*a*c+b*c)*w*d/VmW=0.0024附录 2：问题二求解程序附录 2.1：分析当 v=vm时总淋雨量最小程序clear;sy
15、ms t v;a=1.5;b=0.5;c=0.2;w=0.02/3600; d=1000;u=4;minss=(b*c*d*w*cos(t)+a*b*d*w*sin(t)/v+a*b*w*d/u;h=diff(minss,v) %导数h =-(1/1800*cos(t)+1/240*sin(t)/v2%h=-(1/1800*cos(t)+1/240*sin(t)/v2 恒小于零，原函数为减函数附录 2.2：当 =0 时总淋雨量程序clear;t=0;13v=5;a=1.5;b=0.5;c=0.2;w=0.02/3600; d=1000;u=4;minss=(b*c*d*w*cos(t)+a*b
16、*d*w*sin(t)/v+a*b*w*d/uminss =0.0012附录 2.3：当 =30时总淋雨量clear;t=pi/6;v=5;a=1.5;b=0.5;c=0.2;w=0.02/3600;d=1000;u=4;minss=(b*c*d*w*cos(t)+a*b*d*w*sin(t)/v+a*b*w*d/uminss =0.0016附录 3：问题三中 =30时总淋雨量程序一：min=0.5*0.2*1.732/2*1000/v*0.02/3600+(2-v)/2*1.5*0.5*0.02/3600*1000/v;v=0;v=2;v=5;运行结果：Local optimal solut
17、ion found at iteration: 8Objective value: 0.2405556E-03Variable Value Reduced CostV 2.000000 0.000000Row Slack or Surplus Dual Price1 0.2405556E-03 -1.0000002 0.000000 -0.9213884E-033 3.000000 0.000000附录4：问题四的程序t=sin(t0);v1=0:0.1:t(1);v11=t(1):0.1:5;y1=(5.556e-4.*cos(t(1)+41.67e-4.*sin(t(1)./v1-4.16
18、7e-3/4;y11=(5.556e-4.*cos(t(1)-41.67e-4.*sin(t(1)./v11+4.167e-3/4;plot(v1,y1,r);hold onplot(v11,y11,r);xlabel(v);ylabel(W);v2=0:0.1:t(2);v22=t(2):0.1:5;y2=(5.556e-4.*cos(t(2)+41.67e-4.*sin(t(2)./v2-4.167e-3/4;y22=(5.556e-4.*cos(t(2)-41.67e-4.*sin(t(2)./v22+4.167e-3/4;plot(v2,y2,g-.);hold onplot(v22,
19、y22,g-.);v3=0:0.1:t(3);15v33=t(3):0.1:5;y3=(5.556e-4.*cos(t(3)+41.67e-4.*sin(t(3)./v3-4.167e-3/4;y33=(5.556e-4.*cos(t(3)-41.67e-4.*sin(t(3)./v33+4.167e-3/4;plot(v3,y3,k:);hold onplot(v33,y33,k:);v4=0:0.1:t(4);v44=t(4):0.1:5;y4=(5.556e-4.*cos(t(4)+41.67e-4.*sin(t(4)./v4-4.167e-3/4;y44=(5.556e-4.*cos(
20、t(4)-41.67e-4.*sin(t(4)./v44+4.167e-3/4;plot(v4,y4,b);hold onplot(v44,y44,b);v5=0:0.1:t(5);v55=t(5):0.1:5;y5=(5.556e-4.*cos(t(5)+41.67e-4.*sin(t(5)./v5-4.167e-3/4;y55=(5.556e-4.*cos(t(5)-41.67e-4.*sin(t(5)./v55+4.167e-3/4;plot(v5,y5,c);hold onplot(v55,y55,c);v6=0:0.1:t(6);v66=t(6):0.1:5;y6=(5.556e-4
21、.*cos(t(6)+41.67e-4.*sin(t(6)./v6-4.167e-3/4;y66=(5.556e-4.*cos(t(6)-41.67e-4.*sin(t(6)./v66+4.167e-3/4;plot(v6,y6,y);hold onplot(v66,y66,y);v7=0:0.1:t(7);v77=t(7):0.1:5;y7=(5.556e-4.*cos(t(7)+41.67e-4.*sin(t(7)./v7-4.167e-3/4;y77=(5.556e-4.*cos(t(7)-41.67e-4.*sin(t(7)./v77+4.167e-3/4;plot(v7,y7,m);
22、hold onplot(v77,y77,m);v8=0:0.1:t(8);v88=t(8):0.1:5;y8=(5.556e-4.*cos(t(8)+41.67e-4.*sin(t(8)./v8-4.167e-3/4;y88=(5.556e-4.*cos(t(8)-41.67e-4.*sin(t(8)./v88+4.167e-3/4;16plot(v8,y8,k);hold onplot(v88,y88,k);v9=0:0.1:t(9);v99=t(9):0.1:5;y9=(5.556e-4.*cos(t(9)+41.67e-4.*sin(t(9)./v9-4.167e-3/4;y99=(5.
23、556e-4.*cos(t(9)-41.67e-4.*sin(t(9)./v99+4.167e-3/4;plot(v9,y9,m-.);hold onplot(v99,y99,m-.);v10=0:0.1:t(10);v1010=t(10):0.1:5;y10=(5.556e-4.*cos(t(10)+41.67e-4.*sin(t(10)./v10-4.167e-3/4;y1010=(5.556e-4.*cos(t(10)-41.67e-4.*sin(t(10)./v1010+4.167e-3/4;plot(v10,y10,k);hold onplot(v1010,y1010,k);legend(y1,y11,y2,y22,y3,y33,y4,y44,y5,y55,y6,y66,y7,y77,y8,y88,y9,y99,y10,y1010);t0=0:pi/10:pi/2;
