第8章MATLAB在非線性系統中的應用

2022-12-26 20:30:05 字數 3110 閱讀 1166

1、狀態反饋與極點配置

狀態反饋是指從狀態變數到控制端的反饋,如圖3所示。

設原系統動態方程為:

引入狀態反饋後,系統的動態方程為:

圖3 狀態反饋

狀態反饋極點配置問題就是:通過狀態反饋矩陣k的選取,使閉環系統式(4.2)的極點,即的特徵值恰好處於所希望的一組給定閉環極點的位置上。

線性定常系統可以用狀態反饋任意配置極點的充分必要條件是:該系統必須是完全能控的。所以,在實現極點的任意配置之前,必須判別受控系統的能控性。

例8.1 已知有以下系統

希望將閉環系統的極點配置在s1,2,3,4=-1,-2,-1+j,-1-j

解:參考程式:

a=[-5 8 0 0;-4 7 0 0;0 0 0 4;0 0 -2 6];

b=[4;-2;2; 1];

p=[-1,-2,-1+sqrt(-1),-1-sqrt(-1)];

k=place(a,b,p)

執行結果:

??? error using ==> place

can't place eigenvalues there.

說明:對於不完全可控的系統,解決辦法:將系統分解成可控部分和不可控部分,可控部分可以將極點配置到任意的位置。

2、輸出反饋與極點配置

輸出反饋指從輸出端到狀態變數導數的反饋,如圖4所示。

設原系統動態方程為:

引入輸出反饋後,系統的動態方程為:

圖4 輸出反饋

極點配置演算法:

ackerman演算法

– k=acker (a,b,p)

魯棒演算法

– k=place (a,b,p)

– 函式place( )不適用於含有多重期望極點的配置問題。

例8.2 已知物件模型

如何將閉環系統的極點配置在s1,2,3,4= -1,-2,-1±j?

參考程式:

>> a=[0 1 0 0;0 0 -1 0;0 0 0 1;0 0 11 0];

b=[0;1;0;-1];

eig(a)' %特徵值

p=[-1, -2, -1+sqrt(-1), -1-sqrt(-1)];

k=place (a,b,p) %極點配置

eig(a-b*k)'

執行結果:

ans =

0 0 3.3166 -3.3166

k = -0.4000 -1.0000 -21.4000 -6.0000

ans =

-2.00001.0000 - 1.0000i -1.0000 + 1.0000i -1.0000

可以看出,受控系統的極點位置確定位於0,0,3.3166,-3.3166,即該受控系統是不穩定的。

應用極點配置技術,我們可以將系統的閉環極點配置到某些期望的位置上,從而使得閉環系統得到穩定,並同時得到較好的動態特性。

例8.3 已知物件模型

利用matlab實現將其中的兩個極點配置到。

本題程式參考:

a=[0,1,0,0;0,5,0,0;0,0,-7,0;0,0,0,-8];

b=[1;1;3;4];c=[0,5,0,8];

p=eig(a)'

p1=p

p1(1:2)=[-1,-2]

k1=place(a,b,p1)

eig(a-b*k1)'

執行結果:

p = 0 5 -7 -8

p1 =

0 5 -7 -8

p1 =

-1 -2 -7 -8

k1 =

-0.5000 8.5000 0.0000 0

ans =

-8.0000 -1.0000 -2.0000 -7.0000

如果需要配置三個極點,則例如

p2=p;p2(1:3)=[-1,-2,-3];

k2=place(a,b,p2),eig(a-b*k2)',

如果受控系統並不是完全能控的,那麼我們只能將其中能控的極點配置到指定的位置,稱為部分極點配置問題,而實現部分極點配置的前提條件是:受控系統沒有重極點。

結論:對受控系統中不可控的模態而言,若其是不穩定的,則我們不可能通過狀態反饋的方法將其變成穩定的模態。

例8.4 已知物件模型

(1) 如果我們想將閉環系統的極點配置到-1,-2,-3,利用matlab設計控制器,並繪出閉環系統的階躍響應曲線。(說明:用兩種方法配置極點)

參考程式一:

>> a=[-0.3 0.1 -0.05;1 0.1 0;-1.5 -8.9 -0.05 ];

b=[2;0;4]; c=[1 2 3];

p=[-1, -2, -3]; k=place (a',c',p)'

執行結果:

k = -0.1731

-0.3681

2.2197

參考程式二:

a=[-0.3 0.1 -0.05;1 0.1 0;-1.5 -8.9 -0.05 ];

b=[2;0;4]; c=[1 2 3];

p=[-1, -2, -3]; k=acker (a',c',p)'

執行結果:

k = -0.1731

-0.3681

2.2197

*(2) 如果想將閉環系統的所有極點均配置到-1,怎樣設計控制器?

說明:若希望極點均配置到-1,即多重極點配置,應該使用acker( )函式。

例8.5 以試驗三磁懸浮模型為研究物件,如何通過狀態反饋配置極點,改善系統的效能。

num=1

den=[1 40 20.5 -157552]

[a,b,c,d]=tf2ss(num,den)

g=ss(a,b,c,d)

[y,t,x]=step(g)

p=[-2;-1+sqrt(-1);-1-sqrt(-1)]

k=place(a,b,p)

eig(a-b*k)

cn=[b,a*b,a^2*b]

rank(cn)

a1=a-b*k

g1=ss(a1,b,c,d)

[y,t]=step(g1)

plot(t,y)

第4章非線性規劃張

第四章非線性規劃模型 第一節非線性規劃的例項與基本概念 一 非線性規劃的例項 例1 化學反應的平衡組成 設現有原料由種原子組成,各種原子數量依次為共生成種分子 產品 設生產數量 待求 依次為。設第種分子中含各種原子的數量依次為 所有產品中含第種原子數之和為 由熟知的質量守恆定律有 在一定的溫度 壓力...

第4章非線性規劃張

第四章非線性規劃模型 第一節非線性規劃的例項與基本概念 一 非線性規劃的例項 例1 化學反應的平衡組成 設現有原料由種原子組成,各種原子數量依次為共生成種分子 產品 設生產數量 待求 依次為。設第種分子中含各種原子的數量依次為 所有產品中含第種原子數之和為 由熟知的質量守恆定律有 在一定的溫度 壓力...

第二章線性系統的狀態空間描述

2 併聯實現 1 只含單實極點的情況 設可分解為 式中為n階系統的單實極點,則可化為對角標準型。那麼傳遞函式可展成 式中 取狀態變數 整理後有 即狀態方程為 又有 即輸出方程為 向量 矩陣形式為 對角形動態方程的狀態變數圖為 由於對角形動態方程的狀態變數圖 例2.2.9 已知系統傳遞函式為,試求狀態...