BDF(MD)

打印 上一主题 下一主题

主题 1008|帖子 1008|积分 3028


这句话的语法结构和意思可以分解如下:
语法结构:

  • Forano.d.e. - 这部分看起来像是一个不完备的短语或术语的缩写,大概是指“对于常微分方程(For ordinary differential equations)”,但由于它不完备,这里存在一定的语法错误或不明白性。
  • y′=f(t,y) - 这是一个数学表达式,表现一个常微分方程,其中 y′ 表现未知函数 y 对变量 t 的导数,f(t,y) 是一个已知函数,依靠于 t 和 y。
  • is constructed by approximating - 这部分是被动语态,表现“通过近似来构造”。
  • the derivative of y at the node tn+7 - 这部分指的是在节点 tn+7 处 y 的导数。
  • in the relation y′(tn+7)=f(tn+7,y(tn+7)) - 这部分说明了导数的具体关系,即 y 在 tn+7 处的导数等于 f 在 t=tn+7 和 y=y(tn+7) 处的值。
  • by the derivative P′7(tn+7) of the interpolating polynomial - 这部分说明是通过插值多项式的导数 P′7(tn+7) 来近似 y′(tn+7)。
    意思:
    这句话的意思是,对于一个常微分方程 y′=f(t,y),我们通过近似的方法来构造在节点 tn+7 处 y 的导数。这个近似是通过使用插值多项式 P7(t) 在 tn+7 处的导数 P′7(tn+7) 来实现的。具体来说,我们将 y 在 tn+7 处的导数 y′(tn+7) 近似为 f(tn+7, y(tn+7)),这个值是通过插值多项式的导数 P′7(tn+7) 来计算的。
    总结来说,这句话形貌了一种数值方法,用于近似求解常微分方程在特定节点处的导数。
BDF如何运作

BDF(Backward Differentiation Formula)方法是一类隐式线性多步法,用于求解常微分方程初值题目。BDF7是BDF方法的一种,即七步BDF方法。下面我将详细形貌BDF7方法的运作方式:
根本概念

BDF方法是一类多步法,它们使用过去多个时间步的解来计算当前时间步的解。BDF7方法特别地使用前七个时间步的信息。
BDF7公式

BDF7方法的公式可以写成如下情势:
                                                    a                                           −                                  7                                                            y                                           n                                  −                                  7                                                 +                                       a                                           −                                  6                                                            y                                           n                                  −                                  6                                                 +                                       a                                           −                                  5                                                            y                                           n                                  −                                  5                                                 +                                       a                                           −                                  4                                                            y                                           n                                  −                                  4                                                 +                                       a                                           −                                  3                                                            y                                           n                                  −                                  3                                                 +                                       a                                           −                                  2                                                            y                                           n                                  −                                  2                                                 +                                       a                                           −                                  1                                                            y                                           n                                  −                                  1                                                 +                                       a                               0                                                 y                               n                                      =                            h                            f                            (                                       t                               n                                      ,                                       y                               n                                      )                                  a_{-7} y_{n-7} + a_{-6} y_{n-6} + a_{-5} y_{n-5} + a_{-4} y_{n-4} + a_{-3} y_{n-3} + a_{-2} y_{n-2} + a_{-1} y_{n-1} + a_0 y_n = h f(t_n, y_n)                     a−7​yn−7​+a−6​yn−6​+a−5​yn−5​+a−4​yn−4​+a−3​yn−3​+a−2​yn−2​+a−1​yn−1​+a0​yn​=hf(tn​,yn​)
其中,( y_n ) 是在时间 ( t_n ) 的近似解,( h ) 是时间步长,( f(t, y) ) 是微分方程的右侧函数,( a_{-7}, a_{-6}, …, a_0 ) 是BDF方法的系数。
运作步调


  • 初始化:起首,需要初始条件 ( y_0, y_1, …, y_6 )。这些初始条件可以通过其他数值方法(如Runge-Kutta方法)获得。
  • 创建方程:在每一个时间步 ( t_n ),使用BDF7公式创建关于 ( y_n ) 的方程。这个方程是隐式的,由于它包罗未知的 ( y_n )。
  • 求解非线性方程:由于BDF7公式是隐式的,我们需要解一个非线性方程来找到 ( y_n )。这通常通过牛顿法或拟牛顿法等迭代方法来完成。
  • 更新解:一旦找到 ( y_n ),就可以更新解序列,将 ( y_n ) 添加到解序列中,并准备进入下一个时间步。
  • 重复过程:重复步调2到4,直到到达预定的最终时间或满意其他制止条件。
系数

BDF7方法的系数 ( a_{-7}, a_{-6}, …, a_0 ) 是固定的,它们是由BDF方法的递推关系确定的。这些系数确保了方法的稳定性和收敛性。
特点



  • 隐式性:BDF方法隐式地依靠于当前时间步的解,这使得它们对于刚性方程特别有效。
  • 多步性:BDF方法使用多个时间步的信息,这提高了方法的精度。
  • 稳定性:BDF方法通常比显式方法更稳定,尤其是对于大时间步长。
注意事项

在实际应用中,BDF7方法大概需要与变量时间步长策略相结合,以保持数值稳定性。此外,由于其隐式性质,BDF7方法的每一步都需要求解非线性方程,这大概会增加计算成本。
BDF7方法在许多科学和工程范畴中都有应用,尤其是在求解涉及刚性系统的常微分方程时。
详细

好的,我将实验用更简朴和详细的方式来表明BDF7(七步向后差分公式)是如何运作的。
根本概念

起首,我们需要明白常微分方程初值题目(IVP):
假设我们有一个常微分方程:
[ y’(t) = f(t, y(t)), ]
并且给定初始条件 ( y(t_0) = y_0 )。
我们的目的是找到这个方程在 ( t > t_0 ) 时的解 ( y(t) )。
BDF7的多步性质

BDF7是一种多步方法,意味着在计算当前时间步 ( t_n ) 的解 ( y_n ) 时,它会使用之前时间步的解 ( y_{n-1}, y_{n-2}, …, y_{n-6} )。
运作步调


  • 初始化

    • 使用其他数值方法(如Runge-Kutta方法)来计算前六个时间步的解 ( y_0, y_1, …, y_6 )。

  • 创建BDF7公式

    • BDF7公式可以写成如下情势:
                                                                      α                                     0                                                           y                                     n                                              +                                               α                                     1                                                           y                                                   n                                        −                                        1                                                           +                                               α                                     2                                                           y                                                   n                                        −                                        2                                                           +                                  .                                  .                                  .                                  +                                               α                                     6                                                           y                                                   n                                        −                                        6                                                           =                                  β                                  ⋅                                  f                                  (                                               t                                     n                                              ,                                               y                                     n                                              )                                          \alpha_0 y_n + \alpha_1 y_{n-1} + \alpha_2 y_{n-2} + ... + \alpha_6 y_{n-6} = \beta \cdot f(t_n, y_n)                           α0​yn​+α1​yn−1​+α2​yn−2​+...+α6​yn−6​=β⋅f(tn​,yn​),
    • 其中,(                                                                     α                                        0                                                  ,                                                   α                                        1                                                  ,                                     .                                     .                                     .                                     ,                                                   α                                        6                                                           \alpha_0, \alpha_1, ..., \alpha_6                              α0​,α1​,...,α6​) 是固定的系数,( \beta ) 是与时间步长 ( h ) 干系的系数。

  • 求解非线性方程

    • 由于 ( y_n ) 出如今等式的两边,我们需要解一个非线性方程来找到 ( y_n )。
    • 这通常通过迭代方法完成,例如牛顿法。我们从一个初始猜测开始,然后迭代直到找到足够靠近的解。

  • 更新解

    • 一旦我们找到了 ( y_n ),我们就将它添加到我们的解序列中,准备进入下一个时间步 ( t_{n+1} )。

  • 重复过程

    • 使用新的 ( y_n ) 和旧的 ( y_{n-1}, …, y_{n-6} ) 来计算 ( y_{n+1} )。
    • 重复这个过程,直到到达预定的最终时间或满意其他制止条件。

示例

假设我们有以下BDF7公式:
                                                    y                               n                                      −                                       49                               240                                                 y                                           n                                  −                                  1                                                 +                                       6                               25                                                 y                                           n                                  −                                  2                                                 −                                       867                               1200                                                 y                                           n                                  −                                  3                                                 +                                       301                               240                                                 y                                           n                                  −                                  4                                                 −                                       33                               40                                                 y                                           n                                  −                                  5                                                 +                                       7                               10                                                 y                                           n                                  −                                  6                                                 =                                       7                               10                                      h                            f                            (                                       t                               n                                      ,                                       y                               n                                      )                            .                                  y_n - \frac{49}{240} y_{n-1} + \frac{6}{25} y_{n-2} - \frac{867}{1200} y_{n-3} + \frac{301}{240} y_{n-4} - \frac{33}{40} y_{n-5} + \frac{7}{10} y_{n-6} = \frac{7}{10} h f(t_n, y_n).                     yn​−24049​yn−1​+256​yn−2​−1200867​yn−3​+240301​yn−4​−4033​yn−5​+107​yn−6​=107​hf(tn​,yn​).
在这个公式中,                                                   α                               0                                      =                            1                            ,                                       α                               1                                      =                            −                                       49                               240                                      ,                                       α                               2                                      =                                       6                               25                                      ,                                  \alpha_0 = 1 , \alpha_1 = -\frac{49}{240} , \alpha_2 = \frac{6}{25} ,                     α0​=1,α1​=−24049​,α2​=256​,等等。
总结

BDF7方法的关键在于它使用过去七个时间步的解来预测当前时间步的解。这种方法对于解决刚性方程特别有效,由于它能够保持稳定性。然而,由于其隐式性质,每一步都需要解一个非线性方程,这大概会使计算变得复杂。
chu

这段文字讨论了初值题目(Initial Value Problem, IVP)和数值方法中的七步向后差分公式(Backward Differentiation Formula, BDF)。具体来说,它涉及如何通过拉格朗日插值多项式来近似求解常微分方程(Ordinary Differential Equation, ODE)。
初值题目(IVP)
起首定义了一个初值题目:
                                                                                       {                                                                                                                                        u                                                       ′                                                                      (                                                    t                                                    )                                                    +                                                    A                                                    u                                                    (                                                    t                                                    )                                                    =                                                    f                                                    (                                                    t                                                    )                                                    ,                                                                                                                                              0                                                    <                                                    t                                                    <                                                    T                                                    ,                                                                                                         u                                                    (                                                    0                                                    )                                                    =                                                                       u                                                       0                                                                      .                                                                                                                                                       (1)                                                       \begin{cases} u'(t) + Au(t) = f(t), & 0 < t < T, \ u(0) = u^0. \end{cases} \tag{1}                     {u′(t)+Au(t)=f(t),​0<t<T, u(0)=u0.​(1)
这里 ( u(t) ) 是未知函数,( A ) 和 ( f(t) ) 是已知函数或常数,初始条件是 ( u(0) = u^0 )。
拉格朗日插值多项式
接下来引入了拉格朗日插值多项式 ( P_7 ),用于在节点 ( t_n, t_{n+1}, \ldots, t_{n+7} ) 上对某个函数 ( y ) 举行插值。
常微分方程(ODE)
考虑一个常微分方程:
[
y’ = f(t,y),
]
这个方程表现的是 ( y ) 对时间 ( t ) 的导数等于一个关于 ( t ) 和 ( y ) 的函数 ( f(t,y) )。
七步BDF方法
文本提到使用拉格朗日插值多项式的导数 ( P’7(t{n+7}) ) 来近似计算 ( y ) 在节点 ( t_{n+7} ) 处的导数 ( y’(t_{n+7}) )。这种方法称为七步BDF方法,其情势为:
                                                                                                     ∑                                                       i                                           =                                           0                                                      7                                                                a                                        i                                                                y                                                       n                                           +                                           i                                                                =                                     τ                                     f                                     (                                                   t                                                       n                                           +                                           7                                                                ,                                                   y                                                       n                                           +                                           7                                                                )                                     .                                                                            (2)                                                       \sum_{i=0}^{7} a_i y^{n+i} = \tau f(t_{n+7}, y^{n+7}). \tag{2}                     i=0∑7​ai​yn+i=τf(tn+7​,yn+7).(2)
其中 ( a_i ) 是系数,( \tau ) 是步长,( y^{n+i} ) 表如今第 ( n+i ) 步时的 ( y ) 值。
总结
这段文字主要介绍了如何使用拉格朗日插值多项式和七步BDF方法来近似求解常微分方程。通过这种方式,可以在一系列离散点上逐步逼近解曲线。
这句话的意思是在节点                                              t                                       n                               +                               7                                                 t_{n+7}                  tn+7​处,通过拉格朗日插值多项式的导数                                             P                            7                            ′                                  (                                   t                                       n                               +                               7                                            )                              P_7^{\prime}(t_{n+7})                  P7′​(tn+7​)来近似构造                                   y                              y                  y的导数                                             y                            ′                                  (                                   t                                       n                               +                               7                                            )                         =                         f                         (                                   t                                       n                               +                               7                                            ,                         y                         (                                   t                                       n                               +                               7                                            )                         )                              y'(t_{n+7})=f(t_{n+7},y(t_{n+7}))                  y′(tn+7​)=f(tn+7​,y(tn+7​))的关系。
在这句话中,“in the relation”表现“在……关系中”。它用来引出或强调某个特定的关系表达式或等式。在这个上下文中,它连接了两个部分:一个是关于 ( y ) 的导数的表达式                                                    y                               ′                                      (                                       t                                           n                                  +                                  7                                                 )                            =                            f                            (                                       t                                           n                                  +                                  7                                                 ,                            y                            (                                       t                                           n                                  +                                  7                                                 )                            )                                  y'(t_{n+7}) = f(t_{n+7}, y(t_{n+7}))                     y′(tn+7​)=f(tn+7​,y(tn+7​)),另一个是用于近似的导数 ( P_7’(t_{n+7}) )。因此,“in the relation”在这里的作用是指出这两个表达式之间的关系。
文中的两个表达式分别是 ( y’ = f(t, y) ) 和 ( y’(t_{n+7}) = f(t_{n+7}, y(t_{n+7})) )。它们之间的关系在于后者是对前者在特定节点 ( t_{n+7} ) 处的特例化。也就是说,( y’(t_{n+7}) = f(t_{n+7}, y(t_{n+7})) ) 是对一般情势的微分方程 ( y’ = f(t, y) ) 在时间点 ( t_{n+7} ) 的应用。这种特例化是通过使用插值多项式的导数 ( P_7’(t_{n+7}) ) 来近似实现的。
这段文字形貌的是如何构造一个常微分方程(Ordinary Differential Equation, ODE)的数值解法。下面是对这段话的详细表明:
对于一个给定的常微分方程:
[ y’ = f(t, y), ]
我们的目的是找到函数 ( y(t) ) 的数值解,其中 ( y’ ) 表现 ( y ) 关于 ( t ) 的导数,而 ( f(t, y) ) 是已知的一个函数,它给出了 ( y ) 的导数如何依靠于 ( t ) 和 ( y ) 的值。
在数值解法中,我们通常在一系列离散的时间点上求解这个方程。这里的 ( tn+7 ) 是这些离散时间点中的一个,表如今第 ( n+7 ) 个时间点。
如今,我们要在时间点 ( tn+7 ) 处构造 ( y ) 的导数的近似。根据微分方程,我们有:
[ y’(tn+7) = f(tn+7, y(tn+7)). ]
这是在 ( tn+7 ) 处的准确关系,它告诉我们 ( y ) 在 ( tn+7 ) 处的导数等于 ( f ) 在 ( tn+7 ) 和 ( y(tn+7) ) 处的值。
然而,在实际的数值计算中,我们通常无法直接得到 ( y(tn+7) ) 和 ( y’(tn+7) ) 的准确值。因此,我们需要使用数值方法来近似它们。这里提到的插值多项式 ( P_7(t) ) 是一个用来近似 ( y(t) ) 的多项式,它是基于之前计算得到的 ( y ) 的值构造的。
为了近似 ( y’(tn+7) ),我们使用插值多项式 ( P_7(t) ) 在 ( tn+7 ) 处的导数 ( P_7’(tn+7) )。这个导数 ( P_7’(tn+7) ) 是通过以下步调构造的:

  • 起首,我们使用之前计算得到的 ( y ) 的值来构造拉格朗日插值多项式 ( P_7(t) )。
  • 然后,我们计算这个插值多项式在 ( tn+7 ) 处的导数 ( P_7’(tn+7) )。
  • 最后,我们将 ( P_7’(tn+7) ) 作为 ( y’(tn+7) ) 的近似值。
    这样,我们就得到了在 ( tn+7 ) 处 ( y ) 的导数的近似值,可以用于进一步的数值计算,比如求解常微分方程的数值解。
jutiliz

固然可以。让我们通过一个具体的例子来表明如何使用拉格朗日插值多项式的导数来近似常微分方程在特定节点处的导数。
假设我们有以下常微分方程的初值题目:
y’ = f(t, y) = y + t,
y(0) = 1,
我们想要在 t = tn7 = 0.5 处近似 y’(0.5)。
起首,我们需要一些离散的点来构造插值多项式。假设我们已经通过某种数值方法(比如欧拉法或龙格-库塔法)计算了以下点:
t0 = 0, y0 = 1,
t1 = 0.1, y1 = 1.1,
t2 = 0.2, y2 = 1.22,
t3 = 0.3, y3 = 1.361,
t4 = 0.4, y4 = 1.5264.
如今,我们想要构造一个七阶拉格朗日插值多项式 P7(t) 来近似 y(t) 在这些点上的值,并计算其在 t = 0.5 处的导数。
拉格朗日插值多项式的情势是:
P7(t) = Σ[yi * li(t)],
其中 li(t) 是第 i 个拉格朗日基多项式,定义为:
li(t) = Π[(t - tj) / (ti - tj)], j ≠ i.
对于 i = 0 到 4,我们可以计算出每个 li(t) 的值,然后构造 P7(t)。但是,由于我们只有 5 个点,我们实际上只能构造一个五阶多项式。为了简化,我们假设我们已经有足够的数据点来构造七阶多项式。
接下来,我们需要计算 P7(t) 在 t = 0.5 处的导数 P7’(t)。拉格朗日插值多项式的导数可以表现为:
P7’(t) = Σ[yi * li’(t)],
其中 li’(t) 是 li(t) 的导数,计算为:
li’(t) = Σ[1 / (t - tj)] - Σ[1 / (ti - tj)].
我们可以计算每个 li’(0.5) 的值,然后构造 P7’(0.5)。
以下是一个简化的计算过程:
假设我们已经有了 li(0.5) 和 li’(0.5) 的值,我们可以计算 P7(0.5) 和 P7’(0.5):
P7(0.5) = y0 * l0(0.5) + y1 * l1(0.5) + y2 * l2(0.5) + y3 * l3(0.5) + y4 * l4(0.5),
P7’(0.5) = y0 * l0’(0.5) + y1 * l1’(0.5) + y2 * l2’(0.5) + y3 * l3’(0.5) + y4 * l4’(0.5).
最后,P7’(0.5) 就是 y’(0.5) 的近似值。
请注意,这个例子是简化的,实际上你需要计算所有的基多项式和它们的导数,这大概相当复杂。在实际应用中,通常会使用数值软件来实行这些计算。
jiuye

学习常微分方程(ODEs)及其数值解法,例如BDF(Backward Differentiation Formulas)方法和时间自顺应算法,可以为你在以下范畴的工作打下坚固的底子:

  • 科学计算与工程仿真:许多工程和物理题目都可以转化为常微分方程举行求解,特别是在流体力学、结构分析、热传导等范畴。
  • 金融数学:金融衍生品定价常常涉及到偏微分方程和常微分方程的求解,把握干系技能有助于从事量化分析和风险管理等工作。
  • 数据科学与机器学习:一些高级的优化算法和动力学系统建模也依靠于常微分方程的知识。
  • 软件开发:开发高性能的计算软件或数值模拟工具,特别是那些需要处理复杂数值计算的范畴。
  • 研究与学术岗位:假如你对科学研究感兴趣,可以在高校或研究机构从事干系的教学和研究工作。
  • 控制系统筹划:在自动化控制系统中,常微分方程用于形貌系统的动态举动,明白这些方程对于筹划和优化控制系统至关紧张。
    把握这些技能不但能提拔你在上述范畴的竞争力,还能为你打开更多跨学科合作的时机。
    学习加权和移位的七步BDF方法(WSBDF7)以及干系的数值方法偏微分方程(PDE)求解技术,可以为你打开许多职业时机,尤其是在科学计算工程仿真数据科学学术研究等范畴。以下是一些具体的职业方向和干系的工作时机:

1. 科学与工程计算范畴



  • 数值算法开发工程师

    • 开发高效的数值算法,用于求解偏微分方程(如热方程、流体动力学方程等)。
    • 工作范畴:航空航天、能源、材料科学、气候模拟等。
    • 公司范例:科研机构(如NASA、CERN)、工程公司(如ANSYS、COMSOL)、科技公司(如Intel、NVIDIA)。

  • 高性能计算(HPC)工程师

    • 使用超级计算机或分布式计算系统,优化数值算法的并行计算性能。
    • 工作范畴:气候预报、分子动力学模拟、核反应堆筹划等。
    • 公司范例:国家实验室(如橡树岭国家实验室)、科技公司(如Google、Amazon AWS)。

  • 仿真软件工程师

    • 开发用于工程仿真的软件工具,如有限元分析(FEA)、计算流体动力学(CFD)等。
    • 工作范畴:汽车、航空航天、电子设备筹划等。
    • 公司范例:ANSYS、西门子、达索系统。


2. 数据科学与机器学习



  • 科学计算与机器学习结合

    • 将数值方法应用于机器学习模子的训练和优化,特别是在处理高维数据或复杂模子时。
    • 工作范畴:金融建模、医疗数据分析、自动驾驶等。
    • 公司范例:科技公司(如Google、Facebook)、金融公司(如高盛、摩根大通)。

  • 数据科学家

    • 使用数值方法分析大规模数据集,构建预测模子或优化算法。
    • 工作范畴:市场营销、供应链管理、医疗诊断等。
    • 公司范例:咨询公司(如麦肯锡、波士顿咨询)、科技公司(如Amazon、Microsoft)。


3. 学术研究与教诲



  • 大学传授或研究员

    • 在数学、应用数学、计算科学等范畴从事教学和研究工作。
    • 研究方向:数值分析、偏微分方程、科学计算等。
    • 工作地点:大学、研究所(如中国科学院、MIT、斯坦福)。

  • 博士后研究员

    • 在数值方法、计算数学等范畴举行深入研究,发表高水平论文。
    • 工作地点:国表里着名大学或研究所。


4. 金融与量化分析



  • 量化分析师

    • 使用数值方法开发金融模子,用于期权定价、风险管理、投资策略优化等。
    • 工作范畴:对冲基金、投资银行、保险公司。
    • 公司范例:高盛、摩根士丹利、文艺复兴科技。

  • 金融工程师

    • 筹划和实现金融产物的定价模子和风险管理工具。
    • 工作范畴:衍生品定价、资产设置、信用风险评估。
    • 公司范例:花旗银行、巴克莱银行、黑石团体。


5. 科技公司与创业



  • 算法工程师

    • 开发高效的数值算法,用于解决实际题目(如图像处理、信号处理、优化题目)。
    • 工作范畴:人工智能、计算机视觉、自然语言处理。
    • 公司范例:Google、Facebook、阿里巴巴、腾讯。

  • 科技创业者

    • 使用数值方法和科学计算技术,开发创新的软件产物或服务。
    • 创业范畴:仿真软件、数据分析工具、科学计算平台。


6. 政府与公共部门



  • 科研机构研究员

    • 在政府资助的科研机构从事数值方法和科学计算的研究。
    • 工作范畴:气候模拟、能源政策、国防科技。
    • 机构范例:中国科学院、美国能源部、欧洲核子研究中心(CERN)。

  • 政策分析师

    • 使用数值模子分析政策影响,为政府决定提供支持。
    • 工作范畴:情况保护、公共卫生、经济发展。
    • 机构范例:世界银行、联合国、国家发改委。


7. 跨学科范畴



  • 生物医学工程师

    • 使用数值方法模拟生物系统(如血流动力学、肿瘤生长模子)。
    • 工作范畴:医疗设备筹划、药物研发、生物信息学。
    • 公司范例:强生、辉瑞、西门子医疗。

  • 气候与能源分析师

    • 使用数值模子分析气候变化、能源供需等题目。
    • 工作范畴:可再生能源、情况政策、能源经济。
    • 公司范例:国际能源署(IEA)、世界资源研究所(WRI)。


焦点技能与竞争力


  • 数学底子

    • 偏微分方程、数值分析、线性代数、优化理论。

  • 编程能力

    • 熟练把握Python、MATLAB、C++等编程语言。

  • 算法筹划与优化

    • 能够筹划高效的数值算法,并优化其计算性能。

  • 科学计算工具

    • 熟悉常用的科学计算库(如NumPy、SciPy、FEniCS)。

  • 跨学科能力

    • 能够将数值方法应用于不同范畴(如物理、生物、金融)。


总结

学习WSBDF7方法和干系的数值分析技术,可以为你打开广泛的职业时机,尤其是在科学计算工程仿真数据科学学术研究等范畴。无论你是想进入工业界、金融界,照旧从事学术研究,这些技能都会为你提供强大的竞争力。
BDF(Backward Differentiation Formulas)是一种用于解决常微分方程(ODEs)或多步法的数值方法。时间自顺应算法则是一类根据误差估计自动调整时间步长的算法。两者之间的接洽在于,时间自顺应算法可以应用于BDF方法中来提高其效率和精度。
具体来说,当使用BDF方法时,可以选择不同的多步模式(如一步、两步等),并且每一步都需要计算多个函数值。然而,在某些情况下,大概并不需要所有这些函数值都能到达相同的精度要求。因此,可以通过引入时间自顺应机制来动态地调整每一步的时间步长,以确保在整个积分过程中保持所需的精度水平。
这种方法的优点是可以镌汰不必要的计算量,从而提高团体效率;同时也可以制止因过小的步长而导致的不稳定性和收敛性题目。此外,由于BDF方法自己具有较高的稳定性,因此在接纳时间自顺应策略时通常能够获得更好的结果。
六七部

六步BDF(Backward Differentiation Formula)和七步BDF是用于求解常微分方程(ODEs)和偏微分方程(PDEs)的数值方法。BDF方法是一类隐式多步法,特别适用于求解刚性微分方程。以下是详细表明:
1. BDF方法的根本概念

BDF方法通过使用当前步和前几步的解来近似微分方程的导数。BDF方法的阶数决定了其精度和稳定性。一般来说,BDF方法的阶数越高,精度越高,但稳定性大概会降低。
2. 六步BDF方法

六步BDF方法使用当前步和前六步的解来近似导数。其一般情势为:
                                                    ∑                                           i                                  =                                  0                                          6                                                 a                               i                                                 y                                           n                                  +                                  i                                                 =                            τ                            f                            (                                       t                                           n                                  +                                  6                                                 ,                                       y                                           n                                  +                                  6                                                 )                                  \sum_{i=0}^{6} a_i y^{n+i} = \tau f(t_{n+6}, y^{n+6})                     i=0∑6​ai​yn+i=τf(tn+6​,yn+6)
其中,( y^{n+i} ) 是在时间步 ( t_{n+i} ) 处的解,( \tau ) 是时间步长,( f ) 是微分方程的右端函数。
六步BDF方法的稳定性较好,适用于求解刚性微分方程。其稳定性地区较大,能够处理较大范围的时间步长。
3. 七步BDF方法

七步BDF方法使用当前步和前七步的解来近似导数。其一般情势为:
                                                    ∑                                           i                                  =                                  0                                          7                                                 a                               i                                                 y                                           n                                  +                                  i                                                 =                            τ                            f                            (                                       t                                           n                                  +                                  7                                                 ,                                       y                                           n                                  +                                  7                                                 )                                  \sum_{i=0}^{7} a_i y^{n+i} = \tau f(t_{n+7}, y^{n+7})                     i=0∑7​ai​yn+i=τf(tn+7​,yn+7)
七步BDF方法的阶数更高,理论上精度也更高。然而,七步BDF方法在稳定性方面存在题目,特别是对于刚性方程,它大概不稳定。
4. 稳定性分析



  • 六步BDF方法:六步BDF方法的稳定性较好,适用于求解刚性微分方程。其稳定性地区较大,能够处理较大范围的时间步长。
  • 七步BDF方法:七步BDF方法在稳定性方面存在题目,特别是对于刚性方程,它大概不稳定。因此,直接使用七步BDF方法大概会导致数值解的不稳定。
5. 加权和移位的七步BDF方法(WSBDF7)

为了降服七步BDF方法的不稳定性,研究者提出了加权和移位的七步BDF方法(WSBDF7)。该方法通过将七步BDF方法与其移位版本举行线性组合,构造出一个稳定的数值方法。具体来说,WSBDF7方法的情势为:
                                         ϑ                                       ∑                                           i                                  =                                  0                                          7                                                 a                               i                                                 y                                           n                                  +                                  i                                                 +                            (                            1                            −                            ϑ                            )                                       ∑                                           i                                  =                                  0                                          7                                                             a                                  ~                                          i                                                 y                                           n                                  +                                  i                                                 =                            ϑ                            τ                            f                            (                                       t                                           n                                  +                                  7                                                 ,                                       y                                           n                                  +                                  7                                                 )                            +                            (                            1                            −                            ϑ                            )                            τ                            f                            (                                       t                                           n                                  +                                  6                                                 ,                                       y                                           n                                  +                                  6                                                 )                                  \vartheta \sum_{i=0}^{7} a_i y^{n+i} + (1-\vartheta) \sum_{i=0}^{7} \tilde{a}_i y^{n+i} = \vartheta \tau f(t_{n+7}, y^{n+7}) + (1-\vartheta) \tau f(t_{n+6}, y^{n+6})                     ϑi=0∑7​ai​yn+i+(1−ϑ)i=0∑7​a~i​yn+i=ϑτf(tn+7​,yn+7)+(1−ϑ)τf(tn+6​,yn+6)
其中,( \vartheta ) 是加权参数,( \tilde{a}_i ) 是移位BDF方法的系数。
6. 应用

BDF方法广泛应用于科学计算和工程范畴,特别是在求解刚性微分方程和偏微分方程时。六步BDF方法由于其精良的稳定性,常用于实际计算中。而七步BDF方法虽然精度高,但由于稳定性题目,通常需要与其他方法结合使用,如WSBDF7方法。
7. 总结



  • 六步BDF方法:稳定性好,适用于刚性方程,广泛应用于科学计算。
  • 七步BDF方法:精度高,但稳定性差,通常需要与其他方法结合使用。
  • WSBDF7方法:通过加权和移位技术,构造出一个稳定的七步BDF方法,适用于求解刚性微分方程。
这些方法在数值分析中具有紧张意义,特别是在处理复杂和刚性题目时,选择合适的BDF方法可以显著提高计算效率和精度。
逐句页分析

好的!以下是对每一页内容的更加详细的解读,涵盖了每一页的焦点思想、公式推导、技术细节以及背景知识。

第1页



  • 标题与作者信息

    • 标题为“抛物型方程的加权和移位七步BDF方法”,表明报告的焦点内容是针对抛物型方程(如热方程)的数值解法。
    • 作者是陈明华(兰州大学),合作者包括Georgios Akrivis和Fan Yu。
    • 集会信息:2024年微分方程数值方法及其应用青年学术研讨会,湘潭大学,2024年4月21日。


第2页



  • 报告大纲

    • 引言:介绍WSBDF7方法及其在偏微分方程(PDE)中的应用,讨论稳定性地区。
    • 乘子:介绍多步方法的乘子技术,讨论如何确定乘子。
    • 稳定性与收敛性分析:通过能量技术分析方法的稳定性和收敛性。
    • 数值结果:展示WSBDF7方法在实际题目中的数值表现。


第3页



  • 动机

    • 背景:六步BDF方法通过能量技术被证明是稳定的,但七步BDF方法对抛物型方程不稳定(以致不满意零稳定性)。
    • 目的:构建一种稳定的线性组合,即七步BDF方法及其移位版本(WSBDF7方法),并通过能量技术证明其稳定性。
    • 创新点:确定了新的乘子,使得WSBDF7方法在抛物型方程中稳定。


第4页



  • WSBDF7方法用于偏微分方程

    • 题目形貌:考虑初值题目 ( u’(t) + Au(t) = f(t) ),其中 ( A ) 是椭圆算子(如拉普拉斯算子),( f(t) ) 是外力项。
    • 插值多项式:通过拉格朗日插值多项式 ( P_7 ) 在节点 ( t_n, t_{n+1}, \ldots, t_{n+7} ) 上逼近函数 ( y ) 的导数。
    • 七步BDF方法:通过逼近 ( y’(t_{n+7}) ) 构造七步BDF方法,情势为:
      [
      \sum_{i=0}^{7} a_i y^{n+i} = \tau f(t_{n+7}, y^{n+7}).
      ]


第5页



  • 移位七步BDF方法

    • 构造方法:通过逼近 ( y’(t_{n+6}) ) 构造移位七步BDF方法,情势为:
      [
      \sum_{i=0}^{7} \tilde{a}i y^{n+i} = \tau f(t{n+6}, y^{n+6}).
      ]
    • 加权组合:将七步BDF方法和移位七步BDF方法加权组合,得到WSBDF7方法:
      [
      \vartheta \sum_{i=0}^{7} a_i y^{n+i} + (1-\vartheta) \sum_{i=0}^{7} \tilde{a}i y^{n+i} = \vartheta \tau f(t{n+7}, y^{n+7}) + (1-\vartheta) \tau f(t_{n+6}, y^{n+6}).
      ]


第6页



  • WSBDF7方法的特性多项式

    • 定义特性多项式 ( \alpha ) 和 ( \beta ):
      [
      \alpha(\zeta) := \vartheta a(\zeta) + (1-\vartheta) \tilde{a}(\zeta), \quad \beta(\zeta) := \vartheta \zeta^q + (1-\vartheta) \zeta^{q-1}.
      ]
    • 其中 ( a(\zeta) ) 和 ( \tilde{a}(\zeta) ) 分别是七步BDF方法和移位七步BDF方法的特性多项式。


第7页



  • 稳定性地区

    • BDF方法的稳定性:BDF( q ) 方法是 ( A(\varphi_q) )-稳定的,给出了不同 ( q ) 值对应的 ( \varphi_q )。
    • WSBDF方法的稳定性:WSBDF( q ) 方法也是 ( A(\tilde{\varphi}_q) )-稳定的,给出了不同 ( q ) 值对应的 ( \tilde{\varphi}_q )。
    • 数值计算:通过数值计算得到了 ( \tilde{\varphi}_q ) 的近似值。


第8页



  • Nørsett的稳定性准则

    • 给出了BDF( q ) 方法的 ( A(\varphi_q) )-稳定性准则:
      [
      \tan(\varphi_q) = \min_{x \in D_q} \left( -\sqrt{1-x^2} \cdot \frac{I_q(x)}{R_q(x)} \right).
      ]
    • 讨论了WSBDF( q ) 方法的稳定性角 ( \hat{\varphi}_q ),并证明了 ( \hat{\varphi}q \to \varphi{q-1} ) 当 ( \vartheta \to \infty )。


第9页



  • WSBDF7方法的稳定性地区图

    • 绘制了 ( \vartheta = 1, 3, 10 ) 时WSBDF7方法的稳定性地区图,展示了随着 ( \vartheta ) 增加,稳定性地区逐渐增大。


第10页



  • 多步方法的乘子技术

    • 引用了Dahlquist的引理,讨论了多步方法的A-稳定性和G-稳定性。
    • 引理1指出,假如多项式 ( \rho ) 和 ( \sigma ) 满意一定条件,则存在一个正定对称矩阵 ( G ),使得方法稳定。


第11页



  • 乘子的定义

    • 定义了乘子 ( \mu(\zeta) := \zeta^q - \mu_1 \zeta^{q-1} - \cdots - \mu_q ),并给出了乘子需要满意的条件:

      • 多项式 ( \alpha ) 和 ( \mu ) 没有公因子。
      • 满意A-稳定性条件:
        [
        \text{Re} \frac{\alpha(\zeta)}{\mu(\zeta)} > 0 \quad \text{for} \ |\zeta| > 1.
        ]



第12页



  • 条件(7)和(8)的等价情势

    • 给出了条件(7)和(8)的等价情势,并通过复变函数的方法举行了扼要证明。


第13页



  • WSBDF7方法的乘子

    • 给出了WSBDF7方法的乘子 ( \mu_1 = 1.6, \mu_2 = -1.6, \mu_3 = 1.1, \mu_4 = -0.3, \mu_5 = \mu_6 = \mu_7 = 0 )。
    • 讨论了这些乘子如何满意稳定性条件。


第14页



  • 多项式 ( \alpha ) 和 ( \mu ) 的条件

    • 证明了多项式 ( \alpha ) 和 ( \mu ) 满意条件(7),即它们没有公因子且满意A-稳定性条件。


第15页



  • 多项式 ( \delta ) 和 ( \kappa ) 的条件

    • 证明了多项式 ( \delta ) 和 ( \kappa ) 满意条件(8),即它们没有公因子且满意更严酷的稳定性条件。


第16页



  • 乘子简直定

    • 讨论了确定WSBDF7方法乘子的必要条件,给出了乘子需要满意的不等式。


第17页



  • 乘子的范围

    • 给出了乘子 ( \mu_1, \mu_2, \mu_3, \mu_4 ) 的取值范围:
      [
      1.5561 \leq \mu_1 < 2.3133, \quad -2.2024 < \mu_2 < -1.4259, \quad 0.5394 < \mu_3 < 1.3955, \quad -0.6518 < \mu_4 < -0.0504.
      ]


第18页



  • 稳定性分析

    • 给出了WSBDF7方法的稳定性定理,证明了方法在能量意义下是稳定的。


第19页



  • 稳定性估计

    • 讨论了逼迫项的估计,并给出了稳定性估计的详细推导。


第20页



  • 第二稳定性估计

    • 给出了WSBDF7方法的第二稳定性估计,并举行了扼要证明。


第21页



  • 稳定性证明的细节

    • 详细讨论了稳定性证明的细节,并引入了新的符号。


第22页



  • 稳定性证明的继承

    • 继承讨论了稳定性证明的细节,并给出了干系的不等式。


第23页



  • 第三稳定性估计

    • 给出了WSBDF7方法的第三稳定性估计,并举行了扼要证明。


第24页



  • 误差估计

    • 假设解 ( u ) 足够光滑,并给出了误差估计的命题。


第25页



  • 数值结果

    • 将WSBDF7方法应用于初始和边界值题目,给出了具体的方程和初始条件。


第26页



  • 周期性边界条件的数值结果

    • 给出了周期性边界条件下的离散 ( L^2 )-范数误差和数值收敛阶数。


第27页



  • 应用

    • 讨论了所提出的方法在平均曲率流、梯度流和分数阶方程等范畴的应用。


第28页



  • 致谢

    • 感谢听众的到场。


总结

这份PPT详细介绍了WSBDF7方法的构建、稳定性分析、收敛性分析以及数值结果,展示了该方法在抛物型方程中的应用远景。通过加权和移位技术,成功解决了七步BDF方法的不稳定性题目,并通过能量技术证明了其稳定性和收敛性。

免责声明:如果侵犯了您的权益,请联系站长,我们会及时删除侵权内容,谢谢合作!更多信息从访问主页:qidao123.com:ToB企服之家,中国第一个企服评测及商务社交产业平台。

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?立即注册

x
回复

使用道具 举报

0 个回复

倒序浏览

快速回复

您需要登录后才可以回帖 登录 or 立即注册

本版积分规则

刘俊凯

论坛元老
这个人很懒什么都没写!
快速回复 返回顶部 返回列表