基于winner 滤波平稳降噪效果

最新推荐文章于 2024-05-15 23:15:35 发布

welen_flying

最新推荐文章于 2024-05-15 23:15:35 发布

阅读量2.9k

点赞数

分类专栏：音频算法开发文章标签： winner filter

本文链接：https://blog.csdn.net/weiqiwu1986/article/details/50044267

版权

音频算法开发专栏收录该内容

29 篇文章 33 订阅

订阅专栏

https://en.wikipedia.org/wiki/Wiener_filter

Wiener filter solutions

The Wiener filter problem has solutions for three possible cases: one where a noncausal filter is acceptable (requiring an infinite amount of both past and future data), the case where a causal filter is desired (using an infinite amount of past data), and the finite impulse response (FIR) case where a finite amount of past data is used. The first c

ase is simple to solve but is not suited for real-time applications. Wiener's main accomplishment was solving the case where the causality requirement is in effect, and in an appendix of Wiener's book Levinson gave the FIR solution.

Noncausal solution

G(s) = \frac{S_{x,s}(s)}{S_x(s)}e^{\alpha s}.

Where $S$ are spectra. Provided that $g(t)$ is optimal, then the minimum mean-square error equation reduces to

E(e^2) = R_s(0) - \int_{-\infty}^{\infty}{g(\tau)R_{x,s}(\tau + \alpha)\,d\tau},

and the solution $g(t)$ is the inverse two-sided Laplace transform of $G(s)$ .

Causal solution

G(s) = \frac{H(s)}{S_x^{+}(s)},

where

$H(s)$ consists of the causal part of $\frac{S_{x,s}(s)}{S_x^{-}(s)}e^{\alpha s}$ (that is, that part of this fraction having a positive time solution under the inverse Laplace transform)
$S_x^{+}(s)$ is the causal component of $S_x(s)$ (i.e., the inverse Laplace transform of $S_x^{+}(s)$ is non-zero only for $t \,\ge\, 0$ )
$S_x^{-}(s)$ is the anti-causal component of $S_x(s)$ (i.e., the inverse Laplace transform of $S_x^{-}(s)$ is non-zero only for $t < 0$ )

This general formula is complicated and deserves a more detailed explanation. To write down the solution $G(s)$ in a specific case, one should follow these steps:^[2]

Start with the spectrum $S_x(s)$ in rational form and factor it into causal and anti-causal components:

$S_x(s) = S_x^{+}(s) S_x^{-}(s)$

where $S^{+}$ contains all the zeros and poles in the left half plane (LHP) and $S^{-}$ contains the zeroes and poles in the right half plane (RHP). This is called the Wiener–Hopf factorization.
Divide $S_{x,s}(s)e^{\alpha s}$ by $S_x^{-}(s)$ and write out the result as a partial fraction expansion.
Select only those terms in this expansion having poles in the LHP. Call these terms $H(s)$ .
Divide $H(s)$ by $S_x^{+}(s)$ . The result is the desired filter transfer function $G(s)$ .

原始文件，环境噪音已经很弱了

逐帧实时维纳滤波后

welen_flying

关注

0
点赞
踩
2

收藏

觉得还不错? 一键收藏
0
评论
基于winner 滤波平稳降噪效果

https://en.wikipedia.org/wiki/Wiener_filterWiener filter solutionsThe Wiener filter problem has solutions for three possible cases: one where a noncausal filter is acceptable (requiring an inf
复制链接

扫一扫

专栏目录