云南师范大学dsp数字信号处理试验报告DOC.docxVIP

本科学生实验报告 学号 114090395 姓名 李华福 学院 物电学院 专业、班级 11 电子 实验课程名称 数字信号处理（实验） 教师及职称 李宏宁老师 开课学期 2013 至 2014 学年 下 学期 填报时间 2014 年 5 月 7 日 云南师范大学教务处编印 序号 8 名称 离散系 分析 2014年5月7日 室 数字信号 理 室（ 313） 一． 1． 目的 深刻理解离散 系 的系 函数在分析离散系 的 域特性、 率特性及 定性中的重要作用及意 ，熟 掌握利用 MATLAB分析离散系 的 域响 、 率响 和零极点的方法。掌握利用 DTFT和 DFT确定系 特性的原理和方法。 2． 原理、 流程或装置示意 原理： MATLAB提供了 多可以用于分析 性 不 离散系 的函数， 主要包含有系 函数、系 域响 、系 率响 等分析函数。 （ 1）离散系 的 域响 描述离散 LTI 系 的 入 x[k] 与 出 y[k] 关系的差分方程 ： y[k ] y[ k 1] aN 1 y[ k (N 1)] an y[k N ] (2.2.1) b0 x[ k] b1 x[ k 1] bM 1 x[ k (M 1)] bM x[k M ] 上式两 同取 z 可得离散系 的系 函数 H(z) H (z) Y(z) b0 b1z 1 bM 1z ( M 1) bM z M b( z) X ( z) 1 a1z 1 aN 1z ( N 1) aN z N (2.2.2) a(z) 使用向量 b=[ b0 ,b1, ,bM 1,bM ] 和向量 a [ a0 , a1, , aN 1,aN ] 分别保存分子多 项式和分母多项式的系数，这些系数都是从 z0 按 z 的降幂排列。 离散系 的 位冲激响 h[k] 的 算 [h,k]=impz(b,a) 算系 的 位脉冲响 的相 的 向量，其 中 b 和 a 分 是系 函数 H(z) 的分子多 式和分母多 式的系数矩 ； 也可以 写 h=impz(b,a) 。 [h,k]=impz(b,a,n)算 n 点 位 脉 冲 响 ； 也 可 以 写 h=impz(b,a,n) 。 Impz(b,a) 自 制 位脉冲响 的 形。 2）离散系 响 y[k] 的 算 y=filter(b,a,x) 算系 在 入 x 作用下的零状 响 y。 y=filter(b,a,x ，zi) 算系 在 入 x 和初始状 作用下的完全响 y。zi 是由系 的初始状 filtic 函数 而得到的初始条件： zi=filtic(b,a,Y0),其 中 Y0系的 初 始 状， Y0=[y[-1],y[-2],y[-3], ?] 。 也可利用 lsim 函数求解离散系 的系 响 ， 函数在 系 分 析 中介 ，在此不再重复。 [ 例 2.2.1] 已知系统函数为 H ( z) 1 2z 1 2 ，求 0.4z 1 0.12z 1 离散系统的单位脉冲响应 h[k]; 输入 x[k]=u[k] ，求系统的零状态响应 y[k]; (3) 输入就 x[k]=u[k] ，初始条件 y[-1]=1,y[-2]=2 ，求系统的完全响应 y[k] 。 【解】 1） 计算前 40 个点的单位脉冲响应 h[k] N=40; a=[1,0.4,-0.12];b=[1,2]; y=impz(b,a,N); stem(y) xlabel( k );title( h[k] ) 运行结果如下图 2.1.1 所示。 h[k] 2 1.5 1 0.5 0 -0.5 -1 0 5 10 15 20 25 30 35 40 k 图2.2.1 系统的单位脉冲响应 2） 计算前 100 个点的零状态响应 y[k] N=40; a=[1,0.4,-0.12];b=[1,2]; x=ones(1,N); y=filter(b,a,x) （ 3） 计算前 100 个时刻的完全响应 y[k] Filter(b,a,x,zi) 中的初始值 zi 不是 y[-1]=1,y[-2]=2 ，它可以由 filtic 函数计算。 N=100; a=[1,0.4,-0.12];b=[1,2]; x=ones(1,N); zi=filtic(b,a,[1,2]); y=filter(b,a,x,zi); （ 2）离散系统的系统函数零极点分析 LTI 离散系统的系统函数 H(z) 可以表示为零极点形式，即 Y( z) b0 b1 z 1 bM-1 z ( M 1) bM z M H ( z) 1 a z 1 a z ( N