一个高效求解含时Gross-Pitaevskii方程的算法

Gross-Pitaevskii方程广泛应用于玻色-爱因斯坦凝聚体(Bose-Einstein condensate,BEC)的动力学研究,然而这个方程通常很难解析求解.因此发展相应的高精度数值求解方法非常重要.发展了结合算符劈裂法、Crank-Nicolson算法和四阶精度Numerov算法的高效求解Gross-Pitaevskii方程的新数值计算方法.通过数值计算可以表明,与传统的四阶精度的五点差分法相比,所提出的算法具有高效和消耗内存小的优点.     

求解BBM方程的一个新的高精度线性化差分算法

本文对一类带有齐次边界条件的Benjamin-Bona-Mahony方程的初边值问题进行了数值研究。通过先在时间层外推对问题进行线性化离散处理,然后再利用Richardson外推的思想在空间层进行外推,本文提出了一个理论精度为O(τ~2+h~4)的三层线性差分格式,证明了差分解的存在唯一性。在不能得到问题差分解的最大模估计的情况下,本文综合运用数学归纳法和离散泛函分析方法直接证明了该差分格式的收敛性和稳定性。数值实验表明该格式的精度明显优于已有的线性层差分格式。     

多重网格法求解原始变量形式的Navier-Stokes方程

采用多重网格法 ,求解原始变量形式的Navier Stokes方程 .控制方程在交错网格上离散 ,并采用SIMPLE算法计算 .为加速收敛 ,结合使用适合非线性方程的全近似格式的多重网格法FAS .文章计算了雷诺数Re =10 0 ,10 0 0时不可压粘性流体绕圆柱的流动 .计算结果与实验结果和国外的计算结果基本吻合 ,说明了该方法的准确性和有效性     

一个丢番图方程的求解

设P为素数,P D>1,证明了丢番图方程P2z-PzDm D2m=X2除D=2仅有非负整数解32-3·23 26=72和2 D,P≡1(mod8)时,必有m=2,z=1或m=1以及2|D,D含2hq 1形素因子(这里q|m,q为奇素数,h≥1)之外,推出m=1.     

混沌优化算法求解Kepler方程

利用混沌运动的遍历性、随机性及规律性等特点,对Kepler方程设计了一种混沌优化算法,对其进行求解,并在计算机上予以实现.该算法易于实现且时间复杂度较小,在计算效率上有一定优势.经过一系列数值实验测试,验证了混沌优化算法对Kepler方程求解的可行性,并将计算结果与已有解法所得结果相比较,表明了该算法对此方程的有效性.     

求解非线性发展方程的一个新的辅助椭圆方程

利用一个新的辅助椭圆方程将求解非线性发展方程精确解的问题转化为一个代数方程组进行求解,与已有的辅助椭圆方程法的主要不同是,应用这一新的辅助椭圆方程后降低了平衡次数,减少了所得的代数方程组的个数和方程的项数,从而大大地简化了代数方程组的求解．同时,由于辅助椭圆方程的解中包含了更多的可选参数,从而给出了非线性发展方程的更多形式的解．作为应用,借助于计算机的符号计算,求得了一些非线性发展方程的新的精确周期解．     

求解绝对值方程的光滑牛顿算法

背包问题以及大部分的线性互补问题都可以转化成为绝对值方程组来求解,求解绝对值方程Ax+B|a|=b是较难的问题.将该问题等价为线性互补问题,利用光滑牛顿法算求解该互补问题.当满足一定的条件时,证明了该算法是适定的,更证明了该算法的全局收敛性.利用Matlab软件对200维,500维,800维,和1 000维的情况进行了数值试验.每种情况测试了随机产生的50个可解的例子.精度达到了10-6.800维的用时在10 s左右,1000维的用时在20 s左右.     

求解Sylvester方程的正交迭代算法

对于任意初始矩阵,运用求解Sylvester矩阵方程的正交迭代算法可以在有限步内得到方程的最小二乘解,而且通过选择初始矩阵还可以得到方程的极小范数最小二乘解,这种算法还能用于解决最佳逼近问题,数值例子表明了所提出算法的有效性.     

求解加细方程的快速算法

多尺度加细方程的求解是小波分析与多分辨分析的应用基础。本文应用多分辨分析及矩阵的相关理论,采用迭代方法给出了求解m-尺度加细方程的快速算法(这里m>1;m∈Z)。该算法与己有的方法相比具有存储量小和运算速度快并且求解的范围更广泛的特点。     

一个快速求解哈斯图的算法

在理论证明的基础上，给出了一个哈斯图求解的算法。该算法可直接应用关系图方便、快捷地求取其盖住关系，进而得到哈斯图，从而完善了有关哈斯图求解的理论及算法。     

一个复方程混合边值问题的求解论证

讨论了一个复方程求解混合边值问题的方法，给出了求解问题的边界条件，应用基本定理论证了拟线性抛物型方程混合边值问题的唯一解及其存在性。     

求解mKdV方程呼吸子的一个新方法

给出了一种求mKdV方程呼吸子解的新方法．通过非线性变换把mKdV方程化成一种新的形式，假设其解具有分离变量的形式，在特定的条件下得到了该方程的一个特解——呼吸子解（也称为双子解）.     

方程求根的一个三阶算法

在经典的三阶迭代方法基础上,利用一阶导函数的差商代替函数本身的二阶导数,由此得到求解非线性方程根的一个不带二阶导数的迭代公式,并证明了公式具有三阶收敛性。     

间断Galerkin算法求解三维RANS方程

研究了间断Galerkin(DG)算法在三维可压缩湍流流场模拟中的应用。采用DG算法求解RANS方程和SA湍流模型方程,使用Roe格式计算无粘通量。针对RANS方程中的粘性项,发展了一种适用于DG算法的高效的高阶精度梯度计算方法。该方法以计算坐标系下守恒变量的梯度作为混合有限元法的辅助未知量,可使梯度计算达到与守恒变量相同阶的精度,并且不需要计算数值积分。采用该文描述的算法计算了NACA 64叶栅的三维湍流流场,得到了与实验数据吻合良好的结果,表明DG算法在叶轮机械内部湍流流场计算方面具有广阔的应用前景。     

基于QPSO算法的模糊关系方程求解

由于模糊关系方程在很多领域有着广泛的应用,越来越多的研究人员加入到求其极小解的工作中去,但结果不是很理想.该文用量子行为粒子群算法(QPSO)来求解其极小解,希望在模糊关系方程解的寻优问题上有所进展.将所提出的QPSO算法求解模糊关系方程的方法与遗传算法(GA)的求解方法作比较,发现QPSO优化算法较GA算法更能有效地找出模糊关系方程的近似最优解.     

求解有限元方程的区划分裂迭代算法

研究椭圆边值问题有限元离散方程的形成与求解，在区域分裂的基础上，构造出一类具有并行计算结构的迭代算法，并分析了算法的收敛速度。     

求解运动方程的一种数值算法

针对力学运动方程的数值计算，提出了一种类似欧拉法的新算法，其主要特点是没有算法衰减，与显式及隐式欧拉法相比具有较高的计算精度，且与显式欧拉法有相同的计算工作量，在对其算法稳定性进行讨论的基础上，以对单自由度质量弹簧振动模型的仿真计算为例，通过比较该算法与显式，隐式欧拉法及四阶龙格库塔法，阐明了该算法的这些特点。     

求解Rossler方程的格子Runge-Kutta-Boltzmann算法

提出一种求解Rossler方程的格子Runge-Kutta-Boltzmann算法. 先使用经典Runge-Kutta公式构建格子Boltzmann模型, 得到了高阶截断误差; 再通过Chapmann-Enskog展开和多尺度展开技术, 获得了不同时间尺度的系列偏微分方程和修正Rssler方程. 数值结果表明, 所得模型可用于求解Rossler方程.     

求解有限元方程的区域分裂迭代算法

研究椭圆边值问题有限元离散方程的形成与求解．在区域分裂的基础上，构造出一类具有并行计算结构的迭代算法，并分析了算法的收敛速度     

求解Pe11方程的算法与程序设计

本文给出n~(1/2)的次分数的递推算法与其渐近分数的高精度算法,完成求解Pell方程的Turbo　 C程序设计.