大气波导高度随各种参数变化

通过主函数的运行可以得到大气波导高度随各种气候参数（如气温、压强等）变化的趋势图

function Height=Delte（TakTskRHu）
  z1=6;      %测量参考高度  单位：m
  z0=0.00015;     %空气动力学中表面粗糙度参数  单位：m
  e=inline（‘6.105*exp（25.22*（T-273.2）/T-5.31*log（T/273.2））‘‘T‘）;
  Height=[];
  
  Rib=369*z1*（Tak-Tsk）/（u^2*Tak）;  %总体理查森数
  Gamma=PieceGamma（Rib）;     
  Lmo=（10*z1*Gamma）/Rib;     %Monin-Obukhov长度 
  ea=RH/100*e（Tak）;             %海洋表面周围空气水汽压
  es=e（Tsk）;                    %海洋表面海水水汽压
  
  Na=（77.6/Tak）*（1000+4810*ea/Tak）;    %海洋表面空气折射率
  Ns= （77.6/Tsk）*（1000+4810*es/Tsk）;    %海洋表面海水折射率
  Np=（Na-Ns）;

  if 0%       Height=Np/（-0.125*（log（z1/z0）+5.2*（z1-z0）/Lmo）-5.2*Np/Lmo）;
%       if Height>0||Height/Lmo>1
%          Height=（6.2*Np+0.65*z1）/（-0.125*log（z1/z0））;
%       end
      syms z;
      Height=solve（Np*（1+5.2*z/Lmo）/（z*（log（z1/z0）+5.2*（z1-z0）/Lmo））+0.125）;
      
 elseif Rib<0
      B=log（z1/z0）-PiecePsi（z1Lmo）;
      A=-0.125*B/Np;
      Height=（A^4-18*A^3/Lmo）^（-1/4）;
      
  elseif Rib==0
      Height = Np/（-0.125*log（z1/z0））;
  end
  if Height<0
      Height=[];
  end
      
  %Height=Height;
 
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       4389  2014-06-18 19:16  Main_H.m

     文件        237  2014-04-27 13:37  PieceGamma.m

     文件        390  2014-05-19 22:09  PiecePsi.m

     文件       1152  2014-06-18 18:41  Delte.m

----------- ---------  ---------- -----  ----

                 6168                    4