数字信号处理实验指导书(.docVIP

实验一 常见离散信号的MATLAB产生和图形显示 1、实验目的： （１）熟悉MATLAB应用环境，常用窗口的功能和使用方法。 （２）加深对常用离散时间信号的理解。 （３）掌握简单的绘图命令。 （４）掌握线性卷积的计算机编程方法。 2、实验仪器： PC机一台 MATLAB软件 3、实验原理： （１）单位抽样序列 如果在时间轴上延迟了k个单位，得到即： （２）单位阶跃序列 （３）矩形序列 （４）正弦序列 （５）复正弦序列 X=exp();x1=abs();x2=angle();x3=real();x4=imag();subplot(221);stem(n,x3);subplot(n,x4) （６）指数序列 （７）线性时不变系统的响应为如下的卷积计算式： 4、实验内容及步骤： （１）复习常用离散时间信号的有关内容。 （２）编制程序产生上述6种序列（长度可输入确定，对(4) (5) (6)中的参数可自行选择），并绘出其图形。 （３）已知系统的单位脉冲响应输入信号，试用卷积法求解系统的输出，并绘出、及图形。 5、实验用MATLAB函数介绍 （１）数字信号处理中常用到的绘图指令(只给出函数名，具体调用格式参看help) figure(); plot(); stem(); axis(); grid on; title(); xlabel(); ylabel(); text(); hold on; subplot() （２）离散时间信号产生可能涉及的函数 zeros(); ones(); exp(); sin(); cos(); abs(); angle(); real(); imag(); （３）卷积计算可能涉及的函数 conv(); length() 注：实验过程中也可以使用自己编制的自定义函数，如impseq()、stepseq()等。 6、思考题： MATLAB在处理数据运算时与其他高级语言相比有何不同？ MATLAB是以复数矩阵为基本编程单元的一种程序设计语言，我们可以直接以矩阵和数组为单位，进行数组和矩阵的运算。而其他高级语言没有这个功能。 从第一题的计算中就可以看出，使用矩阵运算比逐点计算的速度要快得多。 另外，MATLAB中有多种数学、工程类的函数可以调用，如这次使用的fft函数，极大的方便了我们进行工程计算。 7、实验报告要求 （１）简述实验目的及原理。 （２）总结实验所得主要结论。 （３）简要回答思考题。 实验二 信号、系统及系统响应 1、实验目的： （１）熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。 （２）熟悉时域离散系统的时域特性。 （３）利用卷积方法观察分析系统的时域特性。 （４）掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。 2、实验仪器： PC机一台 MATLAB软件 3、实验原理： 采样是连续信号数字处理的第一个关键环节。对一个连续信号进行理想采样的过程可用下式表示。 其中为的理想采样，为周期冲激脉冲， 即 ； 由频域卷积定理，得　　※ 上式表明，为的周期延拓，其延拓周期为采样角频率()。采样前后的频谱示意图见课本。只有满足采样定理时，才不会发生频率混叠失真。 在计算机上用高级语言计算很不方便，下面给出用序列的傅里叶变换来计算的方法。 课本中(2.4.7)式，表示序列的傅里叶变换和模拟信号的傅里叶变换之间的关系式。与※式比较，可得，这说明两者之间只在频率度量上差一个常数因子。实验过程中应注意这一差别。 　　为了在数字计算机上观察分析各种序列的频域特性，通常对在上进行M点采样来观察分析。 对长度为N的有限长序列x(n), 有 其中 通常应取得大一些，以便观察谱的细节变化。取模可绘出幅频特性曲线。 一个时域离散线性时不变系统的输入输出关系为 上述卷积运算也可以在频域实现，但要注意右边的相乘是在个频点上的频谱值相乘。 4. 实验内容及步骤 （１）认真复习采样理论、离散信号与系统、线性卷积、序列的傅里叶变换及性质等有关内容，阅读本实验原理与方法。 （２）编制并运行各实验程序，完成下述实验内容： ①分析采样序列的特性。 首先，对连续信号 进行采样，得到采样序列 其中为幅度因子， 为衰减因子，是模拟角频率，为采样间隔。 这些参数可以在实验过程中由键盘输入，也可以在程序开始设定。如可使，，。 分析不同采样频率时所对应的频谱混叠情况。分别取采样频率， 观察的变化， 并做记录(打印曲线)；观察随着采样频率降低频谱混叠是否明显存在，说明原因。 ②时域离散信号、系统和系统响应分析。 先编写单位脉冲序列： 再编写系统的单位