实验二离散系统分析信号分析与处理实验电子教案.pptVIP

实验二离散系统分析一、实验目的深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频域特性以及稳定性中的重要作用及意义，熟练掌握利用MATLAB分析离散系统的时域响应、频响特性和零极点的方法。掌握利用DTFT和DFT确定系统特性的原理和方法。实验二离散系统分析实验原理MATLAB提供了许多可用于分析线性时不变连续系统的函数，主要包含有系统函数、系统时域响应、系统频域响应等分析函数。离散时间LTI系统可用如下的线性常系数差分方程来描述：已知输入信号x[k]以及系统初始状态y[-1]，y[-2]...，就可以求出系统的响应。MATLAB提供了计算差分方程的数值解的函数，可以计算上述差分方程描述的离散系统的单位脉冲响应、零状态响应和完全响应等。1.离散系统的时域响应实验二离散系统分析0102在调用MATLAB函数时，需要利用描述该离散系统的系数函数。对差分方程进行Z变换即可得系统函数：在MATLAB中可使用向量a和向量b分别保存分母多项式和分子多项式的系数：这些系数均从z0按z的降幂排列。1.离散系统的时域响应实验二离散系统分析也可简写为：h=impz(b,a,n)。离散系统的单位冲激响应h[k]的计算1)[h,k]=impz(b,a,n)：计算n点单位脉冲响应h[k]；[h,k]=impz(b,a)：计算系统的单位脉冲响应h[k]和相应的时间向量k；也可简写为：h=impz(b,a)。其中：impz(b,a)：绘制单位脉冲响应h[k]的图形。1.离散系统的时域响应实验二离散系统分析1.离散系统的时域响应实验二离散系统分析离散系统响应y[k]的计算y=filter(b,a,x)：计算系统在输入x作用下的零状态响应y[k]；y=filter(b,a,x,zi)：计算系统在输入x和初始状态作用下的完全响应y[k]。zi是由系统的初始状态经过filtic函数转换而得到的初始条件：zi=filtic(b,a,Y0)，Y0为系统的初始状态，Y0=[y[-1]，y[-2]，y[-3]，...]。实验二离散系统分析已知系统函数为，计算(1)计算前40个点的单位脉冲响应N=40;a=[1,0.4,-0.12];b=[1,2];y=impz(b,a,N);stem(y)xlabel(k);title(h[k])(1)离散系统的单位脉冲响应；(2)输入x[k]＝u[k]，求系统的零状态响应y[k]；(3)输入x[k]＝u[k]，初始条件y[-1]＝1，y[-2]＝2，求系统的完全响应y[k]。实验二离散系统分析(3)计算前100个时刻的完全响应filter(b,a,x,zi)中的初始值zi不是y[-1]＝1，y[-2]＝2，它可以由filtic函数计算。N=100;b=[1,2];a=[1,0.4,-0.12];x=ones(1,N);zi=filtic(b,a,[1,2]);y=filter(b,a,x,zi);(2)计算前100个点的零状态响应N=100;b=[1,2];a=[1,0.4,-0.12];x=ones(1,N);y=filter(b,a,x)离散系统的系统函数零极点分析离散LTI系统的系统函数H(z)可以表示为零极点形式：使用MATLAB提供的roots函数计算离散系统的零极点；使用zplane函数绘制离散系统的零极点分布图。注意：在利用这些函数时，要求H(z)的分子多项式和分母多项式的系数的个数相等，若不等则需要补零。实验二离散系统分析实验二离散系统分析”已知系统函数为计算该系统函数的零极点，并画出系统函数零极点分布图。b=[1,2,0];a=[1,0.4,-0.12];z=roots(b)p=roots(a)zplane(b,a)01020304当离散因果LTI系统的系统函数H(z)的极点全部位于z平面单位圆内时，系统的频率响应可由H(z)求出，即计算系统的n点频率响应H，w为频率点向量。05freqz(b,a)：绘制频率响应曲线。[H,w]=freqz(b,a,n)：H=freqz(b,a,w)：计算系统在指定频率点向量