%----------------------------------------
%-----------------------To notice that the unit of s11a is degree.
%-----------------------len=0.012----------------len is equal to the length of the
sample
%-----------------------------------------
%sr=load('srrrrrrrrrrrrrrr.txt');sa=load('saaaaaaaaaaa.txt');f=sr(:,1);s11r=sr(:,2);s21r=sr(:,3);s11a
=sa(:,2);s21a=sa(:,3);
%------------------------------------------
s11=load('s11.txt');s21=load('s21.txt');f=s11(:,1);s11a=s11(:,3);s11r=s11(:,2);s21a=s21(:,3);s21r
=s21(:,2);
%---------------------------------------------
%s11r=load('s11r.txt');s11a=load('s11a.txt');s21r=load('s21r.txt');s21a=load('s21a.txt');f=s11r(:,
1);s11r=s11r(:,2);s21r=s21r(:,2);s11a=s11a(:,2);s21a=s21a(:,2);
%------------------------------------------------
af11=s11a*pi/180;
af21=s21a*pi/180;
s11=s11r.*(cos(af11)+i*sin(af11));
s21=s21r.*(cos(af21)+i*sin(af21));
s11=transpose(s11);
s21=transpose(s21);
f=transpose(f);
%---------to make sure the size of all matrix is size(1,length(f));
G=1e9;
c=3e8;
epr0=8.854188e-12;
len=0.005
L=length(s11r);
%----------------------------------------------
K=(s11.^2-s21.^2+1)./(2*s11);
gama11=K+(K.^2-1).^0.5;
gama12=K-(K.^2-1).^0.5;
Z11=(1+gama11)./(1-gama11);
Z12=(1+gama12)./(1-gama12);
for j=1:1:L
if real(Z11(j))>0%-- 14-11-2 ÏÂÎç8:53 --%
gama1(j)=gama11(j);
else
gama1(j)=gama12(j);
评论20