信号频谱分析FFT仿真Matlab代码中文解释

以下是一个基于Matlab的信号频谱分析FFT仿真代码，中文解释如下：

% 输入信号
t = linspace(0, 1, 1000); % 生成1000个从0到1的等差数列
x = sin(2*pi*50*t) + sin(2*pi*120*t); % 生成由50Hz和120Hz正弦波组成的信号

% 绘制时域图像
subplot(2, 1, 1); % 将画布分为两行一列，此处是第一个子图
plot(t, x); % 绘制x关于t的时域图像
title('时域图像'); % 设置子图标题
xlabel('时间 (s)'); % 设置x轴标签
ylabel('幅度'); % 设置y轴标签

% 进行FFT变换并绘制频域图像
N = length(x); % 获取信号长度
f = (0:N-1)*(1/N); % 生成频率轴
X = fft(x); % 进行FFT变换
subplot(2, 1, 2); % 将画布分为两行一列，此处是第二个子图
plot(f, abs(X)); % 绘制X关于f的频域图像
title('频域图像'); % 设置子图标题
xlabel('频率 (Hz)'); % 设置x轴标签
ylabel('幅度'); % 设置y轴标签

这个代码生成一个由50Hz和120Hz正弦波组成的信号，并绘制该信号的时域图像和频域图像。其中，时域图像显示信号在时间上的变化，频域图像则显示信号在频率上的分布情况。

代码首先生成一个等差数列t，然后通过两个正弦波的叠加生成信号x。接着，代码使用subplot函数将画布分为两行一列，并绘制时域图像。subplot(2, 1, 1)表示将画布分为两行一列，当前子图是第一个子图。plot(t, x)绘制x关于t的时域图像。title、xlabel和ylabel用于设置子图标题、x轴标签和y轴标签。

接着，代码使用fft函数进行FFT变换，并生成频率轴f。N表示信号长度，f的取值范围为0到N-1，并将其归一化为0到1之间的小数。X为FFT变换后的结果。最后，代码使用subplot函数将画布分为两行一列，并绘制频域图像。subplot(2, 1, 2)表示将画布分为两行一列，当前子图是第二个子图。plot(f, abs(X))绘制X关于f的频域图像。abs函数用于取X的绝对值。title、xlabel和ylabel用于设置子图标题、x轴标签和y轴标签。