中国科学技术大学：《数字信号处理》课程教学资源（实验要求）实验1 信号及系统基本特性分析

中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 实验1 信号及系统基本特性分析 1.1 实验目的 l、学习Matlab编程的基本方法：掌握常用函数用法。 2、了解不同信号的频域特性，理解时域特性与频域特性之间的关联性。 3、掌握典型信号序列的时域和频域基本特性。 4、熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对采样定理的理解。 5、了解离散系统的时域/频域特性及其对输出信号的影响，掌握系统分析方法。 1.2 实验原理 1.2.1 连续时间信号的采样 采样是从连续时间信号到离散时间信号的过渡桥梁，对采样过程的研究不仅可以了采 样前后信号时域和频域特性发生的变化以及信号内容不丢失的条件，而且有助于加深对拉氏 变化、傅氏变换、z变换和序列傅氏变换之间关系的理解。 对一个连续时间信号进行理想采样的过程可以表示为该信号的一个周期冲激脉冲的乘 积，即 元()=x()M) (1-1) 其中x(t)是连续信号x。(t)的理想采样，M(t)是周期冲激脉冲 M()=8u-nT) (1-2) 它也可以用傅立叶级数表示为： (1-3) 其中T为采样周期，2，=2π/T是采样角频率。设X,(S)是连续时间信号x(t)的双边拉 氏变换，即有： X(s)=∫x,()e"dh (1-4) 此时理想采样信号文()的拉氏变换为

中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 1 实验1 信号及系统基本特性分析 1.1 实验目的 1、 学习 Matlab 编程的基本方法；掌握常用函数用法。 2、 了解不同信号的频域特性，理解时域特性与频域特性之间的关联性。 3、 掌握典型信号序列的时域和频域基本特性。 4、 熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对采样定理的理解。 5、 了解离散系统的时域/频域特性及其对输出信号的影响，掌握系统分析方法。 1.2 实验原理 1.2.1 连续时间信号的采样 采样是从连续时间信号到离散时间信号的过渡桥梁，对采样过程的研究不仅可以了采 样前后信号时域和频域特性发生的变化以及信号内容不丢失的条件，而且有助于加深对拉氏 变化、傅氏变换、z 变换和序列傅氏变换之间关系的理解。 对一个连续时间信号进行理想采样的过程可以表示为该信号的一个周期冲激脉冲的乘 积，即 x ˆ (t) x (t)M(t) a  a （1-1） 其中 x ˆ (t) a 是连续信号 x (t) a 的理想采样， M (t) 是周期冲激脉冲 M (t) (t nT) n       （1-2） 它也可以用傅立叶级数表示为：      n jm t s e T M t 1 ( ) （1-3) 其中 T 为采样周期， s  2 /T 是采样角频率。设 X (s) a 是连续时间信号 x (t) a 的双边拉 氏变换，即有：     X s  x t e dt st a a ( ) ( ) （1-4） 此时理想采样信号 x ˆ (t) a 的拉氏变换为

中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 (s)=∫x(0ed (1-5) 号芝-n) 作为拉氏变换的一种特例，信号理想采样的傅立叶变换 .U0=2xU0-a,】 (1-6) 由式(1-5)和式(1-6)可知，信号理想采样后的频谱是原信号频谱的周期延拓，其延拓周期 等于采样频率。根据Shannon取样定理，如果原信号是带限信号，且采样频率高于原信号最 高频率分量的2倍，则采样以后不会发生频谱混淆现象。 在计算机处理时，不采用式(1-6)计算信号的频谱，而是利用序列的傅立叶变换计算信号 的频谱，定义序列n)=x(nT)=x(t=x(t)M(),根据Z变换的定义，可以得到序 列x(n)的Z变换为： K()-m (1-7) 以e°代替上式中的z,就可以得到序列xn)的傅立叶变换 Xe)=艺ne (1-8) 式(1-6)和式(1-8)具有如下关系： X(j)=X(e)lo-or (1-9) 由式(1-9)可知，在分析一个连续时间信号的频谱时，可以通过取样将有关的计算转 化为序列傅立叶变换的计算。 1.2.2有限长序列分析 一般来说，在计算机上不可能，也不必要处理连续的曲线X(),通常，我们只要观 察、分析X(o)在某些频率点上的值。对于长度为N的有限长序列 2

中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 2     X s  x t e dt st a a ( ) ˆ ( ) ˆ （1-5) 作为拉氏变换的一种特例，信号理想采样的傅立叶变换 [ ( )] 1 ( ) ˆ s m a Xa j m T X j        （1-6） 由式(1-5)和式(1-6)可知，信号理想采样后的频谱是原信号频谱的周期延拓，其延拓周期 等于采样频率。根据 Shannon 取样定理，如果原信号是带限信号，且采样频率高于原信号最 高频率分量的 2 倍，则采样以后不会发生频谱混淆现象。 在计算机处理时，不采用式(1-6)计算信号的频谱，而是利用序列的傅立叶变换计算信号 的频谱，定义序列 x(n) x (nT) x ˆ (t) x (t)M(t)  a  a  a ，根据 Z 变换的定义，可以得到序 列 x(n)的 Z 变换为：      n n X (z) x(n)z （1-7） 以 j e 代替上式中的 z，就可以得到序列 x(n)的傅立叶变换      n j n j X e x n e   ( ) ( ) （1-8) 式（1-6）和式（1-8）具有如下关系： T j a X j X e     ( ) ( ) | ˆ （1-9） 由式（1-9）可知，在分析一个连续时间信号的频谱时，可以通过取样将有关的计算转 化为序列傅立叶变换的计算。 1.2.2 有限长序列分析 一般来说，在计算机上不可能，也不必要处理连续的曲线 ( ) j X e ，通常，我们只要观 察、分析 ( ) j X e 在某些频率点上的值。对于长度为 N 的有限长序列                             ( ) 1 ( ) 1 1 ( ) ( ) a s m s jm t a m st m jm t a X s jm T x t e dt T e e dt T x t s s

中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 f(m),0≤n≤N-1 x(n)= (1-10) 0,其他n 般只需要在0-2π之间均匀地取M个频率点，计算这些点上的序列傅立叶变换 x(e)( (1-11) 其中ok=2水/M,k=O,l,M-l。X(e)是一个复函数，它的模就是幅频特性曲线。 1.2.3信号卷积 一个线性时不变离散系统的响应y(n)可以用它的单位冲激响应h(n)和输入信号x(n)的卷 积来表示： y(n)=x(n)*h(n)= 艺mam-m (1-12) =- 根据傅立叶变换和Z变换的性质，与式(1-12)对应应该有 Y(2)=X(H() (1-13) Y(e)=X(e)H(e) (1-14) 式(1-12)告诉我们可以通过对两个序列的移位、相乘、累加计算信号响应：而式(1-14) 告诉我们卷积运算也可以在频域上用乘积实现。 1.3 实验内容 1.3.1 Matlab操作与使用 根据所提供的Matlab操作指南学习Matlab的使用。完成文件操作：矩阵运算：绘图： 图形界面的实现等功能，学会使用Matlab联机帮助查找信息。 1.3.2理想采样信号序列的特性分析 对信号xa(t)=Ae“sin(2o1)(t)进行理想采样，可以得到一个理想的采样信号序列： x()=Ae-“sin(2onT),0≤n<50,其中A为幅度因子，a是衰减因子，2o是频率。T 为采样周期。 3

中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 3        n f n n N x n 0,其他 ( ),0 1 ( ) （1-10） 一般只需要在 0 2 之间均匀地取 M 个频率点，计算这些点上的序列傅立叶变换      1 0 jωk n ( ) ( ) N n j k X e x n e  （1-11） 其中 k  2k / M ，k=0,1,..,M-1。 ( ) j k X e  是一个复函数，它的模就是幅频特性曲线。 1.2.3 信号卷积 一个线性时不变离散系统的响应 y(n)可以用它的单位冲激响应 h(n)和输入信号 x(n)的卷 积来表示：        m y(n) x(n) h(n) x(m)h(n m) (1-12) 根据傅立叶变换和 Z 变换的性质，与式（1-12）对应应该有 Y(z)  X(z)H(z) （1-13） ( ) ( ) ( ) j j j Y e  X e H e （1-14） 式（1-12）告诉我们可以通过对两个序列的移位、相乘、累加计算信号响应；而式（1-14） 告诉我们卷积运算也可以在频域上用乘积实现。 1.3 实验内容 1.3.1 Matlab 操作与使用 根据所提供的 Matlab 操作指南学习 Matlab 的使用。完成文件操作；矩阵运算；绘图； 图形界面的实现等功能，学会使用 Matlab 联机帮助查找信息。 1.3.2 理想采样信号序列的特性分析 对信号 ( ) sin( ) ( ) 0 x t Ae t u t t a    进行理想采样，可以得到一个理想的采样信号序列： ( )  sin(0 ),0   50  x t Ae nT n t a  ，其中 A 为幅度因子，  是衰减因子， 0 是频率。T 为采样周期

中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 产生理想采样信号序列，(m,使A=44.128,a=50√2x,2,=50W2r.q)首 先选用采样频率为1000Hz,T=11000,观察所得理想采样信号的幅频特性，并做记录：(2) 改变采样频率为300Hz,T=1/300,观察所得到的频谱特性曲线的变化，并做记录：（3)进 一步减小采样频率为200Hz,T=1/200,观察频谱“混淆”现象是否明显存在，说明原因， 并记录这时候的幅频特性曲线。 1.3.3典型信号序列的特性分析 1.3.3.1信号序列产生 产生如下基本信号： _(m-p2 (1)高斯序列： xaa(n)= 9,0≤n≤15 0,else e-msin2wfim,0≤n≤15 (2)衰减正弦序列：x66(m)= 0.else n+1,0≤n≤3 (3)三角波序列： xcc(n)=8-n,4≤n≤7 0,else 4-n,0≤n≤3 (4)反三角序列： xa(nm)={n-3,4≤n≤7 0,else 1.3.3.2观赛高斯序列的时域和频域特性 ①固定信号x()中的参数p=8,改变q的值，使q分别等于2,4,8。观察它们的时 域和幅频特性，了解q取不同值的时候，对信号时域特性和幅频特性的影响。②固定q=8, 改变p,使p分别等于8,13,14，观察参数p变化对信号序列时域及幅频特性的影响。注 意p等于多少时，会发生明显的泄漏现象，混淆现象是否也随之出现？记录实验中观察到的 现象，绘制相应的时域序列和幅频特性曲线。 1.3.3.3观赛衰减正弦序列的时域和幅频特性 针对信号x(n):①令a=0.1并且f=0.0625,检查谱峰出现的位置是否正确，注意频

中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 4 产生理想采样信号序列 x (n) a ，使 A  444.128，   50 2,0  50 2 。（1）首 先选用采样频率为 1000Hz，T=1/1000，观察所得理想采样信号的幅频特性，并做记录；（2） 改变采样频率为 300Hz，T=1/300，观察所得到的频谱特性曲线的变化，并做记录；（3）进 一步减小采样频率为 200Hz，T=1/200，观察频谱“混淆”现象是否明显存在，说明原因， 并记录这时候的幅频特性曲线。 1.3.3 典型信号序列的特性分析 1.3.3.1 信号序列产生 产生如下基本信号： （1）高斯序列：           else x n e n q n p a a 0, ( ) ,0 15 2 ( ) （2）衰减正弦序列：        else e fn n x n n b b 0, sin 2 ,0 15 ( )   （3）三角波序列：             else n n n n xcc n 0, 8 ,4 7 1,0 3 ( ) （4）反三角序列：             else n n n n xd d n 0, 3,4 7 4 ,0 3 ( ) 1.3.3.2 观察高斯序列的时域和频域特性 ①固定信号 x (n) aa 中的参数 p=8，改变 q 的值，使 q 分别等于 2，4，8。观察它们的时 域和幅频特性，了解 q 取不同值的时候，对信号时域特性和幅频特性的影响。②固定 q=8， 改变 p，使 p 分别等于 8，13，14，观察参数 p 变化对信号序列时域及幅频特性的影响。注 意 p 等于多少时，会发生明显的泄漏现象，混淆现象是否也随之出现？记录实验中观察到的 现象，绘制相应的时域序列和幅频特性曲线。 1.3.3.3 观察衰减正弦序列的时域和幅频特性 针对信号 x (n) bb ：①令α =0.1 并且 f=0.0625，检查谱峰出现的位置是否正确，注意频

中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 谱的形状，绘制幅频特性曲线。②改变=0.4375，再变化=0.5625，观察这两种情况下，频 谱的形状和谱峰出现的位置，有无混淆和泄漏现象发生？说明产生现象的原因。 1.3.3.4观襄三角波序列和反三角波序列的时域和幅频特性 针对信号xec(n)和xa(n):①用8点FFT分析信号xc(n)和xa(n)的幅频特性，观察 两者的序列形状和频谱曲线有什么异同？（注意：这时候的xa()可以看作是x()经过 圆周移位以后得到的)绘制两者的序列和幅频特性曲线。②在的x(n)和xa(m)末尾补零， 用16点FFT分析这两个信号的幅频特性，观察幅频特性发生了什么变化？两个信号之间的 F℉T频谱还有没有相同之处？这些变化说明了什么？ 1.3.3.5*选作内容 将x,(n)信号的长度N设为63，用MatLab中randn(1,N)函数产生一个噪声信号w(n), 计算将这个噪声信号叠加到x(n)上以后新信号y(n)=x(n)+w(n)的频谱，观察发生的 变化并记录。 在步骤2的基础上，改变参数α和，观察在出现混淆现象和泄漏现象的时候有噪声的 y()信号的频谱有什么变化，是否明显？ 1.3.4 离散信号、系统和系统响应的分析 1.3.4.1信号序列产生 产生如下信号序列： (1)理想采样信号序列：对信号x(t)=Ae“sin(2ot)u(t)进行理想采样，可以得到一个 理想的采样信号序列：x.(t)=Ae“sin(2onT),0≤n<50,其中A为幅度因子，au是衰 减因子，2。是频率。T为采样周期。 1,n=0 (2)单位脉冲序列x(n)=δ(n)= 0,n≠0 1,0≤n<N-1 (3)矩形序列x.(n)=Rw(n)= ,其中N=10 0,其他

中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 5 谱的形状，绘制幅频特性曲线。②改变 f=0.4375，再变化 f=0.5625，观察这两种情况下，频 谱的形状和谱峰出现的位置，有无混淆和泄漏现象发生？说明产生现象的原因。 1.3.3.4 观察三角波序列和反三角波序列的时域和幅频特性 针对信号 x (n) cc 和 x (n) dd ：①用 8 点 FFT 分析信号 x (n) cc 和 x (n) dd 的幅频特性，观察 两者的序列形状和频谱曲线有什么异同？（注意：这时候的 x (n) dd 可以看作是 x (n) cc 经过 圆周移位以后得到的）绘制两者的序列和幅频特性曲线。②在的 x (n) cc 和 x (n) dd 末尾补零， 用 16 点 FFT 分析这两个信号的幅频特性，观察幅频特性发生了什么变化？两个信号之间的 FFT 频谱还有没有相同之处？这些变化说明了什么？ 1.3.3.5 * 选作内容 将 x (n) b 信号的长度 N 设为 63，用 MatLab 中 randn(1,N)函数产生一个噪声信号 w(n)， 计算将这个噪声信号叠加到 x (n) b 上以后新信号 y(n) x (n) w(n)  b  的频谱，观察发生的 变化并记录。 在步骤 2 的基础上，改变参数α 和 f，观察在出现混淆现象和泄漏现象的时候有噪声的 y(n)信号的频谱有什么变化，是否明显？ 1.3.4 离散信号、系统和系统响应的分析 1.3.4.1 信号序列产生 产生如下信号序列： （1）理想采样信号序列：对信号 ( ) sin( ) ( ) 0 x t Ae t u t t a    进行理想采样，可以得到一个 理想的采样信号序列： ( )  sin(0 ),0   50  x t Ae nT n t a  ，其中 A 为幅度因子，  是衰 减因子， 0 是频率。T 为采样周期。 （2）单位脉冲序列        0, 0 1, 0 ( ) ( ) n n xb n  n （3）矩形序列         0,其他 1,0 1 ( ) ( ) n N xc n RN n ，其中 N=10