逆滤波器设计的正则化

信息处理 逆问题
2022-01-31 12:18:33

给定一个2×2矩阵,C,假设我想计算一个过滤矩阵H=C1并且我需要出于实际目的添加正则化(例如,对于音频滤波器,需要正则化以限制滤波器在某些频率下的增益,其中条件数为C很大。)求解的线性方程组是CH=I, 在哪里I是个2×2单位矩阵。

我的方法是使用 Tikhonov 正则化,其中解决方案(在每个频率上)由下式给出H=(CC+βI)1CI, 在哪里β是一个标量正则化参数,* 表示采用复共轭,因为这会最小化目标函数CHI2+βH2. So, with an appropriately chosen β , we can essentially constrain the filter gain (since H effectively corresponds to the max. gain of H ), albeit at the expense of H satisfying CH=I less optimally compared to using β=0

另一种方法可能是将 H 计算HH H=(detC)|detC|2+βadj(C),其中adj(C)对应于 C 的对数\det CC于它的行列式。请注意,对于\beta = 0,这恰好对应于C^{-1}这是我最近遇到但以前没有见过的一种方法。detCβ=0C1

我很确定这两种解决方案不等效,但如果它们是相同的,那么将不胜感激。假设它们不是,我一直在尝试确定采用一种方法相对于另一种方法的好处。一方面,后者只能在C是平方的情况下完成(并且,从计算的角度来看,不能太大,因为它需要显式计算行列式和佐证)。此外,β以不同的方式进入解决方案,因此在一种方法中优化β的值可能比另一种方法更容易。除此之外,我对两者之间的差异没有太多了解,因为我对后一种方法并不熟悉。

作为社区的回应,我正在寻找关于这两种方法是否不同的确认,以及关于何时可能想要使用一种方法而不是另一种方法的评论,包括对后一种方法的更直观的解释,如果可能的话。使用后一种方法的论文/文章/书籍的链接将被计算在内并受到赞赏。谢谢!

1个回答

这种替代方法对我来说似乎很可疑。的二维情况,其特征值是实数和正数,并且你有 你想用正则化防止的是你的逆得到一个太大的特征值,即你想把它限制在最多同时,很好。但是,如果您将相应的特征向量放入该替代公式中,您会得到 Cλ0,λ1λ0λ1

detC=λ0λ1
λ01β1λ11v1
Hv1=detC|detC|2+βadjCv1=λ02λ12λ02λ12+βC1v1=λ02λ1λ02λ12+βv1λ02λ1βv1
所有输入的缩放取决于,而不仅仅是低特征值的。β

我想这个公式在某些情况下可能仍然有效(也许行列式大多是虚构的?),但它肯定不适合一般情况。

使用吉洪诺夫。