数字信号处理matlab实验2-离散系统的时域分析.docVIP

实验2离散系统的时域分析 实验目的：加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。 实验原理：离散系统 其输入、输出关系可用以下差分方程描述： 输入信号分解为冲激信号， 。 记系统单位冲激响应 ， 则系统响应为如下的卷积计算式： 当时，h[n]是有限长度的（n：[0，M]），称系统为FIR系统；反之，称系统为IIR系统。 在MATLAB中，可以用函数y=filter(p,d,x)实现差分方程的仿真，也可以用函数 y=conv(x,h)计算卷积，用y=impz(p,d,N)求系统的冲激响应。 实验内容和要求： 1、以下程序中分别使用conv和filter函数计算h和x的卷积y和y1，运行程序，并分析y和y1是否有差别，为什么要使用x[n]补零后的x1来产生y1；具体分析当h[n]有i个值，x[n]有j个值，使用filter完成卷积功能，需要如何补零？ % Program P2_7 clf; h = [3 2 1 -2 1 0 -4 0 3]; % impulse response x = [1 -2 3 -4 3 2 1]; % input sequence y = conv(h,x); n = 0:14; subplot(2,1,1); stem(n,y); xlabel(Time index n); ylabel(Amplitude); title(Output Obtained by Convolution); grid; x1 = [x zeros(1,8)]; y1 = filter(h,1,x1); subplot(2,1,2); stem(n,y1); xlabel(Time index n); ylabel(Amplitude); title(Output Generated by Filtering); grid; 程序运行结果： 由图可看出，y与y1并无差别。 使用x[n]补零后的x1来产生y1,是因为存在边界效应，只要脉冲响应采样电部分位于输入信号采样值之外，输出就不确定，如： x 1 -2 3 -4 3 2 1 2 h 3 0 -4 0 1 -2 1 2 3 2 4 1 需变换成如下才能确定输出： x 1 -2 3 -4 3 2 1 2 0 0 0 0 0 0 0 0 0 0 0 h 3 0 -4 0 1 -2 1 2 3 2 4 1 此时n＝18。 用conv函数计算能再输入序列后自动补零，而filter函数不能。 分析： (1) h = [1 4 2 3 2 1 -2 1 0 -4 0 3]; % impulse response i=12 x = [1 -2 3 -4 3 2 1 2]; % input sequence j=8 n = 0:17; x1 = [x zeros(1,10)]; %补十个零值 y1 = filter(h,1,x1); stem(n,y1); xlabel(Time index n); ylabel(Amplitude); title(Output Generated by Filtering); grid （1）图 (2) h = [1 4 2 3 2 1 -2 1 0 -4 0 3]; % impulse response i=12 x = [1 -2 3 -4 3 2 1 2]; % input sequence j=8 n = 0:18; x1 = [x zeros(1,11)]; %补11个零值 y1 = filter(h,1,x1); stem(n,y1); xlabel(Time index n); ylabel(Amplitude); title(Output Generated by Filtering); grid （2）图 (3) h = [1 4 2 3 2 1 -2 1 0 -4 0 3]; % impulse response i=12 x = [1 -2 3 -4 3 2 1 2]; % input sequence j=8 n = 0:30; x1 = [x zeros(1,23)]; %补23个零值 y1 = filter(h,1,x1); stem(n,y1); xlabel(Time index n); ylabel(Amplitude); title(Output Generated by Filtering); grid; （3）图 对照（1）、（2）图，当n＝18时两图有区别，(2)图能完全卷积，当补零数