% 创建一个包含两个频率分量的信号
Fs = 1000; % 采样率
t = 0:1/Fs:1-1/Fs; % 时间向量
x = cos(2*pi*50*t) + cos(2*pi*120*t); % 合成信号
% 使用FFT进行频谱分析
N = length(x); % 获取信号长度
X = fft(x, N)/N; % 计算FFT并归一化
f = Fs*(0:(N/2))/N; % 频率向量
% 绘制频谱图
plot(f, abs(X(1:N/2+1))) % 只绘制一半的频谱,因为X是对称的
title('频谱分析')
xlabel('频率 (Hz)')
ylabel('幅度')