小学期数学建模大作业 西安交通大学

Lp，。l，ll（）S(a)1PS(a)2（）a0ax017-3旋转成a角令使0Ll1/a与b2/++kkKCy；KZKk）KKkKKkkk；123/所以y=Cλ+Cλ+k+——a+kk121122取有y1=10+4.3508.^x+1。4192。^x；y2=（12-3.3266。＊x）。＊1。8182。^x；plot(x,y1，'g。',’markersize’，25）on；plot（x,y2，'r.''markersize'legend（Á½¸öʵ¸ù’Ò»¸öʵ¸)4/第三次作业eX=510］;W=［W;（abs（X(i)—X（j））+abs（Y（i）5/bcdefghi0700000用6/A1AB0.30令，=x1B00.70Acx则1x110cxxnn0ii0.30200.31540.32810.3318x1，x10，x30，x50，0.69800.68460.67190.66820.3333x1000.66671AB1133，=A1B223313nx2n3AB12557/5555555553456789108/根据速度平衡原理，对于房顶排水系统水槽中水的容量的变化率=雨水的流入速度-排水管流出的速度。()VtQQ10Q,Q分别是单位时间流入水槽和从水槽流出的雨水量的体积。10表示单位时间里落在水平面上雨水的深度，房顶的面积rtbdbrtbd实际受雨的水平面积cos,房顶上雨水的流速bdcostcos流入水槽的速度应是在铅垂方向的分量Qr1排水管的流出速度应与水槽中水的深度有关。根据能量守恒原理12()mvmghtt2()QAt2()v2水槽中水的体积为0，ah2h2hcos(t)adsincos)V2aa2ah(ah)2h2a2ahVt)ad1)sin2aa2a()2()2()()Vtdhtahth2trtbdcossinAt2()2t)2aht)ht)2tsincosA2ght)rdtht22()2()将表中的数据代入（7）式，用matlab解得9/19.6ht)9600()0.333v1.374038v-0.001449ht)dht)4)dtht0.15()ht0.15()ht()ht22由假设（1）若vt)常数。先讨论水槽的深度趋于一个低于0。075m的稳定值，即h(0)0.05时。h'(t)0将代入上述模型，得到h(t)=899209.148v^2h'(t)0当h（t)<=0。075时，解得由此得到，当m/s时，水不溢出；当v>m/s时，水溢出。下面,用matlab对这一结论进一步分析.不妨取编写a。m文件：v〈=0.00028802functionhp=a(t，h1)hp=(1。374038*0。0002—0。001449*sqrt（h1))/sqrt（0.15*h1—h1^2）；functionhp=b（t，h2)hp=（1.374038＊0.0003-0。001449＊sqrt(h2)）/sqrt（0.15*h2-h2^2);执行matlab程序：>>［t，h1］=ode45('a',［0，60］，0。01）;〉>plot(t，h1）0.040.0350.03m/0.025h0.020.0150.010204060t/s>>［t，h2］=ode45（'b’，[0，60]，0。01）;plot（t,h2)/0.080.060.040.020m/h0204060t/s我国气象部门规定24h降水量在200mm以上（约0。000002m/s）的雨为特大vt)常数m/s即24h降水量在24884.9mm以上的强降雨机率几乎为0，因此,这个公司的承诺是能兑现的。(2）若为周期函数,不妨设为正弦函数,即：v(t)t0.0003sin,0t这表明下雨过程是在60s内发生的vt)t一个短促的强阵雨过程,最大的降雨强度是0。000微分方程：t0.000412sin0.001449h60,0t60dhdt0.15hh20.001449h,t600.15hh2h(0)0.01运用matlab求数值解：编写shuicao.m文件:functionf=shuicao（t,h)f=(0.000412＊sin（pi＊t/60)—0.001449＊sqrt（h))/sqrt（0.15＊h-h^2）;执行matlab程序：〉>［t，h］=ode45(’shuicao',［0,60],0。01);>〉plot（t,h）/m/h020405060槽的水也不会溢出,这个公司的承诺是能兑现的。AABB/P（频率）0.70。20。1P（频率）0.30。40。20.10若XXx1=rand（10000，1)；x2=rand（10000,1)；i=1：10000x1（i)〈0。7T1=0；x1(i）〈0.9//T2=30+2*x3(i）；x2（i)<0。3T3=28；x2(i）<0。7x2(i)<0。9T3=34；T3〈T2+T1第六次作业（2）10至米/为,与x=N=T=M=r=========X==-）++=(—)+1］（％0。，==)/第七次作业(1）常都是极为复杂的,./2aarccos112S22x```````````221/S2abc2x````````````322x042xcca222arccos2xcxc2222matlab程序求解如下：a=39.5;b=9.5；c=1。7;x=0:0。001:10;t=1000；fori=1:10000ift〉Qt=Q;flag=i；endendflagx=flag＊0.001运行结果为：flag=10000，x=10.（巴西里约热内卢的救赎）当a=46。5,b=46。5时，结果为：flag=10000，x=10.（美国自由女神像）当a=26.4，b=7。599三次数据的a与b一定的比例关系，但人观察塑像时与塑像张成的夹角最大的位置相同。/年1年表1表2(/xz3x+=或x=-az(k）+b(k=2，(1）ˆ,ˆ)T则2xX故年//n的4aaaaaa、、、从Q到QD7即即11223344556677/1000000001000100110001010001000000110000001011001000011=1356794aaaaaaa1000000001000100110001010001000000110000001011001000011a39=135679/第九次作业(2）一0表1二+上一级降等到本级的职员人数=下一/三四xi/x2ix3ix4ix5ix6ix7iMn(n)Nn(