用matlab画出一首歌的时域谱和频域谱

在MATLAB中，我们可以使用内置的音频处理函数来绘制歌曲的时间域波形图（时域谱）和频率域谱（频谱图）。以下是基本步骤：

1. **加载音频文件**:
使用`audioread`函数读取音频文件，例如`[y, Fs] = audioread('song.wav');`。`y`是音频数据，`Fs`是采样率。

2. **时域波形图**:
`plot(t, y)`，其中`t = (0:length(y)-1)/Fs`是时间轴，`y`是对应的音频样本值。这将显示音频随时间的变化情况。

3. **计算频域谱**:
使用`fft`函数对信号进行快速傅立叶变换（FFT），得到频域数据`Y = fft(y)；`。通常我们会对结果取实部并归一化，`P = abs(Y(1:end/2+1))/length(y);`。`P`是功率谱密度。

4. **绘制频谱图**:
`plot(frequencies, P)`，其中`frequencies = (0:length(P)-1)*Fs/(length(y)-1)`是频率轴。`P`作为数据点绘制频率分布。

5. **调整视觉效果**:
可能需要添加标题、刻度标签等，并设置合适的颜色和线型。

完整示例代码如下：

[y, Fs] = audioread('song.wav');
t = (0:length(y)-1)/Fs;
plot(t, y);
xlabel('Time (s)');
ylabel('Amplitude');

% 计算频谱
Y = fft(y);
P = abs(Y(1:end/2+1)).^2 / length(y); % 频谱密度（PSD）
frequencies = (0:length(P)-1)*Fs/(length(y)-1);

% 绘制频谱图
subplot(2,1,2);
plot(frequencies, P);
xlabel('Frequency (Hz)');
ylabel('Power/Frequency');
title('Song Spectrogram');