解有界变量凸规划问题的ODE方法

用微分方程的解曲线确定约束优化问题的解即ODE方法已受到人们广泛重视和研究.潘平奇对无约束和带等式约束优化问题提出了很好的ODE方法.该方法的主要优点之一是没有扩大问题的规模.关于带不等式约束的优化问题的ODE方法,尚待研究.另外,虽然问题(1)可以通过标准化处理变成等式约束情形,再用[3]中的ODE方法求解,但这样做会扩大问题规模,因此,本文将在不扩大问题规模的基础上     

一种通用的复模态矩阵摄动法

对于非自伴随系统,提出了一种通用的复模态矩阵摄动法.该法能同时适用于孤立特征值,重特征值及密集(相近)特征值三种复特征值情况.由复特征子空间缩聚技术求解低阶摄动项,高阶摄动项则由逐次逼近过程求得.三个计算实例表明,该通用方法合理可靠,精度高.     

一类具有非线性发生率的戒烟模型动力学分析

对一类具有非线性发生率的PLSQ戒烟模型进行研究.首先,确定了模型的基本再生数,并分析了模型平衡点的存在性.然后,讨论了模型平衡点的局部稳定性和全局稳定性.最后,将发生率具体化,采用微分方程的非标准差分法(NSFD)求解戒烟模型,并将数值结果与ODE23和ODE45的计算结果比较.     

求解特征值问题的 AOR 算法的收敛性

用 AOR 方法求解线性方程组是众所周知的,我们将此方法应用到求解特征值问题方面.考虑下面特征值问题:(A—λI)x=0,(1.1)这里 A 是大型稀疏非奇异对称矩阵.显然,问题(1.1)有下面三条性质:i)其 n 个特征值都是实的,不妨设为λ_1≤λ_2≤…≤λ_n;(1.2)     

Schrdinger特征值问题的Legendre-Galerkin-Chebyshev配置法

提出了求解多维的Schrodinger特征值问题的Legendre-Galerkin-Chebyshev配置法(LGCC).LGCC方法是Legendre-Galerkin方法和Chebyshev配置法的耦合,便于处理变系数和非线性项且保持了Legendre-Galerkin方法的稳定性和高精度.对于拟线性SchrodingerPoisson特征值问题,建立了基于同伦连续法的LGCC方法.数值结果显示了该方法的有效性.     

Dodd-Bullough-Mikhailov方程的精确解

利用(G'/G)法求解了Dodd-Bullough-Mikhailov的精确解,得到了Dodd-Bullough-Mikhailov方程的用双曲函数,三角函数和有理函数表示的三类精确行波解.由于方法中的G为某个二阶常系数线性ODE的通解,故方法具有直接、简洁的优点;更重要的是,方法可用于求得其它许多非线性演化方程的行波解.如果对其中双曲函数表示的行波解中的参数取特殊值,那么可得已有的孤波解.     

关于矩阵An的特征值的研究

从一道线性代数习题出发,举例说明常见教材中关于由矩阵A的特征值确定ψ(A)的特征值的结论不够完备,进而分析问题关键,运用求解特征多项式的方法推导出矩阵多项式的特征值.     

求解特征值互补问题的一类ABS算法

ABS算法是20世纪80年代初,由Abaffy,Broyden和Spedicato完成的用于求解线性方程组的含有三个参量的投影算法,是一类有限次迭代直接法。目前,ABS算法不仅可以求解线性与非线性方程组,还可以求解线性规划和具有线性约束的非线性规划等问题。本文即是利用ABS算法求解特征值互补问题的一种尝试,构造了求解特征值互补问题的ABS算法,证明了求解特征值互补问题的ABS算法的收敛性。数值例子充分验证了求解特征值互补问题的ABS算法的有效性。     

一类Toeplitz线性代数方程组的预处理GMRES方法

本文研究一类来源于分数阶特征值问题的Toeplitz线性代数方程组的求解.构造Strang循环矩阵作为预处理矩阵来求解该Toeplitz线性代数方程组,分析了预处理后系数矩阵的特征值性质.提出求解该线性代数方程组的预处理广义极小残量法（PGMRES）,并给出该算法的计算量.数值算例表明了该方法的有效性.     

一个带固定步长的ODE型算法

提出了一种新的求解无约束优化问题的ODE型方法,其特点是:它在每次迭代时仅求解一个线性方程组系统来获得试探步;若该试探步不被接受,算法就沿着该试探步的方向求得下一个迭代点,其中步长通过固定公式计算得到.这样既避免了传统的ODE型算法中为获得可接受的试探步而重复求解线性方程组系统,又不必执行线搜索,从而减少了计算量.在适当的条件下,还证明了新算法的整体收敛性和局部超线性收敛性.数值试验结果表明:提出的算法是有效的.     

Construction of ODE Curves

This paper deals with a construction problem of free-form curves from data constituted by some approximation points and a boundary value problem for an ordinary differential equation (ODE). The solution of this problem is called an ODE curve. We discretize the problem in a space of B-spline functions. Finally, we analyze a graphical example in order to illustrate the validity and effectiveness of our method.     

A preprocessing method for parameter estimation in ordinary differential equations

Parameter estimation for nonlinear differential equations is notoriously difficult because of poor or even no convergence of the nonlinear fit algorithm due to the lack of appropriate initial parameter values. This paper presents a method to gather such initial values by a simple estimation procedure. The method first determines the tangent slope and coordinates for a given solution of the ordinary differential equation (ODE) at randomly selected points in time. With these values the ODE is transformed into a system of equations, which is linear for linear appearance of the parameters in the ODE. For numerically generated data of the Lorenz attractor good estimates are obtained even at large noise levels. The method can be generalized to nonlinear parameter dependency. This case is illustrated using numerical data for a biological example. The typical problems of the method as well as their possible mitigation are discussed. Since a rigorous failure criterion of the method is missing, its results must be checked with a nonlinear fit algorithm. Therefore the method may serve as a preprocessing algorithm for nonlinear parameter fit algorithms. It can improve the convergence of the fit by providing initial parameter estimates close to optimal ones.     

非线性偏微分方程的约化和精确解

§ 1　IntroductionSeeking the exact solutions of the nonlinear partial differential equation is one of thevery importantsubjectin PDE research.Up to now,many methods offinding the exact so-lutions for NLPDE are constructed,such as inverse scattering transformation(IST) [1 ] ,Liepoint symmetry and similar reductions[2 ,3] ,B cklund[4— 6] and Cole-Hofe transformations,Hirota s bilinear method[7] ,the homogeneous balance method[8,9] ,tanh function method[1 0 ]and so on.In this paper,we giv…     

Une méthode de décomposition en temps avec des schémas dʼintégration réversible pour la résolution de système dʼéquations différentielles ordinaires

We propose a time domain decomposition method that breaks the sequentiality of the integration scheme for systems of ODE. Under the condition of differentiability of the flow, we transform the initial value problem into a well-posed boundary values problem using the symmetrization of the interval of time integration and time-reversible integration scheme. For systems of linear ODE, we explicitly construct the block tridiagonal system satisfied by the solutions at the time sub-intervals extremities. We then propose an iterative algorithm of Schwarz type for updating the interfaces conditions which can extend the method to systems of nonlinear ODE.     

An exact solution of a quasilinear Fisher equation in cylindrical coordinates

Lie symmetry method is applied to analyse Fisher equation in cylindrical coordinates. Symmetry algebra is found and symmetry invariance is used to reduce the equation to a first-order ODE. The first-order ODE is further analysed to obtain exact solution of Fisher equation in explicit form.     

Continuous methods for extreme and interior eigenvalue problems

In this paper, continuous methods are introduced to compute both the extreme and interior eigenvalues and their corresponding eigenvectors for real symmetric matrices. The main idea is to convert the extreme and interior eigenvalue problems into some optimization problems. Then a continuous method which includes both a merit function and an ordinary differential equation (ODE) is introduced for each resulting optimization problem. The convergence of each ODE solution is proved for any starting point. The limit of each ODE solution for any starting point is fully studied. Both the extreme and the interior eigenvalues and their corresponding eigenvectors can be easily obtained under a very mild condition. Promising numerical results are also presented.     

Nonlinear dynamic analysis of a V-shaped microcantilever of an atomic force microscope

This paper is devoted to investigate the nonlinear behaviors of a V-shaped microcantilever of an atomic force microscope (AFM) operating in its two major modes: amplitude modulation and frequency modulation. The nonlinear behavior of the AFM is due to the nonlinear nature of the AFM tip–sample interaction caused by the Van der Waals attraction/repulsion force. Considering the V-shaped microcantilever as a flexible continuous system, the resonant frequencies, mode shapes, governing nonlinear partial and ordinary differential equations (PDE and ODE) of motion, boundary conditions, frequency and time responses, potential function and phase-plane of the system are obtained analytically. The governing PDE is determined by employing the Hamilton principle. Subsequently, the Galerkin method is utilized to gain the governing nonlinear ODE. Afterward, the resulting ODE is analytically solved by means of some perturbation techniques including the method of multiple scales and the Lindsted–Poincare method. In addition, the effects of different parameters including geometrical one on the frequency response of the system are assessed.     

An exact solution of a quasilinear Fisher equation in cylindrical coordinates

Lie symmetry method is applied to analyse Fisher equation in cylindrical coordinates. Symmetry algebra is found and symmetry invariance is used to reduce the equation to a first-order ODE. The first-order ODE is further analysed to obtain exact solution of Fisher equation in explicit form.     

Periodic wave solutions of coupled integrable dispersionless equations by residue harmonic balance

We introduce the residue harmonic balance method to generate periodic solutions for nonlinear evolution equations. A PDE is firstly transformed into an associated ODE by a wave transformation. The higher-order approximations to the angular frequency and periodic solution of the ODE are obtained analytically. To improve the accuracy of approximate solutions, the unbalanced residues appearing in harmonic balance procedure are iteratively considered by introducing an order parameter to keep track of the various orders of approximations and by solving linear equations. Finally, the periodic solutions of PDEs result. The proposed method has the advantage that the periodic solutions are represented by Fourier functions rather than the sophisticated implicit functions as appearing in most methods.     

Steady-state and Hopf bifurcations in the Langford ODE and PDE systems

This paper is concerned with the Langford ODE and PDE systems. For the Langford ODE system, the existence of steady-state solutions is firstly obtained by Lyapunov–Schmidt method, and the stability and bifurcation direction of periodic solutions are established. Then for the Langford PDE system, the steady-state bifurcations from simple and double eigenvalues are intensively studied. The techniques of space decomposition and implicit function theorem are adopted to deal with the case of double eigenvalue. Finally, by the center manifold theory and the normal form method, the direction of Hopf bifurcation and the stability of spatially homogeneous and inhomogeneous periodic solutions for the PDE system are investigated.