fr=S(:,1);%提取频率
s1=S(:,2)+S(:,3)*i;%参数 s1
s2=S(:,4)+S(:,5)*i;%参数 s2
Z_eff=sqrt(((1+s1).^2-s2.^2)./((1-s1).^2-s2.^2));%等效阻抗
gamma=(Z_eff-1)./(Z_eff+1);%反射系数
T=s2./(1-s1.*gamma);
c=3*10^8;%光速
d=140*10^(-3);%端口间的距离
k=fr/c*2*pi*10^6;%波数
m=0;
n_eff=-(log(abs(T))+i*(angle(T)+2*m*pi))./(k*d*i);%等效折射率
ep_eff=n_eff./Z_eff;%等效介电系数
mu_eff=n_eff.*Z_eff;%等效磁导率
figure(1);
plot(fr,real(mu_eff),'color',[0 0 1],'linewidth',2,'linestyle','-');
hold on;
plot(fr,imag(mu_eff),'color',[0 0 1],'linewidth',2,'linestyle','--');
hold off;
xlabel('Frequency (MHz)');ylabel('有效磁导率');
legend('Re(mu_eff)','Im(mu_eff)');
grid on;
兄弟能加个QQ1085694129交流下学习心得吗,也是刚接触超材料和HFSS。抱团取暖