限积分的计算方法

‘壹’ 有限积分法和有限差分法

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

1.2 差分格式
（1）从格式的精度来划分，有一阶格式、二阶格式和高阶格式。
（2）从差分的空间形式来考虑，可分为中心格式和逆风格式。
（3）考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。

目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。

1.3 构造差分的方法
构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。

2. FEM

2.1 概述
有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。

2.2 原理
有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学、土力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。

根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。
（1）从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法；
（2）从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格；
（3）从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。
不同的组合同样构成不同的有限元计算格式。

对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。

有限元插值函数分为两大类，一类只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值；另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状，一维看作长度比，二维看作面积比，三维看作体积比。在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元，常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。

2.3 基本原理与解题步骤
对于有限元方法，其基本思路和解题步骤可归纳为：
(1)建立积分方程，根据变分原理或方程余量与权函数正交化原理，建立与微分方程初边值问题等价的积分表达式，这是有限元法的出发点。
(2)区域单元剖分，根据求解区域的形状及实际问题的物理特点，将区域剖分为若干相互连接、不重叠的单元。区域单元划分是采用有限元方法的前期准备工作，这部分工作量比较大，除了给计算单元和节点进行编号和确定相互之间的关系之外，还要表示节点的位置坐标，同时还需要列出自然边界和本质边界的节点序号和相应的边界值。
(3)确定单元基函数，根据单元中节点数目及对近似解精度的要求，选择满足一定插值条件的插值函数作为单元基函数。有限元方法中的基函数是在单元中选取的，由于各单元具有规则的几何形状，在选取基函数时可遵循一定的法则。
(4)单元分析：将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近；再将近似函数代入积分方程，并对单元区域进行积分，可获得含有待定系数(即单元中各节点的参数值)的代数方程组，称为单元有限元方程。
(5)总体合成：在得出单元有限元方程之后，将区域中所有单元有限元方程按一定法则进行累加，形成总体有限元方程。
(6)边界条件的处理：一般边界条件有三种形式，分为本质边界条件(狄里克雷边界条件)、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件)。对于自然边界条件，一般在积分表达式中可自动得到满足。对于本质边界条件和混合边界条件，需按一定法则对总体有限元方程进行修正满足。
(7)解有限元方程：根据边界条件修正的总体有限元方程组，是含所有待定未知量的封闭方程组，采用适当的数值计算方法求解，可求得各节点的函数值。

3. 有限体积法

有限体积法（FiniteVolumeMethod）又称为控制体积法。其基本思路是：将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积；将待解的微分方程对每一个控制体积积分，便得出一组离散方程。其中的未知数是网格点上的因变量的数值。为了求出控制体积的积分，必须假定值在网格点之间的变化规律，即假设值的分段的分布的分布剖面。从积分区域的选取方法看来，有限体积法属于加权剩余法中的子区域法；从未知解的近似方法看来，有限体积法属于采用局部近似的离散方法。简言之，子区域法属于有限体积发的基本方法。有限体积法的基本思路易于理解，并能得出直接的物理解释。离散方程的物理意义，就是因变量在有限大小的控制体积中的守恒原理，如同微分方程表示因变量在无限小的控制体积中的守恒原理一样。限体积法得出的离散方程，要求因变量的积分守恒对任意一组控制体积都得到满足，对整个计算区域，自然也得到满足。这是有限体积法吸引人的优点。有一些离散方法，例如有限差分法，仅当网格极其细密时，离散方程才满足积分守恒；而有限体积法即使在粗网格情况下，也显示出准确的积分守恒。就离散方法而言，有限体积法可视作有限单元法和有限差分法的中间物。有限单元法必须假定值在网格点之间的变化规律（既插值函数），并将其作为近似解。有限差分法只考虑网格点上的数值而不考虑值在网格点之间如何变化。有限体积法只寻求的结点值，这与有限差分法相类似；但有限体积法在寻求控制体积的积分时，必须假定值在网格点之间的分布，这又与有限单元法相类似。在有限体积法中，插值函数只用于计算控制体积的积分，得出离散方程之后，便可忘掉插值函数；如果需要的话，可以对微分方程中不同的项采取不同的插值函数。

4. 比较分析
有限差分法（FDM）:直观，理论成熟，精度可眩但是不规则区域处理繁琐，虽然网格生成可以使FDM应用于不规则区域，但是对区域的连续性等要求较严。使用FDM的好处在于易于编程，易于并行。
有限元方法（FEM）:适合处理复杂区域，精度可眩缺憾在于内存和计算量巨大。并行不如FDM和FVM直观。不过FEM的并行是当前和将来应用的一个不错的方向。
有限容积法:适于流体计算，可以应用于不规则网格，适于并行。但是精度基本上只能是二阶了。FVM的优势正逐渐显现出来，FVM在应力应变，高频电磁场方面的特殊的优点正在被人重视。

比较一下：
有限容积法和有限差分法：一个区别就是有限容积法的截差是不定的（跟取的相邻点有关，积分方法离散方程），而有限差分就可以直接知道截差（微分方法离散方程）。有限容积法和有限差分法最本质的区别是，前者是根据积分方程推导出来的（即对每个控制体积分），后者直接根据微分方程推导出来，所以前者的精度不但取决于积分时的精度，还取决与对导数处理的精度，一般有限容积法总体的精度为二阶，因为积分的精度限制，当然有限容积法对于守恒型方程导出的离散方程可以保持守恒型；而后者直接由微分方程导出，不涉及积分过程，各种导数的微分借助Taylor展开，直接写出离散方程，当然不一定有守恒性，精度也和有限容积法不一样，一般有限差分法可以使精度更高一些。
当然二者有联系，有时导出的形式一样，但是概念上是不一样的。

至于有限容积法和有限元相比，有限元在复杂区域的适应性对有限容积是毫无优势可言的，至于有限容积的守恒性，物理概念明显的这些特点，有限元是没有的。目前有限容积在精度方面与有限元法有些差距。

有限元方法比有限差分优越的方面主要在能适应不规则区域，但是这只是指的是传统意义上的有限差分，现在发展的一些有限差分已经能适应不规则区域。对于椭圆型方程，如果区域规则，传统有限差分和有限元都能解，在求解效率，这里主要指编程负责度和收敛快慢、内存需要，肯定有限差分有优势。

‘贰’ 二重积分一共有多少种计算方法，分别是什么求归纳

二重积分一共一般有三种计算方法：变限求积分，直角坐标化极坐标，作图构思取最简单的微元。

当f(x,y)在区域D上可积时，其积分值与分割方法无关，可选用平行于坐标轴的两组直线来分割D，这时每个小区域的面积Δσ=Δx·Δy，因此在直角坐标系下，面积元素dσ=dxdy。可以看出二重积分的值是被积函数和积分区域共同确定的。

(2)限积分的计算方法扩展阅读：

对任意取定的x0∈[a，b]，过点（x0，0，0）作垂直于x轴的平面x=x0，该平面与曲顶柱体相交所得截面为区间，z=f（x0，y)为曲边的曲边梯形，由于x0的任意性。

其中y是积分变量在积分过程中视x为常数。上述曲顶柱体可看成平行截面面积S（x）从a到b求定积分的体积。

‘叁’ 积分上下限怎么算

定积分的计算就是把一个函数通过积分公式积分,再把定积分的上下限代入积分后的式子中,用代入上限的值减去代入下限的值.
当上限和下限的值一样时,代入上限的式子和代入下限的式子完全相同,相减就为0了.

‘肆’ 求问，这个变上限定积分是怎么算出来的

分部积分法，不过一般被积变量和上下限的变量会选择不同的表达，比如用t。

以分子为例，原本的定积分被积函数自变量是t，下限是0，上限是x。

令u=x-t，那么当t=0的时候，u=x；当t=x的时候，u=0。

所以当原本下限是t=0的时候，在新的定积分中，就是对应u=x。

在原本上限是t=x的时候，在新的定积分中，就是对应u=0。

所以这并不是什么上下限对调，而是根据u=x-t这个关系式，计算出当t=0和t=x的时候，u对应的值作为新的上下限。而这个“对调”，只是因为u=x-t这个关系的特殊性而已。

定积分

定积分的正式名称是黎曼积分。用黎曼自己的话来说，就是把直角坐标系上的函数的图象用平行于y轴的直线把其分割成无数个矩形，然后把某个区间［a，b］上的矩形累加起来，所得到的就是这个函数的图象在区间［a，b］的面积。实际上，定积分的上下限就是区间的两个端点a，b。

‘伍’ 变上限积分公式到底是怎样的

这里的意思就是
积分下限为a，下限是g(x)
那么对这个变上限积分函数求导，
就用g(x)代替f(t)中的t，
再乘以g(x)对x求导，即g'(x)
所以导数为f[g(x)] *g'(x)

‘陆’ 求极限积分的计算方法

方法如下图所示，
请认真查看，
祝学习愉快，
学业进步！

满意请釆纳！

‘柒’ 上限是无穷大的变限积分如何计算

上限无穷大的变限积分，先不管上下限，先把原函数写出来，然后此时的原函数当变量取无穷大的时候就相当于是取极限为一个定值。

积分下限为a，下限是g(x) 那么对这个变上限积分函数求导， 就用g(x)代替f(t)中的t， 再乘以g(x)对x求导。

即g'(x) 所以导数为f[g(x)]*g'(x)这里的意思就是积分下限为a，下限是g(x)，那么对这个变上限积分函数求导，就用g(x)代替f(t)中的t，再乘以g(x)对x求导，即g'(x)所以导数为f[g(x)] *g'(x)。

积分变限函数是一类重要的函数，它最着名的应用是在牛顿一莱布尼兹公式的证明中．事实上，积分变限函数是产生新函数的重要工具，尤其是它能表示非初等函数，同时能将积分学问题转化为微分学问题。积分变限函数除了能拓展我们对函数概念的理解外，在许多场合都有重要的应用。

(7)限积分的计算方法扩展阅读

反常积分总共就分两类：

1、积分上下限无界。

2、积分区域有界，函数在边界有暇点。

针对第二类，有如下的计算技巧。

∫baf(x)dx∫abf(x)dx，设在(a,b]上，在a处是暇点。

limx→a+f(x)(x−a)δ存在,δ∈(0,1)limx→a+f(x)(x−a)δ存在,δ∈(0,1) ，则积分收敛。

设在[a,b)上，b处是暇点。

limx→b−f(x)(x−b)δ存在,δ∈(0,1)limx→b−f(x)(x−b)δ存在,δ∈(0,1) ，则积分收敛。

‘捌’ 极限与积分的计算

问题出在第2个等号，因为这是n项求和，然后n趋于∞的极限

即无穷多项求和。是不可以交换求和与取极限的顺序的

这道题应该可以用定积分的定义以及夹逼定理来解，具体过程见图片

另：利用定积分定义求极限的问题一般已知的是和式极限，将其化为定积分有一定技巧，需逆向思维，在转化的过程中，关键要将和式中的每项化出一个1/n的因式

基本公式为

lim(n→∞)(1/n)∑f(k/n)=∫(0→1)f(x)dx(k=1→n)

‘玖’ 积分怎么算

计算定积分常用的方法：

拓展资料：

定积分的数学定义：如果函数f(x)在区间[a,b]上连续，用分点xi将区间[a,b]分为n个小区间，在每个小区间[xi-1,xi]上任取一点ri（i=1,2,3„,n），作和式f(r1)+...+f(rn)，当n趋于无穷大时，上述和式无限趋近于某个常数A，这个常数叫做y=f(x)在区间上的定积计做/abf(x)dx即/abf(x)dx＝limn>00[f(r1)+...+f(rn)]，这里，a与b叫做积分下限与积分上限，区间[a,b]叫做积分区间，函数f(x)叫做被积函数，x叫做积分变量，f(x)dx叫做被积式。

几何定义：可以理解为在Oxy坐标平面上，由曲线y=f(x)与直线x=a,x=b以及x轴围成的曲边梯形的面积值。（一种确定的实数值）

‘拾’ 变限积分求导公式是什么

F(x) = ∫(a,x) xf(t) dt

F(x) = x∫(a,x) f(t) dt

F'(x) = ∫(a,x) f(t) dt + x * [x' * f(x) - a' * f(a)]

= (1/x)F(x) + x * [1 * f(x) - 0 * f(a)]，下限a的导数是0，所以整体都会变为0

= (1/x)F(x) + xf(x)

求导注意事项：

（1）区间a可为-∞，b可为+∞；

（2）此定理是变限积分的最重要的性质，掌握此定理需要注意两点：第一，下限为常数，上限为参变量x（不是含x的其他表达式）；第二，被积函数f(x)中只含积分变量t，不含参变量x。

原函数存在定理

若函数f(x)在区间[a,b]上连续，则积分变上限函数就是f(x)在[a,b]上的一个原函数。