4、.075m的穩(wěn)定值,即時(shí)。將代入上述模型,得到h(t)=899209.148v^2當(dāng)h(t)<=0.075時(shí),解得v<=0.00028802由此得到,當(dāng)v<=0.00028802m/s時(shí),水不溢出;當(dāng)v>0.00028802m/s時(shí),水溢出。下面,用matlab對(duì)這一結(jié)論進(jìn)一步分析。不妨取v=0.0002、0.000288分別代入方程(*)求其數(shù)值解,并作出圖形:編寫a.m文件:functionhp=a(t,h1)hp=(1.374038*0.0002-0.001449*sqrt(h1))/sqrt(0.15*h1-h1^2);functionhp=b
5、(t,h2)hp=(1.374038*0.0003-0.001449*sqrt(h2))/sqrt(0.15*h2-h2^2);執(zhí)行matlab程序:>>[t,h1]=ode45('a',[0,60],0.01);>>plot(t,h1)>>[t,h2]=ode45('b',[0,60],0.01);plot(t,h2)我國(guó)氣象部門規(guī)定24h降水量在200mm以上(約0.000002m/s)的雨為特大暴雨。對(duì)于這種情形,v>0.00028802m/s即24h降水量在24884.9mm以上的強(qiáng)降雨機(jī)率幾乎為0,因此,這個(gè)公司的承諾是能兌現(xiàn)的。(2)若為周
6、期函數(shù),不妨設(shè)為正弦函數(shù),即:這表明下雨過(guò)程是在60s內(nèi)發(fā)生的一個(gè)短促的強(qiáng)陣雨過(guò)程,最大的降雨強(qiáng)度是0.00001m/s,由方程(*)得到如下微分方程:運(yùn)用matlab求數(shù)值解:編寫shuicao.m文件:functionf=shuicao(t,h)f=(0.000412*sin(pi*t/60)-0.001449*sqrt(h))/sqrt(0.15*h-h^2);執(zhí)行matlab程序:>>[t,h]=ode45('shuicao',[0,60],0.01);>>plot(t,h)從圖中可以看出,h(t)的最大值不會(huì)超過(guò)0.075m,因此,對(duì)于第(2
7、)種情況,水槽的水也不會(huì)溢出,這個(gè)公司的承諾是能兌現(xiàn)的。5.模型的改進(jìn)基于長(zhǎng)時(shí)間特大暴雨的考慮,可做如下兩種改進(jìn):(1)增大排水管的橫截面積,即增大排水管的直徑d。(2)改變水槽的連接方式,讓水槽往屋檐傾斜一定角度,這相當(dāng)于增加水槽的容水高度。