34、实数域上高效负循环卷积算法解析

实数域上高效负循环卷积算法解析

在信号处理、密码学等众多领域中，卷积运算，特别是负循环卷积，扮演着至关重要的角色。然而，传统的计算方法存在效率低下或精度不足的问题。本文将详细介绍一种高效的实数域负循环卷积算法，深入分析其原理、误差传播，并通过实验验证其性能。

1. 基础知识回顾

在深入探讨负循环卷积算法之前，我们需要回顾一些基本的数学概念。

• 循环卷积与负循环卷积 ：
  设 (f, g \in C^N)，其中 (N \in N)。经典的循环卷积定义为 ((f * g) k := \sum {j = 0}^{N - 1} f_j g_{(k - j) \bmod N})，而负循环卷积则在循环索引每次环绕时添加一个 (-1) 的因子，即 ((f \bar{*} g) k := \sum {j = 0}^{N - 1} (-1)^{\lfloor \frac{k - j}{N} \rfloor} f_j g_{(k - j) \bmod N})。从多项式的角度来看，循环卷积计算的是两个多项式模 (X^N - 1) 的乘积系数，而负循环卷积计算的是两个多项式模 (X^N + 1) 的乘积系数。
• 卷积定理 ：
  卷积定理（CT）表明，卷积向量的傅里叶变换等于它们各自傅里叶变换的逐元素（二元）乘积。即 (F(f * g) = F(f) \odot F(g))，其中 (F(\cdot)) 表示离散傅里叶变换（DFT），(\odot) 表示两个向量的二元乘法。DFT 是一种基变换，定义为 (F(f) k := \su