CN109740204B - 一种基于杜芬振子逆向相变的弱正弦信号快速检测方法 - Google Patents
一种基于杜芬振子逆向相变的弱正弦信号快速检测方法 Download PDFInfo
- Publication number
- CN109740204B CN109740204B CN201811553023.XA CN201811553023A CN109740204B CN 109740204 B CN109740204 B CN 109740204B CN 201811553023 A CN201811553023 A CN 201811553023A CN 109740204 B CN109740204 B CN 109740204B
- Authority
- CN
- China
- Prior art keywords
- phase change
- duffing
- vibrators
- detection
- phase
- 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
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Gyroscopes (AREA)
Abstract
本发明涉及一种基于杜芬振子逆向相变的弱正弦信号快速检测方法,解决的是原有方法计算结果包含因检测阵列振子数有限而引入的量化误差,无法准确限定相变振子数,单次测量误差明显,收敛过程慢,计算量大的技术问题,通过将检测模型中s个振子发生相变时的参考信号理想幅值用QINC(s)、QDEC(s)代替计算,较好解决了该问题。所述QINC(s)、QDEC(s)为搜索逼近得到的参考信号幅值,逼近程度可通过搜索精度的大小调整,s为限定参数;QINC(s)加上搜索精度的量,则检测阵列相变振子数由大于s变为s;QDEC(s)减去搜索精度的量,则检测阵列相变振子数由小于s变为s。
Description
技术领域
本发明涉及信号检测领域,具体涉及一种基于杜芬振子逆向相变的弱正弦信号快速检测方法。
背景技术
杜芬振子是一非线性振动方程,其相轨迹可以随着策动力幅度的增大依次出现同宿轨道、分岔、混沌运动、大周期运动等运动状态,是用于研究混沌的常用模型之一。杜芬振子混沌态和大周期态之间的状态转变具有对周期信号极其敏感,而对噪声有一定免疫能力的特性。杜芬振子弱信号检测方法正是利用此特性进行弱信号检测。
现有的基于杜芬振子逆向相变的弱信号检测方法由于检测阵列振子数有限,相变振子数并不能连续地表示策动力幅度的变化,因此计算结果包含量化误差,无法准确限定相变振子数,单次测量误差明显,收敛过程慢,计算量大。
本发明提供了一种克服上述技术问题的基于杜芬振子逆向相变的弱正弦信号快速检测方法。
发明内容
本发明所要解决的技术问题是现有技术中存在的计算结果包含因相变振子数不能连续表示策动力幅度变化而引起的量化误差,无法准确限定相变振子数,单次测量误差明显,收敛过程慢,计算量大的技术问题。提供一种新的基于杜芬振子逆向相变的弱正弦信号快速检测方法,该基于杜芬振子逆向相变的弱正弦信号快速检测方法具有加速了检测结果的收敛过程,减少了计算量和数据量需求,增强了检测结果的可靠性和检测过程的可控性的特点。
为解决上述技术问题,采用的技术方案如下:
一种基于杜芬振子逆向相变的弱正弦信号快速检测方法,所述弱正弦信号快速检测方法包括:
步骤1,建立弱正弦信号幅度检测的阵列模型;
步骤2,设定检测阵列模型中振子初始稳定大周期态,并输入待测信号;
步骤3,通过快速搜索逼近得到QINC(s)、QDEC(s);所述QINC(s)为参考信号幅值,QINC(s)加上搜索精度的量,则检测阵列相变振子数由大于s变为s;所述QDEC(s)为参考信号幅值,QDEC(s)减去搜索精度的量,则检测阵列相变振子数由小于s变为s;
步骤4,限定N2=N2',将QINC(s)、QDEC(s)带入弱正弦信号幅度检测的计算模型,计算出待测信号的估值为:
其中,N1为检测阵列每组振子总数;N2为第一组发生相变的振子数,N2'为第二组未发生相变的振子数。
本发明的工作原理:现有技术由于检测阵列振子数有限,相变振子数并不能连续地表示策动力幅度的变化,因此计算结果包含“量化误差”,方案中只是要求第一组相变振子数N2和第二组未相变振子数N2'统计上相等,无法准确限定N2、N2'大小,而计算过程是通过多次测量取平均来消除该“量化误差”,因此计算过程对每段数据利用率较低,单次测量误差明显,收敛过程慢,要求更多的数据,计算量大。而且计算过程是考虑到γ1、γ2>>εA这一情况而得到的近似表达式,因此增加了系统噪声,进一步减慢了收敛过程。本发明将原计算模型中角平分线上对应的参考信号幅值进行了转换,转变为可以通过搜索逼近方法计算得到的量,并且给出了一种快速搜索逼近的方法,可以对相变振子数、搜素精度进行严格限制,因此测量结果收敛速度快,需要的数据量变小,检测过程的可控性增强。
上述方案中,为优化,进一步地,步骤3包括:
步骤A1,定义理想状态下,系统策动力临界阈值γc位于相邻的相变振子和未相变振子相位差的平分线上,得出:
其中,
步骤A3,根据快速搜索逼近方法求解出QINC(s)、QDEC(s)。
进一步地,所述步骤4还包括,进行L次测量,将L次测量的平均值作为待测信号的估值:
进一步地,步骤1包括:构建检测阵列,所述检测阵列包括两组1×N1个杜芬振子,第1组第n个杜芬振子和第2组的第n个杜芬振子初始相轨迹点、参考信号初始相位均相同,调整两组1×N1检测振子的参考信号频率等于待测信号频率,使得1×N1个杜芬振子参考信号初始相位等间隔分布于[0,2π],所述间隔为2π/N1,所述检测阵列中第(m,n)杜芬振子的检测方程为:
进一步地,步骤2包括将利用数值仿真方法确定稳定大周期态下步骤1所述杜芬振子相轨迹点(x,y)和系统策动力相位Φ的对应关系,根据所述对应关系设置杜芬振子的初始相轨迹点位置(x0n,y0n)和参考信号初始相位Φn;其中,(x0n,y0n)和Φn分别为第n个振子的初始相轨迹点和初始相位。
QINC(s)、QDEC(s)是表示通过快速搜索使检测阵列满足特定条件而得到的参考信号幅度值,具体通过快速搜索计算过程和限定条件求得。其中,参考信号幅值初始搜索步长InStep一般选择比较大,搜索精度Precision用于控制参考信号幅度搜索最终步长的上限。每次搜索方向LABEL变化,搜索步长都会减小一半,因此能够快速搜索,很快达到搜索精度的要求,比普通的定步长搜索方法要快得多,搜索精度可控。
本发明的有益效果:本发明更加精确、快速地搜索调整参考信号幅值,使检测阵列迅速按照限定的振子数量和要求的精度发生相变,并更加充分利用每段数据进行计算,从而加速了检测结果的收敛过程,减少了计算量和数据量需求,增强了检测结果的可靠性和检测过程的可控性,有利于促进其工程应用。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1,基于杜芬振子逆向相变的弱正弦信号快速检测方法示意图。
图2,QDEC(s)计算流程示意图。
图3,QINC(s)计算流程示意图。
图4,理想状态下系统策动力等效矢量合成及其局部放大示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
实施例1
本实施例提供一种基于杜芬振子逆向相变的弱正弦信号快速检测方法,如图1,所述弱正弦信号快速检测方法包括:
步骤1,建立弱正弦信号幅度检测的阵列模型;
步骤2,设定阵列模型中振子初始稳定大周期态,并输入待测信号r(t);
步骤3,通过快速搜索逼近,获得QINC(s)、QDEC(s);所述QINC(s)为参考信号幅值,QINC(s)加上搜索精度Precision的量,则检测阵列相变振子数由大于s变为s;所述QDEC(s)为参考信号幅值,QDEC(s)减去搜索精度Precision的量,则检测阵列相变振子数由小于s变为s;
步骤4,限定N2=N2',将QINC(s)、QDEC(s)带入弱正弦信号幅度检测的计算模型,计算出待测信号的估值为:
其中,N1为检测阵列每组振子总数;N2为第一组振子发生相变的振子数。
如图4,步骤3包括:
步骤A1,定义理想状态下,系统策动力临界阈值γc位于相邻的相变振子和未相变振子相位差的平分线上,得出:
其中,
步骤A3,根据快速搜索逼近方法求解出QINC(s)、QDEC(s)。
QINC(s)、QDEC(s)表示通过快速搜索使检测阵列满足特定条件而得到的参考信号幅度值,具体快速搜索计算过程和限定条件如图2和图3所示。
具体地,所述步骤4还包括,进行L次测量,将L次测量的平均值作为待测信号的估值:
本实施例中,N1为每组阵列振子数,是根据检测精度要求确定的,越大越好。I可根据实际测量情况调整,2~5为理想取值范围。
具体地,步骤1包括:构建检测阵列,所述检测阵列包括两组1×N1个杜芬振子,第1组第n个杜芬振子和第2组的第n个杜芬振子初始相轨迹点、参考信号初始相位均相同,调整两组1×N1检测振子的参考信号频率等于待测信号频率,使得1×N1个杜芬振子参考信号初始相位等间隔分布于[0,2π],所述间隔为2π/N1,所述检测阵列中第(m,n)杜芬振子的检测方程为:
具体地,步骤2包括将利用数值仿真方法确定稳定大周期态下步骤1所述杜芬振子相轨迹点(x,y)和系统策动力相位Φ的对应关系,根据所述对应关系设置杜芬振子的初始相轨迹点位置(x0n,y0n)和参考信号初始相位Φn;其中,(x0n,y0n)和Φn分别为第n个振子的初始相轨迹点和初始相位。
QINC(s)、QDEC(s)是表示通过快速搜索使检测阵列满足特定条件而得到的参考信号幅度值,具体通过快速搜索计算过程和限定条件求得。其中,参考信号幅值初始搜索步长InStep一般选择比较大,搜索精度Precision用于控制参考信号幅度搜索最终步长的上限。计算过程中,每次搜索方向LABEL变化,搜索步长都缩减为1/2,很快即可收敛到满足限制条件,比普通的定步长搜索方法要快得多,搜索精度可控。因此相同收敛效果情况下比当前方法需要的数据量少,收敛快。
本实施例的其余部分是与“一种基于杜芬振子逆向相变的弱信号检测方法”相同,一种基于杜芬振子逆向相变的弱信号检测方法是直接利用了角平分线上对应的参考信号幅值作为计算用,这个量无法直接测量得到,但是可通过多次独立测量取平均得到其近似,因此原方案只是要求多次测量的相变振子数统计上相同即可,约束较松。本实施例中,将角平分线上对应的参考信号幅值进行了转换,转变为可以通过搜索逼近方法计算得到的量,并且给出了一种快速搜索逼近的方法,可以对相变振子数、搜素精度进行严格限制,因此测量结果收敛速度快,需要的数据量变小,检测过程的可控性增强。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员能够理解本发明,但是本发明不仅限于具体实施方式的范围,对本技术领域的普通技术人员而言,只要各种变化只要在所附的权利要求限定和确定的本发明精神和范围内,一切利用本发明构思的发明创造均在保护之列。
Claims (4)
1.一种基于杜芬振子逆向相变的弱正弦信号快速检测方法,其特征在于:所述弱正弦信号快速检测方法包括:
步骤1,建立弱正弦信号幅度检测的阵列模型;
步骤2,设定检测阵列振子初始稳定大周期态,并输入待测信号r(t);
步骤3,通过快速搜索逼近,得到QINC(s)、QDEC(s);所述QINC(s)为参考信号幅值,QINC(s)加上搜索精度的量,则检测阵列相变振子数由大于s变为s;所述QDEC(s)为参考信号幅值,QDEC(s)减去搜索精度的量,则检测阵列相变振子数由小于s变为s;
步骤4,将QINC(s)、QDEC(s)带入检测的计算模型,计算出待测信号的估值为:
其中,N1为检测阵列每组振子总数;N2为第一组发生相变的振子数;
步骤3包括:
步骤A1,定义理想状态下,系统策动力临界阈值γc位于相邻的相变振子和未相变振子相位差的平分线上,得出:
其中,
步骤A3,根据快速搜索逼近方法求解出QINC(s)、QDEC(s),A为幅度值;N2'为第二组未发生相变的振子数。
4.根据权利要求3所述的基于杜芬振子逆向相变的弱正弦信号快速检测方法,步骤2包括将利用数值仿真方法确定稳定大周期态下步骤1所述杜芬振子相轨迹点(x,y)和系统策动力相位Φ的对应关系,根据所述对应关系设置杜芬振子的初始相轨迹点位置(x0n,y0n)和参考信号初始相位Φn;其中,(x0n,y0n)和Φn分别为第n个振子的初始相轨迹点和初始相位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811553023.XA CN109740204B (zh) | 2018-12-19 | 2018-12-19 | 一种基于杜芬振子逆向相变的弱正弦信号快速检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811553023.XA CN109740204B (zh) | 2018-12-19 | 2018-12-19 | 一种基于杜芬振子逆向相变的弱正弦信号快速检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109740204A CN109740204A (zh) | 2019-05-10 |
CN109740204B true CN109740204B (zh) | 2021-01-26 |
Family
ID=66360461
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811553023.XA Active CN109740204B (zh) | 2018-12-19 | 2018-12-19 | 一种基于杜芬振子逆向相变的弱正弦信号快速检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109740204B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111999556B (zh) * | 2020-07-24 | 2023-06-16 | 南京信息职业技术学院 | 一种基于Duffing振子系统的线谱检测方法、系统及其存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103913222A (zh) * | 2014-04-25 | 2014-07-09 | 重庆邮电大学 | 基于锁相环的Duffing振子弱信号时域检测方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4655871B2 (ja) * | 2005-10-19 | 2011-03-23 | 株式会社日立製作所 | 永久磁石同期電動機の弱め界磁ベクトル制御装置及びモジュール |
CN100571098C (zh) * | 2006-02-24 | 2009-12-16 | 中国科学院上海微系统与信息技术研究所 | 通信系统中低复杂度的极大似然检测方法及装置 |
CN101294845B (zh) * | 2008-05-05 | 2011-11-02 | 西北工业大学 | 转子早期故障的多频率微弱信号检测方法 |
RU2485311C2 (ru) * | 2008-11-22 | 2013-06-20 | Карамай Дриллинг Текнолоджи Рисёрч Институт, Вестерн Дриллинг Инджиниринг Корпорейшн Оф Цнпц | Способ и система для передачи данных о стволе скважины |
CN101982945B (zh) * | 2010-11-12 | 2013-07-24 | 北京邮电大学 | 一种基于混沌理论的频谱感知方法 |
CN103323538B (zh) * | 2013-05-23 | 2016-03-09 | 暨南大学 | 基于杜芬方程Lyapunov指数的超声导波检测方法 |
CN104965123A (zh) * | 2015-07-07 | 2015-10-07 | 哈尔滨电工仪表研究所 | 一种基于混沌振子的电力系统谐波、间谐波检测与估计新方法 |
CN107066740B (zh) * | 2017-04-18 | 2020-04-14 | 中国人民解放军重庆通信学院 | 一种基于杜芬振子逆向相变的弱信号检测方法 |
-
2018
- 2018-12-19 CN CN201811553023.XA patent/CN109740204B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103913222A (zh) * | 2014-04-25 | 2014-07-09 | 重庆邮电大学 | 基于锁相环的Duffing振子弱信号时域检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109740204A (zh) | 2019-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111143989B (zh) | 频率调整量计算方法、模块、系统、存储介质和设备 | |
CN105091970B (zh) | 超声波流量计动态补偿方法 | |
CN107271956A (zh) | 非视距环境中未知起始时间的基于到达时间的定位方法 | |
CN104698845B (zh) | 基于自抗扰控制器的雷达天线伺服跟踪方法与系统 | |
US11199429B2 (en) | Method for measuring a speed of a fluid | |
CN107066740B (zh) | 一种基于杜芬振子逆向相变的弱信号检测方法 | |
CN105928546A (zh) | 编码器的信号处理装置 | |
CN109740204B (zh) | 一种基于杜芬振子逆向相变的弱正弦信号快速检测方法 | |
CN110297250A (zh) | 基于泰勒展开的初始掠射角求解方法、声线弯曲修正方法和设备 | |
CN105245203A (zh) | 高精度低速时钟占空比检测系统及方法 | |
CN103870238A (zh) | 基于fpga的伪随机序列发生器及其生成方法 | |
CN104280047A (zh) | 一种多传感器融合的陀螺漂移滤波系统及方法 | |
US3820022A (en) | Asymmetrical wave digital phase measuring system | |
CN114442078A (zh) | 检测飞行时间的方法、超声波流量计及光学设备 | |
CN105093280A (zh) | 表层模型对地震数据影响的低频与高频成分的分解方法 | |
US3440876A (en) | Flowmeter frequency control system | |
CN111538079A (zh) | 基于全波形反演技术确定地质裂缝柔度参数的方法及装置 | |
CN108732912A (zh) | 被测信号边沿触发的时钟分相法 | |
US20110249782A1 (en) | Clock recovery | |
Luong-Van et al. | Covariance profiling for an adaptive Kalman filter to suppress sensor quantization effects | |
CN109164419A (zh) | 基于交互式多模型的多平台无序量测处理方法 | |
Fortunato et al. | On the representation of bathymetry by unstructured grids | |
US10866079B2 (en) | Position sensing device | |
Vihman et al. | Adaptive Kalman Filter Based Data Aggregation in Fault-Resilient Underwater Sensor Networks | |
CN106052562A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |