【正文】 ]aaa=[R’ Rs’ Rg’]。 Qs39。%Qs是算出來的逐日地表徑流過程endendendend %把地表徑流和地下徑流相加推求地下徑流for i=1:(m+n1)Q(i)=Qg(i)+Qs(i)。endendQs=zeros(1,m+n1)。%Qg是算出來的逐日地下徑流過程 end for i = m+1:(m+n1)Qg(i) = Cg * Qg(i 1)。 %出流系數法推求地下徑流Qg(1)=。 endenddwx=[0,40,80,130,100,80,48,20,10,5,0]。 RG(i) = 0。 RS(i) = R(i) RG(i)。% RG表示地下徑流深 RS(i) = 0。 a(i + 1) = WMM * (1 (1 W(i + 1) / WM) ^ (1 / (1 + B)))。if PE(i) = 0 %當不產流時 R(i) = 0。 else WD(i + 1) = DM。 else WL(i + 1) = LM。 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)。 end end %產流計算完畢if WU(i) + P(i) EU(i) R(i) = UM WU(i + 1) = WU(i) + P(i) EU(i) R(i)。 W(i + 1) = WM。 a(i + 1) = PE(i) + a(i)。 %流域產流計算凈降雨量PE%三層蒸散發(fā)計算，求總蒸發(fā)量和流域凈雨量if PE(i) 0 %當產流時 if PE(i) + a(i) WMM R(i) = PE(i) + W(i) WM + WM * (1 (PE(i) + a(i)) / WMM) ^ (B + 1)。 end endE(i) = EU(i) + EL(i) + ED(i)。 EL(i) = WL(i)。 ED(i) = 0。 WL(i) = C * (EP(i) EU(i)) EU(i) = WU(i) + P(i)。 ED(i) = 0。endif WU(i) + P(i) EP(i) if WL(i) = C * LM EU(i) = WU(i) + P(i)。EL(i) = 0。 for i = 1:m EP(i) = E0(i) * Kc。 %流域三層蒸發(fā)計算，先給出初始 上、下、深層土壤含水量WMM = WM * (1 + B)。 WL(1) = 40。 %逐日降雨深的計算end%求逐日降雨深 W(1) = 110。 %確定暴雨持續(xù)時間R=zeros(1,m)。Kc=。P3=hs(:,6)。P1=hs(:,4)。 %導入常數值。Fc=24。C=。DM=60。UM=20。)。 %DC是確定性系數end：%計算步驟%首先進行三層蒸發(fā)計算，得到凈降雨量，然后進行徑流劃分，然后計算次日到初始土壤含水量%然后根據劃分出的各時段的地表和地下徑流，求取洪水過程 %**************************************************************%% 2004年場次暴雨的洪水預報計算 %% %%**************************************************************%clcclear %清空數據 hs=load(39。Y=Y+(RR(i)Rr)^2。Y=0。for i=1:mRr=Rr+RR(i)/m。)。ylabel(39。天數39。endplot(s ,RR,s,R)。 a(i + 1) = WMM * (1 (1 W(i + 1) / WM) ^ (1 / (1 + B)))。if PE(i) = 0 %當不產流時 R(i) = 0。 else WD(i + 1) = DM。 else WL(i + 1) = LM。 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)。 end end %產流計算完畢if WU(i) + P(i) EU(i) R(i) = UM WU(i + 1) = WU(i) + P(i) EU(i) R(i)。 W(i + 1) = WM。 a(i + 1) = PE(i) + a(i)。 %流域產流計算凈降雨量PE%三層蒸散發(fā)計算，求總蒸發(fā)量和流域凈雨量if PE(i) 0 %當產流時 if PE(i) + a(i) WMM R(i) = PE(i) + W(i) WM + WM * (1 (PE(i) + a(i)) / WMM) ^ (B + 1)。 end endE(i) = EU(i) + EL(i) + ED(i)。 EL(i) = WL(i)。 ED(i) = 0。 WL(i) = C * (EP(i) EU(i)) EU(i) = WU(i) + P(i)。 ED(i) = 0。endif WU(i) + P(i) EP(i) if WL(i) = C * LM EU(i) = WU(i) + P(i)。EL(i) = 0。 for i = 1:m EP(i) = E0(i) * Kc。 %流域三層蒸發(fā)計算，先給出初始 上、下、深層土壤含水量WMM = WM * (1 + B)。 WL(1) = 40。 %逐日降雨深的計算end%求逐日降雨深 W(1) = 110。 EU=zeros(1,m)。 %確定Kc取值m=length(hs)。P4=hs(:,8)。P2=hs(:,6)。E0=hs(:,4)。 %導入常數值。Fc=24。C=。DM=60。UM=20。)。 %：此處僅以1987年為例%**************************************************************%% 1987年日徑流量理論值和實測值的檢驗 %% %%**************************************************************%%計算步驟%首先進行三層蒸發(fā)計算，得到凈降雨量，然后進行徑流劃分，然后計算次日到初始土壤含水量%然后根據劃分出的各時段的地表和地下徑流，求取洪水過程clcclear %清空數據 hs=load(39。 endend %，誤差均在5%以內[y,l]=min(cha)。 if cha(j)=KC(i)=kc(j)。 %直接篩選出1987年誤差最小的年份和差額i=1。)。ylabel(39。Kc值3