实验14 快速傅里叶变换(FFT)
(完美格式版,本人自己完成,所有语句正确,不排除极个别错误,特别适用于山大,勿用冰点等工具下载,否则下载之后的word 格式会让很多部分格式错误,谢谢)
XXXX 学号姓名处XXXX
一、实验目的
1、加深对双线性变换法设计IIR 数字滤波器基本方法的了解。
2、掌握用双线性变换法设计数字低通、高通、带通、带阻滤波器的方法。
3、了解MA TLAB 有关双线性变换法的子函数。
二、实验内容
1、双线性变换法的基本知识
2、用双线性变换法设计IIR 数字低通滤波器
3、用双线性变换法设计IIR 数字高通滤波器
4、用双线性变换法设计IIR 数字带通滤波器
三、实验环境
MA TLAB7.0
四、实验原理
1、实验涉及的MATLAB 子函数
(1)fft
功能:一维快速傅里叶变换(FFT)。 调用格式:
)(x fft y =;利用FFT 算法计算矢量x 的离散傅里叶变换,当x 为矩阵时,y 为矩阵x
每一列的FFT 。当x 的长度为2的幂次方时,则fft 函数采用基2的FFT 算法,否则采用稍慢的混合基算法。
),(n x fft y =;采用n 点FFT 。当x 的长度小于n 时,fft 函数在x 的尾部补零,以构成n
点数据;当x 的长度大于n 时,fft 函数会截断序列x 。当x 为矩阵时,fft 函数按类似的方式处理列长度。
(2)ifft
功能:一维快速傅里叶逆变换(IFFT)。 调用格式:
)
(x ifft y =;用于计算矢量x 的IFFT 。当x 为矩阵时,计算所得的y 为矩阵x 中每一
列的IFFT 。
)
,(n x ifft y =;采用n 点IFFT 。当length(x)n 时,
将x 截断,使length(x)=n 。
(3)fftshift
功能:对fft 的输出进行重新排列,将零频分量移到频谱的中心。 调用格式:
)(x fftshift y =;对fft 的输出进行重新排列,将零频分量移到频谱的中心。当x 为向
量时,)(x fftshift 直接将x 中的左右两半交换而产生y 。