CN116299110A - 一种用于双指数波形式电磁脉冲测量的时域标定方法 - Google Patents

一种用于双指数波形式电磁脉冲测量的时域标定方法 Download PDF

Info

Publication number
CN116299110A
CN116299110A CN202211530339.3A CN202211530339A CN116299110A CN 116299110 A CN116299110 A CN 116299110A CN 202211530339 A CN202211530339 A CN 202211530339A CN 116299110 A CN116299110 A CN 116299110A
Authority
CN
China
Prior art keywords
waveform
pulse source
measured
standard pulse
electromagnetic pulse
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.)
Pending
Application number
CN202211530339.3A
Other languages
English (en)
Inventor
杨静
崔志同
朱志臻
董亚运
邱孟通
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwest Institute of Nuclear Technology
Original Assignee
Northwest Institute of Nuclear Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwest Institute of Nuclear Technology filed Critical Northwest Institute of Nuclear Technology
Priority to CN202211530339.3A priority Critical patent/CN116299110A/zh
Publication of CN116299110A publication Critical patent/CN116299110A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R35/00—Testing or calibrating of apparatus covered by the other groups of this subclass
    • 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
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Abstract

本发明涉及一种用于双指数波形式电磁脉冲测量的时域标定方法，解决标定用标准脉冲源波形选取问题。包括步骤1：获得电磁脉冲测量系统的频响曲线H(jω)，并构建待测波形数据库；2：计算每个待测波形经过电磁脉冲测量系统后，输出波形的峰值；3：获得可选标准脉冲源波形库；4：计算采用每个可选标准脉冲源波形进行标定时，传输系数的理论值；5：计算采用该传输系数对待测波形数据库中每个待测波形进行测量时的误差，并获得每个可选标准脉冲源波形对应的测量结果的误差界，选取最小的误差界对应的可选标准脉冲源波形作为传输系数标定用标准脉冲源波形sz(t)；6：利用sz(t)对电磁脉冲测量系统进行时域标定。

Description

一种用于双指数波形式电磁脉冲测量的时域标定方法
技术领域
本发明涉及用于电磁脉冲测量系统的标定方法，具体涉及一种用于双指数波形式电磁脉冲测量的时域标定方法。
背景技术
电磁脉冲测量系统主要是指脉冲电场和脉冲电流测量系统。测量系统在使用前需进行标定。标定的思路是通过标准源来对电磁脉冲测量系统的某些参数进行标校。对于电磁脉冲测量系统，根据标准源的特征，标定方法可分为时域标定方法和频域标定方法。时域标定的标准源产生脉冲信号，频域标定的标准源产生的是一系列正弦信号。通过分析测量系统对标准信号的响应分析测量系统特性。时域标定常用于电磁脉冲测量系统的动态测量范围、线性区间，传输系数、响应时间等参数的标校，频域标定用于电磁脉冲测量系统的频率响应，带宽等的标校。
电磁脉冲测量系统带宽较宽，在工作频带内频率响应波动小于±3dB。因此当被测信号频率成分在测量系统通带范围内时，可认为测量系统输出波形与被测信号波形相同，仅幅值不同。被测信号与测量系统输出幅值之间的比值即为传输系数。因此，可通过时域标定获得测量系统的传输系数，采用该测量系统测量脉冲信号时，将测量系统输出乘以标定得到的传输系数即可方便地恢复被测信号。目前的时域标定常采用方波标准源开展。若待测信号为双指数形式的电磁脉冲，由于标准脉冲源波形与待测双指数波形式的电磁脉冲的频率成分存在差异，且电磁脉冲测量系统频率响应通带内存在波动，导致采用传输系数恢复待测波形时存在误差，目前还没有对这个误差进行定量分析和优化的方法。
发明内容
本发明的目的是解决由于标准脉冲源波形与待测双指数波形式的电磁脉冲的频率成分存在差异，且电磁脉冲测量系统频率响应通带内存在波动，导致采用传输系数恢复待测波形时存在误差，目前还没有对这个误差进行定量分析和优化方法的问题，提供了一种用于双指数波形式电磁脉冲测量的时域标定方法，该方法的特点在于根据待测双指数波形式电磁脉冲信号范围和电磁脉冲测量系统频率响应优化了传输系数标定用标准脉冲源波形，从而使得由标定脉冲源波形引起的测量误差最小。
本发明的设计思路为：
标准脉冲源波形与待测双指数波形式电磁脉冲频谱差异以及电磁脉冲测量系统频率响应通带内的波动，导致采用标定的传输系数恢复待测信号时存在误差，本发明中根据待测双指数波形式电磁脉冲信号特征及范围和电磁脉冲测量系统的频率响应确定时域标定用标准脉冲源波形，并通过定义和计算不同标准脉冲源波形引起的测量误差界，在误差界最小条件下，确定时域标定用标准脉冲源波形，使得由标定脉冲源波形引起的测量误差最小。
本发明所采用的技术方案是：
一种用于双指数波形式电磁脉冲测量的时域标定方法，其特殊之处在于，包括以下步骤：
步骤1：获得电磁脉冲测量系统的频响曲线H(jω)；以及根据待测双指数波形式电磁脉冲前沿tr的范围[Q，W]和半宽tw的范围[R，F]，构建待测波形数据库，其中，待测波形数据库内包含M个待测波形，j为虚数，ω为角频率；
步骤2：基于步骤1中的频响曲线H(jω)和待测波形数据库，计算每个待测波形经过电磁脉冲测量系统后，输出波形的峰值；
步骤3：基于待测双指数波形式电磁脉冲的前沿tr的范围[Q，W]和半宽tw的范围[R，F]，形成可选标准脉冲源波形库，可选标准脉冲源波形库中包含N个可选标准脉冲源波形；
步骤4：计算采用每个可选标准脉冲源波形对电磁脉冲测量系统进行标定时，传输系数的理论值；
步骤5：基于步骤4中每个传输系数的理论值以及步骤2中每个输出波形的峰值，计算采用该传输系数对待测波形数据库中每个待测波形进行测量时的误差，并获得每个可选标准脉冲源波形对应的测量结果的误差界，选取最小的误差界对应的可选标准脉冲源波形作为传输系数标定用标准脉冲源波形sz(t)；
步骤6：利用步骤5获得的标准脉冲源波形sz(t)对电磁脉冲测量系统进行时域标定。
进一步地，步骤1中，所述根据待测双指数波形式电磁脉冲前沿tr的范围[Q，W]和半宽tw的范围[R，F]，构建待测波形数据库的方式如下：
设定前沿步进步长为U，半宽步进步长为P，则在前沿tr的范围[Q，W]和半宽tw的范围[R，F]内，可以形成((W-Q)/U+1)((F-R)/P+1)个待测波形，((W-Q)/U+1)((F-R)/P+1)个待测波形构成待测波形数据库。
进一步地，所述步骤1中，还包括对待测波形数据库进行预处理，获得极性统一和幅值归一化处理后的离散待测波形数据库，具体为：
通过公式
Figure BDA0003974196940000041
对待测波形数据库中的每个待测波形进行极性统一和幅值归一化，其中，xxl(n)表示预处理前待测波形数据库中，第l个待测波形的采样，n代表时间序列，xl(n)表示预处理后的第l个待测波形，l＝1,2......M。
进一步地，所述步骤2具体包括以下步骤：
2.1计算待测波形数据库中，每一个待测波形经过传输函数为H(jω)的电磁脉冲测量系统后的输出波形；
输出波形的计算公式如下：
Figure BDA0003974196940000042
其中
Figure BDA0003974196940000043
分别表示离散傅里叶变换和离散傅里叶反变换，yl(n)表示第l个待测波形经过传输函数为H(jω)的电磁脉冲测量系统后的输出波形；
2.2基于步骤2.1中获得的M个输出波形，计算每一个输出波形的峰值；
通过公式pl＝||yl(n)||∞计算每一个输出波形的峰值，其中，pl表示第l个输出波形的峰值。
进一步地，所述步骤3中，还包括对可选标准脉冲源波形库进行预处理，获得极性统一和幅值归一化处理后的可选标准脉冲源波形库；
通过公式
Figure BDA0003974196940000044
对可选标准脉冲源波形库中的每个可选标准脉冲源波形进行极性统一和幅值归一化，其中，ssi(n)表示预处理前的第i个可选标准脉冲源波形的采样，si(n)表示预处理后的第i个可选标准脉冲源波形，i＝1,2......N,n代表时间序列。
进一步地，所述步骤4具体包括以下步骤：
4.1计算每一个可选标准脉冲源波形经过传输函数为H(jω)的电磁脉冲测量系统后的输出波形；
计算公式如下：
Figure BDA0003974196940000051
其中，
Figure BDA0003974196940000052
分别表示离散傅里叶变换和离散傅里叶反变换，syi(n)表示第i个可选标准脉冲源波形si(n)经过传输函数为H(jω)的系统后的输出波形；
4.2基于步骤3中，可选标准脉冲源波形库中的每一个可选标准脉冲源波形及步骤4.1中每一个可选标准脉冲源波形经过传输函数为H(jω)的系统后的输出波形，分别计算每一个可选标准脉冲源波形的峰值及其输出波形的峰值；
第i个可选标准脉冲源波形的峰值Si为：Si＝||si(n)||∞；
第i个可选标准脉冲源波形经过传输函数为H(jω)的系统后，输出波形的峰值SYi为：SYi＝||syi(n)||∞；
4.3基于步骤4.2中，每一个可选标准脉冲源波形的峰值及其输出波形的峰值，计算得到每一个可选标准脉冲源波形对应的传输系数理论值；
第i个可选标准脉冲源波形对应的传输系数理论值ki为：
Figure BDA0003974196940000053
进一步地，所述步骤5包括以下步骤：
5.1以每一个可选标准脉冲源波形对传输系数进行标定后，采用电磁脉冲测量系统对待测波形数据库内的每个待测波形进行测量，得到采用该传输系数得到的测量结果的误差界；
计算公式如下：
δi＝1-kipl∞
其中，δi表示采用第i个可选标准脉冲源波形对电磁脉冲测量系统的传输系数进行标定，电磁脉冲测量系统对待测波形库内信号进行测量，得到的测量结果的误差界；
5.2基于步骤5.1中每个标准脉冲源波形的误差界，选取最小的误差界对应的可选标准脉冲源的波形作为传输系数标定用标准脉冲源波形。
本发明的有益效果是：
1、本发明充分利用电磁脉冲测量系统的频率响应和待测双指数波形式电磁脉冲的先验知识(前沿tr的范围[Q，W]和半宽tw的范围[R，F])，从理论上定量给出了由标定用可选标准脉冲源波形引起的测量结果误差界，并选取最小的误差界对应的可选标准脉冲源波形作为传输系数标定用标准脉冲源波形，从而使得由标定脉冲源波形引起的测量误差最小。
2、本发明基于误差界提出的一种用于双指数波形式电磁脉冲测量的时域标定方法，此方法给出了进行传输系数标定时标准脉冲源波形选取方法，为标准脉冲源的选取提供理论依据。
3、本发明提出的方法可集成至标准脉冲源控制程序，依据测量系统频响特性和待测双指数波形式电磁脉冲的范围，调节可选标准脉冲源波形的输出波形，提供高质量的标定结果，为电磁脉冲测量系统标定平台建设提供参考。
4、本发明提出的时域标定方法可推广至其他形式电磁脉冲信号测量。
附图说明
图1是本发明实施例中，获得电磁脉冲电流测量系统的频率响应的试验设置；
图2(a)是本发明实施例中，电磁脉冲电流测量系统的频率与振幅关系图；
图2(b)是本发明实施例中，电磁脉冲电流测量系统的频率与相移关系图；
图3是本发明实施例中，双指数波形式的电磁脉冲波形图；
图4是本发明实施例中，待测波形数据库内的所有待测波形经过电磁脉冲测量系统后输出幅值的分布图，其中P表示幅值；
图5是本发明实施例中，脉冲电流测量系统的时域标定试验装置示意图；
图6本发明用于双指数波形式电磁脉冲测量的时域标定方法实施例的流程图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明提出一种用于双指数波形式电磁脉冲测量的时域标定方法，标定方法主要改进点在于：依据待测双指数波形式电磁脉冲范围和电磁脉冲测量系统的频率响应，优化传输系数时域标定用标准脉冲源波形，如图6所示，具体包括以下步骤：
步骤1：获得电磁脉冲测量系统的频响曲线H(jω)；以及根据待测双指数波形式电磁脉冲前沿tr的范围[Q，W]和半宽tw的范围[R，F]，构建待测波形数据库，其中，待测波形数据库内包含M个待测波形，H表示传输函数，j为虚数，ω为角频率，在本实施例中，tr和tw的单位为ns；
获得电磁脉冲测量系统频响曲线H(jω)的方法如下；
可以通过产品说明书获得电磁脉冲测量系统的频响曲线H(jω)；
或者通过频域标定结果获得电磁脉冲测量系统的频响曲线H(jω)，且所使用的电磁脉冲测量系统相关参数满足测试需求；
通过频域标定结果获得电磁脉冲测量系统的频响曲线H(jω)的方法具体为：如选取的电磁脉冲测量系统为一套电磁脉冲电流测量系统，通过图1的连接方式，对其进行频域标定，通过测量结果计算得到该电磁脉冲测量系统的频响曲线H(jω)，如图2(a)与图2(b)所示，该电磁脉冲测量系统的±3dB频带宽度约为200MHz；
B、如图3所示，根据待测双指数波形式电磁脉冲前沿tr的范围[Q，W]和半宽tw的范围[R，F]，构建待测波形数据库的方式如下：
设定前沿步进步长为U，半宽步进步长为P，则在前沿tr的范围[Q，W]和半宽tw的范围[R，F]内，可以形成((W-Q)/U+1)((F-R)/P+1)个待测波形，((W-Q)/U+1)((F-R)/P+1)个待测波形构成待测波形数据库，即M＝((W-Q)/U+1)((F-R)/P+1)，在本实施例中，U和P的单位为ns；
例如假设待测双指数波形式电磁脉冲前沿tr的范围为[5ns，20ns]，半宽tw的范围为[30ns，200ns]，设定前沿步进步长为0.5ns，半宽步进步长为1ns，在前沿tr的范围为[5ns，20ns]，半宽tw的范围为[30ns，200ns]的前提下，可以形成M＝31×171＝5301个待测波形，即M＝5301；待测波形数据库A是根据用户测试需求生成的；
为方便后续处理，将待测波形数据库内的所有待测波形进行极性统一和幅值归一化预处理，获得预处理后的待测波形数据库，获得预处理后的待测波形数据库方式如下：
通过下列公式计算预处理后的待测波形：
Figure BDA0003974196940000091
其中，xxl(n)表示预处理前待测波形数据库中，第l个待测波形，n代表时间，xl(n)表示预处理后的第l个待测波形，l＝1,2......M；
步骤2：基于步骤1中的频响曲线H(jω)和待测波形数据库，计算每个待测波形经过电磁脉冲测量系统后，输出波形的峰值，结果如图4所示，图4中，水平面放入横轴和纵轴分别为待测波形的前沿和半宽范围，竖直坐标表示峰值；
2.1计算步骤1获得的待测波形数据库中，每一个待测波形经过传输函数为H(jω)的电磁脉冲测量系统后的输出波形，输出波形的数量为M个；
输出波形的计算公式如下：
Figure BDA0003974196940000092
其中
Figure BDA0003974196940000093
分别表示离散傅里叶变换和离散傅里叶反变换，yl(n)表示第l个待测波形经过传输函数为H(jω)的电磁脉冲测量系统后的输出波形；
2.2基于步骤2.1中获得的M个输出波形，计算每一个输出波形的峰值；
通过公式pl＝||yl(n)||∞计算每一个输出波形的峰值，其中，pl表示第l个输出波形的峰值；
步骤3：基于待测双指数波形式电磁脉冲的前沿tr的范围[Q，W]和半宽tw的范围[R，F]，形成可选标准脉冲源波形库，可选标准脉冲源波形库中中包含N个可用于测量系统传输系数标定的标准脉冲源波形；
假设可选标准脉冲源波形有三个，即N＝3，第一个为前沿tr＝5ns，半宽tw＝30ns的双指数波，第二个为前沿tr＝10ns，半宽tw＝100ns的双指数波，第三个为前沿tr＝20ns，半宽tw＝200ns的双指数波；
为方便后续处理，对可选标准脉冲源波形库内的N个可选标准脉冲源波形进行极性统一和幅值归一化预处理，获得由预处理后N个可选标准脉冲源波形组成的可选标准脉冲源波形库；
通过下列公式计算预处理后的可选标准脉冲源波形；
Figure BDA0003974196940000101
其中，ssi(t)表示预处理前的第i个可选标准脉冲源波形的采样，采样需满足奈奎斯特采样定律，si(t)表示预处理后的第i个可选标准脉冲源波形，i＝1,2......N，N＝3；
步骤4：基于步骤3中的可选标准脉冲源波形库，计算采用可选标准脉冲源波形库中，每个可选标准脉冲源波形对电磁脉冲测量系统进行标定时，传输系数的理论值；
4.1如图5所示，计算每个可选标准脉冲源波形经过传输函数为H(jω)的电磁脉冲测量系统后的输出波形；
计算公式如下：
Figure BDA0003974196940000102
其中，
Figure BDA0003974196940000103
分别表示傅里叶变换和傅里叶反变换，syi(n)表示第i个可选标准脉冲源波形si(n)经过传输函数为H(jω)的系统后的输出波形；
4.2基于步骤3可选标准脉冲源波形库中的每一个可选标准脉冲源波形及步骤4.1中每一个可选标准脉冲源波形经过传输函数为H(jω)的系统后的输出波形，分别计算每一个可选标准脉冲源波形的峰值及其输出波形的峰值；
第i个可选标准脉冲源波形的峰值Si为：Si＝||si(n)||∞；
第i个可选标准脉冲源波形经过传输函数为H(jω)的系统后，输出波形的峰值SYi为：SYi＝||syi(n)||∞；
4.3通过公式
Figure BDA0003974196940000111
计算得到每一个可选标准脉冲源波形对应的传输系数的理论值k；
第i个标准脉冲源对应的传输系数ki的计算公式为：
Figure BDA0003974196940000112
由于si(t)是经过归一化预处理后的第i个可选标准脉冲源波形，所以Si＝1；
在本实施例中，计算出来的第一种可选标准脉冲源波形对应的传输系数理论值k1＝1.0781,第二种可选标准脉冲源波形对应的传输系数理论值k2＝1.0520,第三种可选标准脉冲源波形对应的传输系数理论值k3＝1.047；
步骤5：基于步骤4中每个传输系数的理论值以及步骤2中每个输出波形的峰值，计算采用该传输系数对待测波形数据库中每个待测波形进行测量时的误差，并获得每个可选标准脉冲源波形对应的测量结果的误差界，选取最小的误差界对应的可选标准脉冲源波形作为传输系数标定用标准脉冲源波形sz(t)；
5.1以每一个可选标准脉冲源波形对传输系数进行标定后，采用电磁脉冲测量系统对待测波形数据库内的每个待测波形进行测量，得到采用该传输系数得到的测量结果的误差界；
计算公式如下：
δi＝||1-kipl||∞
其中，δi表示采用第i个可选标准脉冲源波形对电磁脉冲测量系统的传输系数进行标定，电磁脉冲测量系统对待测波形库内信号进行测量，得到的测量结果的误差界；
具体的：计算得到三种可选标准脉冲源对应的误差范围依次分别为[-2.9％0％]，[-0.4％2.4％]，[02.8％]，则对应的误差界为δ1＝2.9％，δ2＝2.4％，δ3＝2.8％；
5.2基于步骤5.1中每个标准脉冲源波形的误差界，选取最小的误差界对应的可选标准脉冲源的波形作为传输系数标定用标准脉冲源波形；
其中，误差界最小的为δ2＝2.4％，则sz(t)＝s2(t)，即本发明选取的标准脉冲源波形为前沿tr＝10ns，半宽tw＝100ns的双指数波。
步骤6：利用步骤5获得的传输系数标定用标准脉冲源波形sz(t)，结合测量系统测量范围选定标准脉冲源，对双指数波形式电磁脉冲测量系统进行时域标定；
根据标准脉冲源在示波器输出Vsource和电流传感器输出Vout标定得到时域传输系数kcal
Figure BDA0003974196940000121

Claims (7)

1.一种用于双指数波形式电磁脉冲测量的时域标定方法，其特征在于，包括以下步骤：
步骤1：获得电磁脉冲测量系统的频响曲线H(jω)；以及根据待测双指数波形式电磁脉冲前沿tr的范围[Q，W]和半宽tw的范围[R，F]，构建待测波形数据库，其中，待测波形数据库内包含M个待测波形，j为虚数，ω为角频率；
步骤2：基于步骤1中的频响曲线H(jω)和待测波形数据库，计算每个待测波形经过电磁脉冲测量系统后，输出波形的峰值；
步骤3：基于待测双指数波形式电磁脉冲的前沿tr的范围[Q，W]和半宽tw的范围[R，F]，形成可选标准脉冲源波形库，可选标准脉冲源波形库中包含N个可选标准脉冲源波形；
步骤4：计算采用每个可选标准脉冲源波形对电磁脉冲测量系统进行标定时，传输系数的理论值；
步骤5：基于步骤4中每个传输系数的理论值以及步骤2中每个输出波形的峰值，计算采用该传输系数对待测波形数据库中每个待测波形进行测量时的误差，并获得每个可选标准脉冲源波形对应的测量结果的误差界，选取最小的误差界对应的可选标准脉冲源波形作为传输系数标定用标准脉冲源波形sz(t)；
步骤6：利用步骤5获得的标准脉冲源波形sz(t)对电磁脉冲测量系统进行时域标定。
2.根据权利要求1所述的一种用于双指数波形式电磁脉冲测量的时域标定方法，其特征在于：
步骤1中，所述根据待测双指数波形式电磁脉冲前沿tr的范围[Q，W]和半宽tw的范围[R，F]，构建待测波形数据库的方式如下：
设定前沿步进步长为U，半宽步进步长为P，则在前沿tr的范围[Q，W]和半宽tw的范围[R，F]内，可以形成((W-Q)/U+1)((F-R)/P+1)个待测波形，((W-Q)/U+1)((F-R)/P+1)个待测波形构成待测波形数据库。
3.根据权利要求2所述的一种用于双指数波形式电磁脉冲测量的时域标定方法，其特征在于：
所述步骤1中，还包括对待测波形数据库进行预处理，获得极性统一和幅值归一化处理后的离散待测波形数据库，具体为：
通过公式
Figure FDA0003974196930000021
对待测波形数据库中的每个待测波形进行极性统一和幅值归一化，其中，xxl(n)表示待测波形数据库中，第l个待测波形的采样，n代表时间序列，xl(n)表示预处理后的第l个待测波形，l＝1,2......M。
4.根据权利要求3所述的一种用于双指数波形式电磁脉冲测量的时域标定方法，其特征在于：
所述步骤2具体包括以下步骤：
2.1计算待测波形数据库中，每一个待测波形经过传输函数为H(jω)的电磁脉冲测量系统后的输出波形；
输出波形的计算公式如下：
Figure FDA0003974196930000022
其中
Figure FDA0003974196930000023
分别表示离散傅里叶变换和离散傅里叶反变换，yl(n)表示第l个待测波形经过传输函数为H(jω)的电磁脉冲测量系统后的输出波形；
2.2基于步骤2.1中获得的M个输出波形，计算每一个输出波形的峰值；
通过公式pl＝yl(n)∞计算每一个输出波形的峰值，其中，pl表示第l个输出波形的峰值。
5.根据权利要求4所述的一种用于双指数波形式电磁脉冲测量的时域标定方法，其特征在于：
所述步骤3中，还包括对可选标准脉冲源波形库进行预处理，获得极性统一和幅值归一化处理后的可选标准脉冲源波形库；
通过公式
Figure FDA0003974196930000031
对可选标准脉冲源波形库中的每个可选标准脉冲源波形进行极性统一和幅值归一化，其中，ssi(n)表示预处理前的第i个可选标准脉冲源波形的采样，si(n)表示预处理后的第i个可选标准脉冲源波形，i＝1,2......N，n代表时间序列。
6.根据权利要求5所述的一种用于双指数波形式电磁脉冲测量的时域标定方法，其特征在于：
所述步骤4具体包括以下步骤：
4.1计算每一个可选标准脉冲源波形经过传输函数为H(jω)的电磁脉冲测量系统后的输出波形；
计算公式如下：
Figure FDA0003974196930000032
其中，
Figure FDA0003974196930000033
分别表示离散傅里叶变换和离散傅里叶反变换，syi(n)表示第i个可选标准脉冲源波形si(n)经过传输函数为H(jω)的系统后的输出波形；
4.2基于步骤3中，可选标准脉冲源波形库中的每一个可选标准脉冲源波形及步骤4.1中每一个可选标准脉冲源波形经过传输函数为H(jω)的系统后的输出波形，分别计算每一个可选标准脉冲源波形的峰值及其输出波形的峰值；
第i个可选标准脉冲源波形的峰值Si为：Si＝si(n)∞；
第i个可选标准脉冲源波形经过传输函数为H(jω)的系统后，输出波形的峰值SYi为：SYi＝syi(n)∞；
4.3基于步骤4.2中，每一个可选标准脉冲源波形的峰值及其输出波形的峰值，计算得到每一个可选标准脉冲源波形对应的传输系数理论值；
第i个可选标准脉冲源波形对应的传输系数理论值ki为：
Figure FDA0003974196930000041
7.根据权利要求6所述的一种用于双指数波形式电磁脉冲测量的时域标定方法，其特征在于：
所述步骤5包括以下步骤：
5.1以每一个可选标准脉冲源波形对传输系数进行标定后，采用电磁脉冲测量系统对待测波形数据库内的每个待测波形进行测量，得到采用该传输系数得到的测量结果的误差界；
计算公式如下：
δi＝1-kipl∞
其中，δi表示采用第i个可选标准脉冲源波形对电磁脉冲测量系统的传输系数进行标定，电磁脉冲测量系统对待测波形库内信号进行测量，得到的测量结果的误差界；
5.2基于步骤5.1中每个标准脉冲源波形的误差界，选取最小的误差界对应的可选标准脉冲源的波形作为传输系数标定用标准脉冲源波形。
CN202211530339.3A 2022-11-30 2022-11-30 一种用于双指数波形式电磁脉冲测量的时域标定方法 Pending CN116299110A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211530339.3A CN116299110A (zh) 2022-11-30 2022-11-30 一种用于双指数波形式电磁脉冲测量的时域标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211530339.3A CN116299110A (zh) 2022-11-30 2022-11-30 一种用于双指数波形式电磁脉冲测量的时域标定方法

Publications (1)

Publication Number Publication Date
CN116299110A true CN116299110A (zh) 2023-06-23

Family

ID=86798379

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211530339.3A Pending CN116299110A (zh) 2022-11-30 2022-11-30 一种用于双指数波形式电磁脉冲测量的时域标定方法

Country Status (1)

Country Link
CN (1) CN116299110A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116609715A (zh) * 2023-05-24 2023-08-18 中国人民解放军火箭军工程大学 一种双指数波型脉冲形状参数测量误差估计方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116609715A (zh) * 2023-05-24 2023-08-18 中国人民解放军火箭军工程大学 一种双指数波型脉冲形状参数测量误差估计方法

Similar Documents

Publication Publication Date Title
US6230106B1 (en) Method of characterizing a device under test
US7271575B2 (en) Oscilloscope based return loss analyzer
CN100561162C (zh) 一种虚拟振动台检测信号处理方法及其设备
EP2290392B1 (en) Enhanced impulse response measurement of an FM radar transmitter pulse
US10914782B2 (en) Test system and test method for testing a device under test
US20110238383A1 (en) One-Port De-embedding Using Time Domain Substitution
CN110133564B (zh) 一种电流传感器相位特性的同步测试方法及装置
US8706433B2 (en) Time domain reflectometry step to S-parameter conversion
US10191092B2 (en) Time domain measuring method with calibration in the frequency range
CN116299110A (zh) 一种用于双指数波形式电磁脉冲测量的时域标定方法
EP3051709B1 (en) De-embedding cable effect for waveform monitoring for arbitrary waveform and function generator
US20090125264A1 (en) Nonlinear Measurement System Error Correction
GB2331589A (en) Network analyser measurement method
US10509064B2 (en) Impedance measurement through waveform monitoring
US20110234239A1 (en) Two-Port De-Embedding Using Time Domain Substitution
JP4141961B2 (ja) システム校正されたベクトル・ネットワークアナライザの実効指向性および／または実効ソースポート整合性の測定方法、校正標準セットの作成方法
US20070253562A1 (en) Method and apparatus for measuring characteristics of an audio system using a tapered chirp
CN101308174A (zh) 采用可变带宽参数的小波变换的信号频谱测量分析方法
EP3179259B1 (en) Group delay based averaging
CN115792764A (zh) 一种传输系数法测量电磁脉冲系统误差的分析方法
Blair Error estimates for frequency responses calculated from time-domain measurements
JP3045759B2 (ja) セトリング特性測定方法
KR102036790B1 (ko) 계단-주파수 반사파 계측법을 이용하는 케이블 결함 평가 시스템 및 방법
CN112033518A (zh) 一种基于正弦函数拟合的水听器灵敏度标定方法
CN111555995A (zh) 一种宽带采集系统的相频响应测量及补偿方法

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