发明名称 一种基于模型预测控制的平流层飞艇平面路径跟踪控制方法
摘要 一种基于模型预测控制的平流层飞艇平面路径跟踪控制方法,步骤如下:给定期望跟踪值;制导误差计算:计算期望位置与实际位置之间的距离误差,角度误差;动力学方程纵横向分解,控制器设计只取其横向状态量;求解离散化系统方程:对由以上步骤得到的平流层飞艇横侧向连续系统进行线性化处理,并且也将误差导数和进行线性化处理。然后将飞艇横向状态量和误差当成扩展状态量,并且对扩展连续状态空间方程离进行离散化处理;预测系统未来动态:根据由组合惯导等传感器测量得到的当前状态量预测未来某一段时间的状态量或输出量;构造模型预测控制目标函数:由预测状态量构造目标函数,并用标准QP算法进行求解得到系统输入量。
申请公布号 CN104317300B 申请公布日期 2017.04.05
申请号 CN201410486563.6 申请日期 2014.09.22
申请人 北京航空航天大学 发明人 祝明;余帅先;刘丽莎;闫柯瑜;陈天
分类号 G05D1/10(2006.01)I;G05B13/04(2006.01)I 主分类号 G05D1/10(2006.01)I
代理机构 代理人
主权项 一种基于模型预测控制的平流层飞艇平面路径跟踪控制方法,其特征在于具体步骤如下:步骤一给定期望跟踪值:给定期望平面路径;给定期望俯仰角θ<sub>c</sub>、期望滚转角φ<sub>c</sub>;给定期望速度υ<sub>c</sub>;步骤二制导误差计算:计算期望位置与实际位置之间的距离误差e,角度误差<img file="FDA0001160884340000011.GIF" wi="71" he="70" />所述的制导误差计算,其计算方法如下:1)计算期望路径的方向角Ψ<sub>p</sub>=arctan(a),方向角误差<img file="FDA0001160884340000012.GIF" wi="267" he="79" />其中Ψ为飞艇偏航角,a为路径参数直线斜率;2)计算飞艇当前位置与期望路径垂直距离误差<img file="FDA0001160884340000013.GIF" wi="315" he="159" />p=[x,y]<sup>T</sup>为平流层飞艇的当前位置,其中a路径参数直线斜率,b路径参数直线截距;3)计算期望偏航角为<img file="FDA0001160884340000014.GIF" wi="411" he="130" />Δ>0为控制参数;4)计算方向角误差导数<img file="FDA0001160884340000015.GIF" wi="147" he="79" />r为飞艇偏航角速度;5)计算距离误差导数<img file="FDA0001160884340000016.GIF" wi="477" he="95" />其中u为飞艇前向速度,v为飞艇侧向速度;步骤三动力学方程纵横向分解:将动力学方程按照纵横向进行分解,控制器设计只取其横向状态量;步骤四求解离散化系统方程:对由以上步骤得到的平流层飞艇横侧向连续系统进行线性化处理,并且也将误差导数<img file="FDA0001160884340000017.GIF" wi="30" he="71" />和<img file="FDA0001160884340000018.GIF" wi="50" he="78" />进行线性化处理;然后将飞艇横向状态量和误差当成扩展状态量,并且对扩展连续状态空间方程进行离散化处理,所求解离散化系统方程,其计算方法如下:1)线性化处理定义扩展状态变量<img file="FDA0001160884340000019.GIF" wi="341" he="87" />对扩展后系统方程在参考点<img file="FDA0001160884340000021.GIF" wi="386" he="87" />和u<sub>r</sub>=0进行线性化;线性化后得:<maths num="0001"><math><![CDATA[<mrow><mover><mover><mi>X</mi><mo>~</mo></mover><mo>&CenterDot;</mo></mover><mo>=</mo><msub><mo>&Integral;</mo><mrow><mi>x</mi><mo>,</mo><mi>r</mi></mrow></msub><mover><mi>X</mi><mo>~</mo></mover><mo>+</mo><msub><mo>&Integral;</mo><mrow><mi>u</mi><mo>,</mo><mi>r</mi></mrow></msub><mover><mi>u</mi><mo>~</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001160884340000022.GIF" wi="1398" he="127" /></maths>其中<img file="FDA0001160884340000023.GIF" wi="239" he="79" />和<img file="FDA0001160884340000024.GIF" wi="194" he="77" />为相对于参考点的误差;f<sub>x,r</sub>和f<sub>u,r</sub>为连续系统方程分别对X和u的偏导数;其中误差方程线性化后为<maths num="0002"><math><![CDATA[<mrow><mover><mi>e</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>U</mi><mi> </mi><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><msub><mi>&Psi;</mi><mrow><mi>L</mi><mi>O</mi><mi>S</mi></mrow></msub><mo>)</mo></mrow><mover><mi>&Psi;</mi><mo>&OverBar;</mo></mover><mo>+</mo><mi>v</mi><mi> </mi><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><msub><mi>&Psi;</mi><mrow><mi>L</mi><mi>O</mi><mi>S</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001160884340000025.GIF" wi="1430" he="86" /></maths><maths num="0003"><math><![CDATA[<mrow><mover><mover><mi>&Psi;</mi><mo>&OverBar;</mo></mover><mo>&CenterDot;</mo></mover><mo>=</mo><mi>r</mi></mrow>]]></math><img file="FDA0001160884340000026.GIF" wi="123" he="78" /></maths>2)计算离散系统方程由于采样时间一般都有比较短,我们可以采用前向差分的方法对连续系统进行离散化,设采样周期为T,线性化后的扩展状态空间方程为:<maths num="0004"><math><![CDATA[<mrow><mover><mover><mi>X</mi><mo>~</mo></mover><mo>&CenterDot;</mo></mover><mo>=</mo><mi>A</mi><mover><mi>X</mi><mo>~</mo></mover><mo>+</mo><mi>B</mi><mover><mi>u</mi><mo>~</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001160884340000027.GIF" wi="1438" he="103" /></maths>则离散化之后的离散系统状态空间方程为:<maths num="0005"><math><![CDATA[<mrow><msub><mover><mover><mi>X</mi><mo>~</mo></mover><mo>&CenterDot;</mo></mover><mrow><mi>K</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>A</mi><mi>K</mi></msub><msub><mover><mi>X</mi><mo>~</mo></mover><mi>K</mi></msub><mo>+</mo><msub><mi>B</mi><mi>K</mi></msub><msub><mover><mi>u</mi><mo>~</mo></mover><mi>K</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001160884340000028.GIF" wi="1430" he="118" /></maths>其中A<sub>K</sub>=A*T+I B<sub>K</sub>=B*T;步骤五预测系统未来动态:根据由组合惯导等传感器测量得到的当前状态量预测未来某一段时间的状态量或输出量;步骤六构造模型预测控制目标函数:由预测状态量构造目标函数,并用标准QP算法进行求解得到系统输入量u<sub>k</sub>。
地址 100191 北京市海淀区学院路37号