实验七 信号与系统.doc

学号： 14132501702 信号与系统实验报告 学生姓名 班 级 院 部 物理与电子学院 专 业 电子科学与技术 任课老师 王晓明 指导老师 王晓明 二0一四—— 二0一五 学年第 二 学期 实验项目名称： 离散线性时不变系统的分析 实验成绩： 实验日期： 2015.6.27 实验室： 6404 一、实验目的 深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频域特性以及稳定性中的重要作用及意义，熟练掌握利用MATLAB分析离散系统的时域响应、频响特性和零极点的方法。 二、相关知识点 1．离散时间线性时不变系统的时域响应（单位冲击响应、单位阶跃响应、零输入响应、零状态响应和全响应）。 2．离散时间线性时不变系统的频率特性（幅度谱和相位谱）。 3．离散时间线性时不变系统的零极点及稳定性。 三、实验原理、方法和手段 1. 离散系统的时域响应 离散时间LTI系统可用如下的线性常系数差分方程来描述： 已知输入信号x[k]以及系统初始状态y[-1]，y[-2] ...，就可以求出系统的响应。MATLAB提供了计算差分方程的数值解的函数，可以计算上述差分方程描述的离散系统的单位脉冲响应、零状态响应和完全响应等。 在调用MATLAB函数时，需要利用描述该离散系统的系数函数。对差分方程进行Z变换即可得系统函数： 在MATLAB中可用向量a和向量b分别保存分母多项式和分子多项式的系数： 这些系数均从z0按z的降幂排列。 离散系统的单位冲激响应h[k]的计算 [h,k] = impz(b, a)：计算系统的单位脉冲响应h[k]和相应的时间向量k；也可简写为：h = impz(b, a)。 [h,k] = impz(b, a, n)：计算n点单位脉冲响应h[k]； 也可简写为：h = impz(b, a, n)。 impz(b, a)：绘制单位脉冲响应h[k]的图形。 离散系统零状态响应y[k]的计算 y = filter(b, a, x)：计算系统在输入x作用下的零状态响应y[k]； （3）离散系统全响应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) 离散系统的单位脉冲响应； 2) 输入x[k]＝ u[k]，求系统的零状态响应y[k]； 3) 输入x[k]＝ u[k]，初始条件y[-1]＝ 1， y[-2]＝ 2，求系统的完全响应y[k]。 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]) 2) 计算前100个点的零状态响应N=100; b=[1,2]; a=[1,0.4,-0.12]; x=ones(1,N); y=filter(b,a,x) 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．离散系统的频率响应 当离散因果LTI系统的系统函数H(z)的极点全部位于z平面单位圆内时，系统的频率响应可由H(z)求出，即 [H, w]=freqz(b, a, n)；%计算系统的n点频率响应H，w为频率点向量。 H=freqz(b, a, w)； %计算系统在指定频率点向量w上的频响； freqz(b,a)；%绘制频率响应曲线。 abs(H)；%幅度响应函数，angle(H)；%相位响应。 其中：b和a分别为系统函数H(z)的分子分母系数矩阵； 已知某离散系统的系统函数为： 分析系统的幅频特性。 b=[1,1]; a=[1,-1,0.5]; [H,w]=freqz(b,a); plot(w,abs(H)) xlabel(Freque