MATLAB中完成以下信号的频谱图绘制，并进行反变换与原信号对比：冲激

在MATLAB中，你可以使用`fft`函数来计算信号的频谱图，`plot`函数来绘制图形，以及`ifft`函数来进行反变换。假设你有一个名为`impulse_signal`的冲激序列，下面是简单的步骤：

1. **频谱图绘制**:

% 创建冲激信号
t = 0:0.001:1; % 时间向量
impulse_signal = dirac(t); % 冲击函数

% 计算快速傅里叶变换 (FFT)
spectrogram = fft(impulse_signal);

% 取实部，因为对称性只关心正频率部分
spectrogram = abs(spectrogram);

% 转换频率轴（通常取离散时间信号的一半）
fs = 1/(t(2) - t(1)); % 采样率
frequencies = (0:length(spectrogram)-1)*fs/(length(spectrogram));

% 绘制频谱图
figure;
plot(frequencies, spectrogram);
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('Impulse Signal Spectrum');

2. **反变换并对比**:

% 反变换回时域信号
reconstructed_impulse = ifft(spectrogram);

% 比较原信号和重建信号
figure;
hold on;
stem(t, impulse_signal, 'r', 'Marker', 'o'); % 原信号标记为红色圆点
stem(t, reconstructed_impulse, 'b', 'LineWidth', 2); % 重构信号标记为蓝色线
legend('Original Impulse', 'Reconstructed Impulse');
xlabel('Time (s)');
title('Comparison of Original and Reconstructed Impulse Signal');
hold off;