信号与系统-音乐去噪

【背景知识】

璀璨的中华文明具有非常大的多样性，仅戏曲一项，南方戏曲和北方戏曲格与内容各有侧重。南方戏曲声腔婉转，重于抒情。北方戏曲声腔激情高亢，重于叙事。现有两首受到噪声污染的戏曲diewithoutyou_noise.wav(《梁山伯与祝英台》片段)和commandtaking_noise.wav（《花木兰》片段）（见附件）。在时域观察信号波形较为杂乱，很难找到规律，但在频域观察信号的幅度谱，很容易找到规律，发现噪声频谱。利用所学知识将受到噪声污染的戏曲还原。

【研讨内容】

1. (*)分别利用Matlab函数audioread()读取这两首含噪戏曲，画出其时域波形(横坐标-秒(s))。
2. (*)分别利用Matlab函数fft计算这两首含噪戏曲的频谱，画出其幅度频谱。分析其频谱的大致分布范围，判断频谱中噪声情况。
3. (*)根据其频谱中的噪声位置，判断分别可以利用什么滤波器对其进行去噪？
4. (***)设计响应的滤波器，并画出滤波器的幅度频谱和相位谱；
5. (***)分别利用合适的滤波器对含噪戏曲进行滤波；
6.  (*)分别画出滤波器后的戏曲时域波形和幅度频谱波形，与滤波器前波形相比，有什么变化？
7. (*)利用Matlab函数sound分别试听滤波后的戏曲，并利用audiowrite函数将滤波器后戏曲数据写入新的wav文件中。

一、结果：
（一）梁祝

（二 ）穆桂英挂

1.低通滤波器

2.带阻滤波器-一遍：

3.带阻滤波器-两遍：

二、Matlab代码实现：

（一)梁祝

 （二）穆桂英挂帅

1、低通滤波器：

2.带阻滤波器-一遍

3.带阻滤波器-两遍

 代码类似1.