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）