数字信号处理3个作业-----作业1带阻滤波器过滤1000HZ噪声

1. 题目

       设计一个带阻滤波器，滤除受到1000Hz正弦噪声污染的音频信号。加入正弦噪声，幅度0.1，频率1000Hz。

2.分析

      1000HZ正弦噪声的频谱图是一根(对称两根)谱线,带阻滤波器的阻带覆盖该谱线,然后变换为时域形成一个序列,将该序列与音频信号序列进行卷积即可滤去1000HZ正弦噪声,但也会滤去原信号1000HZ左右频率信号。
       使用切比雪夫1型滤波器来设计滤波器，重点是阻带通带频率的调整。然后使用切比雪夫1型滤波器获得的两个向量参数对音频信号进行卷积滤波，使用移动平均滤波器来对含噪数据进行平滑处理。此示例使用 filter 函数计算沿数据向量的平均值。移动平均值滤波器沿数据移动长度为 windowSize 的窗口，并计算每个窗口中包含的数据的平均值。以下差分方程定义向量 x 的移动平均值滤波器：

3.步骤

1.读取音频文件，格式为wav:
[audio,fs]=audioread(‘audio.wav’);
%取一个通道的信号
audio = audio(:,1); %双通道变单通道

2.产生一个正弦信号A=0.1，f=1000Hz:
noise=0.1sin(10002pi/fstt’);%加噪声 1000HZ
%tt’是(1:n)’，有公式