《信號與系統(tǒng)實(shí)驗(yàn)》課件第17章

一、實(shí)驗(yàn)?zāi)康亩?shí)驗(yàn)原理三、涉及的MATLAB函數(shù)四、實(shí)驗(yàn)內(nèi)容與方法五、實(shí)驗(yàn)要求六、思考題

一、實(shí)驗(yàn)?zāi)康?/p>

(1)掌握離散時間信號Z變換和逆Z變換的實(shí)現(xiàn)方法及編程思想;

(2)掌握系統(tǒng)頻率響應(yīng)函數(shù)幅頻特性相頻特性和系統(tǒng)函數(shù)的零極點(diǎn)圖的繪制方法;

(3)了解函數(shù)ztrans、iztrans、zplane、dimpulse、dstep和freqz的調(diào)用格式及作用;

(4)了解利用零極點(diǎn)圖判斷系統(tǒng)穩(wěn)定性的原理。二、實(shí)驗(yàn)原理

離散系統(tǒng)的分析方法可分為時域解法和變換域解法兩大類。其中離散系統(tǒng)變換域解法只有一種，即Z變換域解法。Z變換域沒有物理性質(zhì)，它只是一種數(shù)學(xué)手段，之所以在離散

系統(tǒng)的分析中引入Z變換的概念，就是要像在連續(xù)系統(tǒng)分析時引入拉氏變換一樣，簡化分析方法和過程，為系統(tǒng)的分析研究提供一條新的途徑。

Z域分析法是把復(fù)指數(shù)信號ejΩk擴(kuò)展為復(fù)指數(shù)信號zk或z=rejΩ，并以zk為基本信號，把輸入信號分解為基本信號與zk之和，則響應(yīng)為基本信號zk的響應(yīng)之和。這種方法的數(shù)學(xué)描述為Z變換及其逆變換。三、涉及的MATLAB函數(shù)

1.變換函數(shù)ztrans

功能：ztrans可以實(shí)現(xiàn)信號f(k)的(單邊)Z變換。

調(diào)用格式：

F=ztrans(f)：實(shí)現(xiàn)函數(shù)f(n)的Z變換，默認(rèn)返回函數(shù)F是關(guān)于z的函數(shù)。

F=ztrans(f，w)：實(shí)現(xiàn)函數(shù)f(n)的Z變換，返回函數(shù)F是關(guān)于w的函數(shù)。

F=ztrans(f，k，w)：實(shí)現(xiàn)函數(shù)f(k)的Z變換，返回函數(shù)F是關(guān)于w的函數(shù)。

2.單邊逆Z變換函數(shù)iztrans

功能：iztrans可以實(shí)現(xiàn)信號F(z)的逆Z變換。

調(diào)用格式：

f=iztrans(F)：實(shí)現(xiàn)函數(shù)F(z)的逆Z變換，默認(rèn)返回函數(shù)f是關(guān)于n的函數(shù)。

f=iztrans(F，k)：實(shí)現(xiàn)函數(shù)F(z)的逆Z變換，返回函數(shù)f是關(guān)于k的函數(shù)。

f=iztrans(F，w，k)：實(shí)現(xiàn)函數(shù)F(w)的逆Z變換，返回函數(shù)f是關(guān)于k的函數(shù)。

3.離散系統(tǒng)頻率響應(yīng)函數(shù)freqz

調(diào)用格式：

［H，w］=freqz(B，A，N)：其中B、A分別是該離散系統(tǒng)系統(tǒng)函數(shù)的分子、分母多項(xiàng)式的系數(shù)向量，N為正整數(shù)，返回向量H則包含了離散系統(tǒng)頻率響應(yīng)H(ejθ)在0~π范圍內(nèi)N個頻率等分點(diǎn)的值，向量θ為0～π范圍內(nèi)的N個頻率等分點(diǎn)。系統(tǒng)默認(rèn)N=512。

［H，w］=freqz(B，A，N，′whole′)：計算離散系統(tǒng)在0~2π范圍內(nèi)N個頻率等分點(diǎn)的頻率響應(yīng)H(ejθ)的值。

在調(diào)用完freqz函數(shù)之后，可以利用函數(shù)abs和angle以及plot命令，繪制出該系統(tǒng)的幅頻特性和相頻特性曲線(事實(shí)上不帶輸出向量的freqz函數(shù)將自動繪制幅頻和相頻曲線)。

4.零極點(diǎn)繪圖函數(shù)zplane

調(diào)用格式：

zplane(Z，P)以單位圓為參考圓繪制Z為零點(diǎn)列向量，P為極點(diǎn)列向量的零極點(diǎn)圖，若有重復(fù)點(diǎn)，在重復(fù)點(diǎn)右上角以數(shù)字標(biāo)出重數(shù)。

zplane(B，A)B，A分別是傳遞函數(shù)H(Z)按Z-1的升冪排列的分子分母系數(shù)行向量，注意當(dāng)B，A同為標(biāo)量時，如B為零點(diǎn)，則A為極點(diǎn)。

5.單位脈沖響應(yīng)繪圖函數(shù)dimpulse

調(diào)用格式：

dimpulse(B，A)繪制傳遞函數(shù)H(Z)的單位脈沖響應(yīng)圖，其中B，A分別是傳遞函數(shù)H(Z)按Z-1的升冪排列的分子分母系數(shù)行向量。

dimpulse(B，A，N)功能同上，其中N為指定的單位脈沖響應(yīng)序列的點(diǎn)數(shù)。

6.單位階躍響應(yīng)繪圖函數(shù)dstep

調(diào)用格式：

dstep(B，A)繪制傳遞函數(shù)H(Z)的單位脈沖響應(yīng)圖，其中B，A分別是傳遞函數(shù)H(Z)按Z-1的升冪排列的分子分母系數(shù)行向量。

dstep(B，A，N)功能同上，其中N為指定的單位階躍響應(yīng)序列的點(diǎn)數(shù)。

7.數(shù)字濾波單位脈沖響應(yīng)函數(shù)impz

調(diào)用格式：

［h，t］=impz(B，A)：B，A分別是傳遞函數(shù)H(Z)按Z-1的升冪排列的分子分母系數(shù)行向量。h為單位脈沖響應(yīng)的樣值，t為采樣序列。

［h，t］=impz(B，A，N)功能同上，其中N為標(biāo)量時指定的單位階躍響應(yīng)序列的點(diǎn)數(shù)，N為矢量時，t=N，為采樣序列。

8.極點(diǎn)留數(shù)分解函數(shù)residuez

調(diào)用格式：

［r，p，k］=residuez(B，A)：B，A分別是傳遞函數(shù)H(Z)按Z-1的升冪排列的分子分母系數(shù)行向量。r為極點(diǎn)對應(yīng)系數(shù)，p為極點(diǎn)，k為有限項(xiàng)對應(yīng)系數(shù)。四、實(shí)驗(yàn)內(nèi)容與方法

1.驗(yàn)證性實(shí)驗(yàn)

1)Z變換

確定信號f1(n)=3nε(n)，f2(n)=cos(2n)ε(n)的Z變換。MATLAB程序：

%確定信號的Z變換

symsnz%聲明符號變量

f1=3^n；

f1_z=ztrans(f1)；

f2=cos(2*n)；

f2_z=ztrans(f2)；運(yùn)行后在命令窗口顯示：

f1=

3^n

f1_z=

1/3*z/(1/3*z-1)

f2=

cos(2*n)

f2_z=

(z+1-2*cos(1)^2)*z/(1+2*z+z^2-4*z*cos(1)^2)

2)Z反變換

已知離散LTI系統(tǒng)的激勵函數(shù)為f(k)=(-1)kε(k)，單位序列響應(yīng)h(k)=

ε(k)，采用變換域分析法確定系統(tǒng)的零狀態(tài)響應(yīng)yf(k)。MATLAB程序：

%由Z反變換求系統(tǒng)零狀態(tài)響應(yīng)

symskz

f=(-1)^k;

f_z=ztrans(f);

h=1/3*(-1)^k+2/3*3^k;

h_z=ztrans(h);

yf_z=f_z*h_z;

yf=iztrans(yf_z)運(yùn)行后在命令窗口顯示：

yf=

1/2*(-1)^n+1/3*(-1)^n*n+1/2*3^n

計算，|z|＞5的反變換。

MATLAB程序：

%由部分分式展開求Z反變換num=［01］;

den=poly(［-5，1，1］);

［r，p，k］=residuez(num，den)

運(yùn)行后在命令窗口顯示：

r=

-0.1389

-0.0278-0.0000i

0.1667+0.0000ip=

-5.0000

1.0000+0.0000i

1.0000-0.0000i

k=

［］

所以反變換結(jié)果為

［-0.1389·(-5)k-0.0278+0.1667·(k+1)］·u(k)

3)離散頻率響應(yīng)函數(shù)

一個離散LTI系統(tǒng)，差分方程為y(k)-0.81y(k-2)=f(k)-f(k-2)，試確定：

(1)系統(tǒng)函數(shù)H(z);

(2)單位序列響應(yīng)h(k)的數(shù)學(xué)表達(dá)式，并畫出波形;

(3)單位階躍響應(yīng)的波形g(k);

(4)繪出頻率響應(yīng)函數(shù)H(ejθ)的幅頻和相頻特性曲線。MATLAB程序：

%(1)求系統(tǒng)函數(shù)H(z)

num=［1，0，-1］;

den=［10-0.81］;

printsys(fliplr(num)，fliplr(den)，′1/z′)

%(2)單位序列響應(yīng)h(k)的數(shù)學(xué)表達(dá)式，并畫出波形subplot(221);

dimpulse(num，den，40);ylabel(′脈沖響應(yīng)′);

%(3)單位階躍響應(yīng)的波形

subplot(222);

dstep(num，den，40);

ylabel(′階躍響應(yīng)′);

%(4)繪出頻率響應(yīng)函數(shù)的幅頻和相頻特性曲線

［h，w］=freqz(num，den，1000，′whole′);

subplot(223);

plot(w/pi，abs(h))；ylabel(′幅頻′);

xlabel(′＼omega/＼pi′);

subplot(224);

plot(w/pi，angle(h))；

ylabel(′相頻′);

xlabel(′＼omega/＼pi′);

運(yùn)行后在命令窗口顯示：

num/den=

-11/z^2+1

-----------------------

-0.811/z^2+1

系統(tǒng)的響應(yīng)與頻率響應(yīng)函數(shù)如圖17.1所示。圖17.1系統(tǒng)的響應(yīng)與頻率響應(yīng)函數(shù)

4)MATLAB繪制離散系統(tǒng)極點(diǎn)圖

采用MATLAB語言編程，繪制離散LTI系統(tǒng)函數(shù)的零極點(diǎn)圖，并從零極點(diǎn)圖判斷系統(tǒng)的穩(wěn)定性。

已知離散系統(tǒng)的H(z)，求零極點(diǎn)圖，并求解h(k)和H(e^jw)。MATLAB程序：

b=［121］;

a=［1-0.5-0.0050.3］;

subplot(3，1，1);

zplane(b，a);

num=［0121］;

den=［1-0.5-0.0050.3］;

h=impz(num，den);

subplot(3，1，2);stem(h);

%xlablel(′k′);

%ylablel(′h(k)′);

［H，w］=freqz(num，den);

subplot(3，1，3);

plot(w/pi，abs(H));

%xlablel(′/omega′);

%ylablel(′abs(H)′);

系統(tǒng)的響應(yīng)與零極點(diǎn)分布如圖17.2所示。圖17.2系統(tǒng)的響應(yīng)與零極點(diǎn)分布5)直線型系統(tǒng)函數(shù)的Z域分布

直線型系統(tǒng)函數(shù)為

H(z)=

試求其零點(diǎn)和極點(diǎn)，并將它轉(zhuǎn)化為二階節(jié)形式。MATLAB程序：

num=［1-0.1-0.3-0.3-0.2］;

den=［10.10.20.20.5］;

［z，p，k］=tf2zp(num，den);m=abs(p);disp(′零點(diǎn)′);disp(z);

disp(′極點(diǎn)′);disp(p);disp(′增益系數(shù)′);disp(k);

sos=zp2sos(z，p，k);disp(′二階節(jié)′);disp(real(sos));zplane(num，den)計算求得零、極點(diǎn)增益系數(shù)和二階節(jié)的系數(shù)分別為：

零點(diǎn)

0.9615-0.5730-0.1443+0.5850i-0.1443-0.5850i

極點(diǎn)

0.5276+0.6997i0.5276-0.6997i-0.5776+0.5635i-0.5776-0.5635i

增益系數(shù)

1二階節(jié)

0.1892-0.0735-0.10431.00001.15520.6511

5.28461.52471.91851.0000-1.05520.7679

系統(tǒng)的零極點(diǎn)分布如圖17.3所示。圖17.3系統(tǒng)的零極點(diǎn)分布系統(tǒng)函數(shù)的二階節(jié)形式為

H(z)=

極點(diǎn)圖見圖17.3。

2.程序設(shè)計實(shí)驗(yàn)

(1)試分別繪制下列系統(tǒng)的零極點(diǎn)圖，并判斷系統(tǒng)的穩(wěn)定性。

(a)H(z)=

(b)H(z)=

(c)H(z)=

(d)H(z)=

(2)