发明名称 一种基于触发式预测控制的飞行器控制方法
摘要 本发明公开一种基于触发式预测控制的飞行器控制方法,其特征在于,首先确定模型预测控制算法每一采样时刻的最优控制序列;然后通过系统的输出和状态变量确定条件触发机制的条件触发方程;最后根据控制需求设置系统的条件阈值,通过比较每一采样时刻条件触发方程的值与阈值之间的大小关系判断系统是否进行条件触发。该方法针对传统的模型预测控制算法在线计算量大的问题,将模型预测控制算法与条件触发机制相结合,以减小其在线计算量,提高模型预测控制算法的在线实施效率。
申请公布号 CN106569503A 申请公布日期 2017.04.19
申请号 CN201610943280.9 申请日期 2016.10.26
申请人 湖北工业大学 发明人 赵熙临;刘斐;明航;吴胧胧;何晶晶
分类号 G05D1/08(2006.01)I 主分类号 G05D1/08(2006.01)I
代理机构 武汉帅丞知识产权代理有限公司 42220 代理人 朱必武;曾祥斌
主权项 一种基于触发式预测控制的飞行器控制方法,其特征在于,首先确定模型预测控制算法每一采样时刻的最优控制序列;然后通过系统的输出和状态变量确定条件触发机制的条件触发方程;最后根据控制需求设置系统的条件阈值,通过比较每一采样时刻条件触发方程的值与阈值之间的大小关系判断系统是否进行条件触发,所述方法包含如下步骤:步骤1:模型预测控制系统的分析是基于研究对象的数学模型,假设以离散形式描述的研究对象数学模型如下所示:x(k+1)=Ax(k)+Bu(k)y=Cx(k)其中,x为状态变量,A为系统矩阵,B为输入矩阵,C为输出矩阵,k为采样时刻,目标函数定义为:<maths num="0001"><math><![CDATA[<mrow><mi>J</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mo>|</mo><mo>|</mo><mi>W</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>Y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><msubsup><mo>|</mo><msub><mi>Q</mi><mi>y</mi></msub><mn>2</mn></msubsup><mo>+</mo><mo>|</mo><mo>|</mo><mi>U</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><msubsup><mo>|</mo><msub><mi>R</mi><mi>y</mi></msub><mn>2</mn></msubsup></mrow>]]></math><img file="FDA0001140463800000011.GIF" wi="622" he="86" /></maths>其中,W(k)是输出的期望值,Q<sub>y</sub>、R<sub>y</sub>是输出和控制加权矩阵,其控制过程表现为针对目标函数的优化,在每一采样时刻系统都会计算出一个最优的输入序列:U(k)=[u(k);u(k+1)…u(k+m‑1)]其中m为系统的控制域,控制过程中,系统只将输入序列第一个元素u(k)作用于对象;步骤2:定义条件触发方程为:<maths num="0002"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><mi>y</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>,</mo><mover><mi>y</mi><mo>~</mo></mover><mo>(</mo><mi>k</mi><mo>)</mo><mo>,</mo><mi>x</mi><mo>(</mo><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo><mo>,</mo><mi>x</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>=</mo><mo>|</mo><mo>|</mo><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>y</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><mo>+</mo><mi>&lambda;</mi><mo>|</mo><mo>|</mo><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup></mrow>]]></math><img file="FDA0001140463800000012.GIF" wi="1086" he="88" /></maths>其中,y(k)表示输出,<img file="FDA0001140463800000013.GIF" wi="93" he="54" />表示预测值,x(k)为状态变量,||·||表示常规欧几里得规范,λ为相邻采样时刻状态变量的差值的权重;步骤3:设置阈值ε,比较ε与条件触发方程P的大小,当P的值大于阀值ε时,触发器被触发;否则,系统将前一时刻计算的输入序列中对应的值作为当前输入值作用于系统,系统的控制变量表示如下:<maths num="0003"><math><![CDATA[<mrow><mover><mi>u</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo></mrow></mtd><mtd><mrow><mi>w</mi><mi>h</mi><mi>e</mi><mi>n</mi></mrow></mtd><mtd><mrow><mi>P</mi><mrow><mo>(</mo><mi>y</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>,</mo><mover><mi>y</mi><mo>~</mo></mover><mo>(</mo><mi>k</mi><mo>)</mo><mo>,</mo><mi>x</mi><mo>(</mo><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo><mo>,</mo><mi>x</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>&gt;</mo><mi>&epsiv;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo></mrow></mtd><mtd><mrow><mi>w</mi><mi>h</mi><mi>e</mi><mi>n</mi></mrow></mtd><mtd><mrow><mi>P</mi><mrow><mo>(</mo><mi>y</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>,</mo><mover><mi>y</mi><mo>~</mo></mover><mo>(</mo><mi>k</mi><mo>)</mo><mo>,</mo><mi>x</mi><mo>(</mo><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo><mo>,</mo><mi>x</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>&le;</mo><mi>&epsiv;</mi></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001140463800000014.GIF" wi="1036" he="173" /></maths>其中,P(·)表示系统的触发方程,<img file="FDA0001140463800000021.GIF" wi="94" he="63" />为输入信号,u(k)为最优控制序列的第一个元素,u(k+j‑1)为输入序列的第j个元素,j表示系统连续不满足触发条件的次数。
地址 430068 湖北省武汉市洪山区南李路28号