【正文】 相對誤差隨著Kc值變化圖xlabel(39。天數(shù)39。)。ylabel(39。流量值39。)。 %坐標軸表示對象標簽Rr=0。for i=1:mRr=Rr+RR(i)/m。 %Rr是實測值的均值endX=0。Y=0。for i=1:mX=X+(R(i)RR(i))^2。Y=Y+(RR(i)Rr)^2。DC=1X/Y。 %DC是確定性系數(shù)end：%計算步驟%首先進行三層蒸發(fā)計算，得到凈降雨量，然后進行徑流劃分，然后計算次日到初始土壤含水量%然后根據(jù)劃分出的各時段的地表和地下徑流，求取洪水過程 %**************************************************************%% 2004年場次暴雨的洪水預報計算 %% %%**************************************************************%clcclear %清空數(shù)據(jù) hs=load(39。39。)。WM=140。UM=20。LM=60。DM=60。B=。C=。 IM=。Fc=24。Cg=。 %導入常數(shù)值。E0=hs(:,3)。P1=hs(:,4)。P2=hs(:,5)。P3=hs(:,6)。P4=hs(:,7)。Kc=。 %確定Kc取值m=length(hs)。 %確定暴雨持續(xù)時間R=zeros(1,m)。for i = 1:m %流域平均降雨量計算P(i) = * P1(i) + * P2(i) + * P3(i) + * P4(i)。 %逐日降雨深的計算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 i = 1:m EP(i) = E0(i) * Kc。%逐日蒸散發(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) + 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 WD(i + 1) = DM。 end endend%此處計算的是降雨補充土壤含水量，由前一天的土壤含水量推求第二天的土壤含水量，補充來源為降雨減去蒸散發(fā)減去徑流量，順序為上、下、深依此補充。if PE(i) = 0 %當不產(chǎn)流時 R(i) = 0。 W(i + 1) = W(i) + PE(i)。 a(i + 1) = WMM * (1 (1 W(i + 1) / WM) ^ (1 / (1 + B)))。end end%以上產(chǎn)流計算完畢f(xié)or i = 1:m %水源劃分 if PE(i) 0 if PE(i) = Fc/8 % Fc穩(wěn)定下滲率 RG(i) = R(i)。% RG表示地下徑流深 RS(i) = 0。 %RS表示地表徑流深 else RG(i) = Fc /8* R(i) / PE(i)。 RS(i) = R(i) RG(i)。 end end if PE(i) = 0 R(i) = 0。 RG(i) = 0。 RS(i) = 0。 endenddwx=[0,40,80,130,100,80,48,20,10,5,0]。n=length(dwx)。 %出流系數(shù)法推求地下徑流Qg(1)=。 for i = 2:m Qg(i) = Cg * Qg(i 1) + (1 Cg) * RG(i) * 553 / (3 * )。%Qg是算出來的逐日地下徑流過程 end for i = m+1:(m+n1)Qg(i) = Cg * Qg(i 1)。end%單位線法求地表徑流for i=1:m for j=1:nqs(i,j)=RS(i)*dwx(j)/10。endendQs=zeros(1,m+n1)。for k=1:(m+n1)for i=1:m for j=1:nif k==i+j1Qs(k)=Qs(k)+qs(i,j)。%Qs是算出來的逐日地表徑流過程endendendend %把地表徑流和地下徑流相加推求地下徑流for i=1:(m+n1)Q(i)=Qg(i)+Qs(i)。%Q表示逐日總的徑流量endhsgc=[Qg39。 Qs39。 Q39。]aaa=[R’ Rs’ Rg’]