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