程序运行结果如下: K =0.4142 0.6104 0.1009 同时得到闭环阶跃响应曲线,如图1-2所示。
图1-2 闭环系统阶跃响应曲线 由图1-1和图1-2知,经最优输出反馈后,闭环系统阶跃响应曲线与经最优状态反 馈后的阶跃响应曲线很接近。
三、离散系统线性二次型最优控制
下面对离散系统线性二次型最优控制进行详细介绍。
其中,A为系统的状态矩阵;B为系统的输出矩 阵;Q为给定的半正定实对称常数矩阵;R为给 定的正定实对称常数矩阵;N代表更一般化性 能指标中交叉乘积项的加权矩阵;K为最优反馈 增益矩阵;S为对应Riccati方程的唯一正定解P (若矩阵A-BK是稳定矩阵,则总有正定解P存 在);E为矩阵A-BK的特征值。
1000 Q= 取 0
,R=1。 用MATLAB函数dlqr()来求解最优控制器,给出程序清 单如下: %求解最优控制器 a=2;b=1;c=1;d=0; Q=[1000,0;0,1]; R=1; A=[a,0;-c*a,1]; B=[b;-c*b]; Kx=dlqr(A,B,Q,R) k1=-Kx(2);k2=Kx(1); axc=[(a-b*k2),b*k1;(-c*a+c*b*k2),(1-c*b*k1)]; bxc=[0;1];cxc=[1,0];dxc=0; dstep(axc,bxc,cxc,dxc,1,100)
•
1.LQG最优控制原理 最优控制原理
x(t ) = Ax(t ) + Bu (t ) + Gw(t ) 假设对象模型的状态方程表示为:
y (t ) = Cx(t ) + v(t )
T
式中,ω(t)和ν(t)为白噪声信号,ω(t)为系统干扰噪声,ν(t)为传感器带来的 量测噪声。假设这些信号为零均值的Gauss过程,它们的协方差矩阵为: