分数傅里叶变换

1、分数傅里叶变换分数傅里叶定义：分数傅里叶变换 的物理意义即做傅里叶变换 便次，其中。不一定要为整数（比傅里叶变换 更加广泛）；通过分数傅里叶变换 之后，图像或信号便会同时拥有 时域与频域两者的特征。1.1 （维基百科）第一种定义：XO-g第二种定义:1 jCOt(!) ；一心 t.21. Ef-27rr°小血/曾吃出001.2从数学上分数傅立叶变换定义了积分形式:A. C. Mcbridc和FH Ke门布1987年给出了 V Namit的分数傅里叶变换 的枳分形式凡具体地说,对信号空间A （曲中的任何信号/，它的分数 傅里叶变换（尸/） （/）可以写成枳分形式（产"）5=，

2、/（）/（p；由（2-5）其积分核是r( )exp | i( cota-2iVcsc(a)4-r cot5 卜-(T)?(喇p = 2(2-6)p=2/l公式中各记号的含义是l-icot（a） pn2H 抗. ”5式中止整数.|）分数傅里叶变换的级次,可取任何实数.(2-7)Wigner分布函数相空间定义的分数傅立叶变换A.W.Lohmann在1993年利用傅里叶变换相当于在 Wigner分布函数相空间中角度为 兀 /2的旋转这一性质，说明分数傅里叶变换在Wigner分布函数空间中相当于角度是 p u /2的旋转，这里，p是分数傅里叶变换的级次。因此，A.WLHimw】in定义第次是p的分数傅

3、里叶变换（户"）（/）为其中，矩阵网川是时-频相平面，炉上角度为（即/2）的旋转矩阵分数傅里叶变换的定义在数学上是等价的。当分数傅里叶变换的哥次p从0连续增长到达1时，分数傅里叶变换的结果相应地从原始信号的纯时间（空间）形式开始逐渐变化成为它的纯频域（谱）形式，哥次p在0到1之间的任何时刻对应的分数傅里叶变换采取了介乎于时（空）域和频域之间的一个过渡域的形式，形成一个既包含时（空）域信息同时也包含频（谱）域信息的混合信号。因此，这样定义的分数傅里叶变换确实是一种时（空）频描述和分析工具分数傅里叶的分类：1 .一维分数傅里叶变换分数傅里叶变换的数学表达式有 积分形式和级数表达式两种等价

4、形式，2 .积分形式2级数表达式形式其中3 .二维分数傅里叶变换其中C为相应常系数。当a=b时，上式就是二维分数傅里叶变换的表达式；当a=b=1时，上式转化为常规二维傅里叶变换 ；当a与b不相等时，我们称这种情况 的二维分数傅里叶变换为不对称分数傅里叶变换。此时在x、y方向实施的变换级次是不同的。分数傅里叶变换的性质1周期性：（k为整数）九3 （ ）=尸""/=尸'/（*） =人（"）2线性：（c1和c2是复常数）尸</（，） +X 3（"） + %/（"）3阶数可加性：尸'/(叨=尸'"/(叨=加(&q

5、uot;)4尺度变换特性：2 八一l-/cot a2 L：exp imr cot a 1 lWcotacos- a 1一 cos2 a J'Afasin 优、i sin a )M 是实数，.FL M 不为零和无穷, a= *疗/2 , af - arctan(A/2 tan a) / = 2(7 7T5时移特性：Ff, /(*一0)=exp (inlr sin a cos a |exp(-/ksiiiiz)i /-Acos a ：6频移特性：尸 exp (ilrrbx f( x)j = exp I -ilrrlr sin a cos a 卜exp | -il7tnb cos a) /；

6、("一sin a)7可逆性：对一个函数进行 P级分数傅里叶变换后，接着进行- P级的分数傅里叶变换， 则可得到原函数：PL /(x) =尸。/(*) = *)分数傅里叶变换的数值算法户二£/(切川(i)基于傅立叶变换矩阵因子哥的离散化算法，利 自 来计算离散的分数傅立叶变换的核矩阵，从而利用FFT来计算离散分数傅立叶。其中W是离散傅立叶变换核矩阵(2)基于正交投影的离散化算法，对连续分数傅立叶变换的特征函数 进行离散化近似和正交投影，得到一组与 Hermite-Gaussian函数形状相似的离散傅立叶变换矩阵的正交化离散Hermite特征向量。然后，仿照连续分数傅立叶变换的

7、核函数谱分解表达式，构造了离散分数傅 立叶变换矩阵 (3)基于chirp分解的离散化算法。根据分数傅立叶变换的表达式，将分数傅立叶变换分解为信号的卷积形式后，直接离散化，利用FFT来计算分数傅立叶变换。2.2.1基于傅立叶变换矩阵1基于chirp分解的离散化算法(FRFT.m和frft22d.m)对应的分数傅里叶函数：21m/xcsc 由 »(工)由rFfl/(w) = cjexp(ijru cot 以1)输入信号f(x)与喟啾信号相乘(j将两个信号分别离散化);2)进彳T FT运算；3)进行尺度变换，系数为csc(f);4)再与喟啾信号相乘；5)最后与常数位相相乘。2基于正交投影的

8、离散化算法(Disfrft.m和cdpei.m)其计算过程如下：(1)计算矩阵s的正交特征向量勺(*)(2)计算连续Hermite*Gaussian函数的取样向最应(上) (#)= /1 I-此Q-37) 将.(司进行平移得到引'N_%(灯一1(2-38)当N为偶数£(*)=2应(左一川)彳一1兰亚三"一1当N为奇数(2-39)（4） £（竹投影到匕&）空间得到吃停）吃伏卜S （吃住）嗔町”/上）（2-40）|"#i mod4口（5）对和（k）进行正交归一化得到/（上）,可以采用了两种不同正交映 射的方法，并根据采用的正交化方法不同而分别

9、命名OPA方法和GSA方 法口6） 再由灯构造U山口同I"*-（2-41）（7）构造"严A'-* Aiuod? I(2-42)（8）最后由少和U得到离散分数傅立叶变换矩阵尸FF = U»U（2-43）离散分数傅立叶算子尸，无论p为任迓实数，均满足特征方程Q-27）,并 且澜足阶数可加性产尸=产:。信号上伏）的离散分数傅立叶变换/任）通过如下公式计算：*<）=尸（2-44）与连续的分数传立叶变换相似，信号工。）也可以通过其逆变换恢复;（2-45）基于正交投影的离散化算法 对连续分数傅立叶变换的特征函数进行离散化近似和正交投影，得到一组与 Hermite

10、-Gaussian函数形状相似的离散傅立叶变换矩阵的正交化离散Hermite特征向量。然后，仿照连续分数傅立叶变换的核函数谱分解表达式，构造了离散分数傅立叶变换矩阵。此算法适合用来计算连续的分数傅里叶变换。基于chirp分解的离散化算法，将分数傅立叶变换分解为信号的卷积形式后，直接离散化，利用FFT来计算分数傅立叶变换。图像的分数傅里叶变换对图像进行分数傅里叶变换分析的目的是确定图像经过分数傅里叶变换后的特性表现，主要包含分数傅里叶变换对图像能量分布和频率分布影响两方面的内容。其中能量分布表现分数傅里叶变换图像的能量聚积性与分数变换阶数的关系，能量聚集性强烈地依赖于其接近于傅里叶变换的程度；频

11、率分布表现在分数傅里叶变换的相位函数包含了图像的纹理频率信息，变换阶数不同，相位函数所含的图像边缘高频信息也不相同。图像经过某种二维离散变换之后的能量分布体现了图像的变换特征。图像分数傅里叶变换域的能量分布特点是：能量向中心区域聚集性。(1)当分数阶次p由小变大时，由相位函数恢复的图像呈现出图像边缘轮廓变得越来越 清晰，这类似于原始图像经历了不同截止频率的高通滤波器。当p较小时对应于截止频率较低的高通滤波器，低频成份浮现出来，图像边缘模糊；当p较大时，对应于截止频率较高的高通滤波器，大部分低频成份被滤掉，图像边缘比较清晰，FRFT逐渐向FT退化。(2)当变换阶数p由小变大时，仅由幅度函数恢复的

12、图像越来越接近原图像的背景，这 类似于原图像经历了不同截止频率的低通滤波器。p较小时，对应于截止频率较高的低通滤波器，高频分量残留较多，能清晰看到原图像的轮廓；p较大时，对应于截止频率较低的低通滤波器，大部分高频分量被滤出只显现原图像背景。(3)当变换阶数p为其它值时，由FRFT相位函数和幅度函数所恢复的图像既包含了原 图像的背景信息又包含了原图像的纹理频率信息。由此可以推论这类似于原图像经历了 FRFT的时频滤波，也即将时频平面旋转某一角度后再进行滤波。假如频域滤波器截止频率 和带宽固定，当旋转的角度不同 (阶数不同)时，时间轴和频率轴上的投影不同，所以频域滤 波器输出的频率成分也不同。 这

13、表现在恢复的图像上即为相位函数和幅度函数包含的频率成 份随阶数而变化。当变换阶数较小时，由图像的FRFT的幅度函数和相位函数恢复的图像都显示出很强的 图像信息，体现出了较强的空域特性；当变换阶数逐渐增大时，图像的 FRFT的幅度函数恢 复的图像所包含的原图像的空域特征逐渐减弱直至消失，相位函数恢复的图像包含原图像的边缘纹理特征逐渐增强。当FRFT的变换阶数增大到一定程度时，其幅度和相位特征越来越接近FT域即频域特征。这些结论有力地体现了FRFT域的空-频双域特征。采用分数傅立叶变换的图像边缘提取方法对图像作连续小级数的分数傅立叶变换，相当于对图像作连续的微小变换，当分数级次很小时，肉眼几乎看不

14、出与原图的区别， 当级次略有增加，图像边缘与原图有了明显的区别， 当继续缓慢增加分数级次时，图像与原图明显不同。 通过分析可以看出， 图像中对比度低的区域随级次变化缓慢，对比度高的区域(即图像边缘)随级数变化快。由此，取不同级数的 分数傅立叶变换后的图像减去原图像，即可得到图像的边缘。 不同的分数级数对应不同的形变，选取不同级数变换后的图像相减，即可提取不同尺度的边缘。图像分数阶Fourier变换的幅度和相位信息假设 F(k,h)是二维图像 f(x, y)的二维 Fourier 变换 F(k,h)=FT2Df(x,y)(8)我们可以把F (k,h )分解成幅度部分和相位部分，即F (k,h 尸

15、 F (k,h j_|P(k,h )=A(k,h ),P(k,h )(9)其中A(k,h尸F (k,h，为幅度函数，P(k,h尸F (k,h，A(k,h )为相位函数,结论:1 .当变换阶数P由小变大时，仅由相位函数恢复的图像， 显现原图像的边缘越来越清晰， 这类似于原图像经历了不同截止频率的高通滤波器。p较小时(0.01)对应于截止频率较低的高通滤波器，低频成份浮现出来，使提取的边缘模糊，如图 4(b)所示；p较大时(0.8),对应于 截止频率较高的高通滤波器，大部分低频成份被滤出，提取的边缘较清晰，如图5(d)所示,此时FRFT基本退化为FT。2 .同理，当变换阶数P由小变大时，仅由幅度函

16、数恢复的图像越来越接近原图像的背景，这类似于原图像经历了不同截止频率的低通滤波器。p较小时(0.01),对应于截止频率较高的低通滤波器，高频分量残留较多，还能清晰看到原图像的轮廓，如图4(a)所示；p较大时(0.8),对应于截止频率较低的低通滤波器，大部分高频分量被滤出，此时仅显现原图像的背景，如图6(c)所示。3 .当变换阶数P为其它任意值时，由FEFT相位函数和幅度函数所恢复的图像既包含了 原图像的背景又包含了原图像的纹理，如图4(b)、图4 (c)、图5 (a)、图5 (b)所示。针对这种情况，我们可以推论这类似于原图像经历了FRFT的时频滤波，也即将时频平面旋转某一角度后再进行滤波。假如频域滤波器截止频率和带宽固定，当旋转的角度不同(阶数不同)时，时间轴和频率轴上的投影不同，所以频域滤波器输出的频率成分也不同。表现在恢复的图像上即为相位函数和幅度函数包含的频率成份随阶数而变。分数傅里叶变换(采用chirp信号的方法)a从0.1-1过程中，相应的分数傅里叶变换。相应的幅值信息(a从0.1-1)相应的相位信息（即：f/|f|,其中f为相应分数傅里叶变化后，输出的信号 ）相应的幅值信息进行逆变换根据相位信息进行相应的逆变换，可以发现图像的边缘信息，消除了图像的背景信息。总结通过分数傅里叶变换和其逆变换可以找到一些图像的主体边缘信息，但是有些图像边缘的采集并不是很好，对