梯形断面正常水深莫洛图

% -- 计算梯形渠道的 正常水深计算
clc; clear all
h0=0.001:0.01:20;  % 无量纲参数的取值范围  β=αQ/b^2.5;
g=9.81; % 重力加速度
 
% plot(lam,entak0,'k-','linewidth',2); hold on
% 迭代一次
% plot(lam,entak1,'r-','linewidth',2); hold on
  
% 当 m=0 时； 矩形
beta=h0.^(5/3)./(1+2*h0).^(2/3) ;
plot(beta,h0,'k-','linewidth',2); hold on


% 当 m=0.25 时； 
m=0.25;
en=m*h0;
a1=(1+en).*en;
a2=(m+2*en.*sqrt(1+m^2));
a3=a1./a2.^0.4;
beta=1/m*a3.^(1/0.6);
plot(beta,h0,'b-','linewidth',2); hold on
 

% 当 m=0.5 时； 
m=0.5;
en=m*h0;
a1=(1+en).*en;
a2=(m+2*en.*sqrt(1+m^2));
a3=a1./a2.^0.4;
beta=1/m*a3.^(1/0.6);
plot(beta,h0,'k-','linewidth',2); hold on

% 当 m=1 时； 
m=1;
en=m*h0;
a1=(1+en).*en;
a2=(m+2*en.*sqrt(1+m^2));
a3=a1./a2.^0.4;
beta=1/m*a3.^(1/0.6);
plot(beta,h0,'b-','