10 有限长脉冲响应滤波器和窗函数

10.1 理想滤波器及其存在的问题

  1. Q: 对于抽样间隔\(\Delta\),理想低通滤波器频谱范围为(),但时域出现了()的问题。
    A: \([-1/2\Delta,1/2\Delta]\),需要无限长的信号
    注:这里的频谱其实都相当于周期函数。
  2. Q: 既然时域无限长,那我们直接频域逐点相乘,出现什么问题?
    A: 注意上面的“注”。这里实际上是“有无穷多个带”的带通滤波了。这些“带”按周期排列。
  3. Q: 某某“通”与某某“阻”滤波器的频谱有何关系?时间函数有何关系?为什么?
    A: 相加之后频谱为常数。相加之后时间函数为常数倍\(\delta\)函数。
  4. Q: 如何计算理想微分滤波器的时间函数?
    A: 回忆频谱为\(i2\pi f\),从而容易计算出时域\(h(n) = \int_{-1/2\Delta}^{1/2\Delta}i2\pi fe^{i2\pi n\Delta f}df\)
  5. Q: 吉布斯现象的原因中为什么有一条“频谱突跳”?
    A: 例如:可以从时间处“截断”,频域卷积\(sinc\)理解。如果没有“突跳”,自己拿“很像\(\delta\)函数的”\(sinc\)卷一下就知道没有Gibbs现象了。
    直观理解:sinc有正有负,中间“主峰”多出的“正”预备着外面“还掉”。但你要是频谱突变,那就可能导致“不用还”,“耍赖”。总之没有突跳就没有此现象。

近似理想滤波器的技术要求

  1. Q: 解说近似理想滤波器技术指标中的关键数值。
    A: 过滤带\([\omega_p,\omega_s]\):其中“不做要求”
    波动\(\delta_p,\delta_s\):“做要求”的频段出现波动,不是严格直线。
  2. Q: 分贝和数字“20”有何关系?
    A: 2是振幅和能量的关系,10对应“分”

10.2 时窗函数

  1. Q: “主瓣”和“旁瓣”的特点和截尾时间函数频谱的畸变形状有何联系?
    A: 主瓣不尖,跃变就不突然。
    旁瓣不可忽略,本来平的地方就出现波动。
  2. Q: 如何理解“主瓣宽度小”和“旁瓣水平低”矛盾?
    A: 实际操作中主瓣高度一定,其宽度小时,旁瓣水平就高。
  3. Q: 钟形时窗\(w_3(t) = e^{-\alpha (t/T)^2}\)
    的频谱近似为(),
    A: \(W_3(f)\approx \frac{\sqrt \pi T}\alpha e^{-\frac {\pi^2 T^2 f^2}{\alpha}}\)
  4. Q: 选取滤波器:先根据需要的()选择(),再根据过滤带宽度(和滤波器长度\(N\)有关)选择()
    A: 阻带衰减 旁瓣峰值,时窗函数,滤波器长度\(N\)
  5. Q: 解说\(h=\hat h w\)的符号含义
    A: 想求实际可用的滤波器时间函数\(h\),我们需要先求出理想的\(\hat h\)(这里要求考虑过渡带中心是截频,以及物理可实现)
  6. Q: 求最大振幅比时窗函数的\(\delta\)参数对求解结果有何影响?
    A: 影响\(cos2\pi f\Delta\)取值范围(\(-1\le x\le a\)),影响优化目标,从而改变切比雪夫多项式的系数。
  7. Q: 简述“最大能量比”和“特征根”的联系
    A: 能量本身是二次,求导后是一次。\(U/V\)求导为0得到\(\partial V/\partial h_k = 常数\cdot \partial U/\partial h_k\),这就有了线性方程。

11 递归滤波器的设计

11.1 递归滤波及其稳定性

  1. Q: 如何理解此处“递归”一词?
    A: 其实是带“记忆”的递归,“空间换时间”,利用之前计算的结果,减少计算量。和算法设计中的带记忆的递归有类似之处。
  2. Q: 如何把\(y(t) = b_0x(t)+\cdots + b_n x(t-n) - a_1y(t-1)-\cdots - a_ny(t-n)\)表示成\(Z\)变换的关系?
    A: 首先注意\(b_kx(t-k)\)对应\(b_kZ^k X(Z)\),其次可以形式上记\(A(Z)=\sum_{k=0}^ma_kZ^k,a_0=1\). 这样容易写出形如\(H(Z) = Y(Z)/X(Z) = B(Z)/A(Z)\)\(H\)就是滤波器的\(Z\)变换。
  3. Q: 误差,稳定和差分方程有何关系?
    A: 可以用差分方程求出误差的通项公式。如果其中出现了模长大于1的底数,就不稳定。
    误差跟\(x\)那边没关系,所以只用考察\(y\)这边(即\(A(Z)\))
    \(A\)根都在单位圆外,倒数都在单位圆内。再回忆差分方程求通项公式时底数是怎么来的,就知道了。
  4. Q: 反向递归滤波中2.对应的结论为什么没有“相反”成“单位圆内”?
    A: 人为定义问题。这里反向滤波时的\(A(Z)\)仍然定义“次数越高,离前线越远”。这点和正向一致。
  5. Q: \(1/(1-\alpha Z)\)\(|\alpha|>1\)时,\(Z\)的根在单位圆内,那怎么找稳定的递归滤波器?那如果分母多项式有的根在单位圆外,有的根在单位圆内呢?
    A: 看\(1/Z\)
    进一步地,如果分母多项式有的根在单位圆外,有的根在单位圆内,则可以串联(相乘)或并联(加减),分两部分处理。一部分正向,一部分反向。
    串联是先算中间信号,再算最终结果。每步都稳定
    并联是分两部分,每部分都稳定
    并联加法的意义?\(H_1-H_2\)变为\(H_1+H_2\)

11.2 模拟滤波器的设计

  1. Q: 为什么不直接设计\(H(\omega)\),而是通过\(|H(\omega)|^2=\)()中的多项式()来设计?
    A: \(\frac 1{1+g^2(\omega)}\)\(g(\omega)\)
    \(g(\omega)\)的性质(圆内很小,圆外很大)很容易构造,例如\(\omega^n\),切比雪夫……等。
  2. Q: 接上,\(|H(\omega)|^2\)怎么求出\(H\)
    A: 复数域上展开,每个\((\omega-\alpha_i)(\omega-\bar \alpha_i)\)对中取一个。
    (注:如果让下半平面无零点、极点,则具有一些特殊性质)