计算天文讲义_19887.pptVIP

12. 代数方程的数值求解 12 . 1 追赶法 12 . 2 迭代法 * * 双曲、抛物型1D隐式 在定常或非定常问题中，显格式可直接由迭代初值或当前时刻值算出新时刻的值，但隐格式则将原方程（组）化为代数方程组 双曲、抛物型2D隐式或椭圆型 需要对方程组进行联合求解。该方程组所含方程的个数等于因变量数与格点数的积，设为N。若采用通常的系数矩阵求逆或高斯消去法，每次求解将涉及N3/3次运算。 寻求有效算法 对一维隐式或多维问题(单个因变量)的交替方向隐式，每个点的差分方程只涉及三个未知量 A为n*n系数矩阵，它具有三对角形状 线性外推 等值外推 固定边界 它可能出现三项，但总可化为两项 例： 非开拓网格的对称边界 作业 作业 且b2=-c2 增广矩阵 其递推关系 令第i-1行已求得，则 消去 此过程由上依次向下的递推过程俗称“追”。 增广矩阵 其解的递推关系为 此过程由下依次向上的回代过程俗称赶”。 追步计算量：2n次除， 2n次乘， 2n次加 赶步计算量：n次乘， n次加 共2n次除， 3n次乘， 3n次加。 内存占用量7N。但若让 分别和 共享内存，则占用量为4N。 计算量是显格式的2倍。而高斯消去法的计算量约 次乘除，内存单元 对涉及m个因变量的偏微分方程组，则