西电数字信号实验第一次上机实验报告.docx

实验一：信号、系统及系统响应一、实验目的(1)熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。(2) 熟悉时域离散系统的时域特性。(3) 利用卷积方法观察分析系统的时域特性。(4) 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。二、实验原理与方法1．连续时间信号的采样采样是从连续时间信号到离散时间信号的过渡桥梁，对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生的变化以及信号内容不丢失的条件，而且有助于加深对拉氏变换、傅氏变换、z变换和序列傅氏变换之间关系的理解。对一个连续时间信号进行理想采样的过程可以表示为该信号和个周期冲激脉冲的乘积，即（1－1）其中是连续信号的理想采样，是周期冲激脉冲（1－2）它也可以用傅立叶级数表示为：（1－3）其中T为采样周期，是采样角频率。设是连续时间信号的双边拉氏变换，即有：（1－4）此时理想采样信号的拉氏变换为（1－5）作为拉氏变换的一种特例，信号理想采样的傅立叶变换（1－6）由式（1－5）和式（1－6）可知，信号理想采样后的频谱是原信号频谱的周期延拓，其延拓周期等于采样频率。根据Shannon采样定理，如果原信号是带限信号，且采样频率高于原信号最高频率分量的2倍，则采样以后不会发生频率混淆现象。在计算机处理时，不采用式（1－6）计算信号的频谱，而是利用序列的傅立叶变换计算信号的频谱，定义序列，根据Z变换的定义，可以得到序列的Z变换为：（1－7）以代替上式中的z，就可以得到序列的傅立叶变换（1－8）式（1－6）和式（1－8）具有如下关系：（1－9）由式（1－9）可知，在分析一个连续时间信号的频谱时，可以通过取样将有关的计算转化为序列傅立叶变换的计算。2．有限长序列分析一般来说，在计算机上不可能，也没有必要处理连续的曲线，通常，我们只要观察、分析在某些频率点上的值。对于长度为N的有限长序列（1－10）一般只需要在之间均匀的取M个频率点，计算这些点上的序列傅立叶变换（1－11）其中，是一个复函数，它的模就是幅频特性曲线。3．信号卷积一个线性时不变离散系统的响应y(n)可以用它的单位冲激响应h(n)和输入信号x(n)的卷积来表示：（1－12）根据傅立叶变换和Z变换的性质，与式（1－12）对应，有（1－13）（1－14）式（1－12）告诉我们可以通过对两个序列的移位、相乘、累加计算信号响应；而式（1－14）告诉我们卷积运算也可以在频域上用乘积实现。三、实验内容及步骤(1) 认真复习采样理论、离散信号与系统、线性卷积、序列的傅里叶变换及性质等有关内容，阅读本实验原理与方法。(2) 编制实验用主程序及相应子程序。①信号产生子程序，用于产生实验中要用到的下列信号序列：a. xa(t)=A*e^-at *sin(Ω0t)u(t)b. 单位脉冲序列：xb(n)=δ(n)c. 矩形序列： xc(n)=RN(n), N=10②系统单位脉冲响应序列产生子程序。本实验要用到两种FIR系统。a. ha(n)=R10(n);b. hb(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3) ③有限长序列线性卷积子程序用于完成两个给定长度的序列的卷积。可以直接调用MATLAB语言中的卷积函数conv。 conv用于两个有限长度序列的卷积，它假定两个序列都从n=0 开始。调用格式如下：y=conv (x, h)四、实验结果记录①分析采样序列的特性。a. 取采样频率fs=1 kHz,，即T=1 ms。b. 改变采样频率，fs=300 Hz，观察|X(e^jω)|的变化，并做记录(打印曲线)；进一步降低采样频率，fs=200 Hz，观察频谱混叠是否明显存在，说明原因，并记录这时的|X(e^jω)|曲线。Matlab程序如下：close all;clear all;clc;A=444.128;a=50*sqrt(2)*pi;m=50*sqrt(2)*pi;fs1=1000;fs2=300;fs3=200; T1=1/fs1;T2=1/fs2;T3=1/fs3;N=100;n=[0:N-1];x1=A*exp(-a*n*T1).*sin(m*n*T1);x2=A*exp(-a*n*T2).*sin(m*n*T2);x3=A*exp(-a*n*T3).*sin(m*n*T3);w=linspace(-pi,pi,10000);X1=x1*exp(-j*n*w);X2=x2*exp(-j*n*w);X3=x3*exp(-j*n*w); figure(1)subplot(131),stem(n,x1,LineWidth,2);xlabel(n);ylabel(x1(n));title(采样频率1KHz x1(n));axis([-1 51 -