Analysis of the factorization method for elliptic differential equations

This paper presents a numerical analysis of the method of factorization for elliptic boundary value problems. A second-order elliptic boundary value problem is transformed into a decoupled system of first-order initial value problems. This elaborates the decoupled system’s occupying operator, the solution of the problem, and the affine part of the operator. The solution of the problem is thus obtained by solving the factorized system using the finite-element method. A comparative study is made with the help of an example that illustrates many of its properties.     

A new computational method for solving two‐dimensional Stratonovich Volterra integral equation

This paper presents a method for computing numerical solutions of two‐dimensional Stratonovich Volterra integral equations using one‐dimensional modification of hat functions and two‐dimensional modification of hat functions. The problem is transformed to a linear system of algebraic equations using the operational matrix associated with one‐dimensional modification of hat functions and two‐dimensional modification of hat functions. The error analysis of the method is given. The method is computationally attractive, and applications are demonstrated by a numerical example. Copyright © 2017 John Wiley & Sons, Ltd.     

A numerically stable least squares solution to the quadratic programming problem

The strictly convex quadratic programming problem is transformed to the least distance problem — finding the solution of minimum norm to the system of linear inequalities. This problem is equivalent to the linear least squares problem on the positive orthant. It is solved using orthogonal transformations, which are memorized as products. Like in the revised simplex method, an auxiliary matrix is used for computations. Compared to the modified-simplex type methods, the presented dual algorithm QPLS requires less storage and solves ill-conditioned problems more precisely. The algorithm is illustrated by some difficult problems.      

有资格限制的指派问题的求解方法

在实际的指派工作中，常会遇到某个人有没有资格去承担某项工作的问题，因此，本建立了有资格限制的指派问题的数学模型。在此数学模型中，将效益矩阵转化为判定矩阵，由此给出了判定此种指派问题是否有解的方法；在有解的情况下，进一步将效益矩阵转化为求解矩阵，从而将有资格限制的指派问题化为传统的指派问题来求解。最后给出了一个数值例子来说明这样的处理方法是有效的。     

Modified Legendre–Gauss–Radau Collocation Method for Optimal Control Problems with Nonsmooth Solutions

A new method is developed for solving optimal control problems whose solutions are nonsmooth. The method developed in this paper employs a modified form of the Legendre–Gauss–Radau orthogonal direct collocation method. This modified Legendre–Gauss–Radau method adds two variables and two constraints at the end of a mesh interval when compared with a previously developed standard Legendre–Gauss–Radau collocation method. The two additional variables are the time at the interface between two mesh intervals and the control at the end of each mesh interval. The two additional constraints are a collocation condition for those differential equations that depend upon the control and an inequality constraint on the control at the endpoint of each mesh interval. The additional constraints modify the search space of the nonlinear programming problem such that an accurate approximation to the location of the nonsmoothness is obtained. The transformed adjoint system of the modified Legendre–Gauss–Radau method is then developed. Using this transformed adjoint system, a method is developed to transform the Lagrange multipliers of the nonlinear programming problem to the costate of the optimal control problem. Furthermore, it is shown that the costate estimate satisfies one of the Weierstrass–Erdmann optimality conditions. Finally, the method developed in this paper is demonstrated on an example whose solution is nonsmooth.

     

Distributionally robust parameter identification of a time-delay dynamical system with stochastic measurements

In this paper, we consider a parameter identification problem involving a time-delay dynamical system, in which the measured data are stochastic variable. However, the probability distribution of this stochastic variable is not available and the only information we have is its first moment. This problem is formulated as a distributionally robust parameter identification problem governed by a time-delay dynamical system. Using duality theory of linear optimization in a probability space, the distributionally robust parameter identification problem, which is a bi-level optimization problem, is transformed into a single-level optimization problem with a semi-infinite constraint. By applying problem transformation and smoothing techniques, the semi-infinite constraint is approximated by a smooth constraint and the convergence of the smooth approximation method is established. Then, the gradients of the cost and constraint functions with respect to time-delay and parameters are derived. On this basis, a gradient-based optimization method for solving the transformed problem is developed. Finally, we present an example, arising in practical fermentation process, to illustrate the applicability of the proposed method.     

Nonmonotone filter DQMM method for the system of nonlinear equations

In this paper, we propose a nonmonotone filter Diagonalized Quasi-Newton Multiplier (DQMM) method for solving system of nonlinear equations. The system of nonlinear equations is transformed into a constrained nonlinear programming problem which is then solved by nonmonotone filter DQMM method. A nonmonotone criterion is used to speed up the convergence progress in some ill-conditioned cases. Under reasonable conditions, we give the global convergence properties. The numerical experiments are reported to show the effectiveness of the proposed algorithm.     

时滞广义时变系统的容许性分析与镇定

研究了时滞广义时变系统的容许性与镇定性问题.首先,基于广义Lyapunov不等式、线性矩阵不等式和受限等价方法,建立时滞广义时变系统的Lyapunov不等式,将时滞广义时变系统的容许性问题转化为求解时滞广义时变系统的Lyapunov不等式问题,得到了系统容许的充分条件.然后,根据充分条件进一步研究了时滞广义时变系统的镇定问题,给出了状态反馈镇定器的设计方法.最后,通过数值算例验证了所得结论的有效性.     

A Generalization of Ritz-Variational Method for Solving a Class of Fractional Optimization Problems

This paper presents an approximate method for solving a class of fractional optimization problems with multiple dependent variables with multi-order fractional derivatives and a group of boundary conditions. The fractional derivatives are in the Caputo sense. In the presented method, first, the given optimization problem is transformed into an equivalent variational equality; then, by applying a special form of polynomial basis functions and approximations, the variational equality is reduced to a simple linear system of algebraic equations. It is demonstrated that the derived linear system has a unique solution. We get an approximate solution for the initial optimization problem by solving the final linear system of equations. The choice of polynomial basis functions provides a method with such flexibility that all initial and boundary conditions of the problem can be easily imposed. We extensively discuss the convergence of the method and, finally, present illustrative test examples to demonstrate the validity and applicability of the new technique.     

Conjugate gradient techniques for the optimal control evolution dam problem

This paper considers the numerical simulation of optimal control evolution dam problem by using conjugate gradient method.The paper considers the free boundary value problem related to time dependent fluid flow in a homogeneous earth rectangular dam.The dam is taken to be sufficiently long that the flow is considered to be two dimensional.On the left and right walls of the dam there is a reservoir of fluid at a level dependent on time.This problem can be transformed into a variational inequality on a fixed domain.The numerical techniques we use are based on a linear finite element method to approximate the state equations and a conjugate gradient algorithm to solve the discrete optimal control problem.This algorithm is based on Armijo's rule in the unconstrained optimization theory.The convergence of the discrete optimal solutions to the continuous optimal solutions,and the convergence of the conjugate gradient algorithm are proved.A numerical example is given to determine the location of the minimum surface     

Symbolic-numeric efficient solution of optimal control problems for multibody systems

This paper presents an efficient symbolic-numerical approach for generating and solving the boundary value problem-differential algebraic equation (BVP-DAE) originating from the variational form of the optimal control problem (OCP). This paper presents the method for the symbolic derivation, by means of symbolic manipulation software (Maple), of the equations of the OCP applied to a generic multibody system. The constrained problem is transformed into a nonconstrained problem, by means of the Lagrange multipliers and penalty functions. From the first variation of the nonconstrained problem a BVP-DAE is obtained, and the finite difference discretization yields a nonlinear systems. For the numerical solution of the nonlinear system a damped Newton scheme is used. The sparse and structured Jacobians is quickly inverted by exploiting the sparsity pattern in the solution strategy. The proposed method is implemented in an object oriented fashion, and coded in C++ language. Efficiency is ensured in core routines by using Lapack and Blas for linear algebra.     

论理想流体与线弹性结构的耦联振动

本文对理想流体与线弹性结构的耦联振动问题作理论分析和数值分析.文中证明了耦联振动的固有频率存在并且均为正实数.将流-固耦合系统分析转化为单一结构物在真空中的自由振动分析后,频率方程中不再含有流体变元.使问题得以大大简化.给出了数值解的收敛性证明,以保证解的可靠性.文中还综合里兹法、边界元和有限元方法,提出一种分析转化后结构的混合算法.利用该算法,只需对现有结构分析程序稍作改进,就可分析那些理想流场与结构的耦合问题.一些数例说明了算法的有效性.     

Existence of solutions to Hamilton-Jacobi functional differential equations

This paper deals with the Cauchy problem for nonlinear first order partial functional differential equations. The unknown function is the functional variable in the equation and the partial derivatives appear in a classical sense. A theorem on the local existence of a generalized solution is proved. The initial problem is transformed into a system of functional integral equations for an unknown function and for their partial derivatives with respect to spatial variables. The existence of solutions of this system is proved by using a method of successive approximations. A method of bicharacteristics and integral inequalities are applied.     

An optimal control problem in the study of liver kinetics

The enzyme activities in the liver are described by a system of ordinary differential equations. A certain substance in the blood is transformed twice by different kinds of enzymes. The mathematical problem is to determine the distribution of the enzymes along the blood flow so that the outflow concentration of the once-transformed form of the substance is as small as possible. This problem has been considered before and solved for particular types of enzyme kinetics. In this paper, we solve the problem for more general types of kinetics (including substrate- inhibition kinetics). The methods used are also different, in that the problem is considered as a problem of optimal control and Pontryagin's maximum principle is applied to derive necessary conditions.     

A new method for solving nonlinear second order partial differential equations

In this paper, a new method for finding the approximate solution of a second order nonlinear partial differential equation is introduced. In this method the problem is transformed to an equivalent optimization problem. Then, by considering it as a distributed parameter control system, the theory of measure is used for obtaining the approximate solution of the original problem.     

On an augmented Lagrangian-based preconditioning of Oseen type problems

The paper deals with a general framework for constructing preconditioners for saddle point matrices, in particular as arising in the discrete linearized Navier-Stokes equations (Oseen’s problem). We utilize the so-called augmented Lagrangian framework, where the original linear system of equations is first transformed to an equivalent one, which latter is then solved by a preconditioned iterative solution method.     

美式看跌期权定价的一种混合数值方法

本文对美式看跌期权的定价提供了一种新的混合数值方法,即快速傅里叶变换法加龙格-库塔法.首先将美式看跌期权价格所满足的Black-Scholes微分方程定解问题转化为一个标准的抛物型初、边值问题,然后通过傅里叶变换,使之转换为一个不带股价变量的常微分方程初值问题,再利用龙格-库塔法对其进行数值求解.数值实验表明,本文算法是一种快速的高精度的算法.     

Stackelberg solutions for fuzzy random two-level linear programming through level sets and fractile criterion optimization

This paper considers Stackelberg solutions for two-level linear programming problems under fuzzy random environments. To deal with the formulated fuzzy random two-level linear programming problem, an α-stochastic two-level linear programming problem is defined through the introduction of α-level sets of fuzzy random variables. Taking into account vagueness of judgments of decision makers, fuzzy goals are introduced and the α-stochastic two-level linear programming problem is transformed into the problem to maximize the satisfaction degree for each fuzzy goal. Through fractile criterion optimization in stochastic programming, the transformed stochastic two-level programming problem can be reduced to a deterministic two-level programming problem. An extended concept of Stackelberg solution is introduced and a numerical example is provided to illustrate the proposed method.     

Classical solutions of hyperbolic functional differential systems

The paper deals with a generalized Cauchy problem for quasi-linear hyperbolic functional differential systems. The unknown function is the functional variable in the system of equations and the partial derivatives appear in the classical sense. A theorem on the local existence of a solution is proved. The initial problem is transformed into a system of functional integral equations for an unknown function and for their partial derivatives with respect to spatial variables. A method of bicharacteristics and integral inequalites are applied.     

Anomalous exponents of self-similar blow-up solutions to an aggregation equation in odd dimensions

We calculate the scaling behavior of the second-kind self-similar blow-up solution of an aggregation equation in odd dimensions. This solution describes the radially symmetric finite-time blowup phenomena and has been observed in numerical simulations of the dynamic problem. The nonlocal equation for the self-similar profile is transformed into a system of ODEs that is solved using a shooting method. The anomalous exponents are then retrieved from this transformed system.