untitled6.m

来自「利用DFT分析模拟信号频谱,并附有实验报告说明和源码」· M 代码 · 共 17 行

N=100;
L=512;
f1=100;
f2=110;
fs=250;
T=1/fs;
ws=2*pi*fs;
t=(0:N-1)*T;
f=cos(2*pi*f1*t)+0.75*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));
xlabel('w');
ylabel('x(t)的幅度谱')
title(['Hamming窗函数'])