用matlab程序写出fft谐波分析

在MATLAB中，Fast Fourier Transform (FFT)是一种常用的数学工具，用于将信号从时间域转换到频率域。下面是使用内置函数`fft()`进行简单的一维信号（如一列数或数组）谐波分析的基本步骤：

% 假设我们有一个长度为N的信号向量x
x = rand(1, N); % 创建一个随机信号作为示例

% 使用fft()函数计算FFT
X = fft(x);

% 计算信号的长度，因为fft返回的是复数数组，通常需要对实部和虚部分别处理
N = length(x);
X_mag = abs(X); % 取绝对值得到幅度谱
X_angle = unwrap(angle(X)); % 获取相位角，unwrap消除角度跳跃

% 将结果分为正负半轴，因为FFT的结果是对称的
X_half = X_mag(1:N/2+1);
frequencies = (0:N/2)*(fs/N); % 创建频率向量，假设采样率fs

% 绘制频谱图
plot(frequencies, X_half)
xlabel('Frequency (Hz)')
ylabel('Amplitude')
title('Single-Sided FFT Spectrum')

% 如果你想看到完整的频率范围和双侧谱，可以这样操作：
% X_full = [X_half, zeros(1, N/2)];
% frequencies_full = [frequencies, frequencies(end:-1:1)];