《ch2_4DFT分析信号频谱》.ppt

N=30, N=30, L=64, = 600/64 N=30, L=128, = 600/128 N=30, L=256, = 600/256 五、DFT参数选取 1. 2. 3. 抽样频率: 抽样间隔: 抽样时间: 抽样点数: * * ch2_4DFT分析信号频谱 问题提出： 如何利用数字方法分析信号的频谱？ 有限长序列 的傅立叶变换DFT. DFT可以直接计算周期序列的DFS 可否利用DFT分析以上四种信号的频谱？ 基本原理 利用信号傅立叶变换具有的信号时域与频域之间的对应关系，建立信号的DFT与四种信号频谱之间的关系。 时域的离散化 时域的周期化 频域周期化 频域离散化 二、四种信号的时域与频域对应关系 FT FS DTFT DFS 抽样 离散化 周期化 三、利用DFT分析连续非周期信号的频谱 DFT实现 假设连续信号持续时间有限，频带有限 讨论： （1） 无限长，其频带有限 加窗 抽样 DFT （2） 有限长，其频带无限 抽样 DFT （3） 无限长，其频带无限 加窗 出现三种现象：混叠、泄漏、栅栏 抽样 DFT 四、混叠现象、泄漏现象、栅栏现象 （1）混叠现象：减小抽样间隔T，抗混滤波 抗混滤波 抽样间隔T 抽样 DFT （2）泄漏现象：选择合适的窗函数 其中： 加窗 DFT 矩形窗 窗函数一： 时域波形 幅度频谱 矩形窗： 矩形窗： 主瓣在处有一个峰值，表示其主要是由直流分量组成。由于矩形窗函数在其两个端点的突然截断，使得频谱中存在许多高频分量。 汉宁窗（hanning)： 窗函数二： 时域波形 幅度频谱 窗函数三： 哈明窗（hamming) 时域波形 幅度频谱 窗函数四： 布拉克曼窗（Blackman） 时域波形 幅度频谱 窗函数五： 凯塞窗（Kaiser） 时域波形 幅度频谱 为了说明时域加窗对连续信号频谱分析的影响，现分析一无穷长的余弦信号的频谱。 加窗 抽样 DFT 加窗 抽样 DFT 已知一连续信号为 若以抽样频率 Hz对该信号进行抽样，试求由DFT分析 其频谱时，能够分辨此两个谱峰所需的最少样本点数。 矩形窗 信号样点数 N=30 信号样点数 N=20 加矩形窗 利用矩形窗计算有限长余弦信号频谱 N=30; %数据的长度 L=512; %DFT的点数 f1=100; f2=120; fs=600; %抽样频率 T=1/fs; %抽样间隔 ws=2*pi*fs; t=(0:N-1)*T; f=cos(2*pi*f1*t)+cos(2*pi*f2*t); F=fftshift(fft(f,L)); w=(-ws/2+(0:L-1)*ws/L)/(2*pi); hd=plot(w,abs(F)); ylabel('幅度谱') 已知一连续信号为 若以抽样频率 Hz对该信号进行抽样，试求由DFT分析其频谱。 矩形窗 N=25 矩形窗 N=50 海明窗 N=25 海明窗 N=50 利用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('幅度谱') （3）栅栏现象：序列后补零，ZFFT （3）栅栏现象：序列后补零，ZFFT 解： ，m=0,1,2,3 ，m=0,1,2,…,7 * *