信号与系统实验二.docxVIP

信号与系统实验报告二离散时间系统的实验分析实验报告学院：计算机与信息工程学院专业：通信工程班级：2012级计算机与信息技术学院设计性实验报告专业：通信工程 年级/班级：2012级 2013—2014学年第二学期课程名称信号与系统指导教师学号姓名实验地点计科楼216实验时间2014-4-15项目名称离散时间系统的时域分析实验类型设计性实验目的1. 通过matlab仿真一些简单的离散时间系统，并研究它们的时域特性。2. 掌握利用matlab工具箱求解LTI 系统的单位冲激响应。实验仪器或设备装有MATLAB软件设备的计算机一台实验原理离散时间系统的时域特性1.1 线性定义对离散时间系统，若y1[ n] 和 y2 [ n]分别是输入序列x1 [ n] 和 x 2 [ n]的响应，若输入x3[n] = x1[n] + x2 [n]的输出响应为y3[n] = y1[n] + y2 [n]该性质对任意输入 x1 [ n] 和 x 2 [ n] 都成立，称为叠加性。若输入x4[n] = αx1[n]的输出响应为y4 [n] = αy1[n]该性质对任意输入 x1 [ n] 和任意常量α 都成立，称为齐次性。若同时满足叠加性和齐次性，即输入x[n] = αx1[n] + βx2 [n]的输出响应为y[n] = αy1[n] + β y 2 [n]对任意常量α 和 β 以及任意输入 x1 [ n] 和 x 2 [ n] 都成立，则称为线性。反之称为非线性。1.2 时不变定义对于离散时间系统，若 y1[ n] 是输入序列 x1 [ n] 的响应，若输入x[n] = x1 [n ?n0 ]的输出响应为y[n] = y1[n ?n0 ]对任意整数 n0 和任意输入 x1 [ n] 及其对应输出都成立，则称为时不变。反之称为时变。1.3 LTI 系统线性时不变(LTI)系统既满足线性特性，又满足时不变特性。1.4 单位冲激响应若输入信号为单位冲激信号δ [n]，离散时间系统的响应称为单位冲激响应h[n]。1.5 单位阶跃响应若输入信号为单位阶跃信号u[n]，离散时间系统的响应称为单位阶跃响应s[n]。2.相关基本 matlab函数2.1 filter可获得系统响应。若系统表达式为=x 为输入向量，可令num=[p0 p1…pn],den=[d0 d2…dm]，则 y=filter(num,den,x)产生的输出向量 y 的长度与输入向量 x 的长度相同，且初始值为0。也可用 y=filter(num,den,x,ic)计算系统输出，其中 ic是初始值。2.2 impzy=impz(num,den,N)可计算 LTI 离散时间系统的单位冲激响应的前 N 个样本。四、 实验内容1、离散时间系统的时域分析1.1 线性与非线性系统1.1.1假定系统为 y[n] ?0.4y[n?1] = 2.24x[n] + 2.49x[n?1]输入三个不同的输入序列 x1[n]、x2[n] 和x[n]，计算并求出相应的三个输出，并判断是否线性: x[n] = αx1[n] + βx2[n]程序代码如下：n=0:40;a=2;b=-3;x1=cos(2*pi*0.1*n);x2=sin(2*pi*0.4*n);x=a*x1+b*x2;num=[2.24 2.49];den=[1 -0.4];y1=filter(num,den,x1);y2=filter(num,den,x2);y=filter(num,den,x);yt=a*y1+b*y2;d=y-yt;%计算差值输出d[n]subplot(3,1,1)stem(n,y);ylabel(振幅);subplot(3,1,2)stem(n,yt);ylabel(振幅);subplot(3,1,3)stem(n,d);ylabel(振幅);title(差信号)运行结果为：代码分析：x1=cos(2*pi*0.1*n);x2=sin(2*pi*0.4*n);x1，x2为频率不同的正弦输入信号，x=a*x1+b*x2;为信号x1和x2的线性叠加。num=[2.24 2.49];den=[1 -0.4];表示系统y[ n] ?0.4 y[ n ?1] = 2.24 x[ n] + 2.49 x[ n ?1]y1=filter(num,den,x1);y2=filter(num,den,x2);y1，y2为使用filter函数得到信号x1与x2通过该系统的系统响应。y=filter(num,den,x);y为叠加后的信号x通过该系统的系统响应。yt=a*y1+b*y2;yt为响应y1和y2的线性叠加由图可知，输出信号y与yt基本相同，差信号d的幅值特别小，基本可看做y=yt，所以该系统为线性系统。1.1.2假定另一个系统