把一段音频导入Matlab，画出时域波形，频谱和时频谱

在MATLAB中，你可以通过以下几个步骤将音频导入并分析其时域、频域和时频特征：

1. **导入音频**：
使用`audioread`函数读取音频文件，它会返回一个向量，包含音频的样本值以及采样率信息。例如：

   [audioData, sampleRate] = audioread('your_audio_file.wav');

2. **绘制时域波形**：
要展示音频的时间序列，可以使用`plot`函数，传入时间轴（一般基于采样率计算）和音频数据：

   timeVector = (0:length(audioData)-1) / sampleRate; % 时间轴
   plot(timeVector, audioData);
   xlabel('Time (seconds)');
   ylabel('Amplitude');
   title('Audio Signal in Time Domain');

3. **绘制频谱**：
可以使用`fft`（快速傅立叶变换）计算频谱，然后用`abs`获取幅度，最后通过`plot`显示结果：

   spectrum = abs(fft(audioData)); % 计算离散傅立叶变换
   freqs = (0:sampleRate/length(audioData):sampleRate-1)' * sampleRate/(length(audioData)/2); % 频率范围
   plot(freqs, 2/sampleRate * spectrum(1:end/2+1)); % 去掉负频率部分并调整幅值
   xlabel('Frequency (Hz)');
   ylabel('Magnitude');
   title('Audio Spectrum');

4. **绘制时频谱（STFT）**：
MATLAB中的`spectrogram`函数非常适合这一步，它可以提供短时傅立叶变换的可视化：

   windowSize = round(sampleRate*0.05); % 5% of the total duration as a short-time window size
   overlap = round(windowSize*0.95); % 95% overlap for adjacent frames
   spectrogram(audioData, hann(windowSize), windowSize, overlap, sampleRate);
   colorbar;
   xlabel('Time (frames)');
   ylabel('Frequency (Hz)');
   title('Short-Time Fourier Transform');