技术笔记测试站

Laplace 变换

Laplace 变换将时域信号 $f(t)$ 映射到复频域 $F(s)$, 是分析线性时不变系统的核心工具:

$$F(s) = \int_{0}^{\infty} f(t) \, e^{-st} \, dt, \quad s = \sigma + j\omega$$

复变量 $s$ 的实部 $\sigma$ 决定包络的指数衰减/增长, 虚部 $\omega$ 决定振荡频率。极点(使 $F(s) \to \infty$ 的 $s$)的位置完全刻画了系统的暂态行为:极点在左半平面 → 稳定衰减;在虚轴上 → 等幅振荡;在右半平面 → 发散。

例子:阻尼正弦

取 $f(t) = e^{-at} \sin(\omega t) \, u(t)$($a > 0$), 其 Laplace 变换为:

$$F(s) = \frac{\omega}{(s + a)^2 + \omega^2}$$

它有一对共轭极点:

$$s = -a \pm j\omega$$

实部 $-a$ 控制衰减速率,虚部 $\omega$ 控制振荡频率。

交互式探索

拖动滑块改变衰减率 $a$ 与频率 $\omega$, 观察三者如何联动:

• 时域波形:衰减快慢与振荡频率;
• s 平面极点:极点位置随参数移动(注意实部恒为 $-a$, 虚部为 $\pm\omega$);
• 3D 复频域曲面:$|F(\sigma + j\Omega)|$, 极点处出现尖峰。

复频域三维曲面

下面的曲面是 $|F(\sigma + j\Omega)|$ 在整个复平面上的取值(松开滑块后更新, 因计算量较大)。两个对称的尖峰正对应极点 $s = -a \pm j\omega$。可以拖动旋转观察。

观察

• 把 $a$ 调到很小(如 $0.1$), 时域振荡衰减很慢, 极点靠近虚轴, 3D 曲面的峰又高又尖——系统接近临界稳定;
• 把 $a$ 调大(如 $1.5$), 时域快速衰减, 极点远离虚轴移向左侧, 峰变矮变宽——系统阻尼强;
• 极点实部始终为负($-a < 0$), 故系统始终稳定;若极点进入右半平面, 系统将发散。

Logistic 映射

Logistic 映射是一个简单到难以置信的一维离散动力系统, 却蕴含着通往混沌的全部复杂性。它的迭代规则为:

$$x_{n+1} = r \, x_n (1 - x_n), \quad x_n \in [0,1], \quad r \in [0,4]$$

其中 $r$ 是增长参数, $x_n$ 是第 $n$ 代种群(归一化后)。虽然只有一个参数和一行递推, 它却能展现出 不动点收敛 → 周期振荡 → 倍周期分岔 → 混沌 的完整光谱。

不动点与稳定性

令 $x_{n+1} = x_n = x^*$, 得到两个不动点: $x^* = 0$ 与 $x^* = 1 - 1/r$(后者当 $r>1$ 时存在)。稳定性由 $|f'(x^*)|$ 决定, 其中 $f(x) = rx(1-x)$:

$$f'(x) = r(1 - 2x)$$

• $x^* = 0$ 在 $r < 1$ 时稳定;
• $x^* = 1 - 1/r$ 在 $1 < r < 3$ 时稳定;
• 当 $r > 3$, 不动点失稳, 系统进入倍周期分岔, 并在 $r \approx 3.5699$ 之后进入混沌。

交互式探索:时间序列与蛛网

拖动下面的滑块改变参数 $r$ 与初值 $x_0$, 实时观察迭代行为。建议依次尝试 $r = 2.5$(单点收敛)、$r = 3.2$(2 周期)、$r = 3.5$(4 周期)、$r = 3.7$(混沌)。

蛛网图说明: 横轴是 $x_n$, 纵轴是 $f(x_n)$。从 $(x_0, 0)$ 出发, 垂直上升到曲线 $y = f(x)$, 再水平移动到对角线 $y = x$, 就完成了一次迭代。蛛网的收敛、环绕、乱舞, 直观对应了系统的稳定、振荡与混沌。

分岔图:通往混沌之路

下图固定展示了 $r$ 从 $2.5$ 到 $4$ 时系统的长期行为(丢弃前 400 次瞬态迭代, 保留后 200 个 $x$ 值)。可以清晰看到倍周期分岔树, 以及混沌区域中突然出现的"周期窗口"。

试着把上方滑块的 $r$ 拖到分岔图的红竖线位置, 对比时间序列与分岔结构——这是理解混沌最直接的方式。混沌的标志性特征之一是对初值的敏感依赖:在 $r = 3.7$ 附近, 把 $x_0$ 改变 $0.01$, 时间序列会在若干步后完全分离。