MATLAB中fft频谱图分析

clear all;clc;close all;
fs=128;
N=128;
t=(0:N-1)/fs;
y1=cos(2*pi*50*t);
y2=cos(2*pi*10*t);
y=y1+y2;
Y=fft(y,N);
% fre=(0:N/2)*fs/N;
fre=linspace(0,fs/2,N/2+1);%N为偶数
% fre=linspace(0,fs/2-fs/N/2,(N+1)/2);%N为奇数

stem(fre,abs(Y(1:N/2+1)),'k');    %N为偶数
% stem(fre,abs(Y(1:(N+1)/2)),'k');  %N为奇数
xlabel('频率（Hz）');
ylabel('幅值');
title('只有正频率刻度');

set(gcf,'color','w');