slab2.m

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

function  y=slab(x,t,d1,d2,n1,n2,n3)
B=x;

k1h=t*2*pi*sqrt(n1^2-B^2);
k2h=t*2*pi*sqrt(B^2-n2^2);
if n2>B
   k2h=t*2*pi*sqrt(n2^2-B^2);
end
k3h=t*2*pi*sqrt(B^2-n3^2);
if n3>B
   k3h=t*2*pi*sqrt(n3^2-B^2);
end
k2D2=d2/t*k2h;
k3D1=d1/t*k3h;
epss=1.e-5;
if abs(k2D2)<epss
   carp1=t/d2;
else
   carp1=k2h*coth(k2D2);
end
if abs(k3D1)<epss
   carp2=t/d1;
else
   carp2=k3h*coth(k3D1);
end
if (n2<=B & n3<=B)
y1=k1h^2*sin(2*k1h)-k1h*carp1*cos(2*k1h)...
   -carp2*(k1h*cos(2*k1h)+carp1*sin(2*k1h));
y=y1;
%y=real(y1);
else
 y1=k1h^2*sin(2*k1h)*sin(k2D2)*sin(k3D1)-k1h*k2h*cos(k2D2)*cos(2*k1h)*sin(k3D1)...
 -k3h*cos(k3D1)*(k1h*cos(2*k1h)*sin(k2D2)+k2h*cos(k2D2)*sin(2*k1h));
y=y1;
%y=real(y1);
end

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


 y1=k1h^2*sin(2*k1h)*sin(k2D2)*sin(k3D1)-k1h*k2h*cos(k2D2)*cos(2*k1h)*sin(k3D1)...
 -k3h*cos(k3D1)*(k1h*cos(2*k1h)*sin(k2D2)+k2h*cos(k2D2)*sin(2*k1h));