趋势分析中坏条件线性方程组的数值解

一、前言解线性方程组 A~Tx=b (1) 是地质数学中经常遇到的问题,其中 A∈R~(n×m),b∈R~m,x∈R~n,m≥n。在趋势分析中最后归结为解超定方程组A~Tx=b的问题。这里直接影响其效果的主要有两方面:一是趋势函数的选择;二是超定方程组解的精确度。在这里,我们仅讨论第二个问题。因为趋势分析是按最小二乘原则进行曲线拟合,而曲线拟合特别是多项式拟合导出的线性方程组,往往条件很坏。对于这样的问题,若用直接法求解,则精确度很差,有时甚至完全失真;若用迭代法求解,收敛速度很慢,有时甚至根本不收敛。故坏条件问题给数值计算带来了很大的困难。因此,讨论坏条件线性方程组的有效     

基于区域分解法的地下水有限元与边界元耦合模型—淄博市王旺庄水源地地下水数值模拟

区域分解法（DDM）是20世纪90年代兴起的一种求解偏微分方程的新方法，方法本身独到的耦合思想和高效的并行计算机理，对于求解复杂的、大型的地下水问题具有相当的优势和广阔的应用前景。本文以淄博市王旺庄水源地地下水流模型为例，应用重叠型区域分解法（DDM）构造了边界单元法（BEM）与有限单元法（FEM）耦合模型，在两种数值方法各自优点的基础上，更形象地再现了实际水文地质原型，有效地消除了人为边界造成的流场失真。