alan.m

来自「ok &ouml nemli bir kaybnkao lasdhlas sdo」· M 代码 · 共 27 行

function  EY=alan(j,B,t,d1,d2,n1,n2,n3)

X=j;
k1h=t*2*pi*sqrt(n1^2-B^2);
k1x=X*2*pi*sqrt(n1^2-B^2);
k2h=t*2*pi*sqrt(B^2-n2^2);
k2x=X*2*pi*sqrt(B^2-n2^2);
k3h=t*2*pi*sqrt(B^2-n3^2);
k3x=X*2*pi*sqrt(B^2-n3^2);
k2D2=d2*2*pi*sqrt(B^2-n2^2);
k3D1=d1*2*pi*sqrt(B^2-n3^2);

f1h=k1h*cos(k1h-k1x)+k2h*coth(k2D2)*sin(k1h-k1x);
f2h=k1h*sinh(k2D2+k2h-k2x)/sinh(k2D2);
f3h=(k1h*cos(2*k1h)+k2h*sin(2*k1h)*coth(k2D2))*sinh(k3D1+k3h+k3x)/sinh(k3D1);
if  (X>=-t & X<=t)
   EY=f1h;
elseif (X>t & X<=(t+d2+0.25))
   EY=f2h;
elseif (X<-t & X>=(-t-d1-0.25))
   EY=f3h;
end

%  X=x/lambda
%  t=h/lambda
%  d1=D1/lambda
%  d2=D2/lambda