線性標準轉換

维基百科,自由的百科全书
跳转至: 导航搜索

在數學的文獻中,線性標準轉換(linear canonical transform, LCT)也稱作"ABCD轉換。在漢米爾頓力學中,線性標準轉換是積分變換的一個代表家族,並且能夠將許多經典的轉換進行廣義化,例如傅立葉變換分數傅立葉變換拉普拉斯變換菲涅爾轉換(Fresnel transform,電磁波在空氣中傳播)、高斯-魏爾斯特拉斯轉換包格曼轉換等等等。此轉換提供了這些最常使用的線性轉換一個統一框架,並且在光學、信號轉換以及系統響應領域中都提供一般化的概念。尤其從系統工程的角度看來,線性標準轉換提供一個強大的光學系統設計和分析的工具。 此轉換有四維變數的線性積分轉換和一個限制條件,因此實際上是一個三維自由度的積分變換的家族。 在群論中,線性標準轉換屬於特殊線性群(SR(2))在時頻域上的一個作用群。


定義[编辑]

從積分變換的定義開始:

 (Tf)(u) = \int \limits_{t_1}^{t_2} K(t, u)\, f(t)\, dt

積分變換是將輸入的訊號f(t)經由核心(Kernel)K(t, u)的作用後對應後的輸出結果即是另一個函數(TF)(u),此時(TF)(u)稱作是f(t)的積分變換。

線性標準轉換(LCT)以一般的線性積分變換關係表示如下:

X_{(a,b,c,d)}(u) = \sqrt{-i} \cdot e^{i  \pi \frac{d}{b} u^{2}} \int_{-\infty}^\infty e^{-i 2 \pi \frac{1}{b} ut}e^{i  \pi \frac{a}{b} t^2} x(t) \; dt \, ,   此時 b ≠ 0,
X_{(a,0,c,d)}(u) = \sqrt{d} \cdot e^{i  \pi cdu^{2}} x(du) \, , 此時 b = 0.

並且需要滿足此條件:

ad-bc=1

此時的X_{(a,b,c,d)}(u)表示的是訊號x(t)經過線性標準變換過後的結果 而我們通常使用O_{F}^{(a,b,c,d)}作為線性標準轉換的操作函數,也就是:

O_{F}^{(a,b,c,d)}(x(t) = X_{(a,b,c,d)}(u))

另外,線性標準轉換(LCT)也可用簡單的2x2的矩陣以及一個行列式限制條件來表示,如同矩陣\begin{bmatrix}
a & b \\
c & d
\end{bmatrix}
,其條件為ad - bc = 1

矩陣形式在時頻分析中所代表的意義是將時頻域上的平行四邊形面積扭曲成另一個平性四邊形面積,而不同的特殊情形的矩陣形式則分別代表不同的幾何轉換


特殊情形[编辑]

線性標準轉換是許多經典轉換的廣義化。透過將矩陣形式帶入不同的變數所得到的特殊情形如下,以下舉出的轉換在時頻分析上均有訊號在時頻域運動的意義。例如:

O^{(cos\phi ,sin\phi ,-sin\phi ,cos\phi )}_F \left\{ f(t) \right\} =  \sqrt{e^{-j\phi}} \cdot O^{\phi }_F \left\{ f(t) \right \}
若寫成矩陣的形式為:
\begin{bmatrix} a & c \\ b & d \end{bmatrix} = \begin{bmatrix} \cos\phi & \sin\phi \\ -\sin\phi & \cos\phi \end{bmatrix}
F(u) = FT{(f(t))} = \int_{-\infty}^{\infty} f(t)\ e^{- j u t}\,dx,   對任意實數 u
可以將LCT整理成為以下形式:
O^{(0,1,-1,0)}_F \left\{ f(t) \right\} = \sqrt{ \frac{-j}{2\pi}} \cdot \int_{-\infty}^{\infty}e^{-jut}\,f(t)\,dt=\sqrt{-j} \cdot FT \left\{ f(t) \right \}
此時即是分數傅立葉變換中,當\phi = \pi /2的特殊情形,代表將時頻域上的面積做\phi = +\pi /2度的旋轉。另外寫成矩陣的形式則為:
\begin{bmatrix} a & b \\ c & d \end{bmatrix} = \begin{bmatrix} 0 & 1 \\ -1 & 0 \end{bmatrix}
  • 傅立葉反變換,傅立葉反變換的定義:
f(t) = IFT{(F(u))} = \int_{-\infty}^{\infty} F(u)\ e^{- j u t}\,dx,   對任意實數 u
可以將LCT整理成為以下形式:
O^{(0,-1,1,0)}_F \left\{ f(t) \right\} = \sqrt{ \frac{j}{2\pi}} \cdot \int_{-\infty}^{\infty}e^{jut}\,f(t)\,dt=\sqrt{j} \cdot FT^{-1} \left\{ f(t) \right \}
此時即是分數傅立葉變換中,當\phi = -\pi /2的特殊情形,也就是代表將時頻域上的面積做\phi = -\pi /2度的旋轉。另外寫成矩陣的形式則為:
\begin{bmatrix} a & b \\ c & d \end{bmatrix} = \begin{bmatrix} 0 & -1 \\ 1 & 0 \end{bmatrix}
  • 菲涅爾轉換(Fresnel transform),用來描述電磁波在空氣中傳播的情形,代表將時頻域上的面積做時間方向上的切變,若將其寫成矩陣的形式為:
\begin{bmatrix} a & b \\ c & d \end{bmatrix} = \begin{bmatrix} 1 & \lambda z/2\pi \\ 0 & 1 \end{bmatrix}
  • 拉普拉斯轉換(Laplace transform),代表將時頻域上的面積對其複數平面做+90度的旋轉,其矩陣如下:
\begin{bmatrix}  a & b \\  c & d \end{bmatrix} = \begin{bmatrix}  0 & i \\  i & 0 \end{bmatrix}
  • 分數拉普拉斯轉換,代表將時頻域上的面積對其複數平面做任意角度的旋轉,其矩陣如下:
\begin{bmatrix}  a & b \\  c & d\end{bmatrix} = \begin{bmatrix}  i \cos \phi& i \sin \phi\\  i \sin \phi& -i \cos \phi\end{bmatrix}
  • 啾啾聲載波(chirp multiplication),其定義是如下:
X_{a,0,c,d}(u)=e^{j\pi \tau u^2}x(u)
寫成矩陣的形式為:
\begin{bmatrix} a & c \\ b & d \end{bmatrix} = \begin{bmatrix} 1 & 0 \\ \tau & 1 \end{bmatrix}

特性[编辑]

  • 加成性:
當我們將兩個參數矩陣不同的LCT組合起來,也就是對同一個訊號做兩次線性標準轉換。這樣的結果相當於將兩矩陣的乘積後的矩陣作為參數的線性標準轉換,並以此新的LCT對原始訊號作用後得到的結果。最常見是維格納分佈(wigner distribution function, WDF)中的加法特性。
當使用O_{F}^{(a,b,c,d)}來表示LCT的作用時:
O_{F}^{(a,b,c,d)}(x(t) = X_{(a,b,c,d)}(u)
若現在有兩種轉換O_{F}^{(a_1,b_1,c_1,d_1)}(x(t) = X_{(a_1,b_1,c_1,d_1)}(u))O_{F}^{(a_2,b_2,c_2,d_2)}(x(t) = X_{(a_2,b_2,c_2,d_2)}(u))兩者依序對x(t)作用
則有第三種LCT其參數為前兩者的矩陣形式的乘積O_{F}^{(a_3,b_3,c_3,d_3)}(x(t) = X_{(a_3,b_3,c_3,d_3)}(u)),將滿足:
O_{F}^{(a_2,b_2,c_2,d_2)}\left\{O_{F}^{(a_1,b_1,c_1,d_1)}[x(t)]\right\}=O_{F}^{(a_3,b_3,c_3,d_3)}[x(t)]
用矩陣來表示為:
\begin{bmatrix} a_3 & b_3 \\ c_3 & d_3 \end{bmatrix} = \begin{bmatrix} a_2 & b_2 \\ c_2 & d_2 \end{bmatrix} \begin{bmatrix} a_1 & b_1 \\ c_1 & d_1 \end{bmatrix}
因此,若是W(u,v)X(u)做WDF,而X(u)x(t)做LCT,可以表示成:
W_{X_{(a,b,c,d)}}(u,v)=W_{x}(du-bv,-cu+av)
W_{X_{(a,b,c,d)}}(au+bv,cu+dv)=W_{x}(u,v)
此兩式互為相同含意,因此可知反向LCT的就是由四個參數代表矩陣的反矩陣
  • 結合律:
當現在有三種參數不同的LCT依序對x(t)作用,則前兩者先作用第三種再作用的結果會與,後二者先作用再由第一種作用的結果相同。

LCT的訊號扭曲運用[编辑]

利用LCT的特性,可以把原本的訊號轉變成另外的形式,此過程可以輕易的用LCT來將一個矩形的區域,變成另外一個平行四邊形區域,此概念是利用線性代數上的基底變換,來把原本的基底換成另外一種型式來呈現,因此要符合線性代數中基底變換的限制,其中心點為(0,0),若無法滿足此條件,則無法使用LCT的線性轉換公式,並且此轉換的行列式為1,不會改變原來區域的面積大小,並且轉換前後皆為平行四邊形,因此可以把一些不好運算的平行四邊形,藉由LCT的方式,轉換到比較容易計算或是比較容易理解的基底,有利於運算的考量。 舉例來說:

在二維的平面上,若已知兩個平面個別為,一個中心點為原點的矩形,四個點分別為:(-1,2),(1,2),(-1,-2),(1,-2),和經過LCT的轉換可以變成中心點為原點的平行四邊形,其四個點為:(1,1),(7,3),(-7,-3),(-1,-1),可以藉由其中兩個點來解聯立方程式,算出a,b,c,d,所代表的數值,

\begin{cases} a(-1)+b(2)=1,\ a(1)+b(2)=7 \\ c(-1)+d(2)=1,\ c(1)+d(2)=3 \end{cases}
\Rightarrow\begin{bmatrix} a & c \\ b & d \end{bmatrix}=\begin{bmatrix} 3 & 1 \\ 2 & 1 \end{bmatrix}

其總面積皆為:8

LCT在各領域的應用[编辑]

LCT在許多的工程領域中都十分地實用,因為LCT能夠用來描述任何的二次相位系統,以及這些二次相位系統的串接組合,以下我們將舉出LCT在多個領域中的應用實例:

  • 電磁波傳播:
電磁波在空氣中傳播的行為是透過菲涅爾轉換(Fresnel transform)來描述,假設一光學系統滿足菲涅爾近似而電磁波由x_i, y_i平面傳播到x, y平面,則式子如下:
 U_0(x,y) = - \frac{j}{\lambda} \frac{e^{jkz}}{z} \int_{-\infty}^\infty \int_{-\infty}^{\infty} e^{j \frac{k}{2z}  [ (x-x_i)^2 +(y-y_i)^2 ] } U_i(x_i,y_i) \; dx_i\; dy_i,
其中
k = 2 \pi / \lambda    : 波數;
\lambda  : 波長;
z  : 傳播距離;
j  : 虛數.
如果我們以LCT的矩陣形式來表示的話則是:
\begin{bmatrix} a & b \\ c & d \end{bmatrix} = \begin{bmatrix} 1 & \lambda z/2\pi \\ 0 & 1 \end{bmatrix}
在時頻域上將訊號作對時間軸的切變,其中\lambda越大,則切變效果越劇烈。
  • 光學系統
  • 球面薄凸透鏡
  • 假設透鏡的折射係數為n,則電磁波經過薄透鏡折射之後的結果,可以用下列式子表示:
 U_0(x,y) = e^{jkn \Delta} e^{-j \frac{k}{2f} [x^2 + y ^2]} U_i(x,y)
  • 如果以LCT的矩陣形式來表示則是:
 \begin{bmatrix}  a & b \\  c & d \end{bmatrix}= \begin{bmatrix}  1 & 0 \\  \frac{-1}{\lambda f} & 1 \end{bmatrix}
在時頻域上將訊號作對頻率軸的切變,但此時與菲涅爾轉換不同的是\lambda越大,切變得效果越小。
  • 球面反射鏡
球面反射鏡的LCT矩陣形式與薄凸透鏡很類似,但此時影響頻率軸上的切變不是焦距而是求面反射鏡的半徑,其矩陣形式如下:
 \begin{bmatrix}  a & b \\  c & d \end{bmatrix}= \begin{bmatrix}  1 & 0 \\  \frac{-1}{\lambda R} & 1 \end{bmatrix}
  • 結合LCT加成性的例子
當電磁波由左端平面輸入,在空氣傳播時使用菲涅爾轉換描述,接著在距離為z1處經過一焦距為f的薄凸透鏡折射後再射入空氣中傳播距離z2後由右端平面輸出,此過程經過三層LCT的操作,依序為菲涅爾轉換、薄凸透鏡折射和菲涅爾轉換。
根據LCT加成性,整個LCT過程相當於作用以下三個矩陣相乘後的矩陣乘積。
\begin{bmatrix} a & b \\ c & d \end{bmatrix} = \begin{bmatrix} 1 & \lambda z_1/2\pi \\ 0 & 1 \end{bmatrix}\begin{bmatrix}  1 & 0 \\  \frac{-1}{\lambda f} & 1 \end{bmatrix} \begin{bmatrix} 1 & \lambda z_/2\pi \\ 0 & 1 \end{bmatrix}

參考[编辑]

Jian-Jiun Ding, Time frequency analysis and wavelet transform class note,the Department of Electrical Engineering, National Taiwan University (NTU), Taipei, Taiwan, 2009.

A. Stern, “Why is the linear canonical transform so little known?” in Proceedings of 5th International Workshop on Information Optics, G. Cristóbal, B. Javidi, and S. Vallmitjana, eds (Springer, 2006), pp. 225–234.

许天周,李炳照,线性正则变换及应用,北京:科学出版社,2013.02