bocengsuanfa.zip_光子晶体_光子晶体滤波

clear % A为分母系数，B为分子系数。 A B for j=1:length(A) a(j)=A(length(A)+1-j); end for j=1:length(B) b(j)=B(length(B)+1-j); end n=length(A); C=zeros(n,n); D=zeros(n,n); R(n)=-(A(n)/B(1)); %最后一层的反射系数 R1(n)=-(B(n)/A(1)); C(n-1,:)=(R(n)/(1-R(n)^2))*b+(1/(1-R(n)^2))*A; %公式3.70 D(n-1,:)=(R(n)/(1-R(n)^2))*a+(1/(1-R(n)^2))*B; C1=zeros(n,n); D1=zeros(n,n); for j=1:n-1 C1(n-1,j)=C(n-1,n-j); %翻转操作 D1(n-1,j)=D(n-1,n-j); end R(n-1)=-(C(n-1,n-1)/D(n-1,1)); R1(n-1)=-(D(n-1,n-1)/C(n-1,1)); %倒数第二层反射系数 i=0; for k=n-2:-1:1 C(k,:)=(R(k+1)/(1-R(k+1)^2))*D1(k+1,:)+(1/(1-R(k+1)^2))*C(k+1,:); D(k,:)=(R(k+1)/(1-R(k+1)^2))*C1(k+1,:)+(1/(1-R(k+1)^2))*D(k+1,:); R1(k)=-(C(k,k)/D(k,1)); R(k)=-(D(k,k)/C(k,1)); %R为最终的反射系数 for j=1:k C1(k,j)=C(k,k+1-j); D1(k,j)=D(k,k+1-j); end end n=zeros(1,length(R)+1); n(1,1)=1; for i=2:length(R)+1 n(1,i)=(n(1,i-1)-R(i-1)*n(1,i-1))/(1+R(i-1)); end n %n为折射率 save E:\研究生课题\mat\D n