Applying the Spectral Method for Modeling Linear Filters: Bessel, Papoulis, and Legendre Filters

技术摘要：应用谱方法对线性滤波器进行建模

问题陈述
本文研究了用于信号处理的线性滤波器（具体为贝塞尔、帕普利斯和勒让德类型）的计算机仿真问题。传统的仿真方法通常需要从连续时间向离散时间转换（例如，通过脉冲不变性或双线性变换）。作者指出，这种转换会引入固有的缺陷，包括频率混叠和频率畸变（频率翘曲）。本研究的目标是开发一种完全在连续时间下运行的建模技术，从而避免这些离散化伪影，这在滤波器设计阶段特别有用。

方法论
所提出的技术利用谱方法来分析线性控制系统。该方法并非直接在时域求解微分方程，也不是转换为离散时间，而是将输入和输出信号表示为选定基函数（在数值实验中特指三角基函数）内展开系数的有序集合。

方法论的关键组成部分包括：

1. 谱表示： 信号 $g(t)$ 和 $x(t)$ 在正交基 $\{q(i, t)\}$ 中展开。从时域信号到其谱图像（系数构成的无限列矩阵）的转换记为 $S$。
2. 二维非平稳传递函数 ($W$)： 滤波器本身不是通过标准的传递函数 $H(s)$ 来建模，而是通过一个二维非平稳传递函数 $W$ 来建模。该矩阵通过代数关系 $X = WG$ 将输入谱系数与输出谱系数联系起来。
3. 直接推导： 该方法的一个核心优势在于，$W$ 可以直接从滤波器的有理传递函数 $H(s)$（或其对应的微分方程）推导得出，而无需先计算脉冲响应。矩阵 $W$ 被表示为矩阵 $P$（代表微分器）或 $P^{-1}$（代表积分器）的函数。
   • 对于由 $H(s) = M(s)/D(s)$ 定义的滤波器，谱矩阵为 $W(P) = D^{-1}(P)M(P)$。
4. 基函数选择： 虽然在谱域中的滤波器描述具有通用性，但数值实验采用了三角基函数。论文提供了矩阵 $P$、$P^{-1}$ 的显式公式，以及用于误差计算所需的时移矩阵 ($S_\tau$) 和区间限制矩阵 ($A_\theta$)。

主要贡献
本文将谱方法的应用扩展到了三种特定的正交多项式线性滤波器家族，此前这些滤波器尚未使用这种特定的连续时间谱技术进行建模：

1. 贝塞尔滤波器： 本文利用反向贝塞尔多项式推导了贝塞尔滤波器（及广义贝塞尔滤波器）的二维非平稳传递函数。文中提供了这些滤波器的显式递推关系和矩阵形式。
2. 帕普利斯滤波器： 详细介绍了帕普利斯滤波器（最优“L”滤波器）的综合过程。这涉及对源自勒让德多项式的多项式进行谱分解。论文概述了寻找必要根并构建相应矩阵 $W$ 的步骤。
3. 勒让德滤波器： 该方法被应用于勒让德滤波器，这类滤波器利用修正后的相关勒让德多项式。文中详细说明了涉及基于这些多项式的方程进行谱分解的综合过程。

此外，这项工作简要展示了如何将该技术扩展到修正贝塞尔、哈尔彭 (Halpern) 和超球面 (Ultraspherical) 滤波器。论文还提供了必要的数学工具，包括相对于三角基函数的初等模块（微分器、积分器、时移）的非平稳传递函数的显式元素。

结果
通过涉及混合噪声的确定性信号滤波的数值实验验证了该技术。

• 设置： 有用信号 $u(t) = \sin(5\pi t)$ 与噪声 $v(t) = 0.1(\cos(85\pi t) + \sin(105\pi t))$ 混合。滤波器在不同的截止频率和不同阶数（$n=2$ 到 $6$）下进行了测试。
• 误差分析： 通过衡量滤波输出与原始有用信号之间的 $L_2$ 误差（考虑相位延迟）来评估滤波质量。
• 发现：
  • 贝塞尔滤波器： 在所有测试类型中准确度最低，这归因于其幅度滚降较慢。误差主要受滤波器固有特性而非谱方法截断误差的影响。
  • 帕普利斯滤波器： 表现出高准确度，其中偶数阶滤波器的表现通常优于奇数阶。误差显示出对截断阶数 $L$ 的依赖性，表明滤波器特性和谱近似误差共同影响了结果。
  • 勒让德滤波器： 通常优于帕普利斯滤波器，特别是当参数 $m=2$ 时。与帕普利斯滤波器类似，较高的阶数需要更大的截断阶数以最小化谱近似误差。
  • 收敛性： 增加截断阶数 $L$（矩阵大小）降低了谱方法误差，证实了谱近似的收敛性。

意义与主张
本文声称提出了一种用于连续时间线性滤波器计算机仿真的新方法。其主要意义在于：

1. 避免离散化： 它消除了向离散时间转换的需求，从而避免了与双线性变换等方法相关的频率混叠和频率畸变。
2. 泛化性： 它成功地将之前应用于巴特沃斯、林克威茨-里利 (Linkwitz–Riley) 和切比雪夫滤波器的数学描述谱形式，扩展到了贝塞尔、帕普利斯和勒让德滤波器。
3. 计算框架： 它为模拟这些特定滤波器家族提供了一个完整的代数框架（矩阵表示），包括推导必要的非平稳传递函数以及在谱域内处理相位延迟和误差计算。

作者总结认为，这种方法在滤波器设计阶段特别有用，并且只要选择合适的基函数（如勒让德多项式），就可以推广到分数阶滤波器。这项工作作为现有线性滤波器理论文献的补充，提供了一种针对标准离散时间仿真方法的连续时间替代方案。