精品（五）代数方程的求解.ppt

* 适用于非对称 矩阵的Bi-conjugate gradient 算法如下: 2倍于CG的计算量，相同的收敛速度，鲁棒性好 * 其他解法 CGSTAB (稳定化的CG) √GMRES (Saad and Shultz, 1986) * 5.3.8 Multigrid methods 大多数迭代法在细网格上可以很快消除误差的高频分量，但低频分量相当顽固。可以在粗网格上消除这些低频分量。 * 典型V循环式多重网格法的粗网格、限制和插值算子 * 两级线性多重网格法步骤 多级多重网格法：继续向更粗的网格限制，直到无更粗的网格为止。在最粗网格上精确求解修正方程。 * 公式描述:线性方程 * 公式描述:非线性方程 * 限制和插值算子: 对于eq (5.63) 1/2*eq(i-1)+ ?*eq(i+1) +eq(i) results in: * Comparison of count for convergence On 2D Poisson equation, k*k grid, n=k2, unknown Method Cost Gaussian elimination O(n3) GS O(n2logn) CG O(n1.5) FFT/Cyclic reduction O(nlogn) Multigrid O(n) optimal * 选择solver MG+SIP or MG +GS ICCG SIP ADI GS GMRES+MG 没有MG时， ICCGSIP ADI GS * 5.4 其他迭代法coupled equations (system of nonlinear equations) Simultaneous approach: All equations are considered part of a single system. Sequential approach: Each equation is solved for its dominant variable, treating the other as known, and one iterates through the equations until the solution of the coupled system is obtained. Iterations performed on each equation are called inner iteration. In order to obtain a solution which satisfy all equations, the coefficient matrices and source vector must be updated after each cycle ad the process repeated. The cycle are called outer iteration. * Sequential solution: Under-Relaxation On the nth iteration the equation for generic variable is Patankar 1980对SIMPLE采用, 稳定求解，但可能降低收敛率 时间相关法就是一种松驰法。 * 5.4.2延迟修正办法deferred-correction approaches 对于高阶差分离散，如果左端项用高阶差分,则计算复杂 如果左端项只保留相邻点的项，远邻点移到右端，则计算可能发散 为克服上述困难，可用延迟修正: 高阶差分移到右端，同时在左右两端加仅涉及相邻点的低阶差分： 用途：可以处理隐式高阶差分、交叉导数项、非正交相等。但不能处理高阶导数项。 * 第5次课阅读提示 傅《计算流体力学》第5章全部。 Peric《书》Chapter 5 全部。 * 第五次课后作业 实践Peric《书》附的代数方程求解程序 （待具体 化） 　PPT（全名：PowerPoint）是美国微软公司出品的办公软件系列重要组件之一（另外还有Excel、Word等）；最常用的就是用来制作幻灯片。MicrosoftOfficePowerPoint是一种演示文稿图形程序，该软件是功能强大的演示文稿制作软件；可协助用户独自或联机创建永恒的视觉效果；它增强了多媒体支持功能，利用该软件制