抛物线断面临界水深莫洛图

% -- 计算  抛物线 形渠道的 临界水深计算
clc; clear all
beta=0.01:0.01:1000;  % 无量纲参数的取值范围  β=αQ/b^2.5;
g=9.81; % 重力加速度
 Q=beta*g^0.5; 
 
 % 当 m=0.5 时； 
m=0.25;
hk025=(27/64*beta.^2/m).^0.25; % 单位 m；
plot(Q,hk025,'b-','linewidth',2); hold on
 
% 当 m=0.5 时； 
m=0.5;
hk05=(27/64*beta.^2/m).^0.25; % 单位 m；
plot(Q,hk05,'k-','linewidth',2); hold on
 
 % 当 m=1 时； 
m=1;
hk10=(27/64*beta.^2/m).^0.25; % 单位 m；
plot(Q,hk10,'b-','linewidth',2); hold on

 % 当 m=1.5 时； 
m=1.5;
hk15=(27/64*beta.^2/m).^0.25; % 单位 m；
plot(Q,hk15,'k-','linewidth',2); hold on

 % 当 m=2 时； 
m=2;
hk20=(27/64*beta.^2/m).^0.25; % 单位 m；
plot(Q,hk20,'b-','linewidth',2); hold on

 % 当 m=2.5 时； 
m=2.5;
hk25=(27/64*beta.^2/m).^0.25; % 单位 m；
plot(Q,hk25,'k-','linewidth',2); hold on
 
 % 当 m=3 时； 
m=3;
hk30=(27/64*beta.^2/m).^0.25; % 单位 m；
plot(Q,hk30,'b-','linewidth',2); hold on

 % 当 m=4 时； 
m=4;
hk40=(27/64*beta.^2/m).^0.25; % 单位 m；
plot(Q,hk40,'k-','linewidth',2); hold on

  axis([0.01,2000,0.1,20])
 

xlabel('流量  Q（m^3/s）');
 ylabel('临界水深 h_k (m) ');
  
% legend('地基加固后2#桩检测值','本文预估计算值（1#桩）', '原状土1#桩检测值','本文计算值（1#桩）') 


text(103,0.9,'p=0.25')
text(103,0.8,'p=0.5')
text(103,0.7,'p=1')
text(103,0.6,'p=1.5')
text(103,0.5,'p=2')
text(103,0.4,'p=2.5')
text(103,0.3,'p=3')
text(103,0.2,'p=4')

title('抛物线形渠道 x^2=2py  临界水深查算图');