伽辽金法

原文链接

伽辽金方法(Galerkin method)是由俄罗斯数学家鲍里斯·格里戈里耶维奇·伽辽金(俄文:Борис Григорьевич Галёркин 英文:Boris Galerkin)发明的一种数值分析方法。应用这种方法可以将求解微分方程问题(通过方程所对应泛函变分原理)简化成为线性方程组的求解问题。而一个高维(多变量)的线性方程组又可以通过线性代数方法简化,从而达到求解微分方程的目的。

伽辽金法采用微分方程对应的弱形式,其原理为通过选取有限多项试函数(又称基函数形函数),将它们叠加,再要求结果在求解域内及边界上的加权积分(权函数为试函数本身)满足原方程,便可以得到一组易于求解的线性代数方程,且自然边界条件能够自动满足。

必须强调指出的是,作为加权余量法的一种试函数选取形式,伽辽金法所得到的只是在原求解域内的一个近似解(仅仅是加权平均满足原方程,并非在每个点上都满足)。

因为伽辽金方法的妙处在于研究它们的抽象方法,所以我们首先给出它们的抽象推导。最后我们再给出应用的例子。

一个问题的弱形式

我们通过一个抽象问题来引入伽辽金方法,将问题表示成在一个希尔伯特空间上的弱形式,也就是,求解使得对于所有

成立。这里,是一个双线性型表达式,是一个上的线性形表达式。

伽辽金离散化

选取一个n 维子空间,然后求解问题在子空间中的投影:求使得对于所有

我们称这个方程为伽辽金方程。注意方程形式没有改变,但是求解域改变了。

伽辽金正交性

这是使得伽辽金方法非常有效的关键性质。因为,我们可以取为原方程的一个试矢量。带入并相减,便得到误差的伽辽金正交性关系

这里是真实解和伽辽金方程的解之间的误差。

矩阵形式

因为伽辽金方法的目标是将问题简化为线性方程组,我们来构造它的矩阵形式,以便利用计算机进行数值求解。

空间中的一组。则显然依次选取这些基矢量作为伽辽金方程的试矢量是充分的,也即:求解使得

用上述基矢量表示出,将其代入上面的方程得到

这样我们就得到了上面这组型的线性方程组,式中

矩阵的对称性

由于矩阵项的定义,伽辽金方程的系数矩阵对称矩阵充要条件是双线性型表达式是对称的。

时间: 2024-10-14 00:59:19

伽辽金法的相关文章

稳态热传导的有限元分析

在分析工程问题时,经常要了解工件内部的温度分布情况,例如发动机的工作温度.金属工件在热处理过程中的温度变化.流体温度分布等.物体内部的温度分布取决于物体内部的热量交换,以及物体与外部介质之间的热量交换,一般认为是与时间相关的.物体内部的热交换采用以下的热传导方程(Fourier方程)来描述,         (6-1) 式中为密度,kg/m3: 为比热容,:为导热系数,:T为温度,℃:t为时间,s:为内热源密度,w/m3. 对于各向同性材料,不同方向上的导热系数相同,热传导方程可写为以下形式,

数值微分方程

F 有限元分析 H Heun方法 K Kansa方法 伽 伽辽金法 刚 刚性方程 努 努梅罗夫方法 奇 奇异边界法 快 快速行进算法 打 打靶法 欧 欧拉方法 正 正则化无网格法 瑞 瑞利-里兹法 蛙 蛙跳积分法 边 边界粒子法 边界节点法 韦 韦尔莱积分法 龙 龙格-库塔法

FDM, FVM, FEM

有限元法.有限差分法和有限体积法的区别 标签: 函数 有限元 插值 差分 格式 有限差分方法(Finite Differential Method)是计算机数值模拟最早采用的方法,至今仍被广泛运用.该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域.有限差分法以泰勒级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组.该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数

扩展欧几里得算法的模板实现

我居然现在还记不住扩欧的板子,我太弱啦! 扩展欧几里得算法解决的是这样的问题: 给定一个不定方程组ax+by=gcd(a,b),求他的一组整数解 先给出实现代码 void exgcd(int a,int b,int &x,int &y) { if(!b) { x=1,y=0;//gcd(a,0)显然等于1*a-0*0=a return a; } int ans=exgcd(b,a%b,x,y); int tem=x; x=y; y-=tem-(a/b)*y; return ans;} 但实

【精品】北京赛车计划冠军定位玩法技巧

車車是一種投資,我們的目標是:細水長流,見好就收,不求日金千金,只求長期穩定!許多人賠本的原因:1. 資金不足,卻大把下注,跟到第4期不出,錢不夠了,心慌了.有人孤注一擲,衝到第5期中了,嚇的半死.有人不敢跟,第5期出號了,氣死, 然後再跟新計劃,沒錢了,郁悶死.這兩種做法都不對,既然是以投資的心態做事,就應該計劃好翻倍的本錢,做到99%的穩賺,狀況不對就要及時止損. 看著連續中,就是不敢跟,最後咬牙跟了,馬上挂了.于是開始哭,我運氣咋這麽差.不買就中,一買就挂.相反,有些人專門等挂,一挂就上,

构建之法——读书笔记(8)

<构建之法>第十&十一章 主要讲述了在软件设计前期的需求分析问题上的方法和实践经验,分为"典型用户和场景"以及"软件设计与实现". 其中第十章大部分内容包含: 用户的分类(典型用户可以包括以下内容: 1. 名字(越自然越好) 2. 年龄(不同年龄和收入的用户有不同的需求) 3. 收入 4. 代表的用户在市场上的比例和重要性(比例大不等同于重要性高,如付费的用户比例较少,但是影响大,所以更重要 5. 使用这个软件的典型场景 6. 使用本软件/服务的

转:梯度下降法(上升法)的几何解释

梯度下降法是机器学习和神经网络学科中我们最早接触的算法之一.但是对于初学者,我们对于这个算法是如何迭代运行的从而达到目的有些迷惑.在这里给出我对这个算法的几何理解,有不对的地方请批评指正! 梯度下降法定义 (维基百科)梯度下降法,基于这样的观察:如果实值函数  在点  处可微且有定义,那么函数 在  点沿着梯度相反的方向  下降最快. 因而,如果                         对于  为一个够小数值时成立,那么 . 考虑到这一点,我们可以从函数  的局部极小值的初始估计  出发

1024. 科学计数法 (20)

科学计数法是科学家用来表示很大或很小的数字的一种方便的方法,其满足正则表达式[+-][1-9]"."[0-9]+E[+-][0-9]+,即数字的整数部分只有1位,小数部分至少有1位,该数字及其指数部分的正负号即使对正数也必定明确给出. 现以科学计数法的格式给出实数A,请编写程序按普通数字表示法输出A,并保证所有有效位都被保留. 输入格式: 每个输入包含1个测试用例,即一个以科学计数法表示的实数A.该数字的存储长度不超过9999字节,且其指数的绝对值不超过9999. 输出格式: 对每个测

五大常用算法之四:回溯法

(转自:http://www.cnblogs.com/steven_oyj/archive/2010/05/22/1741376.html) 1.概念 回溯算法实际上一个类似枚举的搜索尝试过程,主要是在搜索尝试过程中寻找问题的解,当发现已不满足求解条件时,就“回溯”返回,尝试别的路径. 回溯法是一种选优搜索法,按选优条件向前搜索,以达到目标.但当探索到某一步时,发现原先选择并不优或达不到目标,就退回一步重新选择,这种走不通就退回再走的技术为回溯法,而满足回溯条件的某个状态的点称为“回溯点”. 许