水文預(yù)報(bào)課程設(shè)計(jì)(編輯修改稿)
2025-07-21 13:23 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 18 21 24 263 6 9 12 15 18 21 24 273 6 9 12 15 18 21 24 283 將直接徑流、地下徑流出流過(guò)程和總流量過(guò)程繪制出洪水流量過(guò)程線，見(jiàn)下圖9.圖9 洪水流量過(guò)程線第五章 總結(jié)和心得此次課程設(shè)計(jì)，做了很久，期間碰到諸多棘手的問(wèn)題，終于一一解決，得以完成，雖然還有些不盡人意，但總體上還是較好得完成了此次課程設(shè)計(jì)的各項(xiàng)要求。通過(guò)此次課程設(shè)計(jì)，收獲良多。首先是對(duì)于水文預(yù)報(bào)這門(mén)課程，為了完成這次課設(shè)，對(duì)課本進(jìn)行了深入系統(tǒng)的復(fù)習(xí)，尤其是與新安江模型相關(guān)的內(nèi)容，使得對(duì)這門(mén)課程的掌握更加扎實(shí)牢固，理解也更加深入；其二，本次課設(shè)采用的是matlab軟件進(jìn)行數(shù)據(jù)處理，這也是基于數(shù)據(jù)較多，計(jì)算分析比較復(fù)雜的緣故，利用原有的一點(diǎn)兒編程基礎(chǔ)，在這期間進(jìn)一步不斷地學(xué)習(xí)，對(duì)該軟件的使用能力有了很大程度的增強(qiáng)，這對(duì)以后的學(xué)習(xí)工作都是大有裨益的。這次課設(shè)的意義就是在于不斷逼迫自己去學(xué)習(xí)更多的新東西，并把以前學(xué)過(guò)的東西進(jìn)行整合貫通，達(dá)到提升自身水平的效果；其三，在做課程設(shè)計(jì)的過(guò)程中，老師對(duì)我們進(jìn)行了多次耐心認(rèn)真的輔導(dǎo)，同學(xué)們之間也不斷進(jìn)行互相的交流，也得到了很多的樂(lè)趣?？偟膩?lái)說(shuō)，這次課程設(shè)計(jì)受益匪淺。這是畢業(yè)設(shè)計(jì)之前的最后一次課設(shè)，也是由書(shū)本上知識(shí)運(yùn)用到實(shí)踐中的一次嘗試，讓我們意識(shí)到自身知識(shí)的匱乏和有限，長(zhǎng)嘆“書(shū)到用時(shí)方恨少”，同時(shí)這也是一個(gè)有趣的過(guò)程。試想，通過(guò)自己的不懈努力終于把它攻克，把成果展現(xiàn)在自己眼前時(shí)的喜悅，又有什么能比得上？感謝老師，感謝此次課設(shè)。附件：，前三年率定期和后兩年檢驗(yàn)期的計(jì)算徑流量與實(shí)測(cè)徑流量絕對(duì)誤差和相對(duì)誤差計(jì)算程序：%**************************************************************%% 19871989年計(jì)算結(jié)果 %% %%**************************************************************%clcclear ZL=load(39。39。)。%ZL表示導(dǎo)入的數(shù)據(jù)WM=140。UM=20。LM=60。DM=60。B=。C=。 IM=。%流域平均張力水容量WM（mm),上層張力水容量UM（mm),下土層張力水容量LM（mm),深層張力水容量DM(mm),張力水蓄水容量曲線方次B，深層蒸散發(fā)折算系數(shù)C，%不透水面積占全流域面積的比例IMQ=ZL(:,3)。 %日徑流量實(shí)測(cè)值Q（m3/s）E0=ZL(:,4)。 %日蒸發(fā)量蒸發(fā)皿實(shí)測(cè)值（mm）P1=ZL(:,5)。P2=ZL(:,6)。P3=ZL(:,7)。P4=ZL(:,8)。 %P1,P2,P3,P4分別表示流域四個(gè)地區(qū)的日降雨量實(shí)測(cè)值（mm）sumQ=0。 %sumQ表示年徑流量實(shí)測(cè)值（mm）sumR=zeros(4000，1)。 h=length(ZL)。%求出數(shù)據(jù)的天數(shù)R=zeros(h,1) 。 %R表示日徑流量的理論計(jì)算值EU=zeros(1,h)。 for i = 1:h %流域平均降雨量計(jì)算P(i) = * P1(i) + * P2(i) + * P3(i) + * P4(i)。 %逐日降雨深的計(jì)算sumQ = sumQ + Q(i) * 24 * / 553。 %實(shí)際測(cè)得年徑流量的計(jì)算end%求逐日降雨深和實(shí)測(cè)全年徑流量W(1) = 110。 WU(1) = 10。 WL(1) = 40。 WD(1) = 60。 %初始的土壤總，上，下，深層土壤含水量%流域三層蒸發(fā)計(jì)算WMM = WM * (1 + B)。 %初始土壤含水量的確定a(1) = WMM * (1 (1 (W(1) / WM)) ^ (1 / (1 + B)))。 for j=1:4000 Kc(j)=+*j。 %， for i = 1:h EP(i) = E0(i) * Kc(j)。 %流域逐日蒸散發(fā)深if WU(i) + P(i) = EP(i) EU(i) = EP(i)。EL(i) = 0。ED(i) = 0。endif WU(i) + P(i) EP(i) if WL(i) = C * LM EU(i) = WU(i) + P(i)。 EL(i) = (EP(i) EU(i)) * WL(i) / LM。 ED(i) = 0。 elseif WL(i) C * LM amp。 WL(i) = C * (EP(i) EU(i)) EU(i) = WU(i) + P(i)。 EL(i) = (EP(i) EU(i)) * C。 ED(i) = 0。 elseif WL(i) C * (EP(i) EU(i)) EU(i) = WU(i) + P(i)。 EL(i) = WL(i)。 ED(i) = (EP(i) EU(i)) * C EL(i)。 end endE(i) = EU(i) + EL(i) + ED(i)。%求得各層的蒸發(fā)量 PE(i) = P(i) E(i)。 %流域產(chǎn)流計(jì)算凈降雨量PE%三層蒸散發(fā)計(jì)算，求總蒸發(fā)量和流域凈降水量if PE(i) 0 %當(dāng)產(chǎn)流時(shí) if PE(i) + a(i) WMM R(i) = PE(i) + W(i) WM + WM * (1 (PE(i) + a(i)) / WMM) ^ (B + 1)。 W(i + 1) = W(i) + PE(i) R(i)。 a(i + 1) = PE(i) + a(i)。 elseif PE(i) + a(i) = WMM R(i) = PE(i) + W(i) WM。 W(i + 1) = WM。 a(i + 1) = WMM。 end end %產(chǎn)流計(jì)算完畢if WU(i) + P(i) EU(i) R(i) = UM WU(i + 1) = WU(i) + P(i) EU(i) R(i)。 WL(i + 1) = WL(i) EL(i)。 WD(i + 1) = WD(i) ED(i)。 else WU(i + 1) = UM。 if WL(i) EL(i) + (WU(i) + P(i) EU(i) R(i) UM) = LM WL(i + 1) = WL(i) EL(i) + (WU(i) + P(i) EU(i) R(i) UM)。 WD(i + 1) = WD(i) ED(i)。 else WL(i + 1) = LM。 if WD(i) ED(i) +
物理相關(guān)推薦
鄂ICP備17016276號(hào)-1