CN106199182A - 基于随机电报信号噪声调节随机共振检测弱信号的方法 - Google Patents
基于随机电报信号噪声调节随机共振检测弱信号的方法 Download PDFInfo
- Publication number
- CN106199182A CN106199182A CN201610517109.1A CN201610517109A CN106199182A CN 106199182 A CN106199182 A CN 106199182A CN 201610517109 A CN201610517109 A CN 201610517109A CN 106199182 A CN106199182 A CN 106199182A
- Authority
- CN
- China
- Prior art keywords
- random telegraph
- telegraph signals
- signal
- random
- weak signal
- 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
Links
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
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
本发明提供了一种基于随机电报信号噪声调节随机共振检测弱信号的方法,针对具有弱信号的二维非线性系统,通过调节随机电报信号噪声的高低阈值、噪声关联时间和噪声强度,实现二维非线性系统与随机电报信号噪声和周期信号间最强的协作机制,即随机共振。此外,随机电报信号在电子信息工程领域中较为常见且易于操作处理,是一种非高斯型随机噪声,能够使得二维非线性系统在随机电报信号噪声的作用下产生明显的随机共振现象和较高的信噪比值,更容易对系统中的弱信号进行检测。本发明不仅是一种方便快捷的方法,且具有广泛的适用性,能够提高电子信息工程效能。
Description
技术领域
本发明属于电子信号处理领域,涉及一种根据随机共振进行信号检测的方法。
背景技术
利用随机激励下的随机共振原理,提取混杂在复杂系统中的需求信号,特别是微弱信号,是近些年电子信号处理中常用的一种信号检测和获取方法。这种方法能在适当的系统和随机激励条件下将噪声原本对系统有害的能力向需求信号部分转化,使得该部分的信号能力得到增强,从而将混杂在二维非线性系统中的需求信号特别微弱信号提取出来。这一方法在电子工程故障的检测、图像信息加密和识别等方面有着广泛的应用。
目前,二维非线性系统中随机共振的研究主要集中在高斯噪声、白噪声等随机激励下。而随机电报信号噪声能很好地刻画电子元器件中的部分功能的开关随机波动,如:存储器存储陷阱不可靠性、电子设备电荷迁移阈值电压扰动及晶体管器件载流子数涨落与迁移等。随机电报信号噪声在电子信息工程领域中是一种较为常见且易于操作处理的非高斯随机噪声,在电子信息工程领域中有着广泛的应用。
随机电报信号噪声激励下二维非线性系统的随机共振,通过系统功率谱、周期信号和随机噪声功率谱在原驱动频率处的比值即信噪比等的变化反映了二维非线性系统中丰富而复杂的动力学行为特点,并利用随机共振对弱信号进行检测。
发明内容
为了克服现有技术的不足,本发明提供一种基于随机电报信号噪声调节下,利用随机共振机制对二维非线性系统中的弱信号进行检测的方法,针对具有弱信号的二维非线性系统,通过调节随机电报信号噪声的高低阈值、噪声关联时间和噪声强度,实现二维非线性系统、随机电报信号噪声和周期信号之间最强的协作机制,即随机共振,并利用其进行弱信号的检测。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤1,生成随机电报信号噪声;
步骤2,求解二维非线性系统在随机电报信号噪声激励下系统响应的待处理时间序列,待处理时间序列的数据量与随机电报信号噪声序列的数据量相同;
步骤3,对待处理的时间序列进行等间隔采样,获取系统输出的时间序列;
步骤4,对二维非线性系统输出的时间序列进行频谱分析,获取二维非线性系统在随机电报信号噪声激励下的功率谱和信噪比;
步骤5,调节随机电报信号噪声的高低阈值、噪声强度和噪声关联时间,得到不同高低阈值、噪声强度和噪声关联时间下信噪比极大值的集合,并从该集合中选取最大值;获取该最大值所对应的随机电报信号噪声的高低阈值、噪声强度和噪声关联时间,即具有弱信号的二维非线性系统在随机电报信号噪声激励下呈现出最佳随机共振机制的最优随机电报信号噪声调节参数;
步骤6,将最优随机电报信号噪声调节参数下的随机电报信号噪声输入到具有弱信号的二维非线性系统中,对系统输出进行频谱分析,获取功率谱及其最大峰值所对应的弱信号的特征频率,最终完成弱信号的检验。
所述的步骤2采用四阶龙格库塔法求解待处理时间序列。
所述的步骤3采用100个数据点作为采样间隔。
所述的步骤4采用傅里叶变换公式对二维非线性系统输出的时间序列进行频谱分析。
本发明的有益效果是:
1.利用随机共振原理作为对弱信号进行检测的工具,即利用随机噪声、周期信号和二维非线性系统间的协同作用,将外在的随机噪声能量转化为供系统使用的能量,增强了复杂非线性系统中弱信号的展现,从而对其进行检测。
2.随机电报信号噪声是一种只在两个高低阈值间随机波动且等待时间服从指数分布的噪声,是一种非高斯型随机噪声,在电子信息工程领域中较为常见且易于操作处理。二维非线性系统在非高斯随机电报信号噪声的作用下,能够产生更高的信噪比值。
3.采用一种容易实施的非线性系统的随机共振原理对弱信号进行检测,不仅是一种方便快捷的方法,且具有广泛的适用性,能够提高电子信息工程效能。
附图说明
图1是本发明的方法流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
本发明的实现步骤如下:
步骤1,生成随机电报信号噪声。
步骤2,求解二维非线性系统,获取随机电报信号噪声激励下系统响应的待处理时间序列。
利用四阶龙格库塔算法,数值求解二维非线性系统在随机电报信号噪声激励下的系统响应的待处理时间序列,此处待处理时间序列数据量与随机电报信号噪声序列的数据量相同。
步骤3,对待处理的时间序列进行等间隔采样,获取系统输出的时间序列。
此处对待处理的时间序列进行等间隔采样,采样间隔取100个数据。获取系统输出的时间序列。
步骤4,对系统输出的时间序列进行频谱分析,获取功率谱和信噪比。
利用傅里叶变换公式对二维非线性系统输出的时间序列进行频谱分析,获取二维非线性系统在随机电报信号噪声激励下的功率谱,并计算系统的信噪比。
步骤5,调节随机电报信号噪声参数,获取各个参数下信噪比的极大值,并找出其中的最大值对应下的最优随机电报信号噪声调节参数。
在信噪比值具有极大值的情况下,调节随机电报信号的高低阈值、噪声强度和噪声关联时间等参数,确定不同参数调节下信噪比的极大值,并从该极大值的集合中选取最大值。反之获取该信噪比极大值中最大值所对应的随机电报信号参数,即具有弱信号的二维非线性系统在随机电报信号噪声激励下呈现出最佳随机共振机制的最优随机电报信号噪声调节参数。
步骤6,将最优随机电报信号噪声调节参数重新输入到二维非线性系统中,并对系统输出进行频谱分析获取功率谱,从中找出弱信号的特征频率。
根据最优调节参数选定随机电报信号噪声,并将最优调节参数下的随机电报信号噪声输入到具有弱信号的二维非线性系统中,通过该系统输出进行频谱分析并获取功率谱及其最大峰值所对应的弱信号的特征频率,最终完成弱信号的检验。
如图1所示,本发明的实施例包括以下步骤:
步骤1,生成随机电报信号噪声。
随机电报信号噪声数值时间序列的生成是实现弱信号检测的关键点,这里主要通过理论推导和数值算法快速生成。
设随机电报信号噪声为ζ(t),噪声状态只有高低两个阈值H,L,转移率分别为μH,μL,平均等待时间分别为tH,tL,随机电报信号噪声序列的初始时刻t0,初始状态为x0,任意时刻为t。P(H,t|x0,t0)表示初始时间和初始状态条件下任意时刻的状态为高阈值H的概率,P(L,t|x,t0)表示初始时间和初始状态条件下任意时刻的状态为低阈值L的概率。
随机电报信号噪声的稳态均值函数为
稳态关联函数为
<ζ(t)ζ(t′)>=τDexp(-τ|t-t′|).
其中,<ζ(t)ζ(t′)>为稳态关联函数,τ为随机电报信号噪声的关联时间,D为随机电报信号噪声强度,t和t′表示两个不同的任意时刻,exp(·)表示指数。
随机电报信号噪声强度D的表达式为
其中,λ为表示任意时刻的变量。
随机电报信号噪声从高阈值状态H或低阈值状态L跃迁到下一为高阈值状态H的概率PHH,PLH分别为
其中,tn为第n个时刻,tn+1为第n个时刻的下一时刻,即第n+1个时刻,dt为两个时刻间的时间间隔,选取0.001。
为了生成随机电报信号噪声时间序列,设定随机电报信号噪声为零均值噪声时,且噪声强度D、噪声关联时间τ和高低阈值H与L已知时,可以计算出随机电报信号μH和μL的值,从而计算PHH,PLH。通过如下方法生成随机电报信号噪声序列。
首先,设随机电报信号噪声ζ(t)的序列为{xi,i=0,1,2,…,N}初始状态假定为x0=H(或x0=L),同时,利用计算机生成一列在[0,1]区间上均匀分布的随机数{Ri,i=0,1,2,…,N},这里要求随机数个数N不小于106。接着,比较R0是否小于PHH(或R0是否小于PLH)。如果R0小于PHH(或R0小于PLH)取x1为H,如果R0大于等于PHH(或R0大于PLH)取x1为L。其次,比较R1是否小于PHH(或R1是否小于PLH)。如果R1小于PHH(或R1小于PLH)取x2为H,如果R1大于等于PHH(或R1大于PLH)取x2为L。依次做下去,最终将得到ζ(t)的N个数值的序列{xi,i=0,1,2,…,N}。
步骤2,求解二维非线性系统,获取随机电报信号噪声激励系统响应的待处理时间序列。
其中,X(t)为二维非线性系统变量,η为阻尼系数,V为二维非线性系统的势函数,Acos(ωt)为周期信号,A,ω分别为周期信号振幅和频率,ζ(t)为随机电报信号噪声。
利用四阶龙格库塔算法,数值求解该系统在随机电报信号噪声激励下系统响应的待处理时间序列,记为{Xi,i=1,2,…,N},此处待处理时间序列数据量N与随机电报信号噪声序列的数据量相同。
步骤3,对待处理的时间序列进行等间隔采样,获取系统输出的时间序列。
此处对待处理的时间序列进行等间隔采样,采样间隔取100个数据点。获取系统输出的时间序列{Yi,i=1,2,…,[N/100]}。其中,[·]为取整函数。
步骤4,对系统输出的时间序列进行频谱分析,获取功率谱和信噪比。
利用如下计算功率谱的变换公式对系统输出的时间序列进行频谱分析,获取二维非线性系统在随机电报信号噪声激励下的功率谱。
再选取如下公式计算信噪比:
其中,PS(ω)为周期信号处功率谱幅值即功率谱中的最大值,PN(ω)为周期信号处随机电报信号噪声的功率谱幅值,log10(·)为对数函数。
步骤5,调节随机电报信号噪声参数,获取各个参数下信噪比的极大值,并找出其中的最大值对应下的最优随机电报信号噪声调节参数。
在信噪比比值呈现极大值SNRMax的情况下,调节随机电报信号高低阈值H和L、噪声强度D和噪声关联时间τ,寻找信噪比极大值SNRMax的最大值。获取该信噪比极大值SNRMax中最大值所对应的最优随机电报信号噪声调节参数。
步骤6,将最优随机电报信号噪声调节参数重新输入到二维非线性系统中,并对系统输出进行频谱分析获取功率谱,从中找出弱信号的特征频率。
根据最优随机电报信号噪声调节参数选定随机电报信号噪声,并将最优调节参数下的随机电报信号噪声输入到二维非线性系统中,通过变换公式求得系统输出的功率谱,进行功率谱分析,寻找功率谱最大峰值处弱信号的特征频率,完成弱信号的检验。
Claims (4)
1.一种基于随机电报信号噪声调节随机共振检测弱信号的方法,其特征在于包括下述步骤:
步骤1,生成随机电报信号噪声;
步骤2,求解二维非线性系统在随机电报信号噪声激励下系统响应的待处理时间序列,待处理时间序列的数据量与随机电报信号噪声序列的数据量相同;
步骤3,对待处理的时间序列进行等间隔采样,获取系统输出的时间序列;
步骤4,对二维非线性系统输出的时间序列进行频谱分析,获取二维非线性系统在随机电报信号噪声激励下的功率谱和信噪比;
步骤5,调节随机电报信号噪声的高低阈值、噪声强度和噪声关联时间,得到不同高低阈值、噪声强度和噪声关联时间下信噪比极大值的集合,并从该集合中选取最大值;获取该最大值所对应的随机电报信号噪声的高低阈值、噪声强度和噪声关联时间,即具有弱信号的二维非线性系统在随机电报信号噪声激励下呈现出最佳随机共振机制的最优随机电报信号噪声调节参数;
步骤6,将最优随机电报信号噪声调节参数下的随机电报信号噪声输入到具有弱信号的二维非线性系统中,对系统输出进行频谱分析,获取功率谱及其最大峰值所对应的弱信号的特征频率,最终完成弱信号的检验。
2.根据权利要求1所述的基于随机电报信号噪声调节随机共振检测弱信号的方法,
其特征在于:所述的步骤2采用四阶龙格库塔法求解待处理时间序列。
3.根据权利要求1所述的基于随机电报信号噪声调节随机共振检测弱信号的方法,
其特征在于:所述的步骤3采用100个数据点作为采样间隔。
4.根据权利要求1所述的基于随机电报信号噪声调节随机共振检测弱信号的方法,
其特征在于:所述的步骤4采用傅里叶变换公式对二维非线性系统输出的时间序列进行频谱分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610517109.1A CN106199182A (zh) | 2016-07-04 | 2016-07-04 | 基于随机电报信号噪声调节随机共振检测弱信号的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610517109.1A CN106199182A (zh) | 2016-07-04 | 2016-07-04 | 基于随机电报信号噪声调节随机共振检测弱信号的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106199182A true CN106199182A (zh) | 2016-12-07 |
Family
ID=57465943
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610517109.1A Pending CN106199182A (zh) | 2016-07-04 | 2016-07-04 | 基于随机电报信号噪声调节随机共振检测弱信号的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106199182A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110045193A (zh) * | 2019-03-06 | 2019-07-23 | 珠海博雅科技有限公司 | 一种浮栅器件电报噪声测量系统及测量方法 |
CN110133383A (zh) * | 2019-04-15 | 2019-08-16 | 中国科学院上海光学精密机械研究所 | 基于噪声关联的噪声敏感系数测量方法 |
CN112153652A (zh) * | 2020-08-27 | 2020-12-29 | 浙江大学 | 一种基于微弱信号检测的LoRa和WiFi并发传输方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4527116A (en) * | 1981-06-30 | 1985-07-02 | Ecole Superieure D'electricite And Enertec | Process and device for system characterization by spectral analysis |
CN101191804A (zh) * | 2007-12-03 | 2008-06-04 | 中国人民解放军国防科学技术大学 | 自适应随机共振微弱信号检测方法 |
CN101561314A (zh) * | 2009-05-12 | 2009-10-21 | 中国人民解放军国防科学技术大学 | 随机共振-混沌微弱信号检测方法 |
CN101741349A (zh) * | 2009-12-07 | 2010-06-16 | 杭州电子科技大学 | 一种随机共振噪声滤除与微弱周期信号检出电路 |
CN101825665A (zh) * | 2010-05-28 | 2010-09-08 | 桂林电子科技大学 | 随机共振瞬变电磁弱信号检测方法 |
CN103475431A (zh) * | 2013-09-09 | 2013-12-25 | 南京邮电大学 | 一种低信噪比条件下的基于最佳随机共振的频谱感知方法 |
CN103699513A (zh) * | 2013-12-20 | 2014-04-02 | 中国科学技术大学 | 一种基于多尺度噪声调节的随机共振方法 |
CN104113707A (zh) * | 2013-04-22 | 2014-10-22 | 全视科技有限公司 | 用于随机数产生的设备、方法和系统及图像传感器 |
CN105701544A (zh) * | 2016-01-11 | 2016-06-22 | 南京信息工程大学 | 基于遗传算法和频率调制的随机共振微弱信号检测方法 |
-
2016
- 2016-07-04 CN CN201610517109.1A patent/CN106199182A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4527116A (en) * | 1981-06-30 | 1985-07-02 | Ecole Superieure D'electricite And Enertec | Process and device for system characterization by spectral analysis |
CN101191804A (zh) * | 2007-12-03 | 2008-06-04 | 中国人民解放军国防科学技术大学 | 自适应随机共振微弱信号检测方法 |
CN101561314A (zh) * | 2009-05-12 | 2009-10-21 | 中国人民解放军国防科学技术大学 | 随机共振-混沌微弱信号检测方法 |
CN101741349A (zh) * | 2009-12-07 | 2010-06-16 | 杭州电子科技大学 | 一种随机共振噪声滤除与微弱周期信号检出电路 |
CN101825665A (zh) * | 2010-05-28 | 2010-09-08 | 桂林电子科技大学 | 随机共振瞬变电磁弱信号检测方法 |
CN104113707A (zh) * | 2013-04-22 | 2014-10-22 | 全视科技有限公司 | 用于随机数产生的设备、方法和系统及图像传感器 |
CN103475431A (zh) * | 2013-09-09 | 2013-12-25 | 南京邮电大学 | 一种低信噪比条件下的基于最佳随机共振的频谱感知方法 |
CN103699513A (zh) * | 2013-12-20 | 2014-04-02 | 中国科学技术大学 | 一种基于多尺度噪声调节的随机共振方法 |
CN105701544A (zh) * | 2016-01-11 | 2016-06-22 | 南京信息工程大学 | 基于遗传算法和频率调制的随机共振微弱信号检测方法 |
Non-Patent Citations (4)
Title |
---|
付华等: "随机共振瓦斯微弱信号检测方法研究", 《激光技术》 * |
林激: "二值噪声下信号调制二阶线性系统随机共振", 《中国优秀硕士论文全文数据库 基础科学辑》 * |
樊养余等: "基于随机共振的任意大频率微弱信号检测方法研究", 《仪器仪表学报》 * |
许勇等: "二值噪声激励下非线性双稳系统的逻辑随机共振", 《中国电力大会-2013论文摘要集》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110045193A (zh) * | 2019-03-06 | 2019-07-23 | 珠海博雅科技有限公司 | 一种浮栅器件电报噪声测量系统及测量方法 |
CN110133383A (zh) * | 2019-04-15 | 2019-08-16 | 中国科学院上海光学精密机械研究所 | 基于噪声关联的噪声敏感系数测量方法 |
CN110133383B (zh) * | 2019-04-15 | 2022-04-05 | 中国科学院上海光学精密机械研究所 | 基于噪声关联的噪声敏感系数测量方法 |
CN112153652A (zh) * | 2020-08-27 | 2020-12-29 | 浙江大学 | 一种基于微弱信号检测的LoRa和WiFi并发传输方法 |
CN112153652B (zh) * | 2020-08-27 | 2022-05-24 | 浙江大学 | 一种基于微弱信号检测的LoRa和WiFi并发传输方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106199182A (zh) | 基于随机电报信号噪声调节随机共振检测弱信号的方法 | |
Reddy et al. | Fast tracking of power quality disturbance signals using an optimized unscented filter | |
CN103093244B (zh) | 一种基于稀疏分解的雷达信号识别方法 | |
CN105785338A (zh) | 一种频率捷变雷达载频的优化方法 | |
CN109359378A (zh) | 一种分段非线性双稳系统的弱信号检测方法 | |
CN112379859B (zh) | 二值采样的处理方法及装置、对抗样本的生成方法及装置 | |
Peng et al. | Implication of two-coupled differential van der pol duffing oscillator in weak signal detection | |
CN104850752A (zh) | 一种基于自适应随机共振的参数搜索范围确定方法 | |
Li et al. | Research on power quality disturbance identification and classification technology in high noise background | |
CN106197523A (zh) | 基于一阶线性系统的弱信号检测及复原 | |
CN101828911B (zh) | 一种基于奇异谱熵的神经元动作电位特征提取方法 | |
CN107688167B (zh) | 一种多时宽线性调频脉冲压缩信号幅度包络曲线生成方法 | |
CN103944593B (zh) | 基于独立接收单元的双稳态随机共振阵列处理方法及系统 | |
Gao et al. | Adaptive cascaded-bistable stochastic resonance system research and design | |
CN113625101B (zh) | 基于果蝇算法与随机共振的行波信号处理方法 | |
Wei et al. | Weak Signal Detection Based on WT-LSTM | |
CN104991230B (zh) | 基于信号自相关特征的雷达脉冲压缩增益均衡方法 | |
CN108362940A (zh) | 一种半周波畸变信号的动态谐波频率提取方法 | |
CN100464500C (zh) | 时-频混叠信号的几何特征滤波方法 | |
CN107340056A (zh) | 一种匹配随机共振检测方法 | |
Zhu et al. | Aiming to Complex Power Quality Disturbances: A Novel Decomposition and Detection Framework | |
CN108897018B (zh) | 基于混沌振子的弱glonass信号捕获方法 | |
Xiao-hong et al. | Synchronization and circuit experiment simulation of chaotic time-delay systems | |
Shunhu et al. | Travelling wave signal processing method based on fruit fly algorithm and stochastic resonance | |
Hai et al. | Short-term wind speed prediction based on the wavelet transformation and Adaboost neural network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161207 |
|
WD01 | Invention patent application deemed withdrawn after publication |