用matlab完成伤真-歌曲“我的祖国”片段(30S左右)，在时域和频域上表示不同的状态及关系。要求：至少三个频率，截取音频时域图和频域图。

在MATLAB中，可以使用内置的音频处理工具箱来完成对《我的祖国》片段的分析。以下是简要步骤：

1. **导入音频**:
使用`audioread`函数读取音频文件，例如 `[audio, Fs] = audioread('song.mp3');` 这里假设歌曲名为 'song.mp3'，Fs 是采样率。

2. **确定时间范围和帧大小**:
选择30秒的音频片段，`time_start = 0; % 开始位置 (秒)`
`time_end = 30; % 结束位置 (秒)`
`frame_size = round(Fs * time_resolution); % 根据需要设置帧大小`

3. **截取音频段**:
`start_index = round(time_start*Fs);`
`end_index = start_index + frame_size - 1;`
`audio_segment = audio(start_index:end_index);`

4. **时域分析（幅度谱）**:
对音频段进行短时傅立叶变换（STFT），生成时域波形图。例如：

   [t, y] = stft(audio_segment, 'Window', hann(frame_size), 'OverlapLength', frame_size/2);
   plot(t, abs(y(:,1))) % 取第一个频率成分

5. **频域分析（功率谱密度）**:
计算频率响应，如使用`abs`获取功率谱：

   P = abs(y).^2; % 功率谱密度
   freq = linspace(0, Fs/2, size(P, 2)); % 频率轴
   subplot(2,1,2)
   plot(freq, P(:,1)) % 再次取第一个频率成分

6. **选取特定频率点**:
为了显示至少三个频率，你可以选择一些显著的频率点，例如主频峰及其附近的频率：

   selected_freqs = [freq(find(P(:,1) == max(P(:,1)))), freq(find(P(:,2) == max(P(:,2)))), ...]; % 可能需要调整找到的频率

7. **绘制多条线**:
在频域图上画出这些选定频率对应的线：

   hold on
   for i = 1:length(selected_freqs)
       plot(freq, P(:,i), 'LineWidth', 2)
       text(selected_freqs(i), max(P(:,i)), num2str(selected_freqs(i)) + ' Hz');
   end
   hold off

8. **保存图像**:
`saveas(gcf, 'my_song_analysis.png')` 保存时域和频域图。