matlab频域图代码

以下是在Matlab中绘制频域图的示例代码：

% 生成测试信号
fs = 1000; % 采样频率
t = 0:1/fs:1-1/fs;
x = sin(2*pi*50*t) + sin(2*pi*150*t) + sin(2*pi*250*t);

% 计算快速傅里叶变换（FFT）
N = length(x);
X = fft(x);
f = (0:N-1)*(fs/N);

% 绘制幅度谱和相位谱
subplot(2,1,1);
plot(f,abs(X));
title('幅度谱');
xlabel('频率（Hz）');
ylabel('幅度');
subplot(2,1,2);
plot(f,angle(X));
title('相位谱');
xlabel('频率（Hz）');
ylabel('相位（弧度）');

以上代码示例中，首先生成了一个包含三个正弦信号的测试信号，然后使用FFT计算快速傅里叶变换，得到频域信号的幅度谱和相位谱。最后使用Matlab的subplot函数将幅度谱和相位谱分别绘制在两个子图中。

需要注意的是，绘制频域图前需要对信号进行合理的处理，如去除直流分量、补零等，以避免频谱泄漏和频率分辨率不足等问题。