- 1、本文档共16页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实用标准文案
精彩文档
实验报告
课程名称:数字信号处理
实 验 八: FIR数字滤波器设计及应用
班 级:通信1403
学生姓名:强亚倩
学 号:1141210319
指导教师:范杰清
华北电力大学(北京)
一、实验目的
加深理解FIR数字滤波器的时域特性和频域特性,掌握FIR数字滤波器的设计原理与设计方法,以及FIR数字滤波器的应用。
二、实验原理
FIR数字滤波器可以设计成具有线性相位,在数据通信、图像处理、语音信号处理等实际应用领域得到广泛应用。
M阶FIR数字滤波器的系统函数为:
FIR数字滤波器的单位脉冲响应h[k]是长度为M+1的有限长因果序列。当满足对称条件时,该FIR数字滤波器具有线性相位。FIR数字滤波器设计方法主要有窗口法、频率取样法及优化设计法。
MATLAB中提供的常用FIR数字滤波器设计函数有:
fir1 窗函数法设计FIR数字滤波器(低通、高通、带通、
带阻、多频带滤波器)
fir2 频率取样法设计FIR数字滤波器:任意频率响应
firls FIR数字滤波器设计:指定频率响应
firrcos 升余弦型 FIR 数字滤波器设计
intfilt 内插FIR 数字滤波器设计
kaiserord 凯塞(Kaiser)窗函数设计法的阶数估计
firpm Parks-McClellan算法实现FIR数字滤波器优化设计
firpmord Parks-McClellan 数字滤波器的阶数选择
cremez 复系数非线性相位FIR等波纹滤波器设计
1、 窗口法设计FIR数字滤波器
fir1函数可以很容易地实现FIR数字滤波器窗口法设计。
可设计低通、高通、带通、带阻滤波器、多频带滤波器。
b = fir1(M, Wc)
b = fir1(M, Wc, ftype)
b = fir1(M, Wc, window)
b = fir1(M, Wc, ftype, window)
b = fir1(M,Wc,‘ ftype,window)
输出参数:b为FIR数字滤波器的M+1个系数构成的矩阵
(即系统的单位脉冲响应)
输入参数:M为FIR数字滤波器的阶数。 Wc为3dB截频:0 Wc 1, 1 对应数字频率。 ftype指定滤波器类型,当ftype为:
’high’指定一个截频为Wc的高通滤波器;
’stop’指定一个带阻滤波器,其阻带截止频率为Wc=[w1,w2];
’DC-0’ 在多频带滤波器中,使第一个频带0ww1为阻带;
’DC-1’ 在多频带滤波器中,使第一个频带0ww1为通带。
window 指定窗函数,若不指定,默认为哈明窗。
2. 频率取样法设计FIR滤波器
fir2函数可以实现FIR数字滤波器的频率取样法设计。
可设计任意形状频率响应的滤波器。格式如下:
b = fir2(M, f, m)
b = fir2(M, f, m, window)
输出参数:b为FIR数字滤波器的M+1个系数构成的矩阵。
输入参数:M为滤波器的阶数。
f指定归一化的各频带边界频率,从0到1递增, 1对应于fsam/2,即数字频率W=p。
m指定各频带边界频率处的幅度响应, 因此f和m的长度相等,即length(f)=length(a)。
window指定窗函数,若不指定,默认为哈明窗。
三、 实验内容
1.分别使用矩形窗、汉明窗、汉宁窗设计一个阶数 M=9的FIR数字低通滤波器,截频为
(1)画出各种方法设计的数字滤波器的单位脉冲响应。
(2)画出它们的幅频响应,并比较各滤波器的通带纹波和阻带纹波,有何结论?
(3)若输入为 计算各滤波器的输出并画出其波形.
编程
b1=fir1(9,1/3,boxcar(10));
[H1,w]=freqz(b1,1,512);
H1_db=20*log10(abs(H1));
b2=fir1(9,1/3,hamming(10));
[H2,w]=freqz(b2,1,512);
H2_db=20*log10(abs(H2));
b3=fir1(9,1/3,hanning(10));
[H3,w]=freqz(b3,1,512);
H3_db=20*log10(abs(H3));
subplot(4,1,1);
您可能关注的文档
最近下载
- DHT11数字温湿度传感器毕业设计(附程序清单).doc
- 八年级期末复习必背(课件)-2023-2024学年八年级地理上学期期末考点大串讲(人教版).pptx
- 预防接种异常反应的处理课件.ppt VIP
- 苏教版小学数学五年级上册第一单元同步练习.docx VIP
- 数学苏教版五年级(上册)钉子板上的多边形(课件) (1).pptx
- 《个体防护装备安全管理规范AQ 6111—2023》知识培训.pptx
- IACS Rec.47船舶建造及修理质量标准IACS(中文).pdf
- 中国历史练习部分第三册参考答案.pdf
- 2022年6月浙江省普通高校招生选考物理试卷和答案.doc VIP
- 新能源汽车检测维修考证复习题库.docx
文档评论(0)