CN115060978A - 一种基于时域分析法的介电常数估计方法 - Google Patents
一种基于时域分析法的介电常数估计方法 Download PDFInfo
- Publication number
- CN115060978A CN115060978A CN202210744803.2A CN202210744803A CN115060978A CN 115060978 A CN115060978 A CN 115060978A CN 202210744803 A CN202210744803 A CN 202210744803A CN 115060978 A CN115060978 A CN 115060978A
- Authority
- CN
- China
- Prior art keywords
- medium
- dielectric constant
- time
- time domain
- thickness
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
- G01R27/26—Measuring inductance or capacitance; Measuring quality factor, e.g. by using the resonance method; Measuring loss factor; Measuring dielectric constants ; Measuring impedance or related variables
- G01R27/2617—Measuring dielectric properties, e.g. constants
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/02—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring length, width or thickness
- G01B7/06—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring length, width or thickness for measuring thickness
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
- G01R27/04—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant in circuits having distributed constants, e.g. having very long conductors or involving high frequencies
- G01R27/06—Measuring reflection coefficients; Measuring standing-wave ratio
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
- Length-Measuring Devices Using Wave Or Particle Radiation (AREA)
Abstract
本发明提供一种基于时域分析法的介电常数估计方法,包括设置介质厚度测试系统、测量天线驻波、反射系数S11频域数据处理、S11时频变换计算介质厚度、S11时频变换计算介电常数、提高介电常数测量精度,避免了传统时域分析法时域选通带来的误差,误差来源仅为时域分辨精度,该精度可由chirp‑z逆变换或者频域补零进行优化。并且抓住了被测介质的不连续性在时域的物理意义,根据已经得到的厚度和电磁场相关理论计算介电常数的大小。
Description
技术领域
本发明属于微波测量技术领域,是一种基于时域分析法的介电常数估计方法。
背景技术
矢量网络分析仪(VNA)是一种精密的综合性微波测量仪器,它能够测量单端口网络或多端口网络的散射参数(S参数),矢量网络分析仪往往采用扫频的方式测量被测器件(DUT)的S参数(包括幅度特性和相位特性),然后通过内部的运算模块计算它的阻抗参数,导纳参数,传输参数等数据,但是这些数据均为频域数据。很多场景下,例如对器件特性的不连续性进行检测时,往往需要判断不连续点的位置或多径传输的问题,观测其S参数的时域数据可以非常直观的看到DUT的特性。矢量网络分析仪的时域测量技术正是该方法的应用,该技术是基于傅里叶变换理论,采用离散傅里叶变换实现频域数据和时域数据的转换。
常规时域分析法采用时域选通法对不连续点的时域波形进行开窗处理,再转换至频域计算反射系数的模值,再由平面波对理想介质分界平面的垂直入射场景下反射系数与波阻抗的关系进而推算介质的介电常数。再由介电常数结合时域不连续点的时间差计算得到介质的厚度,但是此方法的误差因素较多,首先该方法受测试带宽影响较大,带宽不足会导致时域不连续点波形发生交叠,导致无法进行时域选通。并且在时频变换过程中为了减小旁瓣,振铃效应等影响,会对信号加不同衰减系数的窗函数,关于窗的选择也是误差的来源之一。计算时域不连续点的时间差也会存在时间读取的误差。如果介质是多层介质或者有损介质,则计算误差和计算难度会大大提升,计算精度也会降低。
发明内容
鉴于以上所述现有技术的缺点,本发明的目的在于提供一种基于时域分析法的介质介电常数估计方法,采用矢量网络分析仪测量两个测试系统DUT的S参数作为参照,将两组数据分别扣除天线驻波后做除法运算,然后将数据转换到时域,通过对比分析计算介质的厚度,再由被测器件的不连续性处发生的时间相对位置差值读取电磁波在介质板中传输所耗的时间,根据此时间和已得到的介质厚度计算介电常数。本发明旨在利用新增加的对照组测试系统,结合数字信号处理相关理论,估算介质板的介电常数。
为实现上述发明目的,本发明技术方案如下:
一种基于时域分析法的介电常数估计方法,包括如下步骤:
(3)测量天线驻波;
(4)反射系数S11频域数据处理
将介质厚度测试系统的S11频域数据扣除天线驻波后的数据,除以对照组测试系统的S11频域数据扣除天线驻波的数据,得到S11div;
(5)S11时频变换计算介质厚度
根据步骤(4)得到的S11div,转换到时域,测量时间观测范围终点的时移T0;根据时移大小即可计算出介质厚度;
(6)S11时频变换计算介电常数
(7)提高介电常数测量精度
采用线性调频逆Z变换,提高时间分辨精度,得到更加精确的时间差Δt,从而改善介电常数的精度。
作为优选方式,步骤(1)进一步为:
(1)设置对照组测试系统
作为优选方式,步骤(2)进一步为:
(2)设置介质厚度测试系统
介质厚度测试系统包括矢量网络分析仪、聚束天线、待测介质板、金属板;
聚束天线垂直照射待测介质板,矢量网络分析仪连接聚束天线的校准处,即馈源接收端口,待测介质板背后紧贴金属板,聚束天线垂直于介质板摆放;聚束天线到金属板的距离和对照组测试系统保持一致;通过矢量网络分析仪测试聚束天线校准处的反射系数S11,测试数据记为
作为优选方式,步骤(3)进一步为:
(3)测量天线驻波
通过矢量网络分析仪测试聚束天线在无待测介质板和金属板时的反射系数,即天线驻波记为S11sw。
作为优选方式,步骤(4)进一步为:
(4)反射系数S11频域数据处理
经严格的电磁场理论公式推导得:
根据以上公式将介质厚度测试系统的S11频域数据扣除天线驻波后的数据,除以对照组测试系统的S11频域数据扣除天线驻波的数据得到S11div。
作为优选方式,步骤(5)进一步为:
(5)S11时频变换计算介质厚度
通过matlab将待分析的频域数据S11div转换到时域进行分析:测量的第一个不连续点的时移T0;这个时移的物理意义即为电磁波在真空中传播两倍介质厚度所耗时间;根据时移大小即可计算出介质厚度;即介质厚度d=c*T0/2,c为真空中电磁波的传播速度。
作为优选方式,步骤(6)进一步为:
(6)S11时频变换计算介电常数
用matlab分析其时域波形,波形中第一个反射峰则代表被测材料的第一个不连续处,即空气与介质板的交界面的第一次反射;第二个反射峰代表被测材料的第二个不连续处,即介质板和金属板的介质面第一次反射;前两次反射在时域上存在一定的时间差,该时间差即为电磁波在介质内传播两倍介质厚度所耗时间的两倍;如下公式所示:
d为介质厚度,c为真空中电磁波传播的速度,ε代表待测介质板的介电常数;
根据公式(2)得到待测介质板的介电常数。
作为优选方式,步骤(7)进一步为:
(7)提高介电常数测量精度
本发明的有益效果在于:本发明的方法避免了传统时域分析法时域选通带来的误差,误差来源仅为时间差读取的精度,该精度可由chirp-z逆变换或者频域补零进行优化。并且抓住了被测介质的不连续性在时域的物理意义,根据已经得到的厚度和电磁场相关的理论计算介电常数的大小。
附图说明
图1是本发明的流程图;
图2是本发明对照组测试系统图;
图3是本发明介质厚度测试系统图,图3中d1为聚束天线到待测介质板表面的距离,d2为介质的厚度,图2中d1与图3保持一致,d1+d2表示聚束天线到金属板表面的距离;
1为矢量网络分析仪,2为聚束天线,3为待测介质板,4为金属板。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。
实施例1
本实施例提供一种基于时域分析法的介电常数估计方法,包括如下步骤:
(3)测量天线驻波;
(4)反射系数S11频域数据处理
将介质厚度测试系统的S11频域数据扣除天线驻波后的数据,除以对照组测试系统的S11频域数据扣除天线驻波的数据,得到S11div;
(5)S11时频变换计算介质厚度
根据步骤(4)得到的S11div,转换到时域,测量时间观测范围终点的时移T0;根据时移大小即可计算出介质厚度;
(6)S11时频变换计算介电常数
(7)提高介电常数测量精度
采用线性调频逆Z变换,提高时间分辨精度,得到更加精确的时间差Δt,从而改善介电常数的精度。
实施例2
本实施例提供一种基于时域分析法的介电常数估计方法,包括如下步骤:
(1)设置对照组测试系统
(2)设置介质厚度测试系统
介质厚度测试系统包括矢量网络分析仪、聚束天线、待测介质板、金属板;
矢量网络分析仪连接聚束天线的校准处,即馈源接收端口,待测介质板背后紧贴金属板,聚束天线垂直于介质板摆放;聚束天线到金属板的距离和对照组测试系统保持一致;通过矢量网络分析仪测试聚束天线校准处的反射系数S11,测试数据记为
(3)测量天线驻波
通过矢量网络分析仪测试聚束天线在无待测介质板和金属板时的反射系数,即天线驻波记为S11sw。
(4)反射系数S11频域数据处理
经严格的电磁场理论公式推导得:
根据以上公式将介质厚度测试系统的S11频域数据扣除天线驻波后的数据,除以对照组测试系统的S11频域数据扣除天线驻波的数据得到S11div。
(5)S11时频变换计算介质厚度
通过matlab将待分析的频域数据S11div转换到时域进行分析:测量的第一个不连续点的时移T0;这个时移的物理意义即为电磁波在真空中传播两倍介质厚度所耗时间;根据时移大小即可计算出介质厚度;即介质厚度d=c*T0/2,c为真空中电磁波的传播速度。
(6)S11时频变换计算介电常数
用matlab分析其时域波形,波形中第一个反射峰则代表被测材料的第一个不连续处,即空气与介质板的交界面的第一次反射;第二个反射峰代表被测材料的第二个不连续处,即介质板和金属板的介质面第一次反射;前两次反射在时域上存在一定的时间差,该时间差即为电磁波在介质内传播两倍介质厚度所耗时间的两倍;如下公式所示:
d为介质厚度,c为真空中电磁波传播的速度,ε代表待测介质板的介电常数;
根据公式(2)得到待测介质板的介电常数。
(7)提高介电常数测量精度
很明显在波形起始处存在天线驻波,随后可以看到明显的三个反射峰,第一个反射峰则代表被测材料的第一个不连续处,体现为类似冲激函数的反射峰,即空气与介质板的交界面的第一次反射。第二个反射峰代表被测材料的第二个不连续处,即介质板和金属板的介质面第一次反射,由于金属是全反射,所以第二层反射明显增大。第三个或者后面的反射峰则是介质内多层反射的结果,考虑到传播的损耗,后面层次的反射会越来越小。前两次反射在时域上存在一定的时间差,该时间差即为电磁波在介质内传播两倍介质厚度所耗时间的两倍。
因厚度已经在前面的分析中计算得到,在时域读取该二者不连续处的冲激响应的时间差则可计算得到介质的介电常数。进一步需要解决的问题同样也是时间差读取的精度问题。本发明采用chirp-z逆变换对第一层反射和第二层反射的附近进行局部放大,得到的数据如图5所示:
很明显时域波形图更细致平滑,时间分辨精度更高,读取峰值对应的时间则会更加精确。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
Claims (8)
1.一种基于时域分析法的介电常数估计方法,其特征在于,包括如下步骤:
(3)测量天线驻波;
(4)反射系数S11频域数据处理
将介质厚度测试系统的S11频域数据扣除天线驻波后的数据,除以对照组测试系统的S11频域数据扣除天线驻波的数据,得到S11dic;
(5)S11时频变换计算介质厚度
根据步骤(4)得到的S11div,转换到时域,测量时间观测范围终点的时移T0;根据时移大小即可计算出介质厚度;
(6)S11时频变换计算介电常数
(7)提高介电常数测量精度
采用线性调频逆Z变换,提高时间分辨精度,得到更加精确的时间差Δt,从而改善介电常数的精度。
4.根据权利要求1所述的一种基于时域分析法的介电常数估计方法,其特征在于:步骤(3)进一步为:
(3)测量天线驻波
通过矢量网络分析仪测试聚束天线在无待测介质板和金属板时的反射系数,即天线驻波记为S11sw。
7.根据权利要求1所述的一种基于时域分析法的介电常数估计方法,其特征在于:步骤(6)进一步为:
(6)S11时频变换计算介电常数
用matlab分析其时域波形,波形中第一个反射峰则代表被测材料的第一个不连续处,即空气与介质板的交界面的第一次反射;第二个反射峰代表被测材料的第二个不连续处,即介质板和金属板的介质面第一次反射;前两次反射在时域上存在一定的时间差,该时间差即为电磁波在介质内传播两倍介质厚度所耗时间的两倍;如下公式所示:
d为介质厚度,c为真空中电磁波传播的速度,ε代表待测介质板的介电常数;
根据公式(2)得到待测介质板的介电常数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210744803.2A CN115060978B (zh) | 2022-06-28 | 2022-06-28 | 一种基于时域分析法的介电常数估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210744803.2A CN115060978B (zh) | 2022-06-28 | 2022-06-28 | 一种基于时域分析法的介电常数估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115060978A true CN115060978A (zh) | 2022-09-16 |
CN115060978B CN115060978B (zh) | 2023-05-23 |
Family
ID=83204830
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210744803.2A Active CN115060978B (zh) | 2022-06-28 | 2022-06-28 | 一种基于时域分析法的介电常数估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115060978B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1143858A (zh) * | 1995-03-22 | 1997-02-26 | 株式会社村田制作所 | 介电谐振器以及使用这种介电谐振器的介电谐振器装置 |
US20030102907A1 (en) * | 2001-11-29 | 2003-06-05 | Tasker Paul Juan | Methods and apparatus for time-domain measurement with a high frequency circuit analyzer |
US20100001747A1 (en) * | 2007-02-14 | 2010-01-07 | Joint-Stock Company 'high Tech' | Method and a device for measuring dielectric characteristics of material bodies |
EP3156784A1 (en) * | 2015-10-12 | 2017-04-19 | Vrije Universiteit Brussel | Enhanced characterization of dielectric properties |
CN109669075A (zh) * | 2019-01-25 | 2019-04-23 | 电子科技大学 | 基于开口矩形波导的介质复介电常数无损反射测量方法 |
CN110308334A (zh) * | 2019-05-31 | 2019-10-08 | 西安空间无线电技术研究所 | 一种旋转关节驻波比及插损测试方法 |
RU2713162C1 (ru) * | 2019-07-29 | 2020-02-04 | Акционерное общество «Обнинское научно-производственное предприятие «Технология» им. А.Г.Ромашина» | Способ определения диэлектрической проницаемости материала |
CN111474411A (zh) * | 2020-04-24 | 2020-07-31 | 京信通信技术(广州)有限公司 | 相对介电常数测试系统、方法、装置和存储介质 |
CN112505429A (zh) * | 2020-12-08 | 2021-03-16 | 电子科技大学 | 基于同轴带状线谐振器的复介电常数测试系统及测试方法 |
CN214953147U (zh) * | 2021-05-10 | 2021-11-30 | 中电科西北集团有限公司 | 一种电磁参数测试仪 |
-
2022
- 2022-06-28 CN CN202210744803.2A patent/CN115060978B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1143858A (zh) * | 1995-03-22 | 1997-02-26 | 株式会社村田制作所 | 介电谐振器以及使用这种介电谐振器的介电谐振器装置 |
US20030102907A1 (en) * | 2001-11-29 | 2003-06-05 | Tasker Paul Juan | Methods and apparatus for time-domain measurement with a high frequency circuit analyzer |
US20100001747A1 (en) * | 2007-02-14 | 2010-01-07 | Joint-Stock Company 'high Tech' | Method and a device for measuring dielectric characteristics of material bodies |
EP3156784A1 (en) * | 2015-10-12 | 2017-04-19 | Vrije Universiteit Brussel | Enhanced characterization of dielectric properties |
CN109669075A (zh) * | 2019-01-25 | 2019-04-23 | 电子科技大学 | 基于开口矩形波导的介质复介电常数无损反射测量方法 |
CN110308334A (zh) * | 2019-05-31 | 2019-10-08 | 西安空间无线电技术研究所 | 一种旋转关节驻波比及插损测试方法 |
RU2713162C1 (ru) * | 2019-07-29 | 2020-02-04 | Акционерное общество «Обнинское научно-производственное предприятие «Технология» им. А.Г.Ромашина» | Способ определения диэлектрической проницаемости материала |
CN111474411A (zh) * | 2020-04-24 | 2020-07-31 | 京信通信技术(广州)有限公司 | 相对介电常数测试系统、方法、装置和存储介质 |
CN112505429A (zh) * | 2020-12-08 | 2021-03-16 | 电子科技大学 | 基于同轴带状线谐振器的复介电常数测试系统及测试方法 |
CN214953147U (zh) * | 2021-05-10 | 2021-11-30 | 中电科西北集团有限公司 | 一种电磁参数测试仪 |
Non-Patent Citations (3)
Title |
---|
AMAN 等: "Retrieval of electrical and physical properties of dielectric samples using time-domain multiple reflection method" * |
AMAN 等: "Retrieval of electrical and physical properties of dielectric samples using time-domain multiple reflection method", 《IET MICROWAVES, ANTENNAS & PROPAGATION》 * |
李天凡: "双负介质电磁参数测量方法研究" * |
Also Published As
Publication number | Publication date |
---|---|
CN115060978B (zh) | 2023-05-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109669075B (zh) | 基于开口矩形波导的介质复介电常数无损反射测量方法 | |
Katz | Method to resolve microphone and sample location errors in the two-microphone duct measurement method | |
JP2006208060A (ja) | 伝送遅延評価システムおよび伝送遅延評価方法 | |
US7385932B2 (en) | Wideband frequency domain reflectometry to determine the nature and location of subscriber line faults | |
CN110031838B (zh) | 一种不需要参考发射波形的穿墙雷达墙体参数估计方法 | |
CN113504414A (zh) | 一种等离子体复介电常数瞬态微波透射检测方法及装置 | |
Dhaene et al. | Calibration and normalization of time domain network analyzer measurements | |
CN114265050A (zh) | 一种脉冲压缩和差测距方法 | |
CN113758617B (zh) | 一种基于宽带扫频信号频域计算的应力梯度高效无损检测系统及其检测方法 | |
CN110425972A (zh) | 一种宽工作频段反射式电厚度测试方法 | |
CN115060978B (zh) | 一种基于时域分析法的介电常数估计方法 | |
CN115143875B (zh) | 一种基于时域分析法的介质厚度估计方法 | |
CN109752705A (zh) | 高频水声阵列性能参数测量方法及系统、设备及存储介质 | |
CN109752633B (zh) | 一种对变电站局部放电信号进行定位的方法及系统 | |
US10345339B2 (en) | Group delay based averaging | |
CN112858393A (zh) | 一种基于频率步进原理的tdr测试电导率方法 | |
CN112034406B (zh) | 一种便携式电厚度测试仪相位校准方法 | |
CN111913161B (zh) | 一种提升nlfm波形雷达目标测角精度的方法 | |
CN110361685B (zh) | 一种宽带示波器探头传输特性校准方法和系统 | |
CN104297733B (zh) | 复合背景下时域信号的处理方法及复合信号处理系统 | |
CN109696660B (zh) | 一种检测自由声场传声筒幅值灵敏度与相位精确测量方法 | |
CN114265017B (zh) | 一种基于数字信号处理的相噪测量方法 | |
CN112033518A (zh) | 一种基于正弦函数拟合的水听器灵敏度标定方法 | |
CN116840807B (zh) | 一种基于探地雷达系统的全波反演介电常数估计方法 | |
Scott et al. | Acoustic vector-corrected impedance meter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |