DSP4_快速傅里叶变换讲述.ppt

* 则： 式中： * * CZT 的实际计算方法： 1. 是 点系列，由 所决定： 2. 是双边无穷长序列，由定义所决定： * 3. 是 点序列，由需要所决定。 希望用DFT实现 截短 延拓 * * 点序列 （取前M个点） * 与本章有关的MATLAB文件主要是fft, ifft和 czt.m。fft实现快速傅立叶变换，ifft实现快速傅立叶反变换，czt.m实现线性调频Z变换 fft的调用格式是：X=fft(x), 或 X=fft(x，N) czt.m 调用格式是： X＝czt(x, M, W, A) 。x是待变换的时域信号，其长度设为N，M是变换的长度，W确定变换的步长，A确定变换的起点。若M=N, A=1, 则CZT变成DFT * 例 设模拟信号 ，以 t = 0.01n (n=0: N-1) 进行取样，试用fft函数对其做频谱分析。N分别为： (1) N=45；(2) N=50；(3) N=55；(2) N=60。 程序清单如下 %计算N=45的FFT并绘出其幅频曲线 N=45;n=0:N-1;t=0.01*n; q=n*2*pi/N;