計算機控制系統(tǒng)的描述與分析

動力工程計算機控制系統(tǒng)

能源與動力工程學院

本科專業(yè)學優(yōu)課程

48學時，3學分

7

第四章計算機控制系統(tǒng)的

描述與分析

主講人:黃勇理

7

第四章計算機控制系統(tǒng)的描述與分析

?4.1計算機控制系統(tǒng)的描述與分析

?4.1.1計算機控制系統(tǒng)的信號變換

?4.1.2線性離散系統(tǒng)的數學描述

?4.1.3微分方程向差分方程的近似轉換

?4.1.4差分方程的解法

?4.2Z變換的基本理論

?4.2.1采樣信號的數學描述與采樣定理

?4.2.2變換的定義、定理和性質

?4.2.3求Z變換與Z反變換

計算機控制系統(tǒng)_HYL

第四章計算機控制系統(tǒng)的描述與分析

?4.3脈沖傳諦函數及其捽制系統(tǒng)描述

?4.3.1用z變換法解差分方程

?4.3.2建立z傳遞函數

?4.3.3離散控制系統(tǒng)的結構分析

?4.4線性離散系統(tǒng)性能分析

?4.4.1線性離散系統(tǒng)穩(wěn)定性分析

?4.4.2線件離散系統(tǒng)過渡過程分析

?4.4.3線件離散系統(tǒng)穩(wěn)態(tài)性能分析

?4.4.4線性離散系統(tǒng)動態(tài)件能與極點分布

計算機控制系統(tǒng)_HYL

4.1計算機控制系統(tǒng)的描述與分析

7

4.1.1計算機控制系統(tǒng)的信號變換

?計算機控制系統(tǒng)的信號變換和傳遞流程

、??????一《?.??????.?，、?、?-???、、

??二、，?a.i??一(t，.?-????-??????

―'■>a?■■-，、?，?/、*iz，t-?■?,■?.-.■..?t,-■；b??.、--*aynr<(------,?.1,-.

i.ih-?/riiHili-fve*/TTif、*ii」i=-x=r'37."-n'I，1IM?

”.?*r—1/>*■??*?w(■?，??，y?、?-，.a?.J?▼????■a?i—a?▼?-A-.,—/、.▼?.?v、■(▼“

4.1.1計算機控制系統(tǒng)的信號變換

?信號的四種類型

1.連續(xù)信號：時間連續(xù)、幅值連續(xù)；

2.離散模擬信號：在一系列離散的時刻上，

信號幅值連續(xù)；

3.數字信號：離散時刻上，幅值不連續(xù)，且

數值經過了整量化處理；

4,整量化連續(xù)信號：也稱階梯信號，時間上

連續(xù)，幅值是階梯狀的。

計算機控制系統(tǒng)_HYL

4.1.1計算機控制系統(tǒng)的信號變換

?信號的四獨類型____________________

tt

數字信號整量化連續(xù)信號

4.1.1計算機控制系統(tǒng)的信號變換

(1)采樣開關

?利用采樣開關周期性的瞬時開啟，把連續(xù)信號y(t)轉

變成離散的信號序列y*(t),這個過程稱為采樣。

?通常采樣開關以等時間間隔T開啟，稱T為采樣時間。

?若系統(tǒng)中的多個采樣開關以等周期同時開閉，則稱為

同步采樣；

?若以等周期，但不同時采樣，稱為非同步采樣；

?若多個開關以不同周期采樣，則稱為多速采樣。

計算機控制系統(tǒng)_HYL

4.1.1計算機控制系統(tǒng)的信號變換

?實際采樣過程需要持續(xù)一段時間。

若：T?T,可以認為：

y(kT)=y(kT+AZ)0<M<T

若在采樣周期內有N個巡回監(jiān)測點，應保證：

T>Nx

計算機控金YLIIII""、1

0T2T3T4T5T6Tt

4.1.1計算機控制系統(tǒng)的信號變換

⑵模數轉換(A/D)

?數字計算機只能處理二進制的數值。故采樣得到

的離散模擬信號幅值需要變成二進制數的最小單

位q的整數倍，這個過程就稱為整量化或稱量化。

?若信號y*(t)的幅值并非q的整數倍，則量化過

程存在誤差。從信號傳遞的角度看，量化誤差是

一種噪聲也是信號的損耗。

計算機控制系統(tǒng)_HYL

4.1.1計算機控制系統(tǒng)的信號變換

(3)數模轉換(D/A)與保持器

?控制對象的測量信號轉換成數字信號y(kT)后，與給

定數字信號r(kT)比較，得到偏差數字信號e(kT),輸

入控制器D(z),經運算處理得到輸出數字信號u(kT)。

?對于時間離散的數字信號，模擬執(zhí)行器(機構)通常

是接受不了也執(zhí)行不了的,所以必須將時間離散的數

字信號重構成時間連續(xù)的信號。

計算機控制系統(tǒng)_HYL

4.1.1計算機控制系統(tǒng)的信號變換

⑶數模轉換(D/A)與保持器

?控制量的數字信號u(kT)先經數模轉換，變成離散模

擬信號U*(t),再通過零階保持器將kT時刻的信號值

保持到時刻(k+1)T,得到能被執(zhí)行機構接受的控制量

u(t)o即：

u(kT)=u(kT+X)0<^t<T

?控制量u(t)是時間連續(xù)的幅值整量化連續(xù)信號。通常

D/A具有零階保持器的功能。

計算機控制系統(tǒng)_HYL

4.1.1計算機控制系統(tǒng)的信號變換

零階保持器：

時域特性：

h0(t)=u(t)-u(t-T)

傳遞函數：

11_Ts

-Ts11—e

(S)=—e-二-------------------

0

頻率特性：HSss

幅頻特性:

相頻特性:/Ho(jco)=-coT12

計算機控制系統(tǒng)_HYL

計算機控制系統(tǒng)—

422線性離散系統(tǒng)的數學描述

?計算機處理的是離散數字信號，計算機控制系統(tǒng)是線性

離散系統(tǒng)或近似為線性離散系統(tǒng)。

?線性連續(xù)系統(tǒng)常微分方程：

心0立⑺

"o+…Q+Q(/)

dtnId嚴a〃八,

心?智包攻)

0dtm1dtmA"1dtm

計算機控制系統(tǒng),

4.1.2線性離散系統(tǒng)的數學描述

?數字調節(jié)器可表示成線性離散系統(tǒng)，與線性連續(xù)系統(tǒng)

類似，輸入與輸出之間的關系可用常系數微分方程描

述，即：

y(kT)+a}y(kT-T)+a2y(kT-2。H----卜any(kT-nT)

=bd(kT)+bHkT-T)+b2r(kT-20+…+bmr{kT-mT)

,線性離散系統(tǒng)的分析方法也有古典法、變換法和離散

狀態(tài)空間法。線性之義就是符合疊加原理，因此對線

性離散系統(tǒng)來說，表征其特性的差分方程應滿足疊加

原理。

16

4.1.3微分方程向差分方程的近似轉換

下面以大家所熟知的慣性環(huán)節(jié)為例說明：

K

G(S)=

TS+1

我們可以由差分概念,K微分方程推導出近似的

差分方程。上式作如下變化:

G⑸=^~;=瞿(T.s+l)Y(s)=KR(s)

TQS+1R(S)

作拉普拉斯反變換，得慣性環(huán)節(jié)的微分方程:

T.y(t)+y(t)=Kr(t)

計算機控制系統(tǒng)_HYL

計算機控制系統(tǒng)—

4.1.3微分方程向差分方程的近似轉換

?用y(kT)、=(取)分別代表丫仕)、r(t)的采樣值。用差

分代替微分

=叱J（kT）-y（kT-T）

dt

?將y(kT)、r(kT)及上式代入微分方程得：

yly(kT)-y(kT-T)]+y(kT)^Kr(kT)

+)y(kT)-y(kT-T)^Kjr(kT)

?此即所描述的連續(xù)系統(tǒng)離散化后的差分方程

18

4.1.4差分方程的解法

?單輸入/單輸出(SISO)線性定常系統(tǒng)，對應的差分方

程常用的求解方法有：迭代(遞推)法、古典法和z變

換法二種。

?(1)迭代法：如果已知差分方程和輸入序列，并且給出了輸出

序列的初始值，則可直接利用迭代關系逐步算出所需要的輸出

序列。詳見教材P45頁

?(2)古典法：與微分方程類似，求解差分方程也有古典法，差

分方程的全解也包括兩部分，即對應齊次方程的通解(瞬態(tài)解)

和非齊次方程的一個特解(穩(wěn)態(tài)解)。具體解法見教材P46頁。

?(3)Z變換法：對線性連續(xù)系統(tǒng)，拉氏變換使求解微分方程的微

積分運算簡化成代數運算。同樣，對線性離散系統(tǒng)，下面將要

介紹的Z變換法也可大大簡化差分方程的求解。***

計算機控制系統(tǒng)_HYL

4.2Z變換的基本理論

7

計算機控制系統(tǒng)—

4.2.1采樣信號的數學描述與采樣定理

?(1)、采樣信號的表達式

?把連續(xù)信號變換為脈沖序列的過程稱為采樣。對一個

連續(xù)信號采樣，就是用一系列離散值代替連續(xù)信號，

這一系列的離散值等于采樣時刻連續(xù)信號的幅值，稱

為采樣值，

?亦即：fy(kT)

?如果給采樣器輸入一個連續(xù)函數歹仞，得到的輸出將

是一脈沖序列y^(t)，它可以表示為：

*

y⑴=y(O)^(t)+—+…

+y(5)3(，一初)+…

21

計算機控制系統(tǒng)—

4.2.1采樣信號的數學描述與采樣定理

?上式中，5（0）或8（以7）為DiracDelta函數，也稱為單位

脈沖函數。單位脈沖函數定義為：

「oot=kT

[0t^kT

J二3（，-左T）力=1

?根據單位脈沖的襦選特性，脈沖序列可表示為：

0000

*

y⑺=打)3(t——)=義譜”kT)

k=0k=0

(k<0,y(kT)=0)

22

計算機控制系統(tǒng)

4.2.1采樣信號的數學描述與采樣定理

?（2）、采樣定理

?可以證明，對一個具有有限頻譜（。＜3也辦）的連續(xù)信

號y⑺采樣，當且僅當采樣頻率3s（或采樣周期T）

滿足關系：

CDs—2comaxorTW兀/?max

?此時，采樣信號y*⑺就能無失真地恢復為原來的連

續(xù)信號?、?，這就是著名的“香農（Shanon）采樣定

理”。

23

4.2.2Z變換的定義、定理和性質

?一、z變換的定義

?z變換可由采樣信號的拉氏變換推導出來。

?對于采樣信號：

00

*

y(0=Ev(kT)3(t—kT)

k=a

?其拉斯變換為：

00

*_*Q00________sOt7

L[y(0]=y(s)=J()Zy(kT)5?！猭T)e-dt

k=0

00

=Z>(左(t-kT)e-stdt

k=0

00

=Zy(kT)e-kTS

k=0

計算機控制系統(tǒng)_HYL

計算機控制系統(tǒng)—

4.2.2Z變換的定義、定理和性質

?上式中，含有超越函數05，為了運算和書寫方便，

我們引入一個新的變量Z：

Ts、1

s

z-e或s二一Inz

T

?將y*⑸記為Hz）,則有：

00

必（祝=丫*（5）=＞（口）產

k=0

00

丫⑶=y*=￡y(kT)z.k

k=0

4.2.2Z變換的定義、定理和性質

?由此得出以Z為自變量的函數y（z）,我們稱它為y*⑺

的Z變換，記作：

y（z）=

?Z變換的定義為：

00

y（z）=Ey（kT）z』=X0）+y3+y（2T）z-2++…

k=0

?由上式可以看出，采樣函數V⑺的z變換y（z）與采樣

點的采樣值M股）有關，或者說，級數中4的系數就

是時間序列M股）的值。

?一個函數的Z變換只有在采樣時刻才有意義。

計算機控制系統(tǒng)_HYL

4.2.2Z變換的定義、定理和性質

?二、典型函數(或序列)的z變換

?(1)數字脈沖系列6左

?若將"定義為：]

n1k-n

§,=<

0k^n

n

?則有：Z[8k.n]=z

在求Z變換時，我們作如下規(guī)定：

函數八。在t=kT時刻的采樣值人左7),就是在該時

刻的脈沖強度。而5⑺的脈沖強度是1,于是有：

Z[5(/)]=Z[5d=l

Z[d(kT-nT)]=Z[dk_n]=z-n>0

計算機控制系統(tǒng)_HYL

4.2.2Z變換的定義、定理和性質

?（2）單位階躍函數

?對于單位階躍函數，有：蛇尸1,「三o

*

u（，）=u{kT）=1k>0

00

k1—1—2

U(z)=Z\u{kT)]=Vz=1+z+z+-??

可以看出，上式為無窮遞減等比數列求和。由無窮

等比數列求和公式，得：

t/（z）=l+z-1+z~2H-———

1-zz-1

?思考，單位脈沖序列的Z變換是什么?

計算機控制系統(tǒng)_HYL

計算機控制系統(tǒng)

422Z變換的定義、定理和性質

?（2）單位階躍函數

對于單位數字階躍序列:

1k>0

〃k=q

0左<0

其Z變換結果與U⑴的Z變換結果是一樣的

U(z)=1+z-1+z-2d—=—二-=z

1-zz-1

29

計算機控制系統(tǒng)

4.2.2Z變換的定義、定理和性質

?(3)單位速度函數

f3=t^>0

或：f(kT)=kT左二0,1,2,???

由階躍函數的Z變換

001

U(z)=1+z-1+z-2H——-\^z~k

k=6~l-z

兩邊對Z求導后，再乘以T,整理后可得:

Tz-11

F(z)=Z[kT]=-~~

(1—z)

30

計算機控制系統(tǒng)

4.2.2Z變換的定義、定理和性質

?(4)指數函數

f(t)=e~att>0

或：f左=0,1,2,???

?離散化后指數函數的Z變換為:

F(z)=Z[e~akT]

0000

、、—cikT—k、、/—ctkT\—k

=>ezz)

邑k=b

1z

1-ezz-e—ctT

31

計算機控制系統(tǒng)

4.2.2Z變換的定義、定理和性質

?(5)正余弦函數

/(0=sinco^和f(0=cosco^t>0

f(kT)=sinskT和f(kT)=cosskTk二0,1,2,…

?根據歐拉公式：

f(kT)=eJMkT=cosskT+7sinco左T

001

F(Z)=Zl*kT]=￡U=匚擊廣

_1-(coscoT)z-1(sincoT)z-1

1-(2coscoT)z-1+z-211-(2coscoT)z-1+z-2

?比較上面兩式可得：

z1(sincoT)z-111-(coscoT)z-1

Z[sinco/]=------------------；-----y和Z[cosco/]=------------------；-----5

l-(2coscoT)z-1+z-2l-(2cosmn^^^

計算機控制系統(tǒng)—

4.2.2Z變換的定義、定理和性質

?三、Z變換的性質和定理

?Z變換的性質和定理與拉氏換的性質和定理非常

相似。下面介紹幾種常用的Z變換性質和定理。

?(1)線性性質

?設Z[y/D]=y(z),Z[雙U)]=X(z),且Q、6為常數，則有：

Z[ay(kT)]=aY(z)

Z\bx(kT)]=bX(z)

Z[ay(kT)+Z?x(kT)]=aY(z)+bX(z)

?這個性質說明，z變換是一種線性變換。

33

計算機控制系統(tǒng)—

422Z變換的定義、定理和性質

?(2)平移定理：

?包括滯后平移和超前平移

?滯后定理

?設Z[y/T)]=y(z),且.TV0時,y(kT)=0,則:

Z[y(kT-nT)]=z~nY(z)

?滯后定理說明：時域中滯后某采樣信號丁(左「有〃個采樣周

期的采樣信號y/T—⑺的Z變換，等于原采樣信號的雇

換Y(z)乘以7〃。實際上，7〃代表滯后環(huán)節(jié)，它使采樣信號

延遲了〃個采樣周期，或者說，推遲了時刻才開始采

樣,如后圖所示：

34

計算機控制系統(tǒng)

4.2.2Z變換的定義、定理和性質

?實際上，”代表滯后環(huán)節(jié)，它使采樣信號延遲了〃個采樣

周期，或者說，推遲了〃丁時刻才開始采樣，如圖：

（0+1）T

ZF的滯后特性

35

計算機控制系統(tǒng)—

4.2.2Z變換的定義、定理和性質

?超前定理

?實際的超前環(huán)節(jié)是不存在的，但在運算中可能出現。

?設Z[y%T)]=y(z),則:

Z[y(kT+〃)]=znY(z)-zn￡"〃尸

J=°

?z〃代表超前環(huán)節(jié)，表示輸出信號超前輸入信號〃個采樣周

期。特別:

當n=l時，有：Z[y(kT+T)]=zy(z)-z^(O)

當口=1時，有：Z[y(kT+27)]=Z2Y(Z)-z2y(0)-zy(T)

若火。)=y(r)=??C=。

貝lj：Z[y(kT+nT)]=znY(z)

36

計算機控制系統(tǒng)—

4.2.2Z變換的定義、定理和性質

?(3)初值定理：

?設Z[y/T)]=y(z),則:

y(O)=limy(kt)=lim7(z)

k=0Zf8

?例：求單位階躍序列〃(左T)的初值。

?［解］：已知單位階躍序列的Z變換為：

?故由初值定理得:

1

〃(0)=lim〃(左力=lim------；=1

k=az—>00

37

計算機控制系統(tǒng)—

4.2.2Z變換的定義、定理和性質

?(4)終值定理：

?設Z[y/T)]=y(z),則:

y(oo)=limy(kt)=lim(z-l)y(z)

左=00Zf1

?例：求單位階躍序列〃(左T)的終值。

?[解]：已知單位階躍序列的Z變換為：

Z[〃(W]=J^r

l-z

?故由終值定理得：

u(co)=hmu(kt)=hm(z-l)-^=l

38

計算機控制系統(tǒng)

4.2.2Z變換的定義、定理和性質

?（5）迭值定理：（省略）

?（6）卷積定理：（省略）

39

計算機控制系統(tǒng)—

4.2.3Z變換與Z反變換方法

?一、z變換方法

?z變換方法有三種：

?無窮級數求和法

?部分分式法

?留數計算法

?（1）無窮級數求和法

?是按照定義式直接計算的方法。將z變換定義式:

*00

y（z）=y（S）=fy(kT)z-k

s=—\nz

Tk=0

?展開成無窮級數求和式即可得。前述典型函數的Z變換

都是按定義展成無窮級數，再整理成閉合解析式的。

40

計算機控制系統(tǒng)

4.2.3Z變換與Z反變換方法

?（2）部分分式法

?部分分式法是由函數的拉氏變換求Z變換的方法。

?當已知一個函數y（s）時，首先將y（s）展成部分分式：

?w（z=l,2,n）

ait

口匕⑸的拉氏反變換都是指數函數:匕4）=Ate-

?然后利用指數函數的Z變換式可求得y⑺的夜換。

?將函數的拉氏變換Y（s）展開成部分分式：

n

y（s）=z4

i=lS-p

式中：口.式）為y（s）的一階極點（p尸一4?）；4為各分式的

待定系數。部分分式中的系數4可用初等數學中的待定系數法求

出，或用留數定理求出。

41

計算機控制系統(tǒng)

4.2.3Z變換與Z反變換方法

?由留數定理可知：

4=(s-“)y(s)|…

?對應分式的拉氏反變換為：’

廠[-^]=4/"

?利用指數函數的Z變換公立，有：

z[A-ait]=a——

iez-e'

?由此得到Y(s)對應的y(t)的Z變換：

〃7〃1

y(z)=Z4一廳=Z4]"T

i=iz-ei=i1-ez一

42

計算機控制系統(tǒng)

4.2.3Z變換與Z反變換方法

?例:已知連續(xù)函數了⑴的拉氏變換如下，求丫(Z)。

y(s)=---------------

(S+1)(5+2)

?解：部分分式展開有：

111

Y(s)=---------二――----

6+1J6+2)5+1s+2

?則可得：J1

^(Z)~i~~~——\^2T~~

1-ez1-ez

(1-e)ez

/1-T—1\zi—2T—1\

(1-ez)(\-ez)

43

計算機控制系統(tǒng)

4.2.3Z變換與Z反變換方法

崢例：已知連續(xù)函數/⑴的拉氏變換如下，求Y（z）。

1

y（s）=-.......

S（S+1）

?解：部分分式展開有：

“、1abc

Y(s)=-......=—+-7+-------

s(s+1)sss+1

?求系數a”,c:

d2d11

a=-[sy(s)]=—[----]=-------2二—1

dss=odss+1I(s+1)

5=0

1

1s2y⑸]=oc=[(s+i)y(s)ki=?=1

S+1s=0

計算機控制系統(tǒng)

4.2.3Z變換與Z反變換方法

?因此得：1111

Y(s)=1--------+

5(5+1)75+1

?三部分分別有：

11

￡、T=1⑺Z[1(O]

s1-z-1

1x

[Tz~

L-[-]=tZ[〃=-1x2

s(1-Z)

1-t7「T1_1

]=eZ[e]—"ry__zy

5+11-ez

?三部分相加得：7z-1

y（z）=—1-1+-l\2+i~~T-1

1-Z(1-Z)1=0臺

45

計算機控制系統(tǒng)

4.2.3Z變換與Z反變換方法

?（3）留數計算法

?若已知連續(xù)函數y⑺的拉氏變換y（s）及其全部極點:

“《=12…，辦則可用下列留數計算式求得丫⑶：

m7

r（z）=Z[/（/）]=ZResy（s）一浮

ILz-e人

mf14見T「

=￡\------lim——r（s-pj'Y（s）----7-

總[（%-1）!$*"Iz-eTs]

?式中：/為y（s）的極點口的階數；加為y（s）的彼此不同極點

個數。

46

計算機控制系統(tǒng)

4.2.3Z變換與Z反變換方法

二、z反變換法

?如果已知丫⑵，求對應的時間序列或數值序

歹1上因，則稱此運算為Z反變換，記作：

y(kT)=Z-l[Y(z)]或X^)=Z-1[y(z)]

?常用的Z反變換法也有三種：

?長除法

?部分分式法

留數計算法

47

計算機控制系統(tǒng)—

4.2.3Z變換與Z反變換方法

?(1)長除法

設丫(Z)為多項式分式:

如加+給小+…+圖

丫⑶二nn-\(H>m)

a^z+axzH-------ban

長除法就是用Y(z)的分子直接除以分母，長除的結果是

z的降幕級數:

00

2』

丫(2)=%+M2T+J/2z_+???+”24+???=EykZ

k=0

48

計算機控制系統(tǒng)_HYC、

4.2.3Z變換與Z反變換方法

?(1)長除法

此外，根據Z變換的定義：

00

*__小

y(z)=z)=

k=0

=y(Q)+y(T)z-x+y(2T)z~2+???+y(kT)z-k+…

?對比長除法的結果和定義式可知：

一(。)=%，y(T)=->(2T)=%,…，y(左T)=力，…

?若多項式分式中的分子關于Z的最高階為一,而分母關于z的

最高階為z3,則長除的結果關于Z的最高階顯然只能為2—2,

因止匕作Z反變換后，肯定有興0)=0,興7>0,延時兩拍。

49

計算機控制系統(tǒng)

4.2.3Z變換與Z反變換方法

?（2）部分分式法

?設丫⑵為:

7m7"2-11

bQz+給+…

y（Z）=--------------------------（H>加）

n（z-A）

?將等式兩邊除以zj：再展成部分分式：

y（z）az*2j./4

一n—乙

zn（z-口）aPt

z=l

其中：」（zeW：+L+，

[n（-A）

計算機控制系統(tǒng)

4.2.3Z變換與Z反變換方法

(2)部分分式法

則丫(Z)的z反變換為:

nA7n

y/T)=Z")］二z」E——=EA(PN

i=\Z-Pi\1

觀察并理解下式除以Z的目的……

y(z)boz"—+Az"”?H—(4-

zn—pjiz-Pi

Z=1

其中：4=(zg空二小

n(z_Pi)

_，=1」z="K￡J

計算機控制系統(tǒng)—

4.2.3Z變換與Z反變換方法

?(2)部分分式法

求y（z）==的z反變換。

例:

?'解I：y(-z?+2z?+z+l_z?+2z?+z+l

z-Z3-Z2-8Z+12-(Z-2)2(Z+3)

,AzAz24z

=4x+—二~~7+^—

z-2(z-2)2z+3

2

B^z+B}z+B2Z+B3

(Z-2)2(Z+3)

52

計算機控制系統(tǒng)

4.2.3Z變換與Z反變換方法

?(2)部分分式法

?比較丫⑵分子部分各項系數。

—4。+4+42+43=1

4=+4+34—44—2

<

B2——8Z()—6Al+44-1

B3=12Ao=1

可聯立求解得：A=—,4=—,A

o125022075

所以：y(z)」一2二+”—+〃工

1250z-220(z-2)75z+3

y(kT)=—5(Z:T)--2A+—(/;+l)2k+—(-3/

53

計算機控制系統(tǒng)

4.2.3Z變換與Z反變換方法

■(3)留數計算法

?留數計算法又稱反演積分法。函數Y(z)可以看作是復數

Z平面上的勞倫級數，級數的各項系數可利用積分關系

求出：

1

Z-1[y(z)]=y(左T)=「恒⑵―力

2引*

?積分路徑C應包圍被積式所有極點，據留數定理知

n

歹(左T)=ZRes[y(z)z"，》

i=\

m(11Qi-1

=V\-------lim——p[(z-p^1Y1>

臺[(%-1)"f

54

計算機控制系統(tǒng)—

4.2.3Z變換與Z反變換方法

所以可利用留數計算法求得Z反變換的數學解析式:

zTy（z）]=y%T）

n

="Res[y（z）z"L,

其中：加為y（z）彼此不相等極點的個數；

Pi為y（z）z'T中不相等的極點；

%為極點”相重的個數。

55

4.3脈沖傳遞畫數

及其控制系統(tǒng)描述

7

4.3脈沖傳遞函數及其控制系統(tǒng)描述

?本節(jié)系統(tǒng)地介紹采用z變換求解線性微分方程

的方法。

?同時，類似于線性連續(xù)系統(tǒng)的傳遞函數描述，

引入脈沖傳遞函數（Z傳遞函數）的概念，用它

來描述離散線性系統(tǒng)。

?因此，線性離散系統(tǒng)既可以用差分方程來描

述，也可以用脈沖傳遞函數描述。

計算機控制系統(tǒng)_HYL

4.3脈沖傳遞函數及其控制系統(tǒng)描述

?z變換作為一個數學工具，可直接用來求解線性

差分方程，在給定的初始條件下，可以求得差

分方程的解析解。

?與線性連續(xù)系統(tǒng)相似，脈沖傳遞函數也是用來

分析線性離散系統(tǒng)的有效方法

?首先要了解建立脈沖傳遞函數的方法，與系統(tǒng)

結構相關的關聯形式。

計算機控制系統(tǒng)_HYL

431用Z變換法求解差分方程

?拉普拉斯變換將微分方程轉變?yōu)榇鷶颠\算的

形式來方便微分方程的求解。

?同樣，z變換也可將差分方程轉變成代數運

算的形式，大大地簡化差分方程的求解過

程，方便了離散系統(tǒng)的分析和綜合。

計算機控制系統(tǒng)_HYL

431用Z變換法求解差分方程

?用z變換求解差分方程的步驟如下：

?(1)、首先對差分方程作z變換，主要涉及到z變換

的滯后定理和超前定理；

?(2)、然后利用初始條件或求出的伏0)/(與,一，

并將它們帶入Z變換式；

計算機控制系統(tǒng)_HYL

計算機控制系統(tǒng)—

4.3.1用Z變換法求解差分方程

?(3)、整理所得差分方程的Z變換一般式，即可求出丫(z)：

1m1m-17

b°z+.ZH------卜靡

y(z)一~

aQz+axzHv-an

?(4)、由上式觀察得出y(z)的最高階數zr(q=n—m),從

而可判斷出：

xo)=o,xn=。，…，-q—i)7i=o

?(5)利用長除法、部分分式法或留數法作Z反變換，求出

差分方程的解貝依)。

61

計算機控制系統(tǒng)—

431用Z變換法求解差分方程

?例：求解差分方程(后向差分)

y(kT)-4y(kT—T)+3y(kT-27)=S(kT-2T)

其中B為數字脈沖，BE)[1“二°)

[0(kw0)

?[解]：運用滯后定理對差分方程進行Z變換

r(z)-4z-1r(z)+3Z-2Y(Z)=Z-2-1

Y(z)=——=-----------

l-4z-1+3z-27(z-3)(z-l)

?觀察上式，分子為零階，分母為二階，y(z)的最高階數Z—2,故可推知

系統(tǒng)應有：y(0)=0,(7)=0,用留數法作Z反變換得:

k-\k—1

Z

y(kT)=lim(z-3)-----------+lim(z-1)

—3(z-3)(z-l)z—>l(Z—3)(z—1)

=0.5⑶i-0.5⑴-=-3k-0.562

_____________________6_______

計算機控制系統(tǒng)_HY

431用Z變換法求解差分方程

?可以觀察到，上面的例題中，由y(AT)的解析式：

1

y(kT)=-3k-0.5

6

?當&=0,1,…時，可以計算得y(0)=—1/3,了⑴=0,八2)=1,…

?顯然，其中的9(0)二一1/3與前面分析所得的y(0)=0相矛盾，這表明

的解析式只確定出22時的義助值，前面若干周期的實際結果，由丫⑵

的分析結果決定。

?令采樣周期丁=1,k=2,代入原差分方程：

y(kT)-4y(kT-T)+3y(kT-2T)=3(kT-2T)

差分方程變?yōu)椋?/p>

六2)—4底1)+3>(0)=8(0)

?將/(。)=0,y(l)=。,y(2)=1代入上式，可驗證方程是成立的一

63

計算機控制系統(tǒng)—

431用Z變換法求解差分方程

?例：求解差分方程(前向差分)

y(kT+2T)-4.(-7+T)+3y(kT)=d(kT)(y(kT)=0,kWO)

其中b為數字脈沖，b(5)=P"二°)

[0(kw0)

?[解]：利用超前定理對差分方程進行Z變換

Z俗(左T)]=1Z2Y(Z)-z2y(0)-zy(T)-4[zY(z)-zy(0)]+3y(z)=1

由于已知：y(kT)=0,k<0將：.(0)=0,左=-1代入差分方程

可求得:y(T)=0

]1

代入Z變換式有：Y(z)=

z2-4z+3(z—3)(z—1)

z"lZ'T

y(kT)=lim(z-3)-----------------1-lim(z-1)---------------

z-3(z—3)(z—1)z-i(z_3)(z_i)

1

=0.5(3)J-0.5⑴J_3左_05

6

計算機控制系統(tǒng)—

4.3.2Z傳遞函數

一、z傳遞函數的定義

?Z傳遞函數是分析線性離散系統(tǒng)的重要工具。

?線性連續(xù)系統(tǒng)分析中，我們已經定義了傳遞函數：零初始

條件下系統(tǒng)輸出的拉氏變換與輸入的拉氏變換之比，即：

L[y(t)]丫⑹

G(s)=------=----

加⑺]RG)

?零初始條件：時系統(tǒng)的輸入量外⑺和輸出量j⑺以及它

們的各階導數均為零。

65

計算機控制系統(tǒng)—

4.3.2Z傳遞函數

、Z傳遞函數的定義

?類似地可以給出，線性離散系統(tǒng)z傳遞函數的定義：

?零初始條件下系統(tǒng)輸出脈沖序列的Z變換y(z)與輸入脈

沖序列的Z變換R(z)之比，即

Z[y(kT)]r(z)

Cr(Z)=--------二-----

Z[r(kT)]R(z)

?零初始條件：輸入量序列外(一7),外(一27),...和輸出量序

列V(—7)/(—27),…均為零：

66

計算機控制系統(tǒng)_HYC、

4.3.2Z傳遞函數

'Z傳遞函數的定義

?Z傳遞函數也稱為脈沖傳遞函數(pulsetransferfunction)。Z

傳遞函數G(z)反映了系統(tǒng)的固有特性，它僅取決于描述離

散系統(tǒng)的差分方程。

y(kT)

>(z)

r(kT)

G(s)

R(s)'A(z)Hs)y(kT)=Z-l[Y(z)]

=ZTG(Z>R(Z)]

系統(tǒng)或環(huán)節(jié)的"專遞函數67