CN107329002A - 强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法 - Google Patents
强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法 Download PDFInfo
- Publication number
- CN107329002A CN107329002A CN201710639671.6A CN201710639671A CN107329002A CN 107329002 A CN107329002 A CN 107329002A CN 201710639671 A CN201710639671 A CN 201710639671A CN 107329002 A CN107329002 A CN 107329002A
- Authority
- CN
- China
- Prior art keywords
- frequency
- mrow
- signal
- time domain
- shield
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/0807—Measuring electromagnetic field characteristics characterised by the application
- G01R29/0814—Field measurements related to measuring influence on or from apparatus, components or humans, e.g. in ESD, EMI, EMC, EMP testing, measuring radiation leakage; detecting presence of micro- or radiowave emitters; dosimetry; testing shielding; measurements related to lightning
- G01R29/0835—Testing shielding, e.g. for efficiency
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Shielding Devices Or Components To Electric Or Magnetic Fields (AREA)
Abstract
本发明涉及一种强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法，其特征是：至少包括如下步骤：S1，扫频发射测试屏蔽体的频域屏蔽效能；S2，将强电磁脉冲源时域波形通过傅立叶变换，获得频域谱特性；S3，在脉冲源场强频谱特性基础上，使其相应频点按S1中获得的屏蔽效能进行相应衰减，得到经屏蔽体衰减后的场强频谱特性；S4，将S3中获得的场强频谱特性，进行傅里叶逆变换获得对应脉冲源经屏蔽作用后的时域波形；通过对屏蔽前后时域波形的峰值之比，得到对应脉冲源作用下的时域屏蔽效能。本发明是具有操作简单方便、成本低，能获取真实强电磁脉冲环境下屏蔽体的防护抑制效果，分析屏蔽体对强电磁脉冲的防护能力。
Description
技术领域
本发明属于电磁脉冲防护领域，提出了一种强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法。
背景技术
随着以半导体类器件为核心的各类电子元器件的性能向高频、高速方向发展，使得电子设备和系统对电磁骚扰的敏感度越来越低，电磁毁伤阈值也越来越低。越来越多的电子产品正面临着强电磁脉冲的严重威胁，电磁脉冲可通过孔缝耦合进入电子系统。对于孔口、缝隙耦合通常采取的强电磁脉冲的防护措施为屏蔽，即在空间上将辐射源与敏感设备进行隔离，阻断其传播路径。
屏蔽体对电子设备或者系统的强电磁脉冲防护加固效果，取决于其对脉冲信号瞬时峰值的抑制量级，通常由电磁脉冲时域屏蔽效能来衡量，即通过测试得到电磁脉冲屏蔽前后的时域波形，取其峰值之比计算得出。一般的单位并不具备专用电磁脉冲时域测量系统，并且电磁脉冲时域测试系统非常复杂，强电磁脉冲源种类多(例如核电磁脉冲、高功率微波、雷电等)，测试费用昂贵。因此研究一种屏蔽体的时域屏蔽效能等效测试方法将显得非常有意义。
发明内容
本发明的目的是提供一种操作简单方便、成本低的强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法，以便可以获取真实强电磁脉冲环境下屏蔽体的防护抑制效果，分析屏蔽体对强电磁脉冲的防护能力。
为实现以上目的，本发明通过以下技术方案实现：一种强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法，其特征是：至少包括如下步骤：
S1，扫频发射测试屏蔽体的频域屏蔽效能；
S2，将强电磁脉冲源时域波形通过傅立叶变换，获得频域谱特性；
S3，在脉冲源场强频谱特性基础上，使其相应频点按S1中获得的屏蔽效能进行相应衰减，得到经屏蔽体衰减后的场强频谱特性；
S4，将S3中获得的场强频谱特性，进行傅里叶逆变换获得对应脉冲源经屏蔽作用后的时域波形；通过对屏蔽前后时域波形的峰值之比，得到对应脉冲源作用下的时域屏蔽效能。
所述的步骤S1包含：
S1.1，建立等效测试系统；
S1.2，在没有放置屏蔽体时，计算机控制信号源发射相应频点的信号，频谱仪接收到的信号为P1(f)；
所述的相应频点信号的确定，指的是在9kHz～18GHz频率范围内，根据屏效在该范围内按照对数坐标连续缓慢变化的特点，将总的频率测试范围划分为不同的频段，在每个频段内测试频点按照对数分布；测试的频点总数在500～900之间动态调整；其中，9kHz～1MHz频段内测试频点数范围为40～80；1MHz～20MHz频段内测试频点数范围为40～80；20MHz～300MHz频段内测试频点数范围为150～250；300MHz～3GHz频段内测试频点数范围为150～250；3GHz～18GHz频段内测试频点数范围为100～200。
S1.3，放置屏蔽体时，计算机控制信号源发射与S1.2相对应的频点信号,并保持发射功率不变，频谱仪接收到的信号P2(f)；
S1.4，步骤S1.2和步骤S1.3中测得的信号按照公式
SE(f)＝P1(f)-P2(f)
计算即可得到屏蔽体的频域屏蔽效能SE(f)；其中P1(f)为无屏蔽体时，频谱仪接收到的信号，单位是dBm；P2(f)为放置屏蔽体后，频谱仪接收到的信号，单位是dBm；SE(f)为屏蔽体的屏蔽效能，单位是dB。
所述的步骤S2包含：
S2.1,若脉冲源波形是具有函数表达式的理论波形E1(t)，则通过输入脉冲源波形的参数，经过连续傅立叶变换得到脉冲源的频谱解析表达式E1(f)；
所述的连续傅立叶变换表达式为：
其中f表示频率变量，单位是Hz；t表示时间变量，单位是s；E1(t)表示脉冲源时域波形函数形形，单位是V/m；E1(f)表示脉冲源频域波形函数，单位是V/m/Hz；
若脉冲源信号E1(tn)，为测试获得的波形数据，则首先通过合理确定采样点，再经过离散快速傅里叶变换得到脉冲源的频谱E1(fk)；
所述的离散傅立叶变换为：
其中N表示脉冲源时域波形总点数；tn表示离散时间序列，单位是s；fs表示时域波形采样率，单位是Hz；n表示时域序列；k表示频域序列；表示脉冲源离散时域波形,单位是V/m；E1(fk)表示脉冲源离散频域波形，单位是V/m/Hz；
所述的合理确定采样点，指的是根据测试获得的波形数据得到该波形的上升时间tr，再依据带宽BW和上升时间tr的近似关系式，估算出信号带宽；带宽和上升时间的近似关系式为：
其中BW为信号带宽，单位是Hz；tr为时域信号上升时间，单位是s；
估算出信号的带宽，即对应信号最高频率fq1；根据奈奎斯特采样定律fs2≥2fq1即可，但为提高信号保真度和增加模数转换器分辨率，我们采用3～5倍的最高频率作为采样率fs2；
若实际测试得到的数据的采样率fs1＜fs2，则可以利用线性插值对该数据进行加密处理得E1(tn1)；为了确保结果的准确性，在条件允许的情况下，建议对该信号进行重新测量；
若实际测试得到的数据的采样率2fs2＞fs1＞fs2，则对该离散数据E1(tn)直接进行快速傅立叶变换；
若实际测试得到的数据的采样率fs1＞＞fs2，为减小数据长度及计算量，建议对该数据进一步采样得E1(tn2)。
所述的步骤S3包含：
S3.1，根据步骤S2中得到的源信号最高截止频率fq1和步骤S1中实际测试得到的频域屏蔽效能SE(f)，可以得到在该截止频率下的频域屏蔽效能SE1(f)；
S3.2，将脉冲源场强频谱E1(f)和实际测试得到的频域屏蔽效能SE1(f)通过屏效定义式计算，即可得到屏蔽后的频谱E2(f)；
所述的计算公式为：
其中E1(f)表示脉冲源场强频谱，单位是V/m/Hz；E2(f)表示脉冲源信号经屏蔽作用后的频谱，单位是V/m/Hz；SE1(f)表示实际测试得到的频域屏蔽效能，单位是dB。
所述的步骤S4包含：
S4.1，采用最小相位原理对频域数据E2(f)重构相位；
S4.2，对得到的强电磁脉冲经屏蔽后的信号频谱E2(f)进行线性插值处理，使其频点分布满足反傅立叶变换的条件；
所述的满足反傅里叶变换条件，指的是信号频域均匀分布的频点数Nf应大于等于该信号在时域中满足奈奎斯特采样定理的点数Nt；从经屏蔽体屏蔽后信号的频谱图E2(f)中，可以得到该信号的最大截止频率fq2,根据奈奎斯特采样定律及考虑信号保真度要求，选取采样率fs3＝(3～5)fq2；根据实际测试经验，考虑到屏蔽体对脉冲信号的时延作用，我们预估该信号的时域持续时间T2为脉冲信号持续时间T1的五倍值；利用采样率fs3和信号持续时间T2的关系，得到时域点数Nt；
根据测试经验，一般情况下有Nf＜Nt，则对该信号频域数据E2(f)进行线性插值处理，得到频域数据E2(fk1)；
S4.3，将经过插值后的频谱E2(fk1)进行反傅里叶变换，即可得到屏蔽后的时域波形E2(t)；
S4.4，将屏蔽后时域波形的最大幅度值和强电磁脉冲源时域波形的最大值通过时域屏效计算公式计算，即可得到时域屏蔽效能；
所述的时域屏效计算公式为：
其中E1(t)max表示脉冲源时域波形最大值，单位是V/m；E2(t)max表示脉冲源信号经屏蔽作用后时域波形最大值，单位是V/m；SE(t)表示时域屏蔽效能，单位是dB。
本发明的优点是：由于建立了等效测试系统，可以通过测试屏蔽体对强电磁脉冲相应频点的屏蔽效能，来获取真实强电磁脉冲作用下屏蔽体的时域屏蔽效能，进而分析屏蔽体对强电磁脉冲的防护能力。具有操作简单方便、成本低廉的特点。
附图说明
图1是屏蔽效能扫频测试系统；
图2是频域屏蔽效能转换到时域屏蔽效能框图；
图3是核电磁脉冲的辐射电场标准波形；
图4是核电磁脉冲辐射电场经屏蔽体屏蔽后的时域瞬态响应波形。
图中，1、屏蔽体；2、信号源；3、发射天线；4、接收天线；5、频谱仪；6、计算机；7、交换机；8、光纤发送接收端子；9、外光电转换器；10、内光电转换器。
具体实施方式
以下结合附图，通过详细说明一个具体实例，对本发明做进一步阐述。
实施例1
一种强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法，至少包括如下步骤：
S1，扫频发射测试屏蔽体的频域屏蔽效能；
S2，将强电磁脉冲源时域波形通过傅立叶变换，获得频域谱特性；
S3，在脉冲源场强频谱特性基础上，使其相应频点按S1中获得的屏蔽效能进行相应衰减，得到经屏蔽体衰减后的场强频谱特性；
S4，将S3中获得的场强频谱特性，进行傅里叶逆变换获得对应脉冲源经屏蔽作用后的时域波形；通过对屏蔽前后时域波形的峰值之比，得到对应脉冲源作用下的时域屏蔽效能。
所述的步骤S1包含：
S1.1，建立等效测试系统；
S1.2，在没有放置屏蔽体时，计算机控制信号源发射相应频点的信号，频谱仪接收到的信号为P1(f)；
所述的相应频点信号的确定，指的是在9kHz～18GHz频率范围内，根据屏效在该范围内按照对数坐标连续缓慢变化的特点，将总的频率测试范围划分为不同的频段，在每个频段内测试频点按照对数分布；测试的频点总数在500～900之间动态调整；其中，9kHz～1MHz频段内测试频点数范围为40～80；1MHz～20MHz频段内测试频点数范围为40～80；20MHz～300MHz频段内测试频点数范围为150～250；300MHz～3GHz频段内测试频点数范围为150～250；3GHz～18GHz频段内测试频点数范围为100～200。
S1.3，放置屏蔽体时，计算机控制信号源发射与S1.2相对应的频点信号,并保持发射功率不变，频谱仪接收到的信号P2(f)；
S1.4，步骤S1.2和步骤S1.3中测得的信号按照公式
SE(f)＝P1(f)-P2(f)
计算即可得到屏蔽体的频域屏蔽效能SE(f)；其中P1(f)为无屏蔽体时，频谱仪接收到的信号，单位是dBm；P2(f)为放置屏蔽体后，频谱仪接收到的信号，单位是dBm；SE(f)为屏蔽体的屏蔽效能，单位是dB。
所述的步骤S2包含：
S2.1,若脉冲源波形是具有函数表达式的理论波形E1(t)，则通过输入脉冲源波形的参数，经过连续傅立叶变换得到脉冲源的频谱解析表达式E1(f)；
所述的连续傅立叶变换表达式为：
其中f表示频率变量，单位是Hz；t表示时间变量，单位是s；E1(t)表示脉冲源时域波形函数形形，单位是V/m；E1(f)表示脉冲源频域波形函数，单位是V/m/Hz；
若脉冲源信号E1(tn)，为测试获得的波形数据，则首先通过合理确定采样点，再经过离散快速傅里叶变换得到脉冲源的频谱E1(fk)；
所述的离散傅立叶变换为：
其中N表示脉冲源时域波形总点数；tn表示离散时间序列，单位是s；fs表示时域波形采样率，单位是Hz；n表示时域序列；k表示频域序列；表示脉冲源离散时域波形,单位是V/m；E1(fk)表示脉冲源离散频域波形，单位是V/m/Hz；
所述的合理确定采样点，指的是根据测试获得的波形数据得到该波形的上升时间tr，再依据带宽BW和上升时间tr的近似关系式，估算出信号带宽。带宽和上升时间的近似关系式为：
其中BW为信号带宽，单位是Hz；tr为时域信号上升时间，单位是s；
估算出信号的带宽，即对应信号最高频率fq1；根据奈奎斯特采样定律fs2≥2fq1即可，但为提高信号保真度和增加模数转换器分辨率，我们采用3～5倍的最高频率作为采样率fs2；
若实际测试得到的数据的采样率fs1＜fs2，则可以利用线性插值对该数据进行加密处理得E1(tn1)；为了确保结果的准确性，在条件允许的情况下，建议对该信号进行重新测量；
若实际测试得到的数据的采样率2fs2＞fs1＞fs2，则对该离散数据E1(tn)直接进行快速傅立叶变换；
若实际测试得到的数据的采样率fs1＞＞fs2，为减小数据长度及计算量，建议对该数据进一步采样得E1(tn2)。
所述的步骤S3包含：
S3.1，根据步骤S2中得到的源信号最高截止频率fq1和步骤S1中实际测试得到的频域屏蔽效能SE(f)，可以得到在该截止频率下的频域屏蔽效能SE1(f)；
S3.2，将脉冲源场强频谱E1(f)和实际测试得到的频域屏蔽效能SE1(f)通过屏效定义式计算，即可得到屏蔽后的频谱E2(f)；
所述的计算公式为：
其中E1(f)表示脉冲源场强频谱，单位是V/m/Hz；E2(f)表示脉冲源信号经屏蔽作用后的频谱，单位是V/m/Hz；SE1(f)表示实际测试得到的频域屏蔽效能，单位是dB。
所述的步骤S4包含：
S4.1，采用最小相位原理对频域数据E2(f)重构相位；
S4.2，对得到的强电磁脉冲经屏蔽后的信号频谱E2(f)进行线性插值处理，使其频点分布满足反傅立叶变换的条件；
所述的满足反傅里叶变换条件，指的是信号频域均匀分布的频点数Nf应大于等于该信号在时域中满足奈奎斯特采样定理的点数Nt；从经屏蔽体屏蔽后信号的频谱图E2(f)中，可以得到该信号的最大截止频率fq2,根据奈奎斯特采样定律及考虑信号保真度要求，选取采样率fs3＝(3～5)fq2；根据实际测试经验，考虑到屏蔽体对脉冲信号的时延作用，我们预估该信号的时域持续时间T2为脉冲信号持续时间T1的五倍值；利用采样率fs3和信号持续时间T2的关系，得到时域点数Nt；
根据测试经验，一般情况下有Nf＜Nt，则对该信号频域数据E2(f)进行线性插值处理，得到频域数据E2(fk1)；
S4.3，将经过插值后的频谱E2(fk1)进行反傅里叶变换，即可得到屏蔽后的时域波形E2(t)；
S4.4，将屏蔽后时域波形的最大幅度值和强电磁脉冲源时域波形的最大值通过时域屏效计算公式计算，即可得到时域屏蔽效能；
所述的时域屏效计算公式为：
其中E1(t)max表示脉冲源时域波形最大值，单位是V/m；E2(t)max表示脉冲源信号经屏蔽作用后时域波形最大值，单位是V/m；SE(t)表示时域屏蔽效能，单位是dB。
实施例2
如图2所示，一种强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法，如在核电磁脉冲作用下方舱时域屏蔽效能的等效测试，包含：
S1，扫频发射测试频域屏蔽效能；S1包含：
S1.1，建立等效测试系统，如图1所示。在测试时，信号源2及发射天线3位于屏蔽体1的外侧，接收天线4和频谱仪5位于屏蔽体1的内侧。为了避免计算机6的辐射对测试结果的影响，将计算机6放在屏蔽体1的外侧。发射天线3距屏蔽体1的外侧垂直壁距离为d2,接收天线4距屏蔽体1的内壁垂直距离为d1。外侧的发射天线3与信号源2电连接，信号源2通过发射天线3向屏蔽体内发送信号。信号源2通过交换机7与计算机6电连接，交换机7同时通过光纤发送接收端子8与频谱仪5网络电连接，频谱仪5与屏蔽体内侧的接收天线4电连接；
S1.2，交换机7同时通过光纤发送接收端子8与频谱仪5网络电连接，光纤发送接收端子8的屏蔽体1的外侧有外光电转换器9，屏蔽体1的内侧有内光电转换器10，外光电转换器9光电连接交换机7和光纤发送接收端子8的外端；内光电转换器10光电连接频谱仪5和光纤发送接收端子8的内端；采用光电转换系统实现计算机对频谱仪的程控操作，解决线缆进出方舱导致方舱屏蔽效能降低的问题；
S1.3，在没有放置屏蔽体1时，计算机6控制信号源2发射相应频点的信号，频谱仪接收到信号P1(f)(dBm)。
所述的相应频点信号的确定，指的是在9kHz～18GHz频率范围内，根据屏效在该范围内按照对数坐标连续缓慢变化的特点，将总的频率测试范围划分为不同的频段，在每个频段内测试频点按照对数分布。测试的频点总数可在500～900之间动态调整，频点选取如表一所述；
表一测试频段及测试频点选取
S1.4，放置屏蔽体时，计算机6控制信号源2发射与S1.3相对应的频点信号，并保持发射功率不变，频谱仪5接收到的信号P2(f)(dBm)；
S1.5，步骤S1.3和步骤S1.4中测得的信号按照公式
SE(f)＝P1(f)-P2(f)
计算即可得到频域屏蔽效能SE(f)。
S2,将强电磁脉冲源时域波形通过傅立叶变换，获得频域谱特性；S2包含：
S2.1,核电磁脉冲时域场强波形选用国际电工技术委员会(IEC)标准波形E1(t)(如图3所示)，对其进行连续傅里叶变换，得到脉冲源频谱解析表达式E1(f)。
所述的国际电工技术委员会(IEC)标准波形为：
E1(t)＝E*k*(e-αt-e-βt)
其中，E＝5×104,k＝1.3，α＝4×107，β＝6×108。
所述的连续傅立叶变换为：
所述的脉冲源频谱解析表达式为：
S3，在脉冲源场强频谱特性基础上，使其相应频点按S1中获得的屏蔽效能进行相应衰减，得到经屏蔽体衰减后的场强频谱特性；S3包含：
S3.1，根据步骤S2中得到的源信号最高截止频率fq1为300MHz和步骤S1中实际测试得到的频域屏蔽效能SE(f)，可以得到在该截止频率下的频域屏蔽效能SE1(f)；
S3.2，将脉冲源场强频谱E1(f)和实际测试得到的频域屏蔽效能SE(f)通过屏效定义式，计算即可得到屏蔽后的频谱E2(f)。
所述的计算公式为：
S4，将S3中获得的场强频谱特性，进行傅里叶逆变换(IFFT)获得对应脉冲源经屏蔽作用后的时域波形。通过对屏蔽前后时域波形的峰值之比，得到对应脉冲源作用下的时域屏蔽效能；S4包含：
S4.1，采用最小相位原理对频域数据E2(f)重构相位；
S4.2，对得到的强电磁脉冲经屏蔽后的信号频谱E2(f)进行线性插值处理，使其频点分布满足反傅立叶变换的条件；
所述的满足反傅里叶变换条件，指的是信号频域均匀分布的频点数Nf应大于等于该信号在时域中满足奈奎斯特采样定理的点数Nt。从经屏蔽体屏蔽后信号的频谱图E2(f)中，可以得到该信号的最大截止频率fq2。根据奈奎斯特采样定律及考虑信号保真度要求，选取采样率fs3＝(3～5)fq2。根据实际测试经验，考虑到屏蔽体对脉冲信号的时延作用，我们预估该信号的时域持续时间T2为脉冲信号持续时间T1的五倍值。利用采样率fs3和信号持续时间T2的关系，得到时域点数Nt。
根据测试经验，一般情况下有Nf＜Nt，对该信号频域数据E2(f)进行线性插值处理，得到频域数据E2(fk1)；
S4.3，将经过插值后的频谱E2(fk1)(V/m/Hz)进行反傅里叶变换(IFFT)即可得到屏蔽后的时域波形E2(t)(V/m)(如图4所示)；
S4.4，将强电磁脉冲源时域波形的最大幅度值E1(t)max和屏蔽后时域波形的最大幅度值E2(t)max通过时域屏效计算公式计算，即可得到时域屏蔽效能。时域屏效计算公式如下：
本实施例没有详细叙述的部件和结构属本行业的公知部件和常用结构或常用手段，这里不一一叙述。
Claims (5)
1.一种强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法，其特征是：至少包括如下步骤：
S1，扫频发射测试屏蔽体的频域屏蔽效能；
S2，将强电磁脉冲源时域波形通过傅立叶变换，获得频域谱特性；
S3，在脉冲源场强频谱特性基础上，使其相应频点按S1中获得的屏蔽效能进行相应衰减，得到经屏蔽体衰减后的场强频谱特性；
S4，将S3中获得的场强频谱特性，进行傅里叶逆变换获得对应脉冲源经屏蔽作用后的时域波形；通过对屏蔽前后时域波形的峰值之比，得到对应脉冲源作用下的时域屏蔽效能。
2.根据权利要求1所述的一种强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法，其特征是：所述的步骤S1包含：
S1.1，建立等效测试系统；
S1.2，在没有放置屏蔽体时，计算机控制信号源发射相应频点的信号，频谱仪接收到的信号为P1(f)；
所述的相应频点信号的确定，指的是在9kHz～18GHz频率范围内，根据屏效在该范围内按照对数坐标连续缓慢变化的特点，将总的频率测试范围划分为不同的频段，在每个频段内测试频点按照对数分布；测试的频点总数在500～900之间动态调整；其中，9kHz～1MHz频段内测试频点数范围为40～80；1MHz～20MHz频段内测试频点数范围为40～80；20MHz～300MHz频段内测试频点数范围为150～250；300MHz～3GHz频段内测试频点数范围为150～250；3GHz～18GHz频段内测试频点数范围为100～200。
S1.3，放置屏蔽体时，计算机控制信号源发射与S1.2相对应的频点信号,并保持发射功率不变，频谱仪接收到的信号P2(f)；
S1.4，步骤S1.2和步骤S1.3中测得的信号按照公式
SE(f)＝P1(f)-P2(f)
计算即可得到屏蔽体的频域屏蔽效能SE(f)；其中P1(f)为无屏蔽体时，频谱仪接收到的信号，单位是dBm；P2(f)为放置屏蔽体后，频谱仪接收到的信号，单位是dBm；SE(f)为屏蔽体的屏蔽效能，单位是dB。
3.根据权利要求1所述的一种强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法，其特征是：所述的步骤S2包含：
S2.1,若脉冲源波形是具有函数表达式的理论波形E1(t)，则通过输入脉冲源波形的参数，经过连续傅立叶变换得到脉冲源的频谱解析表达式E1(f)；
所述的连续傅立叶变换表达式为：
<mrow>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mi>&infin;</mi>
</msubsup>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mi>t</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>t</mi>
</mrow>
其中f表示频率变量，单位是Hz；t表示时间变量，单位是s；E1(t)表示脉冲源时域波形函数形形，单位是V/m；E1(f)表示脉冲源频域波形函数，单位是V/m/Hz；
若脉冲源信号E1(tn)，为测试获得的波形数据，则首先通过合理确定采样点，再经过离散快速傅里叶变换得到脉冲源的频谱E1(fk)；
所述的离散傅立叶变换为：
<mrow>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>k</mi>
<mi>n</mi>
<mo>/</mo>
<mi>N</mi>
</mrow>
</msup>
<mo>&rsqb;</mo>
</mrow>
其中N表示脉冲源时域波形总点数；tn表示离散时间序列，单位是s；fs表示时域波形采样率，单位是Hz；n表示时域序列；k表示频域序列；表示脉冲源离散时域波形,单位是V/m；E1(fk)表示脉冲源离散频域波形，单位是V/m/Hz；
所述的合理确定采样点，指的是根据测试获得的波形数据得到该波形的上升时间tr，再依据带宽BW和上升时间tr的近似关系式，估算出信号带宽；带宽和上升时间的近似关系式为：
<mrow>
<mi>B</mi>
<mi>W</mi>
<mo>=</mo>
<mfrac>
<mn>0.35</mn>
<msub>
<mi>t</mi>
<mi>r</mi>
</msub>
</mfrac>
</mrow>
其中BW为信号带宽，单位是Hz；tr为时域信号上升时间，单位是s；
估算出信号的带宽，即对应信号最高频率fq1；根据奈奎斯特采样定律fs2≥2fq1即可，但为提高信号保真度和增加模数转换器分辨率，我们采用3～5倍的最高频率作为采样率fs2；
若实际测试得到的数据的采样率fs1＜fs2，则可以利用线性插值对该数据进行加密处理得E1(tn1)；为了确保结果的准确性，在条件允许的情况下，建议对该信号进行重新测量；
若实际测试得到的数据的采样率2fs2＞fs1＞fs2，则对该离散数据E1(tn)直接进行快速傅立叶变换；
若实际测试得到的数据的采样率fs1＞＞fs2，为减小数据长度及计算量，建议对该数据进一步采样得E1(tn2)。
4.根据权利要求1所述的一种强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法，其特征是：所述的步骤S3包含：
S3.1，根据步骤S2中得到的源信号最高截止频率fq1和步骤S1中实际测试得到的频域屏蔽效能SE(f)，可以得到在该截止频率下的频域屏蔽效能SE1(f)；
S3.2，将脉冲源场强频谱E1(f)和实际测试得到的频域屏蔽效能SE1(f)通过屏效定义式计算，即可得到屏蔽后的频谱E2(f)；
所述的计算公式为：
<mrow>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
</mrow>
<msup>
<mn>10</mn>
<mfrac>
<mrow>
<msub>
<mi>SE</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
</mrow>
<mn>20</mn>
</mfrac>
</msup>
</mfrac>
</mrow>
其中E1(f)表示脉冲源场强频谱，单位是V/m/Hz；E2(f)表示脉冲源信号经屏蔽作用后的频谱，单位是V/m/Hz；SE1(f)表示实际测试得到的频域屏蔽效能，单位是dB。
5.根据权利要求1所述的一种强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法，其特征是：所述的步骤S4包含：
S4.1，采用最小相位原理对频域数据E2(f)重构相位；
S4.2，对得到的强电磁脉冲经屏蔽后的信号频谱E2(f)进行线性插值处理，使其频点分布满足反傅立叶变换的条件；
所述的满足反傅里叶变换条件，指的是信号频域均匀分布的频点数Nf应大于等于该信号在时域中满足奈奎斯特采样定理的点数Nt；从经屏蔽体屏蔽后信号的频谱图E2(f)中，可以得到该信号的最大截止频率fq2,根据奈奎斯特采样定律及考虑信号保真度要求，选取采样率fs3＝(3～5)fq2；根据实际测试经验，考虑到屏蔽体对脉冲信号的时延作用，我们预估该信号的时域持续时间T2为脉冲信号持续时间T1的五倍值；利用采样率fs3和信号持续时间T2的关系，得到时域点数Nt；
根据测试经验，一般情况下有Nf＜Nt，则对该信号频域数据E2(f)进行线性插值处理，得到频域数据E2(fk1)；
S4.3，将经过插值后的频谱E2(fk1)进行反傅里叶变换，即可得到屏蔽后的时域波形E2(t)；
S4.4，将屏蔽后时域波形的最大幅度值和强电磁脉冲源时域波形的最大值通过时域屏效计算公式计算，即可得到时域屏蔽效能；
所述的时域屏效计算公式为：
<mrow>
<mi>S</mi>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>20</mn>
<mi>lg</mi>
<mfrac>
<mrow>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
其中E1(t)max表示脉冲源时域波形最大值，单位是V/m；E2(t)max表示脉冲源信号经屏蔽作用后时域波形最大值，单位是V/m；SE(t)表示时域屏蔽效能，单位是dB。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710639671.6A CN107329002B (zh) | 2017-07-31 | 2017-07-31 | 强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710639671.6A CN107329002B (zh) | 2017-07-31 | 2017-07-31 | 强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107329002A true CN107329002A (zh) | 2017-11-07 |
| CN107329002B CN107329002B (zh) | 2019-08-09 |
Family
ID=60199508
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710639671.6A Active CN107329002B (zh) | 2017-07-31 | 2017-07-31 | 强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107329002B (zh) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109239469A (zh) * | 2018-08-16 | 2019-01-18 | 江苏赛诺格兰医疗科技有限公司 | 检测磁屏蔽材料在低磁场下磁屏蔽效果的装置和方法 |
| CN111983341A (zh) * | 2020-07-15 | 2020-11-24 | 中国工程物理研究院应用电子学研究所 | 一种强电磁脉冲屏蔽效能测试用信号前处理模块及其系统 |
| CN112578192A (zh) * | 2020-12-30 | 2021-03-30 | 北京无线电计量测试研究所 | 一种小屏蔽体的时域脉冲屏蔽效能测试装置和方法 |
| CN113514708A (zh) * | 2021-04-23 | 2021-10-19 | 中国人民解放军军事科学院国防工程研究院工程防护研究所 | 一种大型地下工程整体hemp环境实测系统及测量方法 |
| CN114626268A (zh) * | 2022-03-18 | 2022-06-14 | 北京航空航天大学 | 一种高精度强电磁脉冲传播过程时域计算方法 |
| CN114722882A (zh) * | 2022-05-23 | 2022-07-08 | 北京智芯微电子科技有限公司 | 脉冲电场数据的处理方法与系统 |
| CN115420842A (zh) * | 2022-05-10 | 2022-12-02 | 鲲鹏仪器(大连)有限公司 | 一种高效液相色谱仪高压泵的精准压力控制方法 |
| CN117761393A (zh) * | 2024-02-22 | 2024-03-26 | 南京派格测控科技有限公司 | 一种时域信号的获取方法及装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101813732A (zh) * | 2010-05-10 | 2010-08-25 | 中国人民解放军理工大学 | 近区低频强磁场屏蔽效能测试系统及其测试方法 |
| CN103308798A (zh) * | 2013-05-27 | 2013-09-18 | 华北电力大学 | 一种电磁屏蔽材料屏蔽效能的测试方法 |
| CN106249104A (zh) * | 2016-07-29 | 2016-12-21 | 山东康威通信技术股份有限公司 | 一种通信电缆屏蔽层状态监测的检测装置及方法 |
-
2017
- 2017-07-31 CN CN201710639671.6A patent/CN107329002B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101813732A (zh) * | 2010-05-10 | 2010-08-25 | 中国人民解放军理工大学 | 近区低频强磁场屏蔽效能测试系统及其测试方法 |
| CN103308798A (zh) * | 2013-05-27 | 2013-09-18 | 华北电力大学 | 一种电磁屏蔽材料屏蔽效能的测试方法 |
| CN106249104A (zh) * | 2016-07-29 | 2016-12-21 | 山东康威通信技术股份有限公司 | 一种通信电缆屏蔽层状态监测的检测装置及方法 |
Non-Patent Citations (3)
| Title |
|---|
| 张龙等: ""强电磁脉冲对材料屏蔽效能的影响"", 《高电压技术》 * |
| 邱畅等: ""多芯电缆屏蔽效能的仿真计算与测试"", 《电线电缆》 * |
| 陈翔等: ""基于屏蔽暗箱窗口法的材料电磁脉冲屏蔽效能时域测试"", 《高电压技术》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109239469A (zh) * | 2018-08-16 | 2019-01-18 | 江苏赛诺格兰医疗科技有限公司 | 检测磁屏蔽材料在低磁场下磁屏蔽效果的装置和方法 |
| CN109239469B (zh) * | 2018-08-16 | 2021-11-02 | 江苏赛诺格兰医疗科技有限公司 | 检测磁屏蔽材料在低磁场下磁屏蔽效果的装置和方法 |
| CN111983341A (zh) * | 2020-07-15 | 2020-11-24 | 中国工程物理研究院应用电子学研究所 | 一种强电磁脉冲屏蔽效能测试用信号前处理模块及其系统 |
| CN111983341B (zh) * | 2020-07-15 | 2023-01-20 | 中国工程物理研究院应用电子学研究所 | 一种强电磁脉冲屏蔽效能测试用信号前处理模块及其系统 |
| CN112578192A (zh) * | 2020-12-30 | 2021-03-30 | 北京无线电计量测试研究所 | 一种小屏蔽体的时域脉冲屏蔽效能测试装置和方法 |
| CN113514708A (zh) * | 2021-04-23 | 2021-10-19 | 中国人民解放军军事科学院国防工程研究院工程防护研究所 | 一种大型地下工程整体hemp环境实测系统及测量方法 |
| CN114626268A (zh) * | 2022-03-18 | 2022-06-14 | 北京航空航天大学 | 一种高精度强电磁脉冲传播过程时域计算方法 |
| CN115420842A (zh) * | 2022-05-10 | 2022-12-02 | 鲲鹏仪器(大连)有限公司 | 一种高效液相色谱仪高压泵的精准压力控制方法 |
| CN115420842B (zh) * | 2022-05-10 | 2024-04-12 | 华谱科仪(北京)科技有限公司 | 一种高效液相色谱仪高压泵的精准压力控制方法 |
| CN114722882A (zh) * | 2022-05-23 | 2022-07-08 | 北京智芯微电子科技有限公司 | 脉冲电场数据的处理方法与系统 |
| CN117761393A (zh) * | 2024-02-22 | 2024-03-26 | 南京派格测控科技有限公司 | 一种时域信号的获取方法及装置 |
| CN117761393B (zh) * | 2024-02-22 | 2024-05-07 | 南京派格测控科技有限公司 | 一种时域信号的获取方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107329002B (zh) | 2019-08-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107329002B (zh) | 强电磁脉冲作用下屏蔽体时域屏蔽效能的等效测试方法 | |
| Nucci et al. | Lightning-induced voltages on overhead lines | |
| Krug et al. | The time-domain electromagnetic interference measurement system | |
| CN104459348B (zh) | 基于软件无线电的高功率微波辐射场测量装置及其方法 | |
| KR102203945B1 (ko) | 시설 내부로 침투하는 전자파 펄스의 전계강도 측정방법 | |
| CN108663575A (zh) | 大型屏蔽体低频脉冲强磁场屏蔽效能测试系统及测试方法 | |
| CN109116182B (zh) | 通信电缆接头屏蔽罩屏蔽效能测量的装置及方法 | |
| CN104391199A (zh) | 一种电缆屏蔽层屏蔽效能的测试方法 | |
| CN104849592B (zh) | 一种射电望远镜宽带电磁屏蔽效能检测系统及检测方法 | |
| CN207965113U (zh) | 低频脉冲强磁场环境下屏蔽体磁饱和性能测试装置 | |
| CN101533049A (zh) | 一种测量线缆的电流辐射发射装置 | |
| CN108089075A (zh) | 一种屏蔽线缆的仿真验证测试装置及其方法 | |
| Ceclan et al. | Lightning-inverse reconstruction by remote sensing and numerical-field synthesis | |
| CN112578192A (zh) | 一种小屏蔽体的时域脉冲屏蔽效能测试装置和方法 | |
| CN109839553A (zh) | 一种基于虚拟仪器的emc自动测试系统设计 | |
| CN207516465U (zh) | 一种屏蔽线缆的仿真验证测试装置 | |
| Carbonini | Comparison of analysis of a WTEM cell with standard TEM cells for generating EM fields | |
| CN211826244U (zh) | 一种冲击接地阻抗测试装置 | |
| CN106405341A (zh) | 一种还原局部放电超高频信号的方法 | |
| Putro et al. | Influence of spacer aperture size on PD-induced electromagnetic wave measured with UHF external sensor in 66 kV GIS model | |
| Schulze et al. | Improved transmission-line model for a cable with an attached suppression ferrite | |
| IDITA et al. | Indoor assessment of screening efficiency through absorption and reflection | |
| CN110765548A (zh) | 一种飞机线束防波套的屏蔽衰减测量方法及装置 | |
| CN106053945B (zh) | 一种短时s函数史密斯圆图分析装置及方法 | |
| Bulut et al. | Measurement of shielding effectiveness of different types of wire meshes in a large frequency range |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220926 Address after: Room 2-1-04, Southeast Chunxiao, Beishapo Village, Beilin District, Xi'an City, Shaanxi Province 710048 Patentee after: Xi'an Huanian Electronic Technology Co.,Ltd. Address before: 710071 Xi'an Electronic and Science University, 2 Taibai South Road, Shaanxi, Xi'an Patentee before: XIDIAN University |
|
| TR01 | Transfer of patent right |