CN111487464A - 一种多载波微放电效应频谱特征解析方法 - Google Patents
一种多载波微放电效应频谱特征解析方法 Download PDFInfo
- Publication number
- CN111487464A CN111487464A CN202010354889.9A CN202010354889A CN111487464A CN 111487464 A CN111487464 A CN 111487464A CN 202010354889 A CN202010354889 A CN 202010354889A CN 111487464 A CN111487464 A CN 111487464A
- Authority
- CN
- China
- Prior art keywords
- electron
- carrier
- thin layer
- micro
- discharge
- 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
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/12—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing
- G01R31/1227—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Testing Relating To Insulation (AREA)
Abstract
本发明公开了一种多载波微放电噪声效应的频谱特征解析方法,首先通过施加多载波激励电场信号到两平行导体面形成的缝隙两侧,模拟多载波微放电效应;其次,将微放电效应稳态过程所形成的谐振电子云视为电子密度很高的电子薄层,利用其来表征谐振电子云所有的物理特性,确定其运动速度和形成的电流大小;然后对运动速度进行傅里叶级数展开,确定多载波微放电效应的频谱特征解析式,进而得到多载波微放电产生的电子噪声频率成分。这样得到了激励源频率、谐振腔(缝隙)宽度、激励电压等关键影响因子对微放电效应的作用情况,实现二者数值关系的定量描述,同时推导出多载波微放电效应所产生的电子噪声具体频谱成分,使得从本质上抑制其电子噪声干扰成为可能。
Description
技术领域
本发明属于电子噪声测量技术领域,特别是多载波微放电效应产生的电子噪声抑制技术方面,更为具体地讲,涉及一种多载波微放电噪声效应的频谱特征解析方法。
背景技术
微放电效应指的是由于二次电子倍增效应所发生的一种电子雪崩现象,研究已证实它是射频噪声、微波器件失谐、无源互调干扰产生的关键原因。特别是新形势下高功率、多载波、多制式微波通信系统的广泛应用,多载波微放电效应由于对微波器件及通信系统潜在的巨大威胁,一直是国内外同行关注的热点和研究的重点,例如:西班牙瓦伦西亚理工大学的S.Anza(2014)在他的博士论文中系统给出了多载波微放电效应的理论框架和预测模型,对多载波微放电效应的SEY方程、多载波阈值特征进行了深入研究;浙江大学宋庆庆博士在他的博士论文《多载波微放电过程的概率分析》中,采用随机漫步和BranchingLevy漫步模型,对多载波微放电过程中二次电子横向扩散所需遵循的概率模型进行了推导。
如上所述,现有技术绝大多数都集中于多载波微放电效应的瞬态发生过程,即如何有效地通过各种手段抑制多载波微放电效应的产生。而实际情况是:多载波微放电效应很难得到完全彻底的抑制,那么,如果对多载波微放电的稳态过程进行充分研究,明确其产生的电子噪声频率成分,在通信系统的频带选择时进行有效规避,也能很好地抑制由于多载波微放电效应所导致的射频干扰。
发明内容
本发明的目的在于克服现有技术的不足,提出一种多载波微放电噪声效应的频谱特征解析方法,以明确激励源频率、谐振腔(狭窄缝隙)宽度、激励电压等关键影响因子对微放电效应的作用情况,实现二者数值关系的定量描述,同时推导出多载波微放电效应所产生的电子噪声具体频谱成分,使得从本质上抑制其电子噪声干扰成为可能。
为实现上述发明目的,本发明多载波微放电噪声效应的频谱特征解析方法,其特征在于,包括以下步骤:
(1)、利用射频信号发生器产生多载波激励电场信号并通过耦合器施加到两平行导体面形成的缝隙中(实际中也可为微波器件形成的狭小缝隙),其中,两平行导体面的间距为d、表面积为S、多载波激励电场信号幅值均为E0、载波频率为fl、n为载波数量、l为载波序号、l=1,2,…n,缝隙内会发生微放电效应,达到稳定时,将充满谐振的电子云;
(2)、视谐振电子云全部集中到缝隙正中间,形成一电子薄层(电子密度很高),然后利用电子薄层表征谐振电子云所有的物理特性;
根据电子薄层所受到的电场力、磁场力以及所满足的牛顿-洛伦兹关系,确定其运动速度v(t)为:
其中,m,e为单个电子的质量和电量,M为多载波微放电的阶数,取正奇数(即:M=1,3,5…),t0、V0分别为电子薄层的发射时间和初始运动速度,T为电子薄层的运动周期;
(3)、周期性谐振运动的电子薄层所形成的电流(电子薄层电流)可以表示为:
其中,z表示位置变量,t表示时间变量,n′为电子薄层的总电子数,n′e为电子薄层的总电荷,ze(t)为t时刻的电子薄层位置,δ为狄拉克函数,n′e·δ(z-ze(t))等价为电子薄层电荷密度,通过在微放电效应达到稳定状态时插入正极性探针测出,S为电子薄层面积(即平行导体面的表面积),为电子薄层的运动方向;
其中,k0=2i-1,i为正整数,f0为多载波激励信号的中心频率,且满足f0∈gcd(f1,f2,...fn);
多载波微放电产生的电子噪声频率成分为:k0f0/M。
本发明的目的是这样实现的。
本发明多载波微放电噪声效应的频谱特征解析方法,首先通过施加多载波激励电场信号到两平行导体面形成的缝隙两侧,模拟多载波微放电效应,其次,将微放电效应稳态过程所形成的谐振电子云视为电子密度很高的电子薄层,利用该电子薄层来表征谐振电子云所有的物理特性,确定其运动速度和形成的电流大小,然后对运动速度进行傅里叶级数展开,确定多载波微放电效应的频谱特征解析式,进而得到多载波微放电产生的电子噪声频率成分。这样得到了激励源频率、谐振腔(缝隙)宽度、激励电压等关键影响因子对微放电效应的作用情况,实现二者数值关系的定量描述,同时推导出多载波微放电效应所产生的电子噪声具体频谱成分,使得从本质上抑制其电子噪声干扰成为可能。
附图说明
图1是本发明多载波微放电噪声效应的频谱特征解析方法一种具体实施方式流程图;
图2是本发明中多载波微放电效应稳态分析模型;
图3是本发明中多载波微放电效应电子噪声时域仿真图;
图4是本发明中多载波微放电效应电子噪声频域仿真图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
图1是本发明多载波微放电噪声效应的频谱特征解析方法一种具体实施方式流程图。
在本实施例中,如图1所示,本发明多载波微放电噪声效应的频谱特征解析方法,包括以下步骤:
步骤S1:施加多载波激励电场信号到缝隙,发生微放电效应,使其充满谐振的电子云。
利用射频信号发生器产生多载波激励电场信号并通过耦合器施加到两平行导体面形成的缝隙中(实际中也可为微波器件形成的狭小缝隙),其中,两平行导体面的间距为d、表面积为S、多载波激励电场幅值均为E0、载波频率为fl、n为载波数、l为载波序号、l=1,2,…n,缝隙内会发生微放电效应,达到稳定时,将充满谐振的电子云。该谐振电子云的周期性运动是产生电子噪声的本质原因。
步骤S2:用电子薄层表征谐振电子云所有的物理特征,根据电子行为与力学关系确定电子薄膜层运动速度。
在本发明中,视谐振电子云全部集中到缝隙正中间,形成一电子薄层(电子密度很高),然后利用电子薄层表征谐振电子云所有的物理特性。
本发明首先根据谐振电子云的运动情况,建立一种平行板导体面的多载波微放电效应稳态分析模型,如图2所示。其中,图2(a)为多载波微放电效应达到稳态时,两平行板导体面内充满周期运动的谐振电子云。可将此谐振电子云的全部电子等效至缝隙正中间处(d/2),简化为一电子密度极高、厚度极薄的电子薄层,如图2(b)所示。那么,能够利用此电子薄层完全表征出谐振电子云的所有物理特性。
其中,m为电子质量,t0为电子薄层的起始时间,对公式(2)进行积分结果为:
式中,V0为电子薄层的初始运动速度,当电子薄层谐振时,会在两平行导体面之间往复运动,撞击平行导体面前和平行导体面后的电子薄膜层初始运动速度会反向,那么在一个运动周期内,可将式(3)改写,即确定电子薄层运动速度v(t)为:
其中,m,e为单个电子的质量和电量,M为多载波微放电的阶数,为了满足微放电效应电子的渡越时间为电场半周期的奇数倍这个条件,M取正奇数(即:M=1,3,5…);t0、V0分别为电子薄层的发射时间和初始运动速度,T为电子薄层的运动周期。初始运动速度V0根据平行导体所用的材料以及加载的多载波激励电场幅值E0确定,具体可以查询平行导体所用材料的电场幅值-初始运动速度关系表确定。
由式(4)可知,电子薄层运动速度可以改写为两部分组成,即:v(t)=vosc(t)+vcon(t),vosc(t)为式(4)中黑色实线部分,是速度v(t)的连续分量,它不会产生新的频率分量;vcon(t)为式(4)中黑色虚线部分,是以MT为周期跳变的常数分量,它会产生新的频率成分。因此,影响多载波微放电效应电子噪声的关键在于其速度v(t)的常数分量部分vcon(t)。因此,对分量vcon(t)进行傅里叶级数展开为:
步骤S3:根据牛顿-洛伦兹关系,确定电子薄层周期性谐振运动形成的电流。
如图2所示,上板为A,下板为B,电子薄层的面积与平行板导体面积均为S,运动速度为v(t),且运动方向为z轴的正方向,那么周期性谐振运动的电子薄层所形成的电流(电子薄层电流)可以表示为:
其中,z表示位置变量,t表示时间变量,n′为电子薄层的总电子数,e为单个电子电量,n′e为电子薄层的总电荷,ze(t)为t时刻的电子位置,δ为狄拉克函数。由于狄拉克函数可以表示一个物理量在空间或时间中分布的密度,所以,n′e·δ(z-ze(t))可以等价为面电荷密度,可以通过在微放电效应达到稳定状态时插入正极性探针测出,S为电子薄层面积(即平行导体面的表面积),为电子薄层的运动方向。
步骤S4:联立步骤S3所述的电子薄层电流和步骤S2确定的电子薄层运动速度v(t),同时,对电子薄层的运动速度v(t)进行傅里叶级数展开,即将公式(5)代入公式(6)中,可以确定多载波微放电效应的频谱特征解析式为:
其中,k0=2i-1,i为正整数,f0为多载波激励信号的中心频率,且满足f0∈gcd(f1,f2,...fn)。
分析所述的多载波微放电效应频谱特征表达式可以发现:多载波微放电效应产生的电子噪声频率成分为:k0f0/M,它与多载波微放电的阶数M以及中心频率f0有关。例如:在发生一阶微放电时(M=1),其产生的谐波频率为:k0f0,当发生高阶微放电效时(M≥3),其产生的次谐波频率为:k0f0/M。
仿真
现以最简单多载波微放电情形n=2为例进一步说明,同时,为便于计算忽略单位量级,取:
那么,当发生一阶双载波微放电时(M=1),其频谱特征表达式可简化为:
由式(8)易知:当发生双载波一阶微放电时(M=1),将出现f0,3f0,5f0,7f0…等新的中心频率f0奇数倍的干扰分量。利用PIC-FDTD程序,且设置:f1=1.5GHz,f2=2.0GHz,其仿真结果如图3、4所示。仿真结果与分析式(8)有较好的一致性,由分析式(8)和仿真结果图4可以知道:干扰频率越高,其幅值越小,故衰减越大,所以主要干扰为接近载波频率f1和f2的几个较低阶谐波与次谐波干扰分量,最严重的为一阶、三阶和五阶干扰分量(即:f0,3f0,5f0,其中,3f0与载波频率f1相同,使得频率为f1的载波幅值略高于频率为f2的载波幅值,如图4所示),已被较多实验检测证实。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (1)
1.一种多载波微放电噪声效应的频谱特征解析方法,其特征在于,包括以下步骤:
(1)、利用射频信号发生器产生多载波激励电场信号并通过耦合器施加到两平行导体面形成的缝隙中(实际中也可为微波器件形成的狭小缝隙),其中,两平行导体面的间距为d、表面积为S、多载波激励电场幅值均为E0、载波频率为fl、n为载波数量、l为载波序号、l=1,2,…n,缝隙内会发生微放电效应,达到稳定时,将充满谐振的电子云;
(2)、视谐振电子云全部集中到缝隙正中间,形成一电子薄层(电子密度很高),然后利用电子薄层表征谐振电子云所有的物理特性;
根据电子薄层所受到的电场力、磁场力以及所满足的牛顿-洛伦兹关系,确定其运动速度v(t)为:
其中,m,e为单个电子的质量和电量,M为多载波微放电的阶数,取正奇数(即:M=1,3,5…),t0、V0分别为电子薄层的发射时间和初始运动速度,T为电子薄层的运动周期;
(3)、周期性谐振运动的电子薄层所形成的电流(电子薄层电流)可以表示为:
其中,z表示位置变量,t表示时间变量,n′为电子薄层的总电子数,n′e为电子薄层的总电荷,ze(t)为t时刻的电子薄层位置,δ为狄拉克函数,n′e·δ(z-ze(t))等价为电子薄层电荷密度,通过在微放电效应达到稳定状态时插入正极性探针测出,S为电子薄层面积(即平行导体面的表面积),为电子薄层的运动方向;
其中,k0=2i-1,i为正整数,f0为多载波激励信号的中心频率,且满足f0∈gcd(f1,f2,...fn);
多载波微放电产生的电子噪声频率成分为:k0f0/M。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010354889.9A CN111487464B (zh) | 2020-04-29 | 2020-04-29 | 一种多载波微放电效应频谱特征解析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010354889.9A CN111487464B (zh) | 2020-04-29 | 2020-04-29 | 一种多载波微放电效应频谱特征解析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111487464A true CN111487464A (zh) | 2020-08-04 |
CN111487464B CN111487464B (zh) | 2022-02-01 |
Family
ID=71798708
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010354889.9A Active CN111487464B (zh) | 2020-04-29 | 2020-04-29 | 一种多载波微放电效应频谱特征解析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111487464B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060176045A1 (en) * | 2000-05-12 | 2006-08-10 | Fujitsu Limited | Plasma reactor and method of determining abnormality in plasma reactor |
JP2008075120A (ja) * | 2006-09-20 | 2008-04-03 | Ulvac Japan Ltd | 成膜装置 |
CN101799494A (zh) * | 2010-01-29 | 2010-08-11 | 南京邮电大学 | 基于盲信号处理技术的传导电磁干扰噪声软分离方法 |
CN102801680A (zh) * | 2012-07-20 | 2012-11-28 | 西安空间无线电技术研究所 | 一种确定多载波微放电二次电子数目的方法 |
CN102984109A (zh) * | 2012-10-31 | 2013-03-20 | 西安空间无线电技术研究所 | 一种确定空间微波部件给定相位多载波微放电阈值的方法 |
CN104062565A (zh) * | 2014-06-24 | 2014-09-24 | 西安空间无线电技术研究所 | 一种利用互调分量检测微波部件微放电的方法 |
CN105072076A (zh) * | 2015-08-21 | 2015-11-18 | 西安空间无线电技术研究所 | 一种确定微波部件多载波微放电最坏状态的方法 |
CN109962323A (zh) * | 2017-12-26 | 2019-07-02 | 中国移动通信集团设计院有限公司 | 单激励微放电射频噪声信号抑制方法及装置 |
CN110089038A (zh) * | 2016-12-22 | 2019-08-02 | 松下知识产权经营株式会社 | 噪声抑制装置、噪声抑制方法、以及使用它们的接收装置、接收方法 |
-
2020
- 2020-04-29 CN CN202010354889.9A patent/CN111487464B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060176045A1 (en) * | 2000-05-12 | 2006-08-10 | Fujitsu Limited | Plasma reactor and method of determining abnormality in plasma reactor |
JP2008075120A (ja) * | 2006-09-20 | 2008-04-03 | Ulvac Japan Ltd | 成膜装置 |
CN101799494A (zh) * | 2010-01-29 | 2010-08-11 | 南京邮电大学 | 基于盲信号处理技术的传导电磁干扰噪声软分离方法 |
CN102801680A (zh) * | 2012-07-20 | 2012-11-28 | 西安空间无线电技术研究所 | 一种确定多载波微放电二次电子数目的方法 |
CN102984109A (zh) * | 2012-10-31 | 2013-03-20 | 西安空间无线电技术研究所 | 一种确定空间微波部件给定相位多载波微放电阈值的方法 |
CN104062565A (zh) * | 2014-06-24 | 2014-09-24 | 西安空间无线电技术研究所 | 一种利用互调分量检测微波部件微放电的方法 |
CN105072076A (zh) * | 2015-08-21 | 2015-11-18 | 西安空间无线电技术研究所 | 一种确定微波部件多载波微放电最坏状态的方法 |
CN110089038A (zh) * | 2016-12-22 | 2019-08-02 | 松下知识产权经营株式会社 | 噪声抑制装置、噪声抑制方法、以及使用它们的接收装置、接收方法 |
CN109962323A (zh) * | 2017-12-26 | 2019-07-02 | 中国移动通信集团设计院有限公司 | 单激励微放电射频噪声信号抑制方法及装置 |
Non-Patent Citations (3)
Title |
---|
HUIPING YANG ET AL.: "An Analytical Method to Evaluate the Spectrum of Multicarrier Multipactor Discharge", 《IEEE TRANSACTIONS ON ELECTRON DEVICES》 * |
杨会平: "微放电效应与松散连接对无源互调的影响机理研究", 《中国博士学位论文全文数据库(电子期刊)》 * |
林舒等: "粒子模拟在空间大功率微波器件微放电效应研究中的应用", 《真空电子技术》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111487464B (zh) | 2022-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Anza et al. | Nonstationary statistical theory for multipactor | |
Cooney et al. | Experiments on Korteweg–de Vries solitons in a positive ion–negative ion plasma | |
Trotta et al. | Electron acceleration at quasi-perpendicular shocks in sub-and supercritical regimes: 2D and 3D simulations | |
Thejappa et al. | In situ detection of strong Langmuir turbulence processes in solar type III radio bursts | |
Anza et al. | Long-term multipactor discharge in multicarrier systems | |
Wei et al. | Interferometric and schlieren characterization of the plasmas and shock wave dynamics during laser-triggered discharge in atmospheric air | |
Zeng et al. | Hybrid numerical simulation of decomposition of SF 6 under negative DC partial discharge process | |
Lakhina et al. | Generation of weak double layers and low-frequency electrostatic waves in the solar wind | |
Abbasi et al. | Preventing the recurrence effect in the Vlasov simulation by randomizing phase-point velocities in phase space | |
Zhang et al. | A numerical model of acoustic wave caused by a single positive corona source | |
Xiao et al. | Fundamental theory of townsend discharge | |
Esfahani et al. | Locust: C++ software for simulation of RF detection | |
CN111487464B (zh) | 一种多载波微放电效应频谱特征解析方法 | |
Rasch et al. | Non-resonant multipactor—A statistical model | |
Burke et al. | Experimental study of fluctuations excited by a narrow temperature filament in a magnetized plasma | |
Ayubi et al. | Surface discharge characteristics and influencing factors of polyimide insulations under high‐frequency sinusoidal voltages | |
Hanzelka et al. | Theories of growth and propagation of parallel whistler-mode chorus emissions: A review | |
Zhang et al. | Quasinormal modes of massless scalar field perturbation of a Reissner–Nordström de Sitter black hole with a global monopole | |
Raghukanth et al. | Empirical mode decomposition of earthquake accelerograms | |
Daldorff et al. | Nonlinear beam generated plasma waves as a source for enhanced plasma and ion acoustic lines | |
Thomas et al. | A treecode to simulate dust–plasma interactions | |
Arnaut | Pulse jitter, delay spread, and Doppler shift in mode-stirred Reverberation | |
Mohammadnejad et al. | Formation of ion acoustic rogue waves in warm dense matter | |
Chevalier et al. | Fluid simulation of the collisionless plasma sheath surrounding an electric dipole antenna in the inner magnetosphere | |
Su et al. | A frequency estimation algorithm based on cross information fusion |
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 |