梯形断面临界水深莫洛图

% -- 计算梯形渠道的 临界水深计算
clc; clear all
beta=0.01:0.01:1000;  % 无量纲参数的取值范围  β=α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 时； 矩形
hk00=(1/g*beta.^2).^(1/3);
plot(beta,hk00,'k-','linewidth',2); hold on

% 当 m=0.25 时； 
m=0.25;
lam=m*(1/g*beta.^2).^(1/3);  % 无量纲参数的取值范围  λ=(αq^2/g)^1^/^3.m/b
a1=3.952*lam.*(1+4*lam).^0.2025; 
entak0=((1+a1).^0.5-1)/2; % 无量纲参数值
hk025=entak0/m;
plot(beta,hk025,'b-','linewidth',2); hold on
% 当 m=0.5 时； 
m=0.5;
lam=m*(1/g*beta.^2).^(1/3);  % 无量纲参数的取值范围  λ=(αq^2/g)^1^/^3.m/b
a1=3.952*lam.*(1+4*lam).^0.2025; 
entak0=((1+a1).^0.5-1)/2; % 无量纲参数值
hk05=entak0/m;
plot(beta,hk05,'k-','linewidth',2); hold on

% 当 m=1 时； 
m=1;
lam=m*(1/g*beta.^2).^(1/3);  % 无量纲参数的取