特征分解型离散分数阶Fourier变换
来源:微智科技网
硼PA &AP霞 【本文献信息】樊孝明,阎鹏,王玫.特征分解型离散分数阶Fourier变换[J].电视技术,2012,36(15). ■ 教蠢 i dto屯^qi t cft“ 特征分解型离散分数阶Fourier变换 樊孝明,阎鹏,王玫 (桂林电子科技大学信息与通信与学院,广西桂林541004) 【摘要】分数阶Fourier变换作为Fourier变换的广义形式,广泛应用于科学计算和研究,离散分数阶Fourier变换是其得以应用 的关键。特征分解算法是由可交换对角矩阵得到近似连续Hermite-Gaussian函数的特征向量,再对Hermite—Gaussian函数进行 加权和运算。对一种基于数特征分解的方法进行了改进,并进行计算机仿真。仿真结果表明所得的Hermite—Gaussian函数与连 续函数的近似度更为优异,从而提高了离散分数阶Fourier变换的近似度。 【关键词】分数阶Fourier变换;离散分数阶Fourier变换;特征分解 【中图分类号】TN911 【文献标识码】A Discrete Fractional Fourier Transform Algorithm Based on New Eigendcomposition Algorithm FAN Xiaoming,YAN Peng,WANG Mei (School ofInformation and Communication Engineering,Guilin University ofElectronic Technology,Guangxi Guilin541004,China) 【Abstract】Fractional Fourier transform(FRFr)as a Fourier transform of the generlaized form is widely used in scientiifc computing and engineeirng re- search,and discrete Fourier transform(DFRFF)algorithm is the key for application.Eigendcomposition algorithm is based on the tridiagonal commuting matrix which provides sample approximations of the continuous Hermite—Gaussian like functions,and the DFRFr is interpreted as a weighting summation of Hermite—Gaussian functions.An improved algoritm ihs simulated.The simulation results show that the obtained Hermite-Ganssian function and contin・ UOUS function approximation degree more excellent.SO as to increaSe the discrete Fourier transform fractional order approximate degrees. 【Key words】fractional Fourier transform;discrete fractional Fourier rtansform;eigendcomposition 近年来,分数阶Fourier变换(FRT)在光信号处理、量 子力学、时频域分析、光滤波器等领域得到广泛应用,离散 化是其应用的关键。当前主要的离散算法有直接型算法、 r _二 eJ州 “ ”“ ( ,u)={6( 一M), =n1r L6(t+M), =(2n+1)啊 (2) 线性加权型算法和特征分解型算法。线性加权性算法不 能逼近连续分数阶Fourier变换,采样型算法不能满足分 数阶Fourier变换特有的旋转相加性,因此不是严格的离 散化算法。特征分解型是目前唯一严格定义了DFRFF的 算法,这种算法与连续分数阶Fourier变换十分相近,同时 还有“分数阶”的特有性质:旋转相加性。这样正逆变换 就可以用相同的表达式描述,只是相差一个参数,用一个 (t,“)可以写成以下形式 (£,“):∑e。j Hn(t)Hn(u) ,’(3) (t)为n阶Hermite—Gaussian函数(HGF),公式为 ÷ . (t)=  ̄/2 n! ( £)e (4) h (t)为//,阶Hermite—Gaussian多项式。FRFT在数 统一的程序可以把分数阶Fourier域应用到时域。本文对 学上是普通Fourier变换的广义形式,并且FRFI"是Her- 种特征分解型算法进行改进,仿真结果表明新算法具有 mite—Gaussian函数的加权和,加权系数为输入信号和 一更高的近似性。 Hermite—Gaussian函数内积与相位因子e_j 的乘积。 1 连续分数阶Fourier变换 X。为FRFF算子,时域连续信号 (t)的a阶FRFr定 义为 2 特征分解型DFRFF算法 特征分解型DFRbq'定义为 xo(“)= ( ) ( )d£ (1) 。 VD。vT={kl N=. O2 =f∑e。。l, T ,Ⅳ为奇数 (、 5) 式中:Ka(t,U)为核函数,以2订为周期。 【∑e~', T+e~ l,T ,Ⅳ为偶数 基金项目:国家自然科学基金项目(61172054;61071088) 54 《电视技术》第36卷第15期(总第391期)I投稿网址http://www.VideoE.cn 囊秘教堪idt oqEngi n e cring lp ARTS & U晃 C A与TIO应N用S l器lI ■ &APP式中: : 。 。…‘l, ‘l, ],Ⅳ为奇数L;l, 是 [',o I',l l…l', I v,v],N为偶数 个零交点的k阶近似离散HGF,通过与F的可交换对角 量更接近采样的HGF,故k。要比k 小很多。当k。非零 时,klS+.i} T和S+( ) 有相同的特征值,所以用 + 矩阵S特征向量的标准正交化得来;D是对角矩阵,其对 kT,k≥0来求取F的特征向量。由S+.j} 求取特征向量的相似度要优于由S求出的 角线元素与 中的列向量相对应。 。= (diag 1 (,e~,0≯…,e 一 ,e一),为N’ 偶数 ㈤ 特征向量,近似均方根误差(RMSE)比较小。所得到的 DFT特征向量与连续Fourier变换的特征函数的近似程度 DYF变换矩阵特征值的赋值规则如表1所示 很高。基于 (或 +kT)的DFRFT定义如下 表1 DFRb3"特征值赋值规则 f∑e g g ,N为奇数 =GD。G =J 。 I 4m毫 — I | I e ;k=o 3I“ 每 ≥ 薯毫I 【∑e g,g T+e g,gT,Ⅳ为偶数 4m“ | ,k '2 l4m Il4m [goI g ̄ "'"I g N-2 I g ̄_’ 是 n I| ^ 21 …。 ….。。l…I g l g ],为偶数 N(或 + )的r阶归一化DFT向量。 _ 一一一一. 坚 值得注意的是当Ⅳ为偶数时,最后一个特征值有跳 4 实验仿真 2 1 O 9 8 7 6 跃,这与特征值的多样性是一致的。 实验1:为了求取S+j} 中k的最佳取值,首先计算 S+kT的特征向量,然后对取得的Ⅳ个特征向量分别与连 3 改进型特征分解算法 续HGF进行比较,计算其均方差的总和。分别对Ⅳ取值 为了改进特征分解算法,引入一个新的可交换对角 25,50,75,100,k的取值为0~54。仿真结果表明当k=15 镇阵 .公式为 的时候均方差最小,如图1所示。 0.5 s寺, 夏孚 。 。 … N cos-o 8。—- ̄T__一 s导) 。 。 l \ Ic0s 】 \ cosIN2)7tc osIN1) ̄TNN——2 —【 】 em寺 9 ~ ● H- + ●●●● — (7) O 5 1O l5 20 25 3O 35 40 45 50 的特征向量也是F的特征向量。由于S和 的特 k 征向量都是近似HGF,而S与 是线性无关的,所以 和 图1 N=25时的均方差 的线性组合k mS+|ic: 的特征向量也是近似离散HGF。 l+『 实验2:对函数 (f):l一 ,I t l≤2分别进行连 通过实验,发现: 0,其他 ]< Ⅲ ]> ). 续FTR和基于S, , +15T的DFRT。N=64,取样间隔 [凡]> ( Ⅲ < ( )。 为÷,对变换阶数取0.1~1。然后计算其均方根误差 (RMSE),仿真结果如图2所示。 式中: [ ]和 [rt]是分别由S和 生成的特征向量; h(t)是连续HGF。如果将Js和 进行适当的线性组合, 通过这两个仿真可以发现这种方法求得的特征向量 与连续HGF的近似度要优于以往的方法。 将得到十分近似 ( )的特征向量。由于T的特征向 (下转第63页) 投稿网址http://www.VideoE.cn l《电视技术》第36卷第15期l总第391期) 55 5 4 ■ 我囊 id|o Engi#ct ring IN R TER 磷 硼 储架构,外部存储器SDRAM存储6帧HD图,每一帧图采 用顶底场数据分开的方式存储在SDRAM中,整个MC单 元采用并行多级流水线处理,实现了对外部帧存储单元的 高效存取访问以及运动补偿的高效计算,处理1 Mbyte的 MC操作需要大约520个时钟周期,相比文献[6]节省了 15%。综合结果表明,本设计满足高清视频解码的要求。 参考文献: [1]Document JVT-C167.Draft ITu—T recommendation H.264(a.k.a”H. IEEE International Conference on Solid—State and Integrated Circuit Technology,2010.[S.1.]:IEEE Press,2010:445--447. a1.A novel VLSI architecture of [6] ZHENG Junhao,GAO Wen,WU D,etmotion compensation for multiple standards[J].IEEE Trans.Consumer Electronics,2008,54(2):687-694. C Y,CHEN T C,CHEN T W,et a1.Bandwidth optimized motion [7] TSAIcompensation hardware design for H.264 HDTV decdeor[C]//Proc. 48th Midwest Symposium on Circuits and Systems,2005.[S.1.]:IEEE Press,2005,2:1 199-1202. ◇ 26L”)[S].2002. [2]WANG S Z,LIN T A,LIU T M,et 1.A anew motion compensation design f0r H.264/AVC decoder[C]//Proc.IEEE International Symposium on Circuits and Systems,2005.[S.1.]:IEEE Press,2005:4558-4561. [3]KIM J H,HYUN G H,LEE H J.Cache organizations for H.264/AVC motion compensationf C]//Proe.13tll IEEE International Conference on 作者简介: 曹超(1987一 )。硕士生。研究方向为信号检测与信号处理、专用集 成电路设计与应用; 虞礼贞(1948一 )。教授。研究方向为信号检测与信号处理、专用集成 电路设计与应用; 张延军(1978一 )。副研究员。研究方向为专用数字信号处理器设计、 嵌入式图形处理器: RTCSA,2007.[S.1.]:IEEE Press,2007:534-541. [4]rrU-T H.264建议书[S].2005. [5]LIN C F,CHUNG C C,TSAI Y C,et a1.Bandwidth-efifcient architecture 研究方向为VLSI超大规模集成电路 李广桢(1981一 ).Ic工程师。设计、多标准视频解码器芯片设计。 责任编辑:薛京 收稿日期:2012-02-26 desing ofr motion compensation in H.24/AVC 6decoder[C]//Proc.10th (上接第55页) 0.015 l—・一基于踟阵 }一争一基于鹿阵 — 参考文献: [1]陶然,邓兵,王越.分数阶傅里叶变换及其应用[M].北京:清华大学 出版社,2009. 0.010 ,。 / \ l—e一基于 l5丁矩阵 \ \ / ./ 一一{ \ O.0o5 { , \. \ \ 、 [2]OZAKTAS H M,ARIKAN O.Digital computation of the fractional Fou・ iter trnsfaorm[J].IEEE Trns.Sianagl Processing,1996,44(9): 2141-2150. r [3]CANDAN C,KUTAY M A,OZAKTAS H M.The discrete fractional Fou— rier trnsfaorm[J].IEEE Tram.Signal Processing,2000,48(5): 1329—1337. P /’ — . [4]PEI S C,YEH M H,TSENG C C.Discrete fractional Fourier transfom rO 0.1 O.2 0-3 O-4 O.5 0.6 O.7 0.8 O.9 1.0 阶数 based on orthogonal projcteions[J].IEEE Trans.Signal Processing,1999, 47(5):1335-1348. 图2基于3种矩阵的DFRFT的RMSE [5]TSENG C C.Eigenvalues and eigenvectom of generalized DFT,general— 5 结束语 izod DHT,DCT—IV and DST-IV matirces[J].IEEE Trm.Siganal Pro- cessing,2002,50(4):866-877. 本文对一种基于特征分解的DFRFF进行了改进,其 [6]PEI S C,DING J J.Eigen functions of the offset Fouirer,fractionla Fou— 特征向量为离散近似Hermite—Gaussian函数,绝大多数特 iter,andlinear canonicaltrnsfaorms[J].J.Opt.SOc.Amer.A,Opt.Image 征向量与连续Hermite—Gaussian函数的近似度要优于以 往的方法。同时对此特征值和特征向量的性质做了严谨 Sci.,2003,20(3):522-532. 的研究。更进一步将S与 进行适当的线性组合,用其 线性组合计算出的特征向量与Hermite—Gaussian函数的 ◇ 作者简介: 樊孝明(1971一)。讲师.主要研究方向为扩频通信、超宽带通信等 研究: 近似度更为优异,从而提高了整体DFRFT计算的准确性。 由于特征分解满足分数阶的旋转相加性,特别适合应用于 阎鹏(1983一)。硕士生.主要研究方向为超宽带网络的定位技术。 雯 收稿日期:2012-03-13 LFM信号的处理。 责任编辑:时投稿网址hnp://www.VideoE.clfI《电视技术》第36卷第15期(, ̄391期) 63