[發明專利]一種基于事件驅動的廣義預測自適應補給船航向控制方法有效
| 申請號: | 201610859143.7 | 申請日: | 2016-09-28 |
| 公開(公告)號: | CN106325071B | 公開(公告)日: | 2019-03-05 |
| 發明(設計)人: | 王元慧;王莎莎;付明玉;王成龍;庹玉龍;趙亮博;張放;張博;沈佳穎;葛然 | 申請(專利權)人: | 哈爾濱工程大學 |
| 主分類號: | G05B13/04 | 分類號: | G05B13/04;G05D1/02 |
| 代理公司: | 暫無信息 | 代理人: | 暫無信息 |
| 地址: | 150001 黑龍江省哈爾濱市南崗區*** | 國省代碼: | 黑龍江;23 |
| 權利要求書: | 查看更多 | 說明書: | 查看更多 |
| 摘要: | |||
| 搜索關鍵詞: | 補給船 廣義預測 航向控制 基于事件 預測模型 控制器 自適應 觸發 航向 判定 驅動觸發器 事件觸發器 自適應算法 最小二乘法 驅動 并行階段 參數變化 觸發時刻 觸發條件 低頻運動 舵角控制 函數判斷 航向跟蹤 滑動平均 環境影響 結合事件 離散事件 離散狀態 數學模型 遺忘因子 在線估計 觸發器 控制量 魯棒性 自回歸 遞推 舵角 解算 受控 補給 輸出 | ||
1.一種基于事件驅動的廣義預測自適應補給船航向控制方法,其特征在于:
第一,選用工作在離散狀態下的離散事件觸發器,設計變閾值下的觸發判定函數,在到達觸發時刻時,事件觸發器通過觸發判定函數判斷當前狀態是否滿足觸發條件;
定義時間段為t∈[tk,tk+1),時間觸發器的采樣周期為固定間隔T,事件下一個觸發時刻表示為:
其中:h為觸發周期,0<h<T,γ((h+1)T)為事件觸發器判斷觸發的閾值,為指數遞減函數,參數范圍ε>1,0≤α≤1,ε0≥0;e(hT)為hT時刻的觸發偏差;k∈N為整數
所述觸發判定函數Γ(s(ψ,e))為:
式中:ψ為船舶航向,ψdi為接收船的航向,ψi為補給船的航向,eADRC為控制器前一觸發時刻與當前時刻的偏差,e為觸發偏差,θ為eADRC調節因子,σ、μ為正態分布函數的參數,κ為閾值調節因子,s(ψ,e)表示在t≥0的情況下觸發方程s(t+k+1)=b(s(tk+1))的解,s=[ψΤ eΤ],且觸發條件是:滿足觸發判定函數的上下界之間的值;
第二,通過補給船低頻運動數學模型得出舵角-航向的受控自回歸積分滑動平均模型作為預測模型,采用遺忘因子遞推最小二乘法對預測模型的參數進行在線估計;
第三,結合事件驅動觸發器將需要控制的狀態發送給控制器,控制器通過GPC自適應算法解算后輸出舵角控制增量和控制量,實現補給船快速完成靠近階段和并行階段并保持補給階段的航向跟蹤。
2.根據權利要求1所述的一種基于事件驅動的廣義預測自適應補給船航向控制方法,其特征在于:
補給船低頻運動數學模型的狀態方程為:
式中:狀態向量X=[u v r]T分別為橫向速度、縱向速度和角速度,系數矩陣A(k)為時變矩陣,控制向量U=[Fx Fy Fn]T表示推力器在縱蕩、橫蕩和艏搖三個方向上產生的力和力矩,B為定常矩陣,得出:
式中:m為補給船的質量,Iz為轉動慣性矩,xG為補給船重心G點橫坐標值,其余均為水動力參數;u、v、r分別表示橫向速度、縱向速度、艏搖角速度;
將補給船的低頻運動數學模型離散化后得出:X(k+1)=G(k)X(k)+HU(k),且有
經過離散化和參數的在線估計,表述成航向角ψ與舵角δ的預測模型為:
其中:Δ=1-z-1是差分算子,ψ(k)、δ(k-1)和ξ(k)分別是系統輸出、系統輸入和均值為零、方差為σ2的白噪聲,z表示離散化后的變量算子,A1=-{I+G(k-2)[G(k-2)+I]-1[G(k-1)+I]},A2=G(k-2)[G(k-2)+I]-1[G(k-1)+I],參數表示A1,A2,B0,B1的估計值。
3.根據權利要求2所述的一種基于事件驅動的廣義預測自適應補給船航向控制方法,其特征在于:
由預測模型得到矩陣多項式:
并引入丟番圖方程:
得到預測輸出的向量形式為:
ψ=GΔδ+F(z-1)ψ(k)+H(z-1)Δδ(k-1)
式中:ψ=[ψ(k+1)T,...,ψ(k+N1)T]T,Δδ=[Δδ(k)Τ,...,Δδ(k+Nu-1)Τ]Τ,E′j,Fj,G′j,Hj為j時刻的多項式,
而性能指標函數為:其中:Δδ(k+j)=0,j=Nu,…,N1,Nu為控制時域,N1為最大預測時域,λ(j)是大于零的控制加權系數,ψr(k)有界為期望艏向設定值;
將預測輸出的向量帶入至性能指標函數的向量表達式中,得到未來控制時域內最優控制增量為:
ψ=(GΤG+λI)-1GΤ[ψr-F(z-1)ψ(k)-H(z-1)Δδ(k-1)]
定義(GΤG+λI)-1GΤ的前n行記為PΤ=[P1,…,PN1]、
則舵角控制增量為:Δδ(k)=PΤ[ψr-F(z-1)ψ(k)-H(z-1)Δδ(k-1)],
控制量為:δ(k)=δ(k-1)+Δδ(k)。
該專利技術資料僅供研究查看技術是否侵權等信息,商用須獲得專利權人授權。該專利全部權利屬于哈爾濱工程大學,未經哈爾濱工程大學許可,擅自商用是侵權行為。如果您想購買此專利、獲得商業授權和技術合作,請聯系【客服】
本文鏈接:http://www.szxzyx.cn/pat/books/201610859143.7/1.html,轉載請聲明來源鉆瓜專利網。





