偏微分方程的有限元方法概述.
科学计一，但三者之间相辅相成。
偏微分方程的有限元方法
一 边值问题的变分原理
1 引论 （1）等周问题 在长度一定的所有平面封闭曲线中，求 所围面积为最大的曲线。
dx dy 模型：在条件 ds l 下 s1 ds ds 1 s2 dy dx 求使得泛函 s( x, y ) s x y ds 2 1 ds ds 达到最大的函数 x( s), y ( s) 。
s2 2 2
定义：当求泛函在一个函数集合K中的极小 （或极大）问题，则该问题称为变分问题。 变分问题与微分方程的定解问题有一定的 联系。 （2）初等变分原理 ① 一元二次函数的变分原理 设 J ( x) Lx 2 2 fx ( L 0, L, f为实常数） 考察J(x)的极值情况。 变分原理： 求 x0 R，使 J ( x0 ) min J ( x)
xR
与求解方程 Lx = f 等价。
② 多元二次函数的变分原理
设
n 1 n n J ( x ) J ( x1 , x2 ,, xn ) aij xi x j bi xi 2 i 1 j 1 i 1
求J(x)取极小值的驻点， 其中
a11 a12 a1n a a a 22 2n 对称正定 A 21 an1 an 2 ann T 设 x ( x1 , x2 ,, x n ) b (b1 , b2 ,,b n )T 1 则J(x)可表示为： J ( x ) ( Ax , x ) (b , x ) 2
（2）两点边值问题的变分原理 考察二阶常微分方程边值问题： d du Lu p qu f x (a, b) dx dx u (a) 0 u(b) 0 ① 构造泛函 1 J (u ) ( Lu, u ) ( f , u ) 2 b b 1 b d du 2 p udx qu dx fudx a a 2 a dx dx 1 b ( pu2 qu 2 2 fu)dx 2 a b du dv 引入泛函算子 (u, v) [ p quv]dx a dx dx 1 则 J (u ) (u, u ) ( f , u ) 2
③ 变分原理（变分问题与边值问题的等价性） 设 f C 0 (I ) ， u* C 2 是边值问题
d du Lu p qu f dx dx u (a) 0, u(b) 0 x ( a, b)
的解，则 u* 使 J(u) 达到极小值； 1 反之，若 u C 2 H E 使 J(u) 达到极小值， 则 u* 是边值问题的解。 1 其中 J (u ) (u, u ) ( f , u ) 2 u (a) 0是强制边界条件，u(b) 0 是自然边 界条件，区别这两类边界条件在用有限元方法求 解边值问题时很重要。
（3）虚功原理
d du Lu p qu f x (a, b) 对两点边值问题： dx dx u (a) 0, u(b) 0 1 设 v HE ，以v乘方程两端，沿[a,b]积分， 并利用 v(a) 0, v(b) 0 ，得变分方程 1 (u, v) ( f , v) 0 v H E b du dv 其中 (u, v) [ p quv]dx a dx dx 虚功原理 2 u C 设 * ，则 u*是边值问题解的充要条件是： 1 u* H E ，且满足变分方程： 1 v H 对任意 (u* , v) ( f , v) 0 E 在力学里， (u, v) ( f , v)表示虚功
② 变分问题
与前述二阶常微分方程边值问题相应的变分 问题是
1 求 u H E ，使 J (u ) min1 J (u )
uH E
1 其中 H E [a, b] {u ( x) u ( x) H 1[a, b], u (a) 0} 1 J (u ) (u, u ) ( f , u ) 2 1 b ( pu2 qu 2 2 fu)dx 2 a
由于实际问题的具体特征、复杂性以 及算法自身的适用范围决定了应用中必须 选择、设计适合于自己特定问题的算法， 因而掌握数值方法的思想至关重要。
要在各种可能的求解方法中找到一种 统一地适用于计算材料学领域（或其它领 域）的理想方法，一般是不现实的。 任何模拟方法，都必须在最佳计算速 度和数值精度之间寻找平衡点。
uC0
在数学上，要将某个微分方程的定解问题 转化为一个变分问题求解，必须针对已给的定 解问题构造一个相应的泛函，并证明定解问题 的解与泛函极值问题的解等价。 有限元方法正是利用这种等价性（边值问 题与变分问题的等价性），先将微分方程定解 问题转化为变分问题（或变分方程）的求解问 题，然后再设法近似求解变分问题（或变分方 程）。
平衡原理 求弦的平衡位置归结为求解两点边值问题： Tu( x) f ( x) 0 x l 其中T是弦的张力。 u (0) 0 u (l ) 0 极小位能原理： 弦的平衡位置 (记为 u u ( x))将在满足边 值条件 u(0)=0，u(l)=0 的一切可能位置中，使 位能取极小值。 设弦处于某一位置u=u(x)，可得到其总位能为 1 l J (u ) (Tu2 2uf )dx 2 0 弦的平衡位置 u u ( x)是下列变分问题的解 J (u ) min2 J (u )
变分原理：
设矩阵A对称正定，则下列两个命题等价： n (a) 求 x0 R ，使 J ( x0 ) minn J ( x ) xR 1 其中 J ( x ) ( Ax , x ) (b , x ) 2 (b) x0 是方程 Ax b 的解 上述两个例子表明： 求二次函数的极小值问题和求线性代数 方程（组）的解是等价的。
偏微分方程数值方法
偏微分方程数值方法偏微分方程（Partial Differential Equation, PDE）是数学中的一种重要的方程类型，它描述了一个函数的多个变量的变化关系。
解决偏微分方程的数值方法在科学和工程领域有着广泛的应用。
本文将介绍几种常见的偏微分方程数值方法，并对其进行详细阐述。
1. 差分法（Finite Difference Method）：差分法是最早也是最直接的一种数值方法，它基于连续函数在一些点的导数可以用它的前向、后向或中心的差商来近似的思想。
偏微分方程的差分格式包括向前差分法、向后差分法和中心差分法等。
对于二维的偏微分方程，可以采用网格化的方式将空间离散化，然后利用差分法进行近似求解。
2. 有限元法（Finite Element Method）：有限元法是一种基于原始形式或变分形式对偏微分方程进行离散化的方法。
在有限元法中，将求解域分割成许多小的、简单的几何单元，然后在每个单元上构建近似解函数和试验函数。
通过构建弱形式并应用基本的变分原理，可以得到离散化的方程组，并通过求解这个方程组来得到数值解。
3. 有限差分法（Finite Difference Method）：有限差分法是一种将连续的偏微分方程离散化成差分方程的方法。
它与差分法的主要区别在于有限差分法不需要对求解域进行网格化，而是直接在连续的求解域上进行离散化。
将偏微分方程中的导数通过差商来近似，然后通过求解离散化的差分方程来得到数值解。
4. 有限体积法（Finite Volume Method）：有限体积法是一种将偏微分方程离散化为离散体积元的方法。
在有限体积法中，将求解域划分成离散的控制体积，然后通过对控制体积的积分运算，将偏微分方程转化为离散的代数方程组。
然后通过求解得到的代数方程组，可以得到数值解。
以上介绍的只是几种常见的偏微分方程数值方法，实际上还有很多其他的方法，如边界元法（Boundary Element Method）、谱方法（Spectral Method）、逆问题方法（Inverse Problem Method）等。
偏微分方程的数值方法
偏微分方程的数值方法偏微分方程是描述自然界许多现象的一种数学模型，它包含多个独立变量，并且方程中的未知函数同时取决于这些变量。
偏微分方程的数值方法是一种求解这类方程的途径，它通过将连续的方程转化为离散的方程，从而使得问题成为一个适用于计算机求解的形式。
本文将介绍几种常用的偏微分方程数值方法。
1. 有限差分法 (Finite Difference Method)有限差分法是最常用的偏微分方程数值方法之一、它将连续的偏微分方程转化为离散的差分方程，通过计算差分方程的近似解来获得原方程的数值解。
在有限差分法中，首先将空间域离散化成网格，再将时间域离散化成步长。
通过近似替代偏微分方程中的导数，将方程转化为差分方程。
通过求解差分方程的解，可以得到偏微分方程的数值解。
2. 有限元法 (Finite Element Method)有限元法是另一种常用的偏微分方程数值方法。
它将连续的偏微分方程转化为离散的代数方程，通过求解代数方程来获得原方程的数值解。
在有限元法中，首先将空间域离散化成有限个小区域，称为有限元。
然后通过选取适当的试探函数和权重函数在每个有限元内部进行插值。
通过将插值函数带入原方程，使用变分原理和加权残差法推导出离散的代数方程。
再通过求解代数方程组的解来得到偏微分方程的数值解。
3. 边界元法 (Boundary Element Method)边界元法也是一种常用的偏微分方程数值方法。
它将连续的偏微分方程转化为边界上的积分方程，通过求解积分方程来获得原方程的数值解。
在边界元法中，将问题的物理域分为两个区域：内域和外域。
通过在内域内求解偏微分方程，得到内域的数值解。
然后通过边界条件将内域的解扩展到整个物理域的边界上。
最后将边界上的积分方程转化为代数方程组，并求解之得到最终的数值解。
4. 谱方法 (Spectral Method)谱方法是一种高精度的偏微分方程数值方法，它同时利用了空间域和频率域的特性。
偏微分方程的数值方法
偏微分方程的数值方法偏微分方程（Partial Differential Equations，简称PDEs）是数学中研究的重要分支，广泛应用于物理学、工程学等领域中。
由于一些复杂的PDEs难以找到解析解，因此需要借助数值方法进行求解。
本文将介绍偏微分方程的数值解法，包括有限差分法、有限元法和谱方法等。
一、有限差分法（Finite Difference Method）有限差分法是解偏微分方程最常用的数值方法之一。
它将偏微分方程中的导数用差商来近似，将空间离散成若干个小区间和时间离散成若干个小时间步长。
通过求解离散化后的代数方程，可以得到原偏微分方程的数值解。
以二维的泊松方程为例，偏微分方程可以表示为：∂²u/∂x² + ∂²u/∂y² = f(x, y)其中，u(x, y)为未知函数，f(x, y)为已知函数。
我们可以将空间离散成Nx × Ny个小区间，时间离散成Nt个小时间步长。
利用中心差分法可以近似表示导数，我们可以得到离散化的代数方程组。
二、有限元法（Finite Element Method）有限元法是一种重要的数值解PDEs的方法。
它将求解区域离散化成一系列的单元，再通过插值函数将每个单元上的未知函数近似表达。
然后，利用加权残差方法，将PDEs转化成代数方程组。
在有限元法中，采用形函数来近似未知函数。
将偏微分方程转化为弱形式，通过选取适当的形函数和权函数，可以得到离散化后的代数方程组。
有限元法适用于求解各种各样的偏微分方程，包括静态和动态、线性和非线性、自由边界和固定边界等问题。
三、谱方法（Spectral Method）谱方法是一种基于特殊函数（如正交多项式）的数值方法，用于解PDEs。
谱方法在求解偏微分方程时，利用高阶连续函数拟合初始条件和边界条件，通过调整特殊函数的系数来近似求解解析解。
谱方法具有高精度和快速收敛的特点，适用于各种偏微分方程求解。
数值模拟偏微分方程的三种方法：FDM、FEM及FVM
数值模拟偏微分方程的三种方法：FDM、FEM及FVM偏微分方程数值模拟常用的方法主要有三种：有限差分方法（FDM）、有限元方法（FEM）、有限体积方法（FVM），本文将对这三种方法进行简单的介绍和比较。
有限差分方法有限差分方法(Finite Difference Methods)是数值模拟偏微分方程最早采用的方法，至今仍被广泛运用。
该方法包括区域剖分和差商代替导数两个过程。
具体地，首先将求解区域划分为差分网格，用有限个网格节点代替连续的求解区域。
其次，利用Taylor级数展开等方法将偏微分方程中的导数项在网格节点上用函数值的差商代替来进行离散，从而建立以网格节点上的值为未知量的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。
差商代替导数后的格式称为有限差分格式，从格式的精度来考虑，有一阶格式、二阶格式和高阶格式。
从差分的空间离散形式来考虑，有中心格式和迎风格式。
对于瞬态方程，考虑时间方向的离散，有显格式、隐格式、交替显隐格式等。
目前常见的差分格式，主要是以上几种格式的组合，不同的组合构成不同的差分格式。
差分方法主要适用于结构网格，网格的步长一般根据问题模型和Courant稳定条件来决定。
请输入标题有限元方法(Finite Element Methods)的基础是变分原理和分片多项式插值。
该方法的构造过程包括以下三个步骤。
首先，利用变分原理得到偏微分方程的弱形式(利用泛函分析的知识将求解空间扩大)。
其次，将计算区域划分为有限个互不重叠的单元(三角形、四边形、四面体、六面体等)。
再次，在每个单元内选择合适的节点作为求解函数的插值点，将偏微分方程中的变量改写成由各变量或其导数的节点值与所选用的分片插值基函数组成的线性表达式，得到微分方程的离散形式。
利用插值函数的局部支集性质及数值积分可以得到未知量的代数方程组。
有限元方法有较完善的理论基础，具有求解区域灵活(复杂区域)、单元类型灵活(适于结构网格和非结构网格)、程序代码通用(数值模拟软件多数基于有限元方法)等特点。
《分数阶偏微分方程的几类有限元方法研究》范文
《分数阶偏微分方程的几类有限元方法研究》篇一一、引言分数阶偏微分方程（Fractional Partial Differential Equations，FPDEs）在物理、工程、生物医学等多个领域中有着广泛的应用。
由于分数阶导数引入了额外的非局部性，这类方程的求解比传统的偏微分方程更为复杂。
近年来，有限元方法作为一种有效的数值计算工具，被广泛应用于分数阶偏微分方程的求解中。
本文将介绍几类有限元方法在分数阶偏微分方程中的应用和研究成果。
二、有限元方法概述有限元方法是一种数值计算方法，它将连续的求解区域离散成有限个单元，通过求解一系列代数方程来得到问题的近似解。
在处理分数阶偏微分方程时，有限元方法具有较好的灵活性和适用性。
根据不同的离散策略和近似函数空间，有限元方法可以分为多种类型。
三、几类有限元方法研究1. 常规有限元方法常规有限元方法在处理分数阶偏微分方程时，主要采用线性或高阶多项式作为近似函数空间。
该方法具有较高的计算精度和稳定性，但当问题规模较大时，计算量较大。
针对这一问题，研究者们提出了一些优化策略，如采用稀疏矩阵技术、多尺度分析等来降低计算成本。
2. 谱有限元方法谱有限元方法是一种基于谱逼近的有限元方法，它采用正交多项式作为近似函数空间。
该方法具有较高的计算精度和收敛速度，尤其适用于处理高阶和弱奇异性的问题。
在处理分数阶偏微分方程时，谱有限元方法能够得到较优的近似解。
3. 径向基函数有限元方法径向基函数有限元方法是一种基于径向基函数的有限元方法，它具有良好的局部性和适应性。
在处理分数阶偏微分方程时，该方法能够有效地捕捉到问题的非局部特性，从而得到较为准确的近似解。
此外，径向基函数有限元方法还具有较好的稳定性和收敛性。
四、应用实例分析本部分将通过具体的应用实例来展示上述几类有限元方法在处理分数阶偏微分方程时的应用效果。
以某物理问题为例，分别采用常规有限元方法、谱有限元方法和径向基函数有限元方法进行求解，并对比分析各种方法的计算精度、稳定性和计算成本等方面的表现。
第六章 第一讲 偏微分方程的有限元法
5.1 泛函与变分原理
数 学 物 理 方 法
5.1.2 变分法
研究泛函极值问题的方法可以归为两类：一类叫直接法，
即直接分析所提出的问题；另一类叫间接法，即把问题 转化为求解微分方程．为讨论间接方法，先介绍变分和 泛函的变分．
Beijing Institute of Technology Deng, Junjun
5.1 泛函与变分原理
数 学 物 理 方 法
类似的例子还可以举出很多，例如，闭合曲线围成的
面积，平面曲线绕固定轴而生成的旋转体积或表面积，
等等，它们也都确定了各自的泛函关系。
这里的函数集合，即泛函的定义域，通常包含要求
y(x)满足一定的边界条件，并且具有连续的二阶导数， 这样的y(x)称为可取函数。
极小值充分条件： J y0 变分： j
j 0 0
0
0
j 0 0
0 是 J 在 y0 处的变分，记为
J
16/76
Beijing Institute of Technology Deng, Junjun
dengjunjun@
15/76
5.1 泛函与变分原理
数 学 物 理 方 法
泛函极值问题转化为一般函数的 j() 极值问题，即：
j J y min!
当 =0时泛函取得极小值J(y0)，根据微积分学可知，
泛函在 y0 取得极值的必要条件是
J y0 0 0
变分运算的几条简单法则：
Beijing Institute of Technology Deng, Junjun
dengjunjun@
20/54
5.1 泛函与变分原理
《分数阶偏微分方程的几类有限元方法研究》范文
《分数阶偏微分方程的几类有限元方法研究》篇一一、引言分数阶偏微分方程（Fractional Partial Differential Equations，FPDEs）在物理、工程、生物医学等多个领域有着广泛的应用。
由于分数阶导数能够更好地描述复杂系统的非局部特性，因此研究分数阶偏微分方程的数值解法具有重要的理论和实践意义。
有限元方法作为一种有效的数值计算工具，在处理分数阶偏微分方程问题上表现出了突出的优势。
本文旨在综述分数阶偏微分方程的几类有限元方法的研究进展，探讨各自的优缺点和适用场景。
二、Caputo-Fabrizio 分数阶导数与偏微分方程的引入首先介绍Caputo-Fabrizio 分数阶导数的定义及其性质，这是研究分数阶偏微分方程的基础。
接着介绍不同领域中常见的偏微分方程，以及当其与分数阶导数结合时所形成的FPDEs。
这类方程在描述复杂系统的扩散、传播等过程时具有更高的精度和适应性。
三、有限元方法在分数阶偏微分方程中的应用（一）传统有限元方法传统有限元方法在处理FPDEs时，通过将连续的求解域划分为有限个离散单元，将偏微分方程转化为线性方程组进行求解。
对于FPDEs中的空间分数阶导数和/或时间分数阶导数，可以借助数值积分或差分法进行离散化处理。
本节将详细介绍传统有限元方法在FPDEs中的应用及其优缺点。
（二）局部弱解有限元方法局部弱解有限元方法是一种针对FPDEs的特殊有限元方法，其核心思想是利用弱解形式将FPDEs转化为等价的变分问题。
该方法能够有效地降低求解问题的复杂性，并提高求解精度。
本节将详细介绍局部弱解有限元方法的原理和实现过程，以及其在FPDEs中的应用实例。
（三）混合型有限元方法混合型有限元方法是一种结合了传统有限元方法和其他数值方法的混合型数值方法。
针对FPDEs中的不同类型导数（如空间导数和时间导数），可以灵活地选择不同的数值处理方法，以达到更好的求解效果。
本节将详细介绍混合型有限元方法的原理和实现过程，以及其在FPDEs中的实际应用。
偏微分方程的有限元法
利用有限元法求解弹性力学问题的基本步骤包括建立离散化的数学模型、选择合适的有 限元空间、求解离散化的线性方程组等。
传热学问题
传热学中的偏微分方程
描述热传导、对流、辐射等过程的偏微分方程包括热传导 方程、对流方程等，这些方程描述了温度场的变化规律。
有限元法在传热学中的应用
通过将连续的温度场离散化为有限个单元，有限元法能够 求解复杂的传热学问题，如热传导、对流换热、辐射换热 等。
区域离散
将连续的求解区域离散化为有限 个小的子区域，每个子区域称为
一个有限元。
函数近似
在每个有限元上选择适当的基函数 来近似未知函数，基函数的选择取 决于问题的性质和求解精度要求。
离散化方程
根据微分方程和边界条件，建立离 散化的代数方程组，表示为矩阵形 式。
有限元法的求解过程
线性化
将非线性微分方程转化为线性方程组，以便于求 解。
描述流体运动的偏微分方程包括Navier-Stokes方程、Euler方 程等，这些方程描述了流体的速度、压力、密度等物理量的变
化规律。
有限元法在流体动力学中的应用
通过将连续的流体域离散化为有限个单元，有限元法能够 求解复杂的流体动力学问题，如湍流、非牛顿流体等。
求解方法
利用有限元法求解流体动力学问题的基本步骤包括建立离散化 的数学模型、选择合适的有限元空间、求解离散化的线性方程
组等。
弹性力学问题
弹性力学中的偏微分方程
描述弹性物体变形的偏微分方程包括弹性力学的基本方程、Mindlin-Reissner方程等， 这些方程描述了弹性体的应力、应变等物理量的变化规律。
有限元法在弹性力学中的应用
通过将连续的弹性体离散化为有限个单元，有限元法能够求解复杂的弹性力学问题，如 非线性弹性、复合材料等。
偏微分方程的数值解法
偏微分方程的数值解法偏微分方程（Partial Differential Equation, PDE）是数学和物理学中的重要概念，广泛应用于工程、科学和其他领域。
在很多情况下，准确解析解并不容易获得，因此需要利用数值方法求解偏微分方程。
本文将介绍几种常用的数值解法。
1. 有限差分法（Finite Difference Method）有限差分法是最常见和经典的数值解法之一。
基本思想是将偏微分方程在求解域上进行离散化，然后用差分近似代替微分运算。
通过求解差分方程组得到数值解。
有限差分法适用于边界条件简单且求解域规则的问题。
2. 有限元法（Finite Element Method）有限元法是适用于不规则边界条件和求解域的数值解法。
将求解域划分为多个小区域，并在每个小区域内选择适当的形状函数。
通过将整个域看作这些小区域的组合来逼近原始方程，从而得到一个线性代数方程组。
有限元法具有较高的灵活性和适用性。
3. 有限体积法（Finite Volume Method）有限体积法是一种较新的数值解法，特别适用于物理量守恒问题。
它通过将求解域划分为多个控制体积，并在每个体积内计算守恒量的通量，来建立离散的方程。
通过求解这个方程组得到数值解。
有限体积法在处理守恒律方程和非结构化网格上有很大优势。
4. 局部网格法（Local Grid Method）局部网格法是一种多尺度分析方法，适用于具有高频振荡解的偏微分方程。
它将计算域划分为全局细网格和局部粗网格。
在全局细网格上进行计算，并在局部粗网格上进行局部评估。
通过对不同尺度的解进行耦合，得到更精确的数值解。
5. 谱方法（Spectral Method）谱方法是一种基于傅里叶级数展开的高精度数值解法。
通过选择适当的基函数来近似求解函数，将偏微分方程转化为代数方程。
谱方法在处理平滑解和周期性边界条件的问题上表现出色，但对于非平滑解和不连续解的情况可能会遇到困难。
6. 迭代法（Iterative Method）迭代法是一种通过多次迭代来逐步逼近精确解的求解方法。
偏微分方程的数值方法
偏微分方程的数值方法偏微分方程是包含多个变量的方程，其中包含偏导数，用于描述多变量函数的变化规律。
解决偏微分方程的数值方法是一种近似求解的方式，主要用于那些无法通过解析方法求得精确解的方程。
本文将介绍几种常见的偏微分方程数值方法。
一、有限差分方法（Finite Difference Method）有限差分方法是求解偏微分方程的一种常见数值方法。
其基本思想是将偏微分方程中的各个偏导数用有限差分的形式来近似表示。
将方程中的空间变量和时间变量分别离散化，即将空间和时间分成一系列的网格点，根据差分近似的原理，将方程转化为一系列的代数方程，然后通过迭代计算求解。
常用的有限差分方法包括显式差分法、隐式差分法和Crank-Nicolson差分法。
二、有限元方法（Finite Element Method）有限元方法是求解偏微分方程的一种常见数值方法。
其基本步骤是将求解区域划分为多个小区域（要素），然后根据偏微分方程的特性构造适当的有限元模型，并建立离散化方程，最后通过求解线性代数方程组来获得数值解。
有限元方法具有较高的灵活性和通用性，对各种不规则边界条件和复杂几何形状的求解问题具有很好的适应性。
三、谱方法（Spectral Method）谱方法是求解偏微分方程的一种高精度数值方法。
其基本思想是将待求解的函数表示为一系列基函数的线性组合，而后通过合适的基函数和求解区域内的截断误差最小化，获得函数近似解。
谱方法对于光滑的解具有高精度的逼近性能和收敛性，常用的基函数有Chebyshev多项式、Legendre多项式和傅立叶级数等。
四、边界元方法（Boundary Element Method）边界元方法是求解偏微分方程的另一种常见数值方法。
其基本思想是将区域内的偏微分方程问题转化为对区域边界上的积分方程的求解问题。
通过将边界上的未知函数值和边界上的迹值引入，并应用格林第二定理，将区域内的偏微分方程问题转化为一系列的线性代数方程组，进而获得数值解。