MATLAB 绘制信号单双边频谱代码

close all
clear

T=1.5;
N=10000;
t=linspace(0,T,N);
y=2.*cos(2*pi*300.*t);%频率300Hz单音信号
fs=N/T;

%% 单边
figure
Y = abs(fft(y)/N);%双边频谱
Y = Y(1:floor(N/2+1));%单边取一半
Y(2:end-1) = 2*Y(2:end-1);%非0频幅值要乘2
f = fs*(0:(N/2))/N;%频率刻度
plot(f,abs(Y)) 
xlabel("频率（HZ）")
ylabel("幅度（V）")
grid on

%% 双边,当N为奇数时没有精确零频
figure
Y = abs(fftshift(fft(y))/N);%双边频谱
f = (-N/2:N/2-1)*(fs/N);%频率刻度
plot(f,abs(Y)) 
xlabel("频率（HZ）")
ylabel("幅度（V）")
grid on

单边频谱如图：

双边频谱如图：