西安交通大学 非线性电路实验报告

《西安交通大学 非线性电路实验报告》由会员分享，可在线阅读，更多相关《西安交通大学 非线性电路实验报告（9页珍藏版）》请在装配图网上搜索。

1、Duffing 方程及其在信号检测中的应用李禹锋(西安交通大学电力设备电气绝缘国家重点实验室，陕西西安710049)摘要：在工程领域中，在噪声环境下对信号进行检测一直都是研究的重点课题。混沌理 论表明一类混沌系统在一定条件下对小信号具有参数敏感性，同时对噪声具有免疫力，因此 使得它在信号检测中非常具有发展潜力。为此，本文分析了Duffing方程的动力学特性，研 究了利用Duffing方程来进行微弱信号检测的原理和过程，并在Matlab平台下进行了仿真实 验。结果表明，可以利用Duffing方程在噪声背景下进行信号的检测。关键词：混沌理论；信号检测； Duffing 方程；仿真研究1 引言在噪声

2、背景中检测微弱的有用信号是工程应用中的一个重要内容，前人已经 开展了大量的研究工作。传统的基于线性理论的信号检测方法由于对噪声背景下 的输出信噪比难以提高而存在一定局限性，尤其在对强噪声背景下的微弱信号检 测更是受到了限制。然而很多研究证明，利用“混沌振子对周期小信号具有敏感 依赖性，而对噪声具有免疫性”的特点，从噪声背景中提取微弱的周期信号是一 种行之有效的方法，引起了人们极大的兴趣1。在众多的信号检测中，正弦或余弦信号的检测占有极其重要的地位，在许多 领域中有着极其广泛的应用。本文采用余弦小信号作为检测对象，在 Matlab 平 台下，对Duffing方程及其在信号检测中的应用进行了初步探

3、讨。 2 基于 Duffing 方程的信号检测2.1 Duffing 方程的数学模型及分析Duffing 方程已被证明是混沌系统，大量学者对其进行过许多研究，研究它 的动力学行为可以揭示系统的各种性质。Duffing系统所描述的非线性动力学系 统表现出丰富的非线性动力学特性，目前已成为研究混沌现象的常用模型2。d2xdt 2霍尔姆斯型 Duffing 方程为：+ kdx 一 x(t) + x3 (t) = y cos(t) (1) dt式中，y cos(t)为周期策动力；k为阻尼比；-x(t)+x3(t)为非线性恢复力。其状态方程为：吟=ydtdt=一ky + x 一 x3 + y cos(t

4、) (3)在k固定的情况下，系统状态随Y的变化出现变化，具体分析如下：(1)当策动力Y为0时，计算得到相平面中结点为(0,0)和鞍点为(1,0)。系统在无策动力驱动作用下，周期性地在两个鞍点之一周围运动，围绕哪个鞍点运动则依赖于初始条件。由如图1可见，(x,y)初值分别设置是(1,1)和(2,2)时，(x,y)最 终分别停在(1,0)和(-1,0)不同的鞍点上。图1 y =0时的相轨迹(2)当策动力Y不为0时，当Y处于一定范围内时，系统都处于混沌状态， 如图2所示。随着Y的增大，当大于某一阀值Y d时，系统进入大尺度的周期运动，如图 3 所示。对应的相轨迹将中点、鞍点统统围住，此时系统的 Po

5、incare 映射为不动点。-1.51111111-2-1.5-1-0.500.511.525150050-1图 2 混沌状态图 3 大尺度的周期运动2.2 正弦信号的检测原理 对于一个非线性动力学系统 , 其参数的摄动有时会引起周期解发生本质的变化。因此，将待测信号作为Duffing方程周期策动力的摄动，噪声虽然强烈，但对系统状态的改变无影响, 而一旦带有特定的信号, 即使幅值较小, 也会使系统发生相变。计算机通过辨识系统状态, 可清楚地检测出特定信号是否存在。假设待检信号形式为s = hcos(t)+n(t), h为有用信号的幅值，n(t)为均值为0 的噪声。在对待测信号进行检测之前，调节

6、策动力信号幅度值卩，使尸，这时 系统处于混沌临界状态。然后输入待测信号s,经过暂态过程以后，系统稳定在 某一运动形式上。通过用眼观察得到系统的状态，由此可判断输入是纯粹噪声还 是稍带有周期信号。若输入信号中含有待测周期信号 ,则系统由混沌状态迅速转 入周期状态；若输入的待测信号中只有噪声，则系统状态不改变。因此，由输入 待测信号前后系统状态是否变化可以判断出输入待检测信号中是否含有待测周 期信号。当我们己经检测出输入信号中含有待测周期信号，需要进一步检测该信 号的幅值时，采用以下方法：逐渐减小策动力信号幅度值直到系统恢复到混沌临界状态。我们设这时 的策动力幅值为片，则待测信号中正弦信号幅值h为

7、h=yd-片，该式中Yd为该系 统处于临界状态的策动力幅值，即混沌阀值。3 MATLAB仿真实现与分析3.1 仿真模型的建立d2xdt2前面已经介绍过，霍尔姆斯型Duffing方程的一般形式为：+ kdx 一 x(t) + x3 (t) = y cos(t) (4) dt当少=1时，加入被检测的强噪声背景下微弱周期小信号s，Duffing方程可 写为如下形式：dxdty(5)dt=-ky + x 一 x3 + (y + h)cos(t) + n(t) (6)令阻尼比k=0.5，当e=1时，系统的Matlab仿真框图如下:在仿真模型中，Sine Wave用来产生系统的周期策动力信号，SineWa

8、vel和Random Number用来产生含有周期信号与噪声的待测信号，Fcn为系统的非 线性项， Gain 为放大器， Integrator 和 Integrator1 为积分器， XYGraph 为系统的 相轨迹显示器，Duffing方程中的参数y和h，即为Sine Wave和Sine Wavel中 产生的信号的幅值。根据文献4,这里取y d=0.8245。3.2 仿真实验结果分析3.2.1 只加入纯噪声信号首先，在y =0.5的情况下，在系统中只加入白噪声信号(方差为0.1),得到 加入噪声前后系统的相图。1.510.50-0.5-1-0.51.5(a)未加入噪声10.50-0.5-1-

9、2 -1.5 -1 -0.5 0 0.5 1 1.5 2(b)加入噪声图 5y =0.5 时白噪声信号对系统相图的影响从图5 可以看出，噪声虽然进入到系统中，但系统仍然保持混沌的状态，由 此证明了 Duffing 振子对于噪声具有免疫性。在一定范围内，即使噪声再强， 也不能使系统从混沌状态跃迁至周期状态。其差异仅仅在于：有噪声作用后，轨 道变得比较粗糙。3.2.2 加入待检信号与噪声混合后的信号首先，在尸0.8245的情况下，在系统中加入含有白噪声信号(方差为0.1) 的周期信号0.003cos(t)，得到加入信号后的系统的相图，如图6所示。图 6 加入含有白噪声信号的周期信号时系统的相图由图

10、 6可知，当加入含有白噪声信号的周期信号时，系统进入大尺度的周期 运动，说明此时信号中含有周期信号。为了检测信号中余弦信号幅值的大小，将?逐渐减小，当系统的相图由大尺 度的周期运动变成混沌临界状态时，由2.2可知，此时信号幅值h为h=yd- Y。 本文通过实验测得，h=0.0026。误差产生的原因是难以确定精确的Yd。4 结论通过对Duffing方程及检测原理的介绍.以及Matlab仿真实验和分析，表明 基于Duffing方程设计的余弦信号检测系统能够实现在噪声背景下检测出微弱的 周期小信号。参考文献袁野.基于混沌理论的微弱信号检测与勘探地震学中同相轴的恢复D.吉林大 学， 2006 年.2 刘崇新.非线性电路理论及应用.西安交通大学出版社， 2007年.徐艳春，杨春玲.基于混沌振子的微弱信号检测技术研究J.电气应用，2008, 27(8)：38-41.4聂春燕.基于混沌相平面变化的弱信号检测方法研究J.长春大学学报，1999，9(4)：1-4.