中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 实验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 2k / 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
中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 产生如下系统单位脉冲响应序列,本实验中用到两种FR系统: (1)h(n)=Rio(n) (2)h,(n)=6nm)+2.5δ(n-1)+2.5δ(n-2)+6(n-3) 1.3.4.2离散信号、系统和系统响应的分析 (1)观察信号x(n)和系统h,()的时域和幅频特性:利用线性卷积求信号通过系统以后 的响应。比较系统响应和系统h,()的时域及幅频特性。注意它们之间有无差异,绘出图形。 (2)观察信号x.()和系统h,(n)的时域和幅频特性,利用线性卷积求系统响应。判断响 应序列图形及序列非零值长度是否与理论结果一致,说出一种定性判断响应序列图形正确与 否的方法(提示:x(nm)=h(n)=Ro(n))。利用序列的傅立叶变换数值计算子程序求出 Y(e),观察响应序列的幅频特性。定性判断结果正确与否。改变信号x.(n)的矩形宽度, 使N=5,重复以上动作,观察变化,记录改变参数前后的差异。 (3)将实验步骤2-(2)中的信号换为x(n),其中A=1,a=0.4,2o=2.0734,T=1。 重复实验2-(2)各步骤,改变x()的参数=0.1再重复实验2-(2)各步骤:改变参数 2。=1.2516,重复实验2-(2)各步骤。在实验中观察改变和2。对信号及系统响应的时 域和幅频特性的影响,绘制相应的图形。 1.3.4.3卷积定律的验证 利用式(1-l4)将x(n)和系统hn(nm)的傅氏变换相乘,直接求得Y(e),将得到的 幅频特性曲线和实验2-(3)中得到的曲线进行比较,观察二者有无差异。验证卷积定律。 1.3.5*选作内容 改变信号x(t)中的衰减因子α,先定性估计频谱可能产生的变化,然后观察其频谱的 变化,记录结果,变化是否你所想的一致,这说明了什么? 一个LTI系统的冲激响应为h(n)=(0.9)”u(n),输入序列为x.(n),求系统响应H(eo) 和输出信号y(n)及其频谱Y(e):如果h(n)=x.(n),其结果又如何? 6
中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 6 产生如下系统单位脉冲响应序列,本实验中用到两种 FIR 系统: (1) ( ) ( ) ha n R10 n (2) h (n) (n) 2.5 (n 1) 2.5 (n 2) (n 3) b 1.3.4.2 离散信号、系统和系统响应的分析 (1)观察信号 x (n) b 和系统 h (n) b 的时域和幅频特性;利用线性卷积求信号通过系统以后 的响应。比较系统响应和系统 h (n) b 的时域及幅频特性。注意它们之间有无差异,绘出图形。 (2)观察信号 x (n) c 和系统 h (n) a 的时域和幅频特性,利用线性卷积求系统响应。判断响 应序列图形及序列非零值长度是否与理论结果一致,说出一种定性判断响应序列图形正确与 否的方法(提示: ( ) ( ) ( ) xc n ha n R10 n )。利用序列的傅立叶变换数值计算子程序求出 ( ) k j Y e ,观察响应序列的幅频特性。定性判断结果正确与否。改变信号 x (n) c 的矩形宽度, 使 N=5,重复以上动作,观察变化,记录改变参数前后的差异。 (3)将实验步骤 2-(2)中的信号换为 x (n) a ,其中 A 1, 0.4,0 2.0734,T 1。 重复实验 2-(2)各步骤,改变 x (n) a 的参数 0.1 再重复实验 2-(2)各步骤;改变参数 0 1.2516 ,重复实验 2-(2)各步骤。在实验中观察改变 和0 对信号及系统响应的时 域和幅频特性的影响,绘制相应的图形。 1.3.4.3 卷积定律的验证 利用式(1-14)将 x (n) a 和系统 h (n) a 的傅氏变换相乘,直接求得 ( ) k j Y e ,将得到的 幅频特性曲线和实验 2-(3)中得到的曲线进行比较,观察二者有无差异。验证卷积定律。 1.3.5 *选作内容 改变信号 x (t) a 中的衰减因子α ,先定性估计频谱可能产生的变化,然后观察其频谱的 变化,记录结果,变化是否你所想的一致,这说明了什么? 一个 LTI 系统的冲激响应为 h(n) (0.9) u(n) n ,输入序列为 x (n) c ,求系统响应 ( ) j H e 和输出信号 y(n)及其频谱 ( ) j Y e ;如果 h(n) x (n) c ,其结果又如何?
中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 编写一个程序,将x.()分解为奇偶序列,绘制奇偶序列时域图形并求出它们频谱 Xe(eo)和Xo(eo),同x.(nm)的频谱Xc(eo)进行比较,可以得出什么结论? 针对信号x(t)=c0s6000m+cos1000π分析采样率、信号功率。 1.4 实验报告要求 1、记录实验内容中要求观察、分析、比较的内容及结果并进行分析。 2、总结在上机实验内容中要求比较时域、幅频曲线差异差异部分内容的结果,定性分析它 们正确与否,并简要说明这些结果的含义。 3、在实验报告中附上在实验过程中记录的各个典型信号序列的时域和幅频特性曲线,分析 所得到的结果图形,说明各个信号的参数变化对其时域和幅频特性的影响。 4、总结一下你在用MatLab进行数字信号处理实验项目的时候常用的函数及其功能。 5、总结实验中根据实验现象得到的其他个人结论
中国科学技术大学电子工程与信息科学系 《数字信号处理》课程基本实验 7 编写一个程序,将 x (n) c 分解为奇偶序列,绘制奇偶序列时域图形并求出它们频谱 ( ) j Xe e 和 ( ) j Xo e ,同 x (n) c 的频谱 ( ) j Xc e 进行比较,可以得出什么结论? 针对信号 x(t) cos6000t cos1000t 分析采样率、信号功率。 1.4 实验报告要求 1、 记录实验内容中要求观察、分析、比较的内容及结果并进行分析。 2、 总结在上机实验内容中要求比较时域、幅频曲线差异差异部分内容的结果,定性分析它 们正确与否,并简要说明这些结果的含义。 3、 在实验报告中附上在实验过程中记录的各个典型信号序列的时域和幅频特性曲线,分析 所得到的结果图形,说明各个信号的参数变化对其时域和幅频特性的影响。 4、 总结一下你在用 MatLab 进行数字信号处理实验项目的时候常用的函数及其功能。 5、 总结实验中根据实验现象得到的其他个人结论