实验二离散时间系统的时域分析概论

1、武汉工程大学信号分析与处理实验一专 业：通信02班学生姓名：李瑶华学号：1304200113完成时间：2020年门月16日实验二：离散时间系统的时域分析一、实验目的1. 在时域中仿真离散时间系统，进而理解离散时间系统对输入信号或延迟信 号进行简单运算处理，生成具有所需特性的输出信号的方法。2. 仿真并理解线性与非线性、时变与时不变等离散时间系统。3. 掌握线性时不变系统的冲激响应的计算，并用计算机仿真实现。4. 仿真并理解线性时不变系统的级联、验证线性时不变系统的稳定特性。二、实验设备计算机，MATLAB语言环境。三、实验基础理论1. 系统的线性性质线性性质表现为系统满足线性叠加原理：若某一输

2、入是山个信号的加权和组 成的，则输出就是系统对这X个信号中每一个的响应的相应加权和组成的。设 x(n)和x2(n)分别作为系统的输入序列,其输出分别用(n)和y2)表示，即X (rz) = 71X (z?)J, y2 (/?) = Tx2 0?)J若满足 Tax 00 + 吋2 何=“I 刃()+a2y2 (n)则该系统服从线性叠加原理，或者称该系统为线性系统。2. 系统的时不变特性若系统的变换关系不随时间变化而变化，或者说系统的输出随输入的移位而 相应移位但形状不变，则称该系统为时不变系统(或称为移不变系统)。对时不 变系统，若)=则 Tx(n-m) = y(n-m)3. 系统的因果性系统的

3、因果性即系统的可实现性。如果系统时刻的输出取决于时刻及时刻以前的 输入，而和时刻以后的输入无关，则该系统是可实现的，是因果系统。系统具有 因果性的充分必要条件为力(比)=，刃V 4. 系统的稳定性稳定系统是指有界输入产生有界输出(BIB0)的系统。如果对于输入序列，存在一个不变的正有限值，对于所有值满足lx(n)lMvoo则称该输入序列是有界的。稳定性要求对于每个有界输入存在一个不变的正有限 值，对于所有值，输出序列满足lyS)UKvoo系统稳定的充分必要条件是系统的单位取样响应绝对可和，用公式表示8 /?(/?) I On=5. 系统的冲激响应设系统输入 x(n) = 5(z?),系统输出y

4、S)的初始状态为零，这时系统输出用力)表示，即/?(/?) = T5(77),则称h(n)为系统的单位脉冲响应。0C对于任意输入信号x(n),系统输出二八工兀“-加)为利用系统/H=-OCX0C满足叠加原理得)4)工x(加)5(-7)=工x(加)715(“-讪/?i=-xw=-x00利用系统时不变性质得到)=工兀(加MS一加)=X(n)*h(n)加=一8 上式的运算关系称为卷积运算。6. 卷积的性质1) 交换律 y() = x(ji) * 力(池)=A(n) * x(n)2) 结合律y(n) = x(n) * hx (n) * h2(n) = x(n) * hx (n) * h2 (n)= x

5、(/?)*/?2(/?)*/71(7?)= x(/?)*/iI(n)*/?2(n)3) 分配律 An)*h(n) + h2(n) = x(n)*/?(n) + x(n)*h2(/?)四、实验内容与步骤1. 离散时间系统的仿真M -11)M点因果滑动平滑系统的仿真，时域表达为y) = -xn-k) iV1 k=0通过上述时域平滑系统可实现山若干个正弦信号之和所组成的信号中滤出高 频分量。据此，可以理解M点因果滑动平滑系统。n=0 :100;sl=sin(2*pi* 03 4*n) ; s2=sin (2*pi* 00 2*n) ; x=sl+s2 ;M=input(desired length

6、of che filler =1);num=ones(1,M);y=filter(num,1,x)/M;elf;subplot(2,2,1);stem(nz si);axis(0, 100, -2,2);xla bml ( n1 ) ; ylabel ( si (n) 1 ) ; title (正弦丿了:歹【J 1);subplot(2,2,2);stem(n,s2);axis(0z100,-2,2);xlabel ( n ) ; ylabel ( fs2 (n) ) ; title ( lE 弦丿了：歹U 21 );subplot(2,2,3);stem(nz x);axis(0z100,-

7、2,2);xlabel (n); ylabel Cx(n) f); title (输入信号，);subplot(2,2,4);stem(n,y);axis(0,100,-2,2); xlabel (nf) ; ylabel ( y (n) ) ; title输出信弓)； desired length of the filter =8a)Wk 5 AD ABCX IKXWS Mr, 9 。& v - rs x-ia oc mt. Web 口】Wrwto*-2）线性与非线性离散时间系统的仿真。简单的非线性系统实例：yS） = x（n）2 _x（Z7_i）x（,2 + i）n=0 :100; x=c

8、os (2*pi* 0l *n) ; xl= x 0 0 ; x2= 0 x 0 ; :3= 0 0 x;y=x2 1 *:3; y=y (1:101);subplot (2,1,1); stem (n, x) xlabel ( n 1 ) ; ylabel ( x (n) 1 ) ; title ( 1 张入信 *) subplot (2,1,2): stem (n, y) xlabel ( n 1 ) ; ylabel ( 1 y (n) 1 ) ; title ( F金出信 号)SJ Figure 1文件(F)编辑 查吾:V)括入Q)工具(T)吏直(D)奇口 (W)书助(H) E3O05

9、4)6100输入信号1020 3040506070 * 80 0MclIdI311w输出0信号简单的线性系统实例:y(n) - 0.4y(比-1) + 0J5y(n -2) = 2.2403x() + 2.490&V(比-1) + 2.2403x(比 一 2)n=0:50;p=3;q=-4;xl=cos(2*pi*0 2*n);:2=cos (2*pi*0 5*n) ; x=p*xl+q*x2;b=2.2403 2.4908 2.2403;a=l -0.4 0.75;yl = filter(bz a,xl);y2 = filter(b,a,x2);y=filter(b,azx);yt=p*yl

10、+q*y2;d=y-yt;subplot(3,1,1)stem(nz y);ylabelCyl(n) f)；title ( Talxl(n)+a2x2n) 1);subplot(3,1,2)stem(nz yt);ylabel(1y2(n):title ( *alTxl(n)+a2Tx2(n);subplot(3,1,3);stem(n,d);xlabel(1n);ylabel(1vl(n)-y2(n) 1);title(* d):Q Figure 1文件(F) 起(V隹入(I)工具(T)裒面(D)窗口 (W)扫助(H)50Ta1x1 (n)+a2x2inT亠10152025303540455

11、0a1Tx1(n)ta2Tx2(n) fn f?4 Y J? T-5050%0501015202630n: 1 1 !Q O ? Q ? rX 102c致右二沃1 1 1 1 1 1 1(1 1;.152025303540 4fi3)时变与时不变系统的仿真。时不变系统实例：y(n)-0Ay(fi-1) + 0,75y(n-2) = 2.240女(总)+ 2.490弘(斤-1) + 2.2403x( - 2)n=0:60;m=8;p=2;q=-5;x=p*cos(2*pi*0l*n)+q*cos(2*pi*04*n);xd=zeros(1,m) x;b=22403 24908 22403;a=l

12、 -0.4 0.75;y=filter(b,a,x);yd=filter(b,az xd);d=y-yd(1+m:61+m);subplot(3,1,1)stem(n,y);ylabelCy(n) f)；title (信号 x (“)的响应 y(n) ) ; grid;subplot(3,1,2);stem(n,yd(1:61);ylabel(1y(n-m)1);title ( 1 ，zj x (n-m)的响应 y (n-m) ) ; grid;subplot(3,1,3);stem(n,d);xlabel(1n1);ylabel(波0闰 1 );title(1y(n-m)-y(n-m) 1

13、) ; grid;Q Figure 1文件(F)鯛(E)查看(V) JSA(I)工具(T)克面(D)窗二(W)：(H)口呂61|氐|毬3巴貝”|風| 目|1口信号x(n)的响应y(n)y(n-m)-y(n-m)1i11r|CU易 0 埃：tillIIIII11lli 0102030406060时变系统实例：y(n) = nx(n) + x(n-l)n=-10:10;xl=sin(0 1*pi*n);subplot(2,2,1);stem(nz xl); xlabel(nf);ylabel(fxl1);x2=sin(0 l*pi*(n-1); subplot (2,2,2) ; stem (nz

14、 :2);xlabel ( n1 ) ; ylabel ( :2 ); y=n *xl+:2;subplot(2,2,3);stem(nz y) xlabel(n);ylabel(y*)Q Figure 1 X文件(F)编辑(E)查看(V)指入(I)工具CD虫面(D)窗口 (W)稻助(H)OS A鸟叵壬F晅*： 甬Pi EQnn:,cocTo J(?0:=1 zeros (1, m-1);b=2 3 2;a=l -0.4 0.75;K=0:1:m-l;y=filter(b,az x);stem(Kz y); title C冲激响应J;xlabel(-n);ylabel(h(n);Q Figur

15、e 1 X丈件(F) 歸)兰看(V)孫入(I)(n 宾面Q)窗口 帮已 飞申玄箔物阴詳建 S匚屁冲激响应2)在实际应用中高阶因果线性时不变系统可以用低阶因果线性时不变系统级联 得到，这可简化系统的设计与实现。例如，对于三阶线性时不变系统8y(n)1 Oy(z? -1)+6y(n -2)- y(n - 3) = 16x(n -1)40x(? 2) + 16x(n3)可以用一个一阶和一个二阶系统级联实现。第一级(刃)一0.25必(刃一 1) = 2x(刃)一x(n -1)第二级儿)一)2(n T)+-52(n - 2) = X 1)2% (n一 2)用MATLAB语言编程验证系统的级联。x= 1

16、zeros (1,20);n=0:20;a=8 -10 6 -1;b=0 16 -40 16;y=filter(b,a,x);al=l -0.25;bl=2 -1;a2=l -1 0.5;b2=0 1 -2;yl = filter(bl,alzx);y2 = filter(b2 z a2z yl);d=y-y2;subplot(3,1,1);stem(n,y);ylabel(1yn) 1);title信号通过高阶LTI系统的响应);grid;subplot(3,1/2)stem(nz y2);ylabel(1y2(n) *);title (信号通过低阶LT I系统级联的响应);grid;sub

17、plot(3,1,3)stem(n,d);xlabel ( 1 n * ) ; ylabel (波幅 1 ) ; title ( 1 y (n) -y2 (n) 1 ) ; grid;Q rigure 1一X文件(F)碍旧查看(V)插入 IMCD变面(D)姦口 (Wj3 0 H 怎 c - a 目 QI信号通过高阶LTI系统的响应3. 线性时不变系统的稳定性若一个线性时不变系统的冲激响应是绝对可和，则此系统就是BIBO的稳定 系统。山此，无限冲激响应线性时不变系统稳定的必要条件是，随着输入序列点 的增加，冲激响应衰减到零。用MATLAB语言编程计算一个IIR线性时不变系统 冲激响应的绝对值的和

18、，验证稳定特性。b=0 1 2 1;a=l -0.5 一0005 0.3;subplot (3, 1,1); zplane (b, a) 由分警非母多项:式的系数画出零-极点分布图 xlabelRe* ) ;ylabel ( jlm);axis(-2 2-11);h=impz(b,a);subplot(3,1,2)stem(h);title (单位脉冲响应J ;xlabel(k*);H/w=freqz(b,a)；subplot(3,1,3)plot(w/pi,abs(H);xlabel ( *omega 1 );title (f频率响应);El rigure 1文件(F) 雋(E)- 言看(V

19、)插入 华CD宾面(D)姦口 CM 務旳(H)X0已八轩咲裁溟&*凰|匡)|QI1.:Re 单位脉冲响应X/061015202530k4. 滤波概念实验通过具体的时间系统理解信号滤波概念。如：系统 1 y(n) = 0.5x(n) + 0.27x(n 1) + 0.77x( 2)系统2y(n)O.53y(n 1) + 0A6y(n 2) = 0.45x(n) + 0.5x(n 1) + 0.45 双 n 2)对于输入信号x(n) = cos(20勿256)+ cos(200勿256)0 匡11 回阶跃响应冲激响应求和就是阶跃响应3. 系统级联实验中，三阶线性时不变系统若改用并联实现，乂该如何进

20、行? 8y(n) 10y(n 1) + 6y( 一 2) - y(n -3) = 16x(/7 -1)- 40x(/? 2) + 16x(/1 - 3)11.2vj (n) = x(n) 一 0.25x(71 -1)4.8y2(7?) 一 56儿( 一 1) = x(n) 一 x(n 一 1) + 0.5x(” 一 2)yW = 16x(”) +)+ 儿Matlab实现系统并联的程序x=l zeros (1,20);n=0:20;a=8 -10 6 -1 ;b=0 16 -40 16 ; y=f ilter (b, a, x);al=l 一025;bl=l12;a2=l -1 05;b2=4 8 -56;yl=filter(blz