The Pontryagin maximum principle. Ab ovo usque ad mala

A proof of the Pontryagin maximum principle for a sufficiently general optimal control problem is presented; the proof is based on the implicit function theorem and the theorem on the solvability of a finite-dimensional system of nonlinear equations. The exposition is self-contained: all necessary preliminary facts are proved. These facts are mainly related to the properties of solutions to differential equations with discontinuous right-hand side and are derived as corollaries to the implicit function theorem, which, in turn, is a direct consequence of Newton’s method for solving nonlinear equations.     

常系数线性微分方程组的基解矩阵的一种新求法

在求解常系数线性微分方程组时,关键是基解矩阵的计算.给出了利用哈密顿—凯莱定理计算基解矩阵的一种方法,并通过实例说明了这种方法的特点和在简化计算方面的有效性.     

Distributed optimal control of the viscous Burgers equation via a Legendre pseudo‐spectral approach

This paper presents a computational technique based on the pseudo‐spectral method for the solution of distributed optimal control problem for the viscous Burgers equation. By using pseudo‐spectral method, the problem is converted to a classical optimal control problem governed by a system of ordinary differential equations, which can be solved by well‐developed direct or indirect methods. For solving the resulting optimal control problem, we present an indirect method by deriving and numerically solving the first‐order optimality conditions. Numerical tests involving both unconstrained and constrained control problems are considered. Copyright © 2015 John Wiley & Sons, Ltd.     

Rationalized Haar approach for nonlinear constrained optimal control problems

This paper presents a numerical method for solving nonlinear optimal control problems including state and control inequality constraints. The method is based upon rationalized Haar functions. The differential and integral expressions which arise in the system dynamics, the performance index and the boundary conditions are converted into some algebraic equations which can be solved for the unknown coefficients. Illustrative examples are included to demonstrate the validity and applicability of the technique.     

The Cubic B-Spline Method for a Class of Caputo-Fabrizio Fractional Differential Equations北大核心CSCD

基于分数阶微积分基本定理和三次B样条理论,构造了求解线性Caputo-Fabrizio型分数阶微分方程数值解的三次B样条方法,利用分数阶微积分基本定理将初值问题转化为关于解函数的表达式,再使用三次B样条函数逼近表达式中积分项的被积函数,进而计算了一类Caputo-Fabrizio型分数阶微分方程的数值解.给出了所构造的三次B样条方法的误差估计、收敛性和稳定性的理论证明.数值实验表明,该文数值方法在求解一类Caputo-Fabrizio型分数阶微分方程数值解时具有一定的可行性和有效性,且计算精度和计算效率优于现有的两种数值方法.     

一类微分方程组的积分解法及其在厚板热力弯曲中的应用

提出并证明了一类非齐次线性微分方程组的积分解法,并以求解受激光和横向力联合加载厚板的轴对称弯曲问题作为该方法的应用实例,籍此给出了考虑挤压效应、旋转惯性效应和剪切变形效应时厚板动态弯曲的扰度和转角公式。结果表明运用该定理求解问题具有规范、简明的特点。     

A novel feedforward–feedback suboptimal control of linear time-delay systems

This paper presents a new approach for solving the optimal control problem of linear time-delay systems with a quadratic cost functional. In this approach, a method of successive substitution is employed to convert the original time-delay optimal control problem into a sequence of linear time-invariant ordinary differential equations (ODEs) without delay and advance terms. The obtained optimal control consists of a linear state feedback term and a forward term. The feedback term is determined by solving a matrix Riccati differential equation. The forward term is an infinite sum of adjoint vectors, which can be obtained by solving recursively the above-mentioned sequence of linear non-delay ODEs. A fast-converging iterative algorithm for this purpose is presented which provides a promising possible reduction of computational efforts. Numerical examples demonstrating the efficiency, simplicity and high accuracy of the suggested technique have been included. Simulation results reveal that just a few iterations of the proposed algorithm are required to find an accurate enough feedforward–feedback suboptimal control.     

Existence of solutions and optimal control for reflecting stochastic differential equations with applications to population control theory *

For the d–dimensional reflecting stochastic differential equations (1) with non-smooth boundary and unbounded domain the existence of a strong solution, (weak solution) is obtained under the conditions that the coefficients are less than linear growth and they are non-Lipschitz, (and the diffusion coefficient is non-degenerate, the drift coefficient is bounded and measurable only). Moreover, the Girsanov theorem and the martingale representation theorem with respect to system (1) are also derived. Then by using the Ekeland lemma and the martingale method the existence, necessary and sufficient conditions for an optimal control and an optimal control are obtained. The results are then applied to solve an optimal control problem for a stochastic population model     

多目标规划在条件微分方程组中的应用

从工程实际出发,借助最佳逼近论和总体极值的思想,运用常微分方程组的求解理论,最优化理论与数值方法,为在最优控制中的一类条件微分方程组的求解,开辟了一条新的求解途径,并用多个计算实例,证明了算法的有效性和可行性。     

Global optimization of polynomial-expressed nonlinear optimal control problems with semidefinite programming relaxation

In this paper, we propose a new deterministic global optimization method for solving nonlinear optimal control problems in which the constraint conditions of differential equations and the performance index are expressed as polynomials of the state and control functions. The nonlinear optimal control problem is transformed into a relaxed optimal control problem with linear constraint conditions of differential equations, a linear performance index, and a matrix inequality condition with semidefinite programming relaxation. In the process of introducing the relaxed optimal control problem, we discuss the duality theory of optimal control problems, polynomial expression of the approximated value function, and sum-of-squares representation of a non-negative polynomial. By solving the relaxed optimal control problem, we can obtain the approximated global optimal solutions of the control and state functions based on the degree of relaxation. Finally, the proposed global optimization method is explained, and its efficacy is proved using an example of its application.     

Comparison of approximate symmetry and approximate homotopy symmetry to the Cahn–Hilliard equation

Complete infinite order approximate symmetry and approximate homotopy symmetry classifications of the Cahn–Hilliard equation are performed and the reductions are constructed by an optimal system of one-dimensional subalgebras. Zero order similarity reduced equations are nonlinear ordinary differential equations while higher order similarity solutions can be obtained by solving linear variable coefficient ordinary differential equations. The relationship between two methods for different order are studied and the results show that the approximate homotopy symmetry method is more effective to control the convergence of series solutions than the approximate symmetry one.     

与年龄相关的非线性种群系统的最优收获控制

讨论与年龄相关的一类非线性种群系统的最优收获控,运用泛函分析的Mazru's定理证明了最优收获控制的存在性,利用G-微分和Lions的变分不等式理论,导出了控制为最优的必要条件,得到了由积分-偏微分方程和变分不等式组成的最优性组,由最优性组确定最优控制.     

A Short Note on the Quasilinearization Method for Fractional Differential Equations

Recent literature shows that for certain classes of fractional differential equations the monotone iterative technique fails to guarantee the quadratic convergence of the quasilinearization method. The present work proves the quadratic convergence of the quasilinearization method and the existence and uniqueness of the solution of such a class of fractional differential equations. Our analysis depends upon the classical Kantorovich theorem on Newton's method. Various examples are discussed in order to illustrate our approach.     

An embedding theorem for integral equations

In this paper, an embedding theorem is established for a system of nonlinear integral equations of the Volterra type. The main result is basic in the development of a maximum principle for an optimal control problem in which the state variables are determined as solutions to integral equations.     

Extragradient method for solving an optimal control problem with implicitly specified boundary conditions

An optimal control problem formulated as a system of linear ordinary differential equations with boundary conditions implicitly specified as a solution to a finite-dimensional minimization problem is considered. An extragradient method for solving this problem is proposed, and its convergence is studied.     

First-order strong variation algorithm for optimal control problems involving hyperbolic systems

This paper considers an optimal control problem involving linear, hyperbolic partial differential equations. A first-order strong variational technique is used to obtain an algorithm for solving the optimal control problem iteratively. It is shown that the accumulation points of the sequence of controls generated by the algorithm (if they exist) satisfy a necessary condition for optimality.     

Optimal control of switched systems with time delay

In this note, we develop a computational method for solving an optimal control problem which is governed by a switched dynamical system with time delay. Our approach is to parameterize the switching instants as a new parameter vector to be optimized. Then, we derive the required gradient of the cost function which is obtained via solving a number of delay differential equations forward in time. On this basis, the optimal control problem can be solved as a mathematical programming problem.     

A semi-analytical direct optimal control solution for strongly excited and dissipative Hamiltonian systems

A semi-analytical direct optimal control solution for strongly excited and dissipative Hamiltonian systems is proposed based on the extended Hamiltonian principle, the Hamilton-Jacobi-Bellman (HJB) equation and its variational integral equation, and the finite time element approximation. The differential extended Hamiltonian equations for structural vibration systems are replaced by the variational integral equation, which can preserve intrinsic system structure. The optimal control law dependent on the value function is determined by the HJB equation so as to satisfy the overall optimality principle. The partial differential equation for the value function is converted into the integral equation with variational weighting. Then the successive solution of optimal control with system state is designed. The two variational integral equations are applied to sequential time elements and transformed into the algebraic equations by using the finite time element approximation. The direct optimal control on each time element is obtained respectively by solving the algebraic equations, which is unconstrained by the system state observed. The proposed control algorithm is applicable to linear and nonlinear systems with the quadratic performance index, and takes into account the effects of external excitations measured on control. Numerical examples are given to illustrate the optimal control effectiveness.     

On the Numerical Solution of Differential Games for Neutral-Type Linear Systems

The paper deals with a zero-sum differential game, in which the dynamics of a conflict-controlled system is described by linear functional differential equations of neutral type and the quality index is the sum of two terms: the first term evaluates the history of motion of the system realized up to the terminal time, and the second term is an integral–quadratic evaluation of the corresponding control realizations of the players. To calculate the value and construct optimal control laws in this differential game, we propose an approach based on solving a suitable auxiliary differential game, in which the motion of a conflict-controlled system is described by ordinary differential equations and the quality index evaluates the motion at the terminal time only. To find the value and the saddle point in the auxiliary differential game, we apply the so-called method of upper convex hulls, which leads to an effective solution in the case under consideration due to the specific structure of the quality index and the geometric constraints on the control actions of the players. The efficiency of the approach is illustrated by an example, and the results of numerical simulations are presented. The constructed optimal control laws are compared with the optimal control procedures with finitedimensional approximating guides, which were developed by the authors earlier.     

Solving some problems for systems of linear ordinary differential equations with redundant conditions

Numerical methods are proposed for solving some problems for a system of linear ordinary differential equations in which the basic conditions (which are generally nonlocal ones specified by a Stieltjes integral) are supplemented with redundant (possibly nonlocal) conditions. The system of equations is considered on a finite or infinite interval. The problem of solving the inhomogeneous system of equations and a nonlinear eigenvalue problem are considered. Additionally, the special case of a self-adjoint eigenvalue problem for a Hamiltonian system is addressed. In the general case, these problems have no solutions. A principle for constructing an auxiliary system that replaces the original one and is normally consistent with all specified conditions is proposed. For each problem, a numerical method for solving the corresponding auxiliary problem is described. The method is numerically stable if so is the constructed auxiliary problem.