码迷,mamicode.com
首页 > Windows程序 > 详细

基于winner 滤波平稳降噪效果

时间:2015-11-26 01:18:19      阅读:349      评论:0      收藏:0      [点我收藏+]

标签:

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

技术分享

Where 技术分享 are spectra. Provided that 技术分享 is optimal, then the minimum mean-square error equation reduces to

技术分享

and the solution 技术分享 is the inverse two-sided Laplace transform of 技术分享.

Causal solution

技术分享

where

  • 技术分享 consists of the causal part of 技术分享 (that is, that part of this fraction having a positive time solution under the inverse Laplace transform)
  • 技术分享 is the causal component of 技术分享 (i.e., the inverse Laplace transform of 技术分享 is non-zero only for 技术分享)
  • 技术分享 is the anti-causal component of 技术分享 (i.e., the inverse Laplace transform of 技术分享 is non-zero only for 技术分享)

This general formula is complicated and deserves a more detailed explanation. To write down the solution 技术分享 in a specific case, one should follow these steps:[2]

  1. Start with the spectrum 技术分享 in rational form and factor it into causal and anti-causal components:
    技术分享
    where 技术分享 contains all the zeros and poles in the left half plane (LHP) and 技术分享 contains the zeroes and poles in the right half plane (RHP). This is called the Wiener–Hopf factorization.
  2. Divide 技术分享 by 技术分享 and write out the result as a partial fraction expansion.
  3. Select only those terms in this expansion having poles in the LHP. Call these terms 技术分享.
  4. Divide 技术分享 by 技术分享. The result is the desired filter transfer function 技术分享.

 

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

技术分享

 

逐帧实时维纳滤波后

技术分享

 

基于winner 滤波平稳降噪效果

标签:

原文地址:http://www.cnblogs.com/welen/p/4996310.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!