单参数非线性问题中高阶奇异点的计算

本文考虑计算单参数非线性问题中高阶奇异点的数值方法，基于确定奇异点的一个普适的扩张系统，结合同伦参数的拟弧长延拓，给出了计算各类高阶奇异点的一个统一算法，数值例子表明了算法的有效性．     

非线性动力系统中两鞍-结分岔点间非稳定曲线的确定

采用将伪弧长延拓法与Poincaré映射法相结合的方法,确定非自治动力系统中两鞍-结分岔点间非稳定曲线,并对采用一般延拓法时出现的奇异性进行了证明。该方法引入了一正则化方程,避免了在求解过程中出现的奇异问题,并给出了相应的迭代格式。在曲线的延拓过程中,由于存在两个延拓方向,为保证将曲线延拓出来,给出了一种确定切线方向的方法,该方法在分析非线性振动系统中的双稳态现象等问题是很有效的。     

用修正的平行割线法求解高阶奇异问题

许多实际问题中,方程在解点处的导算子为一高阶奇异算子,如反应扩散系统、优化问题中的歧点等.因此,对于求解高阶奇异问题的研究具有重要的实际意义.利用平行割线法求解高阶奇异问题,得到了渐近收敛速率,最后结合Hilbert空间几何特征,在几乎不增加计算量的前提下,修正了平行割线法,提高了渐近收敛速率.     

小波Galerkin法在非线性分岔问题求解中的应用

通过一个典型的Bratu问题,研究了小波Galerkin法（WGM）在非线性分岔问题求解方面的应用.首先,利用基于Coiflet的小波Galerkin法,对一维和二维Bratu方程进行离散;然后针对单参数问题,推导了追踪解曲线的伪弧长格式和直接计算极值型分岔点的扩展方程;针对双参数问题,推导了追踪稳定边界的伪弧长格式和求解尖点型分岔点的扩展方程.数值结果表明,基于小波Galerkin法的非线性分岔计算不仅具有更高的计算精度,而且能够有效地捕捉双参数分岔问题的折迭线和尖点突变曲面.该算例展示了基于小波Galerkin法的数值分岔计算的具体过程及其求解多参数分岔问题复杂行为的应用潜力.     

加热弹性杆的热过屈曲分析

基于轴线可伸长细杆的过屈曲变形几何理论,建立了两端轴向不可移的均匀加热直杆热弹性过屈曲行为的精确数学模型.这是一个包含杆轴线弧长在内的多未知函数的强非线性一阶常微分方程两点边值问题.采用打靶法和解析延拓法直接数值求解上述非线性边值问题,分别获得了两端横向简支和夹紧杆的热过屈曲状态解,给出了具有不同细长比杆的热过屈曲平衡路径.     

应用第二类Chebyshev小波求解高阶多点边值问题的数值解（英文）

本文研究了一类高阶多点边值问题的数值解法问题.利用第二类Chebyhsev小波及其积分算子矩阵,将线性与非线性高阶常微分方程多点边值问题转化为代数方程组进行求解.通过与现有文献算法结果的比较,说明了该算法求解高阶多点边值问题的准确性与有效性.扩展了高阶多点边值问题的数值求解方法.     

求解微分方程初值问题的一种弧长法

对于连续介质力学问题中导出的微分方程初值问题,常常具有解奇异性,如不连续、Ｓｔｉｆ性质或激波间断·本文通过在相应空间,引入一个或数个弧长参数变量,克服解的奇异性·对于常微分方程组引入弧长参数变量后,奇异性得以消除和削弱,应用一般的解常微分方程组的方法（如Ｒｕｎｇｅ＿Ｋｕｔａ法）求解·对于偏微分方程引入弧长参数变量后,在相应的空间离散成常微分方程组,用解奇异性常微分方程组相同的方法即可求解·本文给出了两个算例     

三维粘弹性分层介质中平稳随机波的传播

研究平稳随机波在粘弹性分层横观各向同性介质中的传播问题．将岩层考虑为分层介质,各层性质不同,岩层位于基岩上面,并且认为基岩比岩层刚很多,在基岩处给出随机激励．在频率和波数域中将控制方程化为常微分方程求解．对常微分方程,应用两点边值问题的精细积分法进行求解．因此,近年来发展的应用于结构随机振动的虚拟激励法可推广于当前分层岩层响应的计算．     

基于凝聚函数的拟牛顿算法求解绝对值方程

绝对值方程Ax-|x|=b是一个不可微的NP-hard问题.在假设矩阵A的奇异值大于1(这里矩阵A的奇异值定义为矩阵ATA特征值的非负平方根)时,给出了求解绝对值方程一个新的光滑化算法.通过引入一种凝聚函数对绝对值方程进行光滑化处理,得到一个非线性方程组;再引入适当的目标函数,进而把绝对值方程化为无约束优化问题,然后利用拟牛顿算法对其进行求解.数值实验结果表明了该方法的正确性和有效性.     

线性常微分方程的全局误差估计和优化求解方法

线性常微分方程初值问题求解在许多应用中起着重要作用.目前,已存在很多的数值方法和求解器用于计算离散网格点上的近似解,但很少有对全局误差(global error)进行估计和优化的方法.本文首先通过将离散数值解插值成为可微函数用来定义方程的残差;再给出残差与近似解的关系定理并推导出全局误差的上界;然后以最小化残差的二范数为目标将方程求解问题转化为优化求解问题;最后通过分析导出矩阵的结构,提出利用共轭梯度法对其进行求解.之后将该方法应用于滤波电路和汽车悬架系统等实际问题.实验分析表明,本文估计方法对线性常微分方程的初值问题的全局误差具有比较好的估计效果,优化求解方法能够在不增加网格点的情形下求解出线性常微分方程在插值解空间中的全局最优解.     

Study of parametric oscillation of an electrostatically actuated microbeam using variational iteration method

This paper deals with the study of parametric oscillation of an electrostatically actuated microbeam using variational iteration method. The paper considers a micro-beam suspended between two conductive micro-plates, subjected to a same actuation voltage. The nonlinear governing differential equation of motion about static equilibrium position using calculus of variation theory and Taylor series expansion has been linearized and implementing a Galerkin based reduced order model a Mathieu type equation has been obtained. By improving variational iteration method combining with method of strained parameters transition curves, separating stable from unstable regions have been obtained. The results of variational iteration method, perturbation and direct numerical integration methods for some cases selected from different regions (stable and unstable regions) have been compared.     

Numerical continuation of solution at a singular point of high codimension for systems of nonlinear algebraic or transcendental equations

Numerical continuation of solution through certain singular points of the curve of the set of solutions to a system of nonlinear algebraic or transcendental equations with a parameter is considered. Bifurcation points of codimension two and three are investigated. Algorithms and computer programs are developed that implement the procedure of discrete parametric continuation of the solution and find all branches at simple bifurcation points of codimension two and three. Corresponding theorems are proved, and each algorithm is rigorously justified. A novel algorithm for the estimation of errors of tangential vectors at simple bifurcation points of a finite codimension m is proposed. The operation of the computer programs is demonstrated by test examples, which allows one to estimate their efficiency and confirm the theoretical results.     

Numerical computation of stability and detection of Hopf bifurcations of steady state solutions of delay differential equations

The characteristic equation of a system of delay differential equations (DDEs) is a nonlinear equation with infinitely many zeros. The stability of a steady state solution of such a DDE system is determined by the number of zeros of this equation with positive real part. We present a numerical algorithm to compute the rightmost, i.e., stability determining, zeros of the characteristic equation. The algorithm is based on the application of subspace iteration on the time integration operator of the system or its variational equations. The computed zeros provide insight into the system’s behaviour, can be used for robust bifurcation detection and for efficient indirect calculation of bifurcation points. This revised version was published online in June 2006 with corrections to the Cover Date.     

An analog of the method of continuation of solution with respect to the parameter for nonlinear operator equations

A process of second order is constructed for the solution of nonlinear operator equations which is an analog of the method of continuation of solution with respect to the parameter. For each value of the parameter the Newton-Kantorovich iteration formula is applied only once in all. The quadratic convergence of the process is ensured by the specification of the parameter by a special formula. The process under consideration enables us to avoid the singular points of the derivative of the nonlinear operator on the left-hand side of the operator equation.Translated from Matematicheskie Zametki, Vol. 23, No. 4, pp. 601–606, April, 1978.     

Homotopy analysis method for solving the MHD flow over a non-linear stretching sheet

The problem of the boundary layer flow of an incompressible viscous fluid over a non-linear stretching sheet is considered. Homotopy analysis method (HAM) is applied in order to obtain analytical solution of the governing nonlinear differential equations. The obtained results are finally compared through the illustrative graphs with the exact solution and an approximate method. The compression shows that the HAM is very capable, easy-to-use and applicable technique for solving differential equations with strong nonlinearity. Moreover, choosing a suitable value of none–zero auxiliary parameter as well as considering enough iteration would even lead us to the exact solution so HAM can be widely used in engineering too.     

A block-Lanczos method for large continuation problems

The computation of solution paths of large-scale continuation problems can be quite challenging because a large amount of computations have to be carried out in an interactive computing environment. The computations involve the solution of a sequence of large nonlinear problems, the detection of turning points and bifurcation points, as well as branch switching at bifurcation points. These tasks can be accomplished by computing the solution of a sequence of large linear systems of equations and by determining a few eigenvalues close to the origin, and associated eigenvectors, of the matrices of these systems. We describe an iterative method that simultaneously solves a linear system of equations and computes a few eigenpairs associated with eigenvalues of small magnitude of the matrix. The computation of the eigenvectors has the effect of preconditioning the linear system, and numerical examples show that the simultaneous computation of the solution and eigenpairs can be faster than only computing the solution. Our iterative method is based on the block-Lanczos algorithm and is applicable to continuation problems with symmetric Jacobian matrices. This revised version was published online in June 2006 with corrections to the Cover Date.     

Numerical stability analysis and computation of Hopf bifurcation points for delay differential equations

We present a numerical technique for the stability analysis and the computation of branches of Hopf bifurcation points in nonlinear systems of delay differential equations with several constant delays. The stability analysis of a steady-state solution is done by a numerical implementation of the argument principle, which allows to compute the number of eigenvalues with positive real part of the characteristic matrix. The technique is also used to detect bifurcations of higher singularity (Hopf and fold bifurcations) during the continuation of a branch of Hopf points. This allows to trace new branches of Hopf points and fold points.     

A maximum principle for fourth order ordinary differential equations

We present a maximum principle for fourth order ordinary differential equations, based on a new approach involving counting of inflection points. We use our results to compute solutions of nonlinear equations describing static displacements of a uniform beam     

Continuation and shooting methods for boundary value problems of Bernstein type

Fixed point continuation methods and shooting methods are combined to produce an effective numerical procedure for solving boundary value problems for nonlinear ordinary differential equations. Typical numerical solution schemes involve an iteration procedure. Continuation methods systematically generate good initial guesses and, when combined with a shooting method and an appropriate update procedure, give systematic means for the numerical solution of nonlinear boundary value problems. This paper concentrates on problems of Bernstein type, which arise naturally in the calculus of variations and in steady-state heat conduction.