透射反射曲线（transmission line ）求解超材料介电常数 磁导率 阻抗 折射率

利用透射反射曲线求解超材料介电常数 磁导率 阻抗 折射率

load s11.txt
load s21.txt


d=4e-3;%样品厚度（m）
m=size（s11）;
n=m（11）;
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
r（p1）=s11（p2）*exp（-i*（s11（p3））*pi/180）;
t（p1）=s21（p2）*exp（-i*s21（p3）*pi/180）;
k=1e9*2*pi*s11（p1）/3e8;%频率单位（Hz）

index（p1）=acos（（1-r（p1）^2+t（p1）^2）/（2*t（p1）））/（k*d）;

if（imag（index（p1））<0）
index（p1）=-index（p1）;
end

%if （p>722）               
%   index（p1）=index（p1）+2*pi/k/d;
%end

impedance（p1）=sqrt（（（1+r（p1））^2-t（p1）^2）/（（1-r（p1））^2-t（p1）^2））;

if（real（impedance（p1））<0）
    impedance（p1）=-impedance（p1）;
end



mue（p1）=index（p1）*impedance（p1）;
epsilon（p1）=index（p1）/impedance（p1）;
M（p1）=s11（p1）;
M（p2）=real（mue（p1））;
M（p3）=imag（mue（p1））;
M（p4）=real（epsilon（p1））;
M（p5）=imag（epsilon（p1））;
M（p6）=real（index（p1））;
M（p7）=imag（index（p1））;
M（p8）=real（impedance（p1））;
M（p9）=imag（impedance（p1））;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      90094  2014-06-16 17:15  s21.txt

     文件        998  2014-06-16 16:26  fu.m

     文件      90094  2014-06-16 17:14  s11.txt

----------- ---------  ---------- -----  ----

               181186                    3