实验五 基于MATLAB的IIR滤波器结构和窗函数
更新时间:2023-08-30 14:18:01 阅读量: 教育文库 文档下载
- 实验五推荐度:
- 相关推荐
dsp实验五,欢迎下载
实验五 基于MATLAB的IIR滤波器结构和窗函数
一、实验目的:
1. 熟悉IIR滤波器的基本结构和函数语言的实现方法; 2. 了解窗函数的性能和特点;
3. 学会利用MATLAB函数绘制窗函数的幅频特性。
二、实验原理:
1. IIR滤波器的结构
IIR滤波器有三种实现结构:直接型、级联型和并联型。一般来说,滤波器是用系统函数的有理函数形式(或直接型结构)描述的,本实验将学习用MATLAB函数进行各种结构形式的转换,会设计数字滤波器,并利用设计的滤波器进行信号的处理。
IIR直接型结构由两个行矢量b和a描述,并由filter函数调用实现,格式为y=filter(b,a,x)。IIR级联型结构的系统函数可以表示为有限个实系数二阶有理分式之积,可调用函数[z,p,k]=tf2zp(b,a)将系统函数以零极点增益表示,其中z表示零点,p表示极点,k表示增益,然后调用用函数sos=zp2sos(z,p,k)将其转换为二阶分式之积。IIR并级联型结构的系统函数通常表示为有限个实系数一阶或二阶有理分式之和,这种形式在MATLAB中通过调用[r,p,k] = residuez(b,a)函数来实现。
模拟滤波器的设计调用函数见课本第147页。 2. 典型的窗函数
(a)矩形窗(Rectangle Window)
w(n) RN(n)
其频率响应和幅度响应分别为:
W(e
j
)
sinN( /2)sin (/2)
e
j
N 12
WR( )
sin(N /2)sin( /2)
,
在matlab中调用w=boxcar(N)函数,其中N为窗函数的长度,而返回值w为一个N阶的向量,它的元素由窗函数的值组成。‘w=boxcar(n)’ 等价于‘w=ones(1,n)’。 (b)三角形窗(Bartlett Window)
2n
N 1,
w(n)
2n 2 ,
N 1
0 n N 12
N 12
n N 1
N 12
其频率响应为:W(e
j
2sin(N /4)2 j ) []e
Nsin( /2)
dsp实验五,欢迎下载
在matlab中调用w=triang(N)函数,N为窗函数的长度。w=triang(N-2)等价于
bartlett(N)。
(c)汉宁(Hanning)窗,又称升余弦窗
w(n)
12
[1 cos(
2n N 1
)]RN(n)
N 12
其频率响应和幅度响应分别为:
W(e
j
) {0.5WR( ) 0.25[WR( W( )e
j a
2 N 1
) WR(
2 N 1
)]}e
j()
W( ) 0.5WR( ) 0.25[WR(
2 N 1
) WR(
2 N 1
)]
在matlab中调用w=hanning(N)函数,N为窗函数的长度
(d)汉明(Hamming)窗,又称改进的升余弦窗
w(n) [0.54 0.46cos(
2n N 1
)]RN(n)
2 N 1
) WR(
2 N 1
)]
其幅度响应为:W( ) 0.54WR( ) 0.23[WR(
在matlab中调用w=hamming(N)函数,N为窗函数的长度
(e)布莱克曼(Blankman)窗,又称二阶升余弦窗
w(n) [0.42 0.5cos(
2n N 1
) 0.08cos(
4n N 1
)]RN(n)
2 N 1
W( ) 0.42WR( ) 0.25[WR(
2 N 1
) WR( 4 N 1
)]
)]
其幅度响应为:
0.04[WR(
4 N 1
) WR(
在matlab中调用w=blackman(N)函数,N为窗函数的长度
(f)凯泽(Kaiser)窗
w(n)
I0( [1 2n/(N 1)])
I0( )
2
,0 n N 1
其中:β是一个可选参数,用来选择主瓣宽度和旁瓣衰减之间的交换关系,一般说来,β越大,过渡带越宽,阻带越小衰减也越大。I0(·)是第一类修正零阶贝塞尔函数。
在matlab中调用w=kaiser(N, beta),函数N为窗函数的长度,beta为窗函数的参数。
三、实验内容:
1.实现下面系统的级联型结构和并联型结构:
dsp实验五,欢迎下载
3
H(z)
1
16z
53
1
z
1
z
23
2
z
2
13
16
。 z
3
解:级联型: clear all
b=[3,5/3,2/3]; a=[1,1/6,1/3,-1/6];
[b,a] = eqtflength(b,a); %制造b,a长度相等; [z,p,k]=tf2zp(b,a); sos=zp2sos(z,p,k); 运行结果为sos =
3 0 0 1 -1/3 0 1 5/9 2/9 1 1/2 1/2 可得级联型系统函数为
31
13z
1
1 1
5912
zz
1
2912
zz
2
H(z)
2
1
并联型:
b=[3,5/3,2/3]; a=[1,1/6,1/3,-1/6]; [b,a] = eqtflength(b,a); [r,p,k]=residuez(b,a); R1=[r(1),r(2)]
P1=[p(1),p(2)]
[b1,a1]=residuez(R1,P1,0) [b2,a2]=residuez(r(3),p(3),0) 运行结果为
r = 1/2 - 1012/1785i ,1/2 + 1012/1785i ,2 p = -1/4 + 506/765i , -1/4 - 506/765i ,1/3 k = []
b1 =1 1 0 a1 =1 1/2 1/2 b2 =2 0 a2 =1 -1/3
其中R1为共轭复数留数所构成的向量,P1为共轭复数极点所构成的向量,用函数
dsp实验五,欢迎下载
[b1,a1]=residues(R1,P1,0)可获得对应的实系数二阶分式的分子、分母多项式系数,b1,a1为有理分式的分子分母多项式系数向量。
可得并联型系统函数为
H(z)
1
213z
1
1
1 z12z
1
1
12
z
2
2 矩形窗
N=16;
w=boxcar(N); nn=0:(N-1); figure stem(nn,w) axis([0 20 0 1]) [H,W]=freqz(w,1); figure
plot(W/pi,abs(H))
四、作业:
(1)画出三角形窗、汉宁窗、汉明窗和布莱克曼窗的时域和频域图,并对比他们的主瓣宽度,得出你的结论。
dsp实验五,欢迎下载
正在阅读:
实验五 基于MATLAB的IIR滤波器结构和窗函数08-30
北矿导师 - 图文10-24
双曲线的简单几何性质(2)典型例题解析01-25
黑龙江省促进中小企业发展条例03-20
CSS-100BE低频低压减载说明书(0SF.457.014) - V1.001-28
2013年物业管理师实战模拟试题三06-03
五年级以上的读书笔记范文5篇.doc05-09
我的小猫作文400字07-16
- exercise2
- 铅锌矿详查地质设计 - 图文
- 厨余垃圾、餐厨垃圾堆肥系统设计方案
- 陈明珠开题报告
- 化工原理精选例题
- 政府形象宣传册营销案例
- 小学一至三年级语文阅读专项练习题
- 2014.民诉 期末考试 复习题
- 巅峰智业 - 做好顶层设计对建设城市的重要意义
- (三起)冀教版三年级英语上册Unit4 Lesson24练习题及答案
- 2017年实心轮胎现状及发展趋势分析(目录)
- 基于GIS的农用地定级技术研究定稿
- 2017-2022年中国医疗保健市场调查与市场前景预测报告(目录) - 图文
- 作业
- OFDM技术仿真(MATLAB代码) - 图文
- Android工程师笔试题及答案
- 生命密码联合密码
- 空间地上权若干法律问题探究
- 江苏学业水平测试《机械基础》模拟试题
- 选课走班实施方案
- 滤波器
- 函数
- 基于
- 结构
- 实验
- MATLAB
- IIR
- 汽车倒车防撞系统毕业论文
- 走近劳动保障 欢迎登陆重庆劳动保障网站
- 电力系统调度自动化第7章 EMS(3)-网络拓扑
- 2017年长春师范大学数学学院333教育综合之外国教育史教程考研题库
- 超级飞侠涂色画
- 酒店与婚庆公司战略合作提案
- 中国瑜伽教练培训产业现状调查及未来五年投资价值分析报告
- 特殊病种申请
- 小功率伺服无刷直流电机驱动器设计
- 金融学形成性考核作业(仅供参考)
- 劳动者权益保障法论文
- Unit 10 Make your choice
- 商业建筑的人流系统规划分析
- 国际品牌酒店财务部服务流程与规范
- 新材料企业对接资本市场的四大障碍
- 总体规划与概念建筑设计咨询服务合同书20101027
- OTC业务员业务市场调查报告
- 犯罪学精讲第四章
- 入党演讲稿(入党积极分子发言)
- 建筑构造形成性考核册答案