多重网格方法(Multigridmethod)

摘要:
多重网格法是一种求解微分方程的方法。Gauss–Seidel方法和Jacobi方法在前面介绍过,现在我们将使用这两种方法作为示例。尽管Gauss–Seidel方法转换速度更快,但对于高维系统而言,它实际上速度较慢。多网格方法的思想是首先粗略处理问题,将原始网格投影到相对简单的新网格上进行计算,然后在快速收敛后通过插值返回到原始系统。很容易发现,上述公式x(t+1)–x=Mx+B^(-1)B–x=Mx+Mx=M-与该公式等价。此外,我们可以使N=B^(-1):x(t+1)=Mx+Nb-不同的迭代方法实际上是A的不同分解方法,反映的公式是取不同的M和N值。

原文链接

多重网格方法是解微分方程的方法。这个方法的好处是在利用迭代法收敛结果的时候速度特别快。并且,不管是否对称,是否线性都无所谓。它的值要思想是在粗糙结果和精细结果之间插值。

前面介绍了Gauss–Seidel方法和Jacobi 方法,现在再用这两个方法来举例。尽管Gauss–Seidel (GS)方法converge更快一些,但其实对于维度很高的系统都很慢。Multigrid(MG)方法的思路是先把问题粗糙化,把原网格投影到一个比较简单的新网格上计算,等到快速收敛以后再经由Interpolation(插值)返回原来的系统。 

对于某个工程数学问题(如泊松方程),可以归纳为线形方程Ax = b,  A为n X n矩阵。那么,最终目的是得到所谓的x = A^(-1)b。定义e(t) = x – x(t),当e(t)为小于某个值的时候,可以认为xconverge到了合适的值。但实际上我们比较的是相邻的值。

把A非奇异分解A = B – C,

Bx – Cx = b

x = B(-1) Cx + B^(-1) b

并分开求解x

Bx(t+1) - Cx(t)= b

x(t+1) = B^(-1) Cx(t) + B^(-1) b                  -(1)

把形如B^(-1) C的矩阵称作迭代矩阵,用M表示。

容易发现,

x(t+1) – x = Mx(t) + B^(-1) b – x = Mx(t) +Mx = M (x(t) – x)   -(2)

以上式子与(1)式等价。

另外可令N = B^(-1):

x(t+1) = Mx(t) + Nb            - (3) 

不同的迭代方法其实就是A的不同分解法,反映到(3)式就是取不同的M和N值。

比如,在GS方法中B = D – L, C = U, 则(3)式为:

x(t+1) = (D - L)^(-1)Ux(t) + (D - L)^(-1)b            - (4)  

分析发现,当n的数值比较大时,以上收敛是极其缓慢的。假设p(i) 是序号为i的原系统有限元基函数(i = 1, 2, … , n),q(i)是粗糙化的网格(i = 1,2, … , m  m<n). 

一种粗糙化方法是构造矩阵H,使p = H * q, H为m X n矩阵。令A’ = HAH^T, x’ = Hx, b’ = Hb. 则 A’ x’ = b’ 是一个m维粗化的网格系统。 

更加具体来说,对于一个k维的问题,如果k小于一个指定的维数,那么直接用jacob等方法解救可以了。否则,粗糙化为更低的维度比如变为原来的二分之一,最后再把维数变回来。把mesh粗糙化的过程叫做Downv-cycle (从k维到k/2维), 反之叫Up v-cycle.(从k维到2k维)。

免责声明:文章转载自《多重网格方法(Multigridmethod)》仅用于学习参考。如对内容有疑问,请及时联系本站处理。

上篇jquery设置控件位置的方法centos7.6下定时监测MySQL进程终止后自动重启的方法下篇

宿迁高防,2C2G15M,22元/月;香港BGP,2C5G5M,25元/月 雨云优惠码:MjYwNzM=

相关文章

Unity里面的自动寻路(二)

        接着我的 上一篇自动寻路文章,这一次我们就来学习一下与自动寻路有关的组件吧。Unity中与自动寻路相关的组件主要有两个:NavMeshAgent (  又称导航网格代理 ),Off Mesh Link( 分离网格链接 )。这两个组件的作用与使用范围是不同的,我们唯一可以确定的是我们必须烘焙地形,产生NavMesh(导航网格)。因为导航网格决...

并发编程概述--C#并发编程经典实例

优秀软件的一个关键特征就是具有并发性。过去的几十年,我们可以进行并发编程,但是难度很大。以前,并发性软件的编写、调试和维护都很难,这导致很多开发人员为图省事放弃了并发编程。新版.NET 中的程序库和语言特征,已经让并发编程变得简单多了。随着Visual Studio 2012 的发布,微软明显降低了并发编程的门槛。以前只有专家才能做并发编程,而今天,每一个...

Direct3D11学习:(九)绘制基本几何体

转载请注明出处:http://www.cnblogs.com/Ray1024 一、概述 Direct3D中很多复杂的几何效果都是由基本的几何体组合而成的,这篇文章中,我们来学习集中常见的基本几何体的绘制方法。 二、准备工作 我们使用一个类来组织这些绘制基本几何体的代码,以方便我们以后的使用。GeometryGenerator是一个工具类,用于生成诸如网格、...

实时数据集成

  企业应用集成 面向服务的体系结构 (SOA) 目前应该是一个很受欢迎的名词,中间件技术人员几乎到了言必称SOA的程度,数据集成当然也不例外,在Oracle openworld2008大会上,就推出了一堆数据集成的专场演讲,其中和SOA结合最紧密的就是实时数据集成 real time data integration。我总结了一下,实时数据集成一般分为两...

cesium风场图【译】

cesium风场图 以下大部分内容来源于 开源博客 的英文翻译,为自己以后查阅方便 。 这是Rayman Ng的特邀帖子,介绍了他基于CesiumJS构建的开源风向图。 风是研究天气和气候的重要元素,它以各种方式影响我们的日常生活。在许多领域,例如气候分析和风电场管理,分析风是至关重要的。可视化它对于能够快速了解​​测量设备收集的数值风数据至关重要。已经有...

grid栅格布局

前面的话 Grid布局方式借鉴了平面装帧设计中的格线系统,将格线运用在屏幕上,而不再是单一的静态页面,可以称之为真正的栅格。本文将详细介绍grid布局 引入 对于Web开发者来说,网页布局一直是个比较重要的问题。但实际上,在网页开发很长的一段时间当中,我们甚至没有一个比较完整的布局模块。总的来说 Web 布局经历了以下四个阶段: 1、table表格布局,通...