ch2_4dft分析信号频谱.pptVIP

DFT分析信号频谱;问题提出：;有限长序列 的傅立叶变换DFT.;可否利用DFT分析以上四种信号的频谱？;二、四种信号的时域与频域对应关系 ;抽样;讨论：;（2） 有限长，其频带无限;（3） 无限长，其频带无限;四、混叠现象、泄漏现象、栅栏现象;（2）泄漏现象：选择合适的窗函数;矩形窗;矩形窗：;矩形窗：;汉宁窗（hanning)：;窗函数三：;窗函数四：;窗函数五：;为了说明时域加窗对连续信号频谱分析的影响，现分析一无穷长的余弦信号的频谱。 ;;加窗;已知一连续信号为 若以抽样频率 Hz对该信号进行抽样，试求由DFT分析 其频谱时，能够分辨此两个谱峰所需的最少样本点数。 ;信号样点数 N=30;利用矩形窗计算有限长余弦信号频谱 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分析其频谱。 ;利用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;N=30, ;五、DFT参数选取;;