求解对流扩散方程的一致四阶紧致格式

目前,许多高精度差分格式,由于未成功地构造与其精度匹配的稳定的边界格式,不得不采用低精度的边界格式.本文针对对流扩散方程证明了存在一致四阶紧致格式,它的边界点的计算格式和内点的计算格式的截断误差主项保持一致,给出了具体内点和边界格式;并分析了此半离散格式的渐近稳定性.数值结果表明该格式是四阶精度;在对流占优情况下,本文边界格式的数值结果比四阶精度的显式差分格式的的数值结果的数值振荡小,取得了不错的效果,理论结果得到了数值验证;驱动方腔数值结果显示,本文对N-S方程的离散格式具有很好的可靠性,适合对复杂流体流动的数值模拟和研究.     

线性传输方程的几种数值格式的比较

本文研究了线性传输方程的数值计算问题.利用Godunov格式、Entropy格式、Ultra-bee格式和Entropy-Ultra-bee格式对线性传输方程进行了数值计算,获得了相应的数值结果.数值实验结果表明Entropy-Ultra-bee格式结合了Entropy格式和Ultra-bee格式的优点,在整个计算区域都有比较高的分辨率,而且没有出现非物理振荡.     

对流扩散方程的数值计算

本文研究了对流扩散方程的一种并行格式.利用一组saul'yev型非对称格式进行二次构造,分别得到了一类并行GE格式和GEL、GER格式;进一步推广,得到绝对稳定的交替分组显式AGE格式,并用数值例子检验AGE格式的数值计算效果.     

具有三阶精度的数值微分紧致格式及其应用

本文从微分代数精度概念出发 ,引入了数值微分的紧致性概念 ,并且构造了在一般意义下的具有三点三阶精度的数值微分格式 ,以及在等距节点这种特殊情况下的计算格式 .最后 ,通过数值实验证实了格式的精度 .     

高分辨率差分格式的数值分析与限制因子的构造

本文利用差分方法余项效应理论,分析比较了一些典型的限制因子.对不同的限制因子,格式的表现明显差异主要是由其数值耗散性、色散性强弱不同所致.在分析比较格式的数值耗散性、色散性之后,本文提出了一种新的限制因子,得到的格式在解的剧烈变化区具有更高的分辨率,在光滑区避免了由于数值色散性较强导致的失真.数值试验表明该格式具有较好的性质.     

解Hamilton-Jacobi方程的MUSCL格式

本文利用单调数值通量和分片线性重构导数的方法构造了一种求HJ方程数值解的有限差分格式:MUSCL格式,并证明该格式具有TVB稳定性.数值实验表明该格式具有二阶精度,能避免产生伪振荡,尤其在类似"角点"的间断处有较好的分辩率.     

修正相场晶体(MPFC)模型的二阶无条件能量稳定数值格式

本文研究修正的相场晶体(MPFC)方程的二阶无条件能量稳定数值格式.首先,利用标量辅助变量(SAV)法和二阶向后欧拉(BDF2)公式,得到了一个数值格式;其次,给出了能量耗散定律,严格证明了该数值格式是质量守恒的,唯一可解的,无条件能量稳定的;最后,通过数值实例验证了格式的精度和稳定性.     

求解Fisher型方程的混合Jacobi-球面调和谱方法

本文提出了两种数值求解单位球内Fisher型方程的混合Jacobi-球面调和谱格式,并分别给出了格式的收敛性及相关的数值结果。     

非线性Schrdinger方程的高精度守恒差分格式

本文首先分析线性Schrodinger方程一种高阶差分格式的构造方法,得到方程的耗散项.在此基础上对三次非线性Schrodinger方程,提出了一种精度为O(r2 h2)的差分格式,证明了该格式保持了连续方程的两个守恒量,且是收敛的与稳定的.并通过数值例子与已有隐格式进行了比较,结果表明,本文格式在计算量类似的情况下,提高了数值精度.     

非结构网格上一类满足局部极值原理的三阶精度有限体积方法

对二维标量双曲型守恒律方程,发展了一类满足局部极值原理的非结构网格有限体积格式.其构造思想是,以单调数值通量为基础,通过应用基于最小二乘法的二次重构和极值限制器,使数值解满足局部极值原理.为保证数值解在光滑区域达到三阶精度,该格式可结合局部光滑探测器使用.本文从理论上分析了格式的稳定性条件,数值实验验证了格式的精度和对间断的分辨能力.     

三维数值流形方法的理论研究

在二维数值流形方法的基础上,对三维数值流形进行了理论研究．研究了三维覆盖位移函数,进行了三维数值流形的力学分析,给出了三维流形单元的刚度矩阵,详细推导了三维数值流形的Hammer积分及剖分规则,系统地研究三维数值流形的理论体系与数值实现方法．作为数值算例,给出了相应的悬臂梁的计算结果,计算结果表明算法的精度和计算效益较高．     

Efficient numerical methods for pricing American options under stochastic volatility

Five numerical methods for pricing American put options under Heston's stochastic volatility model are described and compared. The option prices are obtained as the solution of a two‐dimensional parabolic partial differential inequality. A finite difference discretization on nonuniform grids leading to linear complementarity problems with M‐matrices is proposed. The projected SOR, a projected multigrid method, an operator splitting method, a penalty method, and a componentwise splitting method are considered. The last one is a direct method while all other methods are iterative. The resulting systems of linear equations in the operator splitting method and in the penalty method are solved using a multigrid method. The projected multigrid method and the componentwise splitting method lead to a sequence of linear complementarity problems with one‐dimensional differential operators that are solved using the Brennan and Schwartz algorithm. The numerical experiments compare the accuracy and speed of the considered methods. The accuracies of all methods appear to be similar. Thus, the additional approximations made in the operator splitting method, in the penalty method, and in the componentwise splitting method do not increase the error essentially. The componentwise splitting method is the fastest one. All multigrid‐based methods have similar rapid grid independent convergence rates. They are about two or three times slower that the componentwise splitting method. On the coarsest grid the speed of the projected SOR is comparable with the multigrid methods while on finer grids it is several times slower. ©John Wiley & Sons, Inc. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

A note on estimating parameters in two-parameter Pareto distributions

Pareto distributions are used extensively in modelling income distributions. Estimation of parameters is revisited in two-parameter Pareto distributions. The method of quantile estimates using the elemental estimates and the method of product spacings are applied to the two-parameter Pareto distributions. A comparative study between the maximum likelihood method, the unbiased estimates which are functions of the maximum likelihood method, the minimum mean squared error method, the method of moments, the method of quantile estimation, the method of quantile estimation using the elemental estimates and the method of product spacings is presented.     

KdV方程的特征线混合间断有限元方法

考虑KdV方程的两种特征线性混合间断有限元方法,一种方法建立在标准特征线修正方法的基础上,另一种方法是带有对流项修正的特征线修正方法.利用具有实际物理意义的特征线追踪技巧处理时间导数项和对流项,采用混合间断有限元方法处理三阶导数项,分别证明了两种方法有限元解的唯一性、稳定性和误差估计.     

A partial inexact alternating direction method for structured variational inequalities

In this article, a new method is proposed for solving a class of structured variational inequalities (SVIs). The proposed method is referred to as the partial inexact proximal alternating direction (piPAD) method. In the method, two subproblems are solved independently. One is handled by an inexact proximal point method and the other is solved directly. This feature is the major difference between the proposed method and some existing alternating direction-like methods. The convergence of the piPAD method is proved. Two examples of the modern convex optimization problem arising from engineering and information sciences, which can be reformulated into the encountered SVIs, are presented to demonstrate the applicability of the piPAD method. Also, some preliminary numerical results are reported to validate the feasibility and efficiency of the piPAD method.     

微分连续正则化方法与一维声波方程系数反演问题求解

本文将求解非线性方程组的微分连续法与求解不问题的Ｔｉｋｈｏｎｏｖ正则化方法结合起来，形成了兼有大范围收性与稳定性的“微分连续正则化方法”，给出了方法的收敛性分析，文中最后给出的关于一维波动方程反问题的计算实例表明了方法的有效性。     

求解线性矩问题的一个修正方法

1．引言许多问题都可化为线性矩问题的求解,最典型的例子就是求解第一类线性积分方程的配置法．第一类算子方程在不同国数空间的离散化得到不同形式的线性矩问题．1968年,地质学家GBackusandFGilbert给出了一种求解线性矩问题的方法,用来求解地球物理反问题,后来称之为B－G方法,山从理论上严格论证了其收敛性．问给出了一种求解线形矩问题的光顺方法,门讨论了再生核空间的rG方法,并将其用于信号处理．以上方法的核心思想是在某类函数空间寻求对利一函数的逼近,从而得到线性矩问题的近似解．本交给出了一种求解线性矩问题的修正方…     

条件极值两种求法的不等价条件分析

本文说明了求条件极值的代入法与Ｌａｇｒａｎｇｅ乘数法的条件不等价。分析了这种不等价的原因，得到了求一般条件极值时，代入法与Ｌａｇｒａｎｇｅ乘数法均有效、前者无效而后者有效，以及两种方法均无效的各种不同条件．     

一种求多项式方程根的参数并行加速迭代法

推广了一种在无重根情况下,利用Newton类迭代法对同时求多项式零点的加速的迭代法.讨论了该方法的收敛性和收敛阶;最后给出数值算例表明:计算收敛阶和定理结论是一致的,且本算法具有较大的收敛范围.     

基于支持向量机的武器系统参数费用模型

在武器系统分析中,建立武器参数费用模型时,首先要挑选特征参数,这里采用R ough理论中的知识约简方法选择武器的特征参数;利用支持向量机建立了参数费用模型;给出了实例和解决此问题的支持向量机源程序.通过实例与线性回归法和神经网络法的结果进行了比较,结果表明支持向量机比较精确和简单.