Methodology for the solutions of some reduced Fokker‐Planck equations in high dimensions

In this paper, a new methodology is formulated for solving the reduced Fokker‐Planck (FP) equations in high dimensions based on the idea that the state space of large‐scale nonlinear stochastic dynamic system is split into two subspaces. The FP equation relevant to the nonlinear stochastic dynamic system is then integrated over one of the subspaces. The FP equation for the joint probability density function of the state variables in another subspace is formulated with some techniques. Therefore, the FP equation in high‐dimensional state space is reduced to some FP equations in low‐dimensional state spaces, which are solvable with exponential polynomial closure method. Numerical results are presented and compared with the results from Monte Carlo simulation and those from equivalent linearization to show the effectiveness of the presented solution procedure. It attempts to provide an analytical tool for the probabilistic solutions of the nonlinear stochastic dynamics systems arising from statistical mechanics and other areas of science and engineering.     

Stochastic period-doubling bifurcation analysis of stochastic Bonhoeffer--van der Pol system

In this paper, the Chebyshev polynomial approximation is applied to the problem of stochastic period-doubling bifurcation of a stochastic Bonhoeffer--van der Pol (BVP for short) system with a bounded random parameter. In the analysis, the stochastic BVP system is transformed by the Chebyshev polynomial approximation into an equivalent deterministic system, whose response can be readily obtained by conventional numerical methods. In this way we have explored plenty of stochastic period-doubling bifurcation phenomena of the stochastic BVP system. The numerical simulations show that the behaviour of the stochastic period-doubling bifurcation in the stochastic BVP system is by and large similar to that in the deterministic mean-parameter BVP system, but there are still some featured differences between them. For example, in the stochastic dynamic system the period-doubling bifurcation point diffuses into a critical interval and the location of the critical interval shifts with the variation of intensity of the random parameter. The obtained results show that Chebyshev polynomial approximation is an effective approach to dynamical problems in some typical nonlinear systems with a bounded random parameter of an arch-like probability density function.     

Intrusive Galerkin methods with upwinding for uncertain nonlinear hyperbolic systems

This paper deals with stochastic spectral methods for uncertainty propagation and quantification in nonlinear hyperbolic systems of conservation laws. We consider problems with parametric uncertainty in initial conditions and model coefficients, whose solutions exhibit discontinuities in the spatial as well as in the stochastic variables. The stochastic spectral method relies on multi-resolution schemes where the stochastic domain is discretized using tensor-product stochastic elements supporting local polynomial bases. A Galerkin projection is used to derive a system of deterministic equations for the stochastic modes of the solution. Hyperbolicity of the resulting Galerkin system is analyzed. A finite volume scheme with a Roe-type solver is used for discretization of the spatial and time variables. An original technique is introduced for the fast evaluation of approximate upwind matrices, which is particularly well adapted to local polynomial bases. Efficiency and robustness of the overall method are assessed on the Burgers and Euler equations with shocks.     

The Applications of Order Reduction Methods in Nonlinear Dynamic Systems

Two different order reduction methods of the deterministic and stochastic systems are discussed in this paper. First, the transient proper orthogonal decomposition (T-POD) method is introduced based on the high-dimensional nonlinear dynamic system. The optimal order reduction conditions of the T-POD method are provided by analyzing the rotor-bearing system with pedestal looseness fault at both ends. The efficiency of the T-POD method is verified via comparing with the results of the original system. Second, the polynomial dimensional decomposition (PDD) method is applied to the 2 DOFs spring system considering the uncertain stiffness to study the amplitude-frequency response. The numerical results obtained by the PDD method agree well with the Monte Carlo simulation (MCS) method. The results of the PDD method can approximate to MCS better with the increasing of the polynomial order. Meanwhile, the Uniform-Legendre polynomials can eliminate perturbation of the PDD method to a certain extent via comparing it with the Gaussian-Hermite polynomials.     

Polynomial chaos expansions for optimal control of nonlinear random oscillators

The polynomial chaos decomposition of stochastic variables and processes is implemented in conjunction with optimal polynomial control of nonlinear dynamical systems. The procedure is demonstrated on a base-excited system whereby ground motion is modeled as a stochastic process with a specified correlation function and is approximated by its Karhunen-Loeve expansion. An adaptive scheme for stochastic approximation with polynomial chaos bases is proposed which is based on a displacement-velocity norm and is applied to the identification of phase orbits of nonlinear oscillators. This approximation is then integrated in the design of an optimal polynomial controller, allowing for the efficient estimation of statistics and probability density functions of quantities of interest. Numerical investigations are carried out employing the polynomial chaos expansions and the Lyapunov asymptotic stability condition based control policy. The results reveal that the performance, as gaged by probabilistic quantities of interest, of the controlled oscillators is greatly improved. A comparative study is also presented against the classical stochastic optimal control, whereby statistical linearization based LQG is employed to design the optimal controller. It is remarked that the proposed polynomial chaos expansion is a preferred approach to the optimal control of nonlinear random oscillators.     

随机参数Duffing系统中的随机混沌及其延迟反馈控制

研究了谐和激励下含有界随机参数Duffing系统(简称随机Duffing系统)中的随机混沌及其延迟反馈控制问题.借助Gegenbauer多项式逼近理论，将随机Duffing系统转化为与其等效的确定性非线性系统.这样，随机Duffing系统在谐和激励下的混沌响应及其控制问题就可借等效的确定性非线性系统来研究.分析阐明了随机混沌的主要特点，并采用Wolf算法计算等效确定性非线性系统的最大Lyapunov指数，以判别随机Duffing系统的动力学行为.数值计算表明，恰当选取不同的反馈强度和延迟时间，可分别达到抑制或诱发系统混沌的目的，说明延迟反馈技术对随机混沌控制也是十分有效的. 关键词： 随机Duffing系统 延迟反馈控制 随机混沌 Gegenbauer多项式     

基于Chebyshev多项式逼近的随机 van der Pol系统的倍周期分岔分析

应用 Chebyshev 多项式逼近法研究了谐和激励作用下具有随机参数的随机van der Pol系统 的倍周期分岔现象.随机系统首先被转化成等价的确定性系统，然后通过数值方法求得响应 ，借此探索了随机van der Pol系统丰富的随机倍周期分岔现象.数值模拟显示随机van der Pol 系统存在与确定性系统极为相似的倍周期分岔行为，但受随机因素的影响，又有与之不 同之处.数值结果表明，Chebyshev 多项式逼近是研究非线性系统动力学问题的一种新的有 效方法. 关键词： Chebyshev 多项式 随机van der Pol 系统 倍周期分岔     

基于参数展开的同伦分析法在强非线性随机动力系统中的应用

将基于参数展开的同伦分析法(PE-HAM)进行了推广，使之适用于谐和激励与随机噪声联合作用下的强非线性随机动力系统. 通过构造合适的同伦映射，将对强非线性随机动力系统响应的求解转化为对一组线性随机微分方程的求解. 进一步研究了受到谐和与Gauss白噪声激励的强非线性Duffing振子，由PE-HAM得到了该系统的解过程和稳态概率密度的解析表达式. 数值模拟的结果说明了PE-HAM方法的精确性. 关键词： PE-HAM方法 强非线性随机动力系统 稳态概率密度 解过程 随机激励     

基于Laguerre多项式逼近法的随机双势阱Duffing系统的分岔和混沌研究

应用Laguerre正交多项式逼近法研究了含有随机参数的双势阱Duffing系统的分岔和混沌行为.系统参数为指数分布随机变量的非线性动力系统首先被转化为等价的确定性扩阶系统，然后通过数值方法求得其响应.数值模拟结果的比较表明，含有随机参数的双势阱Duffing系统保持着与确定性系统相类似的倍周期分岔和混沌行为，但是由于随机因素的影响，在局部小区域内随机参数系统的动力学行为会发生突变. 关键词： 双势阱Duffing系统 指数分布概率密度函数 Laguerre多项式逼近 随机分岔     

The underlying linear dynamics of some positive polynomial systems

The conditions of structural dynamical similarity of two special classes of positive polynomial nonlinear systems, the class of quasi-polynomial systems (Brenig (1988) [1]) and that of reaction kinetic networks with mass action kinetics (Horn and Jackson (1972) [2]) are investigated in this Letter. It is shown that both system classes have an underlying reduced linear dynamics. By applying the theory of X-factorable systems (Samardzija et al. (1989) [9]), it can be shown that the reduced linear dynamics is qualitatively similar to the original one within the positive orthant when the original nonlinear system has a unique positive equilibrium point.     

Inversion of Robin coefficient by a spectral stochastic finite element approach

This paper investigates a variational approach to the nonlinear stochastic inverse problem of probabilistically calibrating the Robin coefficient from boundary measurements for the steady-state heat conduction. The problem is formulated into an optimization problem, and mathematical properties relevant to its numerical computations are investigated. The spectral stochastic finite element method using polynomial chaos is utilized for the discretization of the optimization problem, and its convergence is analyzed. The nonlinear conjugate gradient method is derived for the optimization system. Numerical results for several two-dimensional problems are presented to illustrate the accuracy and efficiency of the stochastic finite element method.     

随机系统的概率密度函数形状调节

针对受高斯白噪声激励的非线性随机系统,提出了使状态响应的概率密度函数形状跟踪期望形状的调节方法.首先,确立了非线性随机系统的多项式反馈机制,同时对系统中的非线性部分进行多项式展开;然后,以Fokker-Planck-Kolmogorov方程为工具,导出了与控制增益相关的各阶矩递推方程,并根据跟踪问题的要求,构造了矩逼近优化问题,用梯度搜索法求解该优化问题,获得了调节函数;再依据特征函数与概率密度函数构成Fourier对的关系,对状态响应的概率密度函数进行重构;最后,通过两个例子仿真,验证了本文方法的有效性.     

Analysis of stochastic bifurcation and chaos in stochastic Duffing--van der Pol system via Chebyshev polynomial approximation

The Chebyshev polynomial approximation is applied to investigate the stochastic period-doubling bifurcation and chaos problems of a stochastic Duffing--van der Pol system with bounded random parameter of exponential probability density function subjected to a harmonic excitation. Firstly the stochastic system is reduced into its equivalent deterministic one, and then the responses of stochastic system can be obtained by numerical methods. Nonlinear dynamical behaviour related to stochastic period-doubling bifurcation and chaos in the stochastic system is explored. Numerical simulations show that similar to its counterpart in deterministic nonlinear system of stochastic period-doubling bifurcation and chaos may occur in the stochastic Duffing--van der Pol system even for weak intensity of random parameter. Simply increasing the intensity of the random parameter may result in the period-doubling bifurcation which is absent from the deterministic system.     

永磁同步风力发电机随机分岔现象的全局分析

永磁同步风力发电机在运行过程中不可避免地会受到风能的随机干扰,本文建立了在输入机械转矩存在随机干扰情况下永磁同步风力发电机的数学模型,采用胞映射方法分析了随机干扰强度变化时系统全局结构的演化行为,并通过数值模拟对理论分析进行验证.研究结果表明,随着随机干扰强度的增大,系统中会出现随机内部激变和随机边界激变,即由于随机吸引子与其吸引域内的随机鞍发生碰撞而产生的随机分岔现象和由于随机吸引子与其吸引域边界发生碰撞而产生的随机分岔现象.研究结果揭示了随机干扰对永磁同步风力发电机运行性能影响的机理,为永磁同步风力发电系统的运行和设计提供了理论依据.     

Inference of a nonlinear stochastic model of the cardiorespiratory interaction

We reconstruct a nonlinear stochastic model of the cardiorespiratory interaction in terms of a set of polynomial basis functions representing the nonlinear force governing system oscillations. The strength and direction of coupling and noise intensity are simultaneously inferred from a univariate blood pressure signal. Our new inference technique does not require extensive global optimization, and it is applicable to a wide range of complex dynamical systems subject to noise.     

An Extension of Algorithm on Symbolic Computations of Conserved Densities for High-Dimensional Nonlinear Systems

An improved algorithm for symbolic computations of polynomial-type conservation laws (PCLaws) of a general polynomial nonlinear system is presented. The algorithm is implemented in Maple and can be successfully used for high-dimensional models. Furthermore, the algorithm discards the restriction to evolution equations. The program can also be used to determine the preferences for a given parameterized nonlinear systems. The code is tested on several known nonlinear equations from the soliton theory.     

Sequential reconstruction of driving-forces from nonlinear nonstationary dynamics

This paper describes a functional analysis-based method for the estimation of driving-forces from nonlinear dynamic systems. The driving-forces account for the perturbation inputs induced by the external environment or the secular variations in the internal variables of the system. The proposed algorithm is applicable to the problems for which there is too little or no prior knowledge to build a rigorous mathematical model of the unknown dynamics. We derive the estimator conditioned on the differentiability of the unknown system’s mapping, and smoothness of the driving-force. The proposed algorithm is an adaptive sequential realization of the blind prediction error method, where the basic idea is to predict the observables, and retrieve the driving-force from the prediction error. Our realization of this idea is embodied by predicting the observables one-step into the future using a bank of echo state networks (ESN) in an online fashion, and then extracting the raw estimates from the prediction error and smoothing these estimates in two adaptive filtering stages. The adaptive nature of the algorithm enables to retrieve both slowly and rapidly varying driving-forces accurately, which are illustrated by simulations. Logistic and Moran-Ricker maps are studied in controlled experiments, exemplifying chaotic state and stochastic measurement models. The algorithm is also applied to the estimation of a driving-force from another nonlinear dynamic system that is stochastic in both state and measurement equations. The results are judged by the posterior Cramer-Rao lower bounds. The method is finally put into test on a real-world application; extracting sun’s magnetic flux from the sunspot time series.     

First passage time for a class of one-dimensional stochastic systems

We consider the time evolution of a class of stochastic systems of finite size with polynomial nearest neighbor transition rates. We obtain analytical expressions for the first passage time (FPT) and its moments. We show that the mean FPT, averaged over a uniform initial distribution, shows a simple asymptotoc behavior with the system size and the parameters of the transition rates.     

CUMULANTS OF STOCHASTIC RESPONSE FOR A CLASS OF SPECIAL NONHOLONOMIC SYSTEMS

This paper studies the response cumulants for a kind of special nonholonomic systems under non-Gaussian, delta-correlated excitations. We present a new methodology for formulating the equations governing the evolution of the response cumulants of the stochastic dynamic systems. The response cumulant differential equations (CDEs) derived can be used to calculate the response cumulants for both linear and nonlinear systems. One example is given to illustrate how to use the CDEs for calculating response cumulants.     

Generalized spectral decomposition for stochastic nonlinear problems

We present an extension of the generalized spectral decomposition method for the resolution of nonlinear stochastic problems. The method consists in the construction of a reduced basis approximation of the Galerkin solution and is independent of the stochastic discretization selected (polynomial chaos, stochastic multi-element or multi-wavelets). Two algorithms are proposed for the sequential construction of the successive generalized spectral modes. They involve decoupled resolutions of a series of deterministic and low-dimensional stochastic problems. Compared to the classical Galerkin method, the algorithms allow for significant computational savings and require minor adaptations of the deterministic codes. The methodology is detailed and tested on two model problems, the one-dimensional steady viscous Burgers equation and a two-dimensional nonlinear diffusion problem. These examples demonstrate the effectiveness of the proposed algorithms which exhibit convergence rates with the number of modes essentially dependent on the spectrum of the stochastic solution but independent of the dimension of the stochastic approximation space.