用MATLAB绘制周期方波的频谱图

在MATLAB中，你可以使用`fft`函数（快速傅立叶变换）来绘制周期信号的频谱图，例如周期方波。以下是基本步骤：

1. **生成周期方波**：
首先创建一个特定频率和周期的正弦波。假设你想制作一个频率为`f` Hz、周期为`T`秒的方波，可以这样做：

   T = 1; % 定义周期
   fs = 1000; % 定义采样率
   f = 5; % 方波频率
   t = linspace(0, T, T*fs); % 时间向量
   x = sin(2*pi*f*t); % 生成正弦波

2. **计算频谱**：
使用`fft`函数计算信号的离散傅立叶变换（DFT），这将给出各频率成分的幅度信息：

   X = fft(x);

3. **绘制频谱图**：
`abs(X)`会得到各个频率分量的绝对值，`plot`函数用于绘制：

   X_mag = abs(X)/length(x); % 对结果取归一化
   freq = (0:length(X)-1) * fs / length(t); % 计算频率轴
   plot(freq, X_mag); % 绘制频谱图
   xlabel('Frequency (Hz)');
   ylabel('Magnitude');
   title(['Spectrum of Periodic Square Wave with Frequency ' num2str(f) ' Hz']);

4. **显示图形**：

   grid on;