維納濾波器 （Wiener filter）是由美國數學家諾伯特·維納 （Norbert Wiener）於1940年代提出，並於1949年首次正式發表的一種線性濾波技術。它的設計目標是在最小均方誤差 平穩隨機過程
與理想化的通阻型濾波器（如低通 、高通 、帶通 或帶阻 ）不同，維納濾波器會根據訊號與雜訊的統計特性，動態調整對各頻段的保留或抑制程度。它具有兩種形式：
非因果版本 ：可以利用未來資訊來達成最佳去噪效果，但不適用於即時系統。因果版本 ：只使用當前與過去資料，以適應即時處理需求，但推導上更為複雜。維納濾波器廣泛應用於音訊降噪、影像復原、通訊無線通道均衡等領域；同時也成為多種自適應濾波演算法 如 LMS、RLS 及卡爾曼濾波器的重要理論基石。
維納濾波器的理論基礎可追溯至維納與霍夫（Norbert Wiener and Eberhard Hopf）對隨機訊號估計問題的研究，其推導出的維納–霍夫方程提供了解析解，使得在理想條件下可以精確求出濾波器的脈衝響應或係數。在離散時間域中，維納濾波器可以實作為有限長度（FIR ）或無限長度（IIR ）形式；而若進一步要求因果性，則必須透過頻譜因式分解處理，以確保濾波器可以即時實作。
與傳統設計為低通、高通、帶通等類型的頻率選擇性濾波器不同，維納濾波器不單以頻率為設計依據，而是依據輸入訊號與期望輸出的統計特性，於各頻段進行權重調整。其在頻域上的濾波器傳輸函數通常具有非零相位響應，這表示它可能會改變訊號的波形結構，但可保證輸出在統計意義下最接近目標訊號。
為解決實務中無法事先得知訊號與雜訊統計特性的限制，維納濾波器衍生出一系列自適應濾波方法，像是最小均方（LMS）與遞迴最小平方（RLS）演算法。這些方法可在線估測濾波係數，使濾波器能夠追蹤動態環境變化而維持近似最小均方誤差的性能。
除了用於一維訊號的去雜訊，維納濾波技術也已擴展至多維資料，例如影像處理中的降噪與去模糊。這類應用中，空間域與頻域的維納濾波器皆可使用，並可根據局部統計特性調整濾波強度，達到空間自適應的效果。
維納濾波器提供了一種將統計估計與線性系統設計結合的方式，並開啟了數位訊號處理中以機率模型為基礎的濾波方法發展，其理論影響深遠，至今仍是統計訊號處理與通訊系統中的重要工具。
假設維納濾波器的輸入信號是
s
(
t
)
{\displaystyle s(t)}
n
(
t
)
{\displaystyle n(t)}
x
(
t
)
{\displaystyle x(t)}
g
(
t
)
{\displaystyle g(t)}
x
(
t
)
=
g
(
t
)
∗
(
s
(
t
)
+
n
(
t
)
)
{\displaystyle x(t)=g(t)*(s(t)+n(t))}
s
(
t
)
{\displaystyle s(t)}
n
(
t
)
{\displaystyle n(t)}
x
(
t
)
{\displaystyle x(t)}
s
(
t
)
{\displaystyle s(t)}
g
(
t
)
{\displaystyle g(t)}
誤差是
e
(
t
)
=
s
(
t
+
d
)
−
x
(
t
)
{\displaystyle e(t)=s(t+d)-x(t)}
e
2
(
t
)
=
s
2
(
t
+
d
)
−
2
s
(
t
+
d
)
x
(
t
)
+
x
2
(
t
)
{\displaystyle e^{2}(t)=s^{2}(t+d)-2s(t+d)x(t)+x^{2}(t)}
s
(
t
+
d
)
{\displaystyle s(t+d)}
e
(
t
)
{\displaystyle e(t)}
根據d 的不同，問題名稱可以更換為：
如果
d
>
0
{\displaystyle d>0}
預測
如果
d
=
0
{\displaystyle d=0}
濾波
如果
d
<
0
{\displaystyle d<0}
平滑 將
x
(
t
)
{\displaystyle x(t)}
卷積 積分：
x
(
t
)
=
∫
−
∞
∞
g
(
τ
)
[
s
(
t
−
τ
)
+
n
(
t
−
τ
)
]
d
τ
{\displaystyle x(t)=\int _{-\infty }^{\infty }{g(\tau )\left[s(t-\tau )+n(t-\tau )\right]d\tau }}
計算平方誤差的均值 ，可得
E
(
e
2
)
=
R
s
(
0
)
−
2
∫
−
∞
∞
g
(
τ
)
R
x
s
(
τ
+
d
)
d
τ
+
∫
−
∞
∞
∫
−
∞
∞
g
(
τ
)
g
(
θ
)
R
x
(
τ
−
θ
)
d
τ
d
θ
{\displaystyle E(e^{2})=R_{s}(0)-2\int _{-\infty }^{\infty }{g(\tau )R_{x\,s}(\tau +d)d\tau }+\int _{-\infty }^{\infty }{\int _{-\infty }^{\infty }{g(\tau )g(\theta )R_{x}(\tau -\theta )d\tau }d\theta }}
其中
R
s
{\displaystyle \,\!R_{s}}
s
(
t
)
{\displaystyle s(t)}
自相關 函數
R
x
{\displaystyle \,\!R_{x}}
x
(
t
)
{\displaystyle x(t)}
自相關 函數
R
x
s
{\displaystyle R_{x\,s}}
x
(
t
)
{\displaystyle x(t)}
s
(
t
)
{\displaystyle s(t)}
互相關 函數如果信號
s
(
t
)
{\displaystyle s(t)}
n
(
t
)
{\displaystyle n(t)}
R
x
s
=
R
s
{\displaystyle R_{x\,s}=R_{s}}
R
x
=
R
s
+
R
n
{\displaystyle \,\!R_{x}=R_{s}+R_{n}}
這個的目的是求最優的
g
(
t
)
{\displaystyle g(t)}
E
(
e
2
)
{\displaystyle E(e^{2})}
[ 編輯 ] 維納濾波對於因果系統 與非因果系統 有兩種不同解，如下：
[ 編輯 ]
G
(
s
)
=
S
x
,
s
(
s
)
e
α
s
S
x
(
s
)
{\displaystyle G(s)={\frac {S_{x,s}(s)e^{\alpha s}}{S_{x}(s)}}}
其中
S
{\displaystyle S}
g
(
t
)
{\displaystyle g(t)}
E
(
e
2
)
=
R
s
(
0
)
−
∫
−
∞
∞
g
(
τ
)
R
x
,
s
(
τ
+
d
)
d
τ
{\displaystyle E(e^{2})=R_{s}(0)-\int _{-\infty }^{\infty }{g(\tau )R_{x,s}(\tau +d)d\tau }}
那麼方程的解
g
(
t
)
{\displaystyle g(t)}
G
(
s
)
{\displaystyle G(s)}
G
(
s
)
=
H
(
s
)
S
x
+
(
s
)
{\displaystyle G(s)={\frac {H(s)}{S_{x}^{+}(s)}}}
其中
H
(
s
)
{\displaystyle H(s)}
S
x
,
s
(
s
)
e
α
s
S
x
−
(
s
)
{\displaystyle {\frac {S_{x,s}(s)e^{\alpha s}}{S_{x}^{-}(s)}}}
S
x
+
(
s
)
{\displaystyle S_{x}^{+}(s)}
S
x
(
s
)
{\displaystyle S_{x}(s)}
S
x
−
(
s
)
{\displaystyle S_{x}^{-}(s)}
S
x
(
s
)
{\displaystyle S_{x}(s)}
離散序列FIR維納濾波器框圖。輸入信號 w [n ]與維納濾波器 g [n ]進行卷積，輸出結果與參考信號 s [n ]進行比較，得到濾波誤差 e [n ]. 因果的有限衝激響應（FIR）維納濾波器通過使用輸入和輸出信號的統計信息，而不是使用一些給定的數據矩陣X和輸出矢量Y，來發現最佳的抽頭權重。它將輸入信號的自相關的估計（T）代入輸入矩陣X，將輸出和輸入信號的互相關估計（V）代入的輸出向量Y。
為了獲得維納濾波係數，考慮將信號w[n]輸入到一個N階維納濾波器，係數
{
a
i
}
{\displaystyle \{a_{i}\}}
i
=
0
,
…
,
N
{\displaystyle i\,=\,0,\,\ldots ,\,N}
x
[
n
]
=
∑
i
=
0
N
a
i
w
[
n
−
i
]
{\displaystyle x[n]=\sum _{i=0}^{N}a_{i}w[n-i]}
剩餘誤差來表示e[n]的，被定義為e [n ] = x [n ] − s [n ]（參見相應的框圖）。維納濾波器被設計成最小化均方誤差（MMSE準則），它可以如下簡而言之：
a
i
=
arg
min
E
{
e
2
[
n
]
}
{\displaystyle a_{i}=\arg \min ~E\{e^{2}[n]\}}
E
{
⋅
}
{\displaystyle E\{\cdot \}}
a
i
{\displaystyle a_{i}}
共軛對稱 Toeplitz矩陣 ，而不是對稱 Toeplitz矩陣 。為簡單起見，下面只考慮所有這些量實數的情況。均方誤差（MSE）可改寫為：
E
{
e
2
[
n
]
}
=
E
{
(
x
[
n
]
−
s
[
n
]
)
2
}
=
E
{
x
2
[
n
]
}
+
E
{
s
2
[
n
]
}
−
2
E
{
x
[
n
]
s
[
n
]
}
=
E
{
(
∑
i
=
0
N
a
i
w
[
n
−
i
]
)
2
}
+
E
{
s
2
[
n
]
}
−
2
E
{
∑
i
=
0
N
a
i
w
[
n
−
i
]
s
[
n
]
}
.
{\displaystyle {\begin{array}{rcl}E\{e^{2}[n]\}&=&E\{(x[n]-s[n])^{2}\}\\&=&E\{x^{2}[n]\}+E\{s^{2}[n]\}-2E\{x[n]s[n]\}\\&=&E\{{\big (}\sum _{i=0}^{N}a_{i}w[n-i]{\big )}^{2}\}+E\{s^{2}[n]\}-2E\{\sum _{i=0}^{N}a_{i}w[n-i]s[n]\}.\end{array}}}
矢量
[
a
0
,
…
,
a
N
]
{\displaystyle [a_{0},\,\ldots ,\,a_{N}]}
a
i
{\displaystyle a_{i}}
∂
∂
a
i
E
{
e
2
[
n
]
}
=
2
E
{
(
∑
j
=
0
N
a
j
w
[
n
−
j
]
)
w
[
n
−
i
]
}
−
2
E
{
s
[
n
]
w
[
n
−
i
]
}
i
=
0
,
…
,
N
=
2
∑
j
=
0
N
E
{
w
[
n
−
j
]
w
[
n
−
i
]
}
a
j
−
2
E
{
w
[
n
−
i
]
s
[
n
]
}
.
{\displaystyle {\begin{array}{rcl}{\frac {\partial }{\partial a_{i}}}E\{e^{2}[n]\}&=&2E\{{\big (}\sum _{j=0}^{N}a_{j}w[n-j]{\big )}w[n-i]\}-2E\{s[n]w[n-i]\}\quad i=0,\,\ldots ,\,N\\&=&2\sum _{j=0}^{N}E\{w[n-j]w[n-i]\}a_{j}-2E\{w[n-i]s[n]\}.\end{array}}}
假設，
w
[
n
]
{\displaystyle w[n]}
s
[
n
]
{\displaystyle s[n]}
R
w
[
m
]
{\displaystyle R_{w}[m]}
R
w
s
[
m
]
{\displaystyle R_{ws}[m]}
w
[
n
]
{\displaystyle w[n]}
w
[
n
]
{\displaystyle w[n]}
s
[
n
]
{\displaystyle s[n]}
R
w
[
m
]
=
E
{
w
[
n
]
w
[
n
+
m
]
}
R
w
s
[
m
]
=
E
{
w
[
n
]
s
[
n
+
m
]
}
.
{\displaystyle {\begin{aligned}R_{w}[m]=&E\{w[n]w[n+m]\}\\R_{ws}[m]=&E\{w[n]s[n+m]\}.\end{aligned}}}
因此，注意到
R
w
s
[
−
i
]
=
R
s
w
[
i
]
{\displaystyle R_{ws}[-i]\,=\,R_{sw}[i]}
∂
∂
a
i
E
{
e
2
[
n
]
}
=
2
∑
j
=
0
N
R
w
[
j
−
i
]
a
j
−
2
R
s
w
[
i
]
i
=
0
,
…
,
N
{\displaystyle {\frac {\partial }{\partial a_{i}}}E\{e^{2}[n]\}=2\sum _{j=0}^{N}R_{w}[j-i]a_{j}-2R_{sw}[i]\quad i=0,\,\ldots ,\,N}
令該導數等於零，得到
∑
j
=
0
N
R
w
[
j
−
i
]
a
j
=
R
s
w
[
i
]
i
=
0
,
…
,
N
{\displaystyle \sum _{j=0}^{N}R_{w}[j-i]a_{j}=R_{sw}[i]\quad i=0,\,\ldots ,\,N}
上式可以改寫為矩陣形式
T
a
=
v
⇒
[
R
w
[
0
]
R
w
[
1
]
⋯
R
w
[
N
]
R
w
[
1
]
R
w
[
0
]
⋯
R
w
[
N
−
1
]
⋮
⋮
⋱
⋮
R
w
[
N
]
R
w
[
N
−
1
]
⋯
R
w
[
0
]
]
[
a
0
a
1
⋮
a
N
]
=
[
R
w
s
[
0
]
R
w
s
[
1
]
⋮
R
w
s
[
N
]
]
{\displaystyle {\begin{aligned}&\mathbf {T} \mathbf {a} =\mathbf {v} \\\Rightarrow &{\begin{bmatrix}R_{w}[0]&R_{w}[1]&\cdots &R_{w}[N]\\R_{w}[1]&R_{w}[0]&\cdots &R_{w}[N-1]\\\vdots &\vdots &\ddots &\vdots \\R_{w}[N]&R_{w}[N-1]&\cdots &R_{w}[0]\end{bmatrix}}{\begin{bmatrix}a_{0}\\a_{1}\\\vdots \\a_{N}\end{bmatrix}}={\begin{bmatrix}R_{ws}[0]\\R_{ws}[1]\\\vdots \\R_{ws}[N]\end{bmatrix}}\end{aligned}}}
這些方程被稱作Wiener-Hopf方程。出現在方程中的矩陣T是一個對稱的Toeplitz矩陣。當
R
{\displaystyle R}
a
=
T
−
1
v
{\displaystyle \mathbf {a} \,=\,\mathbf {T} ^{-1}\mathbf {v} }
Levinson-Durbin 算法，因此並不需要對矩陣
T
{\displaystyle \mathbf {T} }
維納濾波器的設計優勢在於：它能以最小化均方誤差（MSE）為目標，達成線性估計的最優效果，這使其在理想統計條件下（訊號與雜訊的功率譜已知）能充分發揮精準去噪與訊號恢復性能；此外，在某些案例中，濾波器係數可以求得封閉形式解，因此能利用簡單的物理元件（RC、LC網絡）加以實現。因其理論明確、可計算性高，維納濾波常成為後續 LMS、RLS 及卡爾曼濾波自適應方法的基礎。
然而，其實務上的限制也不容忽視：首先，在真實環境中精確獲得信號和噪聲的功率頻譜或自相關函數十分困難，尤其在非平穩或時變噪聲條件下更難以估計到位；其次，維納濾波為線性方法，這意味其無法處理非線性現象（如訊號高次諧波、混疊），可能導致濾波後結果出現模糊、振鈴（ringing）等失真現象；第三，設計 FIR 或 IIR 濾波器時需要進行矩陣反演或頻譜因式分解，對計算能力要求較高，尤其在高維或大規模案例中亦不易即時實現。
此外，在應用中也遇到頻率分佈重疊與生態噪聲模式的挑戰，維納濾波解對小幅頻譜估計誤差通常十分敏感，實務效能可能因此明顯下降。而且在一些影像去模糊案例中，若參數估計不精準，更可能導致「過度鈍化」（over-smoothing）或不必要的邊緣損失。
綜上所述，儘管維納濾波在統計理論上具備無可匹敵的最佳特性，但其使用條件限制（已知統計特性、線性假設、計算開銷）與對估計誤差敏感性，使其在較複雜或動態環境中常被替代或與其他自適應、非線性技術（如 LMS、卡爾曼或小波濾波）結合使用，以求更穩健且實用的處理效果。
^ Extrapolation, Interpolation, and Smoothing of Stationary Time Series . New York: Wiley. ISBN 0-262-73005-7 ^ Introduction to Random Signals and Applied Kalman Filtering . 3 ed. New York: John Wiley & Sons. ISBN 0-471-12839-2 ^ 
![{\displaystyle x(t)=\int _{-\infty }^{\infty }{g(\tau )\left[s(t-\tau )+n(t-\tau )\right]d\tau }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b083fe4e243e1e6d1906ff186fa8aaaf91bfae9e)
![{\displaystyle x[n]=\sum _{i=0}^{N}a_{i}w[n-i]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0b2f5cddce4bdf6ab2cad0bd57cc95d149b94708)
![{\displaystyle a_{i}=\arg \min ~E\{e^{2}[n]\}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2a927d41b2e72dfbf4330007de38dfccd4eeadae)
![{\displaystyle {\begin{array}{rcl}E\{e^{2}[n]\}&=&E\{(x[n]-s[n])^{2}\}\\&=&E\{x^{2}[n]\}+E\{s^{2}[n]\}-2E\{x[n]s[n]\}\\&=&E\{{\big (}\sum _{i=0}^{N}a_{i}w[n-i]{\big )}^{2}\}+E\{s^{2}[n]\}-2E\{\sum _{i=0}^{N}a_{i}w[n-i]s[n]\}.\end{array}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d6c92859cecfc7933906b59d6c7383e9656a6829)
![{\displaystyle [a_{0},\,\ldots ,\,a_{N}]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/23b15c3b4639ed827b254c60c133606cec814d33)
![{\displaystyle {\begin{array}{rcl}{\frac {\partial }{\partial a_{i}}}E\{e^{2}[n]\}&=&2E\{{\big (}\sum _{j=0}^{N}a_{j}w[n-j]{\big )}w[n-i]\}-2E\{s[n]w[n-i]\}\quad i=0,\,\ldots ,\,N\\&=&2\sum _{j=0}^{N}E\{w[n-j]w[n-i]\}a_{j}-2E\{w[n-i]s[n]\}.\end{array}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c049f8ef94a1de015893be949e6b25f362020a65)
![{\displaystyle w[n]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2a4e3e5afc2a8c6da9020b8c6b21450959101a18)
![{\displaystyle s[n]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ffe878c033cc70d1aa3ffbd156394791aa8dffde)
![{\displaystyle R_{w}[m]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6fd433a0ac1fe24b6b95d7bad92cea453be5c7f2)
![{\displaystyle R_{ws}[m]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7e7790e74796ce594e25ef54c65a4af6ea895803)
![{\displaystyle {\begin{aligned}R_{w}[m]=&E\{w[n]w[n+m]\}\\R_{ws}[m]=&E\{w[n]s[n+m]\}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b8cfb26024aaa8783a92fe728508710adf7496af)
![{\displaystyle R_{ws}[-i]\,=\,R_{sw}[i]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/29fb90f5d17ceba11cb91e4d4378d5b96c280eea)
![{\displaystyle {\frac {\partial }{\partial a_{i}}}E\{e^{2}[n]\}=2\sum _{j=0}^{N}R_{w}[j-i]a_{j}-2R_{sw}[i]\quad i=0,\,\ldots ,\,N}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3fa993e936efbab3a73f0e2fb25e0873fef01144)
![{\displaystyle \sum _{j=0}^{N}R_{w}[j-i]a_{j}=R_{sw}[i]\quad i=0,\,\ldots ,\,N}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0cbc8f4f8b82e431a3ab67e5be987b32e3fd8e7a)
![{\displaystyle {\begin{aligned}&\mathbf {T} \mathbf {a} =\mathbf {v} \\\Rightarrow &{\begin{bmatrix}R_{w}[0]&R_{w}[1]&\cdots &R_{w}[N]\\R_{w}[1]&R_{w}[0]&\cdots &R_{w}[N-1]\\\vdots &\vdots &\ddots &\vdots \\R_{w}[N]&R_{w}[N-1]&\cdots &R_{w}[0]\end{bmatrix}}{\begin{bmatrix}a_{0}\\a_{1}\\\vdots \\a_{N}\end{bmatrix}}={\begin{bmatrix}R_{ws}[0]\\R_{ws}[1]\\\vdots \\R_{ws}[N]\end{bmatrix}}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6baa68928b85556b809589f5d235d1b362bc0b6c)
