赵雪英10版《数字信号处理》实验讲义要点

《数字信号处理》

实验讲义

信息学院

赵雪英

2013.1

1

前 言

数字信号处理是利用计算机或专用数字处理设备，采用数值计算的方法对信号进行处理的一门学科，它包括数据采集，变换、分析、综合、滤波、估值与识别等加工处理，以便于提取信息和应用。

数字信号处理的主要优点有：

（1） 灵活性好。适合用计算机、可编程器件（如通用单片机、DSP、可编程逻辑器件等）

实现，通过编程很容易改变数字信号处理系统得参数，从而使系统实现各种不同的处理功能。如数字电话系统中采用的时分复用技术。

（2）稳定可靠。 （3）处理精度高。 （4）便于加解密。

（5）便于大规模集成化、小型化。 （6）便于自动化、多功能化。

（7）可实现模拟系统无法实现的复杂处理功能。

数字信号处理原理、实现和应用是本学科研究和发展的三个主要方面。数字信号处理应用非常广泛，涉及语音、雷达、声呐、地震、图像处理、通信系统、系统控制、生物医学工程、机械振动、遥感遥测、航空航天、电力系统、故障检测和自动化仪表等领域。

MATLAB是美国MathWorks公司开发的一种功能极其强大的高技术计算语言和内容极其丰富的软件库。其中的信号处理工具箱是一个内容丰富的信号处理软件库，是学习、应用数字信号处理的一个极好工具。在学习数字信号处理理论的同时，熟练掌握MATLAB的使用，对理工科的学生是非常必要的。

2

目 录

实验一 时域离散信号和系统 ................................................................................................................................. 4 实验二 时域离散信号和系统的频域分析 ............................................................................................................. 6 实验三 离散傅里叶变换及其快速算法 ................................................................................................................. 8 实验四 特殊滤波器 ................................................................................................................................................. 9 实验五 IIR数字滤波器设计 ................................................................................................................................. 10 实验六 FIR数字滤波器设计 ................................................................................................................................ 12 实验七 综合实验-数字滤波器设计 ...................................................................................................................... 14 实验八 时域离散系统的实现 ............................................................................................................................... 15

3

实验一 时域离散信号和系统

一、实验目的

（一）常用时域离散信号的MATLAB表示 （二）应用MATLAB求线性卷积 （三）应用MATLAB求解差分方程 二、实验内容

（一）常用时域离散信号的MATLAB表示

用两个参数向量x和n表示有限长序列x(n),x是x(n)的样值向量，n是位置向量； n与x长度相等。

序列 单位采样序列?(n) 单位阶跃序列u(n) 实指数序列x(n)?an,?n;a?R 复指数序列x(n)?e(??j?0)nMATLAB函数 zeros(m,n) ones(m,n) .^ exp ,?n 正余弦序列x(n)?cos(?0n??),?n cos(或sin)函数 例：应用MATLAB表示正弦序列：x(n)?sin(?n/5),?10?n?10 （二）应用MATLAB计算有限长序列的线性卷积 （1）y?conv(x,h)

例1.3.2 设x(n)?2?(n)??(n?1)?2?(n?2),h(n)??(n)?2?(n?1)??(n?2),求

y(n)?x(n)*h(n)

（2）[y,ny]?convu(h,nh,x,nx)

4

例：h(n)?x(n)?R5(n?2),求:y(n)?x(n)*h(n) （三）应用MATLAB求解差分方程

已知差分方程：y(n)??bix(n?i)??aiy(n?i)，求：

i?0i?1MN（1）零状态解：yn=filter(b,a,x)，其中，b?[b0,b1,...,bM],a?[a0,a1,...,aN] （2）全解：xi=filtic(b,a,ys,xs) 其中，b?[b0,b1,...,bM],a?[a0,a1,...,aN],

ys?[y(?1),y(?2),y(?3),...,y(?N)],xs?[x(?1),x(?2),x(?3),...x(?M)]

yn=filter(b,a,x,xi)

（3）单位脉冲响应h(n)

[h,t]=impz(b,a) [h,t]=impz(b,a,N) impz(b,a)

（4）单位阶跃响应s(n)

[s,t]=stepz(b,a) [s,t]=stepz(b,a,N) stepz(b,a)

（5）判断系统的稳定性。

n????|h(n)|?M???sum(abs(h))

?|zk|?1,k?1,2,...N?z?roots(a)例1.4.1：已知差分方程y(n)?0.8y(n?1)?x(n)，x(n)??(n)，初始条件y(?1)?1，求： （1）差分方程的全解；（2）单位脉冲响应h(n)；（3）单位阶跃响应s(n)； （4）判断系统的稳定性

5

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说教育文库赵雪英10版《数字信号处理》实验讲义要点在线全文阅读。