动力学系统参数辨识问题最优控制解的理论与方法(Ⅰ)--基本概念及确定性系统参数辨识

本文基于动力学系统参数辨识问题最优控制解的概念和确定性动力学系统的最优控制理论,建立了参数辨识研究与最优控制理论的对应关系.将最优控制的数学理论和算法应用于参数辨识问题的研究.依据Hamilton-Jacobi-Bellman (HJB)方程解的理论阐述了动力学系统参数辨识最优控制解的存在唯一性问题,并据此得到了解决确定性系统参数辨识问题的具体算法步骤.     

带边界反馈弦方程的参数辨识——边界控制方法

分布参数系统中的参数辨识问题是一个极具挑战的数学问题.本文研究了带边界反馈弦方程的弹性模量的辨识问题,基于边界控制方法建立了重构未知参数的辨识算法,并且通过数值例子表明了辨识算法的有效性.     

动力学系统参数辨识问题最优控制解的理论与方法（I）—基本概 …

本文基于动力学系统参数辨识问题最优控制解的概念和确定性动力学系统的最优控制理论,建立了参数辨识研究与最优控制理论的对应关系。将最优控制的数学理论和算法应用于参数辨识问题的研究。依据Ｈａｍｉｌｔｏｎ－Ｊａｃｏｂｉ－Ｂｅｌｌｍａｎ（ＨＪＢ）方程解的理论阐述了动力学参数辨识最优控制解的存在唯一性问题,并据此得到了解决确定性系统参数辨识问题的具体算法步骤。     

一类非光滑分布参数系统的可辨识性及最优性条件

变压器温度场参数辨识问题是一种分片光滑的分布参数辨识问题,以流速为辨识参数,针对传质传热的一类分布参数系统参数辨识问题,证明了系统最优参数的存在性和控制参数为最优的必要条件,为变压器温度场的数值模拟研究提供了理论基础.     

动力学系统参数辨识问题最优控制解的理论与方法(Ⅱ)——随机系统参数辨识及应用实例

基于文(Ⅰ)(《应用数学和力学》,1998,20(2))的内容和随机最优控制理论,本文首先介绍了随机动力学系统参数辨识问题最优控制解的概念.然后讨论了建立参数辨识问题HJB方程的过程以及参数辨识的算法.最后给出了一个应用实例:解决动力学系统局部非线性参数辨识问题的方法.     

Theory and Algorithm of Optimal Control Solution to Dynamic Systme Parameters Identification (

基于文(Ⅰ)(《应用数学和力学》,1998,20(2))的内容和随机最优控制理论,本文首先介绍了随机动力学系统参数辨识问题最优控制解的概念。然后讨论了建立参数辨识问题HJB方程的过程以及参数辨识的算法。最后给出了一个应用实例解决动力学系统局部非线性参数辨识问题的方法。     

径流过程的广义Hammersstein模型及其辨识

本文针对一类水文随机模拟系统的参数辨识问题,提出了一种广义 Hammerstein辨识方法,讨论了广义Hammerstein函数的构造问题。文中给出了两个日径流系统辨识的实际例子。     

动力学系统参数辨识问题最优控制解的理论与方法（II）—随机？…

基于（Ｉ）（《应用数学和力学》,１９９８,２０（２））的内容和随机最优控制理论,本首先介绍了随机动力学系统参数辨识问题最优控制解的概念。然后讨论了建立参数辨识问题ＨＪＢ方程的过程以及参数辨识的算法,最后给出了一个应用实例：解决动力学系统局部非线性参数辨识问题的方法。     

几类典型随机非线性系统的辨识

考察实际中常见的三类典型随机非线性系统(即Wiener、Hammerstein和NARX系统)的辨识,首先概述了现有的递推和非递推辨识算法,然后介绍这三类系统的一个统一辨识框架:利用系统所确定的过程的马氏性及混合型,将辨识转化为求函数零点的问题,基于扩张截尾的随机逼近算法,得到了递推、强一致的辨识结果,并给出了数值模拟验证辨识算法收敛到真值.     

发汗冷却控制系统的参数辨识问题

研究发汗冷却控制系统中气动加热热流密度的参数辨识问题．证明了该参数辨识的存在及唯一性,给出了参数辨识所满足的充分必要条件,最后,根据得到的充分必要条件,尝试直接构造极小化序列,进而给出该系统参数辨识的算法．     

On bounds for the mode and median of the generalized hyperbolic and related distributions

Except for certain parameter values, a closed form formula for the mode of the generalized hyperbolic (GH) distribution is not available. In this paper, we exploit results from the literature on modified Bessel functions and their ratios to obtain simple but tight two-sided inequalities for the mode of the GH distribution for general parameter values. As a special case, we deduce tight two-sided inequalities for the mode of the variance-gamma (VG) distribution, and through a similar approach we also obtain tight two-sided inequalities for the mode of the McKay Type I distribution. The analogous problem for the median is more challenging, but we conjecture some monotonicity results for the median of the VG and McKay Type I distributions, from we which we conjecture some tight two-sided inequalities for their medians. Numerical experiments support these conjectures and also lead us to a conjectured tight lower bound for the median of the GH distribution.     

Half Space and Quarter Space Dirichlet Problems for the Partial Differential Equation

A customary, heuristic, method, by which the Poisson integral formula for the Dirichlet problem, for the half space, for Laplace's equation is obtained, involves Green's function, and Kelvin's method of images. Although this heuristic method leads one to guess the correct result, this Poisson formula still has to be verified directly, independently of the method by which it was arrived at, in order to be absolutely certain that a solution of the Dirichlet problem for the half space, for Laplace's equation, has been actually obtained. A similar heuristic method, as seems to be generally known, could be followed in solving the Dirichlet problem, for the half space, for the equation where is a real constant. However, in Part 1, a different, labor-saving, method is used to study Dirichlet problems for the equation. This method is essentially based on what Hadamard called the method of descent. Indeed, it is shown that he who has solved the half space Dirichlet problem for Laplace's equation has already solved the half space Dirichlet problem for the equation In Part 2, the solution formula for the quarter space Dirichlet problem for Laplace's equation is obtained from the Poisson integral formula for the half space Dirichlet problem for Laplace's equation. A representation theorem for harmonic functions in the quarter space is deduced. The method of descent is used, in Part 3, to obtain the solution formula for the quarter space Dirichlet problem for the equation by means of the solution formula for the quarter space Dirichlet problem for Laplace's equation. So that, indeed, it is also shown that he who has solved the quarter space Dirichlet problem for Laplace's equation has already solved the quarter space Dirichlet problem for the " equation" For the sake of completeness and clarity, and for the convenience of the reader, the appendix, at the end of Part 3, contains a detailed proof that the Poisson integral formula solves the half space Dirichlet problem for Laplace's equation. The Bibliography for Parts 1,2, 3 is to be found at the end of Part 1.     

An approximation of stochastic hyperbolic equations: case with Wiener process

In the present paper, the two‐step difference scheme for the Cauchy problem for the stochastic hyperbolic equation is presented. The convergence estimate for the solution of the difference scheme is established. In applications, the convergence estimates for the solution of difference schemes for the numerical solution of four problems for hyperbolic equations are obtained. The theoretical statements for the solution of this difference scheme are supported by the results of the numerical experiment. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimates for norms of matrix-valued and operator-valued functions and some of their applications

A survey is presented of estimates for a norm of matrix-valued and operator-valued functions obtained by the author. These estimates improve the Gel'fand-Shilov estimate for regular functions of matrices and Carleman's estimates for resolvents of matrices and compact operators.From the estimates for resolvents, the well-known result for spectrum perturbations of self-adjoint operators is extended to quasi-Hermitian operators. In addition, the classical Schur and Brown's inequalities for eigenvalues of matrices are improved.From estimates for the exponential function (semigroups), bounds for solution norms of nonlinear differential equations are derived. These bounds give the stability criteria which make it possible to avoid the construction of Lyapunov functions in appropriate situations.     

Superconvergence of mixed finite element semi-discretizations of two time-dependent problems

We will show that some of the superconvergence properties for the mixed finite element method for elliptic problems are preserved in the mixed semi-discretizations for a diffusion equation and for a Maxwell equation in two space dimensions. With the help of mixed elliptic projection we will present estimates global and pointwise in time. The results for the Maxwell equations form an extension of existing results. For both problems, our results imply that post-processing and a posteriori error estimation for the error in the space discretization can be performed in the same way as for the underlying elliptic problem.     

High-order computational methods for option valuation under multifactor models

Many of the different numerical techniques in the partial differential equations framework for solving option pricing problems have employed only standard second-order discretization schemes. A higher-order discretization has the advantage of producing low size matrix systems for computing sufficiently accurate option prices and this paper proposes new computational schemes yielding high-order convergence rates for the solution of multi-factor option problems. These new schemes employ Galerkin finite element discretizations with quadratic basis functions for the approximation of the spatial derivatives in the pricing equations for stochastic volatility and two-asset option problems and time integration of the resulting semi-discrete systems requires the computation of a single matrix exponential. The computations indicate that this combination of high-order finite elements and exponential time integration leads to efficient algorithms for multi-factor problems. Highly accurate European prices are obtained with relatively coarse meshes and high-order convergence rates are also observed for options with the American early exercise feature. Various numerical examples are provided for illustrating the accuracy of the option prices for Heston’s and Bates stochastic volatility models and for two-asset problems under Merton’s jump-diffusion model.     

Nonlocal boundary-value problems for abstract parabolic equations: well-posedness in Bochner spaces

The well-posedness of the nonlocal boundary-value problem for abstract parabolic differential equations in Bochner spaces is established. The first and second order of accuracy difference schemes for the approximate solutions of this problem are considered. The coercive inequalities for the solutions of these difference schemes are established. In applications, the almost coercive stability and coercive stability estimates for the solutions of difference schemes for the approximate solutions of the nonlocal boundary-value problem for parabolic equation are obtained.     

Temperature prediction at critical points in district heating systems

Current methodologies for the optimal operation of district heating systems use model predictive control. Accurate forecasting of the water temperature at critical points is crucial for meeting constraints related to consumers while minimizing the production costs for the heat supplier. A new forecasting methodology based on conditional finite impulse response (cFIR) models is introduced, for which model coefficients are replaced by coefficient functions of the water flux at the supply point and of the time of day, allowing for nonlinear variations of the time delays. Appropriate estimation methods for both are described. Results are given for the test case of the Roskilde district heating system over a period of more than 6 years. The advantages of the proposed forecasting methodology in terms of a higher forecast accuracy, its use for simulation purposes, or alternatively for better understanding transfer functions of district heating systems, are clearly shown.     

Determinants of matrices over semirings

In this paper, the concept of determinants for the matrices over a commutative semiring is introduced, and a development of determinantal identities is presented. This includes a generalization of the Laplace and Binet–Cauchy Theorems, as well as on adjoint matrices. Also, the determinants and the adjoint matrices over a commutative difference-ordered semiring are discussed and some inequalities for the determinants and for the adjoint matrices are obtained. The main results in this paper generalize the corresponding results for matrices over commutative rings, for fuzzy matrices, for lattice matrices and for incline matrices.     

Stochastic service systems,random interval graphs and search algorithms

We consider several stochastic service systems, and study the asymptotic behavior of the moments of various quantities that have application to models for random interval graphs and algorithms for searching for an idle server or for an vacant or occupied waiting station. In some cases the moments turn out to involve Lambert series for the generating functions for the sums of powers of divisors of positive integers. For these cases we are able to obtain complete asymptotic expansions for the moments of the quantities in question. © 2013 Wiley Periodicals, Inc. Random Struct. Alg., 45, 421–442, 2014