离散系统线性二次型最优控制的算法设计.docxVIP

一种离散系统线性二次型最优控制的算法设计

姓名:专业:学号:

一种离散系统线性二次型最优控制的算法设计

一：背景意义

对于许多的控制系统，为得到满意的控制效果，需根据建立的

系统数学模型，选择一个容许的控制规律，在一定的条件下，使得控制系统在完成所要求的控制任务时，使给定的某一性能指标达到最优值，极小值或极大值，以使某一种性能指标为最小，实现最优控制.常用的性能指标有积分型性能指标如最小时间控制和最小能量控制；末值型性能指标如机床工作台移动准确停止控制和复合型性能指标等.线性二次型最优控制是一种常用的最优

控制系统设计方法.这种方法中的性能指标是对象状态与控制输入的二次型函数，在线性系统的约束条件下，选择控制输入使得二次型函数达到最小。

二：模型描述

离散线性定常系统：

，x(k+1)=Ax(k)+Bu(k)

TOC\o1-5\h\zy(k)=Cx(k)Du(k) (1)

式中，x(k)为n维状态向量；u(k)为p维控制向量;A为nxn非奇异阵；B为nxp矩阵，当其可控性矩阵的秩为n,选择完全可控线性离散系统的性能指标为：

1N4「 ]

J=2xT(N)Sx(N)+2瓦xT(k)Qx(k)+uT(k)Ru(k)] (2)

2k=0

Q为nxn维正定或半正定实对称矩阵；R为pxp维正定实对

称矩阵;S为nxn维正定或半正定实对称矩阵

若选择最优反馈距阵为

K(k)二R」BT(AT)」P(k)—Ql (3)

则其对应的最优控制序列和最优性能指标分别表示为公式 (4)和

(5)：

u(k)二-K(k)x(K)二-R—1BT(AT)—1〔P(k)-QX(k)

(4)

(k二0,1,2.…N-1)

J=1xT(0)P(0)x(0)

2 / (5)

其中P(k)二QATP」(k1)BR_1Bt]JA

3、最优控制序列的确定

令N—?8卡则桌统最优控制[Q解为稳态解,系统性能指标变人：

八当￡P⑷创約+/⑷也⑷]

Zk-Q

K(k)变为常数増益矩阵：

⑹(7)(8)K=IR+BPR\bW

⑹

(7)

(8)

P(M变为常数矩阵；

P^Q±AJ[Pl

对应的最优控制用列为：

u\k)=-Kx(k)=-{R+BtPBY[B7PAx(k)

闭环系统的状态方程为:

x(k-1)=Ax(k)+Bu(k)=(/+fi/?B7P)1Ax(k)(9)其戢优性能指标仍为公式(5)

氣实例仿真与结论

某伺服系统动态结构图如图1所/J?

(10)X(A4^1)=?(*)+*?(*)

(10)

Xi)=cx(A)

图I?何服系统动态结构芷

山图1,可得

(u(k)=kyV^k)-k2x(k)

Iv(A)=rx(k)一y⑻+咻一1) (11)

山公式(10)和(11),有

v(t+1)=-cax(k)+v(k)一cbu(k)+r(k+1) (12)

山公式(11)和(12),可得跆阵:

x(k+1)咻+1)_a-caOTx(i)1

x(k+1)咻+1)_

a

-ca

OTx(i)

1咻)

h

-ca

伙)十

卜伙十1)

(13)

令Xc(k)=x(k卜x(8),ve(k)=v(k)-W00),uc(k)=u(k11(°°)

兀伙+1)a

兀伙+1)

a

-ca

：］［常 檢)04)

III公式(11),有：叫⑷=広叫伙)-S) (15)

令鬲⑷=兀⑷,2仏)=Ve(*),W(i)=叫⑹=

一伙2X”)

兀2

一伙2

為伙+1)勺伙+1)

為伙+1)

勺伙+1)

a

-ca

(16)

对于单位阶跃给宦输入，冇:

工仏十1)咻+1)

\-cbky

/iy(k)=cx(k)=[coj

Xik)

吨)

选择性能*f标为；J=

(k)Qx(k)+w*(k)Rw⑹]

100 0

选择0= 0 1 ,R二|】]?

三：程序及仿真图

clear;

a=2;b=0.5;c=1;d=0;Q=[1000;01];R=[1];

A=[a0;-c*a1];B=[b;-c*b];KX=dlqr(A,B,Q,R);

k仁-KX(2);k2=KX(1);axc=[(a-b*k2)b*k1;(c*a+c*b*k2)(1-c*b*k1)];

bxc=[0;1];cxc=[10];dxc=[0];dstep(axc,bxc,cxc,dxc,1,100)

计算出的稳态最优反馈增益矩阵为

K=[3.8785-0.1743]

图2