唐佳实验二详解.doc

《I》 实验报告 名： 唐佳 号：1508103 专 业：通信工程8节 指导老师 杭州电子科技大学 通信工程学院 一、实验目的 1. 通过matlab 仿真一些简单的离散时间系统，并研究它们的时域特性。 2. 掌握利用matlab 工具箱求解LTI 系统的单位冲激响应 实验内容 1. 离散时间系统的时域特性 1.1 线性定义 对离散时间系统，若 y1[ n] 和 y2 [ n]分别是输入序列 x1 [ n] 和 x2 [ 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] + βx 2 [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]。 三、实验过程 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(差信号) 运行结果 假定另一个系统为y[n] = x[n]+ 3.2x[n ? 2]，修改以上程序，通过绘出的图形判断该系统是否线性系统 程序代码如下 clear all; 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=[1 0 3.2]; den=[1 0 0]; 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(差信号); 运行结果 1.2 时变与时不变系统 根据1.1中系统，产生两个不同的输入序列x[n]和x[n-D]，根据输出判断是否时不变系统。 程序代码如下： clear all; n=0:40; a=2;b=-3;D=10; x=cos(2*pi*0.1*n); xd=[zeros(1,D) x]; num=[2.24 2.49]; den=[1 -0.4]; y=filter(num,den,x); yd=filter(num,den,xd); subplot(3,1,1); stem(n,