求解扩散方程的一种高精度隐式差分方法

利用一阶微商和二阶微商的四阶紧致差分逼近公式,推导出了数值求解一维扩散方程的两种新的高精度隐式紧致差分格式,其截断误差分别为O(τ^2 h^4)和O(τ^4 h^4)．通过Fourier分析方法证明了格式O(τ^2 h^4)是无条件稳定的,而格式O(r^4 h^4)是无条件不稳定的．并且由于每一时间层上只用到了3个网格点,所以差分方程可采用追赶法直接进行求解．     

二维非定常对流扩散方程的隐式多重网格方法

提出了数值求解二维非定常对流扩散方程的一种新的加权平均隐格式，利用Fourier分析方法证明了该格式是无条件稳定的．利用多重网格加速技术。提出了基于时间修正的多重网格的全近似格式(FAS)，从而克服了传统迭代法在求解隐格式时收敛速度慢的缺陷，大大加快了迭代收敛速度。提商了问题的求解效率．数值计算结果表明，修正的多重网格FAS格式较传统的多重网格粗网格校正格式(CS)具有更好的收敛效率．并且随着σ和s的增大，它可以将传统迭代法的收敛速度提商几十倍，甚至几百倍．     

基于二维泊松方程六阶紧致格式的多重网格方法

利用六阶紧致差分格式、结合多重网格V循环算法求解了二维泊松方程的Dirichlet边值问题，并用不同的松驰算子与四阶精度格式的多重网格方法进行了比较，计算结果表明，该方法在不明显增加计算量的前提下较四阶精度格式的多重网格方法具有更好的精确度和收敛阶，且ZLGS迭代不论对四阶精度还是对六阶精度格式的多重网格算法，都是一种较其他松弛算子更加有效的“光滑剂”。     

求解一维线性双曲型方程的高精度紧致差分格式

【目的】双曲型方程是一类重要的偏微分方程,由于寻求问题本身的精确解比较困难,数值方法来求解此类方程有极具深远的意义和实际应用价值。【方法】首先对于一维的线性双曲型方程,在空间上采用Kreiss提出的四阶紧致差分公式进行逼近,时间上采用Taylor级数展开及截断误差修正的方法,推导出一个隐式的紧致差分格式。【结果】该格式在时间和空间上都有四阶精度,截断误差为O(τ4+h4)。【结论】采用Fourier方法分析了该格式的稳定性。数值实验证明提出的格式具有较好的稳定性和精确性。
     

三维热传导方程恒稳定的高精度半显式差分方法

提出了数值求解三维热传导方程的一种无条件稳定的高精度半显式差分方法,该方法可以显式计算且计算量小,截断误差为O(τ2+h4).数值算例验证了方法的精确性和可靠性.     

数值求解一维波动方程的四阶紧致差分方法

针对一维波动方程提出了一种有限差分方法.首先,采用泰勒级数展开公式和原方程代入的方法推导出了第一个时间层未知函数值的四阶紧致差分格式.然后,用四阶紧致差分公式近似空间导数项,采用中心差分格式截断误差余项修正的方法处理时间导数项,推导出了第二个时间层以后未知函数的四阶紧致差分格式.该方法时间和空间具有整体四阶精度.利用Fourier方法分析了所提格式的稳定性.由于本文格式在未知时间层仅涉及3个网格点,因此可采用追赶法求解离散化后所得到的线性方程组.最后,用数值算例验证了本文格式的精确性和稳定性.     

饰纹姬蛙Dmrt基因研究过程中异常结果的分析

对DMRT基因的研究需要运用到一系列的分子生物学实验技术,如DNA提取、简并PCR、DNA切胶纯化、克隆、SSCP分析及测序等.在实验的各个环节都可能会出现非预期的异常结果.对这些异常结果的分析研究也是实验的一个重要方面,不仅有助于促进实验的顺利进行,提高整个实验的效率,并且还有可能从中发现一些新的问题.     

三维波动方程的高精度交替方向隐式方法

提出了数值求解三维波动方程的2种精度分别为O(τ2 h4)和O(τ4 h4)的交替方向隐式(ADI)格式,并且通过Fourier方法证明了格式的稳定性.该方法在沿每个空间方向上只涉及三个网格基架点,因此可以重复采用TDMA算法,从而可以大大节省计算时间,数值实验验证了本文方法的精确性和可靠性.