中南大学
《数字信号处理》实验报告
-
目录
实验一 常见离散时间信号的产生和频谱分析........ 3
一、实验目的.................................. 3
二、实验原理.................................. 3
三、 实验内容................................. 6
实验二 数字滤波器的设计.................... 12
一、实验目的................................. 12
二、实验原理................................. 12
三、实验内容................................. 16
(1)熟悉MATLAB应用环境,常用窗口的功能和使用方法;
(2)加深对常用离散时间信号的理解;
(3)掌握简单的绘图命令;
(4)掌握序列傅里叶变换的计算机实现方法,利用序列的傅里叶变换对离散信号进行频域分析。
(1)常用离散时间信号
a)单位抽样序列
如果在时间轴上延迟了k个单位,得到即:
b)单位阶跃序列
c)矩形序列
d)正弦序列
e)实指数序列
f)复指数序列
(2)离散傅里叶变换:
设连续正弦信号为
这一信号的频率为,角频率为,信号的周期为。如果对此连续周期信号进行抽样,其抽样时间间隔为T,抽样后信号以表示,则有,如果令为数字频率,满足,其中是抽样重复频率,简称抽样频率。
为了在数字计算机上观察分析各种序列的频域特性,通常对在上进行M点采样来观察分析。 对长度为N的有限长序列x(n), 有
其中
通常应取得大一些,以便观察谱的细节变化。取模可绘出幅频特性曲线。
(3)用DFT进行普分析的三种误差
三种误差:混叠现象、泄露现象、栅栏效应
a) 混叠现象
当采样频率小于两倍信号(这里指是信号)最大频率时,经过采样就会发生频谱混叠,这使得采样后的信号序列频谱不能真实地反映原信号的频谱。所以在利用DFT分析连续信号的频谱时,必须注意这一问题。避免混叠现象的唯一方法是保证采样速率足够高,使频谱交叠现象不致出现。也就是说,在确定采样频率之前,必须对信号的性质有所了解,一般在采样前,信号通过一个防混叠低通滤波器。
b) 泄漏现象
实际中的信号序列往往很长,为了方便我们往往用截短的序列来近似它们,这样可以使用较短的DFT来对信号进行频谱分析,这种截短等价于给原信号序列乘以一个矩形窗函数。
泄漏是不能与混叠完全分离开的,因为泄漏导致频谱的扩散,从而造成混叠。为了减小泄漏的影响,可以选择适当的窗函数,使频谱的扩散减到最小。
c) 栅栏效应
因为DFT是对单位圆上Z变换的均匀采样,所以他不可能将频谱视为一个连续函数。这样就产生了栅栏效应,就一定意义上看,DFT来观看频谱就好像通过一个尖桩的栅栏来观看一个图景一样,只能在离散点上看到真实频谱,这样就可能发生一些频谱的峰点或谷点被“尖桩的栅栏”所挡住,不能被我们观察到。减小栅栏效应的一个方法就是借助在原序列的末端添补一些零值,从而变动DFT的点数。这一方法实际上是人为地改变了对真实谱采样的点数和位置,相当于搬动了每一根“尖桩栅栏”的位置,从而使得频谱的峰点或者谷点暴露出来。当然,这是每根谱线所对应的频率和原来的不同了。
综上所述,DFT可以用于信号的频谱分析,但必须注意可能产生的误差,在应用过程中要尽可能减少和消除这些误差的影响。
(1)复习常用离散时间信号的有关内容
在时间上依次出现的数值序列,例如,{…,0.5,1,2,-1,0,5,…}。相邻两个数之间的时间间隔可以是相等的,也可以是不等的。在前一情况下,设时间间隔为T秒,则离散信号可用符号x(nT)来表示。在间隔T归一化为1的条件下,T可以省略,即将x(nT)表示为x(n)。x(n)既可表示整个序列, 也可表示离散信号在nT瞬间的值。
(2)用MATLAB编程产生上述任意3种序列(长度可输入确定,对(d) (e) (f)中的参数可自行选择),并绘出其图形;
(3) 混叠现象
对连续信号其中,进行采样,分别取采样频率,观察的变化,并做记录(打印曲线),观察随着采样频率降低频谱混叠是否明显存在,说明原因。
(4)截断效应
给定,截取一定长度的信号,为窗函数,长度为N,。分别取N=6,8,12,计算的N点DFT变换,画出其幅频特性曲线;做2N点DFT变换,分析当N逐渐增大时,分析是否有频谱泄露现象、主瓣的宽度变化?如何减小泄露?
(5)栅栏效应
给定,分别计算在频率区间上的16点、32点、64点等间隔采样,绘制采样的幅频特性图,分析栅栏效应,如何减小栅栏效应?
(1)熟悉用双线性变换法和脉冲响应不变法设计IIR数字滤波器的原理与方 法;
(2)学会调用MATLAB信号处理工具中滤波器设计函数,设计各种IIR滤波 器,学会根据滤波需求确定滤波器指标参数;
(3)掌握用窗函数法设计FIR数字滤波器的原理和方法;
(一)IIR滤波器
模拟滤波器设计
巴特沃兹滤波器的振幅平方函数为
(1)
其传输函数为
(2)
(3)
首先确定技术指标:
i. 通带中允许的最大衰减和通带截止频率;
ii. 阻带允许的最小衰减和阻带起始频率。
由式(8-11)可得:
(4)
(5)
得到
(6)
(7)
再利用上面两式得到
令,
则 (8)
已知,可由式(8)求出滤波器的阶数N。求出的N可能有小数部分一般取大于等于N的最小整数。关于3dB截止频率,有时在技术指标中给出,如果没有给出可以按照式(6)或式(7)求出。
根据以上所述,巴特沃兹滤波器的设计步骤为:
i. 根据要求,由式(8)求出阶数N;
ii. 由式(6)或式(7)求出3dB截止频率;
iii. 由式(3)求出N个极点;
iv. 由式(2)写出传递函数。
双线性变换法和脉冲响应不变法实现将数字频率和模拟频率相互映射,具体见教材相关介绍。
(二)FIR滤波器
设所希望得到的滤波器的理想频率响应为。那么FIR滤波器的设计就在于寻找一个传递函数去逼近。在这种逼近中最直接的一种方法是从单位取样响应序列着手,使逼近理想的单位取样响应。我们知道可以从理想频率响应通过傅里叶反变换来得到,即:
(9)
但是一般来说,这样得到的单位取样响应往往都是无限长序列;而且是非因果的。以一个截止频率为的线性相应位理想低通为例来说明。设低通滤波器的时延为,即:
(10)
则
这是一个以为中心的偶对称的无限长非因果序列。这样一个无限长的序列怎样用一个有限长序列去近似呢?最简单的办法就是直接截取它的一段来代替它。例如把到的一段截取来作为,但是为要保证所得到的是线性相位滤波器。必须满足的对称性,所以时延应该取长度的一半,即
这种直接截取的办法可以形象地想象为,好比是通过一个“窗口”所看到的一段。中表达为和一个“窗口函数”的乘积。在这里,窗口函数就是矩形脉冲函数,即
但是一般来说,窗口函数并不一定是矩形函数,可以在矩形以内还对作一定的加权处理,因此,一般可以表示为
这里就是窗口函数。这种对理想单位取样响应加窗的处理对频率响应会产生以下三点影响:
a) 使理想特性不连续的边沿加宽,形成一过渡带,过渡带的宽度取决于窗口频谱的主瓣宽度。
b) 在过渡带两旁产生肩峰和余振,它们取决于窗口频谱的旁瓣;旁瓣越多,余振也越多;旁瓣相对值越大,肩峰则越强。
c) 增加截取长度N,只能缩小窗口频谱的主瓣宽度而不能改变旁瓣的相对值;旁瓣与主瓣的相对关系只决定于窗口函灵敏的形状。因此增加N,只能相对应减小过渡带宽。而不能改变肩峰值。肩峰值的大小直接决定通带内的平稳和阻带的衰减,对滤波器性能有很大关系。例如矩形窗的情况下,肩峰达8.95%,致使阻带最小衰减只有21分贝,这在工程上往往是不够的。怎样才能改善阻带的衰减特性呢?只能从改善窗口函数的形状上找出路,所以希望的窗口频谱中应该减少旁瓣,使能量集中在主瓣,这样可以减少肩峰和余振,提高阻带的衰减。而且要求主瓣宽度尽量窄,以获得较陡的过渡带,然而这两个要求总不能兼得,往往需要用增加主瓣宽度带换取决瓣的抑制,于是提出了海明窗、凯塞-贝塞尔窗、切比雪夫窗等窗口函数。
(1) 设计IIR数字滤波器,绘制幅频响应特性曲线。
基于脉冲响应和双线性变换设计低通数字滤波器:要求其通带截至频率100Hz,阻带截至频率200Hz,通带衰减小于2dB,阻带衰减大于15dB,采样频率Fs=500HZ。
(2)设计FIR数字滤波器,绘制幅频响应特性曲线。
(a)用矩形窗函数法设计一个线性相位FIR低通滤波器,性能指标:通带截止频率,阻带截止频率,阻带衰减为20dB,通带衰减为3dB。
(b) 当阻带衰减不小于40dB,通带衰减不大于3dB,应选择何种窗函数,画出其幅度响应的曲线。
2012上期090805-001班《数字信号处理》试题(A)
解 答
一、选择题(每小题4分,共20分)
?
1. 序列和?n?
n?sin??????4????n?2?等于
(A) 1 ; (B) ? ; (C) ??n?2? ; (D) u?n?2?
?
解:sin?n??
n??????4???n?2??sin?2??
???4???1
所以正确答案为(A)。
2. 卷积和u?n??????n?2????n?3???等于
(A) ??n?2?; (B) ??n?3?; (C) ??n?2????n?3?; (D) 解:u?n??????n?2????n?3????u?n?2??u?n?3??1
所以正确答案为(D)。
3. 已知系统 T[x(n)]?x(n)?x(?n) ,该系统是
(A) 线性时不变系统; (B) 非线性时不变系统;
(C) 线性时变系统; (D) 非线性时变系统
解:① ? T[ax1(n)?bx2(n)]?a[x1(n)?x1(?n)]?b[x2(n)?x2(?n)] ?aT[x1(n)]?bT[x2(n)]满足叠加原理 ? 是线性系统。
② ?T[x(n?m)]?x(n?m)?x(?n?m)
y(n?m)?x(n?m)?x(?(n?m))?x(n?m)?x(?n?m)T[x(n?m)]?y(n?m)
? 是时变系统。
所以正确答案为(C)。
4. 某线性时不变系统的单位脉冲响应为h?n??0.5nu?n?,该系统是
(A) 非因果稳定系统; (B) 因果稳定系统;
(C) 因果非稳定系统; (D) 非因果非稳定系统
1 1
解:① n?0时,h?n??0.5nu?n??0,故此系统是因果系统;
②
n????h?n???0.5n?0??n?1?2??,故此系统是稳定系统; 1?0.5
所以正确答案为(B)。
5. 已知IIR数字滤波器的系统函数 H(z)?1,该数字滤波器为 1?0.9z?1
(A) 低通数字滤波器; (B) 高通数字滤波器;
(C) 带通数字滤波器; (D) 带阻数字滤波器 解:H(z)?1z? 系统的零点为z?0, 极点为z??0.9?11?0.9zz?0.9
极点在实轴的-0.9处,通带中心在???处。毫无疑问,这是一个高通滤波器。 所以正确答案为(B)。
二、(12分)已知模拟信号x(t)?cos(3000?t)sin(1000?t),若采用FFT算法对其进行
频谱分析,要求频率分辨率?f?20Hz,试确定:
⑴ 信号的采集时间长度T1; ⑵ 信号的采样间隔T2;
⑶ 信号的采样点数N1; ⑷ 频率分辨率提高一倍的采样点数N2。 解:信号的最高频率分量的频率fm?1500?500?2000Hz
⑴ 由分辨率的要求确定信号的采集时间长度: T1?11??0.05s?50ms ?f20
⑵ 信号的采样间隔应满足: T2?11??0.25ms 2fm2?2000
⑶ 采样点数应满足: 2fm2?2?103
N1???200 ?f20
⑷ 频率分辨率提高一倍的采样点数应满足: 2fm2?2?103
N2???400 ?f/210
2
(15分)已知线性时不变系统的差分方程:
y(n)?3
4y(n?1)?1
8y(n?2)?x(n)?1
3x(n?1)
⑴ 求系统的系统函数H(z),并画出零极点分布图; ⑵ 若该系统是因果稳定的,写出H(z)的收敛域; ⑶ 求该因果稳定系统的单位脉冲响应。
解:⑴ 对差分方程两边取z变换:
Y(z)?3?11
4zY(z)?8z?2Y(z)?X(z)?1
3z?1X(z)
系统函数:
1?1z?11?1z?1
H(z)?Y(z)
X(z)??
1?3
4z?1?1
8z?2??1?1??1?1?
?1?2z????1?4z??
零点:z??11
3,0 极点:z?2,1
4
零极点分布图:
⑵ 由于系统是因果稳定系统,故收敛域:z?1
2
⑶ 对系统函数H(z)求z反变换即得单位抽样响应h(n)。 用部分分式法:
1?1z?1??z?1?
H?z?????3??z
?1?1??1?1??1??1?
?1?2z????1?4z????z?2????z?4??
3 三、
1??z???H?z?AA3????1?2 1??1?z?z?z?z?z?????242??4??
?H?z??1?10?A1?Res??z?? ???11z23????1?z???z?????z??122??4?z??2z?1?H?z??1?7?A2?Res??z??? ???11z43????1?z???z?????z??1424?z????4z?1107z?z
?H(z)?? z?z?24
1根据收敛域:z? 得: 2
?10?1?n7?1?n?h(n)????????u?n? ??3?2?3?4???
四、(15分) 已知序列x(n)的z变换为X(z)?
种不同ROC的z反变换。
解:设 1 ,求其可能对应的几(1?2z?1)(1?3z?1)
X(z)?
有 AB? ?1?11?2z1?3z
A?(1?2z?1)X(z)
故 X(z)?z?2z?3??2?3 B?(1?3z)X(z)?1?23? ?1?11?2z1?3z
由于X(z)有两个极点:z?2,z?3。所以X(z)的三个不同ROC分别为:
ROC1:
ROC3:
z?3z?24 ROC2:2?z?3
于是可得X(z)的三个不同的ROC对应的序列分别为:
ROC1:z?3x(n)?(?2n?1?3n?1)u(n)ROC2:2?z?3
x(n)??2n?1u(n)?3n?1u(?n?1) ROC3:
z?2
x(n)?(2n?1?3n?1)u(?n?1)
五、(10分) 已知序列x(n)?R16(n),求x(n)的32点DFT和64点DFT。
解:求x(n)的DTFT:
15
?j16?X?ej????e?j?n
?1?e?j8??ej8??e?j8???j15?sin?8??n?0
1?e?j??e
e?j?2??e2?ej?2?j?
?sin?/2??e2?
?
求x(n)的32点DFT:
X?k??X?ej??
??
2?32
k???
?e?j152??16
k
sin??8?16k??sin?
?1?2??16k??????e?j15sin?k??32
?k
?2?sin????32k???
求x(n)的64点DFT:
X?k??X?ej??
??
2?64
k??15?sin?8?k?
?e?j2?32
k
32?sin?
?1?2??32k???????e?j15sin?k?64
?k
?4?sin????64k???
5
六、(10分)设x(n)是长度为30的因果序列,h(n)是长度为10的因果序列。
X(k)?DFT[x(n)]30
Yc(k)?H(k)X(k)
H(k)?DFT[h(n)]30
yc(n)?IDFT[Yc(k)]30
y(n)?h(n)?x(n)
试确定在哪些点上yc(n)?y(n),并解释为什么? 解:yc(n)?IDFT[Yc(k)]30=x(n)?h(n)?
m???
?y(n?30m)R
?
30
(n)
只有在如上周期延拓序列中无混叠的点上,才满足yc(n)?y(n)。 而目前yc(n)长度为30,y(n)长度为30+10-1=39。 故当0?n?8时存在时域混叠。
?
yc(n)?y(n)=h(n)?x(n),
9?n?29
七、(18分)已知两个4点序列分别为x1(n)??1,2,0,1?,x2(n)??2,2,1,1?,试利用一次
4点DFT运算来同时计算DFT?x1(n)??X1(k)和DFT?x2(n)??X2(k)。 解: 利用两个4点实序列构造一个复序列
w(n)?x1(n)?jx2(n)??1?2j,2?2j,j,1?j? 则
W(k)?DFT?w(n)??
?w(n)W
n?0
3
nk
4?w(0)W40?w(1)W4k?w(2)W42k?w(3)W43k
W4?e
?j
2?4
?e
?j
?
2
??j
W(0)?w(0)W40?w(1)W40?w(2)W40?w(3)W40?4?6j
W(1)?w(0)W40?w(1)W41?w(2)W42?w(3)W43?2W(2)?w(0)W?w(1)W?w(2)W?w(3)W??2W(3)?w(0)W40?w(1)W43?w(2)W46?w(3)W49?2j
04
24
44
64
?W(k)??4?6j,2?,2?,j 2
W((k))4R4(k)??4?6j,2,?2,2j?W((N?k))4R4(k)??4?6j,?2j,?2,2?
?
6
由x1(n)?Re?w(n)?得
X1(k)?DFT?x1(n)??DFT?Re?w(n)???Wep(k) ?1
2??W((k))?
4?W((N?k))4??R4(k) ?1
2?8,2?2j,?4,2?2j?
??4,1?j,?2,1?j?
由x2(n)?Im?w(n)?得
X)?DFT?x?Im?w(n)???12(k2(n)??DFTjWop(k)
?1
2j??W((k))?
4?W((N?k))4??R4(k)
?1
2j?12j,2?2j,0,?2?2j?
??6,1?j,0,1?j?
7
南京信息工程大学数字信号处理实验报告学院:电子与信息工程学院班级:11通信1班学号:XXX姓名:XX指导教师:XX20XX/12/…
数字信号处理第一次实验报告学院:信息工程学院专业:电子信息工程二班学号:***姓名:实验一:系统响应及系统稳定性1.实验目的(1)…
数字信号处理实验报告专业电子信息工程学号111308336姓名张强伟实验一数字滤波器的结构一实验目的1加深对数字滤波器分类与结构的…
北京信息科技大学实验报告课程名称数字信号处理实验项目IIR数字滤波器设计实验仪器计算机MATLAB软件学院信息与通信工程学院班级姓…
西南大学学生实验报告姓名杨剑学号2220xx3220xx058班级1班专业电科实验日期20xx年9月29日实验学时2学时实验名称基…
数字信号处理认识报告20xx90512101文通084-1安森令《数字信号处理》是我们通信工程和电子类专业的一门重要的专业基础课程…
南京信息工程大学数字信号处理实验报告学院:电子与信息工程学院班级:11通信1班学号:XXX姓名:XX指导教师:XX20XX/12/…
北京信息科技大学实验报告课程名称数字信号处理实验项目IIR数字滤波器设计实验仪器计算机MATLAB软件学院信息与通信工程学院班级姓…
目录仿真一LMS算法和RLS算法11自适应滤波的基本原理111自适应最小均方LMS算法112递归最小二乘方RLS算法22仿真实验4…
数字信号处理B课程项目题目组号任课教师实施报告数字音效处理器4020xx年10月18日1报告提纲一课程项目要求设计要求设计一个数字…
机械振动基实验报告机械0814中南大学1机械振动实验报告姓名学号0806081529成绩指导教师隆志力23机械振动实验报告姓名周菁…