从噪声数据学习偏微分方程的复合神经网络

提出一种名为NN-PDE(neural network-partial differential equations)的复合神经网络方法, 用于噪声数据预处理和学习偏微分方程。NN-PDE用一套神经网络负责数据预处理, 另一套网络耦合备选的方程信息, 进而学习潜在的控制方程。两套网络复合为一套网络, 可更加高效地处理噪声数据, 有效减小噪声的影响。使用NN-PDE学习多种物理方程(如Burgers方程、Korteweg-de Vries方程、Kuramoto-Sivashinsky方程和Navier-Stokes方程)的噪声数据, 均可获得准确的控制方程。     

(2+1) 维Broer-Kau-Kupershmidt方程一系列新的精确解

借助于符号计算软件Maple，通过一种构造非线性偏微分方程(组)更一般形式精确解的直接方法即改进的代数方法，求解(2+1) 维 Broer-Kau-Kupershmidt方程，得到该方程的一系列新的精确解，包括多项式解、指数解、有理解、三角函数解、双曲函数解、Jacobi 和 Weierstrass 椭圆函数双周期解. 关键词： 代数方法 (2+1) 维 Broer-Kau-Kupershmidt 方程 精确解 行波解     

完整力学系统的高阶运动微分方程

从质点系的牛顿动力学方程出发,引入系统的高阶速度能量,导出完整力学系统的高阶Lagrange方程、高阶Nielsen方程以及高阶Appell方程,并证明了完整系统三种形式的高阶运动微分方程是等价的.结果表明,完整系统高阶运动微分方程揭示了系统运动状态的改变与力的各阶变化率之间的联系,这是牛顿动力学方程以及传统分析力学方程不能直接反映的.因此,完整系统高阶运动微分方程是对牛顿动力学方程及传统Lagrange方程、Nielsen方程、Appell方程等二阶运动微分方程的进一步补充. 关键词： 高阶速度能量 高阶Lagrange方程 高阶 Nielsen方程 高阶Appell方程     

具有三个任意函数的变系数KdV-MKdV方程的精确类孤子解

利用一个新的变换将变系数KdV-MKdV方程约化为三阶非线性常微分方程(NODE),考虑这个NODE,获得了变系数KdV-MKdV方程的若干精确类孤子解.这种思路也适合于其他的变系数非线性方程,如变系数KP方程、变系数sine-Gordon方程等. 关键词：     

(2+1)维Konopelchenko-Dubrovsky方程新的多孤子解

利用齐次平衡方法,将(2+1)维Konopelchenko-Dubrovsky方程转化为两个变量分离的线性偏微分方程,然后采用三种不同的函数假设,得到相应的常系数微分方程,通过求解特征方程,方便地构造出Konopelchenko-Dubrovsky方程新的多孤子解.     

改进的tanh函数方法与广义变系数KdV和MKdV方程新的精确解

利用改进的tanh函数方法将广义变系数KdV方程和MKdV方程化为一阶变系数非线性常微分方 程组-通过求解这个变系数非线性常微分方程组，获得了广义变系数KdV方程和MKdV方程新的 精确类孤子解、有理形式函数解和三角函数解- 关键词： 改进的tanh函数方法 类孤子解 有理形式函数解 三角函数解     

变质量力学系统的三阶拉格朗日方程

本文从变质量力学系统的三阶D’Alember-Lagrange原理出发，导出了变质量力学系统的Lagrange方程。利用该方程可以使我们得到描述变质量力学系统的运动。此外，变质量力学系统的Lagrange方程也可以使三阶运动微分方程理论得到充实。     

求解含有高阶导数偏微分方程的局部间断Petrov-Galerkin方法

构造一类求解三种类型偏微分方程的间断Petrov-Galerkin方法.求解的方程分别含有二阶、三阶和四阶偏导数,包括Burgers型方程、KdV型方程和双调和型方程.首先将高阶微分方程转化成为与之等价的一阶微分方程组,再将求解双曲守恒律的间断Petrov-Galerkin方法用于求解微分方程组.该方法具有四阶精度且具有间断Petrov-Galerkin方法的优点.数值实验表明该方法可以达到最优收敛阶而且可以模拟复杂波形相互作用,如孤立子的传播及相互碰撞等.     

一个新的对比饱和蒸气压方程

一、前言 目前,大多数饱和蒸气压方程都是在克拉贝隆-克劳修斯方程的基础上,假定△H_v和△Z_v随温度呈一定的变化关系并加以修正再进行积分得出的,例如Riedel方程、Miller方程、Thek-Stiel方程、Gomez Nieto-Thodos方程和徐忠方程等。本文在总结前人工作的基础上,从统计力学的角度,对分子结构及分子间作用力进行适当简化,先用分析的方法导出饱和压力随温度变化关系的基本函数形式,再用最小二乘回归的方法拟合出一套系数,从而得到纯流体从三相点到临界点的饱和蒸气压方程。     

求解含有高阶导数偏微分方程的局部间断Petrov-Galerkin方法（英文）

构造一类求解三种类型偏微分方程的间断Petrov-Galerkin方法.求解的方程分别含有二阶、三阶和四阶偏导数,包括Burgers型方程、KdV型方程和双调和型方程.首先将高阶微分方程转化成为与之等价的一阶微分方程组,再将求解双曲守恒律的间断Petrov-Galerkin方法用于求解微分方程组.该方法具有四阶精度且具有间断Petrov-Galerkin方法的优点.数值实验表明该方法可以达到最优收敛阶而且可以模拟复杂波形相互作用,如孤立子的传播及相互碰撞等.     

A standard and direct method to obtain multiple soliton solutions of the nonlinear partial differential equation

In this paper, we improve some key steps in the homogeneous balance method (HBM), and propose a modified homogeneous balance method (MHBM) for constructing multiple soliton solutions of the nonlinear partial differential equation (PDE) in a unified way. The method is very direct and primary; furthermore, many steps of this method can be performed by computer. Some illustrative equations are investigated by this method and multiple soliton solutions are found.     

广义(3+1)维Zakharov-Kuznetsov方程的对称约化、精确解和守恒律

运用李群分析,得到了广义(3+1)维Zakharov-Kuznetsov(ZK)方程的对称及约化方程,结合齐次平衡原理,试探函数法和指数函数法得到了该方程的群不变解和新精确解,包括冲击波解、孤立波解等.进一步给出了广义(3+1)维ZK方程的伴随方程和守恒律.     

An efficient method for model refinement in diffuse optical tomography

Diffuse optical tomography (DOT) is a non-linear, ill-posed, boundary value and optimization problem which necessitates regularization. Also, Bayesian methods are suitable owing to measurements data are sparse and correlated. In such problems which are solved with iterative methods, for stabilization and better convergence, the solution space must be small. These constraints subject to extensive and overdetermined system of equations which model retrieving criteria specially total least squares (TLS) must to refine model error. Using TLS is limited to linear systems which is not achievable when applying traditional Bayesian methods. This paper presents an efficient method for model refinement using regularized total least squares (RTLS) for treating on linearized DOT problem, having maximum a posteriori (MAP) estimator and Tikhonov regulator. This is done with combination Bayesian and regularization tools as preconditioner matrices, applying them to equations and then using RTLS to the resulting linear equations. The preconditioning matrixes are guided by patient specific information as well as a priori knowledge gained from the training set. Simulation results illustrate that proposed method improves the image reconstruction performance and localize the abnormally well.     

Finite element method for conserved phase fields: Stress-mediated diffusional phase transformation

Phase-field models with conserved phase-field variables result in a 4th order evolution partial differential equation (PDE). When coupled with the usual 2nd order thermo-mechanics equations, such problems require special treatment. In the past, the finite element method (FEM) has been successfully applied to non-conserved phase fields, governed by a 2nd order PDE. For higher order equations, the convergence of the standard Galerkin FEM requires that the interpolation functions belong to a higher continuity class.     

New evolution equations for the joint response-excitation probability density function of stochastic solutions to first-order nonlinear PDEs

By using functional integral methods we determine new evolution equations satisfied by the joint response-excitation probability density function (PDF) associated with the stochastic solution to first-order nonlinear partial differential equations (PDEs). The theory is presented for both fully nonlinear and for quasilinear scalar PDEs subject to random boundary conditions, random initial conditions or random forcing terms. Particular applications are discussed for the classical linear and nonlinear advection equations and for the advection–reaction equation. By using a Fourier–Galerkin spectral method we obtain numerical solutions of the proposed response-excitation PDF equations. These numerical solutions are compared against those obtained by using more conventional statistical approaches such as probabilistic collocation and multi-element probabilistic collocation methods. It is found that the response-excitation approach yields accurate predictions of the statistical properties of the system. In addition, it allows to directly ascertain the tails of probabilistic distributions, thus facilitating the assessment of rare events and associated risks. The computational cost of the response-excitation method is order magnitudes smaller than the one of more conventional statistical approaches if the PDE is subject to high-dimensional random boundary or initial conditions. The question of high-dimensionality for evolution equations involving multidimensional joint response-excitation PDFs is also addressed.     

A search for the simplest chaotic partial differential equation

A numerical search for the simplest chaotic partial differential equation (PDE) suggests that the Kuramoto-Sivashinsky equation is the simplest chaotic PDE with a quadratic or cubic nonlinearity and periodic boundary conditions. We define the simplicity of an equation, enumerate all autonomous equations with a single quadratic or cubic nonlinearity that are simpler than the Kuramoto-Sivashinsky equation, and then test those equations for chaos, but none appear to be chaotic. However, the search finds several chaotic, ill-posed PDEs; the simplest of these, in the discrete approximation of finitely many, coupled ordinary differential equations (ODEs), is a strikingly simple, chaotic, circulant ODE system.     

基于稀疏贝叶斯学习的复杂网络拓扑估计

提出了一种噪声环境下复杂网络拓扑估计方法, 仅利用含噪时间序列估计未知结构混沌系统的动力学方程和参数, 以及由混沌系统组成的复杂网络的拓扑结构、节点动力学方程、所有参数、 节点间耦合方向和耦合强度.通过采用动力学方程的统一形式, 将动力系统方程结构和参数估计看成线性回归问题的系数估计, 该估计问题利用贝叶斯压缩传感的信号重建算法求解, 含噪信号的模型重建使用相关向量机方法,即通过稀疏贝叶斯学习求解稀疏欠定线性方程得到上面提到的可估计对象.以单个Lorenz系统及由200个 Lorenz系统组成的无标度网络为例说明方法的有效性. 仿真结果表明,提出的方法对噪声有很强的鲁棒性,收敛速度快,稳态误差极小, 克服了最小二乘估计方法收敛速度慢、 稳态误差大以及压缩传感估计方法对噪声鲁棒性不强的缺点.     

A Novel Low-Dimensional Method for Analytically Solving Partial Differential Equations

This paper is concerned with a low-dimensional dynamical system model for analytically solving partial differential equations (PDEs). The model proposed is based on a posterior optimal truncated weighted residue (POT-WR) method, by which an infinite dimensional PDE is optimally truncated and analytically solved in required condition of accuracy. To end that, a POT-WR condition for PDE under consideration is used as a dynamically optimal control criterion with the solving process. A set of bases needs to be constructed without any reference database in order to establish a space to describe low-dimensional dynamical system that is required. The Lagrangian multiplier is introduced to release the constraints due to the Galerkin projection, and a penalty function is also employed to remove the orthogonal constraints. According to the extreme principle, a set of ordinary differential equations is thus obtained by taking the variational operation of the generalized optimal function. A conjugate gradient algorithm by FORTRAN code is developed to solve the ordinary differential equations. The two examples of one-dimensional heat transfer equation and nonlinear Burgers’ equation show that the analytical results on the method proposed are good agreement with the numerical simulations and analytical solutions in references, and the dominant characteristics of the dynamics are well captured in case of few bases used only.     

基于稀疏表示和特征加权的离格双耳声源定位*

基于头相关传递函数数据库的传统双耳声源定位方法的定位角度往往被限定在头相关传递函数数据库的离散测量点上。当头相关传递函数数据库的测量方位角间隔较大时,这类算法的性能会显著下降,这就是典型的离格问题。该文提出了基于加权宽带稀疏贝叶斯学习的离格双耳声源定位算法。首先该算法建立离格双耳信号的稀疏表示模型,然后利用双耳相干与扩散能量比特征对各个频点进行加权以降低噪声和混响的影响,最后通过加权宽带稀疏贝叶斯学习方法估计离格声源的方位角。实验结果表明,该算法在各种复杂的声学环境下都有着较高的定位精度和鲁棒性,特别是提高了离格条件下的声源定位性能。     

Double-Parameter Solutions of Projective Riccati Equations and Their Applications

The purpose of the present paper is twofold. First, the projective Riccati equations （PREs for short） are resolved by means of a linearized theorem, which was known in the literature. Based on the signs and values of coeffcients of PREs, the solutions with two arbitrary parameters of PREs can be expressed by the hyperbolic functions, the trigonometric functions, and the rational functions respectively, at the same time the relation between the components of each solution to PREs is also implemented. Second, more new travelling wave solutions for some nonlinear PDEs, such as the Burgers equation, the mKdV equation, the NLS^＋ equation, new Hamilton amplitude equation, and so on, are obtained by using Sub-ODE method, in which PREs are taken as the Sub-ODEs. The key idea of this method is that the travelling wave solutions of nonlinear PDE can be expressed by a polynomial in two variables, which are the components of each solution to PREs, provided that the homogeneous balance between the higher order derivatives and nonlinear terms in the equation is considered.