数字信号处理课程设计 数字滤波器的设计及实现.docVIP

数字滤波器的设计及实现 摘要：通过MATLAB软件的信号处理工具箱中的滤波器设计各种IIR和FIR数字滤波器，并绘制滤波器的幅频特性、相频特性。通过观察滤波器的输入、输出信号的时域波形及其频谱，建立了数字滤波的概念。 关键词：MATLAB；IIR数字滤波器；FIR数字滤波器 （一）设计目的 1、熟悉IIR数字滤波器和FIR数字滤波器的设计原理和方法； 2、学会调用MATLAB信号处理工具箱中的滤波器设计函数设计各种IIR和FIR数字滤波器，学会根据滤波要求确定滤波器指标参数； 3、掌握用IIR和FIR数字滤波器的MATLAB实现方法，并能绘制滤波器的幅频特性、相频特性； 4、通过观察滤波器的输入、输出信号的时域波形及其频谱，建立数字滤波的概念。 （二）设计要求 用MATLAB软件设计IIR数字滤波器和FIR数字滤波器，并绘制滤波器的幅频特性、相频特性。 （三）设计原理 数字滤波器是将输入数字序列通过一定的运算后转变为输出数字序列的数字信号处理器。数字滤波器的输入、输出均为数字信号, 通过一定运算系改变输入信号所含频率成分的相对例或消除某些频率成分。与模拟滤波器相比, 数字滤波器的主要优点是：（1）精度和稳定性高; (2)系统函数容易改变, 因而灵活性高; （3）不存在阻抗匹配问题; ( 4)便于大规模集成; ( 5)可以实现多维滤波。它不仅能实现模拟处理的大部分功能, 而且还能完成模拟处理由于成本、可靠性等原因而无法具体实现的功能。 所谓抑制载波单频调制信号，就是两个正弦信号相乘，它有2个频率成分：和频+，差频-，这两个频率成分关于载波频率对称。所以，1路抑制载波单频调幅信号的频谱图是关于载波频率对称的两根谱线。显然，当调制频率和（或）载波频率不同时，可以得到包含不同频率成分的单频调幅信号。 （四）设计内容 1、调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st，观察st的时域波形和幅频特性曲线； 根据题目的要求编写代码如下： function st=mstg %产生信号序列st，并显示st的时域波形和频谱 %st=mstg返回三路调幅信号相加形成的混合信号，长度N=800 N=800; %信号长度N为800 Fs=10000;T=1/Fs;Tp=N*T; %采样频率Fs=10kHz，Tp为采样时间 t=0:T:(N-1)*T;k=0:N-1;f=k/Tp; fc1=Fs/10; %第1路调幅信号载波频率fc1=1000Hz fm1=fc1/10; %第1路调幅信号的调制信号频率fm1=100Hz fc2=Fs/20; %第2路调幅信号载波频率fc2=500Hz fm2=fc2/10; %第2路调幅信号的调制信号频率fm2=50Hz fc3=Fs/40; %第3路调幅信号载波频率fc3=250Hz fm3=fc3/10; %第3路调幅信号的调制信号频率fm3=25Hz xt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); %产生第1路调幅信号 xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); %产生第2路调幅信号 xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); %产生第3路调幅信号 st=xt1+xt2+xt3; %三路信号相加，得到复合信号 fxt=fft(st,N); %计算信号st的频谱 %以下为绘图命令 subplot(2,1,1); plot(t,st);grid;xlabel(t/s);ylabel(s(t)); axis([0,Tp,min(st),max(st)]);title((a)s(t)的波形) subplot(2,1,2); stem(f,abs(fxt)/max(abs(fxt)),.);grid;title((b)s(t)的频谱) axis([0,Fs/8,0,1.2]); xlabel(f/Hz);ylabel(幅度); 运行的波形如下： 图一三路调幅信号st的时域波形和幅频特性曲线 通过观察st的幅频特性曲线，分别确定可以分离st中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、带通滤波器、高通滤波器）的通带截止频率和阻带截止频率；假定要求滤波器的通带最大衰减为0.1dB，阻带最小衰减为60dB，编程调用MATLAB滤波器设计函数分别设计这三个数字滤波器，并绘图显示其幅频特性曲线。 根据题目编写代码如下： % 低通 巴特沃