极点配置自适应控制器的设计
具有抑制噪声性质的极点配置自校正控制器设计
摘要：本文首先介绍了极点配置自校正控制器设计中的一些不足之处，然后提
出了一种极点配置自校正控制算法.它具有以下三个特点:（1）除了配置系统的闭环极点外，还可抑制噪声对系统输出的干扰；（2）与其他极点配置自校正控制算法相比，本算法所需的计算量大为减少；（3）理论上可得到算法的收敛性结果.该算法已在实际系统中得到了成功的应用。
关键字：极点配置 自校正 仿真 一、前言
极点配置自校正控制器在最近几年里得到了很大的发展，出现了许多自校正算法，但现有的各种自校正控制算法还存在一些不够完善之处。首先，由于这种控制器的控制目的是配置系统的闭环极点，所以，在控制器的设计过程中，对如何减少噪声对系统输出的干扰这个问题没有特别加以考虑；第二，计算量比较大，难以实时控制；第三，理论上不够完善，如现有的各种算法都没有收敛性证明。所以说，极点配置自校正控制器的理论和应用都不太成熟，有待于进一步发展。
二、算法
设所研究的系统由下面方程描述。
(1)
)e (t )C (q )u (t )(q q (t))(-1-1-k 1-+=B y q A 其中(t)y 为输出量，(t)u 为控制量，(t)e 为零均值白噪声且与t 时刻以前（不包括
t 时刻）的控制(t)
u 无关。
c
c b
b a a n n n n n n q
c q c q C q
b q b b q B q
a q
a A -1
-11--1
-101
--1
-11
-1)()(1)q (+???++=+???++=+???++=
)(1
-q A ，)(1
-q C 为稳定多项式，0
0≠b ，1-q 为延后一步算子。
对上述系统如采用如下控制律
)
()
()
()()
F(q )
G(q
-(t)1
1111
--1t y q F q
G t y u r --+
= （2)
其中(t)r y 为参考输入，)(1-q G ，)(1-q F ，)(1-1q G ，)(q -11F 为四个适当次数的多项。则可得闭环方程
(4)
(t )u (t )u u (t )(3) (t )(t )y (t )2121+=+=y y
其中
(5)
(t))
)G(q
(q q )(q )()
()()
(F )(q q (t)1
-1
-k
-1
-1
-1
-1-11
-1-1
-k
1r y B F q A q F q G q B y +?
=
为跟踪参考输入部分；
(6) (t))
)G(q (q q ))F(q A(q
)
)F(q C(q
(t)1
-1
-k
-1
-1
--1
-1
2e B y +=
为噪声所引起的扰动部分
(7)
(t))
)G(q
(q q )(q )()
()()
(F )
A(q
(t)1
-1
-k
-1
-1
-1
-11-1
-1-11r y B F q A q G q F q u +?
=
为跟踪参考输入所需的控制量
(8) (t))
)G(q (q q ))F(q A(q
)
)C(q
G(q
-(t)1
-1
-k
-1
-1
--1
-1
2e B u ?+=
为抑制噪声干扰所需的控制量。
为对于参考输入实现极点配置，应使下式成立。
(9) (t))
T(q )
)Z(q
(q q (t)1
--1
-1
-k
1r y B y =
其中)(1-q T ，)(1-q Z 由设计者给定，它们决定了系统对于参考输入的极点和部分零点。这里我们假设系统的开环零点都是希望的闭环零点。如此假设不成立时，可采用如下两种方法：1）假设原系统为最小相位系统，我们可取B T T '=。由式（9）可知，此时闭环传递函数为)(q )/T'(q -1-1-Z q k ，从而实现了零极点配置。2）如原系统是非最小相位系统，可将B 分解为稳定部分+B 和不稳定部分-B ，取
+
=B T T '，此时闭环传递函数为'/--T Z B q k
，从而实现了部分零极点配置。
为减少噪声对系统输出的干扰，应使指标函数
})(t )u ~
()k)(t {(2
222λ++=y E J (10) 最小。其中λ~
为某一常数。
不加证明地我们选取如下的控制器参数，可达到上述两个目的。
)()(1
-1
-q G q G =
(11)
)C(q )(q F )(q )(-1
-1
-1
1
-λ+=B q F (12)
))F(q (q )(-1
-1
1
-1T q F =
(13)
)B (q
))A(q )()C(q
(q )(-1
-1
-1
-11-1+=λZ q G (14)
其中)(1-q F ，)(1-q G 由下面方程解得。
)C (q )(q )()(-1-1-1-1-=+G q q F q A k (15) deg 1-)(1-k q F =，deg k)-n 1,-(n max )(c a 1-=q G ，02/~
b λλ= 于是，我们可以得到下面的显式自校正控制算法。
1）在时刻t ，用递推增广最小二乘法辨识系统(1)的参数)(1-q A ，)(1
-q B ，
)(1-q C 得到估计值)(q ?-1t A ，)(q ?-1t B ，)(q ?-1t
C 。 2）由方程)(q ?)q (?q )(q F )(q ?-11
1-1-k -1t -11C G A =+ (16) 解得)(q -1t F ，)(q ?-1t G ，其中deg 1-k )(q -1=t F ，1)-n k,-(n max ?deg a
c 1-=q G t ，然后计算 )(q C ?)(q B ?)(q F )(q ?-1t -1t -1t -1λ+=t F (17) )(q F ?)T(q )(q -1
t
-1
-1
1=t F (18) ))(q B ?)(q A ?)((q C ?)Z(q )(q -1t
-1t -1t -1-11+=λt G (19) 式(17)，(19)中的λ的选取准则是使)B(q )(q -1-1+A λ成为一个稳定多项式。对于最小相位系统我们可取0=λ。
3）记由方程(t))
(q F ?)(q ?(t))(q F ?)(q G ?-(t)1-1t -1t
11-t -1t r y G y u +=得到的(t)u 为1M ，则取 t 时刻的控制量为
(t )u 11M )s g n (M {
M ?=
M
|M |M |M |12≤>如果如果
其中M 为一考虑工程上对控制量所加的限制而给定的正数。此控制器具有“结
合型”性质，即兼有极点配置自校正控制器和最优型自校正控制器的特点。对于参考输入，系统具有所希望的闭环极点，这有利于利用工程师们的实际经验设计出一个满意的跟踪系统。而对于噪声，控制器使(10)式的指标函数J 达到了最小, 从而减小了噪声对输出的影响，提高了跟踪精度。
本算法的另一优点是计算量小。一般的极点配置自校正控制器在每一采样周期内要解一个形如))C(q T(q ))G(q B(q )(q )(-1-1-1-1-11-=+F q A 的多项式方程，这相当于解一组k n n b a ++维的线性方程组。当原系统的阶数较高时，其计算量是很大的。在本算法中，这个方程被简化为)C(q )(q )()(-1-11-1-=+G q q F q A k ，我们只需用)(1-q A 除以)(1-q C ，得到k-1阶商式即为)q (1-F ，余式为)(1--q G q k ，另外方程
(12)，(13)，(14)也只需经过几次多项式乘法即可得到。所以，同其他极点配置自校正控制算法相比，本算法所需的计算量大为减少，从而提高了控制器的实时控制能力。
在本算法中，我们假设延时k 是已知的。如果k 事先并不知道时，我们可按k=1的情形设计控制器。这样做的结果有可能使系统失去抑制噪声的能力。但系统 的跟踪性质并不会受到影响。也就是说，当时延k 事先不知道时，控制器同样可达到极点配置的目的。
由方程(3)和方程(15)~(11)可知，系统的理想闭环方程为
(20) (t))
B(q )A(q )
F(q (t))
T(q )
)Z(q
(q q (t)1
-1--1
1
--1
-1
-k
e y B y r ++
=
λ
设在相同的参考输入和噪声干扰情况下，理想系统(20)的输出为(t)*y ，按上述自校正控制算法得到的输出为(t)y ，则有以下的收敛性定理
定理 对系统(1)，如采用本文所述的自校正控制算法，且满足 1）(t)r y 有界；
2）1/2-)(q 1/-1C 严格正实；
3）0/lim 1->∞
→t P t t ，t P 为递推增广最小二乘法中的增益矩阵；
4）λ的选择使)B(q )(q -1-1+A λ为稳定多项式。且1=ρ，ρ为递推增广最小二乘法中的遗忘因子；
5）对充分大的t 有|M u t <|则有
0y (t ))-(t )
(y lim *=∞
→t 三、仿真结果
设被控对象为如下(本例子中加入两人有色噪声)：
2)
-0.6065y(k 1)-1.6065y(k -(k)+y
1)-(k 0.5(k)4)-0.0902u(k 3)-(k 0.1065ξξ+++=u 式中，(k)ξ为方差为0.01的白噪声。
期望传递函数分母多项式为 2
--1
-1m 4966.01.3205z
-1)(z z
A +=
取初值0.001(0)?10)0(6==θ、I P ，遗忘因子1=ρ；期望输出(k)r y 为幅值
10的方波信号。极点配置自校正控制间接算法也可分两种情况进行仿真，仿真结果如下如所示：
(1)考虑系统零点不被对消的情况(仿真程序中V al=0.0)
(a)控制效果
(b)对象参数估计结果
以上(a)、(b)两图为极点配置间接自校正控制零点不对消的情况。仿真程序见附录
(2)考虑系统零点被对消的情况(仿真程序中取V al=0.9)
(c)控制效果
(d)对象参数估计结果
可见本算法在减少噪声干扰方面的效果是非常明显的。
四、参考文献：
[1]一种具有抑制噪声性质的极点配置自校正控制器郑勤中国科学技术大学
[2]自校正控制器袁著祉中北大学
[3]抑制未知确定性扰动的极点配置自适应控制刘贺平孙一康
北京科技大学自动化信息工程学院北京100083
[4]自校正控制器控制规律的快速解算陈慧清梁慧冰
[5]一种简单的多变量极点配置自适应控制方案吕民张成乾
[6]极点配置自校正PID控制算法及其应用武省陈述平王研张宏宇
基于极点配置的控制器设计与仿真
计算机控制理论与设计作业 题目：基于极点配置方法的直流调速系统的控制器设计
摘要 本文目的是用极点配置方法对连续的被控对象设计控制器。基本思路是对连续系统进行数学建模，将连续模型进行离散化，针对离散的被控对象，用极点配置的方法分别在用状态方程和传递函数两种描述方法下设计前馈和反馈控制器，并用MATLAB仿真。文中具体以直流调速系统作为研究对象，对直流调速系统的组成和结构进行了分析，把各个部分进行数学建模，求出其传递函数，组成系统结构框图，利用自控原理的知识对结构图化简，求出被控对象的传递函数和状态方程，进一步得将其离散化。第一种是通过极点配置设计方法的原理，用状态方程设计被控对象的控制律，因为直流调速系统存在噪声，实际状态不可测，故选择了全阶的观测器，又因为采样时间小于计算延时，所以选择了预报观测器。利用所学知识对此闭环系统设计前馈和反馈控制器[1]。第二种利用传统的离散传递函数，从代数多项式的角度进行复合控制器的设计，在保证系统稳定的情况下，分析系统的可实现性，稳定性，静态指标，动态指标，抗干扰等方面性能研究前馈反馈相结合控制器设计。重点是保证被控对象的不稳定的零极点不能被抵消。最后利用MATLAB的Simulink进行仿真，观察系统的输出的y和u和收敛性，并加入扰动看其抗干扰性能，得出结论。 经研究分析，对于直流调速系统，基于极点配置设计的前馈反馈相结合的控制器，具有良好的稳定性能和抗干扰性能。运行结果符合实际情况。 关键词：极点配置；状态方程；直流调速系统；代数多项式；Matlab；
1绪论 1.1论文的背景及意义 在工业生产和日常生活中，自动控制系统分为确定性系统和不确定性系统两类，确定性系统是指系统的结构和参数是确定的，确定的输入下，输出也确定的一类系统。确定性系统相对于不确定性系统而言的。在确定的系统中所用的变量都可用确切的函数关系来描述，系统的运动特性可以完全确定。以确定性系统为研究对象的控制理论称为确定性控制理论。本文以直流调速系统为研究对象，利用极点配置的设计方法，包括利用状态空间模型和传递函数模型分别描述线性系统，采用闭环极点为指标的控制器设计的理论和方法，设计出前馈和反馈控制器，组建闭环控制系统，用Matlab进行仿真可以逼真地还原出实际系统。 1.2 论文的主要内容 本文直流电机的调速系统的模型作为研究对象，利用线性系统极点配置的设计方法，设计前馈反馈控制器。论文研究的主要内容: (1)阅读学习国内外期刊文献，研究了极点配置的基本原理和Matlab的实现方法。 (2)系统的说明直流电机的系统结构和工作原理并分析，建立直流调速系统的数学模型，将其进行离散化，并讨论其传递函数与状态方程之间的关系。 (3)分析极点配置控制器的设计原理，利用状态方程设计控制器。 （4）将被控对象的传递函数离散化，利用传递函数模型设计控制器。 (4)在MATLAB中建立闭环直流调速系统的模型，根据闭环极点配置的设计步骤编写程序，用Simulink搭建仿真系统，对闭环直流调速系统的输出进行仿真分析。 (5)对仿真结果分析。将仿真结果与实际直流调速系统的阶跃响应的各项参数相比较，得出结论。
倒立摆状态空间极点配置控制实验实验报告
《现代控制理论》实验报告 状态空间极点配置控制实验 一、实验原理 经典控制理论的研究对象主要是单输入单输出的系统，控制器设计时一般需要有关被控对象的较精确模型，现代控制理论主要是依据现代数学工具，将经典控制理论的概念扩展到多输入多输出系统。极点配置法通过设计状态反馈控制器将多变量系统的闭环系统极点配置在期望的位置上，从而使系统满足瞬态和稳态性能指标。 １．状态空间分析 对于控制系统X = AX + Bu 选择控制信号为：u = ?KX 式中：X 为状态向量（ n 维）u 控制向量（纯量） A n × n维常数矩阵 B n ×1维常数矩阵 求解上式，得到 x(t) = (A ? BK)x(t) 方程的解为： x(t) = e( A?BK )t x(0) 状态反馈闭环控制原理图如下所示： 从图中可以看出，如果系统状态完全可控，K 选择适当，对于任意的初始状态，当t趋于无穷时，都可以使x(t)趋于0。 ２．极点配置的设计步骤 1) 检验系统的可控性条件。 2) 从矩阵 A 的特征多项式 来确定 a1, a2,……,an的值。 3) 确定使状态方程变为可控标准型的变换矩阵 T：T = MW 其中 M 为可控性矩阵， 4) 利用所期望的特征值，写出期望的多项式 5) 需要的状态反馈增益矩阵 K 由以下方程确定： 二、实验内容 针对直线型一级倒立摆系统应用极点配置法设计控制器，进行极点配置并用Matlab进行仿真实验。 三、实验步骤及结果 1.根据直线一级倒立摆的状态空间模型，以小车加速度作为输 入的系统状态方程为： 可以取1 l 。则得到系统的状态方程为： 于是有：
直线一级倒立摆的极点配置转化为： 对于如上所述的系统，设计控制器，要求系统具有较短的调整时间（约 3 秒）和合适的阻尼（阻尼比? = 0.5）。 2.采用四种不同的方法计算反馈矩阵 K。 方法一：按极点配置步骤进行计算。 1) 检验系统可控性，由系统可控性分析可以得到，系统的状态完全可控性矩阵的秩等于系统的状态维数(4)，系统的输出完全可控性矩阵的秩等于系统输出向量y 的维数(2)，所以系统可控。 倒立摆极点配置原理图 2) 计算特征值 根据要求，并留有一定的裕量（设调整时间为 2 秒），我们选取期望的闭环极点s =μi (i = 1,2,3,4) ，其中： 其中，μ 3,μ 4 使一对具有的主导闭环极点，μ 1 ,μ 2 位于 主导闭环极点的左边，因此其影响较小，因此期望的特征方程为： 因此可以得到： 由系统的特征方程： 因此有 系统的反馈增益矩阵为： 3) 确定使状态方程变为可控标准型的变换矩阵 T：T = MW 式中： M = 0 1.0000 0 0 1.0000 0 0 0 0 0.7500 0 5.5125 0.7500 0 5.5125 0 W = 0 -7.3500 -0.0000 1.0000 -7.3500 -0.0000 1.0000 0 -0.0000 1.0000 0 0 1.0000 0 0 0 于是可以得到： T = -7.3500 -0.0000 1.0000 0 0 -7.3500 -0.0000 1.0000 0 -0.0000 0.7500 0 -0.0000 0 -0.0000 0.7500 T’= -7.3500 0 0 -0.0000 -0.0000 -7.3500 -0.0000 0 1.0000 -0.0000 0.7500 -0.0000 0 1.0000 0 0.7500
实验 6 极点配置与全维状态观测器的设计(优.选)
实验 6 极点配置与全维状态观测器的设计 一、实验目的 1. 加深对状态反馈作用的理解。 2. 学习和掌握状态观测器的设计方法。 二、实验原理 在MATLAB 中，可以使用acker 和place 函数来进行极点配置，函数的使用方法如下：K = acker(A,B,P) A，B为系统系数矩阵，P为配置极点，K为反馈增益矩阵。 K = place(A,B,P) A，B为系统系数矩阵，P为配置极点，K为反馈增益矩阵。 [K,PREC,MESSAGE] = place(A,B,P) A，B为系统系数矩阵，P为配置极点，K为反馈增益矩阵，PREC 为特征值，MESSAGE 为配置中的出错信息。 三、实验内容 1.已知系统 （1）判断系统稳定性，说明原因。 （2）若不稳定，进行极点配置，期望极点：-1，-2，-3，求出状态反馈矩阵k。 （3）讨论状态反馈与输出反馈的关系，说明状态反馈为何能进行极点配置？ （4）使用状态反馈进行零极点配置的前提条件是什么？ 1. （1） （2） 代码： a=[-2 -1 1;1 0 1;-1 0 1]; b=[1,1,1]'; p=[-1,-2,-3]'; K=acker(a,b,p) K = -1 2 4 （3）讨论状态反馈与输出反馈的关系, 说明状态反馈为何能进行极点配置?
在经典控制理论中,一般只考虑由系统的输出变量来构成反馈律，即输出反馈。在现代控制理论的状态空间分析方法中,多考虑采用状态变量来构成反馈律，即状态反馈。从状态空间模型输出方程可以看出，输出反馈可视为状态反馈的一个特例。状态反馈可以提供更多的补偿信息，只要状态进行简单的计算再反馈，就可以获得优良的控制性能。 (4)使用状态反馈配置极点的前提是系统的状态是完全可控的。 2.已知系统 设计全维状态观测器，使观测器的极点配置在12+j，12-j 。 （1）给出原系统的状态曲线。 （2）给出观测器的状态曲线并加以对比。（观测器的初始状态可以任意选取）观察实验结果，思考以下问题： （1）说明反馈控制闭环期望极点和观测器极点的选取原则。 （2）说明观测器的引入对系统性能的影响。 （1）A=[0 1;-3 -4]; B=[0;1]; C=[2 0]; D=[]; G=ss(A,B,C,D); x=0:0.001:5; U=0*(x<0)+1*(x>0)+1*(x==0); X0=[0 1]'; T=0:0.001:5; lsim(G,U,T,X0);
7状态空间设计法极点配置观测器解析
第7章线性定常离散时间状态空间设计法 7.1引言 7.2状态反馈配置极点 7.3状态估值和状态观测器 7.4利用状态估值构成状态反馈以配置极点 7.5扰动调节 7.6无差调节
7.1 引言 一个被控对象： (1)()()()() ():1,():1,:,:,:x k Fx k Gu k y k Cx k x k n u k m F n n G n m C r n +=+?? =?????? 7.1 当设计控制器对其控制时，需要考虑如下各因素： ● 扰动，比如负载扰动 ● 测量噪声 ● 给定输入的指令信号 ● 输出 如图7.1所示。 给d L (k )扰动 图7.1 控制系统示意图 根据工程背景的不同，控制问题可分为调节问题和跟踪问题，跟踪问题也称为伺服问题。 调节问题的设计目标是使输出迅速而平稳地运行于某一平衡状态。包括指令变化时的动态过程，和负载扰动下的动态过程。但是这二者往往是矛盾的，需要折衷考虑。 伺服问题的设计目标是对指令信号的快速动态跟踪。 本章研究基于离散时间状态空间模型的设计方法。 7.2研究通过状态变量的反馈对闭环系统的全部特征值任意配置——稳定性与快速线。 7.3考虑当被控对象模型的状态无法直接测量时，如何使用状态观测器对状态进行重构。 7.4讨论使用重构状态进行状态反馈时闭环系统的特征值。 7.5简单地讨论扰动调节问题。 7.6状态空间设计时的无差调节问题。
7.2 状态反馈配置极点 工程被控对象如式7.1，考虑状态反馈 ()()()u k v k Lx k =+ 7.2 如图7.2所示。式7.2带入式7.1，得 (1)()()()() ()()()x k Fx k Gu k y k Cx k u k v k Lx k +=+?? =??=+? 7.3 整理得 ()(1)()() ()()x k F GL x k Gv k y k Cx k +=++?? =? 7.4 (k ) v (k ) 图7.2 状态反馈任意配置闭环系统的极点 闭环系统的特征方程为 []det ()0zI F GL -+= 7.5 问题是在什么情况下式7.5的特征根是可以任意配置的？即任给工程上期望的n 个特征根λ1, λ2, ..., λn ，有 []1det ()()0n i i zI F GL z λ=-+=-=∏ 7.6 定理：状态反馈配置极点
控制器极点配置方法
控制器极点配置方法 如果已知系统的模型或传递函数，通过引入某种控制器，使得闭环系统的极点可以移动到指定的位置，从而使系统的动态性能得到改善。这种方法称为极点配置法。 例6－12 有一控制系统如图6－38，其中，要求设计一个控制器，使系统稳定。 图6-38 解：（1）校正前，闭环系统的极点： > 0 因而控制系统不稳定。 （2）在控制对象前串联一个一阶惯性环节，c>0，则闭环系统极点： 显然，当，时，系统可以稳定。但此对参数c 的选择依赖于 a 、b 。因而，可 选择控制器，c 、d ，则有特征方程： 当，时，系统稳定。 本例由于原开环系统不稳定，因而不能通过简单的零极点相消方式进行控制器的设计，其原因在于控制器的参数在具体实现中无法那么准确，从而可能导致校正后的系统仍不稳定。 例6－13 已知一单位反馈控制系统的开环传递函数：
要求设计一串联校正装置Gc(s) ，使校正后系统的静态速度误差系统，闭环主导极点在 处。 解：首先，通过校正前系统的根轨迹可以发现，如图6－39所示，其主导极点为： 。 图6-39 为使主导极点向左偏移，宜采用超前校正装置。 （2）令超前校正装置，可采用待定系数法确定相关参数： 又
其中、、、为待定系数。 进一步可得： 即 将代入式子可以得到：，，，。进一步可得超前校正装置的传递函数： 校正后系统的根轨迹如图6-39所示。 该校正装置与例6－7中由超前装置获取的校正装置结果基本相同，说明结果是正确的。 在matlab中，亦有相应的命令可进行极点配置，主要有三个算法可实现极点配置算法：Bass-Gura算法、Ackermann 算法和鲁棒极点配置算法。这些算法均以状态空间进行表征，通过设定期望极点位置，获取状态反馈矩阵K。下面通过示例介绍其中的一种算法。 例6－14 考虑给定的系统，其状态方程模型如下：
基于输出反馈的区域极点配置
第22卷第2期南 京 理 工 大 学 学 报Vol.22No.21998年4月 Journal of Nanjing University of Science and Technology Apr.1998 基于输出反馈的区域极点配置 X 王子栋X X 郭 治 (南京理工大学信息学院,南京210094)摘要 该文研究输出反馈情形下线性定常连续及离散系统区域极点配置的统一代数刻划问题,即利用完全参数化方法,设计输出反馈控制器,使闭环极点配置于指定圆形区域内。文中导出了期望输出反馈控制器存在的充要条件,并进一步给出了这类控制器的全部参数化刻划。最后,得到了若干有益的推论,包括线性离散及连续系统稳定化控制器的统一代数表示等。 关键词 线性系统,输出反馈,极点配置,参数法,代数刻划 分类号 TP 202.1,T P 214.1 众所周知,线性定常系统的稳态及动态特性直接受其极点所在位置的影响,因而极点配置问题一直是控制理论研究中基本而重要的课题之一,其在工程实践中也具有明显的应用背景,如飞行控制系统的设计以及柔性结构的振动控制等[1]。迄今为止,精确极点的配置问题已得到了很好的研究。在过去的十年中,区域极点的配置问题也开始受到充分的注意,涌现出一批成果[2][3]。 目前,区域极点配置的相关文献中的大部分均是针对某性能指标给出具体的设计方法,且均集中于状态反馈情形,缺乏一定的通用性。本文对连续及离散线性定常系统使用统一的代数方法,给出了配置闭环极点至给定圆形区域的输出反馈控制器的全部参数化刻划,为区域极点配置问题提供了一条具有理论意义及应用价值的新途径。 1 问题的描述 考虑线性定常连续系统x a (t )=A x (t )+B u (t ),y (t )=Cx (t )及线性定常离散系统x (k +1)=A x (k )+Bu (k ),y (k )=Cx (k ),其中x ∈R n 为状态,u ∈R m 为控制输入,y ∈R p 为测量输出,A 、B 、C 为适维已知常数阵。(A ,B )及(A ,C )分别为可控和可观的。 考虑圆形区域D (A ,r ),其中在连续时间情形D (A ,r )表示圆心在A +j 0(A <0)处、半径为r (r <-A )的圆,在离散时间情形D (A ,r )表示单位圆内圆心位于A +j 0、半径为r 的圆。这里均考虑复平面。 X X XX 王子栋 男 32岁 副教授 国家自然科学基金及高校博士学科点专项科研基金资助项目 本文于1997年1月14日收到
单级倒立摆系统的极点配置与状态观测器设计
单级倒立摆系统的极点配置与状态观测器设计 14122156 杨郁佳 （1）倒立摆的运动方程并将其线性化 选取小车的位移z ，及其速度z g 、摆的角位置θ及其角速度θg 作为状态变量，即T x z z θθ??=??? ?g g 则系统的状态空间模型为 01000100000010()1000mg M M x u M m g Ml Ml x ????????????-????=+????????+-????????????g []1000y x = 设M=2kg ，m=0.2kg ，g=9.81m/2 s ，则单级倒立摆系统的状态方程为 (1010) 01010 01020.500013030 011040.54x x x x u x x x x ??????????????????-????????=+????????????????-???????????? []12100034x x y x x ???? ??=?????? （2）状态反馈系统的极点配置。 首先，使用MATLAB ，判断系统的能控性矩阵是否为满秩。 MATLAB 程序如下：
A=[0 1 0 0; 0 0 -1 0; 0 0 0 1; 0 0 11 0]; B=[0; 0.5; 0; -0.5]; C=[1 0 0 0]; D=0; rct=rank(ctrb(A,B)) [z,p,k]=ss2zp(A,B,C,D) MATLAB程序执行结果如下： 系统能控，系统的极点为 1=0 λ 2=0 λ 3=3.3166 λ 4=-3.3166 λ 可以通过状态反馈来任意配置极点，将极点配置在 1=-3 λ* 2=-4 λ* 3=-5 λ* 4=-6 λ*
基于极点配置方法的直流电机转速控制系统设计
摘要 建模、控制与优化是控制理论要解决的主要问题。在这些问题中，广泛采用了现代数学方法，使得控制理论的研究不断深入，取得了丰硕的成果。建模是控制理论中所要解决的第一个问题。控制理论中的建模方法主要有两种，一是经验建模，二是根据物理规律建模。所研究的对象主要是动态模型，一般用微分方程或差分方程来描述。设计控制系统是控制理论的核心内容。在线性系统中，我们所用到的数学工具是拓扑、线性群。在非线性系统中，我们用到了微分几何。可以说微分几何是非线性控制理论的数学基础。优化是控制的一个基本目的，而最优控制则是现代控制理论的一个重要组成部分。例如庞特里亚金的极大值原理、贝尔曼的动态规划，都是关于优化和最优控制问题的。 本报告首先介绍了直流电动机的物理模型, 并测量计算了它的具体参数。然后根据牛顿第二定律和回路电压法分别列写运动平衡方程式和电机电枢回路方程式，从而通过一些数学变换抽象出了以电压为输入、转速为输出、电流和转速为状态变量的数学模型。通过对抽象出来的模型进行性能分析，确定需要使用状态观测器来修正系统。继而借助MATLAB软件对转速环进行了状态反馈控制器的设计，使系统的阶跃响应达到了设计指标。 关键词：建模控制理论设计控制系统直流电动机转速状态反馈控制器
1 系统的物理模型、参数及设计要求 -------------------- 4 1.1 系统模型 ------------------------------------- 4 1.2 系统参数 ------------------------------------- 5 1.3 设计要求 ------------------------------------- 5 2 系统模型的建立------------------------------------ 6 2.1 模型抽象 ------------------------------------- 6 2.2 所建模型的性能分析 --------------------------- 7 3 系统状态观测器的设计----------------------------- 11 3.1 期望配置的极点的确定以及状态观测器的设计----- 11 3.1.1 第一组极点配置-------------------------- 11 3.1.2 第二组极点配置-------------------------- 11 3.2 状态观测器的设计 ---------------------------- 12 3.2.1 第一组极点------------------------------ 12 3.2.2 第二组极点------------------------------ 14 3.3 状态观测器的仿真图 -------------------------- 16 3.4 原系统加了状态观测器后的仿真结果图及分析----- 17 3.4.1 第一组极点------------------------------ 17 3.4.2 第二组极点------------------------------ 18 4 状态观测器极点配置与PID方法的比较 --------------- 20 4.1 直流电机转速、电流PID控制的设计------------- 20 4.2 两种方法的比较 ------------------------------ 21
控制系统的极点配置设计法
控制系统的极点配置设计法 一、极点配置原理 1.性能指标要求 2.极点选择区域 主导极点： 2 11 1 cos tan ξ βξ ξ -- - == 图3.22 系统在S平面上满足 时域性能指标的范围 n s t ζω 4 = ；当Δ=0.02时，。 n s t ζω 3 = 当Δ=0.05时,
3.其它极点配置原则 系统传递函数极点在s 平面上的分布如图（a ）所示。极点s 3距虚轴距离不小于共轭复数极点s 1、s 2距虚轴距离的5倍，即n s s ξω5Re 5Re 13=≥（此处ξ，n ω对应于极点s 1、s 2） ；同时，极点s 1、s 2的附近不存在系统的零点。由以上条件可算出与极点s 3所对应的过渡过程分量的调整时间为 135 1 451s n s t t =?≤ ξω 式中1s t 是极点s 1、s 2所对应过渡过程的调整时间。 图（b ）表示图（a ）所示的单位阶跃响应函数的分量。由图可知，由共轭复数极点s 1、s 2确定的分量在该系统的单位阶跃响应函数中起主导作用，即主导极点。因为它衰减得最慢。其它远离虚轴的极点s 3、s 4、s 5 所对应的单位阶跃响应衰减较快，它们仅在极短时间内产生一定的影响。因此，对系统过渡过程进行近似分析时。可以忽略这些分量对系统过渡过程的影响。 n x o (t) （a ） （b ） 系统极点的位置与阶跃响应的关系
二、极点配置实例 磁悬浮轴承控制系统设计 1.1磁悬浮轴承系统工作原理 图1是一个主动控制的磁悬浮轴承系统原理图。主要由被悬浮转子、传感器、控制器和执行器(包括电磁铁和功率放大器)四大部分组成。设电磁铁绕组上的电流为I0，它对转子产生的吸力F和转子的重力mg相平衡，转子处于悬浮的平衡位置，这个位置称为参考位置。 （a）（b） 图1 磁悬浮轴承系统的工作原理 Fig.1 The magnetic suspension bearing system principle drawing 假设在参考位置上，转子受到一个向下的扰动，转子就会偏离其参考位置向下运动，此时传感器检测出转子偏离其参考位置的位移，控制器将这一位移信号变换成控制信号，功率放大器又将该控制信号变换成控制电流I0+i，控制电流由I0增加到I0+i，因此，电磁铁的吸力变大了，从而驱动转子返回到原来的平衡位置。反之，当转子受到一个向上的扰动并向上运动，此时控制器使得功率放大器的输出电流由I0，减小到I0-i，电磁铁的吸力变小了，转子也能返回到原来的平衡位置。因此，不论转子受到向上或向下的扰动，都能回到平衡状态。这就是主动磁轴承系统的工作原理。即传感器检测出转子偏移参考点的位移，作为控制器的微处理器将检测到的位移信号变换成控制信号，然后功率放大器将这一控制信号转换成控制电流，控制电流在执行磁铁中产生磁力从而使转子维持其悬浮位置不变。悬浮系统的刚
倒立摆系统的状态空间极点配置控制设计
摘要：为实现多输入、多输出、高度非线不稳定的倒立摆系统平衡稳定控制，将倒立摆系统的非线性模型进行近似线性化处理，获得系统在平衡点附近的线性化模型。利用牛顿—欧拉方法建立直线型一级倒立摆系统的数学模型。在分析的基础上，基于状态反馈控制中极点配置法对直线型倒立摆系统设计控制器。由MATLAB仿真表明采用的控制策略是有效的，设计的控制器对直线型一级倒立摆系统的平衡稳定性效果好，提高了系统的干扰能力。 关键词：倒立摆、极点配置、MATLAB仿真 引言：倒立摆是进行控制理论研究的典型试验平台，由于倒立摆本身所具有的高阶次、不稳定、非线性和强耦合性，许多现代控制理论的研究人员一直将他视为典型的研究对象，不断从中发掘出新的控制策略和控制方法。控制器的设计是倒立摆系统的核心内容，因为倒立摆是一个绝对不稳定的系统，为使其保持稳定并且可以承受一定的干扰，基于极点配置法给直线型一级倒立摆系统设计控制器 1．数学模型的建立 倒立摆系统其本身是自不稳定的系统，实验建模存在着一定的困难。在忽略掉一些次要的因素之后，倒立摆系统就是一典型的运动的刚体系统，可以在惯性坐标系中应用经典力学理论建立系统动力学方程。下面采用牛顿-欧拉方法建立直线型一级倒立摆系统的数学模型。 1.1微分方程的数学模型 在忽略了空气阻力和各种摩擦力之后，可将直线一级倒立摆系统抽象成小车和匀质杆组成的系统，如图1所示：
图1：直线一级倒立摆模型 设系统的相关参数定义如下： M：小车质量 m：摆杆质量 b：小车摩擦系数 l：摆杆转动轴心到杆质心的长度 I：摆杆质量 F：加在小车上的力 x：小车位置 Φ：摆杆与垂直方向上方向的夹角 θ：摆杆与垂直方向下方向的夹角（摆杆的初始位置为竖直向下） 如下图2所示为小车和摆杆的受力分析图。其中，N和P为小车与摆杆相互作用力的水平和垂直方向的分量。
基于MATLAB的状态观测器设计
基于MATLAB 的状态观测器设计 预备知识： 极点配置 基于状态反馈的极点配置法就是通过状态反馈将系统的闭环极点配置到期望的极点位置上，从而使系统特性满足要求。 1. 极点配置原理 假设原系统的状态空间模型为： ???=+=Cx y Bu Ax x 若系统是完全可控的，则可引入状态反馈调节器，且： Kx u input -= 这时，闭环系统的状态空间模型为： ???=+-=Cx y Bu x )BK A (x 2. 极点配置的MATLAB 函数 在MATLAB 控制工具箱中，直接用于系统极点配置的函数有acker()和place()。调用格式为： K=acker(A,C,P) 用于单输入单输出系统 其中：A ，B 为系统矩阵，P 为期望极点向量，K 为反馈增益向量。 K=place(A,B,P) (K,prec,message)=place(A,B,P) place()用于单输入或多输入系统。Prec 为实际极点偏离期望极点位置的误差；message 是当系统某一非零极点偏离期望位置大于10%时给出的警告信息。 3. 极点配置步骤： （1）获得系统闭环的状态空间方程； （2）根据系统性能要求，确定系统期望极点分布P ； （3）利用MATLAB 极点配置设计函数求取系统反馈增益K ； （4）检验系统性能。 已知系统模型 如何从系统的输入输出数据得到系统状态？
初始状态：由能观性，从输入输出数据确定。 不足：初始状态不精确，模型不确定。 思路：构造一个系统，输出逼近系统状态 称为是的重构状态或状态估计值。实现系统状态重构的系统称为状态观 测器。 观测器设计 状态估计的开环处理： 但是存在模型不确定性和扰动！初始状态未知！ 应用反馈校正思想来实现状态重构。 通过误差来校正系统：状态误差，输出误差。 基于观测器的控制器设计 系统模型 若系统状态不能直接测量， 可以用观测器来估计系统的状态。 L是观测器增益矩阵，对偏差的加权。 真实状态和估计状态的误差向量 误差的动态行为：
控制系统的极点配置实验报告
课程名称： 控制理论乙 指导老师： 姚唯 成绩： 实验名称： 控制系统的极点配置 实验类型： 同组学生姓名： 郁明非 一、实验目的和要求（必填） 二、实验内容和原理（必填） 三、主要仪器设备（必填） 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析（必填） 七、讨论、心得 一、实验目的和要求 1．掌握全状态反馈系统的极点配置方法 2．在Simulink 仿真环境中，研究极点配置对系统特性的影响 二、实验内容和原理 （一）实验内容 1．一被控对象，其传递函数为 ) 3)(2)(1(10 )(+++= s s s s G 设计反馈控制器u=-kx ，使闭环系统的极点为3221j +-=μ，3222j --=μ，103-=μ。 2． 在Simulink 仿真环境下，用基本环节组成经过极点配置后的系统，通过图形观察环节，观察系统的各点响应。 （二）实验原理 对一给定控制系统如果其状态完全可控，则可进行任意极点配置即通过设计反馈増益K 使闭环系统具有期望的极点。极点配置有二种方法：第一种方法是采用变换矩阵T ，使系统具有期望的极点，从而求出矩阵K ；第二种方法基于Caylay-Hamilton 理论，通过矩阵多项式φ（a），可求出K （这种方法称为Ackermann 公式）。在MATLAB 中，利用控制系统工具箱函数place 和acker 进行极点配置设计。 三、主要仪器设备 一台PC 电脑，matlab 仿真软件，simulink 仿真环境 四、实验源代码及实验结果
function jidianpeizhi num=[10]; den=[1,6,11,6]; [A,B,C,D]=tf2ss(num,den); J=[-2-j*2*sqrt(3),-2+j*2*sqrt(3),-10]; K=place(A,B,J); K sys=ss(A-B*K,[0;0;0],eye(3),0); t=0::4; X=initial(sys,[1;0;0],t); x1=[1,0,0]*X'; x2=[0,1,0]*X'; x3=[0,0,1]*X'; subplot(3,1,1); plot(t,x2); grid on; title('Reponse to initial condition'); ylabel('x1'); subplot(3,1,2); plot(t,x2); grid on; ylabel('x2'); subplot(3,1,3); plot(t,x3); grid on; ylabel('x3'); xlabel('t(sec)');
状态反馈与极点配置报告
自 动 控 制 原 理 （课程设计）
一、题目 用MATLAB创建用户界面，并完成以下功能： （1）由用户输入被控系统的状态空间模型、闭环系统希望的一组极点； （2）显示未综合系统的单位阶跃响应曲线； （3）显示采用一般设计方法得到的状态反馈矩阵参数； （4）显示闭环反馈系统的单位阶跃响应曲线； （5）将该子系统嵌入到寒假作业中程序中。 分别对固定阶次和任意阶次的被控系统进行设计。分别给出设计实例。 二、运行结果 界面：如图 由用户输入被控系统的状态空间模型、闭环系统希望的一组极点 例如，输入 010 001 034 A ?? ?? =?? ?? -- ?? ， 1 B ?? ?? =?? ?? ?? ，[] 2000 C=，0 D=，闭环系统 希望的一组极点:22j -+、22j --、5 -如图所示：
被控系统的单位阶跃响应曲线 闭环系统的单位阶跃响应曲线
状态反馈矩阵显示 三、讨论 该闭环控制系统的状态反馈与极点配置设计系统可用于任意阶次的控制系统。在此之前，我还做了一个固定阶次的控制系统状态反馈与极点配置的Matlab 控制台程序（见附录二）。 该系统的利用状态反馈进行极点任意配置所采用的方法为一般方法，其步骤如下： ①判断受控系统是否完全能控； ②由给定的闭环极点要求确定希望的闭环特征多项式的n个系数 ~ i a； ③确定原受控系统的特征多项式系数i a； ④确定系统状态反馈矩阵 ~ ~~ ~ [,,,] 12n f f f F=的诸元素~~1 1i i i f a a - =- -； ⑤确定原受控系统化为能控标准形的变换阵的逆1 P-， ⑥确定受控系统完成闭环极点配置任务的状态反馈阵 ~ 1 F F P-=。 四、参考文献 [1]黄家英.《自动控制原理》.高等教育出版社,2010.5 [2]唐向红，郑雪峰.《MATLAB及在电子信息类》.电子工业出版社,2009.6 [3]吴大正，高西全.《MATLAB新编教程》.机械工业出版社，2008.4 五、附录 function varargout = tufeiqiang(varargin) %TUFEIQIANG M-file for tufeiqiang.fig % TUFEIQIANG, by itself, creates a new TUFEIQIANG or raises the existing % singleton*. % % H = TUFEIQIANG returns the handle to a new TUFEIQIANG or the handle to % the existing singleton*. % % TUFEIQIANG('Property','Value',...) creates a new TUFEIQIANG using
极点配置自适应控制器的设计
具有抑制噪声性质的极点配置自校正控制器设计 摘要：本文首先介绍了极点配置自校正控制器设计中的一些不足之处，然后提 出了一种极点配置自校正控制算法.它具有以下三个特点:（1）除了配置系统的闭环极点外，还可抑制噪声对系统输出的干扰；（2）与其他极点配置自校正控制算法相比，本算法所需的计算量大为减少；（3）理论上可得到算法的收敛性结果.该算法已在实际系统中得到了成功的应用。 关键字：极点配置 自校正 仿真 一、前言 极点配置自校正控制器在最近几年里得到了很大的发展，出现了许多自校正算法，但现有的各种自校正控制算法还存在一些不够完善之处。首先，由于这种控制器的控制目的是配置系统的闭环极点，所以，在控制器的设计过程中，对如何减少噪声对系统输出的干扰这个问题没有特别加以考虑；第二，计算量比较大，难以实时控制；第三，理论上不够完善，如现有的各种算法都没有收敛性证明。所以说，极点配置自校正控制器的理论和应用都不太成熟，有待于进一步发展。 二、算法 设所研究的系统由下面方程描述。 (1) )e (t )C (q )u (t )(q q (t))(-1-1-k 1-+=B y q A 其中(t)y 为输出量，(t)u 为控制量，(t)e 为零均值白噪声且与t 时刻以前（不包括 t 时刻）的控制(t) u 无关。 c c b b a a n n n n n n q c q c q C q b q b b q B q a q a A -1 -11--1 -101 --1 -11 -1)()(1)q (+???++=+???++=+???++= )(1 -q A ，)(1 -q C 为稳定多项式，0 0≠b ，1-q 为延后一步算子。 对上述系统如采用如下控制律 ) () () ()() F(q ) G(q -(t)1 1111 --1t y q F q G t y u r --+ = （2) 其中(t)r y 为参考输入，)(1-q G ，)(1-q F ，)(1-1q G ，)(q -11F 为四个适当次数的多项。则可得闭环方程 (4) (t )u (t )u u (t )(3) (t )(t )y (t )2121+=+=y y
鲁棒控制简介
当今的自动控制技术都是基于反馈的思想。反馈理论的要素包括三个部分：测量、比较和执行。测量关心的变量，与期望值相比较，用这个误差纠正调节控制系统的响应。 这个理论应用于自动控制的关键是，做出正确的测量和比较后，如何利用误差才能更好地纠正系统（即控制器的设计）。 鲁棒控制（Robust Control）方面的研究始于20世纪50年代。在过去的20年中，鲁棒控制一直是国际自控界的研究热点。所谓“鲁棒性”，是指控制系统在一定（结构，大小）的参数摄动下，维持某些性能的特性。根据对性能的不同定义，可分为稳定鲁棒性和性能鲁棒性。以闭环系统的鲁棒性作为目标设计得到的固定控制器称为鲁棒控制器。 由于工作状况变动、外部干扰以及建模误差的缘故，实际工业过程的精确模型很难得到，而系统的各种故障也将导致模型的不确定性，因此可以说模型的不确定性在控制系统中广泛存在。如何设计一个固定的控制器，使具有不确定性的对象满足控制品质，也就是鲁棒控制，成为国内外科研人员的研究课题。 主要的鲁棒控制理论有：（1）Kharitonov区间理论；（2）H∞控制理论（IMPORTANT）；（3）结构奇异值理论（μ理论）等等。 H∞控制理论 H∞控制理论是20世纪80年代开始兴起的一门新的现代控制理论。H∞控制理论是为了改变近代控制理论过于数学化的倾向以适应工程实际的需要而诞生的,其设计思想的真髓是对系统的频域特性进行整形(Loopshaping),而这种通过调整系统频率域特性来获得预期特性的方法,正是工程技术人员所熟悉的技术手段,也是经典控制理论的根本。 1981年Zames首次用明确的数学语言描述了H∞优化控制理论,他提出用传递函数阵的H∞范数来记述优化指标。1984年加拿大学者Fracis和Zames用古典的函数插值理论提出了H∞设计问题的最初解法,同时基于算子理论等现代数学工具,这种解法很快被推广到一般的多变量系统,而英国学者Glover则将H∞设计问题归纳为函数逼近问题,并用Hankel算子理论给出这个问题的解析解。Glover 的解法被Doyle在状态空间上进行了整理并归纳为H∞控制问题,至此H∞控制理论体系已初步形成。 在这一阶段提出了H∞设计问题的解法,所用的数学工具非常繁琐,并不像问题本身那样具有明确的工程意义。直到1988年Doyle等人在全美控制年会上发表了著名的DGKF论文,证明了H∞设计问题的解可以通过适当的代数Riccati方程得到。DGKF的论文标志着H∞控制理论的成熟。迄今为止,H∞设计方法主要是DGKF等人的解法。不仅如此,这些设计理论的开发者还同美国的The Math Works公司合作,开发了MA TLAB中鲁棒控制软件工具箱(Robust Control Toolbox),使H∞控制理论真正成为实用的工程设计理论。 研究——现代鲁棒控制（有界摄动） 鲁棒控制的早期研究，主要针对单变量系统（SISO）的在微小摄动下的不确定性，具有代表性的是Zames提出的微分灵敏度分析。然而，实际工业过程中故障导致系统中参数的变化，这种变化是有界摄动而不是无穷小摄动。因此产生了以讨论参数在有界摄动下系统性能保持和控制为内容的现代鲁棒控制。 现代鲁棒控制是一个着重控制算法可靠性研究的控制器设计方法。其设计目
系统稳定性分析 、利用MATLAB 实现极点配置、设计状态观测器
实验报告 实验名称系统稳定性分析、利用MATLAB实现极点配置、设计状态观测器系专业班 姓名学号授课老师 预定时间实验时间实验台号 一、目的要求 掌握系统稳定性的概念。学会使用MATLAB确定线性定常系统和非线性定常系统的稳定性。 掌握状态反馈和输出反馈的概念及性质。 掌握利用状态反馈进行极点配置的方法。学会用MATLAB求解状态反馈矩阵。 掌握状态观测器的设计方法。学会用MATLAB设计状态观测器。 熟悉分离定理，学会设计带有状态观测器的状态反馈系统。 二、原理简述 函数eig()的调用格式为V=eig(A)返回方阵A的特征值。 函数roots()的调用格式为roots(den)，其中den为多项式的系数行向量。计算多项式方程的解。 函数pole()的调用格式为pole(G)，其中G为系统的LTI对象。计算系统传递函数的极点。 函数zpkdata()的调用格式为[z,p,k]=zpkdata(G,’v’)，其中G为系统LTI对象。返回系统的零点、极点和增益。 函数pzmap()的调用格式为pzmap(G)，其中G为LTI对象。绘制系统的零点和极点。 对于线性定常连续系统x Ax，若A是非奇异矩阵，则原点是其唯一的平衡状态。统在原点处大范围渐近稳定的充分条件是：存在李氏函数v(x)x T px，且v(x)正定，v(x)负定。 如果SISO线性定常系统完全能控，则可通过适当的状态反馈,将闭环系统极点配置到 任意期望的位置。 MATLAB提供的函数acker()是用Ackermann公式求解状态反馈阵K。 MATLAB提供的函数place()也可求出状态反馈阵K。 如果线性定常系统完全能观测，则可构造全维（基本）观测器。全维（基本） 状态观测器的状态方程为观测器的反馈矩阵L为 其中为系统的能观测矩阵。 其中为期望的状态观测器的极点。观测器设计是极点配置的对偶问题，故可利用函数acker()和place()进行求解。
基于极点配置的横滚自动驾驶仪设计
基于极点配置的横滚自动驾驶仪设计 文章基于尾舵控制的弹体，采用极点配置设计方法，建立了控制参数与开环穿越频率，阻尼比及零极点比值的解析关系。在综合考虑执行机构带宽及相位裕度的约束的前提下，完成横滚自动驾驶仪设计。仿真结果表明系统时域响应品质良好，具备良好的稳定裕度，文章提出的设计方法有效。 标签：极点配置；横滚自动驾驶仪；穿越频率 同时四代机为实现隐身的目的，导弹需内埋弹射发射，在初始弹射段，导弹位于载机机腹下，载机流场复杂，面临严重的干扰力矩，这都对导弹的横滚控制的稳定性、快速性及无静差控制提出了严格的要求[1]。 文章针对尾舵控制导弹，设计了一种无静差快响应两回路横滚控制自动驾驶仪。利用简化的三阶横滚自动驾驶仪，采用极点配置方法，综合考虑执行机构带宽及稳定裕度约束的前提下，完成横滚自动驾驶仪的设计，实现横滚无静差控制，给出导弹典型空域的设计结果及线性化仿真结果。 1 横滚自动驾驶仪设计 忽略执行机构、速率陀螺及结构滤波器的动态特性，横滚自动驾驶仪为三阶系统，其简化的结构框图如图1所示： 自动驾驶仪待设计参数包含内回路角速度反馈增益Ig，PI控制器参数Ki，Kg。 从（6）式可以看出，通过理想极点配置方法，横滚控制增益大小与弹体气动特性c1、c2、系统开环截止频率？棕cr、实零点与实极点大小比值？浊及阻尼比？孜有关。c1、c2为弹体本体特性，决定于导弹的气动外形。通过相关公式可知？浊为略大于1的数值，对于横滚控制来说，阻尼？孜通常取为0.9左右，因此最终横滚控制增益的计算取决于系统开环截止频率的选取。 对于简化的三阶横滚自动驾驶仪开环穿越频率选取因遵循以下两个原则：（1）稳定性要求：相位裕度大于60°；（2）快速性要求：系统开环截止频率决定了横滚控制的快速性，工程上通常将开环截止频带取执行机构等效带宽的1/3左右。 2 设计仿真 取某典型飞行状态，选取期望的设计参数，计算相应的控制增益，完成横滚自动驾驶仪设计及线性化仿真。飞行状态包含飞行高度（H=10000）、飞行马赫数（M=3.5），设计约束包含舵机带宽（180 rad/s）、系统开环相位裕度（>60°），期望参数包含开环穿越频率？棕cr（180 rad/s）、阻尼比（？孜=0.9）比例值？浊