fig161.m

利用matlab编程实现含负折射材料的表面波和波导特性

function [s]=fig161(e12,u12)
%clear
f1=10.5e9;
f2=11.5e9;
w1=2*pi*f1;
w2=2*pi*f2;
theta=30*pi/180;
c=3e8;
k01=w1/c;
k02=w2/c;
kx1=k01*sin(theta);
kx2=k02*sin(theta);
k0z1=k01*cos(theta);
k0z2=k02*cos(theta);
e11=2;u11=1;

if e12<0&u12<0
k1z1=-sqrt(w1^2*e11*u11/c^2-kx1^2);
k1z2=-sqrt(w2^2*e12*u12/c^2-kx2^2);
else
k1z1=sqrt(w1^2*e11*u11/c^2-kx1^2);
k1z2=sqrt(w2^2*e12*u12/c^2-kx2^2);
end
t=0;
w3=(w1+w2)/2;
lamda=2*pi*c/w3;
for k=1:400
    x(1)=-5*lamda;
    x(k+1)=x(k)+0.05*lamda;
    for  l=1:400
    z(1)=-5*lamda;
    z(l+1)=z(l)+0.05*lamda;
    if z(l)<0
    kai1(k,l)=kx1*x(k)+k0z1*z(l);
    kai2(k,l)=kx2*x(k)+k0z2*z(l);
    s(k,l)=kx1*cos(kai1(k,l))^2/(w1*1)+kx2*cos(kai2(k,l))^2/(w2*1)+(kx1/(w1*1)+kx2/(w2*1))*cos(kai1(k,l))*cos(kai2(k,l));
    else 
    kai1(k,l)=kx1*x(k)+k1z1*z(l);
    kai2(k,l)=kx2*x(k)+k1z2*z(l); 
    s(k,l)=kx1*cos(kai1(k,l))^2/(w1*u11)+kx2*cos(kai2(k,l))^2/(w2*u12)+(kx1/(w1*u11)+kx2/(w2*u12))*cos(kai1(k,l))*cos(kai2(k,l));
%     else  
%     kai1(k,l)=kx1*x(k)+k0z1*z(l);
%     kai2(k,l)=kx2*x(k)+k0z2*z(l);  
    end
    
   
    
    end
    
end

% hold on
% mesh (s(1:400,1:400)); figure(gcf);
% hold off