应用MATLAB平台给语音信号叠加噪声，噪声类型是：(1）白噪声（2）单频噪声（3）多频噪声或者其它噪声，画出时域和频域图形，并给出听觉上的比较。用MATLAB代码实现

在MATLAB中，你可以通过以下几个步骤来处理语音信号并叠加不同类型的噪声：

1. **导入语音信号**:

   [speech, fs] = audioread('your_speech_file.wav'); % 替换为实际语音文件路径和采样频率

2. **创建噪声样本**:
- **白噪声**: 使用`awgn()`函数生成等方差的噪声，假设SNR（信噪比）已知：

     snr = 10; % 设定信噪比值
     noise_white = awgn(speech, snr, 'measured');

- **单频噪声**: 可以使用`sin()`函数生成正弦波噪声：

     freq = 1000; % 设定噪声频率
     noise_mono = speech * sin(2*pi*freq*(0:length(speech)-1)/fs);

- **多频噪声**: 可以使用傅里叶变换添加随机频率成分：

     noise_multi = speech + randn(size(speech)) * abs(fft(speech))/length(speech); % 随机相位

3. **叠加噪声到语音**:

   noisy_speech_white = speech + noise_white;
   noisy_speech_mono = speech + noise_mono;
   noisy_speech_multi = speech + noise_multi;

4. **绘制时域和频域图**:
- 时域图可以使用`plot()`或` stem() `，比如` plot(t, noisy_speech)`。
- 频域图使用` spectrogram()` 或 `fft() `，如 `specgram(noisy_speech)`。

5. **听觉比较**:
创建音频播放器来分别播放原始语音、叠加了不同噪声类型的语音，可以用`sound`或`audioplayer`函数。

以下是完整的例子代码片段：

% ... (之前的代码)

t = (0:length(speech)-1)./fs;

% 白噪声示例
noisy_speech_white = speech + noise_white;
sound(noisy_speech_white, fs);

% 单频噪声示例
noisy_speech_mono = speech + noise_mono;
sound(noisy_speech_mono, fs);

% 频谱分析
figure, specgram(noisy_speech_white, 256, [], [], fs);
figure, specgram(noisy_speech_mono, 256, [], [], fs);
figure, specgram(noisy_speech_multi, 256, [], [], fs);

% 听觉比较
[~, ~, playhandle1] = sound(speech, fs);
pause(.1);
stop(playhandle1);
[~, ~, playhandle2] = sound(noisy_speech_white, fs);
pause(.1);
stop(playhandle2);
[~, ~, playhandle3] = sound(noisy_speech_mono, fs);