发明名称 基于非线性H<sub>∞</sub>逆优化输出反馈控制器的船舶航向运动控制方法
摘要 本发明公开一种基于非线性H<sub>∞</sub>逆优化输出反馈控制器的船舶航向运动控制方法,通过设计船舶航向运动及海浪干扰滤波器,结合船载罗经实时采集的船舶航向信息,估计出船舶舵作用产生的航向运动信息,以及海浪扰动信息,再结合估计得到的航向运动信息利用基于船舶航向局部优化H<sub>∞</sub>控制器而设计的船舶航向运动H<sub>∞</sub>逆优化控制器,得到最优控制舵角,最终使船舶沿预定航向运动;并设计了控制效果监控器实时记录船舶的航向角与控制舵角,判断船舶是否能够跟踪船舶航向参考模型,若跟踪效果不好,则及时调整控制器。
申请公布号 CN103645736B 申请公布日期 2016.05.25
申请号 CN201310652115.4 申请日期 2013.12.05
申请人 哈尔滨工程大学 发明人 彭秀艳;赵新华;胡忠辉;王显峰
分类号 G05D1/02(2006.01)I 主分类号 G05D1/02(2006.01)I
代理机构 代理人
主权项 一种基于非线性H<sub>∞</sub>逆优化输出反馈控制器的船舶航向运动控制方法,其特征在于,包括如下步骤:步骤一:执行器接收到指令航向,开始进行船舶航向运动控制;步骤二:利用船载罗经实时测量舰船航向角,将航向角信息传输至船载计算机,利用船舶航向运动及海浪干扰滤波器估计出船舶舵作用产生的船舶航向角ψ<sub>L</sub>的估计值<img file="FDA0000932079710000011.GIF" wi="94" he="70" />船舶舵作用产生的船舶航向角速率r<sub>L</sub>的估计值<img file="FDA0000932079710000012.GIF" wi="71" he="70" />由海浪扰动产生的高频航向角ξ<sub>H</sub>的估计<img file="FDA0000932079710000013.GIF" wi="93" he="84" />海浪扰动产生的高频航向角速率ψ<sub>H</sub>的估计值<img file="FDA0000932079710000014.GIF" wi="101" he="71" />及海浪低频扰动量n的估计值<img file="FDA0000932079710000015.GIF" wi="64" he="62" />通过对高低频航向运动参数的辨识分离,从受海浪扰动的测量信号中重构了低频运动信号,实现海浪滤波;所涉及的船舶航向运动及海浪干扰滤波器表达式为<img file="FDA0000932079710000016.GIF" wi="1042" he="447" />式中,<img file="FDA0000932079710000017.GIF" wi="476" he="159" />ω<sub>n</sub>为遭遇频率,ζ<sub>n</sub>表示附加阻尼,且有ζ<sub>n</sub>∈[0.01,0.1];<img file="FDA0000932079710000018.GIF" wi="51" he="66" />为滤波误差,且有<img file="FDA0000932079710000019.GIF" wi="245" he="70" />ψ=ψ<sub>L</sub>+ψ<sub>H</sub>=ψ<sub>L</sub>+C<sub>w</sub>ξ;<img file="FDA00009320797100000110.GIF" wi="324" he="151" />其中<img file="FDA00009320797100000111.GIF" wi="445" he="111" /><img file="FDA00009320797100000112.GIF" wi="298" he="159" />n<sub>3</sub>、n<sub>1</sub>为常数,对航向不稳定的船舶有n<sub>1</sub><0,航向稳定n<sub>1</sub>>0;<img file="FDA00009320797100000113.GIF" wi="165" he="127" />u=δ,T为追随性时间指数,K为舵增益系数,δ为控制舵角;<img file="FDA00009320797100000114.GIF" wi="244" he="166" />k<sub>1</sub>∈R<sup>2</sup>,k<sub>2</sub>、k<sub>3</sub>、k<sub>4</sub>∈R为滤波增益系数,且满足k<sub>1</sub>=‑2(ζ<sub>i</sub>‑ζ<sub>n</sub>)ω<sub>c</sub>/ω<sub>n</sub>,k<sub>2</sub>=2(ζ<sub>i</sub>‑ζ<sub>n</sub>)ω<sub>n</sub>,k<sub>3</sub>=ω<sub>c</sub>,ζ<sub>i</sub>为陷波滤波阻尼、且ζ<sub>i</sub>>ζ<sub>n</sub>,ω<sub>c</sub>为截止频率、且ω<sub>c</sub>>ω<sub>n</sub>;T<sub>n</sub>为偏移时间常数,T<sub>n</sub>>>1且<img file="FDA00009320797100000115.GIF" wi="487" he="78" />步骤三:将步骤二中估计得到的<img file="FDA00009320797100000116.GIF" wi="191" he="71" />传输至基于船舶航向运动局部优化H<sub>∞</sub>控制器而设计的船舶航向运动H<sub>∞</sub>逆优化控制器,形成最优控制舵角指令,进而作用于舵机,产生实际舵角;控制效果监控部分实时记录船舶的航向角与控制舵角;所涉及的船舶航向运动局部优化H<sub>∞</sub>控制器为<img file="FDA0000932079710000021.GIF" wi="1133" he="86" />其中<img file="FDA0000932079710000022.GIF" wi="1333" he="239" />u=δ,<img file="FDA0000932079710000023.GIF" wi="301" he="230" /><img file="FDA0000932079710000024.GIF" wi="494" he="246" />ψ<sub>d</sub>为理想航向角;所涉及的船舶航向运动H<sub>∞</sub>逆优化控制器为u<sub>l</sub>=‑R<sup>‑1</sup>B<sup>T</sup>Px              (3)式中,u<sub>l</sub>为代价函数最优的控制器;R为第一加权系数矩阵,P=P<sup>T</sup>>0;所涉及的代价函数表达式为PA+A<sup>T</sup>P+P(γ<sup>‑2</sup>DD<sup>T</sup>‑BR<sup>‑1</sup>B<sup>T</sup>)P+Q=0    (4)式中,Q为第二加权系数矩阵,γ为给定的最优扰动抑制指标,且有γ≥γ<sup>*</sup>,γ<sup>*</sup>为非最优扰动抑制指标;步骤四:重复步骤二、步骤三,直至船舶按照指令航向运动。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室