窗长,DFT点数对dft影响
%x(t)=cos(2pi#f1t)+cos(2pif2t).f1=100HZ,f2=120HZ,fsam=600HZ.
clc;
clear all;
N=1024;
fs=600;
f1=100;f2=120;
t=(0:N-1)/fs;%t处出错了,应该用【0:N-1】乘取样间隔而不是(0:N-1)fs/N,这是频域表示
x=cos(2pif1t)+cos(2pif2t);
Xm=fftshift(fft(x));
%Xm=(fft(x));
w=(0:N-1)fs/N
w=w-fs/2;
plot(w,abs(Xm));
title(‘窗长1024,1024点DFT’)
figure(2);
N=30;
%N=1024;
L=512;
%L=1024;
f1=100;
f2=120;
fs=600;
T=1/fs;
ws=2pifs;
t=(0:N-1)T;
x=cos(2pif1t)+cos(2pif2t);
X=fftshift(fft(x,L));
w=(-ws/2+(0:L-1)ws/L)/(2pi);
plot(w,abs(X));
%第二个程序中涉及到N,L,N*T记过就是信号长度,即矩形窗的长度,可以认为此程序中矩形窗长为30,在fs=1/T确定的情况下,N越多
% 截取的时间越长,描述的信号频谱有效主瓣宽度越窄,信号细节越丰富,fs=1/T,表示频率分辨率。L是做DFT的点数,L越大显示出的
% 图形越接近于街区信号的原本频率图,图形越平滑