temp1.m

用传输矩阵的方法来计算光子晶体带隙

a=1;
b=1;
c=1;
d=1;
i=sqrt(-1);
r=[0.1:0.1:10];
k=2*pi/2;
z=(a*besselh(1+1/2,1,k*r)+b*besselh(1+1/2,2,k*r)).*sqrt(pi/2/k./r);
P1=legendre(1,1/sqrt(2));
P1=P1(2);
P2=1/sqrt(2)*P1-2*legendre(0,1/sqrt(2));
M=1/sqrt(2)*(c*P1*z-i*d*P2*