(数字信号处理实验一.docxVIP

一、实验目的1．? 熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对采样定理的理解。2．? 熟悉离散信号何系统的时域特性。3．? 熟悉线性卷积的计算编程方法：利用卷积的方法，观察、分析系统响应的时域特性。4．? 掌握序列傅氏变换的计算机实现方法，利用序列的傅氏变换对离散信号、系统及系统响应进行频域分析。二、实验原理（一）连续时间信号的采样采样是从连续信号到离散时间信号的过渡桥梁，对采样过程的研究不仅可以了解采样前后信号时域何频域特性发生的变化以及信号内容不丢失的条件，而且有助于加深对拉氏变换、傅氏变换、Z变换和序列傅氏变换之间关系的理解。对一个连续时间信号进行理想采样的过程可以表示为信号与一个周期冲激脉冲的乘积，即：（1-1）其中，是连续信号的理想采样，是周期冲激脉冲（1-2）此时，采样信号的拉氏变换可以表示为：,其中是信号的拉氏变换作为拉氏变换的一种特例，信号理想采样的傅立叶变换为：（1-3）由（1-3）式可知，信号理想采样后的频谱式原来信号频谱的周期延拓，其延拓周期等于采样频率。根据香农定理，如果原信号是带限信号，且采样频率高于原信号最高频率的2倍，则采样后的离散序列不会发生频谱混叠现象。在计算机处理时，不采用时计算信号的频谱，而是利用序列的傅立叶变换计算信号的频谱，定义序列：（1-4）根据Z变换的定义，可以得到序列的Z变换为：（1-5）以代替上式中的Z，就可以得到序列的傅立叶变换：（1-6）式（1-3）和式（1-6）具有如下关系：（1-7）由式可知，在分析一个连续时间信号的频谱时，可以通过取样将有关的计算转化为序列傅立叶变换的计算。（二）有限长序列的分析一般来说，在计算机上不可能、也不必要处理连续的曲线，通常，我们的做饭时只要观察、分析在某些频率点上的值。对于长度为N的有限长序列：（1-8）一般只需要在之间均匀的取M个频率点，计算这些点上的序列傅立叶变换：（1-9）其中，，。是一个复函数，它的模就是幅频特性曲线。（三）信号卷积一个线性时不变离散系统的响应可以用它的单位冲激响应和输入信号的卷积来表示：（1-10）根据傅立叶变换和 Z变换的性质，应该有：（1-11）（1-12）上两式告诉我们：可以通过对两个序列的移位、相乘和累加计算信号响应；卷积运算可以在频域通过乘积实现。三、实验内容及步骤(一) 编制设计用主程序及相应子程序1．信号产生子程序，包括：（1）理想采样信号序列：对信号进行理想采样，可以得到一个理想的采样信号序列，其中为幅度因子，是衰减因子，是频率，为采样周期。根据设计内容的需要，这些参量请设计为在程序运行过程中输入。（2）单位脉冲序列（3）矩形序列, 其中取：N=102．系统单位冲激响应序列产生子程序，本实验中用到两种FIR系统：（1）（2）3．有限长序列线性卷积子程序，用于计算两个给定长度（分别为M和N）的序列的卷积，输出序列长度为L=M+N-1。M文件：1．function [xa] = lixiangcaiyangxinhaoxulie(A,a,w0,T)globalxa;A=input(A=),a=input(a=),w0=input(w0=),T=input(T=);n=0:50;xa=A*exp(-a*n*T).*sin(w0*n*T);end2.function [xb]=danweimaichongxulie()globalxb;n=0:50;xb=[1 zeros(1,50)];3.function [xc] = juxingxulie(N)globalxc;N=input(N=);xc=[ones(1,N) zeros(1,41)];end4.function [ha] = FIR1()global ha;ha=[ones(1,10) zeros(1,41)];end5.function [hb] = FIR2()globalhb;n=1:50;hb=zeros(1,50);hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1;end（二）上机设计内容在编制以上各部分程序以后，编制主程序调用各个功能模块实现对信号、系统和系统响应的时域和频域分析，完成以下设计内容。1．分析理想采样信号序列的特性产生理想采样信号序列，使：（1）首先选用采样频率为1000Hz，T=1/1000，观察所得理想采样信号的幅频特性，在折叠频率以内和给定的理想幅频特性无明显差异，并作记录。lixiangcaiyangxinhaoxulieA=44.128A = 44.1280a=50*sqrt(2.0)*pia = 222.1441w0=50*sqrt(2.0)*piw0 = 222.1441T=0.001subplot(1,1,1);stem(n,xa);title(理想采样信号序列)