CN111665719A - Supply ship synchronous control algorithm with timeliness and stability - Google Patents

Supply ship synchronous control algorithm with timeliness and stability Download PDF

Info

Publication number
CN111665719A
CN111665719A CN202010530817.5A CN202010530817A CN111665719A CN 111665719 A CN111665719 A CN 111665719A CN 202010530817 A CN202010530817 A CN 202010530817A CN 111665719 A CN111665719 A CN 111665719A
Authority
CN
China
Prior art keywords
formula
tender
control
vessel
matrix
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.)
Granted
Application number
CN202010530817.5A
Other languages
Chinese (zh)
Other versions
CN111665719B (en
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.)
Dalian Maritime University
Original Assignee
Dalian Maritime 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 Dalian Maritime University filed Critical Dalian Maritime University
Priority to CN202010530817.5A priority Critical patent/CN111665719B/en
Publication of CN111665719A publication Critical patent/CN111665719A/en
Application granted granted Critical
Publication of CN111665719B publication Critical patent/CN111665719B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • 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

Abstract

The invention provides a tender vessel synchronous control algorithm with timeliness and stability, which comprises a terminal cost function module, a tender vessel model prediction control module and a neural dynamic optimization module. The terminal cost function module ensures the closed loop stability of the supply ship control system by using the terminal cost function; the Model prediction Control module converts the synchronous Control problem of the tender boat into a tracking Control problem following a pilot, and designs a tender boat Control algorithm based on a Model Predictive Control (MPC) algorithm, so that the anti-interference capability of the tender boat during the execution of a tender task is improved, the speeds of two boats can be synchronized during the execution of the tender task, and the tender efficiency is improved; the neural dynamic optimization module solves the problems of large calculation amount and low calculation speed in the traditional model predictive control algorithm by using a parallel calculation framework of a recurrent neural network, and can enable a tender boat to better cope with emergency situations during tender.

Description

Supply ship synchronous control algorithm with timeliness and stability
Technical Field
The invention relates to the technical field of automatic control, in particular to a synchronous control algorithm for a tender boat with timeliness and stability.
Background
Chinese patent CN 201610859143.7 discloses a generalized prediction adaptive tender vessel course control method based on event driving, and the invention provides an event driving-based generalized prediction adaptive tender vessel course control method, firstly, a discrete event trigger working in a discrete state is selected, a trigger judgment function under a variable threshold value is designed, and when the trigger time is reached, the event trigger judges whether the current state meets the trigger condition or not through the trigger judgment function; secondly, a controlled autoregressive integral sliding average model of a rudder angle-course is obtained through a supply ship low-frequency motion mathematical model and is used as a prediction model, and parameters of the prediction model are estimated on line by adopting a forgetting factor recursion least square method; thirdly, the state to be controlled is sent to the controller by combining with the event-driven trigger, and the controller outputs rudder angle control increment and control quantity after resolving through a GPC adaptive algorithm, so that the tender vessel can quickly finish the approaching stage and the parallel stage and keep course tracking of the tender stage. The invention mainly researches the course control problem of the supply ship, but in the course of sailing supply of the supply ship, only the course tracking effect in the supply stage is not enough, the speed synchronization problem of the supply ship and the supplied ship and the position relation between the two ships are considered, and the timeliness and the stability of the algorithm are increased, so that the supply ship can efficiently and safely complete the supply task.
Disclosure of Invention
According to the technical problems provided by the invention, a supply ship synchronous control algorithm with timeliness and stability is provided. The invention mainly utilizes a tender boat synchronous control algorithm with timeliness and stability, and is characterized by comprising the following steps: the system comprises a terminal cost function module, a tender ship model prediction control module and a neural dynamic optimization module.
For the tracking control problem of the supply ship, a three-degree-of-freedom ship motion model is adopted:
Figure BDA0002535343600000011
wherein η ═ x, y, ψ]TX represents the lateral displacement of the vessel, y represents the longitudinal displacement of the vessel, ψ represents the yaw angle; u, v, r]TU represents the forward speed of the vessel, v represents the lateral speed of the vessel, and r represents the yaw rate of the vessel; j (ψ) represents a rotation matrix,
Figure BDA0002535343600000021
the tender vessel is typically operated at low speed, and therefore the dynamic model of vessel motion is summarized as:
Figure BDA0002535343600000022
in the formula (I), the compound is shown in the specification,
Figure BDA0002535343600000023
a matrix of the inertia is represented and,
Figure BDA0002535343600000024
a damping matrix is represented which,
Figure BDA0002535343600000025
representing the Coriolis and central centripetal force matrix, τT=[τuτvτr]TIndicating control force and moment, X(·),Y(·),N(·)Is a hydrodynamic parameter.
The leader ship serves as a reference of a following ship, and a mathematical model of the leader ship is defined as:
Figure BDA0002535343600000026
coordinate conversion is performed on equations (1) and (3) by the following equation (4):
z1=xcosψ+ysinψ
z2=-xsinψ+ycosψ (4)
z3=ψ
offset (x) of displacemento,yo) In combination with equation (4) above, equation (4) can be converted to:
z1r=(xr-xo)cosψr+(yr-yo)sinψr
z2r=-(xr-xo)sinψr+(yr-yo)cosψr
z3r=ψr(5)
defining a synchronous control error system for the replenishment of a dynamically positioned vessel as:
zie=zi-ziri=1,2,3
ue=u-ur(6)
ve=v-vr
re=r-rr
the formula (6) is differentiated and linearized to obtain:
Figure BDA0002535343600000031
in the formula, xt=yt=[z1ez2ez3euevere]T
Figure BDA0002535343600000032
Wherein the content of the first and second substances,
Figure BDA0002535343600000033
d11=-Xu,d22=-Yv,d33=-Nr.
according to the actual requirement, selecting proper sampling time, the equation (8) can be converted into discrete form:
x(k+1)=Ax(k)+Bτ(k) (9)
in the formula, x, y ∈ R6×1,A∈R6×6,B∈R6×3K denotes sampling time points, and A and B are discrete system parameter matrixes.
In the invention, the stability of the control system is ensured by designing a secondary terminal cost function, and in the terminal cost function module, a linear state feedback law is firstly designed for the control system:
τ=Kx (10)
in the formula, K represents a vector for controlling the gain. The state matrix P of the terminal cost is obtained by solving the lyapunov equation (11)):
(Ak+κI)TP+P(Ak+κI)=-(Qlq+KTRlqK) (11)
in the formula, Ak=A+BK,QlqAnd RlqA weight matrix representing the output and input, κ being selected by equation (12):
κ<-λmax(Ak) (12)
thereby ensuring AkThe real part of all eigenvalues of + kappa I is negative, where lambda ismax(Ak) Is represented by AkMaximum eigenvalue of the real part.
Further, in order to improve the robustness and the anti-interference capability of the system, a tender vessel control algorithm is designed based on a model prediction control algorithm. Defining a prediction time domain as N in the supply ship model prediction control modulepControl time domain as NcIn the prediction calculation process, let τ (k + N)c)=τ(k+Nc+1)=...=τ(k+Nc+Np) Then, the future state information of the system can be calculated sequentially by a certain prediction model (x (k +1) ═ ax (k) + B τ (k)):
Figure BDA0002535343600000041
where x (k + i | k) i is 1,2pIndicates future state information x (k + i) i, which is predicted from the state information at the sampling time point k, is 1,2p
According to the above formula (13), the following definitions are made:
Figure BDA0002535343600000042
rewriting the above formula (13) into a matrix form as shown in formula (15):
X=Fx(k)+Φ (15)
in the formula (I), the compound is shown in the specification,
Figure BDA0002535343600000043
then, the synchronous multivariable control problem of the supply ship is converted into an optimal control problem with terminal cost:
Figure BDA0002535343600000044
the constraint conditions are as follows:
min≤≤max(17)
in the formula, QmpcAnd RmpcIs an output weight matrix and an input weight matrix,
Figure BDA0002535343600000051
PXpis the terminal cost, XpIs the terminal state of X and P is the state matrix of the terminal cost function determined by equation (11).
Furthermore, in order to solve the problems of large calculation amount and low calculation speed in the traditional model predictive control algorithm, in the neural dynamic optimization system module, the parallel calculation architecture of the recurrent neural network is utilized to accelerate the calculation speed, so that the proposed algorithm has timeliness. In the neural dynamics optimization system module, bringing formula (15) into formula (16) results in:
Figure BDA0002535343600000052
by definition
Figure BDA0002535343600000053
Equation (18) is simplified to:
J=TH+W+ξ (20)
the objective function (equation (20)) is solved iteratively by establishing a neural dynamic optimization system (equation (21)), so that the calculation speed is accelerated:
Figure BDA0002535343600000054
in the formula (I), the compound is shown in the specification,
Figure BDA0002535343600000055
is the convergence rate, and Ω () is taken to be ▽ J, PΩIs the projection operator, and the projection operator,
Figure BDA0002535343600000061
compared with the prior art, the invention has the following advantages:
(1) the method is based on a model prediction control algorithm, has robustness, and can excellently complete the replenishment task under the condition that external interference exists in the ship;
(2) the method solves the problems of large calculation amount and low calculation speed in the traditional model predictive control algorithm by utilizing the parallel calculation architecture of the recurrent neural network, can ensure that the proposed algorithm has timeliness, and ensures that the supply vessel adopting the method can better cope with changeable external environments;
(3) the invention ensures the closed loop stability of the control system by using the terminal cost function technology, so that the proposed algorithm has stability;
(4) the technology provided by the invention is used for controlling the supply ship, so that the supply ship can automatically follow the leader ship from different initial positions under the condition of interference and achieve the expected target of synchronous speed, the working strength of a rudder during supply can be effectively reduced, the supply ship can reach a supply state more quickly, the anti-interference capability of the supply ship is increased, the speed of the two ships can be ensured to be synchronous during supply, the influence of waves on supply is reduced, and the supply ship can finish the supply task efficiently.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a block diagram of an MPC design with timeliness and stability of the present invention.
FIG. 2 is a performance diagram of synchronous control at different starting points according to the present invention.
FIG. 3 is a schematic diagram of the velocity synchronization process at different initial points according to the present invention.
FIG. 4 is a graphical representation of the response of the control input to various initiation points of the present invention.
Detailed Description
In order to make the technical solutions of the present invention better understood, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be noted that the terms "first," "second," and the like in the description and claims of the present invention and in the drawings described above are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used is interchangeable under appropriate circumstances such that the embodiments of the invention described herein are capable of operation in sequences other than those illustrated or described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
As shown in fig. 1-4, the present invention provides a tender vessel synchronization control algorithm with timeliness and stability, which includes: the system comprises a terminal cost function module, a tender ship model prediction control module and a neural dynamic optimization module.
In the application, (1) in a terminal cost function module, the closed loop stability of a control system is guaranteed by designing a secondary terminal cost function; (2) in a model prediction control module of the supply ship, converting a synchronous control problem of the supply ship into a tracking control problem following a pilot, and solving the problem based on a model prediction control algorithm, so that the anti-interference capability of the supply ship is improved, the speed of two ships can be synchronized when a supply task is carried out, and the supply efficiency and the safety are improved; (3) in the neural dynamic optimization module, a parallel computing architecture of a recurrent neural network is used for quickly solving an objective function in the tender vessel model prediction control module, so that the ship can quickly react to a changeable external environment.
For the tracking control problem of the supply ship, a three-degree-of-freedom ship motion model is adopted:
Figure BDA0002535343600000071
wherein η ═ x, y, ψ]TX represents the lateral displacement of the vessel, y represents the longitudinal displacement of the vessel, ψ represents the yaw angle;υ=[u,v,r]Tu represents the forward speed of the vessel, v represents the lateral speed of the vessel, and r represents the yaw rate of the vessel; j (ψ) represents a rotation matrix,
Figure BDA0002535343600000081
the tender vessel is typically operated at low speed, and therefore the dynamic model of vessel motion is summarized as:
Figure BDA0002535343600000082
in the formula (I), the compound is shown in the specification,
Figure BDA0002535343600000083
a matrix of the inertia is represented and,
Figure BDA0002535343600000084
a damping matrix is represented which,
Figure BDA0002535343600000085
representing the Coriolis and central centripetal force matrix, τT=[τuτvτr]TIndicating control force and moment, X(·),Y(·),N(·)Is a hydrodynamic parameter.
The leader ship serves as a reference of a following ship, and a mathematical model of the leader ship is defined as:
Figure BDA0002535343600000086
coordinate conversion is performed on equations (1) and (3) by the following equation (4):
Figure BDA0002535343600000087
offset (x) of displacemento,yo) In combination with equation (4) above, equation (4) can be converted to:
z1r=(xr-xo)cosψr+(yr-yo)sinψr
z2r=-(xr-xo)sinψr+(yr-yo)cosψr
z3r=ψr(5)
defining a synchronous control error system for the replenishment of a dynamically positioned vessel as:
zie=zi-ziri=1,2,3
ue=u-ur(6)
ve=v-vr
re=r-rr
the formula (6) is differentiated and linearized to obtain:
Figure BDA0002535343600000091
in the formula, xt=yt=[z1ez2ez3euevere]T
Figure BDA0002535343600000092
Wherein
Figure BDA0002535343600000093
d11=-Xu,d22=-Yv,d33=-Nr
According to the actual requirement, selecting proper sampling time, the equation (8) can be converted into discrete form:
x(k+1)=Ax(k)+Bτ(k) (9)
in the formula, x, y ∈ R6×1,A∈R6×6,B∈R6×3K denotes sampling time points, and A and B are discrete system parameter matrixes.
In the invention, the stability of the control system is ensured by designing a secondary terminal cost function, and in the terminal cost function module, a linear state feedback law is firstly designed for the control system:
τ=Kx (10)
in the formula, K represents a vector for controlling the gain. The state matrix P of the terminal cost is obtained by solving the lyapunov equation (11)):
(Ak+κI)TP+P(Ak+κI)=-(Qlq+KTRlqK) (11)
in the formula, Ak=A+BK,QlqAnd RlqA weight matrix representing the output and input, κ being selected by equation (12):
κ<-λmax(Ak) (12)
thereby ensuring AkThe real part of all eigenvalues of + kappa I is negative, where lambda ismax(Ak) Is represented by AkMaximum eigenvalue of the real part.
Further, in order to improve the robustness and the anti-interference capability of the system, a tender vessel control algorithm is designed based on a model prediction control algorithm. Defining a prediction time domain as N in the supply ship model prediction control modulepControl time domain as NcIn the prediction calculation process, let τ (k + N)c)=τ(k+Nc+1)=...=τ(k+Nc+Np) Then, the future state information of the system can be calculated sequentially by a certain prediction model (x (k +1) ═ ax (k) + B τ (k)):
Figure BDA0002535343600000101
where x (k + i | k) i is 1,2pIndicates future state information x (k + i) i, which is predicted from the state information at the sampling time point k, is 1,2p
According to the above formula (13), the following definitions are made:
Figure BDA0002535343600000102
rewriting the above formula (13) into a matrix form as shown in formula (15):
X=Fx(k)+Φ (15)
in the formula (I), the compound is shown in the specification,
Figure BDA0002535343600000103
then, the synchronous multivariable control problem of the supply ship is converted into an optimal control problem with terminal cost:
Figure BDA0002535343600000104
the constraint conditions are as follows:
min≤≤max(17)
in the formula, QmpcAnd RmpcIs an output weight matrix and an input weight matrix,
Figure BDA0002535343600000111
PXpis the terminal cost, XpRepresents the terminal state of X, and P is the state matrix of the terminal cost function determined by equation (11).
Furthermore, in order to solve the problems of large calculation amount and low calculation speed in the traditional model predictive control algorithm, in the neural dynamic optimization system module, the parallel calculation architecture of the recurrent neural network is utilized to accelerate the calculation speed, so that the proposed algorithm has timeliness. In the neural dynamics optimization system module, bringing formula (15) into formula (16) results in:
Figure BDA0002535343600000112
by definition
Figure BDA0002535343600000113
Equation (18) is simplified to:
J=TH+W+ξ (20)
the objective function (equation (20)) is solved iteratively by establishing a neural dynamic optimization system (equation (21)), so that the calculation speed is accelerated:
Figure BDA0002535343600000114
in the formula (I), the compound is shown in the specification,
Figure BDA0002535343600000115
is the convergence rate, and Ω () is taken to be ▽ J, PΩIs the projection operator, and the projection operator,
Figure BDA0002535343600000121
example 1
The embodiment shows the synchronous control performance of the dynamically positioned ship when the novel replenishment ship synchronous control algorithm with timeliness and stability provided by the invention is adopted to finish the tracking replenishment task from different starting points. Fig. 2 shows the tracking performance of a dynamic positioning vessel when the replenishment task is performed from a different starting point, fig. 3 shows the synchronization process of the speed, and fig. 4 shows the response of the control input. In this embodiment, it can be seen from fig. 1 that although the starting point is at a different starting position, the following ship can still automatically follow the leader ship for replenishment, and it can be seen from fig. 2 that the synchronization of the speeds is automatically achieved in a steady state.
In addition, in the embodiment, the calculation time is compared, the comparison result is shown in the following table, and according to the data in the table, the novel supply ship synchronous control algorithm with timeliness and stability provided by the invention can be obtained, so that the calculation time of the traditional MPC algorithm can be effectively reduced, and the timeliness of the designed algorithm can be improved.
TABLE 1 calculated time comparison
Figure BDA0002535343600000122
The above-mentioned serial numbers of the embodiments of the present invention are merely for description and do not represent the merits of the embodiments.
In the above embodiments of the present invention, the descriptions of the respective embodiments have respective emphasis, and for parts that are not described in detail in a certain embodiment, reference may be made to related descriptions of other embodiments.
In the embodiments provided in the present application, it should be understood that the disclosed technology can be implemented in other ways. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units may be a logical division, and in actual implementation, there may be another division, for example, multiple units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, units or modules, and may be in an electrical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.

Claims (4)

1. A tender ship synchronous control algorithm with timeliness and stability is characterized in that: the system comprises a terminal state matrix module, a model prediction control module and a neural dynamic optimization module;
for the tracking control problem of the supply ship, a three-degree-of-freedom ship motion model is adopted:
Figure FDA0002535343590000011
wherein η ═ x, y, ψ]TX represents the lateral displacement of the vessel, y represents the longitudinal displacement of the vessel, ψ represents the yaw angle; u, v, r]TU represents the forward speed of the vessel, v represents the lateral speed of the vessel, and r represents the yaw rate of the vessel; j (ψ) represents a rotation matrix,
Figure FDA0002535343590000012
the tender vessel is typically operated at low speed, and therefore the dynamic model of vessel motion is summarized as:
Figure FDA0002535343590000013
in the formula (I), the compound is shown in the specification,
Figure FDA0002535343590000014
a matrix of the inertia is represented and,
Figure FDA0002535343590000015
a damping matrix is represented which,
Figure FDA0002535343590000016
representing the Coriolis and central centripetal force matrix, τT=[τuτvτr]TRepresenting control force and moment; wherein X(·),Y(·),N(·)All represent hydrodynamic parameters;
the leader ship serves as a reference of a following ship, and a mathematical model of the leader ship is defined as:
Figure FDA0002535343590000017
Figure FDA0002535343590000018
the formula (1) is converted into by coordinate conversion of the formula (1) and the formula (3) by the following formula (4):
Figure FDA0002535343590000021
offset (x) of displacemento,yo) In combination with equation (4) above, equation (4) can be converted to:
z1r=(xr-xo)cosψr+(yr-yo)sinψr
z2r=-(xr-xo)sinψr+(yr-yo)cosψr
z3r=ψr; (5)
defining a synchronous control error system for the replenishment of a dynamically positioned vessel as:
Figure FDA0002535343590000022
the formula (6) is obtained by differentiation and linearization:
Figure FDA0002535343590000023
in the formula, xt=yt=[z1ez2ez3euevere]T
Figure FDA0002535343590000024
Wherein the content of the first and second substances,
Figure FDA0002535343590000025
d11=-Xu,d22=-Yv,d33=-Nr
by actually giving the appropriate sampling time, said equation (8) is converted into a discrete form:
x(k+1)=Ax(k)+Bτ(k); (9)
in the formula, x, y ∈ R6×1,A∈R6×6,B∈R6×3K denotes sampling time points, and A and B denote a discretized system parameter matrix.
2. The tender vessel synchronization control algorithm with timeliness and stability of claim 1, wherein:
the closed loop stability of the control system is ensured through a secondary terminal cost function, in the terminal cost function module, a linear state feedback law is firstly designed for the control system:
τ=Kx; (10)
wherein K represents a vector of control gains; the state matrix P of the terminal cost function is obtained by solving the Lyapunov equation, namely the formula (11):
(Ak+κI)TP+P(Ak+κI)=-(Qlq+KTRlqK); (11)
in the formula, Ak=A+BK,QlqAnd RlqA weight matrix representing the output and input, κ being selected by equation (12):
κ<-λmax(Ak); (12)
Akthe real part of all eigenvalues of + kappa I are negative values, where lambda ismax(Ak) Is represented by AkMaximum eigenvalue of the real part.
3. The tender vessel synchronization control algorithm with timeliness and stability of claim 1, wherein:
defining a prediction time domain as N in the supply ship model prediction control modulepControl time domain as NcIn the prediction calculation process, it is assumed that:
τ(k+Nc)=τ(k+Nc+1)=...=τ(k+Nc+Np);
the future state information of the system can be calculated sequentially by a certain prediction model (x (k +1) ═ ax (k) + B τ (k)):
Figure FDA0002535343590000031
where x (k + i | k) i is 1,2pIndicates future state information x (k + i) i, which is predicted from the state information at the sampling time point k, is 1,2p
According to the above formula (13), the following definitions are made:
Figure FDA0002535343590000041
rewriting the above formula (13) into a matrix form as shown in formula (15):
X=Fx(k)+Φ; (15)
in the formula (I), the compound is shown in the specification,
Figure FDA0002535343590000042
then, the synchronous multivariable control problem of the supply ship is converted into an optimal control problem with terminal cost:
Figure FDA0002535343590000043
the constraint conditions are as follows:
min≤≤max; (17)
in the formula, QmpcAnd RmpcIs an output weight matrix and an input weight matrix,
Figure FDA0002535343590000044
is the terminal cost, XpIs the terminal state of X and P is the state matrix of the terminal cost function determined by equation (11).
4. The tender vessel synchronization control algorithm with timeliness and stability of claim 1, wherein:
in the neural dynamic optimization system module, the proposed algorithm has timeliness through a parallel computing architecture of a recurrent neural network; in the neural dynamics optimization system module, bringing formula (15) into formula (16) results in:
Figure FDA0002535343590000051
by definition
Figure FDA0002535343590000052
Equation (18) is simplified to:
J=TH+W+ξ; (20)
the method comprises the following steps of (1) iteratively solving an objective function, namely the formula (20), by establishing a neural dynamic optimization system, namely the formula (21), so that the calculation speed is accelerated:
Figure FDA0002535343590000053
where, theta is the convergence speed,Ω()is taken as
Figure FDA0002535343590000054
PΩIs the projection operator, and the projection operator,
Figure FDA0002535343590000055
CN202010530817.5A 2020-06-11 2020-06-11 Supply ship synchronous control algorithm with timeliness and stability Active CN111665719B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010530817.5A CN111665719B (en) 2020-06-11 2020-06-11 Supply ship synchronous control algorithm with timeliness and stability

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010530817.5A CN111665719B (en) 2020-06-11 2020-06-11 Supply ship synchronous control algorithm with timeliness and stability

Publications (2)

Publication Number Publication Date
CN111665719A true CN111665719A (en) 2020-09-15
CN111665719B CN111665719B (en) 2022-11-22

Family

ID=72386574

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010530817.5A Active CN111665719B (en) 2020-06-11 2020-06-11 Supply ship synchronous control algorithm with timeliness and stability

Country Status (1)

Country Link
CN (1) CN111665719B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112068440A (en) * 2020-09-18 2020-12-11 江苏科技大学 Model prediction control-based AUV (autonomous Underwater vehicle) recovery butt joint power positioning control method

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107544500A (en) * 2017-09-18 2018-01-05 哈尔滨工程大学 A kind of unmanned boat berthing action trail planing method for considering constraint
CN108828955A (en) * 2018-08-16 2018-11-16 大连海事大学 Accurate Track In Track control method based on finite time extended state observer
CN109116731A (en) * 2018-07-27 2019-01-01 北京理工大学 A kind of data-driven method for optimally controlling switching linear system
CN109521797A (en) * 2018-12-28 2019-03-26 西北工业大学 The optimization mowing type formation control method of unmanned boat guidance submarine navigation device group
CN109656143A (en) * 2019-02-25 2019-04-19 曲阜师范大学 A kind of sea drives the adaptive tracking control method and system of ship entirely
CN110161853A (en) * 2019-05-09 2019-08-23 大连海事大学 A kind of novel ship craft integrated automated driving system with real-time
CN110879531A (en) * 2019-11-22 2020-03-13 北京理工大学重庆创新中心 Data-driven self-adaptive optimization control method and medium for random disturbance system

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107544500A (en) * 2017-09-18 2018-01-05 哈尔滨工程大学 A kind of unmanned boat berthing action trail planing method for considering constraint
CN109116731A (en) * 2018-07-27 2019-01-01 北京理工大学 A kind of data-driven method for optimally controlling switching linear system
CN108828955A (en) * 2018-08-16 2018-11-16 大连海事大学 Accurate Track In Track control method based on finite time extended state observer
CN109521797A (en) * 2018-12-28 2019-03-26 西北工业大学 The optimization mowing type formation control method of unmanned boat guidance submarine navigation device group
CN109656143A (en) * 2019-02-25 2019-04-19 曲阜师范大学 A kind of sea drives the adaptive tracking control method and system of ship entirely
CN110161853A (en) * 2019-05-09 2019-08-23 大连海事大学 A kind of novel ship craft integrated automated driving system with real-time
CN110879531A (en) * 2019-11-22 2020-03-13 北京理工大学重庆创新中心 Data-driven self-adaptive optimization control method and medium for random disturbance system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘程: "《船舶路径跟踪与减横摇综合控制研究》", 《万方学位论文》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112068440A (en) * 2020-09-18 2020-12-11 江苏科技大学 Model prediction control-based AUV (autonomous Underwater vehicle) recovery butt joint power positioning control method

Also Published As

Publication number Publication date
CN111665719B (en) 2022-11-22

Similar Documents

Publication Publication Date Title
Witkowska et al. Adaptive dynamic control allocation for dynamic positioning of marine vessel based on backstepping method and sequential quadratic programming
Tannuri et al. Dynamic positioning systems: An experimental analysis of sliding mode control
CN108319140B (en) Redefined output type model-free self-adaptive course control method and system
CN111007854B (en) Under-actuated ship trajectory tracking control system
Yang et al. A recurrent neural network based fuzzy sliding mode control for 4-DOF ROV movements
Zhang et al. Anti-disturbance control for dynamic positioning system of ships with disturbances
CN110162039A (en) A kind of novel integrated ship path trace and rollstabilization optimal control method
Shen et al. Dynamic surface control for tracking of unmanned surface vessel with prescribed performance and asymmetric time-varying full state constraints
CN112034858B (en) Model-free adaptive heading control method fusing weak observation high-order output data
CN110161853A (en) A kind of novel ship craft integrated automated driving system with real-time
CN111665719B (en) Supply ship synchronous control algorithm with timeliness and stability
CN115113524B (en) ASV multiport event trigger path tracking control method based on intervention LVS guidance
Meng et al. Distributed leaderless formation control for multiple autonomous underwater vehicles based on adaptive nonsingular terminal sliding mode
McNinch et al. Sliding mode setpoint control of an underactuated surface vessel: simulation and experiment
Fu et al. Formation control for water-jet USV based on bio-inspired method
Jia et al. Distributed dynamic rendezvous control of the AUV-USV joint system with practical disturbance compensations using model predictive control
Børhaug et al. Formation control of underactuated marine vehicles with communication constraints
Witkowska et al. Adaptive backstepping tracking control for an over–actuated DP marine vessel with inertia uncertainties
CN112578797A (en) Ship formation control method based on distributed guidance under switching topology
CN110109357B (en) Semi-global self-adaptive control method for non-standard type non-linear aircraft
CN116449829A (en) Ship formation system switching control method based on two-degree-of-freedom LQR
Zhu et al. Adaptive neural network fixed-time sliding mode control for trajectory tracking of underwater vehicle
Xiang et al. Coordinated 3D path following for autonomous underwater vehicles via classic PID controller
Sun et al. Distributed anti-windup NN-sliding mode formation control of multi-ships with minimum cost
Cui et al. Trajectory re-planning and tracking control of unmanned underwater vehicles on dynamic model

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
GR01 Patent grant
GR01 Patent grant