[go: up one dir, main page]

CN115903802A - H∞ heading control method based on nonlinear variable parameter unmanned ship model - Google Patents

H∞ heading control method based on nonlinear variable parameter unmanned ship model Download PDF

Info

Publication number
CN115903802A
CN115903802A CN202211421529.1A CN202211421529A CN115903802A CN 115903802 A CN115903802 A CN 115903802A CN 202211421529 A CN202211421529 A CN 202211421529A CN 115903802 A CN115903802 A CN 115903802A
Authority
CN
China
Prior art keywords
heading
variable parameter
model
nonlinear
unmanned ship
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211421529.1A
Other languages
Chinese (zh)
Inventor
黄宴委
黄鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fuzhou University
Original Assignee
Fuzhou University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fuzhou University filed Critical Fuzhou University
Priority to CN202211421529.1A priority Critical patent/CN115903802A/en
Publication of CN115903802A publication Critical patent/CN115903802A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Feedback Control In General (AREA)

Abstract

The invention provides an H infinity heading control method based on a nonlinear variable parameter unmanned ship model, which considers the complex and changeable impact of water flow on a ship, obtains the circulation force borne by the ship by using a low aspect ratio wing theory and further establishes the nonlinear variable parameter heading control model. And (3) deducing a nonlinear matrix inequality condition meeting the H infinity robust stability of the system by a parameter-related Lyapunov function aiming at the established nonlinear variable parameter heading control model, and solving the nonlinear matrix inequality condition by an SOS (sequence of oriented services) tool box in MATLAB to obtain a nonlinear variable parameter H infinity robust controller.

Description

基于非线性变参数无人船模型的H∞艏向控制方法H∞ heading control method based on nonlinear variable parameter unmanned ship model

技术领域Technical Field

本发明属于无人船艏向控制技术领域,尤其涉及一种基于非线性变参数无人船模型的H∞艏向控制方法。The invention belongs to the technical field of unmanned ship heading control, and in particular relates to an H∞ heading control method based on a nonlinear variable parameter unmanned ship model.

背景技术Background Art

随着海洋与河域相关商业、军事的发展,无人船由于其具有体积小、速度高、成本低和无人员上网风险等优点,成为执行搜救、侦察、监测等任务的重要平台。艏向控制是无人船实现自主、可靠和安全航行的关键技术之一。然而,在无人船航行的过程中,无人船的状态往往是时变的,如速度时变,这使得无人船收到的水动力不断变化,破坏了系统的稳定结构。水动力变化主要体现在无人船模型参数上的不确定,这种现象导致无人船系统具有强时变非线性。With the development of ocean and river-related commerce and military, unmanned ships have become an important platform for performing search and rescue, reconnaissance, monitoring and other tasks due to their advantages such as small size, high speed, low cost and no risk of personnel being online. Bow control is one of the key technologies for unmanned ships to achieve autonomous, reliable and safe navigation. However, during the navigation of unmanned ships, the state of the unmanned ships is often time-varying, such as time-varying speed, which makes the hydrodynamic force received by the unmanned ships constantly changing, destroying the stable structure of the system. The change in hydrodynamic force is mainly reflected in the uncertainty of the parameters of the unmanned ship model. This phenomenon causes the unmanned ship system to have strong time-varying nonlinearity.

在目前的无人船艏向控制领域中,所采用的无人船模型主要有Fossen、Norrbin和Nomoto等,为了便于模型参数辨识和系统控制,这些模型都被简化为简单的形式,甚至仅表现出线性时不变特性。然而,在无人船实际航行的过程中,存在各种时变因素如航速、风浪、水流等引起模型参数不确定,会破坏无人船系统的稳定结构,且无人船模型本身固有非线性特性,导致线性时不变特性模型难以表征无人船实际情况,从而进一步增大无人船艏向控制的难度。另一方面,针对具有非线性时变特性的无人船模型设计艏向控制器以抑制外界干扰和参数不确定性,可以采用滑模控制、反步法和H∞鲁棒控制等。然而,滑模控制方法有着切换系统时存在抖振的致命缺点,反步法需要关于船舶航向控制系统的模型及其变化参数的精确信息,这在实际应用中十分困难。H∞鲁棒控制方法通过抑制干扰到期望输出之间传递函数的最大增益,从而达到抑制扰动和不确定性的目的,具有强鲁棒性、控制效果好和不依赖扰动模型的优点。在设计H∞鲁棒控制器时,需要将模型代入至线性变参数系统的控制框架内,以此描述被控模型的非线性时变特性,但其本质是在模型时变参数特定工作点的周围对模型进行线性化处理,因此线性变参数系统的控制框架不能反映被控模型的完整动态特性,尤其是非线性特性。In the current field of unmanned ship heading control, the unmanned ship models used mainly include Fossen, Norrbin and Nomoto. In order to facilitate model parameter identification and system control, these models are simplified to simple forms, and even only show linear time-invariant characteristics. However, in the actual navigation process of unmanned ships, there are various time-varying factors such as speed, wind and waves, and water flow, which cause model parameter uncertainty, which will destroy the stable structure of the unmanned ship system. In addition, the unmanned ship model itself has inherent nonlinear characteristics, which makes it difficult for the linear time-invariant characteristic model to characterize the actual situation of the unmanned ship, thereby further increasing the difficulty of unmanned ship heading control. On the other hand, the heading controller is designed for the unmanned ship model with nonlinear time-varying characteristics to suppress external interference and parameter uncertainty. Sliding mode control, backstepping method and H∞ robust control can be used. However, the sliding mode control method has the fatal disadvantage of chattering when switching the system, and the backstepping method requires accurate information about the model of the ship heading control system and its changing parameters, which is very difficult in practical applications. The H∞ robust control method suppresses disturbances and uncertainties by suppressing the maximum gain of the transfer function from disturbance to desired output. It has the advantages of strong robustness, good control effect and independence from disturbance model. When designing an H∞ robust controller, the model needs to be substituted into the control framework of a linear variable parameter system to describe the nonlinear time-varying characteristics of the controlled model. However, its essence is to linearize the model around a specific operating point of the model's time-varying parameters. Therefore, the control framework of a linear variable parameter system cannot reflect the complete dynamic characteristics of the controlled model, especially the nonlinear characteristics.

发明内容Summary of the invention

为解决现有技术存在的缺陷和不足的问题,本发明提出一种基于非线性变参数无人船模型的H∞艏向控制方法,针对一种小型欠驱动无人船,建立非线性变参数艏向控制模型,并设计非线性变参数H∞鲁棒控制器使系统获得干扰抑制与鲁棒性能,实现艏向的快速、高精度跟踪。In order to solve the defects and shortcomings of the prior art, the present invention proposes an H∞ heading control method based on a nonlinear variable parameter unmanned ship model. Aiming at a small under-actuated unmanned ship, a nonlinear variable parameter heading control model is established, and a nonlinear variable parameter H∞ robust controller is designed to enable the system to obtain interference suppression and robust performance, thereby achieving fast and high-precision tracking of the heading.

考虑水流对船舶的冲击复杂多变,利用低展弦比机翼理论得到船舶所受环流力,进而建立非线性变参数艏向控制模型。针对所建立的非线性变参数艏向控制模型,通过参数相关Lyapunov函数,推导出满足系统H∞鲁棒稳定的非线性矩阵不等式条件,并通过MATLAB中的SOS工具箱求解非线性矩阵不等式条件得到非线性变参数H∞鲁棒控制器。非线性变参数艏向控制模型是基于低展弦比机翼理论建立的,水动力项具有较为明确的物理意义,可以充分描述船舶受到的水动力,且体现了无人船的非线性时变特性。非线性变参数H∞鲁棒控制器可以在保证系统稳定的前提下很好的抑制外界干扰,并且通过包含了系统的非线性状态和时变参数,从而克服系统固有的非线性和时变参数引起的参数摄动现象,从而实现对给定无人船艏向角跟踪,保证快速、平滑和高精度的艏向跟踪性能。本发明使无人船具有新的非线性变参数艏向控制模型,同时为该模型的艏向控制问题提出了非线性变参数H∞鲁棒控制解决方案。Considering the complex and changeable impact of water flow on ships, the circulation force on the ship is obtained by using the low aspect ratio wing theory, and then a nonlinear variable parameter bow control model is established. According to the established nonlinear variable parameter bow control model, the nonlinear matrix inequality conditions that satisfy the H∞ robust stability of the system are derived through the parameter-related Lyapunov function, and the nonlinear variable parameter H∞ robust controller is obtained by solving the nonlinear matrix inequality conditions through the SOS toolbox in MATLAB. The nonlinear variable parameter bow control model is established based on the low aspect ratio wing theory. The hydrodynamic term has a relatively clear physical meaning, which can fully describe the hydrodynamic force on the ship and reflect the nonlinear time-varying characteristics of the unmanned ship. The nonlinear variable parameter H∞ robust controller can well suppress external interference while ensuring the stability of the system, and by including the nonlinear state and time-varying parameters of the system, it overcomes the parameter perturbation phenomenon caused by the inherent nonlinearity and time-varying parameters of the system, thereby realizing the tracking of the given unmanned ship's bow angle and ensuring fast, smooth and high-precision bow tracking performance. The present invention enables the unmanned ship to have a new nonlinear variable parameter heading control model, and at the same time proposes a nonlinear variable parameter H∞ robust control solution for the heading control problem of the model.

本发明解决其技术问题采用的技术方案是:The technical solution adopted by the present invention to solve the technical problem is:

一种基于非线性变参数无人船模型的H∞艏向控制方法,其特征在于,用于小型欠驱动无人船的控制系统,,基于非线性变参数艏向控制模型和非线性变参数H∞鲁棒控制器,具体包括以下步骤:A H∞ heading control method based on a nonlinear variable parameter unmanned ship model, characterized in that the control system for a small underactuated unmanned ship is based on a nonlinear variable parameter heading control model and a nonlinear variable parameter H∞ robust controller, and specifically includes the following steps:

步骤S1:建立Fossen框架下的船舶动力学模型,在建立的过程中忽略垂荡、横摇和纵摇运动,将六自由度无人船简化为与纵荡、横荡与艏摇运动相关的三自由度动力学模型;Step S1: Establish a ship dynamics model under the Fossen framework, ignore heave, roll and pitch motions during the establishment process, and simplify the six-degree-of-freedom unmanned ship into a three-degree-of-freedom dynamics model related to pitch, roll and bow motions;

步骤S2:根据低展弦比机翼理论,将建立的三自由度动力学模型中的水动力阻尼矩阵分解为船舶所受环流力和横流阻力矩阵,得到基于低展弦比机翼的三自由度无人船模型;Step S2: According to the low aspect ratio wing theory, the hydrodynamic damping matrix in the established three-degree-of-freedom dynamic model is decomposed into the circulation force and cross-flow resistance matrix of the ship, and a three-degree-of-freedom unmanned ship model based on the low aspect ratio wing is obtained;

步骤S3:根据小型欠驱动无人船的实际物理特性,忽略横荡运动与横流阻力,将基于低展弦比机翼的三自由度无人船模型分解得到操纵动力学模型;Step S3: According to the actual physical characteristics of the small underactuated unmanned ship, ignoring the swaying motion and cross-flow resistance, the three-degree-of-freedom unmanned ship model based on the low-aspect-ratio wing is decomposed to obtain a maneuvering dynamics model;

步骤S4:根据艏向角与艏向角速度的关系,对操纵动力学模型引入艏向角状态变量,并以艏向角误差作为反馈,得到非线性变参数艏向误差模型;Step S4: According to the relationship between the heading angle and the heading angular velocity, the heading angle state variable is introduced into the maneuvering dynamics model, and the heading angle error is used as feedback to obtain a nonlinear variable parameter heading error model;

步骤S5:将假定的非线性状态反馈控制器代入非线性变参数艏向误差模型,建立非线性变参数艏向控制闭环系统;Step S5: Substituting the assumed nonlinear state feedback controller into the nonlinear variable parameter heading error model to establish a nonlinear variable parameter heading control closed-loop system;

步骤S6:构建状态和参数相关的李雅普诺夫函数,对非线性变参数艏向控制闭环系统进行H∞鲁棒性能和稳定性证明,并推导出系统具有H∞鲁棒稳定性的条件;Step S6: construct the state and parameter-related Lyapunov function, prove the H∞ robust performance and stability of the nonlinear variable parameter heading control closed-loop system, and derive the conditions for the system to have H∞ robust stability;

步骤S7:可以利用Schur补、SOS等引理,将H∞鲁棒稳定性条件转换为多项式线性矩阵不等式;Step S7: The H∞ robust stability condition can be converted into a polynomial linear matrix inequality by using the Schur complement, SOS and other lemmas;

步骤S8:通过MATLAB软件中的SOS工具箱对多项式线性矩阵不等式进行求解,得到步骤S5假定的非线性状态反馈控制器,作为非线性变参数H∞鲁棒控制器;Step S8: Solve the polynomial linear matrix inequality by using the SOS toolbox in MATLAB software to obtain the nonlinear state feedback controller assumed in step S5 as a nonlinear variable parameter H∞ robust controller;

步骤S9:无人船系统根据测量的艏向角和艏向角速度作为反馈,分析得到艏向角误差输入到H∞鲁棒控制器中,以实现对给定艏向角的跟踪、外部扰动和参数摄动的抑制。Step S9: The unmanned ship system uses the measured heading angle and heading angular velocity as feedback, analyzes the heading angle error and inputs it into the H∞ robust controller to achieve tracking of the given heading angle, suppression of external disturbances and parameter perturbations.

进一步地,步骤S1得到的三自由度动力学模型为:Furthermore, the three-degree-of-freedom dynamic model obtained in step S1 is:

Figure BDA0003941397160000031
Figure BDA0003941397160000031

其中,v=[u,v,r]T是无人船状态矢量,u,v,r分别为纵荡速度,横荡速度和艏摇速度;M∈R3×3为惯性矩阵;C(v)∈R3×3为科里奥利向心力矩阵;D(v)∈R3×3为水动力阻尼矩阵;τ=[τuvr]T为推进器输出力矩,τuvr分别为纵荡、横荡和艏摇运动的力矩;τw=[τuwvwrw]T为外界干扰,τuwvwrw分别为纵荡、横荡和艏摇运动的外界干扰。Among them, v = [u, v, r] T is the state vector of the unmanned ship, u, v, r are the surge velocity, sway velocity and bow pitch velocity respectively; M∈R 3×3 is the inertia matrix; C(v)∈R 3×3 is the Coriolis centripetal force matrix; D(v)∈R 3×3 is the hydrodynamic damping matrix; τ = [τ uvr ] T is the thruster output torque, τ uvr are the torques of surge, sway and bow pitch motion respectively; τ w = [τ uwvwrw ] T is the external disturbance, τ uwvwrw are the external disturbances of surge, sway and bow pitch motion respectively.

进一步地,步骤S2具体为:根据低展弦比机翼理论,将船舶模型中的水动力阻尼矩阵用船舶所受环流、横流力代替,即令:Further, step S2 is specifically as follows: according to the low aspect ratio wing theory, the hydrodynamic damping matrix in the ship model is replaced by the circulation and cross flow forces on the ship, that is,

D(v)=FLD(v)+Fcf(v)D(v)=F LD (v)+F cf (v)

其中,Fcf(v)为船舶所受横流阻力;FLD(v)=FL(v)+FD(v)为船舶所受环流力,FL(v)=[XL,YL,NL]T为船舶所受环流升力,FD(v)=[XD,YD,ND]T为船舶所受环流阻力;XL和XD分别为纵荡运动的环流升力和阻力,YL和YD分别为横荡运动的环流升力和阻力,NL和ND分别艏摇运动的环流升力矩和阻力矩。Wherein, F cf (v) is the transverse flow resistance of the ship; F LD (v) = F L (v) + F D (v) is the circulating force of the ship, F L (v) = [X L , Y L , N L ] T is the circulating lift of the ship, and F D (v) = [X D , Y D , N D ] T is the circulating resistance of the ship; XL and X D are the circulating lift and drag of surge motion respectively, Y L and Y D are the circulating lift and drag of surge motion respectively, and NL and ND are the circulating lift moment and drag moment of yaw motion respectively.

进一步地,步骤S3的操纵动力学模型具体为:Furthermore, the manipulation dynamics model of step S3 is specifically:

Figure BDA0003941397160000032
Figure BDA0003941397160000032

其中,

Figure BDA0003941397160000033
为系统参数,m为船舶质量,xg为重心到船舶坐标原点的x轴距离,
Figure BDA0003941397160000034
为附加质量系数,N(u,v,r)为艏摇环流力矩系数。in,
Figure BDA0003941397160000033
is the system parameter, m is the ship mass, xg is the x-axis distance from the center of gravity to the origin of the ship coordinate system,
Figure BDA0003941397160000034
is the additional mass coefficient, and N (u,v,r) is the yaw circulation moment coefficient.

进一步地,在步骤S4中,所述非线性变参数艏向误差模型为:Further, in step S4, the nonlinear variable parameter heading error model is:

Figure BDA0003941397160000035
Figure BDA0003941397160000035

其中,x=[r,ψe]T为状态变量;uτ=c4τr为执行器输出;w=c4τrw为总外部干扰;z=ψe为被控输出,纵荡速度u为时变参数记作u(t);

Figure BDA0003941397160000041
B1=B2=[1,0]T,C=[0,1]为非线性变参数模型系统矩阵。Among them, x = [r, ψ e ] T is the state variable; u τ = c 4 τ r is the actuator output; w = c 4 τ rw is the total external disturbance; z = ψ e is the controlled output, and the sway speed u is the time-varying parameter denoted as u(t);
Figure BDA0003941397160000041
B 1 =B 2 =[1,0] T , C =[0,1] is the nonlinear variable parameter model system matrix.

进一步地,在步骤S5中,所述假定的非线性状态反馈控制器具体为:对于无人船艏向控制系统,若存在一个依赖时变参数的对称正定多项式矩阵P(u(t)),一个多项式矩阵

Figure BDA0003941397160000042
使得多项式线性矩阵不等式:Further, in step S5, the assumed nonlinear state feedback controller is specifically: for the unmanned ship heading control system, if there is a symmetric positive definite polynomial matrix P(u(t)) that depends on time-varying parameters, a polynomial matrix
Figure BDA0003941397160000042
This makes the polynomial linear matrix inequality:

Figure BDA0003941397160000043
Figure BDA0003941397160000043

Figure BDA0003941397160000044
Figure BDA0003941397160000044

成立,则该控制系统渐近稳定,且对于外部扰动w的H∞范数小于γ,γ为H∞性能指标,此时闭环系统的状态反馈控制器形式为:If it holds, the control system is asymptotically stable, and the H∞ norm of the external disturbance w is less than γ, where γ is the H∞ performance index. At this time, the state feedback controller of the closed-loop system is in the form of:

Figure BDA0003941397160000045
Figure BDA0003941397160000045

进一步地,在步骤S8中,对于所述多项式线性矩阵不等式,将其转化为平方和形式通过MATLAB软件中的SOS工具箱求解,得到一个对称正定多项式矩阵P(u(t))和一个多项式矩阵

Figure BDA0003941397160000046
从而设计出控制律uτ。Further, in step S8, the polynomial linear matrix inequality is converted into a square sum form and solved by the SOS toolbox in MATLAB software to obtain a symmetric positive definite polynomial matrix P(u(t)) and a polynomial matrix
Figure BDA0003941397160000046
Thus the control law u τ is designed.

与现有技术相比,本发明及其优选方案具有以下有益效果:Compared with the prior art, the present invention and its preferred embodiment have the following beneficial effects:

1、提供的非线性变参数艏向控制模型在物理层面建立了船舶水动力机理模型,代替了以往的水动力经验模型,具有较为明确的物理意义;1. The nonlinear variable parameter bow control model provided establishes a ship hydrodynamic mechanism model at the physical level, replacing the previous hydrodynamic empirical model and having a clearer physical meaning;

2、提供的非线性变参数艏向控制模型相比于其他艏向控制模型,不仅具有简洁的形式,而且充分还原了无人船系统的非线性时变特性,有助于设计能精准调节艏向的控制器;2. Compared with other heading control models, the nonlinear variable parameter heading control model provided not only has a simple form, but also fully restores the nonlinear time-varying characteristics of the unmanned ship system, which is helpful for designing a controller that can accurately adjust the heading;

3、提供的非线性变参数H∞鲁棒控制器通过纳入系统状态和时变参数及其导数,不需对系统线性化处理,从而不再隐藏系统的非线性时变特性,降低了控制器设计的保守性,提高了无人船艏向控制系统的动态性能和鲁棒性;3. The nonlinear variable parameter H∞ robust controller provided by the system incorporates the system state and time-varying parameters and their derivatives, and does not need to linearize the system, thus no longer hiding the nonlinear time-varying characteristics of the system, reducing the conservatism of the controller design and improving the dynamic performance and robustness of the unmanned ship bow control system;

4、将非线性变参数H∞鲁棒控制器的设计转化为多项式线性矩阵不等式条件的求解问题,可以利用MATLAB软件的SOS工具箱进行求解,降低了计算难度,易于工程实现。4. The design of nonlinear variable parameter H∞ robust controller is transformed into the problem of solving polynomial linear matrix inequality conditions, which can be solved using the SOS toolbox of MATLAB software, reducing the computational difficulty and facilitating engineering implementation.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

下面结合附图和具体实施方式对本发明进一步详细的说明:The present invention is further described in detail below with reference to the accompanying drawings and specific embodiments:

图1是本发明实施例中基于非线性变参数无人船模型的H∞艏向控制方法结构图;1 is a structural diagram of an H∞ heading control method based on a nonlinear variable parameter unmanned ship model in an embodiment of the present invention;

图2为本发明实施例中非线性变参数无人船模型建立流程图;FIG2 is a flow chart of establishing a nonlinear variable parameter unmanned ship model according to an embodiment of the present invention;

图3为本发明实施例中非线性变参数H∞鲁棒控制器设计流程图;FIG3 is a flow chart of the design of a nonlinear variable parameter H∞ robust controller according to an embodiment of the present invention;

图4是本发明实施例艏向跟踪仿真图;FIG4 is a simulation diagram of heading tracking according to an embodiment of the present invention;

图5是本发明实施例跟踪艏向角实际误差仿真图。FIG5 is a simulation diagram of the actual error of tracking the heading angle according to an embodiment of the present invention.

具体实施方式DETAILED DESCRIPTION

为让本专利的特征和优点能更明显易懂,下文特举实施例,作详细说明如下:In order to make the features and advantages of this patent more obvious and easy to understand, the following embodiments are specifically described in detail as follows:

应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本说明书使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。It should be noted that the following detailed descriptions are illustrative and are intended to provide further explanation of the present application. Unless otherwise specified, all technical and scientific terms used in this specification have the same meanings as those commonly understood by those skilled in the art to which the present application belongs.

需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present application. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should be understood that when the terms "comprise" and/or "include" are used in this specification, it indicates the presence of features, steps, operations, devices, components and/or combinations thereof.

请参照图1,本发明提供一种基于非线性变参数无人船模型的H∞艏向控制方法,根据低展弦比机翼理论得到船舶所受环流、横流力,建立非线性变参数艏向控制模型,推导出多项式线性矩阵不等式,从而利用MATLAB软件中的SOS工具箱求解,最终得到非线性变参数H∞鲁棒控制器。该方法在建立充分描述无人船非线性时变特性模型的基础上,设计了一种非线性变参数H∞鲁棒控制器,给无人船艏向控制系统提供了良好的稳定性和较强的鲁棒性。具体包括如下步骤:Please refer to Figure 1. The present invention provides an H∞ bow control method based on a nonlinear variable parameter unmanned ship model. According to the low aspect ratio wing theory, the circulation and cross flow forces on the ship are obtained, a nonlinear variable parameter bow control model is established, and a polynomial linear matrix inequality is derived, which is then solved using the SOS toolbox in the MATLAB software to finally obtain a nonlinear variable parameter H∞ robust controller. This method designs a nonlinear variable parameter H∞ robust controller based on the establishment of a model that fully describes the nonlinear time-varying characteristics of the unmanned ship, which provides good stability and strong robustness for the unmanned ship bow control system. Specifically, it includes the following steps:

步骤1:如图2所示,根据Fossen框架的船舶理论,以建立纵荡、横荡和艏摇运动的三自由度动力学模型:Step 1: As shown in Figure 2, according to the ship theory of the Fossen framework, a three-degree-of-freedom dynamic model of surge, sway and pitch motion is established:

Figure BDA0003941397160000051
Figure BDA0003941397160000051

其中,v=[u,v,r]T是无人船状态矢量,u,v,r分别为纵荡速度,横荡速度和艏摇速度;M∈R3×3为惯性矩阵;C(v)∈R3×3为科里奥利向心力矩阵;D(v)∈R3×3为水动力阻尼矩阵;τ=[τuvr]T为推进器输出力矩,τuvr分别为纵荡、横荡和艏摇运动的力矩;τw=[τuwvwrw]T为外界干扰,τuwvwrw分别为纵荡、横荡和艏摇运动的外界干扰。Among them, v = [u, v, r] T is the state vector of the unmanned ship, u, v, r are the surge velocity, sway velocity and bow pitch velocity respectively; M∈R 3×3 is the inertia matrix; C(v)∈R 3×3 is the Coriolis centripetal force matrix; D(v)∈R 3×3 is the hydrodynamic damping matrix; τ = [τ uvr ] T is the thruster output torque, τ uvr are the torques of surge, sway and bow pitch motion respectively; τ w = [τ uwvwrw ] T is the external disturbance, τ uwvwrw are the external disturbances of surge, sway and bow pitch motion respectively.

步骤2:如图2所示,将船舶视为一种低展弦比机翼,根据低展弦比机翼理论,将船舶模型中的水动力阻尼矩阵用船舶所受环流、横流力代替,即令:Step 2: As shown in Figure 2, the ship is regarded as a low aspect ratio wing. According to the low aspect ratio wing theory, the hydrodynamic damping matrix in the ship model is replaced by the circulation and cross flow forces on the ship, that is,

D(v)=FLD(v)+Fcf(v)D(v)=F LD (v)+F cf (v)

其中,Fcf(v)为船舶所受横流阻力;FLD(v)=FL(v)+FD(v)为船舶所受环流力,FL(v)=[XL,YL,NL]T为船舶所受环流升力,FD(v)=[XD,YD,ND]T为船舶所受环流阻力;XL和XD分别为纵荡运动的环流升力和阻力,YL和YD分别为横荡运动的环流升力和阻力,NL和ND分别艏摇运动的环流升力矩和阻力矩。Wherein, F cf (v) is the transverse flow resistance of the ship; F LD (v) = F L (v) + F D (v) is the circulating force of the ship, F L (v) = [X L , Y L , N L ] T is the circulating lift of the ship, and F D (v) = [X D , Y D , N D ] T is the circulating resistance of the ship; XL and X D are the circulating lift and drag of surge motion respectively, Y L and Y D are the circulating lift and drag of surge motion respectively, and NL and ND are the circulating lift moment and drag moment of yaw motion respectively.

步骤3:如图2所示,将船舶所受环流、横流力在船舶不同运动方向上进行受力分析得到基于低展弦比机翼的无人船动力学模型,并根据小型无人船欠驱动的特性将基于低展弦比机翼的无人船动力学模型分解,得到操纵动力学模型:Step 3: As shown in Figure 2, the circulation and cross-flow forces on the ship are analyzed in different directions of movement to obtain the unmanned ship dynamics model based on the low aspect ratio wing. The unmanned ship dynamics model based on the low aspect ratio wing is decomposed according to the under-actuated characteristics of the small unmanned ship to obtain the maneuvering dynamics model:

Figure BDA0003941397160000061
Figure BDA0003941397160000061

其中,

Figure BDA0003941397160000062
为系统参数,m为船舶质量,xg为重心到船舶坐标原点的x轴距离,
Figure BDA0003941397160000063
为附加质量系数,N(u,v,r)为艏摇环流力矩系数。in,
Figure BDA0003941397160000062
is the system parameter, m is the ship mass, xg is the x-axis distance from the center of gravity to the origin of the ship coordinate system,
Figure BDA0003941397160000063
is the additional mass coefficient, and N (u,v,r) is the yaw circulation moment coefficient.

步骤4:如图3所示,根据艏向角与艏向角速度的物理关系,在操纵动力学模型内引入艏向角变量,得到非线性变参数艏向控制模型,并考虑在一个艏向采样步长内给定艏向角不变,根据艏向角误差ψe=ψd-ψ,建立非线性变参数艏向误差模型:Step 4: As shown in Figure 3, according to the physical relationship between the heading angle and the heading angular velocity, the heading angle variable is introduced into the maneuvering dynamics model to obtain a nonlinear variable parameter heading control model. Considering that the heading angle is given unchanged within a heading sampling step, a nonlinear variable parameter heading error model is established according to the heading angle error ψ e = ψ d - ψ:

Figure BDA0003941397160000064
Figure BDA0003941397160000064

其中,x=[r,ψe]T为状态变量;uτ=c4τr为执行器输出;w=c4τrw为总外部干扰;z=ψe为被控输出,纵荡速度u为时变参数记作u(t);

Figure BDA0003941397160000065
B1=B2=[1,0]T,C=[0,1]为非线性变参数模型系统矩阵。Among them, x = [r, ψ e ] T is the state variable; u τ = c 4 τ r is the actuator output; w = c 4 τ rw is the total external disturbance; z = ψ e is the controlled output, and the sway speed u is the time-varying parameter denoted as u(t);
Figure BDA0003941397160000065
B 1 =B 2 =[1,0] T , C =[0,1] is the nonlinear variable parameter model system matrix.

步骤5:如图3所示,针对非线性变参数艏向误差模型,构造李雅普诺夫函数V(x,t)=xTP-1(u(t))x。根据系统H∞稳定条件:(1)V(x,u,t)>0;

Figure BDA0003941397160000066
Figure BDA0003941397160000067
推导出无人船艏向控制系统H∞鲁棒稳定的控制器求解条件为:对于无人船艏向控制系统,若存在一个依赖时变参数的对称正定多项式矩阵P(u(t)),一个多项式矩阵
Figure BDA0003941397160000071
使得多项式线性矩阵不等式:Step 5: As shown in Figure 3, for the nonlinear variable parameter heading error model, construct the Lyapunov function V(x,t)=x T P -1 (u(t))x. According to the system H∞ stability condition: (1) V(x,u,t)>0;
Figure BDA0003941397160000066
Figure BDA0003941397160000067
The controller solution condition for the H∞ robust stability of the unmanned ship bow control system is derived as follows: For the unmanned ship bow control system, if there exists a symmetric positive definite polynomial matrix P(u(t)) that depends on time-varying parameters, a polynomial matrix
Figure BDA0003941397160000071
This makes the polynomial linear matrix inequality:

Figure BDA0003941397160000072
Figure BDA0003941397160000072

Figure BDA0003941397160000073
Figure BDA0003941397160000073

成立,则该控制系统渐近稳定,且对于外部扰动w的H∞范数小于γ,γ为H∞性能指标,此时闭环系统的状态反馈控制器形式为:If it holds, the control system is asymptotically stable, and the H∞ norm of the external disturbance w is less than γ, where γ is the H∞ performance index. At this time, the state feedback controller of the closed-loop system is in the form of:

Figure BDA0003941397160000074
Figure BDA0003941397160000074

步骤6:如图3所示,将多项式线性矩阵不等式转化为平方和形式通过MATLAB软件中的SOS工具箱求解,得到一个对称正定多项式矩阵P(u(t))和一个多项式矩阵

Figure BDA0003941397160000075
从而设计出具有H∞鲁棒性能的控制律uτ。Step 6: As shown in Figure 3, the polynomial linear matrix inequality is converted into a square sum form and solved by the SOS toolbox in MATLAB software to obtain a symmetric positive definite polynomial matrix P(u(t)) and a polynomial matrix
Figure BDA0003941397160000075
Thus, a control law u τ with H∞ robust performance is designed.

步骤7:如图1、图3所示,根据需求设置系统给定艏向角ψd,无人船的艏向角速度r和艏向角误差ψe作为状态反馈,从而搭建无人船艏向控制仿真系统。将求解得到的控制律代入uτ仿真系统,并对参数进行优化以达到预期的控制性能,最终实现对给定艏向角的跟踪,仿真结果如图4、图5所示。Step 7: As shown in Figures 1 and 3, the system is set to give a given heading angle ψ d , the unmanned ship's heading angular velocity r and the heading angle error ψ e as state feedback according to the requirements, so as to build an unmanned ship heading control simulation system. Substitute the solved control law into the u τ simulation system, and optimize the parameters to achieve the expected control performance, and finally achieve tracking of the given heading angle. The simulation results are shown in Figures 4 and 5.

本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present application is described with reference to the flowchart and/or block diagram of the method, device (system) and computer program product according to the embodiment of the present application. It should be understood that each process and/or box in the flowchart and/or block diagram, and the combination of the process and/or box in the flowchart and/or block diagram can be realized by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device produce a device for realizing the function specified in one process or multiple processes in the flowchart and/or one box or multiple boxes in the block diagram.

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including an instruction device that implements the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions may also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, whereby the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

以上所述,仅是本发明的较佳实施例而已,并非是对本发明作其它形式的限制,任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型,仍属于本发明技术方案的保护范围。The above is only a preferred embodiment of the present invention, and does not limit the present invention in other forms. Any technician familiar with the profession may use the above disclosed technical content to change or modify it into an equivalent embodiment with equivalent changes. However, any simple modification, equivalent change and modification made to the above embodiment according to the technical essence of the present invention without departing from the technical solution of the present invention still belongs to the protection scope of the technical solution of the present invention.

本专利不局限于上述最佳实施方式,任何人在本专利的启示下都可以得出其它各种形式的基于非线性变参数无人船模型的H∞艏向控制方法,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本专利的涵盖范围。This patent is not limited to the above-mentioned optimal implementation mode. Anyone can derive other various forms of H∞ bow control methods based on nonlinear variable parameter unmanned ship models under the inspiration of this patent. All equal changes and modifications made according to the scope of the patent application of the present invention should be covered by this patent.

Claims (7)

1. An H infinity heading control method based on a nonlinear variable parameter unmanned ship model is characterized in that the H infinity heading control method is used for a control system of a small under-actuated unmanned ship, is based on the nonlinear variable parameter heading control model and a nonlinear variable parameter H infinity controller, and specifically comprises the following steps:
step S1: building a ship dynamic model under a Fossen frame, neglecting heaving, rolling and pitching motions in the building process, and simplifying the six-freedom-degree unmanned ship into a three-freedom-degree dynamic model related to the heaving, rolling and pitching motions;
step S2: decomposing a hydrodynamic damping matrix in the established three-degree-of-freedom dynamic model into a circulating flow force and a cross flow resistance matrix borne by a ship according to a low-aspect-ratio wing theory to obtain a low-aspect-ratio wing-based three-degree-of-freedom unmanned ship model;
and step S3: according to the actual physical characteristics of the small under-actuated unmanned ship, ignoring the swaying motion and the cross flow resistance, decomposing a three-degree-of-freedom unmanned ship model based on the low-aspect-ratio wing to obtain an operation dynamic model;
and step S4: according to the relation between the heading angle and the heading angular velocity, introducing a heading angle state variable into the control dynamic model, and taking a heading angle error as feedback to obtain a nonlinear variable parameter heading error model;
step S5: substituting the assumed nonlinear state feedback controller into a nonlinear variable parameter heading error model to establish a nonlinear variable parameter heading control closed-loop system;
step S6: constructing a Lyapunov function related to states and parameters, carrying out H-infinity robust performance and stability verification on a nonlinear variable parameter heading control closed-loop system, and deducing the condition that the system has H-infinity robust stability;
step S7: converting the H infinity robust stability condition into a polynomial linear matrix inequality;
step S8: solving a polynomial linear matrix inequality through an SOS tool box in MATLAB software to obtain a nonlinear state feedback controller assumed in the step S5 as a nonlinear variable parameter H infinity robust controller;
step S9: the unmanned ship system takes the measured heading angle and heading angular velocity as feedback, analyzes the obtained heading angle error and inputs the error into an H infinity controller so as to realize the tracking of the given heading angle, the suppression of external disturbance and parameter perturbation.
2. The H infinity heading control method based on nonlinear variable parameter unmanned ship model according to claim 1 wherein: the three-degree-of-freedom dynamic model obtained in the step S1 is as follows:
Figure FDA0003941397150000011
wherein v = [ u, v, r] T The unmanned ship state vector is shown, and u, v and r are respectively a longitudinal oscillation speed, a transverse oscillation speed and a heading speed; m belongs to R 3×3 Is an inertia matrix; c (v) is epsilon to R 3×3 Is a coriolis centripetal force matrix; d (v) epsilon R 3×3 Is a hydrodynamic damping matrix; τ = [ τ ] uvr ] T For the output torque of the propeller, τ uvr Moments of surge, sway and yaw motions, respectively; tau is w =[τ uwvwrw ] T For external interference, τ uwvwrw External disturbances of surge, sway and yaw motions, respectively.
3. The H infinity heading control method based on the nonlinear variable parameter unmanned ship model according to claim 2, characterized in that: the step S2 specifically comprises the following steps: according to the low aspect ratio wing theory, the hydrodynamic damping matrix in the ship model is replaced by circulation and cross flow force borne by the ship, namely:
D(v)=F LD (v)+F cf (v)
wherein, F cf (v) The ship is subjected to cross flow resistance; f LD (v)=F L (v)+F D (v) Is the circulation force to which the ship is subjected, F L (v)=[X L ,Y L ,N L ] T Is the circulation lift force borne by the ship, F D (v)=[X D ,Y D ,N D ] T The circulation resistance of the ship; x L And X D Lift and drag of the surging motion, Y, respectively L And Y D Circulating lift and drag, N, respectively, of the oscillatory motion L And N D And respectively the circular flow lifting moment and the resisting moment of the yawing motion.
4. The H infinity heading control method based on nonlinear variable parameter unmanned ship model according to claim 3 wherein: the operation dynamic model of the step S3 is specifically as follows:
Figure FDA0003941397150000021
wherein,
Figure FDA0003941397150000022
is the system parameter, m is the vessel mass, x g Is the x-axis distance from the center of gravity to the origin of coordinates of the vessel>
Figure FDA0003941397150000023
To add a mass coefficient, N (u,v,r) Is the coefficient of the heading circulation moment.
5. The H infinity heading control method based on nonlinear variable parameter unmanned ship model according to claim 4 wherein: in step S4, the nonlinear variable parameter heading error model is:
Figure FDA0003941397150000024
wherein x = [ r, ψ) e ] T Is a state variable; u. u τ =c 4 τ r Outputting for the actuator; w = c 4 τ rw Is the total external interference; z = ψ e For controlled output, the surging speed u is a time-varying parameter and is recorded as u (t);
Figure FDA0003941397150000025
B 1 =B 2 =[1,0] T ,C=[0,1]is a nonlinear variable parameter model system matrix.
6. The H infinity heading control method based on nonlinear variable parameter unmanned ship model according to claim 5 wherein: in step S5, the assumed nonlinear state feedback controller is specifically: for the unmanned ship heading control system, if the unmanned ship heading control system existsA symmetric positive definite polynomial matrix P (u (t)) dependent on time-varying parameters, a polynomial matrix
Figure FDA0003941397150000026
Such that the polynomial linear matrix inequality:
Figure FDA0003941397150000031
Figure FDA0003941397150000032
if yes, the control system is asymptotically stable, the H infinity norm for the external disturbance w is less than gamma, gamma is the H infinity performance index, and the state feedback controller of the closed-loop system has the form:
Figure FDA0003941397150000033
7. the H infinity heading control method based on nonlinear variable parameter unmanned ship model according to claim 6 wherein: in step S8, the polynomial linear matrix inequality is converted into a sum of squares form and solved by SOS toolbox in MATLAB software to obtain a symmetric positive definite polynomial matrix P (u (t)) and a polynomial matrix
Figure FDA0003941397150000034
Thereby designing a control law u τ 。/>
CN202211421529.1A 2022-11-14 2022-11-14 H∞ heading control method based on nonlinear variable parameter unmanned ship model Pending CN115903802A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211421529.1A CN115903802A (en) 2022-11-14 2022-11-14 H∞ heading control method based on nonlinear variable parameter unmanned ship model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211421529.1A CN115903802A (en) 2022-11-14 2022-11-14 H∞ heading control method based on nonlinear variable parameter unmanned ship model

Publications (1)

Publication Number Publication Date
CN115903802A true CN115903802A (en) 2023-04-04

Family

ID=86493347

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211421529.1A Pending CN115903802A (en) 2022-11-14 2022-11-14 H∞ heading control method based on nonlinear variable parameter unmanned ship model

Country Status (1)

Country Link
CN (1) CN115903802A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118795916A (en) * 2024-09-13 2024-10-18 中国海洋大学 Solution method for bow-swing attitude controller of underwater robot

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118795916A (en) * 2024-09-13 2024-10-18 中国海洋大学 Solution method for bow-swing attitude controller of underwater robot
CN118795916B (en) * 2024-09-13 2025-01-28 中国海洋大学 Solution method for bow-swing attitude controller of underwater robot

Similar Documents

Publication Publication Date Title
CN108008628B (en) A Preset Performance Control Method for Uncertain Underactuated Unmanned Vehicle System
Liu et al. Disturbance observer-based adaptive intelligent control of marine vessel with position and heading constraint condition related to desired output
CN109656143B (en) Self-adaptive tracking control method and system for sea surface full-drive ship
CN108803632A (en) Unmanned surface vehicle total state constrained trajectory tracking and controlling method based on saturation compensation technology
CN111736617A (en) A speed observer-based trajectory tracking control method for benthic underwater robot with preset performance
CN110262494B (en) Collaborative learning and formation control method for isomorphic multi-unmanned ship system
Mu et al. Tracking control of podded propulsion unmanned surface vehicle with unknown dynamics and disturbance under input saturation
CN106054884B (en) L1 self-adaptive ship dynamic positioning Double Loop Control System neural network based
CN110865539A (en) Unmanned ship tracking error constraint control method under random interference
CN106154835A (en) A kind of underwater research vehicle TSM control method estimated based on time delay
CN109917657A (en) Anti-jamming control method, device and electronic equipment for dynamic positioning ship
Liu et al. Finite-time prescribed performance trajectory tracking control for underactuated autonomous underwater vehicles based on a tan-type barrier Lyapunov function
CN109828586B (en) Unmanned ship robust H-infinity course control method based on nonlinear uncertainty
CN115903802A (en) H∞ heading control method based on nonlinear variable parameter unmanned ship model
Wang et al. State-constrained control strategy for safe navigation trajectory tracking of hovercraft based on improved barrier Lyapunov function
Wang et al. Output-feedback control for cooperative diving of saucer-type underwater gliders based on a fuzzy observer and event-triggered communication
Gao et al. Disturbance observer based finite-time trajectory tracking control of unmanned surface vehicles with unknown dead-zones
CN105700356A (en) UUV positioning-to-line recovery method based on robustness constraint model prediction control
CN111596670A (en) Automatic berthing method of under-actuated surface craft based on stabilization control
CN116482977A (en) Novel nonsingular rapid terminal sliding mode control method and system for underwater robot
Dong et al. Model predictive control for ship path tracking with disturbances
CN116224990A (en) A track tracking method for surface ships
Wang et al. Adaptive Integral Sliding Mode Control for Attitude Tracking of Underwater Robots With Large Range Pitch Variations in Confined Spaces
Li et al. Underactuated Autonomous Underwater Vehicle Trajectory Tracking Control in Three Dimensions Based on Fractional Order Sliding Mode
Liu et al. Robust event-triggered formation control for underactuated ships under the practice constraint

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination