fs=1000;%设定
采样频率N=100;
n=0:N-1;
t=n/fs;
f0=100;%设定正弦信号频率
%生成正弦信号
x=sin(2*pi*f0*t);
figure(1);
subplot(211);
plot(t,x);%作正弦信号的时域波形
xlabel('t');
ylabel('y');
title('正弦信号y=2*pi*100t时域波形');
grid;
%进行FFT变换并做
频谱图y=fft(x,N);%进行fft变换
mag=abs(y);%求幅值
f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
figure(1);
subplot(212);
plot(f,mag);%做频谱图
xlabel('频率(Hz)');
ylabel('幅值');
title('正弦信号y=2*pi*100t幅频谱图N=100');
grid;