線性標準轉換

在數學的，線性標準轉換(Linear Canonical Transform, LCT)也稱作線性正則變換、ABCD轉換、廣義Fresnel變換等。在漢米爾頓力學中，線性標準轉換是積分變換的一個代表家族，並且能夠將許多經典的轉換進行廣義化，例如傅立葉變換、分數傅立葉變換、拉普拉斯變換、菲涅爾轉換(電磁波在空氣中傳播)、、等等。此轉換提供了這些最常使用的線性轉換一個統一框架，並且在光學、訊號轉換以及系統響應領域中都提供一般化的概念。尤其從系統工程的角度看來，線性標準轉換提供一個強大的光學系統設計和分析的工具。此轉換有四維變數的線性積分轉換和一個限制條件，因此實際上是一個三維自由度的積分變換的家族。在群論中，線性標準轉換屬於特殊線性群(SR(2))在時頻域上的一個作用群。定義.從積分變換的定義開始:formula_1積分變換是將輸入的訊號f(t)經由核心(Kernel)K(t, u)的作用後對應後的輸出結果即是另一個函式formula_2，此時formula_2稱作是formula_4的積分變換。線性標準轉換(LCT)以一般的線性積分變換關係表示如下:並且需要滿足此條件:formula_5此時的formula_6表示的是訊號formula_7經過線性標準變換過後的結果而我們通常使用formula_8作為線性標準轉換的操作函式，也就是:formula_9另外，線性標準轉換(LCT)也可用簡單的2x2的矩陣以及一個行列式限制條件來表示，如同矩陣formula_10，其條件為formula_11矩陣形式在時頻分析中所代表的意義是將時頻域上的平行四邊形面積扭曲成另一個平性四邊形面積，而不同的特殊情形的矩陣形式則分別代表不同的幾何轉換特殊情形.線性標準轉換是許多經典轉換的廣義化。透過將矩陣形式帶入不同的變數所得到的特殊情形如下，以下舉出的轉換在時頻分析上均有訊號在時頻域運動的意義。例如：縮放.縮放, formula_12，是一種在時間軸和頻率軸進行縮小或放大操作的轉換，其中時間軸和頻率軸的操作是相反的，若放大時間軸就會縮小頻率軸，反之亦然：formula_13傅立葉轉換.傅立葉變換，傅立葉變換的定義:formula_14, 對任意實數 u可以將 LCT 整理成為以下形式:formula_15此時即是分數傅立葉變換中，當formula_16的特殊情形，代表將時頻域上的面積做formula_17度的旋轉。另外寫成矩陣的形式則為：formula_18分數傅立葉變換.分數傅立葉變換，代表將時頻域上的面積做formula_19度的旋轉，其積分變換形式為:formula_20若寫成矩陣的形式為：formula_21傅立葉反變換.傅立葉反變換，傅立葉反變換的定義:formula_22,   對任意實數 u可以將LCT整理成為以下形式:formula_23此時即是分數傅立葉變換中，當formula_24的特殊情形，也就是代表將時頻域上的面積做formula_24度的旋轉。另外寫成矩陣的形式則為：formula_26菲涅爾轉換.菲涅爾轉換，用來描述電磁波在空氣中傳播的情形，代表將時頻域上的面積做時間方向上的切變，若將其寫成矩陣的形式為：formula_27其中 formula_28 是距離，formula_29 是波長。拉普拉斯轉換.拉普拉斯轉換，代表將時頻域上的面積對其複數平面做+90度的旋轉，其矩陣如下:formula_30分數拉普拉斯轉換.分數拉普拉斯轉換，代表將時頻域上的面積對其複數平面做任意角度的旋轉，其矩陣如下:formula_31其中拉普拉斯轉換是分數拉普拉斯轉換中 formula_32 的特例，而拉普拉斯反轉換則是 formula_32 的情形。啾啾聲載波.啾啾聲載波(chirp multiplication)，其定義是如下:formula_34寫成矩陣的形式為：formula_35當我們將兩個引數矩陣不同的LCT組合起來，也就是對同一個訊號做兩次線性標準轉換。這樣的結果相當於將兩矩陣的乘積後的矩陣作為引數的線性標準轉換，並以此新的LCT對原始訊號作用後得到的結果。最常見是維格納分佈(wigner distribution function, WDF)中的加法特性。當使用formula_8來表示LCT的作用時:formula_9若現在有兩種轉換formula_38和formula_39兩者依序對formula_7作用則有第三種LCT其引數為前兩者的矩陣形式的乘積formula_41，將滿足:formula_42用矩陣來表示為：formula_43因此，若是formula_44是formula_45做WDF，而formula_45是formula_7做LCT，可以表示成：formula_48formula_49此兩式互為相同含意，因此可知反向LCT的就是由四個引數代表矩陣的反矩陣。當現在有三種引數不同的LCT依序對formula_7作用，則前兩者先作用第三種再作用的結果會與，後二者先作用再由第一種作用的結果相同。反線性正則轉換.引數為formula_51的線性標準轉換其反線性標準轉換引數為formula_52formula_53LCT的訊號扭曲運用.利用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，所代表的數值，formula_54formula_55其總面積皆為：8LCT在各領域的應用.LCT在許多的工程領域中都十分地實用，因為LCT能夠用來描述任何的二次相位系統，以及這些二次相位系統的串接組合，以下我們將舉出LCT在多個領域中的應用例項:對一訊號做線性正則轉換相當於對時間-頻率域的分佈做扭曲一訊號  formula_56  及其線性正則變換 formula_57 的維格納分佈分別表示為 formula_58 和 formula_59 ，它們在時頻域的分佈有以下關係:formula_60formula_61利用線性正則轉換改變訊號在時頻域的分佈來設計濾波器。formula_62其中 formula_63 表示轉換後的域的濾波器電磁波在空氣中傳播的行為是透過菲涅爾轉換來描述，假設一光學系統滿足菲涅爾近似而電磁波由formula_64平面傳播到formula_65平面，則式子如下:formula_66其中如果我們以LCT的矩陣形式來表示的話則是:formula_67在時頻域上將訊號作對時間軸的切變，其中formula_29越大，則切變效果越劇烈。* 假設透鏡的折射係數為n，則電磁波經過薄透鏡折射之後的結果，可以用下列式子表示:formula_69* 如果以LCT的矩陣形式來表示則是:formula_70在時頻域上將訊號作對頻率軸的切變，但此時與菲涅爾轉換不同的是formula_29越大，切變得效果越小。球面反射鏡的LCT矩陣形式與薄凸透鏡很類似，但此時影響頻率軸上的切變不是焦距而是求面反射鏡的半徑，其矩陣形式如下:formula_72優點:在計算上只需要用到2X2的矩陣運算，避免了複雜的數學積分運算與物理理論缺點:得到的結果只有在近軸(靠近光傳播方向)準確性高當電磁波由左端平面輸入，在空氣傳播時使用菲涅爾轉換描述，接著在距離為z1處經過一焦距為f的薄凸透鏡折射後再射入空氣中傳播距離z2後由右端平面輸出，此過程經過三層LCT的操作，依序為菲涅爾轉換、薄凸透鏡折射和菲涅爾轉換。根據LCT加成性，整個LCT過程相當於作用以下三個矩陣相乘後的矩陣乘積。formula_73基本性質.在這個部份，我們列出了一些 LCT 的基本性質：給定一個二維的行向量 formula_74 ，我們列出在某些特定輸入之下 LCT 所展現的特殊性質。範例.考慮一個包含兩個雷達天線的系統，兩個天線面對面放置，其中一個為訊號發出端另一個為接收端，兩個天線之間的距離為 formula_75 ，訊號在兩個天線之間傳遞。首先，對天線 A (發出端) 來說，LCT 矩陣可以寫成以下形式：formula_76接著，對天線 B (接收端) 來說，LCT 矩陣可以寫成以下形式：formula_77最後，對於在空氣中傳遞的訊號來說，LCT 矩陣可以寫成以下形式：formula_78將三個部份合在一起，整個系統的 LCT 可以整合成：formula_79和粒子物理的關聯.研究指出費米子的某些性質可能可以在粒子物理的標準模型中，和線性標準轉換的自旋表示中建立一些關係。在這個方法中，粒子的電荷，弱超荷和弱同位旋被表示成某些用 LCT 的自旋表示相關的克里福代數定義的運運算元的線性組合。參考.Jian-Jiun Ding, Time frequency analysis and wavelet transform class note,the Department of Electrical Engineering, National Taiwan University (NTU), Taipei, Taiwan, 2009.Jian-Jiun Ding, Time frequency analysis and wavelet transform class note,the Department of Electrical Engineering, National Taiwan University (NTU), Taipei, Taiwan, 2021.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