信號(hào)與系統(tǒng)實(shí)驗(yàn)10

1、西北工業(yè)大學(xué)信號(hào)與系統(tǒng)實(shí)驗(yàn)報(bào)告第十章 Z變換實(shí)驗(yàn)報(bào)告學(xué) 院： 軟件與微電子學(xué)院 學(xué)號(hào)： 姓名： 專 業(yè)： 軟件工程 實(shí)驗(yàn)時(shí)間： 2016.11 實(shí)驗(yàn)地點(diǎn)： 毅字樓311 指導(dǎo)教師： 陳勇 西北工業(yè)大學(xué)2016 年 11 月一、 實(shí)驗(yàn)?zāi)康亩?實(shí)驗(yàn)要求 三、 實(shí)驗(yàn)設(shè)備（環(huán)境）四、 實(shí)驗(yàn)內(nèi)容與步驟五、 實(shí)驗(yàn)結(jié)果§10.1 作離散時(shí)間零極點(diǎn)圖1用plotdpz畫出(10.3)或的零極點(diǎn)圖。代碼：%畫零極點(diǎn)圖,num表示分子多項(xiàng)式，den表示分母多項(xiàng)式function splane(num,den)p = roots(den) % 求極點(diǎn)q = roots(num) % 求零點(diǎn)p = p&#

2、39;q = q'x = max(abs(p q); % 確定實(shí)軸范圍 x = x+1;y = x; % 確定虛軸范圍 plot(-x x,0 0,':');hold on; % 畫實(shí)軸plot(0 0,-y y,':');hold on; % 畫虛軸plot(real(p),imag(p),'*');hold on; % 畫極點(diǎn)plot(real(q),imag(q),'o');hold on; % 畫零點(diǎn)title('zero-pole plot');xlabel('實(shí)部');ylabe

3、l('虛部')axis(-x x -y y); % 確定顯示范圍主程序：clear;clc;b=1 -1 0;a=1 3 2;splane(b,a)運(yùn)行結(jié)果：2用plotdpz畫出滿足差分方程的濾波器的零極點(diǎn)圖。代碼：clear;clc;b=1;a=1 1 0.5;splane(b,a)運(yùn)行結(jié)果：3用plotdpz畫出滿足差分方程 所描述濾波器的零極點(diǎn)圖。代碼：clear;clc;b=1 0.5;a=1 -1.25 0.75 -0.125;splane(b,a)運(yùn)行結(jié)果：§10.2離散時(shí)間頻率響應(yīng)的幾何解釋1定義b1和a1包含按filter要求格式的分子和分母多項(xiàng)式的

4、系數(shù)，按10.1節(jié)所介紹的方法求出并畫出的零極點(diǎn)圖。解：clc;clear;close all;b1=1 0 0;a1=1 -0.9 0.81;plotdpz(b1,a1);grid on;title('H1(z的零極點(diǎn)分布');2定義omega=0:511*pi/256和unitcirc=exp(j*omega)得到在單位圓上512等分點(diǎn)，在這些點(diǎn)上將要對(duì)頻率響應(yīng)求值。定義polevectors1是一個(gè)2×512的矩陣，其中每一行包含一些復(fù)數(shù)，這些復(fù)數(shù)是由unitcirc的相應(yīng)列減去一個(gè)極點(diǎn)位置得到的。如果ps1是一個(gè)列向量，它包含了極點(diǎn)的位置，那么就能用下列命令來

5、完成>> polevectors1=ones(2,1)*unitcirc-ps1*ones(1,512);在MATLAB中總是愿意用矩陣和向量運(yùn)算，因?yàn)镸ATLAB幾乎總是計(jì)算矩陣運(yùn)算比用for循環(huán)構(gòu)成的運(yùn)算要快。用abs和atan2定義polelength1和poleangle1作為polevectors1中每一元素的幅值和相位。解：ps1=roots(a1);omega=0:511*pi/256;unitcirc=exp(j*omega);polevectors1=zeros(2,512);polevectors1=ones(2,1)*unitcirc-ps1*ones(1,5

6、12);polelength1=abs(polevectors1);poleangle1=atan2(imag(polevectors1),real(polevectors1);3類似于polevecters1定義zerovectors1，使得它是2×512的矩陣，其中包含從零點(diǎn)位置到unitcirc元素的向量。定義zerolength1和zeroangle1分別是這些向量的幅值和相位。解：zerovectors1=ones(2,512);zs1=roots(b1);zerovectors1=ones(2,1)*unitcirc-zs1*ones(1,512);zerolength1

7、=abs(zerovectors1);zeroangle1=atan2(imag(zerovectors1),real(zerovectors1);4畫出polelength1和zerolength1對(duì)omega的圖。根據(jù)這些圖，你預(yù)計(jì)在哪里有最大值和最小值？解：figure;subplot(1,2,1);plot(omega,polelength1);grid on;subplot(1,2,2);plot(omega,zerolength1);grid on;5用polelength1和zerolength1計(jì)算，并將這個(gè)結(jié)果存入geonH1mag中。用poleangle1和zeroangle1計(jì)算，并將這個(gè)結(jié)果存入gemoH1phase中。在定義geonH1mag和gemoH1phase時(shí)，你或許會(huì)發(fā)現(xiàn)prod和sum是很有用的。畫出從幾何上到處的幅值和相位，并將它和經(jīng)由下面計(jì)算所得結(jié)果比較：>> H1=freqz(b1,a1,512,'whole');你在4中估計(jì)的峰值頻率正確嗎？解：geon