数字信号实验一 信号、系统及系统响应.docVIP

实验一 信号、系统及系统响应 一、实验目的 1.熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对采样定理的理解。 2.熟悉离散信号和系统的时域特性。 3.熟悉现行卷积的计算编程方法，利用卷积的方法，观察、分析系统响应的时域特性。 4.掌握序列傅氏变换的计算机实现方法，利用序列的傅氏变换对离散信号、系统及系统响应进行频域分析。 二、实验原理 (一) 连续时间信号的采样 采样时从连续时间信号到离散时间信号的过渡桥梁，对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生的变化以及信号内容不丢失的条件，而且有助于加深对拉氏变换、傅氏变换、Z变换和序列傅氏变换之间的关系。 对一个连续时间信号进行理想采样的过程可以表示为该信号的一个周期冲激脉冲的乘积，即 其中，是连续信号的理想采样，M（t）是周期冲激脉冲 根据Shannon取样定理，如果原信号时带限信号，且采样频率高于原信号最高频率分量2倍，则采样以后不会发生频谱混淆现象。 在计算机处理时是利用序列的傅里叶变换计算信号的频谱，在分析一个连续时间信号的频谱时，可以通过取样将有关的计算转化为序列福利叶变换的计算。 (二) 有限长序列分析 一般来说，在计算机上不可能，也不必要处理连续的曲线X(ejω)，通常，我们只要观察、分析X(ejω)在某些频率点上的值。对于长度为N的有限长序列： 一般只需要在0-2π之间均匀地取M个频率点，计算这些点上的序列傅里叶变换： 其中=2πk/M，k=0，1……，M-1。X()是一个复函数，它的模就是幅频特性曲线。 (三) 信号卷积 一个线性时不变离散系统的响应y(n)可以用它的单位冲击响应h(n)和输入信号x(n)的卷积来表示： y(n)=x(n)*h(n)= (1) 根据傅里叶变换和Z变换的性质，与上式对应应该有： Y(z)=X(z)H(z) (2) Y()= X()H() (3) 式(1)告诉我们可以通过对两个序列的移位、相乘、累加计算信号响应；而(3)式告诉我们卷积运算也可以在频域上用乘积实现。 三、实验内容及步骤 (一) 编制实验用主程序及相应子程序 1.信号产生子程序，包括： (1)理想采样信号序列：对信号进行理想采样，可以得到一个理想的采样信号序列：，其中，A为幅度因子，α是衰减因子，是频率，T为采样周期。 (2)单位脉冲序列： (3)矩形序列：，其中N=10 2.系统单位脉冲响应序列产生子程序，本实验中用到两种FIR系统： (1) (2) 3.有限长序列线性卷积子程序，用于计算两个给定长度(分别是M和N)的序列的卷积，输出序列长度为L=N+M-1。 (二) 上机实验内容 1.分析理想采样信号序列的特性。 1.理想采样信号序列 n=0:50; A=444.128;a=50*sqrt(2)*pi; w0=50*sqrt(2)*pi; T=1/1000; x=A*exp(-a*n*T).*sin(w0*n*T); x1=fft(x); subplot(3,1,1) stem(n,x1) clear T T=1/300; x=A*exp(-a*n*T).*sin(w0*n*T); x2=fft(x); subplot(3,1,2) stem(n,x2) clear T T=1/200; x=A*exp(-a*n*T).*sin(w0*n*T); x3=fft(x); subplot(3,1,3) stem(n,x3) 2.离散信号，系统和系统响应的分析 （1） n=1:50; hb=zeros(1,50); hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1; close all;subplot(3,1,1);stem(hb); m=1:50;T=0.001; A=444.128;a=50*sqrt(2.0)*pi; w0=50*sqrt(2.0)*pi; x=A*exp(-a*m*T).*sin(w0*m*T); subplot(3,1,2);stem(x); y=conv(x,hb); subplot(3,1,3);stem(y); （2） n=0:50;k=-25:25; xc=[ones(1,10) zeros(1,41)]; subplot(3,1,1);stem(n,xc); axis([0 50 0 1]); X=xc*(exp(-j*pi/25)).^(n*k); magX=abs(X); sub