數(shù)學(xué)模型的相互轉(zhuǎn)換

數(shù)學(xué)模型的相互轉(zhuǎn)換第1頁，課件共30頁，創(chuàng)作于2023年2月2.1連續(xù)系統(tǒng)數(shù)學(xué)模型

2.1.1表達形式描述控制系統(tǒng)的主要模型有微分方程、狀態(tài)空間表達式等形式的時域描述法和用傳遞函數(shù)描述的頻域描述法。即對于一個連續(xù)的控制系統(tǒng)，數(shù)字仿真常用的數(shù)學(xué)模型一般有3種表示方式：◆直接用微分方程描述；◆用傳遞函數(shù)描述；

多項式形式零極點形式◆狀態(tài)方程描述；這三種描述方式是可以相互轉(zhuǎn)換的。

第2頁，課件共30頁，創(chuàng)作于2023年2月(1).微分方程設(shè)連續(xù)系統(tǒng)的輸出量為y(t)，輸入量為u(t)，采用微分方程的形式來表示的系統(tǒng)數(shù)學(xué)模型一般式可描述如下：（2.1）上式中，為常數(shù)。(2).傳遞函數(shù)第3頁，課件共30頁，創(chuàng)作于2023年2月對式（2.1）等號兩邊逐項進行拉氏變換，并考慮系統(tǒng)輸出、輸入及其各階導(dǎo)數(shù)的初值均為零，可得到

（2.2）式中，-系統(tǒng)輸出的拉氏變換；

-系統(tǒng)輸入的拉氏變換；可得系統(tǒng)的傳遞函數(shù)為：(2.3)第4頁，課件共30頁，創(chuàng)作于2023年2月微分方程或傳遞函數(shù)是用系統(tǒng)的輸入、輸出之間的關(guān)系來描述系統(tǒng)的，表示了系統(tǒng)的外部特征，所以稱其為外部模型。用微分方程表示的系統(tǒng)可以是非線性或線性系統(tǒng)，而對于傳遞函數(shù)表示的系統(tǒng)，只適用于單輸入-單輸出的線性定常系統(tǒng)，所以傳遞函數(shù)的模型表示有一定的局限性。

(3).狀態(tài)空間表達式狀態(tài)空間表達式可以由兩個途徑獲得，由微分方程或系統(tǒng)結(jié)構(gòu)方框圖導(dǎo)出，這里對微分方程推導(dǎo)作簡單說明。設(shè)系統(tǒng)由不含輸入量導(dǎo)數(shù)項的n階微分方程表示：(2.4)第5頁，課件共30頁，創(chuàng)作于2023年2月定義n個狀態(tài)變量為，且令寫出各個狀態(tài)變量的一階微分方程形式第6頁，課件共30頁，創(chuàng)作于2023年2月將上述n個一階微分方程寫成矩陣向量形式為(2.5)上式稱為狀態(tài)空間表達式，其中A、B、C為系數(shù)矩陣，x為狀態(tài)變量。第7頁，課件共30頁，創(chuàng)作于2023年2月2.1.2.數(shù)學(xué)模型的相互轉(zhuǎn)換由于要解決的控制問題所需的數(shù)學(xué)模型與所給定的已知數(shù)學(xué)模型往往是不一致的，不同的應(yīng)用場合需要對控制系統(tǒng)的數(shù)學(xué)模型進行轉(zhuǎn)換

（1）微分方程轉(zhuǎn)化為傳遞函數(shù)和狀態(tài)空間表達式例2.1已知某控制系統(tǒng)的微分方程為將其分別表示為傳遞函數(shù)、一階微分方程組和狀態(tài)空間描述。第8頁，課件共30頁，創(chuàng)作于2023年2月解：①將給定系統(tǒng)微分方程的兩端取拉氏變換，并令初始值為零，則可用以下傳遞函數(shù)表示根據(jù)傳遞函數(shù)定義有第9頁，課件共30頁，創(chuàng)作于2023年2月③按照狀態(tài)空間描述，將各變量和系數(shù)矩陣表達為②由于是二階導(dǎo)數(shù)，可以引入兩個狀態(tài)變量，將給定的二階微分方程寫成一階微分方程組形式第10頁，課件共30頁，創(chuàng)作于2023年2月（2）傳遞函數(shù)轉(zhuǎn)換成狀態(tài)空間表達式轉(zhuǎn)換采用的方法是狀態(tài)變量圖法，用基本模擬單元替代系統(tǒng)的傳遞函數(shù)得到的圖形式系統(tǒng)結(jié)構(gòu)圖，在系統(tǒng)結(jié)構(gòu)圖上標(biāo)上狀態(tài)變量的圖形是狀態(tài)變量圖，然后再求出狀態(tài)空間表達式。

圖2-1積分器的系統(tǒng)結(jié)構(gòu)圖和狀態(tài)變量圖由狀態(tài)變量圖根據(jù)積分器的輸入、輸出關(guān)系寫出：輸出方程狀態(tài)方程對于初始條件為零的積分器

第11頁，課件共30頁，創(chuàng)作于2023年2月對于帶反饋的積分器，其傳遞函數(shù)為圖2-2帶反饋積分器的狀態(tài)變量圖由積分器輸入、輸出關(guān)系得到第12頁，課件共30頁，創(chuàng)作于2023年2月從上面得到由系統(tǒng)結(jié)構(gòu)圖到狀態(tài)變量圖并到處狀態(tài)空間表達式的步驟如下：根據(jù)系統(tǒng)的傳遞函數(shù)，畫出系統(tǒng)結(jié)構(gòu)圖，n階系統(tǒng)有n個積分器；把積分器輸出處定為狀態(tài)變量x，積分器輸入處為狀態(tài)變量微分，并把狀態(tài)變量x，和狀態(tài)變量微分分別標(biāo)在積分器輸入和輸出處，得到狀態(tài)變量圖；根據(jù)積分器輸入、輸出的方程寫出系統(tǒng)的狀態(tài)方程和輸出方程。對于高階、復(fù)雜系統(tǒng)采用級聯(lián)法、并聯(lián)法和串聯(lián)法得到代表實際系統(tǒng)傳遞函數(shù)的系統(tǒng)結(jié)構(gòu)圖及相應(yīng)的狀態(tài)變量圖，依據(jù)同樣方法求得狀態(tài)空間表達式。第13頁，課件共30頁，創(chuàng)作于2023年2月對于一個三階系統(tǒng)的傳遞函數(shù)，

(2.6)（1）級聯(lián)法

采用如下步驟進行：◆用傳遞函數(shù)的最高階次除以傳遞函數(shù)分子分母多項式得到◆利用信號流圖法畫出該系統(tǒng)的信號流圖及系統(tǒng)狀態(tài)變量圖

第14頁，課件共30頁，創(chuàng)作于2023年2月圖2-3系統(tǒng)信號流圖圖2-4系統(tǒng)狀態(tài)變量圖第15頁，課件共30頁，創(chuàng)作于2023年2月根據(jù)積分器輸入、輸出關(guān)系得到如下方程◆寫成矩陣表達式

A、B、C為系數(shù)矩陣第16頁，課件共30頁，創(chuàng)作于2023年2月其中，分析系數(shù)矩陣A、B、C可見：系數(shù)矩陣A是一個方陣，以I表示行號，J表示列號，最末一行元素和傳遞函數(shù)分母多項式系數(shù)按s0升冪排列的負值一一對應(yīng)，其余各行的元素在J=I+1時為1，其他全部為0；系數(shù)矩陣B是一個單列矩陣，最后一行元素為1，其余為零；系數(shù)矩陣C是一個單行矩陣，各列元素與傳遞函數(shù)分子多項式系數(shù)按s0升冪排列值相同。第17頁，課件共30頁，創(chuàng)作于2023年2月推廣到n階方程，系數(shù)矩陣

A、B、C分別為

這種形式的矩陣稱為可控標(biāo)準(zhǔn)型狀態(tài)表達式。(2.7)第18頁，課件共30頁，創(chuàng)作于2023年2月例2.2已知某控制系統(tǒng)的傳遞函數(shù)為利用式（2.7）將系統(tǒng)模型轉(zhuǎn)換為狀態(tài)空間描述和一階微分方程組描述。

解：①這是一個三階系統(tǒng)，將給定的系統(tǒng)傳遞函數(shù)按照狀態(tài)空間描述的系數(shù)矩陣A、B、C的對應(yīng)關(guān)系，依據(jù)式（2.7）可得各系數(shù)矩陣如下第19頁，課件共30頁，創(chuàng)作于2023年2月組合為狀態(tài)空間描述有②將上述矩陣展開即可得到系統(tǒng)模型的一階微分方程組表示形式第20頁，課件共30頁，創(chuàng)作于2023年2月（Ⅱ）并聯(lián)法并聯(lián)法的思路是把高階系統(tǒng)的傳遞函數(shù)轉(zhuǎn)變成若干個一階環(huán)節(jié)傳遞函數(shù)之和，如將式（2-6）表示成下式形式然后，對各個一階環(huán)節(jié)的傳遞函數(shù)畫出系統(tǒng)結(jié)構(gòu)圖，并標(biāo)上狀態(tài)變量后，得到如圖2-5所示的狀態(tài)變量圖。圖2.5并聯(lián)法系統(tǒng)狀態(tài)變量圖第21頁，課件共30頁，創(chuàng)作于2023年2月寫成矩陣表達式由圖2-5可到狀態(tài)空間表達式為第22頁，課件共30頁，創(chuàng)作于2023年2月其中利用級聯(lián)法和并聯(lián)法得到的狀態(tài)空間表達式的系數(shù)矩陣A、B、C是不同的。可以用特征方程判斷特征值是否相同，確定幾個狀態(tài)空間表達式是否屬于同一個外部模型，因為同一系統(tǒng)有相同的特征值。

第23頁，課件共30頁，創(chuàng)作于2023年2月（3）狀態(tài)空間表達式轉(zhuǎn)換成傳遞函數(shù)已知系統(tǒng)狀態(tài)空間表達式①對上式兩邊取拉氏變換②在零初始條件下整理上式并在等式兩邊乘以單位矩陣（2.9）

（2.10）第24頁，課件共30頁，創(chuàng)作于2023年2月③將X(s)代入式（2.10）得到

可得（2.11）

在系數(shù)矩陣已知的情況下，根據(jù)式（2.11）就可以求出狀態(tài)空間表達式所對應(yīng)的外部模型的傳遞函數(shù)。

在MATLAB控制系統(tǒng)工具箱中提供了大量的控制系統(tǒng)模型相互轉(zhuǎn)換的函數(shù)，如表2-1所示。第25頁，課件共30頁，創(chuàng)作于2023年2月表2-1數(shù)學(xué)模型轉(zhuǎn)換函數(shù)及其功能函數(shù)名函數(shù)功能ss2tf將系統(tǒng)狀態(tài)空間模型轉(zhuǎn)換為傳遞函數(shù)模型ss2zp將系統(tǒng)狀態(tài)空間模型轉(zhuǎn)換為零極點增益模型tf2ss將系統(tǒng)傳遞函數(shù)模型轉(zhuǎn)換為狀態(tài)空間模型tf2zp將系統(tǒng)傳遞函數(shù)模型轉(zhuǎn)換為零極點增益模型zp2ss將系統(tǒng)零極點增益模型換為狀態(tài)空間模型zp2tf零極點增益模型換為傳遞函數(shù)模型第26頁，課件共30頁，創(chuàng)作于2023年2月[z,p,k]=tf2zp(num,den)[A,B,C,D]=tf2ss(num,den)[num,den]=zp2tf(z,p,k)[A,B,C,D]=zp2ss(z,p,k)[num,den]=ss2tf(A,B,C,D)[z,p,k]=ss2zp(A,B,C,D)第27頁，課件共30頁，創(chuàng)作于2023年2月對線性定常系統(tǒng)，式中s的系數(shù)均為常數(shù)，且a1不等于零，這時系統(tǒng)在MATLAB中可以方便地由分子和分母系數(shù)構(gòu)成的兩個向量唯一地確定出來，這兩個向量分別用num和den表示。

num=[b1,b2,…,bm,bm+1] den=[a1,a2,…,an,an+1]

注意：它們都是按s的降冪進行排列的。MATLAB表示數(shù)學(xué)模型舉例連續(xù)系統(tǒng)的傳遞函數(shù)模型第28頁，課件共30頁，創(chuàng)作于2023年2月在MATLAB中零極點增益模型用[z,p,K]矢量組表示。即：z=[z1,z2,…,zm]p=[p1,p2,...,pn]K=[k]函數(shù)tf2zp()可以用來求傳遞函數(shù)的零極點和增益。零極點增益模型K為系統(tǒng)增益，zi為零點，pj為極點零極點模型實際上是傳遞函數(shù)模型的另一種表現(xiàn)形式，其原理是分別對原系統(tǒng)傳遞函數(shù)的分子、分母進行分解因式處理，以獲得系統(tǒng)的零點和極點的表示形式。第29頁，課件共30頁，創(chuàng)作于2023年2月用法舉例：1）已知系統(tǒng)狀態(tài)空間模型為：》A=[01;-1-2];B=[0;1];》C=[