混沌系统概述

混沌系统概述

目录

  1. Lorenz 系统（洛伦兹系统）
  2. Rössler 系统（罗斯勒系统）
  3. Chen 系统（陈氏系统）
  4. Chua 电路（蔡氏电路）
  5. Duffing 振子（杜芬振子）
  6. Logistic 映射（逻辑斯蒂映射）
  7. Hénon 映射（埃农映射）
  8. Mackey-Glass 系统
  9. Van der Pol 振子
  10. Ikeda 映射（池田映射）

1. Lorenz 系统（洛伦兹系统）

简介

Lorenz 系统由美国气象学家 Edward Lorenz 于 1963 年提出，是混沌理论中最著名的连续动力系统。Lorenz 在研究大气对流模型时发现了这个系统，它展示了对初始条件极其敏感的”蝴蝶效应”。

系统方程

微分方程组：

{x˙=σ(y−x)y˙=x(ρ−z)−yz˙=xy−βz\begin{cases} \dot{x} = \sigma(y - x) \\ \dot{y} = x(\rho - z) - y \\ \dot{z} = xy - \beta z \end{cases}

典型参数： σ=10\sigma = 10, ρ=28\rho = 28, β=8/3\beta = 8/3

混沌吸引子： Lorenz 吸引子形状像一只蝴蝶，具有两个对称的”翅膀”。

2. Rössler 系统（罗斯勒系统）

简介

Rössler 系统由德国化学家 Otto Rössler 于 1976 年提出。与 Lorenz 系统相比，Rössler 系统只有一项非线性项，但其混沌行为同样丰富。该系统产生著名的”螺旋型混沌”。

系统方程

微分方程组：

{x˙=−y−zy˙=x+ayz˙=b+z(x−c)\begin{cases} \dot{x} = -y - z \\ \dot{y} = x + ay \\ \dot{z} = b + z(x - c) \end{cases}

典型参数： a=0.2a = 0.2, b=0.2b = 0.2, c=5.7c = 5.7

3. Chen 系统（陈氏系统）

简介

陈氏系统由陈关荣教授于 1999 年提出，是一个典型的反控制混沌系统。它通过简单的状态反馈控制器从线性系统中产生混沌，是混沌化（chaotification）理论的重要成果。

系统方程

微分方程组：

{x˙=a(y−x)y˙=(c−a)x−xz+cyz˙=xy−bz\begin{cases} \dot{x} = a(y - x) \\ \dot{y} = (c - a)x - xz + cy \\ \dot{z} = xy - bz \end{cases}

典型参数： a=35a = 35, b=3b = 3, c=28c = 28

4. Chua 电路（蔡氏电路）

简介

蔡氏电路由蔡少棠（Leon Chua）于 1983 年提出。它是第一个能够通过物理实验验证混沌行为的电路系统，也是最简单的混沌电路之一。蔡氏电路能够产生著名的”双涡卷”吸引子。

系统方程

无量纲化微分方程组：

{x˙=α[y−x−f(x)]y˙=x−y+zz˙=−βy\begin{cases} \dot{x} = \alpha[y - x - f(x)] \\ \dot{y} = x - y + z \\ \dot{z} = -\beta y \end{cases}

其中非线性函数：

f(x)=m1x+12(m0−m1)(∣x+1∣−∣x−1∣)f(x) = m_1 x + \frac{1}{2}(m_0 - m_1)(|x+1| - |x-1|)

典型参数： α=15.6\alpha = 15.6, β=28\beta = 28, m0=−1.143m_0 = -1.143, m1=−0.714m_1 = -0.714

5. Duffing 振子（杜芬振子）

简介

Duffing 振子是一个描述受迫阻尼谐振子的非线性二阶微分方程。在特定参数下，Duffing 振子表现出混沌行为。

系统方程

二阶微分方程：

x¨+δx˙−x+x3=γcos⁡(ωt)\ddot{x} + \delta \dot{x} - x + x^3 = \gamma \cos(\omega t)

一阶微分方程组形式：

{x˙=yy˙=x−x3−δy+γcos⁡(ωt)\begin{cases} \dot{x} = y \\ \dot{y} = x - x^3 - \delta y + \gamma \cos(\omega t) \end{cases}

典型参数： δ=0.25\delta = 0.25, γ=0.3\gamma = 0.3, ω=1.0\omega = 1.0

6. Logistic 映射（逻辑斯蒂映射）

简介

Logistic 映射是最经典的一维离散混沌映射，常用于描述种群增长模型。虽然形式简单，但它能展示出丰富的动力学行为，包括倍周期分岔通向混沌。

系统方程

差分方程：

xn+1=rxn(1−xn)x_{n+1} = rx_n(1 - x_n)

参数范围： xn∈[0,1]x_n \in [0, 1], r∈[0,4]r \in [0, 4]

混沌区域： 当 r∈[3.57,4]r \in [3.57, 4] 时系统进入混沌状态

7. Hénon 映射（埃农映射）

简介

Hénon 映射由法国天文学家 Michel Hénon 于 1976 年提出，是对 Lorenz 系统 Poincaré 截面的简化二维离散模型。它是研究耗散系统中奇怪吸引子的重要模型。

系统方程

差分方程组：

{xn+1=1−axn2+ynyn+1=bxn\begin{cases} x_{n+1} = 1 - ax_n^2 + y_n \\ y_{n+1} = bx_n \end{cases}

典型参数： a=1.4a = 1.4, b=0.3b = 0.3

8. Mackey-Glass 系统

简介

Mackey-Glass 系统由 Mackey 和 Glass 于 1977 年提出，是描述生理系统中（如血红细胞生成）时滞反馈调节的模型。该系统的特点在于含有时间延迟项，能够产生高维混沌。

系统方程

时滞微分方程：

x˙(t)=ax(t−τ)1+xc(t−τ)−bx(t)\dot{x}(t) = \frac{ax(t-\tau)}{1 + x^c(t-\tau)} - bx(t)

典型参数： a=0.2a = 0.2, b=0.1b = 0.1, c=10c = 10, τ=17\tau = 17

9. Van der Pol 振子

简介

Van der Pol 振子由荷兰电气工程师 Balthasar van der Pol 于 1920 年提出。它是一个具有非线性阻尼的自激振荡器，在受迫激励下可以表现出混沌行为。

系统方程

二阶微分方程（受迫振动）：

x¨−μ(1−x2)x˙+x=Fcos⁡(ωt)\ddot{x} - \mu(1 - x^2)\dot{x} + x = F\cos(\omega t)

一阶微分方程组形式：

{x˙=yy˙=μ(1−x2)y−x+Fcos⁡(ωt)\begin{cases} \dot{x} = y \\ \dot{y} = \mu(1 - x^2)y - x + F\cos(\omega t) \end{cases}

典型参数： μ=5\mu = 5, F=5F = 5, ω=2.466\omega = 2.466

10. Ikeda 映射（池田映射）

简介

Ikeda 映射由 Kensuke Ikeda 于 1979 年提出，用于描述光学环形腔中的光传输动力学。该映射由复变量描述，可以转化为二维实映射。

系统方程

复变量形式：

zn+1=A+Bznei(∣zn∣2+C)z_{n+1} = A + Bz_n e^{i(|z_n|^2 + C)}

实变量二维映射形式：

{xn+1=A+B(xncos⁡ϕ−ynsin⁡ϕ)yn+1=B(xnsin⁡ϕ+yncos⁡ϕ)\begin{cases} x_{n+1} = A + B\big(x_n\cos\phi - y_n\sin\phi\big) \\ y_{n+1} = B\big(x_n\sin\phi + y_n\cos\phi\big) \end{cases}

其中：ϕ=C−D1+xn2+yn2\phi = C - \dfrac{D}{1 + x_n^2 + y_n^2}

典型参数： A=1.0A = 1.0, B=0.9B = 0.9, C=0.4C = 0.4, D=6.0D = 6.0

混沌系统的数值求解方法

连续混沌系统（微分方程）的仿真通常采用以下数值方法：

方法精度特点
Euler 法一阶最简单，但精度低
Runge-Kutta 2 (中点法)二阶速度和精度折中
Runge-Kutta 4 (RK4)四阶最常用，精度高
Dormand-Prince (ode45)四五阶MATLAB 默认求解器

离散映射系统直接迭代即可，无需数值积分。

推荐：在仿真网站中使用 RK4 方法 求解连续混沌系统，离散系统采用直接迭代。

← 返回混沌系统仿真