《数字信号处理》课程基本实验
实验1  信号及系统基本特性分析
1.1  实验目的
1、 学习Matlab编程的基本方法；掌握常用函数用法。
2、 了解不同信号的频域特性，理解时域特性与频域特性之间的关联性。
3、 掌握典型信号序列的时域和频域基本特性。
4、 熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对采样定理的理解。
5、 了解离散系统的时域/频域特性及其对输出信号的影响，掌握系统分析方法。
1.2  实验原理
1.2.1  连续时间信号的采样
采样是从连续时间信号到离散时间信号的过渡桥梁，对采样过程的研究不仅可以了采样前后信号时域和频域特性发生的变化以及信号内容不丢失的条件，而且有助于加深对拉氏变化、傅氏变换、z变换和序列傅氏变换之间关系的理解。
对一个连续时间信号进行理想采样的过程可以表示为该信号的一个周期冲激脉冲的乘积，即
            （1-1）
其中是连续信号的理想采样，是周期冲激脉冲
            （1-2）
它也可以用傅立叶级数表示为：
                （1-3)
其中T为采样周期，是采样角频率。设是连续时间信号的双边拉氏变换，即有：
                （1-4）
此时理想采样信号的拉氏变换为
                （1-5)
作为拉氏变换的一种特例，信号理想采样的傅立叶变换
            （1-6）
由式(1-5)和式(1-6)可知，信号理想采样后的频谱是原信号频谱的周期延拓，其延拓周期等于采样频率。根据Shannon取样定理，如果原信号是带限信号，且采样频率高于原信号最高频率分量的2倍，则采样以后不会发生频谱混淆现象。
在计算机处理时，不采用式(1-6)计算信号的频谱，而是利用序列的傅立叶变换计算信号的频谱，定义序列，根据Z变换的定义，可以得到序列x(n)的Z变换为：
                    （1-7）
以代替上式中的z，就可以得到序列x(n)的傅立叶变换
                （1-8)
式（1-6）和式（1-8）具有如下关系：
                （1-9）
短时傅里叶变换matlab程序
由式（1-9）可知，在分析一个连续时间信号的频谱时，可以通过取样将有关的计算转化为序列傅立叶变换的计算。
1.2.2  有限长序列分析
一般来说，在计算机上不可能，也不必要处理连续的曲线，通常，我们只要观察、分析在某些频率点上的值。对于长度为N的有限长序列
            （1-10）
一般只需要在之间均匀地取M个频率点，计算这些点上的序列傅立叶变换
            （1-11）
其中，k=0,1,..,M-1。是一个复函数，它的模就是幅频特性曲线。
1.2.3  信号卷积
一个线性时不变离散系统的响应y(n)可以用它的单位冲激响应h(n)和输入信号x(n)的卷积来表示：
        (1-12)
根据傅立叶变换和Z变换的性质，与式（1-12）对应应该有
                （1-13）
            （1-14）
式（1-12）告诉我们可以通过对两个序列的移位、相乘、累加计算信号响应；而式（1-14）告诉我们卷积运算也可以在频域上用乘积实现。
1.3  实验内容
1.3.1  Matlab操作与使用
根据所提供的Matlab操作指南学习Matlab的使用。完成文件操作；矩阵运算；绘图；图形界面的实现等功能，学会使用Matlab联机帮助查信息。
1.3.2  理想采样信号序列的特性分析
对信号进行理想采样，可以得到一个理想的采样信号序列：，其中A为幅度因子，是衰减因子，是频率。T为采样周期。
产生理想采样信号序列，使，。（1）首先选用采样频率为1000Hz，T=1/1000，观察所得理想采样信号的幅频特性，并做记录；（2）改变采样频率为300Hz，T=1/300，观察所得到的频谱特性曲线的变化，并做记录；（3）进一步减小采样频率为200Hz，T=1/200，观察频谱“混淆”现象是否明显存在，说明原因，并记录这时候的幅频特性曲线。
1.3.3  典型信号序列的特性分析
1.3.3.1 信号序列产生
产生如下基本信号：
（1）高斯序列：   
（2）衰减正弦序列：
（3）三角波序列： 
（4）反三角序列： 
1.3.3.2 观察高斯序列的时域和频域特性
固定信号中的参数p=8，改变q的值，使q分别等于2，4，8。观察它们的时域和幅频特性，了解q取不同值的时候，对信号时域特性和幅频特性的影响。固定q=8，改变p，使p分别等于8，13，14，观察参数p变化对信号序列时域及幅频特性的影响。注意p等于多少时，会发生明显的泄漏现象，混淆现象是否也随之出现？记录实验中观察到的现象，绘制相应的时域序列和幅频特性曲线。
1.3.3.3 观察衰减正弦序列的时域和幅频特性
针对信号：令α=0.1并且f=0.0625，检查谱峰出现的位置是否正确，注意频谱的形状，绘制幅频特性曲线。改变f=0.4375，再变化f=0.5625，观察这两种情况下，频谱的形状和谱峰出现的位置，有无混淆和泄漏现象发生？说明产生现象的原因。
1.3.3.4 观察三角波序列和反三角波序列的时域和幅频特性
针对信号和：用8点FFT分析信号和的幅频特性，观察两者的序列形状和频谱曲线有什么异同？（注意：这时候的可以看作是经过圆周移位以后得到的）绘制两者的序列和幅频特性曲线。在的和末尾补零，用16点FFT分析这两个信号的幅频特性，观察幅频特性发生了什么变化？两个信号之间的FFT频谱还有没有相同之处？这些变化说明了什么？
1.3.3.5 * 选作内容
将信号的长度N设为63，用MatLab中randn(1,N)函数产生一个噪声信号w(n)，计算将这个噪声信号叠加到上以后新信号的频谱，观察发生的变化并记录。
在步骤2的基础上，改变参数α和f，观察在出现混淆现象和泄漏现象的时候有噪声的y(n)信号的频谱有什么变化，是否明显？
1.3.4  离散信号、系统和系统响应的分析
1.3.4.1 信号序列产生
产生如下信号序列：
（1）理想采样信号序列：对信号进行理想采样，可以得到一个理想的采样信号序列：，其中A为幅度因子，是衰减因子，是频率。T为采样周期。
（2）单位脉冲序列 
（3）矩形序列    ，其中N=10
产生如下系统单位脉冲响应序列，本实验中用到两种FIR系统：
（1）   
（2）
1.3.4.2 离散信号、系统和系统响应的分析
（1）观察信号和系统的时域和幅频特性；利用线性卷积求信号通过系统以后的响应。比较系统响应和系统的时域及幅频特性。注意它们之间有无差异，绘出图形。
（2）观察信号和系统的时域和幅频特性，利用线性卷积求系统响应。判断响应序列图形及序列非零值长度是否与理论结果一致，说出一种定性判断响应序列图形正确与否的方法（提示：）。利用序列的傅立叶变换数值计算子程序求出，观察响应序列的幅频特性。定性判断结果正确与否。改变信号的矩形宽度，使N=5，重复以上动作，观察变化，记录改变参数前后的差异。

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。