专业知识3期精要.doc

音乐信号滤波去噪 ——使用双线性变换法设计的并联型椭圆滤波器 学生姓名 指导老师： 摘 要 本文简单、直观地介绍了椭圆低通滤波器的基本理论和设计思想，阐述了设计椭圆低通滤波器的具体步骤，并对仿真结果进行分析和处理。本课程设计主要内容是利用双线性变换法设计一个椭圆IIR滤波器，对一段音乐信号（含燥声）进行滤波去噪处理，对比滤波前后的波形和频谱图分析滤波器的性能。本课程设计仿真平台为MATLAB。最后。可知。（1）滤波器指标必须符合工程实际。（2）设计完后应检查其频率响应曲线是否满足指标。 （3）处理结果和分析结论应该一致，而且应符合理论。 （4）独立完成课程设计并按要求编写课程设计报告书。IIR系统。 无限长单位冲激响应（IIR）h(n)是无限长； 系统函数H(z)在有限z平面（0|z|∞）； 结构上存在着输出到输入的反馈，也就是结构上是递归型的。 IIR滤波器采用递归型结构，即结构上带有反馈环路。H（z）可以直接Ⅰ型、Ⅱ型、级联型、并联型四种，都具有反馈回路。IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等，这些典型的滤波器各有特点。有现成的设计数据或图表可查在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。2.3 双线性变换法 脉冲响应不变法的主要缺点是频谱交叠产生的混淆，这是从S平面到Z平面的标准变换z＝esT的多值对应关系导致的,为了克服这一缺点，设想变换分为两步：第一步：将整个S平面压缩到S1平面的一条横带里； 第二步：通过标准变换关系将此横带变换到整个Z平面上去。由此建立S平面与Z平面一一对应的单值关系，如图2-1所示。 图 2-2 双线性变换法中的复平面映射 为了将s平面的jΩ轴压缩到s1平面jΩ轴上的- 一段上，可通过以下的正切变换实现： 这里C是待定常数，下面会讲到用不同的方法确定C，可使模拟滤波器的频率特性与数字滤波器的频率特性在不同频率点有对应关系。 经过这样的频率变换，当Ω1由 时,Ω由 ,即映射了整个jΩ轴。 将这一关系解析延拓至整个s平面，则得到s平面 平面的映射关系： 再将s1平面通过标准变换关系映射到z平面，即令 ，通常取C=2/T ，最后得S平面与Z平面的单值映射关系： 双线性换法的主要优点是S平面与Z平面一一单值对应，S平面的虚轴(整个jΩ)对应于Z平面单位圆的一周，S平面的Ω=0处对应于Z平面的ω=0处， 对应 即数字滤波器的频率响应终 止于折迭频率处，所以双线性变换不存在混迭效应。 3设计步骤 3.1 设计流程图 音乐信号滤波去噪——使用双线性变化法设计的并联型椭圆滤波器的设计流程如图3-1所示。 图3-1 设计流程图 3.2 拉弦音乐的采集 先在windows平台下用麦克风录制一段语音，然后加入一段音乐背景。然后在Matlab软件平台下，利用函数wavread对语音信号进行采样，记住采样频率和采样点数。Wavread[x,fs,nbits]=wavread(file);功能说明：采样值放在向量x中，s表示采样频率（Hz）bits表示采样位数。 采集完成后在信号中加入一频率为2200Hz的单频噪声，设计的任务即为从含噪信号中滤除单频噪声，还原原始信号。 [x,fs,bits]=wavread(e:\201054080226.wav); % 输入参数为文件的全路径和文件名， %%输出的第一个参数是每个样本的值， %%fs是生成该波形文件时的采样率， %%bits是波形文件每样本的编码位数。 sound(x,fs,bits); % 按指定的采样率和每样本编码位数回放 N=length(x); % 计算信号x的长度 fn=2300; % 单频噪声频率，此参数可改 t=0:1/fs:(N-1)/fs; % 计算时间范围，样本数除以采样频率x x=x(:,1); y=x+0.1*sin(fn*2*pi*t); sound(y,fs,bits); %可以明显听出有尖锐的单频啸叫声 3.3语音信号的频谱分析 在MATLAB中编辑m函数，使用wavread函数读取采集的声音文件（.wav)将它赋值给某一向量，再对其进行采样后使用