線性非時變系統理論 俗稱LTI 系統理論 ,源自應用數學 ,有在核磁共振頻譜學 、地震學 、電路 、信號處理 和控制理論 等技術領域的直接運用。該理論研究的是線性 的非時變系統 對任意輸入信號的響應。
此條目
需要補充更多來源 。
(2015年5月24日 )
LTI 系統通常僅關注系統在時間軸上的行為,但是類似的理論也可以擴展到空間維度。例如應用到圖像處理 和場論 時,系統的輸入也可在空間維度上變化。如果系統具有類似LTI 的行為,則這類系統也被稱為線性平移不變系統 。在離散(即採樣)系統中,對應的術語是線性時平移不變系統 。由電阻、電容、電感組成的電路是LTI系統的一個常見例子。[ 1]
顧名思義,線性非時變系統必須同時滿足線性 和非時變性 :
線性 ,指系統的輸入和輸出之間的關係是一個線性映射 :如果輸入
x
1
(
t
)
{\displaystyle x_{1}(t)\,}
產生響應
y
1
(
t
)
{\displaystyle y_{1}(t)\,}
,而輸入
x
2
(
t
)
{\displaystyle x_{2}(t)\,}
產生響應
y
2
(
t
)
{\displaystyle y_{2}(t)\,}
,那麼放縮 和加和 輸入
a
1
x
1
(
t
)
+
a
2
x
2
(
t
)
{\displaystyle a_{1}x_{1}(t)+a_{2}x_{2}(t)\,}
產生放縮、加和的響應
a
1
y
1
(
t
)
+
a
2
y
2
(
t
)
{\displaystyle a_{1}y_{1}(t)+a_{2}y_{2}(t)\,}
,其中
a
1
{\displaystyle a_{1}}
和
a
2
{\displaystyle a_{2}}
為實 純量 。此性質可以拓展到任意項,於是對於實數
c
1
,
c
2
,
…
,
c
k
{\displaystyle c_{1},c_{2},\ldots ,c_{k}}
,
輸入
∑
k
c
k
x
k
(
t
)
{\displaystyle \sum _{k}c_{k}\,x_{k}(t)}
產生輸出
∑
k
c
k
y
k
(
t
)
.
{\displaystyle \sum _{k}c_{k}\,y_{k}(t).\,}
特別地,
輸入
∫
−
∞
∞
c
ω
x
ω
(
t
)
d
ω
{\displaystyle \int _{-\infty }^{\infty }c_{\omega }\,x_{\omega }(t)\,\operatorname {d} \omega }
產生輸出
∫
−
∞
∞
c
ω
y
ω
(
t
)
d
ω
{\displaystyle \int _{-\infty }^{\infty }c_{\omega }\,y_{\omega }(t)\,\operatorname {d} \omega \,}
Eq.1
其中,
c
ω
{\displaystyle c_{\omega }}
和
x
ω
{\displaystyle x_{\omega }}
是純量,而輸入在序號為
ω
{\displaystyle \omega }
的連續統 內變化。因此,如果輸入函數可以由一個連續統的輸入函數像上面展示的那樣,「線性」組合而成,則對應的輸出函數,可以通過相應連續統的輸出函數以相同的方式縮放 和求和 得到。
非時變性 ,指如果將系統的輸入信號延遲
τ
{\displaystyle \tau }
秒而得到的輸出,除了這
τ
{\displaystyle \tau }
秒延時以外,是完全相同的,則稱這樣的系統是「非時變」的。即對於具有非時變性的系統,若系統輸入
x
(
t
)
{\displaystyle x(t)}
,對應的輸出為
y
(
t
)
{\displaystyle y(t)}
,則輸入為
x
(
t
+
τ
)
{\displaystyle x(t+\tau )}
時,系統的輸出為
y
(
t
+
τ
)
{\displaystyle y(t+\tau )}
。
LTI系統的理論的基本結論是任何LTI系統都可以完全用一個單一方程來表示,該方程稱為系統的衝激響應 。系統的輸出可以簡單表示為輸入信號與系統的衝激響應的卷積 。這種分析方法通常稱為時域 觀點。相同的結果對於離散時間線性移位不變系統也成立,其中信號為離散時間取樣信號,並且卷積對序列定義。
時域和頻域之間的關係
同理,任何LTI系統的特徵可由頻域 的系統傳遞函數 刻畫,它是系統衝激響應的拉普拉斯變換 (在離散時間系統的情況下為Z變換)。由於這些變換的性質,該系統在頻域的輸出是傳遞函數與輸入的變換的乘積。換句話說,時域中的卷積相當於頻域中的乘法。
對於所有的LTI系統中,本徵函數 和所用變換的基函數,是復 指數函數 。這即是說,如果一個系統的輸入是復波形
A
e
s
t
{\displaystyle Ae^{st}}
,復振幅為
A
{\displaystyle A}
,復頻率為
s
{\displaystyle s}
,輸出將是輸入的復常數倍,表示為新復振幅
B
{\displaystyle B}
的式子
B
e
s
t
{\displaystyle Be^{st}}
。比值
B
/
A
{\displaystyle B/A}
是頻率
s
{\displaystyle s}
的傳遞函數。
因為是正弦 的復指數與復共軛頻率的總和,如果輸入到該系統是一個正弦波,則系統的輸出也將是一個正弦波。輸出的正弦波可以具有不同振幅 和不同相位 ,但是輸出的穩態頻率一定與輸入相同。LTI系統不能產生頻率成分中沒有的輸入。
LTI系統理論適用於描述許多重要的系統。相對於時間變化的和/或非線性的系統,LTI系統通常是「容易」分析的。任何可以被模擬為常係數線性齊次微分方程 系統都是LTI系統。由電阻器 ,電感 和電容器 組成的電路 (RLC電路)是這類系統的一個常見實例。理想的彈簧 - 質量 - 阻尼系統也是LTI系統,並且在數學上與某個RLC電路等效。
LTI系統概念都是連續時間和離散時間(線性移位不變)的情況下相似。在圖像處理中,時間變量被替換為2空間變量,時間不變性的概念被替換為二維移不變性。當分析濾波器組s和MIMO系統中,常常是有用考慮的信號向量 。
線性系統不是非時變可以用其他方法來解決,如格林函數 方法。同樣的方法時,必須使用問題的初始條件是不為空。
輸入信號為x(t),輸出信號為y(t)的線性非時變系統的行為可以用卷積積分描述:[ 2]
y
(
t
)
=
x
(
t
)
∗
h
(
t
)
{\displaystyle y(t)=x(t)*h(t)\,}
=
d
e
f
∫
−
∞
∞
x
(
t
−
τ
)
⋅
h
(
τ
)
d
τ
{\displaystyle {}\quad {\stackrel {\mathrm {def} }{=}}\ \int _{-\infty }^{\infty }x(t-\tau )\cdot h(\tau )\,\operatorname {d} \tau }
=
∫
−
∞
∞
x
(
τ
)
⋅
h
(
t
−
τ
)
d
τ
,
{\displaystyle {}\quad =\int _{-\infty }^{\infty }x(\tau )\cdot h(t-\tau )\,\operatorname {d} \tau ,}
(使用交換律 )
其中
h
(
t
)
{\displaystyle \textstyle h(t)}
為當輸入信號
x
(
τ
)
=
δ
(
τ
)
{\displaystyle \textstyle x(\tau )=\delta (\tau )}
時系統的衝激 響應。因此
y
(
t
)
{\displaystyle \textstyle y(t)}
與輸入函數
x
(
τ
)
{\displaystyle \textstyle x(\tau )}
的加權平均成正比。權重函數為
h
(
−
τ
)
{\displaystyle \textstyle h(-\tau )}
,就是平移了
t
{\displaystyle \textstyle t}
的量。隨着
t
{\displaystyle \textstyle t}
改變,權重函數會突出輸入函數的不同部分。當對所有非負
τ
{\displaystyle \textstyle \tau }
,
h
(
τ
)
{\displaystyle \textstyle h(\tau )}
均為零時,
y
(
t
)
{\displaystyle \textstyle y(t)}
只由時間
t
{\displaystyle \textstyle t}
之前的
x
{\displaystyle \textstyle x}
值決定,而系統稱為因果系統 。
要理解為何LTI系統的輸出可以用卷積產生,就令記號
{
x
(
u
−
τ
)
;
u
}
{\displaystyle \textstyle \{x(u-\tau );\ u\}}
表示變量
u
{\displaystyle \textstyle u}
和常量
τ
{\displaystyle \textstyle \tau }
的函數
x
(
u
−
τ
)
{\displaystyle \textstyle x(u-\tau )}
。用簡潔的記號
{
x
}
{\displaystyle \textstyle \{x\}\,}
表示
{
x
(
u
)
;
u
}
{\displaystyle \textstyle \{x(u);\ u\}}
。那麼就會有一個從輸入函數
{
x
}
,
{\displaystyle \textstyle \{x\},}
轉換到
{
y
}
{\displaystyle \textstyle \{y\}}
的連續時間系統。在一般情況下,輸出的每一個值可以對應輸入的每一個值。這個概念表示為:
y
(
t
)
=
def
O
t
{
x
}
,
{\displaystyle y(t)\ {\stackrel {\text{def}}{=}}\ O_{t}\{x\},}
其中
O
t
{\displaystyle \textstyle O_{t}}
為對時間
t
{\displaystyle \textstyle t}
的變換算子。在典型的系統中,
y
(
t
)
{\displaystyle \textstyle y(t)}
很大程度上取決於
t
{\displaystyle \textstyle t}
臨近時間的
x
{\displaystyle \textstyle x}
的值。除非變換本身隨着
t
{\displaystyle \textstyle t}
變化,否則輸出函數就是常數,系統也沒有意義。
對一個線性系統,
O
{\displaystyle \textstyle O}
必須滿足Eq.1 :
O
t
{
∫
−
∞
∞
c
τ
x
τ
(
u
)
d
τ
;
u
}
=
∫
−
∞
∞
c
τ
y
τ
(
t
)
⏟
O
t
{
x
τ
}
d
τ
.
{\displaystyle O_{t}\left\{\int _{-\infty }^{\infty }c_{\tau }\ x_{\tau }(u)\,\operatorname {d} \tau ;\ u\right\}=\int _{-\infty }^{\infty }c_{\tau }\ \underbrace {y_{\tau }(t)} _{O_{t}\{x_{\tau }\}}\,\operatorname {d} \tau .\,}
Eq.2
而非時變系統的要求是:
O
t
{
x
(
u
−
τ
)
;
u
}
=
y
(
t
−
τ
)
=
def
O
t
−
τ
{
x
}
.
{\displaystyle {\begin{aligned}O_{t}\{x(u-\tau );\ u\}\ &{\stackrel {\quad }{=}}\ y(t-\tau )\\&{\stackrel {\text{def}}{=}}\ O_{t-\tau }\{x\}.\,\end{aligned}}}
Eq.3
在這種記號下,我們可以把衝激響應 寫成
h
(
t
)
=
def
O
t
{
δ
(
u
)
;
u
}
{\displaystyle \textstyle h(t)\ {\stackrel {\text{def}}{=}}\ O_{t}\{\delta (u);\ u\}}
。
同樣:
h
(
t
−
τ
)
{\displaystyle h(t-\tau )\,}
=
def
O
t
−
τ
{
δ
(
u
)
;
u
}
{\displaystyle {}{\stackrel {\text{def}}{=}}\ O_{t-\tau }\{\delta (u);\ u\}}
=
O
t
{
δ
(
u
−
τ
)
;
u
}
.
{\displaystyle {}=O_{t}\{\delta (u-\tau );\ u\}.\,}
(使用Eq.3 )
將此結果代入卷積積分:
x
(
t
)
∗
h
(
t
)
=
∫
−
∞
∞
x
(
τ
)
⋅
h
(
t
−
τ
)
d
τ
=
∫
−
∞
∞
x
(
τ
)
⋅
O
t
{
δ
(
u
−
τ
)
;
u
}
d
τ
,
{\displaystyle {\begin{aligned}x(t)*h(t)&=\int _{-\infty }^{\infty }x(\tau )\cdot h(t-\tau )\,\operatorname {d} \tau \\&=\int _{-\infty }^{\infty }x(\tau )\cdot O_{t}\{\delta (u-\tau );\ u\}\,\operatorname {d} \tau ,\,\end{aligned}}}
該形式為
c
τ
=
x
(
τ
)
{\displaystyle \textstyle c_{\tau }=x(\tau )}
且
x
τ
(
u
)
=
δ
(
u
−
τ
)
{\displaystyle \textstyle x_{\tau }(u)=\delta (u-\tau )}
情形下Eq.2 等式右側的形式。
那麼Eq.2 允許這個延拓:
x
(
t
)
∗
h
(
t
)
=
O
t
{
∫
−
∞
∞
x
(
τ
)
⋅
δ
(
u
−
τ
)
d
τ
;
u
}
=
O
t
{
x
(
u
)
;
u
}
=
def
y
(
t
)
.
{\displaystyle {\begin{aligned}x(t)*h(t)&=O_{t}\left\{\int _{-\infty }^{\infty }x(\tau )\cdot \delta (u-\tau )\,\operatorname {d} \tau ;\ u\right\}\\&=O_{t}\left\{x(u);\ u\right\}\\&\ {\stackrel {\text{def}}{=}}\ y(t).\,\end{aligned}}}
綜上所述,輸入函數
{
x
}
{\displaystyle \textstyle \{x\}}
可以用Eq.1 中描述的時移衝激函數的連續統的「線性」組合來表示。系統的線性特性允許系統由相應的以相同方式組合的衝激響應 的連續統來表示系統的響應。而非時變特性允許用卷積積分來表示這種組合。
上述數學運算可以用一個簡單的圖形模擬。[ 3]
本徵函數 是算子輸出為經過放縮的相同函數的函數。即,
H
f
=
λ
f
{\displaystyle {\mathcal {H}}f=\lambda f}
,
其中f 是本徵函數而
λ
{\displaystyle \lambda }
是特徵值 (一個常數)。
指數函數
A
e
s
t
{\displaystyle Ae^{st}}
(其中
A
,
s
∈
C
{\displaystyle A,s\in \mathbb {C} }
)是線性 非時變 算子的本徵函數 。可以用一個簡單的證明來說明這個概念。假設輸入是
x
(
t
)
=
A
e
s
t
{\displaystyle x(t)=Ae^{st}}
。系統衝激響應
h
(
t
)
{\displaystyle h(t)}
的輸出就是
∫
−
∞
∞
h
(
t
−
τ
)
A
e
s
τ
d
τ
{\displaystyle \int _{-\infty }^{\infty }h(t-\tau )Ae^{s\tau }\,\operatorname {d} \tau }
由卷積 的交換性質,上式等價於
∫
−
∞
∞
h
(
τ
)
A
e
s
(
t
−
τ
)
d
τ
⏞
H
f
=
∫
−
∞
∞
h
(
τ
)
A
e
s
t
e
−
s
τ
d
τ
=
A
e
s
t
∫
−
∞
∞
h
(
τ
)
e
−
s
τ
d
τ
=
A
e
s
t
⏟
Input
⏞
f
H
(
s
)
⏟
Scalar
⏞
λ
,
{\displaystyle {\begin{aligned}\overbrace {\int _{-\infty }^{\infty }h(\tau )\,Ae^{s(t-\tau )}\,\operatorname {d} \tau } ^{{\mathcal {H}}f}&=\int _{-\infty }^{\infty }h(\tau )\,Ae^{st}e^{-s\tau }\,\operatorname {d} \tau &=Ae^{st}\int _{-\infty }^{\infty }h(\tau )\,e^{-s\tau }\,\operatorname {d} \tau \\&=\overbrace {\underbrace {Ae^{st}} _{\text{Input}}} ^{f}\overbrace {\underbrace {H(s)} _{\text{Scalar}}} ^{\lambda },\end{aligned}}}
其中純量
H
(
s
)
=
def
∫
−
∞
∞
h
(
t
)
e
−
s
t
d
t
{\displaystyle H(s)\ {\stackrel {\text{def}}{=}}\ \int _{-\infty }^{\infty }h(t)e^{-st}\,\operatorname {d} t}
只與參數s 有關。
因此,系統的響應是一個縮放的輸入。特別地,對任意I
A
,
s
∈
C
{\displaystyle A,s\in \mathbb {C} }
,系統輸出為輸入
A
e
s
t
{\displaystyle Ae^{st}}
和常量
H
(
s
)
{\displaystyle H(s)}
的乘積。因此,
A
e
s
t
{\displaystyle Ae^{st}}
是LTI系統的本徵函數 ,對應的特徵向量 為
H
(
s
)
{\displaystyle H(s)}
。
因果性和穩定性是描述系統的兩個重要性質。如果獨立變量是時間,那麼因果性是必須的,但並不是所有系統的獨立變量都是時間。例如,一個處理靜止圖像的系統不需要具備因果性。非因果系統可以建立,並可以在許多情況下發揮作用。即使是非實數 系統也可以構建,並且在很多場合也是非常有用的。
如果系統輸出只與當前以及過去的輸入有關,那麼該系統就是因果系統。因果性的充分必要條件是
h
(
t
)
=
0
∀
t
<
0
,
{\displaystyle h(t)=0\quad \forall t<0,}
其中
h
(
t
)
{\displaystyle h(t)}
是衝激響應。由於拉普拉斯變換的逆變換不唯一,所以通常不能根據拉普拉斯變換確定系統的因果性。只有在確定了系統的收斂域 之後才能確定該系統的因果性。
幾乎所有的連續時間系統都能找到與之對應的離散時間系統。
在許多情況下,離散時間(DT)系統實際上是較大的連續時間(CT)系統的一部分。例如,數字錄音系統記錄模擬聲音、數碼化、或許對數碼訊號進行處理、然後重放模擬信號。
正式場合下所研究的離散時間信號幾乎總是連續時間信號的均勻採樣。如果
x
(
t
)
{\displaystyle x(t)}
是一個連續時間信號,那麼模數轉換器 將把它轉換成離散時間信號
x
[
n
]
{\displaystyle x[n]}
,
x
[
n
]
=
x
(
n
T
)
{\displaystyle x[n]=x(nT)}
,
其中T 是採樣周期 。為了保證離散信號能夠忠實地表示輸入信號,非常重要的一點就是需要限制輸入信號的頻率範圍。根據採樣定理 ,離散時間信號所包括的最大頻率範圍是
1
/
(
2
T
)
{\displaystyle 1/(2T)}
。其它頻率都成為這個範圍的混疊 信號。
我們從一個衝激響應是二維函數的時變系統開始來看看非時變這個條件是如何將系統降到一維的。例如,假設輸入信號是
x
[
n
]
{\displaystyle x[n]}
,其中n是整數,即
n
∈
Z
{\displaystyle n\in \mathbb {Z} }
。線性算子
H
{\displaystyle {\mathcal {H}}}
表示系統在輸入信號上的操作,對於這個index set來說合適的算子是一個二維函數
h
[
n
1
,
n
2
]
where
n
1
,
n
2
∈
Z
{\displaystyle h[n_{1},n_{2}]{\mbox{ where }}n_{1},n_{2}\in \mathbb {Z} }
。
由於
H
{\displaystyle {\mathcal {H}}}
是一個線性算子,系統在輸入信號
x
[
n
]
{\displaystyle x[n]}
上的作用就是下面累加 和所表示的線性變換
y
[
n
1
]
=
∑
n
2
=
−
∞
∞
h
[
n
1
,
n
2
]
x
[
n
2
]
,
{\displaystyle y[n_{1}]=\sum _{n_{2}=-\infty }^{\infty }h[n_{1},n_{2}]\,x[n_{2}],}
如果線性算子
H
{\displaystyle {\mathcal {H}}}
也是非時變 的,那麼
h
[
n
1
,
n
2
]
=
h
[
n
1
+
m
,
n
2
+
m
]
∀
m
∈
Z
{\displaystyle h[n_{1},n_{2}]=h[n_{1}+m,n_{2}+m]\qquad \forall \,m\in \mathbb {Z} }
。
如果取
m
=
−
n
2
,
{\displaystyle m=-n_{2},\,}
那麼
h
[
n
1
,
n
2
]
=
h
[
n
1
−
n
2
,
0
]
.
{\displaystyle h[n_{1},n_{2}]=h[n_{1}-n_{2},0].\,}
為了簡化通常我們丟棄
h
[
n
1
,
n
2
]
{\displaystyle h[n_{1},n_{2}]}
的第二個參數零,這樣重疊積分現在變成了濾波中常見的卷積 和
y
[
n
1
]
=
∑
n
2
=
−
∞
∞
h
[
n
1
−
n
2
]
x
[
n
2
]
=
(
h
∗
x
)
[
n
1
]
{\displaystyle y[n_{1}]=\sum _{n_{2}=-\infty }^{\infty }h[n_{1}-n_{2}]\,x[n_{2}]=(h*x)[n_{1}]}
。
這樣,卷積 和表示一個線性 非時變系統 在任意輸入函數上所起的作用,對於類似的有限維參數,參見輪換矩陣
因果性和穩定性是系統的重要特性。與連續時間系統不同,我們可以實現非因果的離散時間系統。通過在系統中加入延時就很容易將非因果有限衝激響應 系統變成因果系統。甚至可以構建非因果的無限衝激響應 系統(參見Vaidyanathan and Chen, 1995)。我們也可以構建不穩定的系統,這種系統在很多場合都很有用,甚至也可以構建在很多情況下非常有用的non-real 系統。
如果系統的輸出只與當前以及過去的輸入有關,那麼系統就是因果系統。因果性的必要且充分條件是
h
[
n
]
=
0
∀
n
<
0
,
{\displaystyle h[n]=0\ \forall n<0,}
其中
h
[
n
]
{\displaystyle h[n]}
是衝激響應。由於逆變換不是唯一的,所以通常很難從Z變換確定系統的因果性。如果收斂域 確定,系統的因果性也就隨之確定。