Comsol超材料S参数反演等效参数 负折射率超材料等效折射率、阻抗、介电常数与磁导率求解

Comsol超材料S参数反演等效参数 负折射率超材料等效折射率、阻抗、介电常数与磁导率求解

在实验室鼓捣超材料的时候，最让人头秃的就是等效参数反演。就像拆开收音机却看不懂电路板上的元件参数，明明测到了S参数，怎么倒推出材料的等效折射率这些核心指标？今儿咱们用COMSOL和MATLAB配合着实战一把。

Comsol超材料S参数反演等效参数 负折射率超材料等效折射率、阻抗、介电常数与磁导率求解

先打开COMSOL建个典型开口环结构，扫频范围设到8-12GHz。边界条件记得用周期性端口，这里有个坑：波端口要离结构至少1/4波长，不然近场干扰会让S参数飘得亲妈都不认识。跑完仿真直接导出S11和S21的实部虚部，存成txt备用。

拿到原始数据先做个预处理：

raw_data = importdata('s_params.txt'); freq = raw_data(:,1); % 第一列是频率 S11 = raw_data(:,2) + 1i*raw_data(:,3); S21 = raw_data(:,4) + 1i*raw_data(:,5);

这里有个骚操作——用传输矩阵法处理相位模糊问题。传统NRW方法在谐振区会翻车，咱们改进下：

Z = sqrt( ((1+S11).^2 - S21.^2)./((1-S11).^2 - S21.^2) ); % 阻抗计算 n = 1/(1i*k*d) .* log( (S21./(1 - S11.*(Z-1)./(Z+1))) ); % 折射率计算

注意这里的波数k=2pifreq/3e8，d是超材料厚度。算出来的n可能会出现负值的实部，这就是负折射率的体现。但直接这么算会有跳点，得加个相位展开：

n_real = real(n); n_imag = imag(n); n_unwrap = unwrap(angle(exp(1i*(n_real + 1i*n_imag)))); % 相位解缠绕 n_corrected = abs(n).*exp(1i*n_unwrap);

处理完的数据画个图验证下，正常应该看到在10GHz附近实部穿到负区间。这时候别急着高兴，记得用K-K关系验证因果性：

[eps_real, mu_real] = deal(real(epsilon), real(mu)); if any(eps_real <0 & mu_real <0) warning('双负区域出现，检查是否真实负折射') end

最后导出参数时，建议用相对值格式保存，避免科学计数法丢失精度。碰到阻抗突变点别慌，可能是Fabry-Perot谐振，这时候需要结合场分布图判断。等效参数反演就像拼乐高，数据清洗占七成功夫，剩下三成靠对物理机制的直觉把握。