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