左手材料参数

load S11.txt
load S21.txt
load arg_S11.txt
load arg_S21.txt
load f.txt;

d=0.00025;  %样品厚度
%m=size（S11）;
n=101;
r=zeros（n1）;
t=zeros（n1）;
index=zeros（n1）;
impedance=zeros（n1）;
mue=zeros（n1）;
epsilon=zeros（n1）;
M=zeros（n9）;


for p=1:n

%if（arg_S11（p1）<0）
%     arg_S11（p1）=arg_S11（p1）+2*pi;
%end

%if（arg_S21（p1）<0）
%     arg_S21（p1）=arg_S21（p1）+2*pi;
%end
r（p1）=S11（p1）*exp（-1i*arg_S11（p1））;
t（p1）=S21（p1）*exp（-1i*arg_S21（p1））;
%k=1e11*2*pi*S11（p1）/3e8;                    % 频率单位
k=2*pi*1e10*f（p1）/（3e8）;

impedance（p1）=sqrt（（（1+r（p1））^2-t（p1）^2）/（（1-r（p1））^2-t（p1）^2））;
if   real（impedance（p1））<0               % passive material 
   impedance（p1）=-impedance（p1）;
end



%e（p1）=t（p1）/（1-r（p1）*（impedance（p1）-1）/（impedance（p1）