一种材料电磁参数测试的方法及用于材料电磁参数测试的多值性问题解决办法
【专利摘要】本发明提出了一种基于自由空间法进行材料电磁参数测试的方法和一种用于材料电磁参数测试的多值性问题解决办法。一种基于自由空间法进行材料电磁参数测试的方法，包括：步骤(一)、初始n值的计算，步骤(二)、后续所有n值的计算。本发明的对材料电磁参数测试中多值性问题的解决方法对待测样品的厚度没有特殊要求；对测试的起始频率无特殊要求；避免了群时延法对所有频点进行的繁琐的迭代求解；可适用于微波毫米波尤其适用于太赫兹波段自由空间法材料电磁参数测试。
【专利说明】
一种材料电磁参数测试的方法及用于材料电磁参数测试的多 值性问题解决办法
技术领域
[0001] 本发明涉及测试技术领域，特别涉及一种基于自由空间法进行材料电磁参数测试 的方法，还涉及一种用于材料电磁参数测试的多值性问题解决办法。
【背景技术】
[0002] 微波介质材料作为电磁波传输媒质已广泛地应用于微波的各个领域中，如微波通 信、卫星通信、导弹制导、电子对抗、雷达导航、遥感、遥测等系统已大量使用微波介质材料。 介质材料的电磁参数是指材料的复介电常数和复磁导率，它们表征了材料与电磁场的相互 作用。电磁参数测试随着介质材料的广泛应用具有越来越重要的意义。自由空间法、传输反 射法作为典型的材料介电常数测试方法，与其他方法相比具有测试方便、可进行宽带测量、 适用于电特和磁性符合材料无损测试等优点。
[0003] 以典型的NRW传输/反射法原理为例，其基本计算公式如下：
[0007] T = e-jYd (4)
[0008] 其中，Sn、S21均为散射参数，可由矢量网络分析仪获取，％为待测样品的相对空间 波阻抗，γ为样品的传播常数，d为待测样品厚度。联立以上四式，求解反射系数和透射系 数，得
[0012]其中，式(5)正负号选择| Γ |<1，根据空间波阻抗和传播常数与介电常数、磁导率
、丫=231八，式(3)、（4)可化为
[0015]联立以上两式，于是得到相对介电常数和磁导率计算公式：
[0018] 求解式（10)、（11)即可得样品的电磁参数。
[0019] 根据计算磁导率的公式（10)知，要计算样品的磁导率必须求解传输系数T的对数， 而T是复数，复数有以下性质，
[0022] n = int(dAg) (14)
[0023] 可见，虽然1/T的值是可以确定的，但是ln(l/T)却是周期性不确定的，虚部相差2n τ，η的取值由波导波长和样品的厚度d决定，如何确定η的取值既是材料电磁参数测试中 的多值性问题。
[0024] 常用的虚部补偿法是一种比较简单的求解多值性问题的解决办法，空气的传播 常数为：
[0026]样品中的传播常数为
[0028] 由式（15)、（16)可知，传播常数γ的虚数部分与测试频率f之间应该是线性相关且 是递增的函数，图2虚线所示传播常数γ的虚数部分随着测试频率f的增大而呈周期性变 化，因此可以通过线性递增函数原理进行虚部补偿，补偿之后传播常数γ的虚数部分与测 试频率f的关系图如图2中实线所示。
[0029] 虚部补偿法确定η的方法，当测量频率升高时（fm<fm+1)，若有imag( γ (fm+1))< imag( γ (fm))，则知传播常数γ的虚数部分发生了周期性的变化，因此，频率fm+1所对应的η 值应该相应的加1(即η = η+1)。测量频率尽量采用扫频的方式进行测量，使得传播常数的虚 数部分在η的范围内变化。该方法相对于群时延法比较简单，但是该方法的缺点在于要对测 量的起始频率进行估算，而初始频率的确定式为
，在不知道介质材料的电 磁参数的情况下很难估算起始频率。
[0030] 自由空间法进行材料电磁参数测试的多值性问题解决方案，目前普遍采用的方法 主要有以下不足：
[0031] (1)限定待测样品厚度，一般使待测样品厚度低于波导波长的一半，由式（14)可 见，满足此条件后，η的取值只有0,因此不存在多值性问题了。该方法的缺点一是需要估计 样品的波导波长，在样品电磁特性未知的情况下很难做到，二是不适合高频测量，高频条件 下，波长过小，样品厚度难以加工到半波长。
[0032] (2)相位群时延法的原理是η从0依次递增计算样品介电常数，当电磁波通过待测 样品时，其群时延有两种计算方法，一是由当前计算的介电常数和磁导率，根据群时延与 介电常数、磁导率、厚度以及工作频率的关系计算群时延;二是由矢量网络分析仪直接读取 电磁波通过样品的群时延。当两种方法得到的群时延一致时，则表示当前的η为正确值。该 方法的缺点是每个频点都要进行多次迭代，比较繁琐，而且判断两种方法计算的群时延是 否一致需要设定误差范围，不同频点的误差差距过大会严重影响判断结果，从而影响计算 精度。
[0033] (3)虚部补偿法，实际测试的传播常数虚部随着测试频率的增加而呈现周期性变 化，根据传播常数虚部随着频率应该是递增的规律，每当传播常数虚部突然降低时，则对η 进行加1。该方法相对于群时延法比较简单，但是该方法的缺点在于要对测量的起始频率进 行估算，而初始频率的确定式为
，在不知道介质材料的电磁参数的情况下 很难估算起始频率。
【发明内容】
[0034] 为解决上述现有技术中的不足，本发明的目的在于提出一种基于自由空间法进行 材料电磁参数测试的方法和一种用于材料电磁参数测试的多值性问题解决办法，先采用相 位群时延法求解η的初始值，再采用"虚部补偿"法求解后续所有η值。
[0035] 本发明的技术方案是这样实现的：
[0036] -种基于自由空间法进行材料电磁参数测试的方法，包括以下步骤：
[0037] 步骤(一）、初始η值的计算：
[0038] 首先采用相位-群时延法计算测试带宽的第一个频点fQ对应的η，在工作频点fo时， η从零开始依次循环，加1递增，在当前循环下，η是已知的，根据式（10)~（13)计算当前介电 常数ε#Ρ磁导率μΓ，
[0043]并按照式（17)计算电磁波通过样品的理论群时延数值，
[0045]或者，群时延通过散射参数校准之后，待测样品两端的S21的相位Φ计算得到，计 算公式如式(18)，
[0047]步骤(二）、后续所有η值的计算：
[0048] 初始η值计算后，按照虚部补偿法进行计算后续所有的η值，虚部补偿以周期性变 化的S21相位为η取值的判断依据，具体为：
[0049] 从频率第二个点m彡2开始，当测量频率升高时，即fm<fm+1，若有phase(S21(f m))> phase (S21 (fm+1))，则知S21的相位部分发生了周期性的变化，此时，频率匕+1所对应的η值相 应的加1，即n = n+l。
[0050]可选地，根据式（17)、（18)计算群时延，当二者结果有较大偏差时，增加 η的数值， 继续计算样品电磁参数并按照式（17)计算群时延，直到与式（18)的计算结果满足误差精度 时，则表示当前循环下η值正确。
[0051 ]可选地，测量频率采用扫频的方式进行测量，使得S21的相位在± 2π的范围内变 化。
[0052]本发明还提出了一种用于材料电磁参数测试的多值性问题解决办法，基于上述 的材料电磁参数测试的方法，包括以下步骤：
[0053] (1)首先，对测试夹具仅进行二端口散射参数校准，经过校准后，使得矢量网络分 析仪读取的S21参数即为电磁波通过待测样品的S21;
[0054] (2)设定周期性起始值η = 0，频率fm的第一个频点数记为m=0;
[0055] (3)根据式（17)、（18)分别计算理论群时延1。31和测试群时延1?^
[0056] (4)设置误差项error，判断是否满足条件卜c；ai-T：me3a | /卜。ai | terror，若是则表示 当前η值即为初始η值，否则进行n = n+l，并返回步骤(3)重新计算理论群时延测试群 时延Tmea ;
[0057] (5)令 m=m+l;
[0058] (6)判断m是否是最后的频点，若是则执行步骤(8)跳出循环，否则执行步骤(7);
[0059] (7)判断当前S21的相位phase(S21(fm))是否大于前一个频点对应的S21的相位 phase(S21(fm-1))，若是则取n(m) =n(m-l)，并返回步骤(5)，若否则取n(m) =n(m-l)+l并返 回步骤(5);
[0060] (8)完成多值性问题的求解，输出所有η值。
[0061 ] 可选地，设置误差项error取1 %。
[0062]本发明的有益效果是：
[0063] (1)对待测样品的厚度没有特殊要求；
[0064] (2)对测试的起始频率无特殊要求；
[0065] (3)避免了群时延法对所有频点进行的繁琐的迭代求解；
[0066] (4)可适用于微波毫米波尤其适用于太赫兹波段自由空间法材料电磁参数测试。
【附图说明】
[0067] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以 根据这些附图获得其他的附图。
[0068]图1为测试场景示意图；
[0069]图2为现有的虚部补偿法原理示意图；
[0070]图3为本发明的基于自由空间法进行材料电磁参数测试的方法电磁波穿过待测样 品原理不意图；
[0071 ]图4为本发明相位补偿前后变换情况示意图；
[0072]图5为本发明对材料电磁参数测试中多值性问题的解决方法的流程图。
【具体实施方式】
[0073] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0074] 多值性问题是自由空间法、传输反射法求解电磁参数必须要解决的首要问题之 一，电磁参数求解过程中，需要求解相位信息，相位具有2ηπ的周期性，因此存在多值性问 题。关于多值性问题目前主要有三种解决方案：
[0075] 一是限制待测样品的厚度，一般使样品厚度低于波导波长的一半，满足此条件后， η的值只有0,因此就不存在多值性问题了。该方法常用于微波波段传输反射法或自由空间 法材料电磁参数测试，但不再适用于高频测量，以太赫兹频段为例，波导波长最高也在毫米 量级，样品太薄难以加工且大大增加了厚度不确定度，严重影响测试精度；
[0076] 二是通过预估待测样品的介电常数和磁导率，从而计算出待测样品的波导波长， 再根据相关公式计算η的取值。显然，当待测材料的电磁特性是未知的，该方法是不适用的；
[0077] 三是通过调整传播常数虚部的方法来解决多值性问题，该方法最大的难题是需要 限制样品厚度或限定测试起始频率以保证η的初值为零。
[0078] 本发明的目的在于提出一种基于自由空间法进行材料电磁参数测试的方法和一 种对材料电磁参数测试中多值性问题的解决方法，不增加任何硬件成本、对材料厚度无特 殊要求、对工作起始频率无特殊要求、尤其适用于超宽带ΤΗζ材料电磁参数测试。
[0079]电磁通过待测样品情景如图3所示，在结合虚部补偿法和群时延法优点的基础上， 本发明材料电磁参数测试的方法包括以下两步：
[0080] -、初始η值的计算。
[0081]本发明首先采用相位-群时延法计算测试带宽的第一个频点fo对应的η，即η的初 始值。在工作频点f〇时，η从零开始依次循环，加1递增，在当前循环下，η是已知的，因此可以 根据式（10)~（13)计算当前介电常数^和磁导率y r，并按照下式（17)计算电磁波通过样品 的理论群时延数值
[0083]此外，群时延还可以通过散射参数校准之后，待测样品两端的S21的相位Φ计算 得到，具体计算公式如式(18)，
[0085]根据式（17)、（18)可见群时延有两种计算方法，实际电磁波通过固定厚度的样品 的群时延应该是唯一的，因此两者的计算结果应该是一致的，当二者结果有较大偏差时，增 加 η的数值，继续计算样品电磁参数并按照式（17)计算群时延，直到与式（18)的计算结果满 足一定的误差精度时，则表示当前循环下η值正确，从而解决了初始频点fo的多值性问题 了。因为只有一个频点需要迭代比较求解，因此这一步骤的计算速度很快而且精确。
[0086] 二、后续所有η值的计算。
[0087] 初始η值计算后，可以按照虚部补偿法进行计算后续所有的η值了，由于相位群时 延法迭代比较繁琐，宽带频率范围内误差判断难以一致，因此后续η值的计算不再采用此方 法了，而是采用改进的虚部补偿法，虚部补偿以周期性变化的S21相位为η取值的判断依据。 散射参数经过全二端口网络校准后，S21的相位随着工作频率的升高呈现处周期性的变化 趋势，如图4虚线所示，经过类似于虚部补偿的相位补偿后，S21的相位将变化成为图4所示 的实线部分。
[0088] 具体相位补偿法确定η的方法是，从频率第二个点开始(m多2)，当测量频率升高时 (fm<fm+l)，若有口1^86(321(匕)）>?11&86(321(匕 +1))，则知321的相位部分发生了周期性的 变化，此时，频率^+1所对应的η值应该相应的加1 (即n = n+l)。测量频率尽量采用扫频的方 式进行测量，使得S21的相位在± 2π的范围内变化。该方法相对于群时延法比较简单，由于 起始点对应的η值已经由上一步骤计算得出，因此不必再按式
限定初始频 率了，避免了在不知道介质材料的电磁参数情况下估算起始频率。
[0089] 利用本发明的上述方法，本发明还提出了一种用于材料电磁参数测试中多值性问 题的解决方法，如图5所示，包括以下步骤：
[0090] (1)首先，对测试夹具仅进行二端口散射参数校准，经过校准后，使得矢量网络分 析仪读取的S21参数即为电磁波通过待测样品的S21;
[0091 ] (2)设定周期性起始值η = 0，频率fm的第一个频点数记为m=0;
[0092] (3)根据式（17)、（18)分别计算理论群时延Tcal和测试群时延i mea;
[0093] (4)设置误差项error，一般取1%，判断是否满足条件卜cai-T mea|/卜cal|彡error， 若是则表示当前η值即为初始η值，否则进行n = n+l，并返回步骤三重新计算理论群时延Tcal 和测试群时延
[0094] (5)令 m=m+l;
[0095] (6)判断m是否是最后的频点，若是则执行步骤8跳出循环，否则执行步骤7;
[0096] (7)判断当前S21的相位phase(S21(fm))是否大于前一个频点对应的S21的相位 phase(S21(fm-1))，若是则取n(m) =n(m-l)，并返回步骤5,若否则取n(m)=n(m-l)+l并返回 步骤5
[0097] (8)完成多值性问题的求解，输出所有η值。
[0098] 本发明的对材料电磁参数测试中多值性问题的解决方法对待测样品的厚度没有 特殊要求;对测试的起始频率无特殊要求;避免了群时延法对所有频点进行的繁琐的迭代 求解;可适用于微波毫米波尤其适用于太赫兹波段自由空间法材料电磁参数测试。
[0099] 以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精 神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种基于自由空间法进行材料电磁参数测试的方法，其特征在于，包括W下步骤： 步骤(一）、初始η值的计算： 首先采用相位-群时延法计算测试带宽的第一个频点f日对应的η，在工作频点f日时，η从 零开始依次循环，加1递增，在当前循环下，η是已知的，根据式（10)~（13)计算当前介电常 数Er和磁导率化，或者，群时延通过散射参数校准之后，待测样品两端的S21的相位Φ计算得到，计算公 式如式(18)，(18) 步骤(二）、后续所有η值的计算： 初始η值计算后，按照虚部补偿法进行计算后续所有的η值，虚部补偿W周期性变化的 S21相位为η取值的判断依据，具体为： 从频率第二个点m>2开始，当测量频率升高时，即若有phase(S21(fm))> 地ase (S21 (f))，则知S21的相位部分发生了周期性的变化，此时，频率f所对应的η值相 应的加1，即η = η+1。2. 如权利要求1所述的一种基于自由空间法进行材料电磁参数测试的方法，其特征在 于，根据式（17)、（18)计算群时延，当二者结果有较大偏差时，增加 η的数值，继续计算样品 电磁参数并按照式（17)计算群时延，直到与式（18)的计算结果满足误差精度时，则表示当 前循环下η值正确。3. 如权利要求1所述的一种基于自由空间法进行材料电磁参数测试的方法，其特征在 于，测量频率采用扫频的方式进行测量，使得S21的相位在±2π的范围内变化。4. 一种用于材料电磁参数测试的多值性问题解决办法，其特征在于，基于权利要求1至 3任一项所述的材料电磁参数测试的方法，包括W下步骤： (1) 首先，对测试夹具仅进行二端口散射参数校准，经过校准后，使得矢量网络分析仪 读取的S21参数即为电磁波通过待测样品的S21; (2) 设定周期性起始值η = 0,频率fm的第一个频点数记为m=0; (3) 根据式（17)、（18)分别计算理论群时延τ。。讯测试群时延Tmea; (4) 设置误差项error,判断是否满足条件I Teai-Tmeal/I Teal I《error,若是则表示当前η 值即为初始η值，否则进行η = η+1，并返回步骤（3)重新计算理论群时延Teal和测试群时延 Tmea; (5) 令m=m+l; (6) 判断m是否是最后的频点，若是则执行步骤(8)跳出循环，否则执行步骤(7); (7) 判断当前S21的相位phase(S21(fm))是否大于前一个频点对应的S21的相位phase (S21(fm-i))，若是则取n(m)=n(m-l)，并返回步骤(5)，若否则取n(m)=n(m-l) + l并返回步 骤(5); (8) 完成多值性问题的求解，输出所有η值。5.如权利要求4所述的用于材料电磁参数测试的多值性问题解决办法，其特征在于，设 置误差项error取1 %。
【文档编号】G01R33/12GK106093810SQ201610353303
【公开日】2016年11月9日
【申请日】2016年5月19日
【发明人】胡大海, 赵锐, 杜刘革, 王亚海
【申请人】中电科仪器仪表有限公司