信號與系統(tǒng)上機實驗報告.doc

姓 名： 學 號： 學 院：Q1:0=n=31 x1(n)=sin(pi*n/4)*cos(pi*n/4); x2(n)=cos(pi*n/4)*cos(pi*n/4) x3(n)=sin(pi*n/4)*cos(pi*n/8) 分別畫出圖形,求出其周期。代碼：n=0:31;x1=sin(pi*n/4).*cos(pi*n/4); x2=cos(pi*n/4).*cos(pi*n/4); x3=sin(pi*n/4).*cos(pi*n/8);stem(n,x1)stem(n,x2)stem(n,x3)結果：由圖形可知周期T1=4由圖形可知周期T2=4由圖形可知周期T3=16Q2:當0=n=5時，h(n)=n; 其他h(n)=0; x(n)=h(n);求y(n)=x(n)*h(n); 用stem函數(shù)畫出y(n). 代碼：n=0:5;h=n;x=h;y=conv(x,h);stem(y)結果：Q3：（a）.定義用向量a1和b1描述差分方程y(n)-0.8y(n-1)=2x(n)-x(n-2)表征的因果LTI系統(tǒng)，(b).用在（a）中的系數(shù)向量，利用freqz定義H1是在0和pi之間4個等份頻率上頻率響應的值，omega1是這些頻率值。(c).用在（a）中的系數(shù)向量，利用freqz定義H2是在0和2*pi之間4個等份頻率上頻率響應的值，omega2是這些頻率值。代碼： n=4;a1=5,0,-4b1=10,0,-5H1,W1 = freqz(b1,a1,n)H2,W2 = freqz(b1,a1,n,whole)結果：a1 = 5 0 -4b1 = 10 0 -5H1 = 5.0000 1.7073 - 0.3659i 1.6667 1.7073 + 0.3659iW1 = 0 0.7854 1.5708 2.3562H2 = 5.0000 1.6667 5.0000 1.6667W2 = 0 1.5708 3.1416 4.7124Q4: X1(n)=u(n)-u(n-8); 其周期N1=8, X2(n)=u(n)-u(n-8); 其周期N2=16, X3(n)=u(n)-u(n-8); 其周期N3=32, (1) 畫出這些周期信號在0=n=1&i=17&i=33&i=49&i=1&i=33&i0.9 (1). 畫出H(z)的零極點圖。(2). 定義omega=0:511*pi/256和unitcirc=exp(j*omega)得到單位圓上512個等分點。定義ps為極點的列向量，zs為零點的列向量，我們定義polevector=ones(2,1)*unitcirc-ps*ones(1,512) zerovector=ones(2,1)*unitcirc-zs*ones(1,512) 令：polelength和poleangle分別為polevector的輻值和相位角，zerolength和zeroangle分別為zerovector的輻值和相位角。畫出polelength和zerolength對于omega的圖，畫出poleangle和zeroangle對于omega的圖。（3）用polelength和zerolength求出abs(H(exp(jw) 用poleangle和zeroangle求出angle(H(exp(jw). 并分別畫出其圖形。（1） 代碼：a=1 0 0;b=1 -0.9 0.81;zplane(b,a);結果：(2) （3）代碼：w=0:511*pi/256;unitcirc=exp(1j*w);b=100,0,0;a=100,-9,81;zplane(b,a);p=roots(a);z=roots(b);polevector=ones(2,1)*unitcirc-p*ones(1,512);zerovector=ones(2,1)*unitcirc-z*ones(1,512);polelength=abs(polevector);poleangle=angle(polevector);zerolength=abs(zerovector);zeroangle=angle(zerovector);plot(w,polelength);figure;plot(w,poleangl