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