实验一、数字信号处理基础.docVIP

实验一：数字信号处理基础 一、实验目的 1、熟悉并掌握离散系统的差分方程表示法； 2、加深对冲激响应和卷积分析方法的理解； 3、加深对离散信号的DFT的理解； 4、熟悉离散系统的频率响应分析方法、加深对零、极点分布的概念理解； 5、掌握Matlab常用函数的使用方法。 二、实验原理 1. LTI系统的表示 在时域中，离散时间系统对输入信号或者延迟信号进行运算处理，生成具有所需特性的输出信号，具体框图如下: 其输入、输出关系可用以下差分方程描述： 输入信号分解为冲激信号, 记系统单位冲激响应，则系统响应为如下的卷积计算式： 当时，h[n]是有限长度的(n=1:M)，称系统为FIR 系统；反之，称系统为IIR系统。 2. 离散傅里叶变换(DFT)的定义 N点序列的DFT和IDFT变换定义式如下： , , 利用旋转因子具有周期性，可以得到快速算法（FFT）。 3. 离散时间系统的变换域分析方法 离散系统的时域方程为 其变换域分析方法如下： 频域: 系统的频率响应为: Z域: 系统的转移函数为: 分解因式: ， 其中和称为零、极点。 三、预习要求 1、在MATLAB 中，熟悉利用函数y = filter(b,a,x)实现差分方程的使用方法； 2、在MATLAB 中，熟悉用函数y=conv(x,h)计算卷积的使用方法； 3、在MATLAB 中，熟悉用函数y=impz(b,a)求系统冲激响应的使用方法。 4、在MATLAB 中，熟悉用函数y=stepz(b,a)求系统冲激响应的使用方法。 5、在MATLAB中，函数和计算N点序列的DFT正、反变换的使用方法。 6、在MATLAB中，熟悉函数tf2zp、zplane、freqz、residuez、zp2sos的使用。 函数的使用方法： 1、filter:这是一个一维数字滤波器；也可以用作卷积运算，结果与conv函数一样。 y = filter(A,B,X)，x为待处理序列，H=B/A为一个系统函数，y的结果是H与X相乘后作傅氏反变换的结果，如果A=1，则y为两个序列的卷积。 2、conv:求卷积函数，y=conv(x,h)，y是x与h卷积的结果。 3、impz:求系统的冲激响应。 [h,t]=impz(b,a):b、a分别为系统的传递函数的分子和分母系数向量。返回系统(b,a)的冲激响应h和相应的时间轴向量。 4、stepz:求系统的单位阶跃响应。 [s,n]=stepz(b,a):求解离散系统的单位阶跃响应，其中b、a为向量，n为时间序列 [s,n]=stepz(b,a,N)：求解离散系统的单位阶跃响应，N为采样点数 stepz(b,a)：在当前窗口，用stem(n,s)绘出图形。 5、求离散时间序列的傅立叶变换和反变换：X=fft(x,N)，x=ifft(X,N) 其中，X为x的DFT变换，而x则是X的傅立叶反变换，N为计算的点数，可以能过上面两个函数实现DFT。 6、tf2ss：由传递函数到状态空间的转换。 [A,B,C,D]=tf2ss(NUM,DEN),NUM,DEM为向量 H(s)=NUM(s)/DEN(s)到x=Ax+Bu,y=Cx+Du 7、zplane(b,a)：绘制由行向量b和a构成的系统函数的零极点分布图； zplane(z,p)：绘制由列向量z确定的零点、列向量p确定的极点构成的零极点分布图。 8、freqz：求离散时间系统的频率响应函数； [h,w]=freqz(b,a,n)。可以得到数字滤波器的n点复频响应值，这n个点均匀地分布这[0,pi]上，并将这n个频点的频率记录在w中，相应的频响值记录在h中，n缺省时取512点。 [h,f]=freqz(b,a,n,Fs)，用于对H(jw)在[0,Fs/2]上等间隔采样n点，采样点频率及相应频响值分别记录在f和h中，由用户指定Fs(以Hz为单位)的值。 9、[r,p,c]=residuez(b,a)：把b(z)/a(z)展开成部分分式；? [b,a]=residuez(r,p,c):根据部分分式的r、p、c数组，返回有理多项式。 10、zp2sos:变系统零极点形式为二阶分割形式 [SOS,G] = ZP2SOS(Z,P,K)； SOS为矩阵形式，G为增益。 ……………………………………………………………………………………………….. 四、实验内容 1. Matlab函数conv和filter的使用 1.1 实验要求 以下程序中分别使用conv和filter函数计算h和x的卷积y和y1，运行程序，并分析y和y1是否有差别，为什么要使用x[n]补零后的x1来产生y1；具体分析当h[n]有M个值，x[n]有M个值，使用filter完成卷积功能，需要如何补零？ % Program P2_7 clf; h = [3