CN104931912B - 一种矢量网络分析仪的时域补偿方法 - Google Patents
一种矢量网络分析仪的时域补偿方法 Download PDFInfo
- Publication number
- CN104931912B CN104931912B CN201510323996.4A CN201510323996A CN104931912B CN 104931912 B CN104931912 B CN 104931912B CN 201510323996 A CN201510323996 A CN 201510323996A CN 104931912 B CN104931912 B CN 104931912B
- Authority
- CN
- China
- Prior art keywords
- time domain
- data
- compensation
- response
- network analyzer
- 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.)
- Active
Links
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开了一种矢量网络分析仪的时域补偿方法,包括:通过傅里叶逆变换,将频域测得的散射参数变换到时域,分别得到反射参数和传输参数的时域数据;每个网络在时域数据中有不同的位置,根据各个网络在时域数据中的位置,构造时域选通函数;利用选通函数,对反射参数和传输参数的时域数据进行选通,得到各个网络的时域选通数据;将时域选通后的时域数据通过傅里叶变换,得到频域选通数据;构造补偿因子,利用遮蔽补偿公式,得到没有遮蔽误差的反射参数和传输参数。利用本发明,可以无误差地观测到微波电路中局部网络的散射参数。此外,本发明的方法简单直观,降低了操作人员的工作难度,而且具有更好的补偿结果。
Description
技术领域
本发明属于信号分析技术领域,尤其涉及一种矢量网络分析仪的时域补偿方法。
背景技术
矢量网络分析仪是微波工程中用于测量微波电路网络散射参数的电子测量仪器,而散射参数是微波工程中最常用的描述微波电路网络性能的参数指标。在微波电路中,一个电路系统往往是由多个微波网络级联而成,相互级联的网络之间相互耦合、相互影响,利用矢量网络分析仪只能测量微波电路系统整体的散射参数,而无法测量其中部分网络的散射参数。工作中,经常需要了解被测电路中部分网络的散射参数,但是这些散射参数一般会受到电路中其他部分的影响。为了消除各种级联的误差网络对测试结果造成的影响,人们提出了许多方法和理论,其一般过程是:首先,构造误差网络的模型;然后测量并计算得到误差网络模型中的各项参数;最后,在测量的结果中消除误差网络的影响,得到所需要网络的参数。
但是在实际应用中,不易构造一个未知误差网络的网络模型并确定模型中的各个参数;工程应用中该类方法一般采用经验模型,但经验模型很难达到较高的精度;这类方法不够直观,过程复杂,对操作人员要求较高。
发明内容
本发明的目的在于提供一种矢量网络分析仪的时域补偿方法,旨在解决现有矢量网络分析仪存在的不易构造一个未知网络的网络模型并确定模型中的各个参数,经验模型的精度较低,不够直观,过程复杂,对操作人员要求较高的问题。
本发明是这样实现的,一种矢量网络分析仪的时域补偿方法,所述矢量网络分析仪的时域补偿方法包括以下步骤:
步骤一,利用傅里叶逆变换,将频域测得的散射参数变换到时域,得到时域冲击响应,反射参数对应时域反射响应,传输参数对应时域传输响应;
步骤二,由于被测电路中的各个不连续点的物理位置是不同的,各个不连续点造成的时域发射响应和传输响应中的脉冲会分布在时间轴上的不同位置,根据反射响应和传输响应中前两个脉冲在时间轴上的位置,分别构造四个时域选通函数;
步骤三,利用步骤二中选通函数,对步骤一中时域的反射响应和传输响应进行选通,分别提取出反射响应和传输响应中的前两个脉冲;
步骤四,将时域选通后的时域脉冲分别通过傅里叶变换,得到频域选通数据;
步骤五,频域选通数据中包含了一些被测电路的信息,利用得到的信息,根据公式构造补偿因子Fcf(i);
步骤六,利用遮蔽补偿公式,得到没有遮蔽误差的反射参数FS11(i)和传输参数FS21(i)。
进一步,步骤五,利用下面的公式,构造补偿因子Fcf(i):
(i=1,2…)
G1(i)~G4(i)是在上步中得到的频域选通数据;
R(i)是比率因子;
Fcf(i)补偿因子。
进一步,步骤六,利用下面的遮蔽补偿公式,得到没有遮蔽误差的反射参数FS11(i)和传输参数FS21(i):
(i=1,2…)。
进一步,在步骤一之前需要进行以下步骤:
首先根据需要设置矢量网络分析仪的测量参数,得到被测网络整体的散射参数,包括反射参数和传输参数;
其次对散射参数进行预处理,对数据序列进行补零,进行快速傅里叶变换;根据傅里叶逆变换后时域脉冲的分布情况和分辨率的要求,选择不同的窗函数对傅里叶变换前的数据进行处理。
进一步,在步骤六之后需要消除数据预处理,舍弃添加零点位置处的数据,得到的反射参数和传输参数除以窗函数,消除加窗产生的影响,得到最终的补偿后的反射参数FS11_C(i)和传输参数FS21_C(i):
(i=1,2…)
W(i)数据预处理中的窗函数。
本发明提供的矢量网络分析仪的时域补偿方法,通过对数据加窗,可以减小快速傅里叶变换后的截断效应;通过补零,可以进行快速傅里叶变换,加快了数据的运算速度;预处理的目的是加快傅里叶逆变换和消除傅里叶逆变换后的振铃现象。为了加快傅里叶逆变换的速度,可以对数据序列进行补零,从而可以进行快速傅里叶变换;为了消除傅里叶变换后的振铃现象,根据傅里叶逆变换后时域脉冲的分布情况和分辨率的要求,可以选择不同的窗函数对傅里叶变换前的数据进行处理。本发明可以测量相互级联的微波电路中各部分网络的时域测量参数,通过对各个网络的时域参数分别进行选通操作,可以实现对被测微波电路中局部网络的散射参数的测量,并且消除电路中其他部分带来的误差。相比于传统方法,本发明的方法简单直观,降低操作人员的工作难度,而且具有更好的补偿结果。
附图说明
图1是本发明实施例提供的矢量网络分析仪的时域补偿方法流程图;
图2是本发明实施例提供的被测电路连接示意图;
图3是本发明实施例提供的作为参照实验的没有C2的参照电路示意图;
图4是本发明实施例提供的被测电路的反射参数处理结果比较示意图;
图5是本发明实施例提供的被测电路的传输参数处理结果比较示意图;
图6是本发明实施例提供的被测电路的S11参数处理结果比较示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
由于相互级联的网络在物理位置上一般有先后顺序,通过将矢量网络分析仪的频域测量结果转换到时域,可以观察到各个级联网络的时域测量结果在时间轴上是分离的,本发明利用这一特性,通过本发明的补偿算法,提出一种补偿微波电路中部分网络误差的算法。
下面结合附图1对本发明的补偿过程做详细的说明:
本发明实施例的矢量网络分析仪的时域补偿方法包括以下步骤:
步骤一,利用傅里叶逆变换,将频域测得的散射参数变换到时域,得到时域冲击响应,反射参数对应时域反射响应,传输参数对应时域传输响应;
步骤二,由于被测电路中的各个不连续点的物理位置是不同的,各个不连续点造成的时域发射响应和传输响应中的脉冲会分布在时间轴上的不同位置,根据反射响应和传输响应中前两个脉冲在时间轴上的位置,分别构造四个时域选通函数;
步骤三,利用步骤二中选通函数,对步骤一中时域的反射响应和传输响应进行选通,分别提取出反射响应和传输响应中的前两个脉冲;
步骤四,将时域选通后的时域脉冲分别通过傅里叶变换,得到频域选通数据;
步骤五,频域选通数据中包含了一些被测电路的信息,利用信息,根据公式构造补偿因子Fcf(i);
步骤六,将时域选通后的时域数据通过傅里叶变换,得到频域选通数据;
步骤七,利用下面的公式,构造补偿因子Fcf(i):
(i=1,2…)
G1(i)~G4(i)是在上步中得到的频域选通数据;
R(i)是比率因子;
Fcf(i)补偿因子;
步骤八,利用下面的遮蔽补偿公式,得到没有遮蔽误差的反射参数FS11(i)和传输参数FS21(i):
(i=1,2…);
步骤九,消除数据预处理带来的影响,舍弃添加零点位置处的数据;对上一步得到的反射参数和传输参数除以窗函数,消除加窗产生的影响,得到最终的补偿后的反射参数FS11_C(i)和传输参数FS21_C(i):
(i=1,2…)
W(i)数据预处理中的窗函数。
通过以下的实验对本发明的应用效果做详细的说明:
实验1:
起始频率为100MHz,频率间隔为100MHz,截止频率为8GHz,对如下图所示的有干扰网络进行测量。在实际网络中,假设电容C2是一个干扰,会对其他网络造成遮蔽误差;本实验的目的是使用本发明中的补偿算法,将C2电容的遮蔽补偿掉,与消除电容C2的理想网络的测试结果进行对比,验证本发明的使用效果。
根据测量参数,首先测量如图2所示的有干扰网络C2的散射参数,然后利用本发明算法,对干扰电路电容C2造成的遮蔽效应进行屏蔽,然后利用本发明算法进行补偿,并将补偿后的数据与作为参照实验的没有C2的参照电路的散射参数测量结果进行对照,通过观察反射参数和传输参数的曲线,可以发现没有C2的参照电路的散射参数测量结果和经过本发明算法补偿的数据具有很好的一致性,证明通过本发明提出的算法,可以有效的消除遮蔽效应带来的误差,如图4和图5所示。
实验2:
起始频率为100MHz,频率间隔为100MHz,截止频率为8GHz,对测试板进行测量,测试板的左右两侧焊接有电阻,可以加大反射,形成不连续点。在本实验中,通过在测试板的左侧焊接不同的电阻,然后利用本发明算法对实验板左侧的电阻造成的遮蔽进行补偿,将补偿后的数据进行对比。
第一次测试时,在实验板左侧位置焊接22欧姆电阻,右侧焊接10欧姆电阻。左侧的22欧姆电阻会对右侧10欧姆的电阻造成遮蔽效应,造成误差。利用本发明中的补偿算法对遮蔽效应进行补偿,并将补偿后的数据保存。
第二次测试时,将实验板左侧22欧姆的电阻换为任意其他阻值电阻,右侧电阻不改变,同样,左侧的电阻会对右侧10欧姆电阻造成遮蔽,利用本发明中的补偿算法进行补偿。
将两次补偿得到的结果进行对比,如图6所示,反射参数的第一次补偿数据与第二次补偿数据非常吻合,发现在实验板左侧任意阻值的电阻对后续电路造成的遮蔽,通过本发明的补偿可以消除其造成的遮蔽误差,验证了本发明的补偿算法可以有效的消除遮蔽效应带来的误差,证明了本发明的有效性和实用性。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种矢量网络分析仪的时域补偿方法,其特征在于,所述矢量网络分析仪的时域补偿方法包括以下步骤:
步骤一,利用傅里叶逆变换,将频域测得的散射参数变换到时域,得到时域冲击响应,反射参数对应时域反射响应,传输参数对应时域传输响应;
步骤二,根据反射响应和传输响应中前两个脉冲在时间轴上的位置,分别构造四个时域选通函数;
步骤三,利用步骤二中选通函数,对步骤一中时域的反射响应和传输响应进行选通,分别提取出反射响应和传输响应中的前两个脉冲;
步骤四,将时域选通后的时域脉冲分别通过傅里叶变换,得到频域选通数据;
步骤五,频域选通数据中包含了被测电路的信息,利用得到的信息,根据公式构造补偿因子Fcf(i);利用下面的公式,构造补偿因子Fcf(i):
(i=1,2…)
G1(i)~G4(i)是在上步中得到的频域选通数据;
R(i)是比率因子;
Fcf(i)补偿因子;
步骤六,利用遮蔽补偿公式,得到没有遮蔽误差的反射参数FS11(i)和传输参数FS21(i)。
2.如权利要求1所述的矢量网络分析仪的时域补偿方法,其特征在于,步骤六,利用下面的遮蔽补偿公式,得到没有遮蔽误差的反射参数FS11(i)和传输参数FS21(i):
(i=1,2…)。
3.如权利要求1所述的矢量网络分析仪的时域补偿方法,其特征在于,在步骤一之前需要进行以下步骤:
首先根据需要设置矢量网络分析仪的测量参数,得到被测网络整体的散射参数,包括反射参数和传输参数;
其次对散射参数进行预处理,对数据序列进行补零,进行快速傅里叶变换;根据傅里叶逆变换后时域脉冲的分布情况和分辨率的要求,选择不同的窗函数对傅里叶变换前的数据进行处理。
4.权利要求3所述的矢量网络分析仪的时域补偿方法,其特征在于,在步骤六之后需要消除数据预处理,舍弃添加零点位置处的数据,得到的反射参数和传输参数除以窗函数,消除加窗产生的影响,得到最终的补偿后的反射参数FS11_C(i)和传输参数FS21_C(i):
(i=1,2…)
W(i)数据预处理中的窗函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510323996.4A CN104931912B (zh) | 2015-06-12 | 2015-06-12 | 一种矢量网络分析仪的时域补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510323996.4A CN104931912B (zh) | 2015-06-12 | 2015-06-12 | 一种矢量网络分析仪的时域补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104931912A CN104931912A (zh) | 2015-09-23 |
CN104931912B true CN104931912B (zh) | 2018-07-24 |
Family
ID=54119165
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510323996.4A Active CN104931912B (zh) | 2015-06-12 | 2015-06-12 | 一种矢量网络分析仪的时域补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104931912B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107273634A (zh) * | 2017-07-03 | 2017-10-20 | 湖南工业大学 | 一种实时在线验证复杂交通控制方法 |
CN108387808A (zh) * | 2018-01-31 | 2018-08-10 | 杨彦青 | 一种智能电力变压器运行监测系统 |
CN108364381A (zh) * | 2018-01-31 | 2018-08-03 | 湖南城市学院 | 一种智能火情逃生门禁系统 |
CN108809165A (zh) * | 2018-06-29 | 2018-11-13 | 广东水利电力职业技术学院(广东省水利电力技工学校) | 一种交流伺服驱动器系统及控制方法 |
CN109077723A (zh) * | 2018-07-05 | 2018-12-25 | 青岛大学附属医院 | 基于远程控制的神经内科护理中的信息采集系统及方法 |
CN109032889A (zh) * | 2018-07-11 | 2018-12-18 | 广东水利电力职业技术学院(广东省水利电力技工学校) | 一种新型液冷式服务器系统及管理方法、计算机程序 |
CN109193670A (zh) * | 2018-08-16 | 2019-01-11 | 广东阿尔派智能电网有限公司 | 一种自供电电路的保护装置及方法 |
CN109243594A (zh) * | 2018-08-28 | 2019-01-18 | 青岛大学附属医院 | 一种小儿外科刀口清洁护理装置及控制方法 |
CN109528175A (zh) * | 2018-08-30 | 2019-03-29 | 吉林大学 | 一种基于互联网的妇科肿瘤护理中监控预警系统及预警方法 |
CN109173097A (zh) * | 2018-09-17 | 2019-01-11 | 陈永廷 | 一种肿瘤内科药物介入治疗装置及控制方法 |
CN109507078A (zh) * | 2018-11-01 | 2019-03-22 | 迈博瑞生物膜技术(南通)有限公司 | 无需外接电源的加压测试针头滤器泡点值的装置及方法 |
CN109541586A (zh) * | 2018-11-08 | 2019-03-29 | 山东交通学院 | 一种探地雷达快速探测道路工程混凝土垫层厚度分析方法 |
CN109800466A (zh) * | 2018-12-21 | 2019-05-24 | 内江师范学院 | 一种机场沥青混凝土道面变形过程的rbf模拟方法 |
CN109787923A (zh) * | 2019-02-13 | 2019-05-21 | 中联云港数据科技股份有限公司 | 一种基于大数据分析的idc+cdn+混合云的产品线布局方法 |
CN110361685B (zh) * | 2019-07-01 | 2021-07-30 | 北京无线电计量测试研究所 | 一种宽带示波器探头传输特性校准方法和系统 |
CN111624410B (zh) * | 2020-06-29 | 2023-03-28 | 普联技术有限公司 | 一种散射参数模型的获取方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1688890A (zh) * | 2002-11-05 | 2005-10-26 | 罗德施瓦兹两合股份有限公司 | 利用非正弦测量信号测量被测多端口设备散射参数的方法 |
CN1909527A (zh) * | 2006-08-18 | 2007-02-07 | 清华大学 | 带时域包络加权的正交频分复用系统的频域信道估计方法 |
CN103399286A (zh) * | 2013-07-16 | 2013-11-20 | 中国电子科技集团公司第四十一研究所 | 一种多特性阻抗网络的测量校准方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001069778A2 (de) * | 2000-03-16 | 2001-09-20 | Rohde & Schwarz Gmbh & Co. Kg | Verfahren zur bestimmung von parametern eines n-tors |
US7124049B2 (en) * | 2005-01-03 | 2006-10-17 | Agilent Technologies, Inc. | Method for implementing TRL calibration in VNA |
-
2015
- 2015-06-12 CN CN201510323996.4A patent/CN104931912B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1688890A (zh) * | 2002-11-05 | 2005-10-26 | 罗德施瓦兹两合股份有限公司 | 利用非正弦测量信号测量被测多端口设备散射参数的方法 |
CN1909527A (zh) * | 2006-08-18 | 2007-02-07 | 清华大学 | 带时域包络加权的正交频分复用系统的频域信道估计方法 |
CN103399286A (zh) * | 2013-07-16 | 2013-11-20 | 中国电子科技集团公司第四十一研究所 | 一种多特性阻抗网络的测量校准方法 |
Non-Patent Citations (3)
Title |
---|
"矢量网络分析仪中时域测试功能的应用";白莉萍;《科技信息》;20090525;全文 * |
"矢量网络分析仪时域功能分析";徐杰 等;《电子测量技术》;20080915;第31卷(第9期);第35-37页 * |
"矢量网络分析仪时域功能及应用";刘丽 等;《中国测试》;20120330;第38卷(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104931912A (zh) | 2015-09-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104931912B (zh) | 一种矢量网络分析仪的时域补偿方法 | |
Karimzadeh et al. | Method for determining utility and consumer harmonic contributions based on complex independent component analysis | |
CN105866543B (zh) | 一种消除基波、谐波对间谐波检测干扰的间谐波检测方法 | |
CN108469560A (zh) | 一种基于快速s变换时频空间模型的电磁干扰客观复杂度评估方法 | |
US20110238383A1 (en) | One-Port De-embedding Using Time Domain Substitution | |
CN105354397B (zh) | 一种电机驱动系统共模电磁干扰滤波器的设计方法 | |
Mandrekar et al. | Causality enforcement in transient simulation of passive networks through delay extraction | |
CN108387808A (zh) | 一种智能电力变压器运行监测系统 | |
CN203643514U (zh) | 用于测量散射参数的耦合组件和数据测定装置 | |
Cho et al. | An efficient extrapolation method of band-limited S-parameters for extracting causal impulse responses | |
Zhang et al. | Micro-motion frequency estimation of radar targets with complicated translations | |
CN105277819A (zh) | 抗噪性评估装置和评估抗噪性的方法 | |
Gazda et al. | Harmonic balance surrogate-based immunity modeling of a nonlinear analog circuit | |
Liu et al. | Application of pulse compression technique in fault detection and localization of leaky coaxial cable | |
US20110234239A1 (en) | Two-Port De-Embedding Using Time Domain Substitution | |
Weber et al. | Cable fault characterization by time-domain analysis from S-parameter measurement and sparse inverse Chirp-Z transform | |
Tan et al. | Harmonic analysis based on time domain mutual-multiplication window | |
Rämmal et al. | Modified multi-load method for nonlinear source characterisation | |
Araneo et al. | Two-port equivalent of PCB discontinuities in the wavelet domain | |
Meng et al. | NMR log data de-noising method based on a variable order wavelet packet domain adaptive filtering | |
Lin | Development of interharmonics identification using enhanced‐fft algorithm | |
Maier et al. | Frequency-domain cable model with a real representation of joints and smooth failure changes and a diagnostic investigation with line resonance analysis | |
Schutt-Ainé et al. | Comparing fast convolution and model order reduction methods for S-parameter simulation | |
Bahena et al. | Impact of neglecting the metal losses on the extraction of the relative permittivity from PCB transmission line measurements | |
Razavi et al. | Measurement-based circuit modeling of conducted and radiated automotive immunity test setups |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |