马斯京根法matlab程序,（一）马斯京根流量演算法

(１)x值试算表：

经过对取不同的x值作图可知，当x＝０.１时，直线较为接近作直线求的K＝△s／△Q’=7.5，其中x是流量比重因子，k是河段发生稳定流的传播时间;根据求的x和k求的流量系数Co,C1,C2.

Co=(0.5△t-Kx)/(K-Kx+0.5△t); C1=(0.5△t+Kx)/(K-Kx+0.5△t);

C2=(K-0.5△t-Kx)/(K-Kx+0.5△t);求的Co=0.091; C1＝0.273; C2=0.636;

(2)流量演算表(包括整河段流量演算表和分河段流量演算表)

利用Q2= CoI2+ C1I1+

C2O1进行整河段流量演算

n=K/△t=2.5; Ll=L/n=44.8; l=(1-2x)L=89.6;

Xl=0.5-0.5*n(1-2x)=-0.5

查表得到Pm,n,进行分河段流量演算.

(3)将两种洪水过程线和实测洪水过程线;

(二)流域产汇流计算

(1)瞬时单位线和河段单位线计算表

地面径流Qd=Q-Qg;q1=Q1*10/h1;q2=(Q2-h2q1/10)*10/h1;

利用h=(∑q△t)/F对单位线修正.

Mh(1)、Mh(2)为地面净雨过程的一阶和二阶原点矩;

Mq(1)、Mq(2)为地面径流过程的一阶和二阶原点矩;

计算公式：Mh(1)= ∑hiti/∑hi;

Mh(2)= ∑hi(ti)/∑hi;

Mq(1)