数字信号处理ch24DFT分析信号频谱课件.ppt

DFT分析信号频谱 混叠现象、泄漏现象、栅栏现象 2. 泄漏现象：选择合适的窗函数 凯塞窗（Kaiser） 窗函数五： 时域波形 幅度频谱 DFT分析信号频谱 窗函数类型 主瓣宽度 旁瓣峰值衰耗(dB) 矩形 4p / N -13 Hanning 8p / N -31 Hamming 8p / N -41 Blackman 12p / N -57 Kaiser（ ） 10p / N -57 86 . 5 = b 混叠现象、泄漏现象、栅栏现象 2. 泄漏现象：选择合适的窗函数 常用窗函数特性 例：为了说明时域加窗对连续信号频谱分析的影响，现分析一无穷长的余弦信号的频谱。 加窗 抽样 DFT 加窗 抽样 DFT 例：已知一连续信号为 若以抽样频率 Hz对该信号进行抽样，试求由DFT分析 其频谱时，能够分辨此两个谱峰所需的最少样本点数。 矩形窗 例：已知一连续信号为 若以抽样频率 Hz对该信号进行抽样，试求由DFT分析 其频谱时，能够分辨此两个谱峰所需的最少样本点数。 利用矩形窗计算有限长余弦信号频谱 N=30; %数据的长度 L=512; %DFT的点数 f1=100; f2=120; fs=600; %抽样频率 T=1/fs; %抽样间隔 ws=2*pi*fs; t=(0:N-1)*T; x=cos(2*pi*f1*t)+cos(2*pi*f2*t); X=fftshift(fft(x,L)); w=(-ws/2+(0:L-1)*ws/L)/(2*pi); plot(w,abs(X)); ylabel(‘幅度谱’); 例：已知一连续信号为 若以抽样频率 Hz对该信号进行抽样，试求由DFT分析 其频谱时，能够分辨此两个谱峰所需的最少样本点数。 信号样点数 N=30 信号样点数 N=20 加矩形窗 例：已知一连续信号为 若以抽样频率 Hz对该信号进行抽样，试求由DFT分析其频谱。 利用Hamming窗计算有限长余弦信号频谱 N=50; %数据的长度 L=512; %DFT的点数 f1=100;f2=150; fs=600; %抽样频率 T=1/fs; %抽样间隔 ws=2*pi*fs; t=(0:N-1)*T; f=cos(2*pi*f1*t)+0.15*cos(2*pi*f2*t); wh=(hamming(N))'; f=f.*wh; F=fftshift(fft(f,L)); w=(-ws/2+(0:L-1)*ws/L)/(2*pi); plot(w,abs(F)); ylabel('幅度谱') 例：已知一连续信号为 若以抽样频率 Hz对该信号进行抽样，试求由DFT分析其频谱。 矩形窗 N=25 矩形窗 N=50 海明窗 N=25 海明窗 N=50 DFT分析信号频谱 混叠现象、泄漏现象、栅栏现象 3. 栅栏现象：序列后补零，ZFFT DFT分析信号频谱 混叠现象、泄漏现象、栅栏现象 3. 栅栏现象：序列后补零，ZFFT * 数字信号处理(Digital Signal Processing) 离散傅里叶变换(DFT) 问题的提出 有限长序列的傅里叶分析 离散傅里叶变换的性质 利用DFT计算线性卷积 利用DFT分析信号的频谱 DFT分析信号频谱 利用DFT分析信号频谱 问题的提出 四种信号频谱之间的关系 利用DFT分析连续非周期信号频谱 混叠现象、泄漏现象、栅栏现象 DFT参数选取 工程实际应用 DFT分析信号频谱 1.连续时间非周期信号 图1 连续非周期信号及其频谱 问题的提出 DFT分析信号频谱 2.连续时间周期信号 问题的提出 图2 连续周期信号及其频谱 DFT分析信号频谱 3.离散时间非周期信号 问题的提出 图3 离散非周期信号及其频谱 DFT分析信号频谱 问题的提出 4.离散时间周期信号 图4 离散周期信号及其频谱 DFT分析信号频谱 问题的提出 如何利用数字方法分析信号的频谱？ DFT分析信号频谱 问题的提出 有限长序列xN[k]的傅立叶变换DFT DFT可以直接计算周期序列的DFS DFT分析信号频谱 问题的提出 可否利用DFT分析以上四种信号的频谱？ 基本原理 利用信号傅立叶变换具有的信号时域与频域之间的对应关系，建立信号的DFT与四种信号频谱之间的关系。 时域的离散