線性非移變系統理論

線性非移變系統理論建立在圖訊號處理中對線性非移變系統的定義上，致力於承襲線性非時變系統理論在操作上的便利性。 對於一個N點的圖 G {\displaystyle G} （可為有向，但通常有限、不重邊），可定義圖位移 S {\displaystyle S} ，為一線性映射從圖訊號時域映射到圖訊號時域，可表為一N階方陣 。 圖位移是一個抽象定義，並沒有特別指對 G {\displaystyle G} 使用哪種特定方法構造出來的為圖位移。 比較被使用的圖位移有：連接矩陣A、拉普拉斯矩陣L、正規化拉普拉斯矩陣LN。 線性非移變系統是圖訊號處理中，有最簡單結構的系統。其性質可最大限度地和數位訊號處理中的線性非時變系統對照，並發展出與之類似的運算，達到想要的訊號處理效果。 在一 N {\displaystyle N} 點的圖上，給定已定義的圖位移 S {\displaystyle S} ，若線性映射 H : C N → C N {\displaystyle H:\mathbb {C} ^{N}\rightarrow \mathbb {C} ^{N}} 滿足 H S = S H {\displaystyle HS=SH} ，則稱其為線性非移變系統（LSI）。 不受圖位移 S {\displaystyle S} 的作用前後順序影響。 可表示為圖圖位移 S {\displaystyle S} 的多項式，也就是說存在 h {\displaystyle h} 多項式，使得 H = h ( S ) {\displaystyle H=h(S)} 。更嚴謹地，若考慮的是N-1階以下的多項式，則 h {\displaystyle h} 唯一存在。 根據性質2.，再以給定好圖位移的情況下，對所有線性非移變系統 H {\displaystyle H} ，可唯一定義其傳遞函數 h {\displaystyle h} ， h {\displaystyle h} 為最高次項小於N-1階的多項式，且 H = h ( S ) {\displaystyle H=h(S)} 。 詳見圖論傅立葉轉換。 摺積在圖訊號時域中因為沒有顯著的前後關係，無法定義一個良好的計算式。但因同一張圖上的圖訊號在頻域中之頻率值相同，摺積的推廣定義可使用摺積原本的性質：頻域中的峰值相乘，也就是 s ^ o u t ( λ k ) = s ^ 1 ( λ k ) s ^ 2 ( λ k ) {\displaystyle {\hat {s}}_{out}(\lambda _{k})={\hat {s}}_{1}(\lambda _{k}){\hat {s}}_{2}(\lambda _{k})} 其中 s 1 {\displaystyle s_{1}} 、 s 2 {\displaystyle s_{2}} 為輸入訊號， s o u t {\displaystyle s_{out}} 為輸出訊號， { λ k } k = 0... N − 1 {\displaystyle \{\lambda _{k}\}_{k=0...N-1}} 為頻域之特徵頻率。對上式作逆圖論傅立業轉換，得 s o u t ( i ) = ( s 1 ∗ s 2 ) ( i ) = ∑ k = 0 N − 1 s ^ 1 ( λ k ) s ^ 2 ( λ k ) u k ( i ) {\displaystyle s_{out}(i)=(s_{1}*s_{2})(i)=\sum _{k=0}^{N-1}{\hat {s}}_{1}(\lambda _{k}){\hat {s}}_{2}(\lambda _{k})u_{k}(i)} 其中 { u k } k = 0... N − 1 {\displaystyle \{u_{k}\}_{k=0...N-1}} 為圖論傅立葉基底。將此式作為摺積之推廣定義。 上式有另一個使用矩陣的表示法 s o u t = V D s ^ 2 V − 1 s 1 {\displaystyle s_{out}=VD_{{\hat {s}}_{2}}V^{-1}s_{1}} 其中 V − 1 {\displaystyle V^{-1}} 為圖論傅立葉轉換矩陣， V {\displaystyle V} 即為逆傅立葉轉換矩陣， D s ^ 2 {\displaystyle D_{{\hat {s}}_{2}}} 表示一對角矩陣，其主對角線的元素為 s ^ 2 {\displaystyle {\hat {s}}_{2}} 平移在圖訊號時域中亦因為沒有好的順序關係，無法被良好的直接定義。 但因一般數位訊號處理的平移運算等價於將訊號與一個單位脈衝訊號摺積，其推廣定義可以利用此性質 s o u t = s i n ∗ δ n = V D s ^ i n V − 1 δ n = V D s ^ i n u n {\displaystyle s_{out}=s_{in}*\delta _{n}=VD_{{\hat {s}}_{in}}V^{-1}\delta _{n}=VD_{{\hat {s}}_{in}}u_{n}} ，其中 δ n ( k ) = { 1 , if  n = k 0 , if  n ≠ k {\displaystyle \delta _{n}(k)={\begin{cases}1,&{\text{if }}n=k\\0,&{\text{if }}n\neq k\end{cases}}} 濾波是讓輸入訊號經過一系統，產生輸出訊號的過程。在時域中的表現形式為 s o u t = H s i n {\displaystyle s_{out}=Hs_{in}} 若 H {\displaystyle H} 是線性非移變系統，存在多項式傳遞函數 h {\displaystyle h} 使得 h = h ( S ) {\displaystyle h=h(S)} ，故 H s i n = h ( S ) s i n = h ( V Λ V − 1 ) s i n = V h ( Λ ) V − 1 s i n = V D s ^ H V − 1 s i n {\displaystyle Hs_{in}=h(S)s_{in}=h(V\Lambda V^{-1})s_{in}=Vh(\Lambda )V^{-1}s_{in}=VD_{{\hat {s}}_{H}}V^{-1}s_{in}} 與摺積的定義比較，可知其實可將「對 s i n {\displaystyle s_{in}} 使用 H {\displaystyle H} 濾波」看作「將 s i n {\displaystyle s_{in}} 與另一個信號 s H {\displaystyle s_{H}} 作摺積」，其中 s H {\displaystyle s_{H}} 符合 s ^ H = ( h ( λ 1 ) , h ( λ 2 ) , . . . h ( λ N − 1 ) ) T {\displaystyle {\hat {s}}_{H}=(h(\lambda _{1}),h(\lambda _{2}),...h(\lambda _{N-1}))^{T}} ， 此時稱 s H {\displaystyle s_{H}} 為 H {\displaystyle H} 的脈衝響應， s ^ H {\displaystyle {\hat {s}}_{H}} 為 H {\displaystyle H} 的頻率響應。 A. Sandryhaila and Jose M. F. Moura, Discrete Signal Processing on Graphs, [[arxiv:1210.4752|http://arxiv.org/abs/1210.4752 （页面存档备份，存于互联网档案馆）]] David I Shuman, Sunil K. Narang, Pascal Frossard, Antonio Ortega, Pierre Vandergheynst, The Emerging Field of Signal Processing on Graphs: Extending High-Dimensional Data Analysis to Networks and Other Irregular Domains, [[arxiv:1211.0053|http://arxiv.org/abs/1211.0053 （页面存档备份，存于互联网档案馆）]] Sandryhaila, Aliaksei; Moura, Jose M. F. Discrete Signal Processing on Graphs. IEEE Transactions on Signal Processing: 1644–1656. [2016-06-30]. ISSN 1053-587X. doi:10.1109/TSP.2013.2238935. （原始内容存档于2018-12-12）.  Shuman, David I.; Narang, Sunil K.; Frossard, Pascal; Ortega, Antonio; Vandergheynst, Pierre. The Emerging Field of Signal Processing on Graphs: Extending High-Dimensional Data Analysis to Networks and Other Irregular Domains. IEEE Signal Processing Magazine: 83–98. [2016-06-30]. ISSN 1053-5888. doi:10.1109/MSP.2012.2235192. （原始内容存档于2020-01-11）.