【正文】 WL(i)。 WL(i) = C * (EP(i) EU(i)) EU(i) = WU(i) + P(i)。endif WU(i) + P(i) EP(i) if WL(i) = C * LM EU(i) = WU(i) + P(i)。 for i = 1:m EP(i) = E0(i) * Kc。 WL(1) = 40。 %確定暴雨持續(xù)時(shí)間R=zeros(1,m)。P3=hs(:,6)。 %導(dǎo)入常數(shù)值。C=。UM=20。 %DC是確定性系數(shù)end：%計(jì)算步驟%首先進(jìn)行三層蒸發(fā)計(jì)算，得到凈降雨量，然后進(jìn)行徑流劃分，然后計(jì)算次日到初始土壤含水量%然后根據(jù)劃分出的各時(shí)段的地表和地下徑流，求取洪水過程 %**************************************************************%% 2004年場(chǎng)次暴雨的洪水預(yù)報(bào)計(jì)算 %% %%**************************************************************%clcclear %清空數(shù)據(jù) hs=load(39。Y=0。)。天數(shù)39。 a(i + 1) = WMM * (1 (1 W(i + 1) / WM) ^ (1 / (1 + B)))。 else WD(i + 1) = DM。 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)。 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)。 a(i + 1) = PE(i) + a(i)。 end endE(i) = EU(i) + EL(i) + ED(i)。 ED(i) = 0。 ED(i) = 0。EL(i) = 0。 %流域三層蒸發(fā)計(jì)算，先給出初始 上、下、深層土壤含水量WMM = WM * (1 + B)。 %逐日降雨深的計(jì)算end%求逐日降雨深 W(1) = 110。 %確定Kc取值m=length(hs)。P2=hs(:,6)。 %導(dǎo)入常數(shù)值。C=。UM=20。 %：此處僅以1987年為例%**************************************************************%% 1987年日徑流量理論值和實(shí)測(cè)值的檢驗(yàn) %% %%**************************************************************%%計(jì)算步驟%首先進(jìn)行三層蒸發(fā)計(jì)算，得到凈降雨量，然后進(jìn)行徑流劃分，然后計(jì)算次日到初始土壤含水量%然后根據(jù)劃分出的各時(shí)段的地表和地下徑流，求取洪水過程clcclear %清空數(shù)據(jù) hs=load(39。 if cha(j)=KC(i)=kc(j)。)。Kc值39。for j=1:4000kc(j)=+*j。 W(i + 1) = W(i) + PE(i)。 if WD(i) ED(i) + WL(i) EL(i) + (WU(i) + P(i) EU(i) R(i) UM) LM = DM WD(i + 1) = WD(i) ED(i) + WL(i) EL(i) + (WU(i) + P(i) EU(i) R(i) UM) LM。 else WU(i + 1) = UM。 a(i + 1) = WMM。 W(i + 1) = W(i) + PE(i) R(i)。 ED(i) = (EP(i) EU(i)) * C EL(i)。 EL(i) = (EP(i) EU(i)) * C。 EL(i) = (EP(i) EU(i)) * WL(i) / LM。 %流域逐日蒸散發(fā)深if WU(i) + P(i) = EP(i) EU(i) = EP(i)。 %初始的土壤總，上，下，深層土壤含水量%流域三層蒸發(fā)計(jì)算WMM = WM * (1 + B)。 %實(shí)際測(cè)得年徑流量的計(jì)算end%求逐日降雨深和實(shí)測(cè)全年徑流量W(1) = 110。%求出數(shù)據(jù)的天數(shù)R=zeros(h,1) 。P4=ZL(:,8)。 %日徑流量實(shí)測(cè)值Q（m3/s）E0=ZL(:,4)。B=。%ZL表示導(dǎo)入的數(shù)據(jù)WM=140。試想，通過自己的不懈努力終于把它攻克，把成果展現(xiàn)在自己眼前時(shí)的喜悅，又有什么能比得上？感謝老師，感謝此次課設(shè)。首先是對(duì)于水文預(yù)報(bào)這門課程，為了完成這次課設(shè)，對(duì)課本進(jìn)行了深入系統(tǒng)的復(fù)習(xí)，尤其是與新安江模型相關(guān)的內(nèi)容，使得對(duì)這門課程的掌握更加扎實(shí)牢固，理解也更加深入；其二，本次課設(shè)采用的是matlab軟件進(jìn)行數(shù)據(jù)處理，這也是基于數(shù)據(jù)較多，計(jì)算分析比較復(fù)雜的緣故，利用原有的一點(diǎn)兒編程基礎(chǔ)，在這期間進(jìn)一步不斷地學(xué)習(xí)，對(duì)該軟件的使用能力有了很大程度的增強(qiáng)，這對(duì)以后的學(xué)習(xí)工作都是大有裨益的。運(yùn)用單位線法和出流系數(shù)法分別計(jì)算直接徑流出流量和地下徑流出流量，兩者之和即為總的流量。（2）各年先對(duì)誤差均已經(jīng)很小，但不能全部達(dá)到5%的范圍之內(nèi)，可能是因?yàn)檎{(diào)試不夠，或者是因?yàn)槟Ｐ图僭O(shè)與市級(jí)的情況不盡相同。根據(jù)相對(duì)誤差規(guī)律得以下三種方法，同時(shí)說明：①通過人為觀察、比較，人工選取所給Kc的最優(yōu)值；②計(jì)算19871989年三年的相對(duì)誤差累計(jì)值，相對(duì)誤差累計(jì)值在一定程度上反映了Kc對(duì)計(jì)算徑流量造成的影響，相對(duì)誤差累計(jì)值越小越好；③計(jì)算19871989年三年相對(duì)誤差的均方差，它在一定程度上說明了各相對(duì)誤差偏離平均相對(duì)誤差的程度，也反映了Kc取值造成誤差的穩(wěn)定程度，相對(duì)誤差的均方差越小越好。雨量站分布較均勻，有10年逐日降水資料和時(shí)段降水資料；該水文站具有10年以上水位、流量資料；流域?qū)偕絽^(qū)性小流域且受到地形、地貌等下墊面條件影響，洪水陡漲緩落，匯流時(shí)間一般2—3h，有時(shí)更短；一次洪水總歷時(shí)2~5d。流域內(nèi)地勢(shì)平緩，土壤主要有黃壤和砂壤，具有明顯的腐殖層，淀積層和母質(zhì)土等層次結(jié)構(gòu)，透水性好。 定性分析該流域集雨面積856km2。7）IM:不透水面積占全流域面積的比例，它的值可由大比例尺的地形圖，通過地理信息系統(tǒng)現(xiàn)代技術(shù)量測(cè)出來，也可用歷史上干旱期小洪水資料來分析。 W[i+1]=WU[i+1]+WL[i+1]+WD[i+1]；式中：i表示第i天。 WD[i+1]=WD[i]ED[i]+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)。三層蓄水量變化的具體計(jì)算如下：1）WU[i]+P[i]EU[i]R[i]=UM，WU[i+1]=WU[i]+P[i]EU[i]R[i]。流域坡地上的降雨產(chǎn)流量因產(chǎn)流過程的條件和運(yùn)動(dòng)路徑不同，受流域的調(diào)蓄作用不同，各徑流成分在流量過程線上的反應(yīng)是不一樣的。蓄滿產(chǎn)流，即任一地點(diǎn)上，土壤含水量達(dá)田間持水量前，降雨量全部補(bǔ)充土壤含水量，不產(chǎn)流；當(dāng)土壤蓄滿后，其后續(xù)降雨量全部產(chǎn)生徑流。2)當(dāng)WU+PEP, WL≥CWLM, EU=WU+P,EL=(EPEU)WL/WLM,ED=0。(2) 流域內(nèi)有四個(gè)雨量站，、。流域上游有寶口水文站，流域面積553km2，%。暴雨成因主要是鋒面雨和臺(tái)風(fēng)雨，常受熱帶風(fēng)暴影響。（3）采用2004年暴雨數(shù)據(jù)進(jìn)行預(yù)報(bào)。根據(jù)已給的設(shè)計(jì)暴雨資料和任務(wù)一率定的Kc，將流域作為整體進(jìn)行如下計(jì)算：（1）次洪產(chǎn)流量計(jì)算，劃分水源；（2）直接徑流匯流，地下徑流匯流的計(jì)算。流域地處粵東沿海的西部，海