chirp信号的分数阶fourier域的带通滤波 matlab仿真,LFM信号时频分布与分数阶域滤波仿真研究...

1LFM线性调频(LFM)信号是脉冲压缩体制雷达中广泛应用的大时宽带宽积信号,它的突出优点是匹配滤波器对回波信号的多普勒频移不敏感。LFM信号又称chirp信号,在雷达、声纳中得到了广泛应用,但由于其非平稳性,在时域频域都有较大的时宽带宽,如果采用传统处理平稳信号的滤波方法进行滤波效果不理想。本文利用FRFT进行分数阶Fourier域的LFM信号滤波,并得到了良好效果。本文给出chirp信号为128点,频率从0到0.3Hz,经过高斯幅度调制后如图一所示,加入高斯白噪声后如图二所示图一图二2时频分布1、Wigner-VilleWigner-Ville分布是时频分析的有力工具,如公式(1)(2)所示=(1)=(2)这两个表达式的等价性通过用频谱表示信号,并代人方程(1)是很容易验证的。根据Lohmann的理论研究推倒出FRFT的Wigner-Ville分布(WVD)是原信号的WVD的坐标轴旋转形式。对上文中的chirp信号进行WVD可以看出时频聚集性(近似为一条直线),如图三和图四分别为该chirp信号的WVD的等高线图和三维图。图三图四2、Wigner-HoughHough变换是Hough在1962年提出的形状匹配技术,它科将被检测图像的参数曲线在参数空间中凝聚起来形成与相应曲线对应的参数峰点,从而得到图像中各个曲线的参数。对于离散的有限图像来说,Hough变换的思想是将所有可能的线条参数组成的参数空间量化为有限的参数表。考虑如下的直线极坐标参数方程:。图像I的每一点,Hough变换会对应平面上的一个正弦波,其幅值等于点的密度。所以对I中的所有点,Hough变换对应了一系列的正弦波,而且会在平面上相互交叉。也就是说,Hough变换在图像I上沿着直线进行积分且积分值与参数相关。当对信号的Wigner-Ville分布进行Hough变换时,可以得到一种新的变换,称作Wigner-Hough变换,简写为W-H变换。例如考虑如下的信号:式中为白噪声。其W-H变换为(3)W-H变换与Wigner-Ville分布相比具有两个优点:抑制噪声以及交叉项。对上文的chirp信号的WVD进行W-H变换后得到等高线图和三维图如图五图六所示。通过W-H变换可以看到LFM信号更好的时频聚集性。图五图六3FRFTFRFT(分数阶Fourier变换,FractionalFourierTransform)为经典Fourier变换的分数阶推广,把单纯的时域频域变换引入分数阶域中,并形成了新的时频分析的有力工具,FRFT公式如(4)所示。(4)分数阶Fourier变换核如公式(5)所示(5)实质上是一组调频率为的chirp信号,其初始频率为,复包络为,分数阶Fourier域由该组完备正交基所表征,通过改变旋转角度,便可以得到不同调频率的基。而由chirp信号的解析式,可知调频率为,如果取,则该chirp信号会在FRFT的某个基上聚焦,由此可得,即为chirp信号与FRFT的某个基正交。利用FRFT建立(p,u)(p为FRFT的阶数,满足)的二维平面,其中p从0到4满足FRFT的一个周期,在此二维平面上进行二维搜索,得到输入信号的极值点和极值点坐标,Fmax=4.8347,Pmax=1.1800,u=83,提取在Pmax=1.1800时输入信号关于u的FRFT的u域分布,如图七所示。在u=83处发现峰值信号,即为有用chirp信号的1.1800阶FRFT。在u域以u=83为中心,选择窄带通滤波器进行滤波,得到滤波后信号如图八所示。图七图八对滤波后的信号进行-1.1800阶FRFT,则回复