背景噪声分析
1.频谱的概念
一个复杂的振动信号,可以看成是由许多简谐分量叠加而成;这些简谐分量及其各自的振幅、频率和初相,这就叫做复杂振动的频谱。
指数衰减周期震荡信号合成微地震子波,公式为:
A
(
t
)
=
A
0
e
−
h
(
t
−
t
0
)
s
i
n
[
w
(
t
−
t
0
)
A(t)=A_0e^{-h(t-t_0)}{sin}[w(t-t_0)
A(t)=A0e−h(t−t0)sin[w(t−t0)]
,
t
>
t
0
,t>t_0
,t>t0
其中
A
0
A_0
A0为信号的初始振幅,
h
h
h为阻尼因子,
w
w
w为角频率,
t
0
t_0
t0为信号的初始相位。
通过对
h
h
h和
w
w
w等参数选取合适的参数,可以得到如下所示的为地震子波。
我们取单道的时域图
我们取多道的时域图
我们取100道数据的时域图
![](https://i-blog.csdnimg.cn/blog_migrate/4239a245bbc8598f4e2383f664425627.png)
- 混沌性
在杜芬方程中微小的摄动回忆起系统的周期解从本质上发生巨大的变化。Duffing振子方程是非线性振动系统,具有非线性恢复力,在外部激励的作用下,由于外加周期摄动力的频率与系统本身的频率相互作用,系统将发生振荡,表现出极其丰富的动力学行为。
相态图的变化
Duffing方程
x
¨
+
k
x
˙
+
α
x
+
β
x
3
=
γ
c
o
s
(
w
t
)
\ddot{x}+k\dot{x}+\alpha{x}+\beta{x}^3 = \gamma cos(wt)
x¨+kx˙+αx+βx3=γcos(wt)
其中
k
k
k为阻尼比,
γ
\gamma
γ
w
w
w分别代表周期策动力的幅值和频率,
α
\alpha
α
β
\beta
β分别为 -1 , 1。
故
x
¨
+
k
x
˙
−
x
+
x
3
=
γ
c
o
s
(
w
t
)
\ddot{x}+k\dot{x}-{x}+{x}^3 = \gamma cos(wt)
x¨+kx˙−x+x3=γcos(wt)
我们取k = 0.5 随着
γ
\gamma
γ的变化,系统状态也会发生变化如下图所示
(1)当
γ
\gamma
γ = 0时
① 当 γ \gamma γ = 0时,初值为(1,1),这个时候时域图与相轨图如下所示
![](https://i-blog.csdnimg.cn/blog_migrate/675ed5f645e43adef6e6ed0b555ab218.png)
![](https://i-blog.csdnimg.cn/blog_migrate/b6287527c651db49839a2164572b134d.png)
②当 γ \gamma γ = 0时,初值为(-1,-1),这个时候时域图与相轨图如下
![](https://i-blog.csdnimg.cn/blog_migrate/d33faeb45912b9df9dda0d900108e064.png)
![](https://i-blog.csdnimg.cn/blog_migrate/6fdc705fc0a1991295804147f39bab4e.png)
(2) γ > 0 时 \gamma>0时 γ>0时
①当 γ > 0 \gamma>0 γ>0时, γ = 0.2 \gamma=0.2 γ=0.2, w = 1 w=1 w=1,初值为(0, 0)
![](https://i-blog.csdnimg.cn/blog_migrate/43413407e5d2357e11adef63c44d3042.png)
![](https://i-blog.csdnimg.cn/blog_migrate/223466bfbaad5ae010eb5ee9fce1fe8d.png)
②当 γ > 0 \gamma>0 γ>0时, γ = 0.3289 \gamma=0.3289 γ=0.3289, w = 1 w=1 w=1,初值为(0, 0)
![](https://i-blog.csdnimg.cn/blog_migrate/0c46f4a2fee89e4b685c663dde7107f4.png)
![](https://i-blog.csdnimg.cn/blog_migrate/dde400ee9fdc64e7a7510e8eaf7512b1.png)
③当 γ > 0 \gamma>0 γ>0时, γ = 0.733 \gamma=0.733 γ=0.733, w = 1 w=1 w=1,初值为(0, 0)
![](https://i-blog.csdnimg.cn/blog_migrate/3fc33a02e2a718b9e6dcf8ab2db05829.png)
![](https://i-blog.csdnimg.cn/blog_migrate/74ba31e5aefacb1007d161f1e29a0b6c.png)
④当 γ > 0 \gamma>0 γ>0时, γ = 0.734 \gamma=0.734 γ=0.734, w = 1 w=1 w=1,初值为(0, 0)
![](https://i-blog.csdnimg.cn/blog_migrate/13119f1af8e257ff54061fe4ed4af0a3.png)
![](https://i-blog.csdnimg.cn/blog_migrate/e342c4055ac8cdd1c7cb164bd527e955.png)
根据上几个图片我们不难发现当 γ \gamma γ取到一定的值的时候我们发现相轨图突然由混沌态转为周期态,而 γ \gamma γ的取值变化是非常细微的。本周学会了对混沌性的duffing方程中的图像进行绘制,下周继续深度的学习混沌性。