數(shù)字信號處理系統(tǒng)響應(yīng)及系統(tǒng)穩(wěn)定性實(shí)驗(yàn)報(bào)告

1、word上海電力學(xué)院實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)課程名稱：數(shù)字信號處理實(shí)驗(yàn)工程名稱：系統(tǒng)響應(yīng)及系統(tǒng)穩(wěn)定性班 級：2022071姓 名：王濤 學(xué) 號：20221938成 績：_ 實(shí) 驗(yàn) 時(shí) 間：第三周 一、 實(shí)驗(yàn)?zāi)康?掌握求系統(tǒng)響應(yīng)的方法2掌握時(shí)域離散系統(tǒng)的時(shí)域特性3分析、觀察及檢驗(yàn)系統(tǒng)的穩(wěn)定性二、 實(shí)驗(yàn)原理在時(shí)域中，本實(shí)驗(yàn)僅在時(shí)域求解。在計(jì)算機(jī)上用遞推法求差分方程的解，最簡單的方法是采用MATLAB語言的工具箱函數(shù)filter函數(shù)。也可以用MATLAB語言的工具箱函數(shù)conv函數(shù)計(jì)算輸入信號和系統(tǒng)的單位脈沖響應(yīng)的線性卷積，求出系統(tǒng)的響。三、 實(shí)驗(yàn)步驟習(xí)題一第14題2：1. 新建.m文件ex114-2;2. 程

2、序代碼：xn=ones(1,12) 2 1 1 0.5 1 1 1 1.5 ones(1,12);hn=0.2 0.2 0.2 0.2 0.2;yn=conv(hn,xn);n=0:length(yn)-1;subplot(2,2,1);stem(n,yn,.)title(a)y(n)波形);xlabel(n);ylabel(y(n)3. 運(yùn)行結(jié)果：習(xí)題一第15題：1. 新建.m文件ex115;2. 程序代碼：xn=1,2,3,4,2,1,zeros(1,10);B=1,0,2;A=1,0.5;yn=filter(B,A,xn)n=0:length(yn)-1;subplot(3,2,1);s

3、tem(n,yn,.);axis(0,15,-2,8)title(系統(tǒng)的零狀態(tài)響應(yīng));xlabel(n);ylabel(y(n)3. 運(yùn)行程序：習(xí)題一第16題：1. 新建.m文件ex116;2. 程序代碼：B1=1;A1=1,-0.6,0.08;B2=2,0,-1;A2=1,-0.7,0.1;xn=1,zeros(1,30);xi=filtic(B1,A1,0);hn1=filter(B1,A1,xn,xi);n=0:length(hn1)-1;subplot(3,2,1);stem(n,hn1,.)title(a)系統(tǒng)1的系統(tǒng)單位脈沖響應(yīng));xlabel(n);ylabel(h(n)xn=o

4、nes(1,30);sn1=filter(B1,A1,xn,xi);n=0:length(sn1)-1;subplot(3,2,2);stem(n,sn1,.)title(b)系統(tǒng)1的系統(tǒng)單位階躍響應(yīng));axis(0,30,0,3);xlabel(n);ylabel(s(n)xn=1,zeros(1,30);xi=filtic(B2,A2,0);hn2=filter(B2,A2,xn,xi);n=0:length(hn2)-1;subplot(3,2,5);stem(n,hn2,.)title(c)系統(tǒng)2的系統(tǒng)單位脈沖響應(yīng));axis(0,30,-1,2);xlabel(n);ylabel(h

5、(n)xn=ones(1,30);sn2=filter(B2,A2,xn,xi);n=0:length(sn2)-1;subplot(3,2,6);stem(n,sn2,.)title(d)系統(tǒng)2的系統(tǒng)單位階躍響應(yīng));xlabel(n);ylabel(s(n)3. 運(yùn)行結(jié)果：實(shí)驗(yàn)10.11. 新建.m文件ex10.1；2. 程序代碼：B=0.05,0.05;A=1,-0.9;xn=1,zeros(1,57);xi=filtic(B,A,0);hn=filter(B,A,xn,xi);n=0:length(hn)-1;subplot(2,2,1);stem(n,hn,.)title(a)系統(tǒng)單位

6、脈沖響應(yīng));xlabel(n);ylabel(h(n)x1n=1 1 1 1 1 1 1 1 zeros(1,50);x2n=ones(1,128);y1n=filter(B,A,x1n);subplot(2,2,2);stem(n,y1n,.)title(b)系統(tǒng)對R8(n)的響應(yīng)y1(n);axis(0,60,0,0.7);xlabel(n);ylabel(y1(n)y2n=filter(B,A,x2n);n=0:length(y2n)-1;subplot(2,2,4);stem(n,y2n,.)title(c)系統(tǒng)對u(n)的響應(yīng)y2(n);axis(0,130,0,1.2);xlabe

7、l(n);ylabel(y2(n)x1n=1 1 1 1 1 1 1 1;h1n=ones(1,10) zeros(1,10);h2n=1 2.5 2.5 1 zeros(1,10);y21n=conv(h1n,x1n);y22n=conv(h2n,x1n);n=0:length(h1n)-1;subplot(2,2,1);stem(n,h1n,.);axis(0,20,0,1.2);xlabel(n);ylabel(h1(n)title(d)系統(tǒng)單位脈沖響應(yīng)h1(n)n=0:length(y21n)-1;subplot(2,2,2);stem(n,y21n,.);xlabel(n);ylab

8、el(y21(n)title(e)h1(n)與R8(n)的卷積y21(n)n=0:length(h2n)-1;subplot(2,2,3);stem(n,h2n,.);xlabel(n);ylabel(h2(n)title(f)系統(tǒng)單位脈沖響應(yīng)h2(n)n=0:length(y22n)-1;subplot(2,2,4);stem(n,y22n,.);xlabel(n);ylabel(y22(n)title(g)h2(n)與R8(n)的卷積y22(n)A=1,-0.9;B=0.05,0.05;un=ones(1,256);n=0:255;xsin=sin(0.014*n)+sin(0.4*n);A=1,-1.8237,0.9801;B=1/100.49,0,-1/100.49;y31n=filter(B,A,un);y32n=filter(B,A,xsin);subplot(2,1,1);stem(n,y31n,.);title(h)諧振器對u(n)的響應(yīng)y31(n);axis(0,260,-0.06,0.06);xlabel(n);ylabel(y31(n)subplot(2,1,2);stem(n,y32n,.);title(i)諧振器對正弦信號的響應(yīng)y32(n);axis(0,260,-1,1);xlabel(n);ylabel(y31(n)3. 運(yùn)行結(jié)果：四、 實(shí)驗(yàn)心得及體