2021-8-2 混沌性

背景噪声分析

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)=A0eh(tt0)sin[w(tt0)] , 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道数据的时域图

  1. 混沌性

在杜芬方程中微小的摄动回忆起系统的周期解从本质上发生巨大的变化。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),这个时候时域图与相轨图如下所示

②当 γ \gamma γ = 0时,初值为(-1,-1),这个时候时域图与相轨图如下

(2) γ > 0 时 \gamma>0时 γ>0

①当 γ > 0 \gamma>0 γ>0时, γ = 0.2 \gamma=0.2 γ=0.2, w = 1 w=1 w=1,初值为(0, 0)

②当 γ > 0 \gamma>0 γ>0时, γ = 0.3289 \gamma=0.3289 γ=0.3289, w = 1 w=1 w=1,初值为(0, 0)

③当 γ > 0 \gamma>0 γ>0时, γ = 0.733 \gamma=0.733 γ=0.733, w = 1 w=1 w=1,初值为(0, 0)

④当 γ > 0 \gamma>0 γ>0时, γ = 0.734 \gamma=0.734 γ=0.734, w = 1 w=1 w=1,初值为(0, 0)

根据上几个图片我们不难发现当 γ \gamma γ取到一定的值的时候我们发现相轨图突然由混沌态转为周期态,而 γ \gamma γ的取值变化是非常细微的。本周学会了对混沌性的duffing方程中的图像进行绘制,下周继续深度的学习混沌性。

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值