CN103616817B - 一种基于非线性切换系统的伺服转台控制器 - Google Patents

Abstract

一种基于非线性切换系统的伺服转台控制器，其特征在于：将多模型切换系统引入伺服转台的精确建模和控制中，采用包括线性子系统和非线性子系统的全局在线滚动优化计算，对任意一个子控制器应用的预测函数控制器，对式所示的线性模型2，4直接应用局部预测函数控制器u₂(k)，u₄(k)；对其中的非线性模型1，3应用非线性预测函数控制器u₁(k)，u₃(k)。对其他型号具体转台还可根据非线性因素的具体情况，增加切换系统的模态和基于子系统模型的控制器，其操作方法与已表述部分相同。

Description

一种基于非线性切换系统的伺服转台控制器

技术领域

本发明涉及空间探测、卫星通信、国防工业等领域，尤其是一种基于非线性切换系统的伺服转台控制器。

背景技术

伺服转台在空间探测、卫星通信、国防工业等领域有着广泛的应用，是精密跟踪雷达、射电天文望远镜、惯导测试系统、火炮、云台等多种设备的重要组成装置，其技术水平的高低直接影响整个设备性能的优劣，是国家军事科技能力的体现点之一。对其展开的基于非线性切换系统的建模与控制器设计，具有重要意义和良好的应用发展前景，由此产生的技术方法和方案结论具有深刻的基础研究意义和应用价值。

过去对伺服转台的研究多停留在机理建模获得的线性系统模型，因为此类模型的局限性，应用的控制算法多为PID甚至是PI控制，无法保证精确的控制效果。目前，以伺服转台为基座和支撑的装备，朝大/微型两极化、数字化、精密化、低速平稳化等方向发展，在机械形变、低速爬行、摩擦、死区、齿轮空隙、风力扰动等方面呈现出更强的非线性因素和多种工作模态，也对系统控制策略提出了更高的要求。

伺服转台在不同的工作状态、环境和应用途径下，非线性因素的突出体现形式有所不同。在厘米波毫米波射电望远镜中，为克服场源遥远对灵敏度和抗干扰性的影响，天线口径、伺服转台的机械结构、驱动能力均往大型化发展，天线在过顶时和水平时伺服转台存在很大的机械结构形变，对系统性能产生不良影响；天线在换向、低速跟踪目标时，因粘滞摩擦、库伦摩擦、齿轮空隙、电机力矩波动等非线性因素的作用，使得整个系统的性能只有采用非线性模型才能准确描述，并需要在控制器中对其进行补偿；而大量早期工作和本课题组的科研结果表明：惯量适中的伺服转台在某些平稳运行的工况下，线性模型可以对系统的结构进行精确建模，并运用基于该模型的LQG控制器及H∞控制器。

切换系统是混杂动态系统的一种重要类型，可以准确描述许多实际模型；对连续系统采用切换控制器进行控制，可以获得比传统的反馈控制器更好的性能。传统控制在伺服转台的建模和控制中存在模型失配和粗犷设计、控制器结构过于简单、参数试凑等问题，为解决这些问题并应对复杂伺服系统在相关行业飞速发展中而出现的多工况多模态的情况，将多模型切换系统引入伺服转台的精确建模和控制中。建立伺服转台的非线性切换模型，设计切换Hammerstein系统预测控制器，在低速跟踪运行阶段、加速/减速运行阶段、过顶机械形变等情况下保持良好的跟踪效果。

目前现有的用于伺服转台控制器多为PID或PI控制器。常见天线伺服转台，控制目标为在风扰作用下对输入目标值进行位置跟踪，其伺服转台的负载为天线机构，均为刚性结构，控制器选用的是PI控制器。

PID控制器的可调参数较少，缺乏明确的工程意义，而且无法应对在复杂特殊工况下伺服转台的精确控制。现在存在一些对PID控制的改进，通常为分段参数的模糊PID或基于时间切换控制器的方案，均可在本专利的研究中，纳入切换控制器，对伺服转台的性能进一步提高和技术推广。

发明内容

针对上述现有技术的缺点，本发明的目的是提供基于非线性切换系统的伺服转台控制器。

一种基于非线性切换系统的伺服转台控制器，将多模型切换系统引入伺服转台的精确建模和控制中，采用包括线性子系统和非线性子系统的全局在线滚动优化计算，切换系统控制器施加于伺服转台的控制量，依切换律选择子控制器的输出。对线性子系统应用状态反馈控制器或线性预测控制器；对非线性子系统，构造基于Hammerstein模型和预测函数控制的非线性预测控制器，对实际输入解算方法进行改进。对子系统的控制采用的是预测控制，具体为基于约束Hammerstein模型的预测控制函数器。预测函数控制(PredictiveFunctional Control)的控制作用是若干已知基函数的组合，本专利中子控制器采用多项式基函数的形式，对一阶指数形式的参考轨迹进行滚动优化，通过最小化误差指标，求取中间变量。再采用改进的两步法非线性预测控制策略，将非线性函数和线性子系统都纳入目标函数，采用遗传算法(GA)的非线性规划法在线优化求解控制量。

附图说明

图1为伺服转台非线性切换控制器示意图；

图2为子系统非线性预测函数控制器结构示意图。

具体实施方式

下面结合附图对本发明作进一步的详细说明。

传统控制在伺服转台的建模和控制中存在模型失配和粗犷设计、参数试凑、控制器结构过于简单等问题，为解决这些问题并应对复杂伺服系统在相关行业飞速发展中而出现的多工况多模态的情况，将多模型切换系统引入伺服转台的精确建模和控制中，以获得良好的跟踪效果。

对伺服专台的连续系统采用切换控制器进行控制，可以获得比传统的PID控制器或其他连续控制其更好的性能。当转台以一定数值范围内的速度匀速运行，则控制器根据正转或反转的工作情况，基于相应的线性系统模型，应用线性预测控制器；当转台运行于超低速、急加速/急停等工况下，切换系统基于事件的条件将非线性控制器切入转台控制中，应用非线性约束Hammerstein模型预测控制器，以应对转台的非线性因素，保证复杂工况下的精确控制品质。

对子系统的控制采用的是预测控制，具体为基于约束Hammerstein模型的预测控制函数器。预测函数控制(Predictive Functional Control)的控制作用是若干已知基函数的组合，本专利中采用多项式基函数的形式，对一阶指数形式的参考轨迹进行滚动优化，通过最小化误差指标，求取中间变量。再采用改进的两步法非线性预测控制策略，将非线性函数和线性子系统都纳入目标函数，采用遗传算法(GA)的非线性规划法现在优化求解控制量。

预测函数控制器的控制输入公式可简单表示为：

$u\left(k\right)=k_0*\[c\left(k\right)-y_p\left(k\right)\]+\underset{j=1}{\overset{\max (N_c,N_c)}{\Σ}}{k}_j*\[c\left(k\right)-e_j\left(k\right)\]+k_m*X_m\left(k\right)---\left(4\right)$

其中参数k₀，k_j，k_m由结构参数和状态变量共同计算所得，c(k)为设定值，X_m(k)为状态值，y_p(k)和e_j(k)分别为预测输出和偏差信号，u(k)为当前控制量，具体为k₀为子控制器中设定值与预测输出之间偏差信号的权值系数，k_j为控制周期内局部偏差部分的权值系数，k_m为状态量的权值系数。以上参数可根据工程经验和仿真调试结果设定，因其为系数，因此没有量纲，以及公式中c(k)与y_p(k)分别为设定值和预测输出，具体物理意义与控制系统输出量有关，一般为伺服转台的输出转速(rad/s)或角位移(rad)。e_j(k)为某一模态下控制步长内第k步的偏差信号，X_m(k)为第k步状态输出，两者的物理意义和量纲定义同上。

对某非线性子系统控制器的基本原理图见图1，此方法原理简单、计算方便，但需要注意优化过程中高阶方程求解和根的分布是否符合约束等问题，所以目前常用的方法是针对奇数阶方程的实根求解或进行简化运算。本专利申请所采用的技术手段提出采用包括线性子系统和非线性子系统的全局在线滚动优化计算，解决非线性代数方程无实根或结算误差的问题。这样可以有效对内部结构类型为偶数次幂多项式的情况进行处理，同时解决了原两步法中不可避免的结算误差问题，提高了非线性因素较强的工况下，伺服转台控制器可以达到的精度。

对任意一个子控制器应用如上所述的预测函数控制器，对式所示的线性模型2，4直接应用局部预测函数控制器u₂(k)，u₄(k)；对其中的非线性模型1，3应用如图1所示的非线性预测函数控制器u₁(k)，u₃(k)，

$u_{\mathrm{\&sigma;}}\left(k\right)=\left\{\begin{array}{cccccc}{u}_1\left(k\right)& \mathrm{\&sigma;}=1,& s_1:& u\left(k\right)\&GreaterEqual;0\left|y\left(k\right)\right|<{\mathrm{\&epsiv;}}_1& or& \left|\mathrm{\&Delta;}y\left(k\right)\right|>{\mathrm{\&gamma;}}_1\\ u_2\left(k\right)& \mathrm{\&sigma;}=2,& s_2:& u\left(k\right)\&GreaterEqual;0\left|y\left(k\right)\right|\&GreaterEqual;{\mathrm{\&epsiv;}}_1& and& \left|\mathrm{\&Delta;}y\left(k\right)\right|\&le;{\mathrm{\&gamma;}}_1\\ u_3\left(k\right)& \mathrm{\&sigma;}=3,& s_3:& u\left(k\right)<0\left|y\left(k\right)\right|<{\mathrm{\&epsiv;}}_2& or& \left|\mathrm{\&Delta;}y\left(k\right)\right|>{\mathrm{\&gamma;}}_2\\ u_4\left(k\right)& \mathrm{\&sigma;}=4,& s_4:& u\left(k\right)>0\left|y\left(k\right)\right|\&GreaterEqual;{\mathrm{\&epsiv;}}_2& and& \left|\mathrm{\&Delta;}y\left(k\right)\right|\&le;{\mathrm{\&gamma;}}_2\end{array}\right.---\left(5\right)$

切换控制器结构如上式(5)。当伺服转台工作在匀速以正常范围内的速度旋转时，控制器为局部预测函数控制器；当伺服转台工作于绝对值低于某值的低速运转状态、突然加速或负载过顶机械形变大时，切换至基于Hammerstein非线性模型的预测控制器输出，达到特殊工况下的精确控制。在具体某转台的建模实施上可根据具体的非线性特点建立模态更多的切换系统模型及相应子系统控制算法；操作手法及辨识算法与前述情况相同。式中a_σ1，a_σ2，a_σ3和b_σ0，b_σ1为离散线性系统脉冲传递函数的系数，

即z^-1满足z^-1x(k)＝x(k-1)。

非线性子模块x(k)＝f_σ1u(k)+f_σ2u²(k)，σ＝1，2，3，4f_σ1，f_σ2为非线性多项式的系数，a_σ1，a_σ2，a_σ3，b_σ0，b_σ1，f_σ1，f_σ2均为建模的模型系数。首先对待控制的伺服转台实施系统辨识实验，记录现场数据后，应用本发明中提出的基于PSO的辨识算法获得。无量纲。式中s_i表示第i个模态，受状态信息驱动。式中|y(k)|＜ε₁表示伺服转台正向运动时，由线性子模型切换到非线性子模型的一种切换条件。ε₁为切换条件的阈值，对速度跟踪的伺服转台来说，为过零状态附近的一个较小的速度值，如正向旋转最大速度的1/10(rad/s)；ε₂为反向运动时，切换条件的阈值，意义同上。伺服转台正向运动和反向运动时摩擦作用方式不同，因此ε₁与ε₂大小可以不相同，具体取值应视伺服转台的实际情况而定，驱动机构及转台惯量等都对其数值有影响。|Δy(k)|＞γ₁为伺服转台正向运动时，由线性子模型切换到非线性子模型的另一种切换条件。γ₁，γ₂为正/反向运动时角加速度的阈值，单位为rad/s²。

上述实施例只为说明本发明的技术构思及特点，其目的在于让熟悉此领域技术的人士能够了解本发明的内容并加以实施，并不能以此限制本发明的保护范围，凡根据本发明精神实质所作的等效变化或修饰效果相似，都应涵盖在本发明的保护范围内。

Claims (2)

1.一种基于非线性切换系统的伺服转台控制器，其特征在于：将多模型切换系统引入伺服转台的精确建模和控制中，采用包括线性子系统和非线性子系统的全局在线滚动优化计算，对任意一个子控制器应用的预测函数控制器，对式所示的线性模型2，4直接应用局部预测函数控制器u₂(k)，u₄(k)，其中k₀为子控制器中设定值与预测输出之间偏差信号的权值系数，k_j为控制周期内局部偏差部分的权值系数，k_m为状态量的权值系数，c(k)与y_p(k)分别为设定值和预测输出，e_j(k)为某一模态下控制步长内第k步的偏差信号，X_m(k)为第k步状态输出，a_σ1，a_σ2，a_σ3和b_σ0，b_σ1为离散线性系统脉冲传递函数的系数，f_σ1，f_σ2为非线性多项式的系数，a_σ1，a_σ2，a_σ3，b_σ0，b_σ1，f_σ1，f_σ2均为建模的模型系数；对其中的非线性模型1，3应用非线性预测函数控制器u₁(k)，u₃(k)，

$u_{\mathrm{\&sigma;}}\left(k\right)=\left\{\begin{array}{cccccc}{u}_1\left(k\right)& \mathrm{\&sigma;}=1,& s_1:& u\left(k\right)\&GreaterEqual;0,\left|y\left(k\right)\right|<{\mathrm{\&epsiv;}}_1& or& \left|\mathrm{\&Delta;}y\left(k\right)\right|>{\mathrm{\&gamma;}}_1\\ u_2\left(k\right)& \mathrm{\&sigma;}=2,& s_2:& u\left(k\right)\&GreaterEqual;0,\left|y\left(k\right)\right|\&GreaterEqual;{\mathrm{\&epsiv;}}_1& and& \left|\mathrm{\&Delta;}y\left(k\right)\right|\&le;{\mathrm{\&gamma;}}_1\\ u_3\left(k\right)& \mathrm{\&sigma;}=3,& s_3:& u\left(k\right)<0,\left|y\left(k\right)\right|<{\mathrm{\&epsiv;}}_2& or& \left|\mathrm{\&Delta;}y\left(k\right)\right|>{\mathrm{\&gamma;}}_2\\ u_4\left(k\right)& \mathrm{\&sigma;}=4,& s_4:& u\left(k\right)>0,\left|y\left(k\right)\right|\&GreaterEqual;{\mathrm{\&epsiv;}}_2& and& \left|\mathrm{\&Delta;}y\left(k\right)\right|\&le;{\mathrm{\&gamma;}}_2\end{array}\right.$

切换控制器结构如上式，当伺服转台工作在匀速以正常范围内的速度旋转时，控制器为局部预测函数控制器；当伺服转台工作于绝对值低于某值得低速运转状态、突然加速或负载过顶机械形变大时，切换至基于Hammerstein非线性模型的预测控制器输出，达到特殊工况下的精确控制，式中s_i表示第i个模态，i为整数；|y(k)|＜ε₁表示伺服转台正向运动时，由线性子模型切换到非线性子模型的一种切换条件，ε₁为切换条件的阈值，ε₂为反向运动时，切换条件的阈值；|Δy(k)|＞γ₁为伺服转台正向运动时，由线性子模型切换到非线性子模型的另一种切换条件，γ₁，γ₂为正/反向运动时角加速度的阈值。

2.如权利要求1所述的基于非线性切换系统的伺服转台控制器，其特征在于：具体模型的个数和子控制器个数可适当增减。