11t.m

来自「包含50个数字信号处理方面的matlab程序」· M 代码 · 共 13 行

clear;
clc;
N=512;n=0:N-1;t=0.01*n;
x=sin(2*pi*5*t)+sin(2*pi*15*t)+sin(2*pi*30*t);
X=fft(x,N);
magx=abs(X);
k=[0:1:N-1];w=2*pi/N*k;
plot(k/N*100,magx);
title('FFT N=512')
xlabel('频率（单位：Hz)');
ylabel('|X|');grid
axis([0,100,0,300])