un

guest
1 / ?
back to lessons

e^{i2πf}沿单位圆运动

一个复指数e^{iθ}位于复平面上的单位圆上。随着θ的增加,点逆时针旋转。

对于在整数时间n = 0, 1, 2, …采样的数字滤波器,特征函数e^{i2πfn}在每个采样时刻沿圆周运动一个角度2πf。

频率作为旋转速度:f表示每个采样发生满额volution的比例。

- f = 0:没有旋转;点保持在(1, 0)位置

- f = 1/4:每步旋转1/4圆

- f = 1/2:每步旋转1/2圆(纽奎斯特频率)

- f = 1:每步旋转满圆——与f = 0不可区分

这个最后一点包含了几何上的全部抖动故事。

频率作为角度:单位圆与抖动

为什么是单位圆

单位圆是所有|z| = 1的集合。计算Z-变换H(z)的结果在单位圆上——设置z = e^{i2πf}——给出频率响应H(f)。单位圆是离散时间稳定性和频率分析的交汇点。

角度与频率

每个频率f对应于一个角度θ = 2πf弧度/采样。所有不同的频率的范围跨越一个完整的革命:f ∈ [0, 1)或等价地θ ∈ [0, 2π)。

在纽奎斯特频率f = 1/2时,每个采样前进正好π弧度——半个革命。

在单位圆上,频率f = 1/6对应于角度2π/6 = π/3弧度。一个滤波器在其传递函数中有一个零点z = e^{i2π/3}(即在f = 1/3)。解释在输入中包含频率f = 1/3时,滤波器输出发生了什么几何变化。为什么将零点放在单位圆上产生完全消除?

抖音的几何图像

单位圆的周长为2π。一个完整的旋转对应于频率f = 1(每个样本一个完整周期)。在采样信号中的不同频率恰好占据一个旋转。

当f = 1/2 + δ时,发生什么?每个样本的旋转 = 2π(1/2 + δ) = π + 2πδ。经过k个样本后,角度 = k(π + 2πδ)。但角度π + 2πδ与角度-π + 2πδ几何上相同,后者对应于频率f = 1/2 - δ的旋转。

抖音是圆周上的模数算术。 高于纠缠频率的频率会回绕。圆圈对过去的方向没有记忆。

采样定理说:保持在半圆[0, π)内。确保信号永远不会达到另一半。反向抖音滤波器执行此边界。

计算抖音的几何方法

在采样率f_s下,频率f的别名出现在|f − round(f / f_s) · f_s|处——这是一个对f_s的最近整数的距离,表示为一个分数。

对于f_s = 1(归一化):别名f = 1 − f,对于f ∈ (1/2, 1)。这是关于纠缠点f = 1/2的反射。

几何上:f & 1 − f在单位圆上处于镜像位置,同样距离π轴远。

一个采样频率为1000 Hz的信号包含一个700 Hz的音调。纠缠频率为500 Hz。使用几何上的反射论证——抖音频率在纠缠点上单位圆的反射处计算别名频率。然后在一个单位圆上绘制700 Hz和其别名的位置,该圆的全周代表1000 Hz采样周期。

幅度响应作为距离乘积

对于传递函数H(z)具有零点z_1, z_2, …和极点p_1, p_2, …:

|H(f)| = (∏ |e^{i2πf} − z_k|) / (∏ |e^{i2πf} − p_k|)

这是图形方法,直接从极零图读取频率响应。

规则

- 在单位圆上的零点在该频率上产生完全消除。

- 在单位圆附近的极点在响应中产生峰值。

- 在单位圆附近(但不在它上面)的零点产生低峰,不是消除。

- 在单位圆内的极点使滤波器稳定。

Z平面几何表示了整个滤波器行为的可视化。工程师在计算系数之前绘制极零图。

一个二阶滤波器的两个零在z = ±j处(在单位圆上频率为f = 1/4和f = 3/4),还有一个极点在z = 0.7处(在实轴上,单位圆内)。不计算任何系数:描述其频率响应形状。哪些频率通过?哪些被消除?响应峰值在哪里?用几何方法为每个主张提供证据。