MATLAB画正弦信号幅度谱相位谱

MATLAB是一种强大的数学软件，常用于科学计算和工程分析。绘制正弦信号的幅度谱和相位谱通常涉及频域分析，可以使用内置函数如`fft()`（快速傅立叶变换）来进行。

1. **幅度谱**（Magnitude Spectrum）：正弦信号通过FFT会分解为频率分量，幅度谱表示每个频率分量的大小。你可以按照以下步骤操作：

   % 创建一个正弦信号
   t = linspace(0, 2*pi, 1000); % 时间范围
   x = sin(t);

   % 计算幅度谱
   X_mag = abs(fft(x)); % 使用abs()获取复数的模值

   % 将结果转换为频率轴（以Hz为例）
   fs = 1/(t(2)-t(1)); % 采样率
   f = (0:length(X_mag)-1)*fs / length(X_mag); % 频率向量

   % 绘制幅度谱
   plot(f, X_mag);
   xlabel('Frequency (Hz)');
   ylabel('Magnitude');
   title('Amplitude Spectrum of a Sinusoidal Signal');

2. **相位谱**（Phase Spectrum）：相位谱则包含每个频率分量相对于原信号相位的信息。计算相位需要对FFT结果取角（unwrap如果存在负角度），例如：

   X_phase = angle(fft(x));
   X_phaseunwrap = unwrap(X_phase); % 解缠绕相位

   % 绘制相位谱
   plot(f, X_phaseunwrap);
   xlabel('Frequency (Hz)');
   ylabel('Phase (radians)');
   title('Phase Spectrum of a Sinusoidal Signal');