简介面向超材料研究的CST参数反演MATLAB脚本资源适合电磁仿真工程师和科研人员使用重点解决从S参数实验数据逆向推导超材料结构参数的问题。压缩包内仅有1个m文件大小约1KB为get_S_Parameter.m脚本可加载S参数数据、设定反演目标函数并调用优化算法输出与实测数据匹配的结构参数同时支持与CST Studio Suite联动实现仿真流程自动化。资源围绕S参数散射特性、参数反演原理、MáSMATLAB-CST接口等核心知识展开涵盖超材料单元建模、电场磁场分布计算、谐振频率与阻抗匹配分析等关键环节帮助读者建立从CST仿真到材料参数提取的完整技术链条。借助该脚本可减少手动反复调试结构尺寸的工作量加速超材料单元优化与实验验证进程尤其适合周期性结构设计、负折射与波束控制等研究场景。已有940人学习适合正在开展超材料设计与验证、需要快速搭建反演流程的中高级研究者参考。1. 超材料S参数反演是什么从CST仿真到等效材料的最后一公里CST超材料仿真最容易被低估的一步不是建模也不是扫参而是仿真结束后怎么把模型变成等效介电常数和等效磁导率。很多人在CST里跑完单元结构盯着S11、S21两根曲线只能说出“这里有个谐振”但说不出这个单元在某个频段到底是负折射、磁谐振还是电谐振更没法把结果交给下一步电路设计。get_S_Parameter负责把复数S参数正确取出来超材料CST参数反演则用NRW这类算法把S参数换算成等效电磁参数两者合起来才是一条完整链路。这篇文章面向做吸波体、频率选择表面、超表面透镜和天线罩周期结构的人从CST仿真设置讲到反演代码再把最常见的坑逐个点名。2. CST超材料仿真建模单元尺寸、unit cell边界与S参数获取2.1 单元模型设置周期、基板与频率范围的先决参数超材料单元建模的第一个决定是周期尺寸。周期不是随便给的它决定了工作频带内会不会出现高阶Floquet模。经验法则是单元周期P小于自由空间波长的二分之一在正入射条件下高阶模截止S参数才是单模的反演公式才成立。以10 GHz附近工作为例周期取4 mm比较稳妥如果目标是20 GHz周期要缩到2 mm左右。建模前先把周期、基板厚度、金属层厚度这三个量定死后续所有扫参都围绕它们展开。基板材料的选择直接影响反演结果的损耗虚部。我常用RO4350B介电常数3.48、损耗角正切0.0037、厚度0.508 mm铜厚0.035 mm。金属谐振结构用Lossy Metal建模线宽不要小于0.1 mm太细的线在频域求解器里会显著增加网格量而且实际加工也做不出来。频率范围建议比目标工作频带外扩30%比如做8到12 GHz反演时远离谐振的频点能用来约束折射率分支这是后面分支选择的关键。求解器推荐用频域求解器Frequency Domain Solver配合同一频带内的离散扫频。时域求解器也能跑周期结构但谐振结构Q值高时时域信号衰减慢收敛时间会拉长频域求解器在每个频点单独求解对窄带高Q结构更稳。CST里单元结构的仿真频率点数用自适应扫频时不用手动设太多后处理时再插值加密。设置参数可参考下面的表参数推荐值说明单元周期P λ₀/2防止高阶Floquet模在工作频带内传播基板厚度0.5 ~ 1.0 mm太薄反演厚度误差敏感太厚等效媒质假设失效铜厚0.035 mm1oz覆铜板标准值频率范围目标频带外扩30%给分支约束留出非谐振区数据求解器Frequency Domain高Q结构比时域收敛更可控2.2 unit cell边界与Floquet端口让S参数在周期结构下成立超材料仿真必须用周期性边界这是很多新手第一次翻车的地方。如果X、Y方向用理想电壁和理想磁壁组合只能模拟特定极化、特定入射角的对称模式斜入射和交叉极化信息全丢。正确做法是X、Y方向都设unit cell边界Z方向设Floquet端口。正入射时Floquet端口只需要Mode 1斜入射时需要把端口模式数开到足够覆盖入射角否则S参数不完整。Floquet端口到样品表面的距离也要留够。端口面若紧贴基板近场高次模会影响S参数相位我一般让端口面距样品表面至少4到5 mm再在CST端口设置里用Deembed功能把参考面移到样品表面。这一步特别重要它直接决定S21相位的零点位置。反演时如果发现折射率曲线整体带有线性漂移先回来检查Deembed距离而不是怀疑算法有问题。端口极化方向也需要确认。Floquet端口默认有两个正交极化模式正入射计算时只激励一个极化另一个极化模式的结果如果不使用就不要在频带内开启它。开启多余的模式会引入虚假的耦合路径导致S11曲线出现不应有的毛刺。仿真完成后正确的结果是S11在谐振点附近深陷远离谐振时接近0 dBS21远离谐振时接近0 dB且相位随频率近似线性变化。看到这样的趋势再往下走才有意义。2.3 get_S_Parameter从后处理结果里把复数S11/S21导出来仿真结束后导航树里进入1D Results展开S-Parameters就能看到S1,1和S2,1。界面上默认显示幅度dB要看相位就切到Phase标签。想同时读两个频率点的值按住Shift键用鼠标在曲线上点两个Marker两个点的频率和幅值会同时显示在图的右上角这个操作在调谐振点位置时非常顺手。真正给反演用的是复数S参数不是图上那根曲线。导出方式有两种一种是在导航树里右键选中S参数Export Text导出带频率、实部、虚部的文件另一种更推荐在后处理模板Result Templates里创建Touchstone导出任务让CST把S参数写成标准的.s2p文件。.s2p文件第一行会注明格式是RI实部/虚部还是MA幅度/相位反演脚本读取时必须区分这是第3章的重点。用cst扫参做完参数扫描后每个样本的S参数会按样本名分好组通过Result Templates批量导出每个样本一个.s2p文件。文件名里带上扫描参数值比如P4.0_d0.5.s2p后续批处理脚本一眼就能看出几何参数。导出前确认端口阻抗是50欧姆频率单位是GHz这两个信息写在.s2p文件头里脚本读取时要处理单位换算否则算k0时频率差三个数量级反演结果完全不可用。3. S参数预处理从Touchstone文件到干净的反演输入3.1 Touchstone文件里有什么格式判断与频点对齐CST导出的.s2p文件本质是纯文本表格每行一个频点后面依次是S11实部、S11虚部、S21实部、S21虚部。文件头若干行以感叹号开头的注释紧接着是# GHz S RI R 50这样的格式行。这里的RI表示数据是实部虚部格式如果是# GHz S MA R 50则表示幅度相位格式。很多反演脚本读MA格式时当RI解析结果S参数变成一堆毫无意义的数反演出来全是噪声。还有一个容易忽略的点是频点对齐。扫参时不同几何参数的样本可能设置了不同的频点数量或者用了自适应扫频导致频点不均匀。反演前必须把所有样本的S参数插值到同一组频点上否则后面批处理时数组长度不一致代码直接报错。CST自适应扫频的频点在谐振区更密插值到均匀频点后谐振细节会丢失我一般保留原频点做反演只有需要对比多个样本时才统一插值。3.2 用Python解析S2P并按需去嵌给出一段能直接用的Python读取函数兼容RI和MA两种格式import numpy as np def load_touchstone(path): 读取CST导出的.s2p返回频率(Hz)、S11、S21复数数组 f, s11_re, s11_im, s21_re, s21_im [], [], [], [], [] is_mag_angle False with open(path, encodingutf-8, errorsignore) as fh: for line in fh: line line.strip() if not line or line.startswith(!): continue if line.startswith(#): if MA in line.upper(): is_mag_angle True continue parts line.split() if len(parts) 5: continue f.append(float(parts[0])) if is_mag_angle: # 幅度(dB或线性)/相位(度) mag 10 ** (float(parts[1]) / 20) ang np.deg2rad(float(parts[2])) s11_re.append(mag * np.cos(ang)) s11_im.append(mag * np.sin(ang)) mag 10 ** (float(parts[3]) / 20) ang np.deg2rad(float(parts[4])) s21_re.append(mag * np.cos(ang)) s21_im.append(mag * np.sin(ang)) else: # RI格式直接读实部虚部 s11_re.append(float(parts[1])) s11_im.append(float(parts[2])) s21_re.append(float(parts[3])) s21_im.append(float(parts[4])) f np.asarray(f) * 1e9 # GHz - Hz S11 np.asarray(s11_re) 1j * np.asarray(s11_im) S21 np.asarray(s21_re) 1j * np.asarray(s21_im) return f, S11, S21读文件时先把格式行的关键字MA或RI抓出来再决定每列怎么解析。频率单位统一换算成Hz这一步不及时做后面计算k0时全乱。CST如果导出的是线性幅度加角度注意幅度列若是dB要转回线性值很多CST版本的Touchstone幅度列默认写dB直接当线性值用会导致|S21|远超1反演出的折射率虚部完全失去物理意义。预处理还应该包含可选的端口去嵌。如果CST里没做Deembed端口参考面仍在端口物理位置S参数相位里含有一段空气传播路径。去嵌公式是S11乘上exp(2jk₀L)S21乘上exp(jk₀L)L是端口到样品表面的空气距离。这段运算放在预处理里既不用重新仿一遍又能多次微调L值看相位曲线的敏感性。3.3 预处理阶段最容易埋雷的三个细节第一个细节是相位解包裹的时机。反演需要连续的S21相位但atan函数返回的主值区间是[-π, π]相位跨过±180°时会跳变。Python的np.angle得到的就是跳变相位必须用np.unwrap解包裹才能用于反演。解包裹放在反演脚本里做预处理阶段只保证数据格式正确。第二个细节是频点密度。谐振区的S21相位变化非常快如果频点间隔过大相邻频点相位差超过π解包裹算法会判断错误反演出错误分支。做8到12 GHz仿真时我习惯至少取161个频点谐振附近相位变化再剧烈也能覆盖住。CST频域求解器的离散扫频里直接指定频点数量不要用默认的稀疏频点。第三个细节是保留原始文件。预处理脚本每次从.s2p原始文件读入不要覆盖式修改也不要只保存处理后的中间结果。分支选择不对时你大概率要回到原始S参数重新看相位原始文件一旦被改过前面几个小时的排查都白费。4. 参数反演核心算法NRW公式、分支处理与最小实现4.1 从S到ε/μ的四步换算与公式含义标题里的MáS泛指这类“CST出S参数、脚本做反演”的自动化方案它不是非用不可的黑匣子核心就是Nicolson-Ross-Weir方法也就是常说的NRW反演。NRW把超材料单元等效成一块均匀介质板用S11和S21反推这层“等效介质”的折射率n和波阻抗z再通过ε n/z、μ n·z得到等效介电常数和等效磁导率。换算过程分四步。第一步定义中间量V1 S21 S11、V2 S21 - S11进而算出X (1 - V1·V2)/(V1 - V2)。第二步用X求反射系数Γ X ± √(X² - 1)根号外符号的选择原则是保证|Γ| ≤ 1。第三步算传输系数T (V1 - Γ)/(1 - V1·Γ)然后由Γ求归一化阻抗z (1 Γ)/(1 - Γ)。第四步用传播相位求折射率n其中关键是要对S21的相位解包裹并处理对数函数的多值性问题。这套公式的物理前提是均匀化。单元周期远小于波长、样品厚度内场分布近似均匀时反演结果才有意义。基板厚度0.508 mm、工作频率10 GHz时等效介质厚度只有约λ₀/60均匀化误差可以接受如果单元电尺寸大到接近半个波长反演出的ε、μ会随厚度剧烈震荡那时要回头缩小周期而不是继续调算法。4.2 一份能直接跑的反演脚本下面是完整的NRW反演函数输入频率、S11、S21和样品厚度输出ε、μ、n、z四条曲线import numpy as np def nrw_retrieval(freq, S11, S21, d): NRW参数反演适用于CST导出的正入射S参数 freq: 频率, Hz S11, S21: 复数S参数数组 d: 样品厚度, m 返回: freq, eps_r, mu_r, n, z S11 np.asarray(S11, dtypecomplex) S21 np.asarray(S21, dtypecomplex) V1 S21 S11 V2 S21 - S11 denom V1 - V2 denom[abs(denom) 1e-12] 1e-12 # 防止除零 X (1 - V1 * V2) / denom # 两个根都算出来选|Gamma|1的作为物理根 G1 X np.sqrt(X**2 - 1) G2 X - np.sqrt(X**2 - 1) Gamma np.where(np.abs(G1) 1.0 1e-12, G1, G2) T (V1 - Gamma) / (1 - V1 * Gamma) # 阻抗和折射率 z (1 Gamma) / (1 - Gamma) phase np.unwrap(np.angle(T)) # 低频相位基线校准假定低频趋近普通介质而非负折射 phase phase - phase[0] t_abs np.abs(T) k0 2 * np.pi * freq / 3e8 n (-phase 1j * np.log(t_abs)) / (k0 * d) eps n / z mu n * z return freq, eps, mu, n, z这段代码的逻辑是先构造V1、V2让X变成一个中间量对Gamma的两个根用np.where按模长选择这一步直接处理了分支选择中最常见的问题T算出来后np.unwrap把S21相位展平再减去第一个频点的相位作为基线让折射率曲线在低频端回到正常介质的位置。模长部分取对数得到衰减项对应折射率虚部。参数上需要重点说明两点。第一厚度d是基板厚度而不是CST里unit cell的Z向总高度很多超材料单元在基板上下还有空气层如果直接把总高度拿来用n会系统性偏大。第二CST的时谐因子采用e^{jωt}无源材料的折射率虚部反演出来是负值对应的ε″、μ″也是负值这是约定问题不是错误和文献里用e^{-iωt}的符号习惯相反时肉眼看着像“损耗为负”实际上并没有违反无源判据。4.3 厚度、相位基线与分支选择反演可靠性的三个闸门反演结果可不可信第一个闸门是厚度。NRW公式里n与kd成反比厚度差10%n就差10%谐振频率位置也会跟着漂移。厚度取值我按基板厚度算铜层和印刷结构对厚度的贡献在10 GHz频段可以忽略。如果结构是金属-介质-金属三层厚度取两层介质加中间层的总厚度不能只取单层基板。第二个闸门是相位基线。phase - phase[0]这个操作把低频相位归零隐含假设是第一个频点远离谐振且相速接近自由空间。如果仿真频带从8 GHz起而谐振点在8.2 GHz基线校准就会失真。稳妥做法是让最低仿真频率比谐振频率低至少20%给基线校准留出“正常介质区”。第三个闸门是n的连续性。解包裹能解决2π跳变但NRW在某些频点仍会误入非物理分支表现是n在某一点突然跳变到另一个连续曲线。处理方式是先看n曲线是否平滑如果只有一个孤立的跳变点手动修正计算该点前后n的差值若差值超过k₀d分之π就整体减去或加上2π/(k₀d)。如果全频带都在无规律震荡多半是S参数本身的问题回第3章查相位解包裹和频点密度。5. S参数反演避坑五个高频问题的现象、原因与对策5.1 S21相位锯齿跳变解包裹后依然不平滑现象反演出的ε曲线呈高频锯齿状完全没有文献里那种平滑趋势n曲线上还有周期性尖峰。原因最常见的两个来源。一是频点太粗谐振区相邻频点相位差超过πnp.unwrap把真实变化误判成跳变二是.s2p文件幅度列是dB读取时没转线性|T|出现大于1的假象折射率虚部跟着震荡。解决把CST扫频点数加到160以上确认Touchstone格式行的MA/RI标记幅度用10**(dB/20)转线性值。改完重跑反演锯齿通常直接消失。5.2 反演出的损耗虚部符号和文献相反现象介质基板明明有损耗反演的ε″却是负值μ″在谐振区也是负值看着像材料的能量在增加。原因时谐因子约定。CST材料库和场求解器用e^{jωt}在这个约定下无源介质的介电常数虚部为负很多论文和教科书用e^{-iωt}虚部为正。两边符号习惯不同反演结果自然“相反”。这不是算法错误也不是材料异常。解决用能量判据代替符号判断。无源材料要求Re(n)0且Im(n)0CST约定下。如果反演结果的实部为正、谐振形态正常只是虚部符号和某篇论文相反按目标平台的约定统一取共轭即可不要为了“凑正号”乱改分支。5.3 厚度d取值不同结果差出一倍现象同一组S参数d取0.5 mm和1.0 mm低频折射率差一倍谐振频率位置也明显偏移。原因NRW的n与k₀d成反比d是全局比例因子。d取错n、n″同步缩放ε、μ也跟着偏。常见的取错方式是直接把CST里unit cell在Z方向的总高度当厚度这个值往往包含了上下空气层比实际样品厚度大好几倍。解决以基板物理厚度为准金属层在10 GHz以下厚度贡献小于1%可忽略。定义成脚本常量后不要每次手改多个样本对比时用同一个d值。5.4 曲线在某个频点“劈叉”或突变现象ε曲线整体平滑突然在某一个或几个频点跳到另一条曲线上回跳后趋势接不上S11模值在这些频点恰好接近0或接近1。原因NRW在|Γ|→0阻抗接近自由空间或|T|→0全反射时公式病态。S11接近0时X的计算变成0/0型的数值不稳定T接近0时log(T)的模长趋于负无穷虚部失去物理意义。另一个可能是该频点附近高阶Floquet模开始传播S参数不再只包含主模信息。解决先看S11模值确认异常频点是否处于匹配状态若是匹配点手动剔除该频点附近的几个数据再插值。若异常出现在高频端且单元周期接近半波长缩小周期尺寸后重新仿真。5.5 扫参到某些尺寸时仿真发散现象参数扫描到某个几何尺寸组合时频域求解器提示不收敛或发散S参数曲线在部分频点突然出现毛刺反演结果完全不可用。原因谐振结构出现极细的缝隙或尖锐拐角网格剖分不够高Q谐振在频域求解器里需要更多迭代次数薄金属片厚度小于网格步长时边缘场分布失真。解决频域求解器的自适应网格属性打开收敛目标设到-40 dB对薄金属结构改用Lossy Metal的thin sheet模型而不是按实体厚度剖分。毛刺若只出现在个别频点把这些频点从反演数据里剔除不要硬着头皮解包裹。6. 把反演流程自动化批处理、结果验证与表面参数进阶6.1 批量导出与自动反演参数扫描跑完CST的Result Templates已经为每个样本导出了.s2p文件剩下的就是让Python脚本批量处理。我一般用文件名解析几何参数再逐个调用前面三个函数读取、预处理、反演结果存成.npz文件import glob import numpy as np for path in sorted(glob.glob(./s2p/*.s2p)): freq, S11, S21 load_touchstone(path) freq, S11, S21 preprocess(freq, S11, S21, d_port0.0) freq, eps, mu, n, z nrw_retrieval(freq, S11, S21, d0.000508) tag path.split(/)[-1].replace(.s2p, ) np.savez(f./out/{tag}.npz, freqfreq, epseps, mumu, nn, zz) print(f{tag}: Re(eps) min{eps.real.min():.2f})这个脚本的输出是每个样本一个.npz文件后续只读.npz画图或做优化不再重复解析原始文件。打印语句输出每个样本的介电常数最小值扫参时一眼就能看出哪个尺寸组合出现了负折射行为。6.2 用等效参数回代验证反演结果不能只看曲线好不好看必须回代验证。做法是在CST里新建一个均匀介质板厚度与反演时用的d一致材料色散用反演得到的ε(f)、μ(f)表格导入边界和端口设置保持不变跑一次透射仿真。把回代的S21幅度和相位与原始超材料单元的S21放在同一张图里对比谐振点幅度偏差小于1 dB、相位趋势一致说明均匀化假设成立反演可信。对比项原始单元S21回代等效板S21谐振频率9.82 GHz9.78 GHz谐振深度-18.3 dB-17.9 dB带外相位斜率-2.1 rad/GHz-2.2 rad/GHz偏差大时先查厚度取值再查相位基线最后才怀疑公式本身。回代这一步也帮我判断过好几次分支选错曲线形态像谐振位置对不上基本都是n分支整体偏移了一个2π/(k₀d)的常数。6.3 薄超表面优先用表面参数如果单元厚度远小于工作波长比如d小于λ₀/20体积等效参数会变得对厚度极其敏感ε、μ的数值大得离谱物理意义弱化。这时更适合提取表面参数比如表面阻抗Z_s或表面极化率公式从传输线模型出发把超表面当作一个零厚度的阻抗边界S参数直接换算成Z_s。判断标准很简单厚度方向上的场分布如果不是近似均匀就别硬用NRW。把手头数据先做一遍回代验证验证不过就切到表面参数路线这是效率最高的工作顺序。我自己每次反演前固定检查三件事S21相位解包裹是否干净、S11模值有没有接近0或1的频点、厚度d是不是写成了常量而不是从文件名里再解析一遍。三次检查全过反演结果才敢往报告里放。希望帮到你。本文还有配套的精品资源点击获取 
长河君 · 2026/9/23 22:04:09 · 阅读 1.2k CST超材料仿真S参数反演：从S11/S21到等效电磁参数
简介面向超材料研究的CST参数反演MATLAB脚本资源适合电磁仿真工程师和科研人员使用重点解决从S参数实验数据逆向推导超材料结构参数的问题。压缩包内仅有1个m文件大小约1KB为get_S_Parameter.m脚本可加载S参数数据、设定反演目标函数并调用优化算法输出与实测数据匹配的结构参数同时支持与CST Studio Suite联动实现仿真流程自动化。资源围绕S参数散射特性、参数反演原理、MáSMATLAB-CST接口等核心知识展开涵盖超材料单元建模、电场磁场分布计算、谐振频率与阻抗匹配分析等关键环节帮助读者建立从CST仿真到材料参数提取的完整技术链条。借助该脚本可减少手动反复调试结构尺寸的工作量加速超材料单元优化与实验验证进程尤其适合周期性结构设计、负折射与波束控制等研究场景。已有940人学习适合正在开展超材料设计与验证、需要快速搭建反演流程的中高级研究者参考。1. 超材料S参数反演是什么从CST仿真到等效材料的最后一公里CST超材料仿真最容易被低估的一步不是建模也不是扫参而是仿真结束后怎么把模型变成等效介电常数和等效磁导率。很多人在CST里跑完单元结构盯着S11、S21两根曲线只能说出“这里有个谐振”但说不出这个单元在某个频段到底是负折射、磁谐振还是电谐振更没法把结果交给下一步电路设计。get_S_Parameter负责把复数S参数正确取出来超材料CST参数反演则用NRW这类算法把S参数换算成等效电磁参数两者合起来才是一条完整链路。这篇文章面向做吸波体、频率选择表面、超表面透镜和天线罩周期结构的人从CST仿真设置讲到反演代码再把最常见的坑逐个点名。2. CST超材料仿真建模单元尺寸、unit cell边界与S参数获取2.1 单元模型设置周期、基板与频率范围的先决参数超材料单元建模的第一个决定是周期尺寸。周期不是随便给的它决定了工作频带内会不会出现高阶Floquet模。经验法则是单元周期P小于自由空间波长的二分之一在正入射条件下高阶模截止S参数才是单模的反演公式才成立。以10 GHz附近工作为例周期取4 mm比较稳妥如果目标是20 GHz周期要缩到2 mm左右。建模前先把周期、基板厚度、金属层厚度这三个量定死后续所有扫参都围绕它们展开。基板材料的选择直接影响反演结果的损耗虚部。我常用RO4350B介电常数3.48、损耗角正切0.0037、厚度0.508 mm铜厚0.035 mm。金属谐振结构用Lossy Metal建模线宽不要小于0.1 mm太细的线在频域求解器里会显著增加网格量而且实际加工也做不出来。频率范围建议比目标工作频带外扩30%比如做8到12 GHz反演时远离谐振的频点能用来约束折射率分支这是后面分支选择的关键。求解器推荐用频域求解器Frequency Domain Solver配合同一频带内的离散扫频。时域求解器也能跑周期结构但谐振结构Q值高时时域信号衰减慢收敛时间会拉长频域求解器在每个频点单独求解对窄带高Q结构更稳。CST里单元结构的仿真频率点数用自适应扫频时不用手动设太多后处理时再插值加密。设置参数可参考下面的表参数推荐值说明单元周期P λ₀/2防止高阶Floquet模在工作频带内传播基板厚度0.5 ~ 1.0 mm太薄反演厚度误差敏感太厚等效媒质假设失效铜厚0.035 mm1oz覆铜板标准值频率范围目标频带外扩30%给分支约束留出非谐振区数据求解器Frequency Domain高Q结构比时域收敛更可控2.2 unit cell边界与Floquet端口让S参数在周期结构下成立超材料仿真必须用周期性边界这是很多新手第一次翻车的地方。如果X、Y方向用理想电壁和理想磁壁组合只能模拟特定极化、特定入射角的对称模式斜入射和交叉极化信息全丢。正确做法是X、Y方向都设unit cell边界Z方向设Floquet端口。正入射时Floquet端口只需要Mode 1斜入射时需要把端口模式数开到足够覆盖入射角否则S参数不完整。Floquet端口到样品表面的距离也要留够。端口面若紧贴基板近场高次模会影响S参数相位我一般让端口面距样品表面至少4到5 mm再在CST端口设置里用Deembed功能把参考面移到样品表面。这一步特别重要它直接决定S21相位的零点位置。反演时如果发现折射率曲线整体带有线性漂移先回来检查Deembed距离而不是怀疑算法有问题。端口极化方向也需要确认。Floquet端口默认有两个正交极化模式正入射计算时只激励一个极化另一个极化模式的结果如果不使用就不要在频带内开启它。开启多余的模式会引入虚假的耦合路径导致S11曲线出现不应有的毛刺。仿真完成后正确的结果是S11在谐振点附近深陷远离谐振时接近0 dBS21远离谐振时接近0 dB且相位随频率近似线性变化。看到这样的趋势再往下走才有意义。2.3 get_S_Parameter从后处理结果里把复数S11/S21导出来仿真结束后导航树里进入1D Results展开S-Parameters就能看到S1,1和S2,1。界面上默认显示幅度dB要看相位就切到Phase标签。想同时读两个频率点的值按住Shift键用鼠标在曲线上点两个Marker两个点的频率和幅值会同时显示在图的右上角这个操作在调谐振点位置时非常顺手。真正给反演用的是复数S参数不是图上那根曲线。导出方式有两种一种是在导航树里右键选中S参数Export Text导出带频率、实部、虚部的文件另一种更推荐在后处理模板Result Templates里创建Touchstone导出任务让CST把S参数写成标准的.s2p文件。.s2p文件第一行会注明格式是RI实部/虚部还是MA幅度/相位反演脚本读取时必须区分这是第3章的重点。用cst扫参做完参数扫描后每个样本的S参数会按样本名分好组通过Result Templates批量导出每个样本一个.s2p文件。文件名里带上扫描参数值比如P4.0_d0.5.s2p后续批处理脚本一眼就能看出几何参数。导出前确认端口阻抗是50欧姆频率单位是GHz这两个信息写在.s2p文件头里脚本读取时要处理单位换算否则算k0时频率差三个数量级反演结果完全不可用。3. S参数预处理从Touchstone文件到干净的反演输入3.1 Touchstone文件里有什么格式判断与频点对齐CST导出的.s2p文件本质是纯文本表格每行一个频点后面依次是S11实部、S11虚部、S21实部、S21虚部。文件头若干行以感叹号开头的注释紧接着是# GHz S RI R 50这样的格式行。这里的RI表示数据是实部虚部格式如果是# GHz S MA R 50则表示幅度相位格式。很多反演脚本读MA格式时当RI解析结果S参数变成一堆毫无意义的数反演出来全是噪声。还有一个容易忽略的点是频点对齐。扫参时不同几何参数的样本可能设置了不同的频点数量或者用了自适应扫频导致频点不均匀。反演前必须把所有样本的S参数插值到同一组频点上否则后面批处理时数组长度不一致代码直接报错。CST自适应扫频的频点在谐振区更密插值到均匀频点后谐振细节会丢失我一般保留原频点做反演只有需要对比多个样本时才统一插值。3.2 用Python解析S2P并按需去嵌给出一段能直接用的Python读取函数兼容RI和MA两种格式import numpy as np def load_touchstone(path): 读取CST导出的.s2p返回频率(Hz)、S11、S21复数数组 f, s11_re, s11_im, s21_re, s21_im [], [], [], [], [] is_mag_angle False with open(path, encodingutf-8, errorsignore) as fh: for line in fh: line line.strip() if not line or line.startswith(!): continue if line.startswith(#): if MA in line.upper(): is_mag_angle True continue parts line.split() if len(parts) 5: continue f.append(float(parts[0])) if is_mag_angle: # 幅度(dB或线性)/相位(度) mag 10 ** (float(parts[1]) / 20) ang np.deg2rad(float(parts[2])) s11_re.append(mag * np.cos(ang)) s11_im.append(mag * np.sin(ang)) mag 10 ** (float(parts[3]) / 20) ang np.deg2rad(float(parts[4])) s21_re.append(mag * np.cos(ang)) s21_im.append(mag * np.sin(ang)) else: # RI格式直接读实部虚部 s11_re.append(float(parts[1])) s11_im.append(float(parts[2])) s21_re.append(float(parts[3])) s21_im.append(float(parts[4])) f np.asarray(f) * 1e9 # GHz - Hz S11 np.asarray(s11_re) 1j * np.asarray(s11_im) S21 np.asarray(s21_re) 1j * np.asarray(s21_im) return f, S11, S21读文件时先把格式行的关键字MA或RI抓出来再决定每列怎么解析。频率单位统一换算成Hz这一步不及时做后面计算k0时全乱。CST如果导出的是线性幅度加角度注意幅度列若是dB要转回线性值很多CST版本的Touchstone幅度列默认写dB直接当线性值用会导致|S21|远超1反演出的折射率虚部完全失去物理意义。预处理还应该包含可选的端口去嵌。如果CST里没做Deembed端口参考面仍在端口物理位置S参数相位里含有一段空气传播路径。去嵌公式是S11乘上exp(2jk₀L)S21乘上exp(jk₀L)L是端口到样品表面的空气距离。这段运算放在预处理里既不用重新仿一遍又能多次微调L值看相位曲线的敏感性。3.3 预处理阶段最容易埋雷的三个细节第一个细节是相位解包裹的时机。反演需要连续的S21相位但atan函数返回的主值区间是[-π, π]相位跨过±180°时会跳变。Python的np.angle得到的就是跳变相位必须用np.unwrap解包裹才能用于反演。解包裹放在反演脚本里做预处理阶段只保证数据格式正确。第二个细节是频点密度。谐振区的S21相位变化非常快如果频点间隔过大相邻频点相位差超过π解包裹算法会判断错误反演出错误分支。做8到12 GHz仿真时我习惯至少取161个频点谐振附近相位变化再剧烈也能覆盖住。CST频域求解器的离散扫频里直接指定频点数量不要用默认的稀疏频点。第三个细节是保留原始文件。预处理脚本每次从.s2p原始文件读入不要覆盖式修改也不要只保存处理后的中间结果。分支选择不对时你大概率要回到原始S参数重新看相位原始文件一旦被改过前面几个小时的排查都白费。4. 参数反演核心算法NRW公式、分支处理与最小实现4.1 从S到ε/μ的四步换算与公式含义标题里的MáS泛指这类“CST出S参数、脚本做反演”的自动化方案它不是非用不可的黑匣子核心就是Nicolson-Ross-Weir方法也就是常说的NRW反演。NRW把超材料单元等效成一块均匀介质板用S11和S21反推这层“等效介质”的折射率n和波阻抗z再通过ε n/z、μ n·z得到等效介电常数和等效磁导率。换算过程分四步。第一步定义中间量V1 S21 S11、V2 S21 - S11进而算出X (1 - V1·V2)/(V1 - V2)。第二步用X求反射系数Γ X ± √(X² - 1)根号外符号的选择原则是保证|Γ| ≤ 1。第三步算传输系数T (V1 - Γ)/(1 - V1·Γ)然后由Γ求归一化阻抗z (1 Γ)/(1 - Γ)。第四步用传播相位求折射率n其中关键是要对S21的相位解包裹并处理对数函数的多值性问题。这套公式的物理前提是均匀化。单元周期远小于波长、样品厚度内场分布近似均匀时反演结果才有意义。基板厚度0.508 mm、工作频率10 GHz时等效介质厚度只有约λ₀/60均匀化误差可以接受如果单元电尺寸大到接近半个波长反演出的ε、μ会随厚度剧烈震荡那时要回头缩小周期而不是继续调算法。4.2 一份能直接跑的反演脚本下面是完整的NRW反演函数输入频率、S11、S21和样品厚度输出ε、μ、n、z四条曲线import numpy as np def nrw_retrieval(freq, S11, S21, d): NRW参数反演适用于CST导出的正入射S参数 freq: 频率, Hz S11, S21: 复数S参数数组 d: 样品厚度, m 返回: freq, eps_r, mu_r, n, z S11 np.asarray(S11, dtypecomplex) S21 np.asarray(S21, dtypecomplex) V1 S21 S11 V2 S21 - S11 denom V1 - V2 denom[abs(denom) 1e-12] 1e-12 # 防止除零 X (1 - V1 * V2) / denom # 两个根都算出来选|Gamma|1的作为物理根 G1 X np.sqrt(X**2 - 1) G2 X - np.sqrt(X**2 - 1) Gamma np.where(np.abs(G1) 1.0 1e-12, G1, G2) T (V1 - Gamma) / (1 - V1 * Gamma) # 阻抗和折射率 z (1 Gamma) / (1 - Gamma) phase np.unwrap(np.angle(T)) # 低频相位基线校准假定低频趋近普通介质而非负折射 phase phase - phase[0] t_abs np.abs(T) k0 2 * np.pi * freq / 3e8 n (-phase 1j * np.log(t_abs)) / (k0 * d) eps n / z mu n * z return freq, eps, mu, n, z这段代码的逻辑是先构造V1、V2让X变成一个中间量对Gamma的两个根用np.where按模长选择这一步直接处理了分支选择中最常见的问题T算出来后np.unwrap把S21相位展平再减去第一个频点的相位作为基线让折射率曲线在低频端回到正常介质的位置。模长部分取对数得到衰减项对应折射率虚部。参数上需要重点说明两点。第一厚度d是基板厚度而不是CST里unit cell的Z向总高度很多超材料单元在基板上下还有空气层如果直接把总高度拿来用n会系统性偏大。第二CST的时谐因子采用e^{jωt}无源材料的折射率虚部反演出来是负值对应的ε″、μ″也是负值这是约定问题不是错误和文献里用e^{-iωt}的符号习惯相反时肉眼看着像“损耗为负”实际上并没有违反无源判据。4.3 厚度、相位基线与分支选择反演可靠性的三个闸门反演结果可不可信第一个闸门是厚度。NRW公式里n与kd成反比厚度差10%n就差10%谐振频率位置也会跟着漂移。厚度取值我按基板厚度算铜层和印刷结构对厚度的贡献在10 GHz频段可以忽略。如果结构是金属-介质-金属三层厚度取两层介质加中间层的总厚度不能只取单层基板。第二个闸门是相位基线。phase - phase[0]这个操作把低频相位归零隐含假设是第一个频点远离谐振且相速接近自由空间。如果仿真频带从8 GHz起而谐振点在8.2 GHz基线校准就会失真。稳妥做法是让最低仿真频率比谐振频率低至少20%给基线校准留出“正常介质区”。第三个闸门是n的连续性。解包裹能解决2π跳变但NRW在某些频点仍会误入非物理分支表现是n在某一点突然跳变到另一个连续曲线。处理方式是先看n曲线是否平滑如果只有一个孤立的跳变点手动修正计算该点前后n的差值若差值超过k₀d分之π就整体减去或加上2π/(k₀d)。如果全频带都在无规律震荡多半是S参数本身的问题回第3章查相位解包裹和频点密度。5. S参数反演避坑五个高频问题的现象、原因与对策5.1 S21相位锯齿跳变解包裹后依然不平滑现象反演出的ε曲线呈高频锯齿状完全没有文献里那种平滑趋势n曲线上还有周期性尖峰。原因最常见的两个来源。一是频点太粗谐振区相邻频点相位差超过πnp.unwrap把真实变化误判成跳变二是.s2p文件幅度列是dB读取时没转线性|T|出现大于1的假象折射率虚部跟着震荡。解决把CST扫频点数加到160以上确认Touchstone格式行的MA/RI标记幅度用10**(dB/20)转线性值。改完重跑反演锯齿通常直接消失。5.2 反演出的损耗虚部符号和文献相反现象介质基板明明有损耗反演的ε″却是负值μ″在谐振区也是负值看着像材料的能量在增加。原因时谐因子约定。CST材料库和场求解器用e^{jωt}在这个约定下无源介质的介电常数虚部为负很多论文和教科书用e^{-iωt}虚部为正。两边符号习惯不同反演结果自然“相反”。这不是算法错误也不是材料异常。解决用能量判据代替符号判断。无源材料要求Re(n)0且Im(n)0CST约定下。如果反演结果的实部为正、谐振形态正常只是虚部符号和某篇论文相反按目标平台的约定统一取共轭即可不要为了“凑正号”乱改分支。5.3 厚度d取值不同结果差出一倍现象同一组S参数d取0.5 mm和1.0 mm低频折射率差一倍谐振频率位置也明显偏移。原因NRW的n与k₀d成反比d是全局比例因子。d取错n、n″同步缩放ε、μ也跟着偏。常见的取错方式是直接把CST里unit cell在Z方向的总高度当厚度这个值往往包含了上下空气层比实际样品厚度大好几倍。解决以基板物理厚度为准金属层在10 GHz以下厚度贡献小于1%可忽略。定义成脚本常量后不要每次手改多个样本对比时用同一个d值。5.4 曲线在某个频点“劈叉”或突变现象ε曲线整体平滑突然在某一个或几个频点跳到另一条曲线上回跳后趋势接不上S11模值在这些频点恰好接近0或接近1。原因NRW在|Γ|→0阻抗接近自由空间或|T|→0全反射时公式病态。S11接近0时X的计算变成0/0型的数值不稳定T接近0时log(T)的模长趋于负无穷虚部失去物理意义。另一个可能是该频点附近高阶Floquet模开始传播S参数不再只包含主模信息。解决先看S11模值确认异常频点是否处于匹配状态若是匹配点手动剔除该频点附近的几个数据再插值。若异常出现在高频端且单元周期接近半波长缩小周期尺寸后重新仿真。5.5 扫参到某些尺寸时仿真发散现象参数扫描到某个几何尺寸组合时频域求解器提示不收敛或发散S参数曲线在部分频点突然出现毛刺反演结果完全不可用。原因谐振结构出现极细的缝隙或尖锐拐角网格剖分不够高Q谐振在频域求解器里需要更多迭代次数薄金属片厚度小于网格步长时边缘场分布失真。解决频域求解器的自适应网格属性打开收敛目标设到-40 dB对薄金属结构改用Lossy Metal的thin sheet模型而不是按实体厚度剖分。毛刺若只出现在个别频点把这些频点从反演数据里剔除不要硬着头皮解包裹。6. 把反演流程自动化批处理、结果验证与表面参数进阶6.1 批量导出与自动反演参数扫描跑完CST的Result Templates已经为每个样本导出了.s2p文件剩下的就是让Python脚本批量处理。我一般用文件名解析几何参数再逐个调用前面三个函数读取、预处理、反演结果存成.npz文件import glob import numpy as np for path in sorted(glob.glob(./s2p/*.s2p)): freq, S11, S21 load_touchstone(path) freq, S11, S21 preprocess(freq, S11, S21, d_port0.0) freq, eps, mu, n, z nrw_retrieval(freq, S11, S21, d0.000508) tag path.split(/)[-1].replace(.s2p, ) np.savez(f./out/{tag}.npz, freqfreq, epseps, mumu, nn, zz) print(f{tag}: Re(eps) min{eps.real.min():.2f})这个脚本的输出是每个样本一个.npz文件后续只读.npz画图或做优化不再重复解析原始文件。打印语句输出每个样本的介电常数最小值扫参时一眼就能看出哪个尺寸组合出现了负折射行为。6.2 用等效参数回代验证反演结果不能只看曲线好不好看必须回代验证。做法是在CST里新建一个均匀介质板厚度与反演时用的d一致材料色散用反演得到的ε(f)、μ(f)表格导入边界和端口设置保持不变跑一次透射仿真。把回代的S21幅度和相位与原始超材料单元的S21放在同一张图里对比谐振点幅度偏差小于1 dB、相位趋势一致说明均匀化假设成立反演可信。对比项原始单元S21回代等效板S21谐振频率9.82 GHz9.78 GHz谐振深度-18.3 dB-17.9 dB带外相位斜率-2.1 rad/GHz-2.2 rad/GHz偏差大时先查厚度取值再查相位基线最后才怀疑公式本身。回代这一步也帮我判断过好几次分支选错曲线形态像谐振位置对不上基本都是n分支整体偏移了一个2π/(k₀d)的常数。6.3 薄超表面优先用表面参数如果单元厚度远小于工作波长比如d小于λ₀/20体积等效参数会变得对厚度极其敏感ε、μ的数值大得离谱物理意义弱化。这时更适合提取表面参数比如表面阻抗Z_s或表面极化率公式从传输线模型出发把超表面当作一个零厚度的阻抗边界S参数直接换算成Z_s。判断标准很简单厚度方向上的场分布如果不是近似均匀就别硬用NRW。把手头数据先做一遍回代验证验证不过就切到表面参数路线这是效率最高的工作顺序。我自己每次反演前固定检查三件事S21相位解包裹是否干净、S11模值有没有接近0或1的频点、厚度d是不是写成了常量而不是从文件名里再解析一遍。三次检查全过反演结果才敢往报告里放。希望帮到你。本文还有配套的精品资源点击获取