MATLAB课程设计要点.docx
- 文档编号:24684021
- 上传时间:2023-05-31
- 格式:DOCX
- 页数:15
- 大小:60.98KB
MATLAB课程设计要点.docx
《MATLAB课程设计要点.docx》由会员分享,可在线阅读,更多相关《MATLAB课程设计要点.docx(15页珍藏版)》请在冰豆网上搜索。
MATLAB课程设计要点
郑州轻工业学院
课程设计任务书
题目基于双线性Z变换法的IIR数字带阻滤波器设计
专业、班级电子信息工程学号姓名
主要内容、基本要求、主要参考资料等:
主要内容:
利用MATLAB软件读取一段含有噪声的.wav格式的语音信号,然后基于FFT对该信号进行频谱分析;基于含噪语音信号的频谱确定滤波器的参数,利用双线性Z变换法设计IIR数字带阻滤波器,并利用所设计的滤波器对信号进行滤波处理。
比较滤波前后语音信号的时域波形及频谱,分析滤波前后的语音变化。
基本要求:
1、基于含噪语音信号的频谱确定滤波器的参数;
2、基于双线性Z变换法分别设计一个巴特沃斯IIR数字带阻滤波器和切比雪夫IIR数字带阻滤波器;
3、掌握利用wavread函数读取、播放.wav格式语音信号的方法;
4、对语音信号进行滤波,绘制滤波前后信号的时域波形及频谱;
5、回放语音信号,分析滤波前后的语音变化。
主要参考资料:
1、从玉良.数字信号处理原理及其MATLAB实现[M].北京:
电子工业出版社.2009.7
2、胡广书.数字信号处理理论、算法与实现[M].北京:
清华大学出版社.2003,8
完成期限:
2010.6.22—2010.6.26
指导教师签名:
课程负责人签名:
2013年6月22日
郑州轻工业学院
课程设计说明书
题目:
基于双线性Z变换法的IIR数字带阻滤波器设计
姓名:
院(系):
电子信息工程
专业班级:
电信工程
学号:
指导教师:
蔡超峰
成绩:
时间:
2015年6月29日至2015年7月3日
目录
1引言2
2课程设计要求4
2.1课程设计题目4
2.2课程设计目的4
2.3设计要求4
2.4数字滤波器的基本概念5
2.5数字滤波器设计的基本步骤5
2.5.1确定指标5
2.5.2模型逼近5
2.5.3实现性能分析和计算机仿真5
3IIR数滤波器设计的原理与方法字6
3.1IIR数字滤波器设计的原理6
3.2IIR数字滤波器设计的基本方法7
4IIR带阻滤波器的MATLAB设计8
4.1IIR带阻滤波器的设计流程8
5IIR带阻滤波器的仿真结果及波形12
6IIR带阻滤波器的仿真结果分析14
总结16
参考文献18
1.引言
在现代通信系统中,由于信号中经常混有各种复杂成分,所以很多信号的处理和分析都是基于滤波器而进行的。
数字滤波器可以理解为是一个计算程序或算法,将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列,并在转化过程中,使信号按预定的形式变化。
数字滤波器有多种分类,根据数字滤波器冲激响应的时域特征,可将数字滤波器分为两种,即无限长冲激响应(IIR)滤波器和有限长冲激响应(FIR)滤波器。
但是,传统的数字滤波器的设计使用繁琐的公式计算,改变参数后需要重新计算,从而在设计滤波器尤其是高阶滤波器时工作量很大。
利用MATLAB信号处理箱(SignalProcessingToolbox)可以快速有效地实现数字滤波器的设计与仿真。
本文设计一个IIR数字带阻滤波器。
IIR数字滤波器具有无限宽的冲激响应,与模拟滤波器相匹配,所以IIR滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。
其设计方法主要有经典设计法、直接设计法和最大平滑滤波器设计法。
在对滤波器实际设计时,整个过程的运算量是很大的。
设计完成后对已设计的滤波器的频率响应要进行校核,要得到幅频、相频响应特性,运算量也是很大的。
平时所要设计的数字滤波器,阶数和类型并不一定是完全给定的,很多时候要根据设计要求和滤波效果不断地调整,以达到设计的最优化。
在这种情况下,滤波器设计就要进行大量复杂的运算,单纯的靠公式计算和编制简单的程序很难在短时间内完成。
利用MATLAB强大的计算功能进行计算机辅助设计,可以快速有效地设计数字滤波器,大大地简化了计算量。
滤波器的功能是让一定频率范围内的信号通过,而将此频率范围之外的信号加以抑制或使其急剧衰减。
当干扰信号与有用信号不在同一频率范围之内,可使用滤波器有效的抑制干扰。
用LC网络组成的无源滤波器在低频范围内有体积重量大,价格昂贵和衰减大等缺点,而用集成运放和RC网络组成的有源滤波器则比较适用于低频,此外,它还具有一定的增益,且因输入与输出之间有良好的隔离而便于级联。
由于大多数反映生理信息的光电信号具有频率低、幅度小、易受干扰等特点,因而RC有源滤波器普遍应用于光电弱信号检测电路中
数字滤波器的实用型式很多,大略可分为有限冲激响应型和无限冲激响应型两类,可用硬件和软件两种方式实现。
在硬件实现方式中,它由加法器、乘法器等单元所组成,这与电阻器、电感器和电容器所构成的模拟滤波器完全不同。
数字信号处理系统很容易用数字集成电路制成,显示出体积小、稳定性高、可程控等优点。
数字滤波器也可以用软件实现。
软件实现方法是借助于通用数字计算机按滤波器的设计算法编出程序进行数字滤波计算。
带阻滤波器是指能通过某一频率范围内的频率分量、但将其他范围的频率分量衰减到极低水平的滤波器,与带阻滤波器的概念相对。
一个模拟带阻滤波器的例子是电阻-电感-电容电路(RLCcircuit)。
这些滤波器也可以用低通滤波器同高通滤波器组合来产生.
2.课程设计要求
2.1课程设计题目
数字滤波系统的MATLAB实现与分析
——双线性变换法设计数字带阻IIR滤波器
2.2课程设计目的
通过本次综合设计,运用已学的课程知识,根据题目要求进行软件仿真设计和调试,对《数字信号处理》课程中涉及的滤波器原理和设计方法有一定的感性认识和实践操作能力,从而加深对本课程知识点的理解,使学生应用知识能力、设计能力、调试能力以及报告撰写能力等方面有显著提高。
学会MATLAB的使用,掌握MATLAB的程序设计方法;掌握在Windows环境下语音信号采集的方法;掌握数字信号处理的基本概念、基本理论和基本方法;掌握MATLAB设计FIR和IIR数字滤波器的方法;学会用MATLAB对信号进行分析和处理。
2.3设计要求
本设计的主要内容包括掌握数字信号处理的基础知识和利用MATLAB软件进行编程,利用几种常见的数字滤波器对信号进行处理,并对过滤结构进行分析。
具体要求包括以下几个方面:
(1)滤波器与数字信号处理的理论知识学习和分析
通过查找、学习相关文献资料,掌握滤波器与数字信号处理的基本理论知识,掌握滤波器设计的方法和基础知识,掌握滤波器过滤信号的基本过程和原理,并对相应信号进行时域和频域分析,理解信号的频谱变化的特点。
(2)利用MATLAB对上述过程进行仿真
利用MATLAB软件进行相应的设计,实现滤波器对信号的过滤处理。
(3)结果分析
利用滤波器处理相应信号,对过滤结果进行分析,对比处理前后的信号频率的衰变,说明信号过滤的情况好坏,并分析该类型的滤波器的优缺点和适用场合。
2.4数字滤波器的基本概念
滤波器的种类很多,从功能上可分为低通、高通、带通和带阻滤波器,每一种又有模拟滤波器和数字滤波器两种形式。
如果滤波器的输人和输出都是离散时间信号,则该滤波器的冲击响应也必然是离散的,这种滤波器称之为数字滤波器。
数字滤波器是一种用来过滤时间离散信号的数字系统,通过对抽样数据进行数学处理来达到频域滤波的目的。
数字滤波器也是具有一定传输选择特性的数字信号处理装置,其输入、输出均为数字信号,实质上是一个由有限精度算法实现的线性时不变离散系统。
它的基木工作原理是利用离散系统特性对系统输人信号进行加工和变换,改变输入序列的频谱或信号波形,让有用频率的信号分量通过,抑制无用的信号分量输出。
2.5数字滤波器设计的基本步骤
数字滤波器的设计要经过如下三个步骤:
2.5.1确定指标
在设计一个滤波器之前,必须首先根据工程实际的需要确定滤波器的技术指标。
在很多实际应用中,数字滤波器常常被用来实现选频操作。
因此,指标的形式一般在频域中给出幅度和相位响应。
2.5.2模型逼近
一旦确定了技术指标,就可利用已学习过的基本原理和关系式,提出一个滤波器模型来逼近给定的指标体系。
这是滤波器设计所要研究的主要问题。
2.5.3实现性能分析和计算机仿真
以上两步的结果得到的滤波器,通常是以差分方程、系统函数或脉冲响应描述的。
根据这个描述就可以分析其频率特性和相位特性、以验证设计结果是否满足指标要求;或者利用计算机仿真实现设计的滤波器,再分析滤波结果来判断。
3.IIR数字滤波器设计的原理与方法
3.1IIR数字滤波器设计的原理
IIR数字滤波器具有无限持续时间冲激响应,需要用递归模型来实现,其系统函数为:
(2.1)
设计IIR滤波器的任务就是寻求一个因果、物理上可实现的系统函数H(z),使其频率响应
满足所希望得到的频域指标,即符合给定的通带截止频率、阻带截止、通带衰减和阻带衰减.。
利用冲激响应不变法设计数字滤波器时可能会导致频域混叠现象,为了克服这一问题,需要找到由s平面到z平面的另外的映射关系,这种关系应保证:
1)s平面的整个jΩ轴仅映射为z平面单位圆上的一周;
2)若G(s)是稳定的,由G(s)映射得到的H(z)也应该是稳定的;
3)这种映射是可逆的,既能由G(s)得到H(z),也能由H(z)得到G(s);
4)如果G(j0)=1,那么
=1。
双线性Z变换满足以上4个条件的映射关系,其变换公式为
(2.2)
双线性Z变换的基本思路:
首先将整个S平面压缩到一条从-π/Ts变换到2π/Ts的横带里,然后通过标准的变换关系
将横带变换到整个Z平面上去,这样就得到了S平面与Z平面间的一一对应的单值关系。
3.2IIR数字滤波器设计的基本方法
IIR数字滤波器的设计方法有两类,一类是借助于模拟滤波器的设计方法设计出模拟滤波器,利用冲激响应不变法或双线性变换法转换成数字滤波器,再用硬件或软件实现;另一类是直接在频域或时域中进行设计,设计时需要计算机作辅助工具。
随着MATLAB软件尤其是MATLAB的信号处理工作箱的不断完善,不仅数字滤波器的计算机辅助设计有了可能,而且还可以使设计达到最优化。
IIR数字滤波器设计的基本步骤如下:
(1)根据任务,确定性能指标。
在设计一个滤波器之前,首先根据工程实际的需要确定滤波器的技术指标如:
边界频率:
ωp,ωs,ωc;阻带最小衰减As和通带最大衰减Rp;
(2)将数字滤波器的技术指标转换成模拟滤波器指标。
利用冲激响应不变法与双线性变换法进行频率间的转换,主要是边界频率Wp与Ws的转换。
(3)用模拟滤波器设计方法得到模拟滤波器的传输函数Ha(s);可借助巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等,这些滤波器都有严格的设计公式、现成曲线和图表供设计人员使用。
4.IIR带通滤波器的MATLAB设计
4.1IIR带阻滤波器的设计流程
本文设计的IIR带阻滤波器是从低通变换过来的,利用的是双线性变换以及切比雪夫II滤波器的原型,其具体的设计流程为上图所示。
首先根据题目要求确定带阻滤波器的技术指标,先要进行频率的预畸变,并且归一化频率,再设计出切比雪夫II模拟低通滤波器,并求出其阶数等相关参数。
其次利用双线性变换法设计数字带通滤波器,,再调用函数进行双线性变换,并求出分子、分母的系数向量。
最后通过画图求出其幅频响应、相频响应、幅度特性曲线与零极点,并画出波形图。
最后进行验证,看所设计的滤波器能否达到要求的指标,若能达到,则说明该滤波器设计符合要求。
3.2IIR带阻滤波器的设计步骤
(1)根据设计流程,首先确定所要设计的数字带阻滤波器的相关指标:
① 通带截止频率wp1=0.4π,wp2=0.6π,通带最大衰减Rp=2dB;
② 阻带截止频率ws1=0.2π和ws2=0.8π,阻带最小衰减Rs=30dB;
③ 取样间隔T=0.1s。
其实现程序如下(程序中pi代表π):
clear;clc;closeall;
[S,Fs,Bit]=wavread('C:
\Users\yyy\Desktop\cc\SHE.wav');%读取音频信号
[S1,Fs,Bit]=wavread('C:
\Users\yyy\Desktop\cc\Bandstop.wav');%读取音频信号(含有噪声)
sound(S,Fs);
sound(S1,Fs);
X=abs(fft(S));X1=abs(fft(S1));
xlabel=0:
Fs/(length(S)-1):
Fs/2;
subplot(2,1,1);plot(xlabel,X(1:
150001));title('原始信号频谱图');%分析音频信号,找出干扰频率段为8000Hz
subplot(2,1,2);plot(xlabel,X1(1:
150001));title('加噪信号频谱图');
%对于带阻滤波器:
fp=5000;
fs=6000;
Rp=3;
Rs=40;
%求数字频率
wp=2*pi*fp/Fs;
ws=2*pi*fs/Fs;
%预畸变处理
Wp=2*Fs*tan(wp/2);
Ws=2*Fs*tan(ws/2);
%开始设计模拟低通滤波器
[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s');
[B,A]=butter(N,Wn,'s');
[b,a]=bilinear(B,A,Fs);
[h,w]=freqz(b,a,1000,'whole');
h=h(1:
501);
w=w(1:
501);
mag=abs(h);
db=10*log10((mag+eps)/max(mag));
figure;subplot(1,2,1);plot(w/pi,mag);
title('幅频响应');
axis([0,1,0,1.1]);
subplot(1,2,2);plot(w/pi,db);
title('幅频响应(dB)');
axis([0,1,-40,1]);
set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1]);
set(gca,'YTickMode','manual','YTick',[-40,-20,-3,0]);grid;
S2=filter(b,a,S1);
sound(S2,Fs);
X2=abs(fft(S2));
figure;plot(xlabel,X2(1:
150001));title('滤波以后的频谱')
5.IIR带阻滤波器频谱
6.IIR带阻滤波器的仿真结果分析
根据上述的仿真波形,可以看出:
(1)由ChebyshevI型、ChebyshevII型低通原型变换成带通模型的幅频特性可以用分贝形式表示。
对于ChebyshevI型的幅频响应在通带内为等纹波衰减、阻带内为单调减小的,且通带(归一化)截止频率在波形上显示为:
0.4与0.6,阻带(归一化)截止频率在波形上显示为:
0.2与0.8,与设计要求基本一致。
对于ChebyshevII型的幅频响应来说,在通带内单调减小、阻带内为等纹波衰减的,而通带、阻带(归一化)截止频率为:
0.4与0.6,0.2与0.8,均满足设计要求。
(2)此次设计的IIR带阻滤波器的阶数可由函数[N,OmegaC]=cheb1ord(Wp,Ws,Rp,Rs,'s')
或[N,OmegaC]=cheb2ord(Wp,Ws,Rp,Rs,'s')求得,N=3。
则从ChebyshevI型带阻滤波器的幅频响应曲线上看,当N=3时,在
(归一化后为0--1)范围内通带波动只有两个波谷,阻带则单调减小,这与理论内容一致。
对于ChebyshevII型带阻滤波器来说,其幅频响应曲线的通带波动在
(归一化后为0--1)范围内阻带波动只有两个波峰,通带则单调减小,这与理论内容一致。
因此,对于滤波器阶数这一指标来说,满足要求。
(3)对于分贝化后的幅度特性,ChebyshevI型带阻滤波器的通带最大衰减大约为2dB,阻带衰减最小为30dB,而ChebyshevII型带阻滤波器的阻带最小衰减大约为30dB,通带最小衰减大约为2dB,这与带阻滤波器的技术指标相同。
因此,在衰减性能上,所设计的滤波器达到了要求。
(4)对于相频响应来说,ChebyshevI型带阻滤波器在通带内的相频响应曲线接近为一条平滑曲线,在阻带内则存在衰减、畸变。
ChebyshevII型带阻滤波器在通带截止频率以内的相频响应曲线几乎为一条直线,但在通带截止频率处有很大角度的转折,使相频特性在通带内产生了的畸变,而在阻带内存在很大的畸变。
因此,对于双线性变换法,其相位特性得不到满足,必须用其它方法加以校正或采用其它方法重新设计滤波器,使之满足线性相位。
(5)群延迟是衡量一个滤波器或是整个系统性能指标的一个重要参数。
对于利用这两种低通滤波器原型模型设计的这两种带通滤波器,带阻滤波器的延迟是同等带宽的低通滤波器延迟的两倍。
这个结果是由低通特性向带通特性变换引起的,n阶低通滤波器传递函数总是变换为2n阶带通滤波器传递函数,其滤波器的群延迟特性的中心频率与理论上的很接近。
总结
此次Matlab课程设计是在一定的理论基础之上进行的,在先修课程《信号与系统》与《数字信号处理》中,大量有Matlab设计方面的知识,所以做起来还比较容易,而且经过了自己的亲身实践,学到了许多实践方面的知识。
首先,在信号滤波系统中,有时因为模拟滤波器阶数太高,硬件占用空间太大为某些仪器的实现设置了障碍,而对于一些窄带情况下的低通滤波器用模拟手段往往很难实现。
在这些情况下,数字滤波器将会是一个很好的解决办法。
MATLAB信号处理工具箱提供了丰富而简便的设计、实现FIR和IIR的方法,使原来繁琐的程序设计简化成函数的调用,特别是滤波器的表达方式和滤波器之间的相互转换显得十分简便。
其次,IIR数字滤波器的设计和模拟滤波器的设计有着紧密的关系。
通常要先设计出适当的模拟滤波器,再通过一定的频带变换把它转换成为所需的数字IIR滤波器。
此外,任何数字信号处理系统中也还不可避免地用到模拟滤波器,比如A/D变换器前的抗混叠滤波器和D/A变换器后的平滑滤波器,因此模拟滤波器设计也是很重要的。
最后,在比较设计滤波器的方法上应该明确其技术指标以及某些参数的实际意义。
比如本文用双线性变换法设计数字带阻滤波器时,必须先将频率归一化,并且进行频率预畸变,然后设计模拟滤波器,再利用频率变换法将模拟低通变为模拟带通,最后经过双线性变换法将模拟带通变换为数字带通。
如果不进行频率预畸变,那么设计出来的带通滤波器的幅频特性与相频特性将会产生很严重的畸变(如图5、图6所示,为设计时没有进行频率预畸变的特性),使设计的结果不满足给定的要求,在实际中会造成很严重的危害。
在做本次课程设计的过程中,我深深地感受到了自己所学到知识的有限,明白了只学好课本上的知识是不够的,要通过图书馆和互联网等各种渠道来扩充自己的知识。
在实验过程中我们曾经遇到过问题。
但是从中我们学习到了如何对待遇到的困难,进一步培养了我们一丝不苟的科学态度和不厌其烦的耐心。
所有的这些心得会对我以后的学习和工作有帮助作用,忠心感谢学校给我们提供这次实验机会。
参考文献
1、从玉良.数字信号处理原理及其MATLAB实现[M].北京:
电子工业出版社.2009.7
2、胡广书.数字信号处理理论、算法与实现[M].北京:
清华大学出版社.2003,8
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 课程设计 要点