1、PAGE 信號(hào)、系統(tǒng)與信號(hào)處理實(shí)驗(yàn) = 1 * ROMAN I實(shí)驗(yàn)報(bào)告 實(shí)驗(yàn)名稱：離散時(shí)間系統(tǒng)的時(shí)域分析 姓 名： 學(xué) 號(hào)： 專 業(yè)： 通信工程 實(shí)驗(yàn)時(shí)間 杭州電子科技大學(xué)通信工程學(xué)院一、實(shí)驗(yàn)?zāi)康?.通過(guò) matlab 仿真一些簡(jiǎn)單的離散時(shí)間系統(tǒng)，并研究它們的時(shí)域特性。 2.掌握利用 matlab 工具箱求解 LTI 系統(tǒng)的單位沖激響應(yīng)。實(shí)驗(yàn)內(nèi)容1、離散時(shí)間系統(tǒng)的時(shí)域分析 1.1 線性與非線性系統(tǒng) 假定系統(tǒng)為 yn-0.4yn-1=2.24xn+2.49xn-1 (2.9) 輸入三個(gè)不同的輸入序列 x1n、x2n和，計(jì)算并求出相應(yīng)的三個(gè)輸出，并判斷是否線性。xn=ax1n+bx2n clear

2、 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=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;%計(jì)算差值輸出 dn subplot(3,1,1) stem(n,y); ylabel(振幅); subplot(3,1,2) stem(n,yt); ylabel(振幅); subplot(3,1,3) stem(n,d)

3、; ylabel(振幅); title(差信號(hào)) (1)假定另一個(gè)系統(tǒng)為 yn=xn+3.2xn-2，修改以上程序，通過(guò)繪出的圖形判斷該系統(tǒng)是否線性系統(tǒng)。 1.2 時(shí)變與時(shí)不變系統(tǒng) 根據(jù)(2.9)的系統(tǒng)，產(chǎn)生兩個(gè)不同的輸入序列 xn和 xn-D，根據(jù)輸出判斷是否時(shí)不變系統(tǒng)。 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); d=y-yd(1+D:41+D);%計(jì)算差值輸dn

4、 subplot(3,1,1) stem(n,y); ylabel(振幅); title(輸出 yn);grid; subplot(3,1,2) stem(n,yd(1:41); ylabel(振幅); title(由于延時(shí)輸入 xn,num2str(D),的輸出);grid; subplot(3,1,3) stem(n,d); ylabel(振幅); title(差信號(hào));grid; (1)yn和 ydn有什么關(guān)系？延時(shí)是由哪個(gè)參數(shù)控制的？ (2)該系統(tǒng)是時(shí)不變系統(tǒng)嗎？ 2、線性時(shí)不變系統(tǒng)的單位沖激響應(yīng) 根據(jù)(2.9)的系統(tǒng)，計(jì)算并繪出單位沖激響應(yīng)。 clear all; N=40;num=

5、2.24 2.49; den=1 -0.4; y=impz(num,den,N); stem(y); xlabel(時(shí)間序號(hào));ylabel(振幅);title(沖激響應(yīng));grid; (1)修改以上程序，產(chǎn)生如下 LTI 系統(tǒng)的單位沖激響應(yīng)前 45 個(gè)樣本： yn+0.71yn-1-0.46yn-2-0.62yn-3=0.9xn-0.45xn-1+0.35xn-2+0.002xn-3 三、實(shí)驗(yàn)過(guò)程及實(shí)驗(yàn)結(jié)果 1.1clear 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=2.24 2.49;

6、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;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(差信號(hào))線性僅在極小的范圍內(nèi)有誤差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;

7、num=1 0 3.2;den=1;y1=filter(num,den,x1);y2=filter(num,den,x2);y=filter(num,den,x);yt=a*y1+b*y2;d=y-yt;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(差信號(hào))線性僅在極小的范圍內(nèi)有誤差1.2clear all;n=0:40;a=2;b=-3;D=10;x=cos(2*pi*0.1*n);xd=zeros(1,D) x;n

8、um=2.24 2.49;den=1 -0.4;y=filter(num,den,x);yd=filter(num,den,xd);d=y-yd(1+D:41+D);subplot(3,1,1)stem(n,y);ylabel(振幅);title(輸出yn);grid;subplot(3,1,2)stem(n,yd(1:41);ylabel(振幅);title(由于延時(shí)輸入xn-,num2str(D),的輸出);grid;subplot(3,1,3)stem(n,d);ylabel(振幅);title(差信號(hào));grid;（1）ydn由yn向右平移10個(gè)單位得到，延時(shí)由參數(shù)D控制。（2）該系統(tǒng)

9、是時(shí)不變系統(tǒng)。yn經(jīng)過(guò)時(shí)移后得到y(tǒng)dn,差信號(hào)為0，所以是時(shí)不變系統(tǒng)。2clear all;N=45;num=0.9 -0.45 0.35 0.002;den=1 0.71 -0.46 -0.62;y=impz(num,den,N);stem(y);xlabel(時(shí)間信號(hào));ylabel(振幅);title(沖激響應(yīng));grid;N=40;num=0.9 -0.45 0.35 0.002;den=1 0.71 -0.46 -0.62;y2=impz(num,den,N);x1=1,zeros(1,N-1); y1=filter(num,den,x1);y3=y2;m=y1-y3stem(m);xlabel