第三章：渐进方法.ppt

数学物理方法概论 之 渐进方法 主讲教师 白璐联系电话mail blu 第三章渐近方法 本章渐进方法着重介绍数学物理中的近似方法 内容包括积分的渐近展开分析与常微分方程的渐进解法两大部分 通过本章的学习目的是为提高数学分析的能力和将理论应用于解决实际问题的本领 该方法在力学 大气科学 物理海洋 光学 声学等研究领域具有广泛的应用 渐近计算是数学计算的近似方法之一 它是解析方法在一定条件下的发展 其与数值方法相结合可以提高计算的精确程度及计算速度 特别在非线性问题的处理中渐近方法具有重要的地位 1 量级符号 2 渐近展开 3 渐近展开式的运算 4 积分的渐近展开式 5 最陡下降法 6 驻定相位法 7 常微分方程的渐近解 第三章渐近方法 由于某些特殊函数具有积分表示式 如果这些函数是微分方程的解 就可以得到一种以它们的拉普拉斯变换或傅立叶变换的积分表达式表达的解 因此求解积分的渐近展开式的问题在解析函数理论中就起特别重要的作用 它可以使我们得到积分解另一种表达 称此为渐近方法 比较函数趋于某个极限时的性质常定义 例 3渐近方法 3 1量级符号 1 同量级 例 称函数f x 至多与g x 同阶 3渐近方法 3 1量级符号 2 量级最多为 也可以说若存在某个常数A 使对定义域D某个内点x0的邻域V内的所有x 满足 例 的意义是说f x 有界 而的意义是说f x 趋于零 3渐近方法 3 1量级符号 3 量级小于 也可以说若存在任一 定义域D内点x0总有一的邻域存在 使得所有 满足 称函数f x 是函数g x 的高阶小量 3 2渐近展开 下面给出渐近展开的定义和它的一些性质 讨论在扩充的复平面上进行 一 渐近序列 设 是定义在区间D上的连续函数序列 是D中的一固定点 若对每一个固定的n 有 则称为点的渐进序列 渐近序列可以是有限项也可以是无限项的 例如 是对零点的渐近序列 3渐近方法 是对于无穷的渐近序列 二 渐近展式 设是一个给定的函数 而是点的一个渐近序列 如果对每个固定的整数n 有那么称此为在点的渐近展式 记为注意 渐近展式与函数的级数展式不同 对确定的z值 渐近展式的项数无限增多时 所得级数一般是发散的 但若满足渐近展式的定义式 则当时 取确定的项数n会得到对函数非常好的近似 3 2渐近展开 3渐近方法 例1 求当时的积分值 即求时的渐近展式 解 余项 3 2渐近展开 3渐近方法 因此 取展开式的前n项 略去余项 当时 其误差量级小于所取的最后一项 符合渐近展式的定义 可记为 3 2渐近展开 3渐近方法 注意 这个级数对于有限的x值均不收敛 但是 取确定的项数 会得到对函数很好的近似 如果仅用一项 给出的相对误差为1 x 结果粗略一些 但已经足够用了 三 展开式系数 当时 的渐近展式的系数为 证明略 3 2渐近展开 3渐近方法 四 展开式的构成 设在区域D中有定义 若有定义且不为零 则是时 的一个直到N项的渐近展开式 证明 首先证明是一个渐近序列 由的定义得 3 2渐近展开 3渐近方法 所以 又因为 故存在一个的邻域使z在其中时 所以 由此 各个都由这种方式定义得 3 2渐近展开 3渐近方法 五 唯一性 设是在D中 的一个已知渐近序列 若是当时 直到N项的一个渐近展式 则此展式是唯一的 注意 这个定理只表示用同一个已知渐近序列表示的展开式的唯一性 但是可能有多个不同的渐近序列对应同一个函数的渐近展式 它们可以不同 而且可以是收敛的也可是发散的 反过来 一个已知的渐近展式可以表示不止一种函数 的一个渐近幂级数展式 记为 六 幂函数的展式 则 是D中 时 3 2渐近展开 3渐近方法 其中一种重要的特殊情形是在D中 当时 如果 则在D中 当时 3 3渐近展式的运算 若在D中 当时 直到N项有则 和 1 加法 2 乘法 3渐近方法 本节讨论渐近展开式的普通运算 由于实际应用中 展式多用幂函数 以下均以幂函数作为渐近序列 3 除法 即除法为两个函数渐近展开式分别保留到N项相除 推论 3 3渐近展式的运算 3渐近方法 4 积分 当时 若则 其中积分沿从到的一条直线路径 推论 当时 若则 3 3渐近展式的运算 3渐近方法 5 求导 当时 若 且当时 在D中存在并有 则在D中渐近展开式满足可逐项积分的条件时 有 推论1 在D中 当有 且在D中 3 3渐近展式的运算 3渐近方法 存在并有 若在D中 渐近幂级数满足逐项积分的条件 则 推论2 对 当时有且存在于相同的区域 当时 有则 对于解析函数 若在区域当时有则在中 当有 3 3渐近展式的运算 3渐近方法 根据渐近展式的定义和相关运算法则 就可以讨论在解析函数理论中常用的积分的渐近展式 获得积分渐近展式的方法有两种把被积函数的一部分展开为级数 然后形式上逐项积分 重复地进行分布积分 一 逐项积分法 瓦特森引理 设 3 4积分的渐近展式 3渐近方法 式对Re z 0成立 因为在此定义域两边都解析且在实轴上它们一致 可应用瓦特森引理得到其积分的渐近展开式 做变量代换 令 解 令则 例 求当 的 函数的渐近展式 3渐近方法 则对给定的值上述变换给出两个解s u 和 u 其中 3 4积分的渐近展式 即 且 两个解分别位于最大值s 1的两边其中 于是 3渐近方法 3 4积分的渐近展式 可以证明 且因当时 故在有界 3渐近方法 3 4积分的渐近展式 则可得与 的关系 剩下要证明的是其中对小的有一个 麦克劳林展开式 再做代换 令 它在处是解析的 因为当 时 有 即 与的邻域有两个分支 根据复变函数理论 若 解析 且 则在 的邻域存在解析的反函数 现在 在 邻域解析 且 在 点不等于零 故在 3渐近方法 3 4积分的渐近展式 另一支是 注意到 则对足够小的有 故 令 的邻域存在解析的反函数 式中 是处的留数 容易算出等 将最后的表示式带入被积式 并在形式上逐项积分 则由瓦特森引理 在时 有 3渐近方法 3 4积分的渐近展式 对 的条件下得 式中 二 分部积分法 形式进行分部积分 在 且当时 可以看出 所得积分是前面考虑过的形式 故可重复同一过程 3渐近方法 3 4积分的渐近展式 存在 且 当时 首先证明在和 的一定假设条件下 式中第一项是 的积分渐近形式 设 1 对 连续且有界 同时 2 对 为实函数且连续 3 4 对所有正 且当时 5 对 存在 则对 3渐近方法 3 4积分的渐近展式 例 在 条件下 求误差函数 的渐近展开式 令 这表明 对 的两边是解析的 且在实轴上相合 现在的积分 因为在此区域中方程 和定理 的假设相符 重复地应用此定理 对 可得 3渐近方法 3 4积分的渐近展式 如果 则应将方法修改 但对这种情形 可以采用 下面两节的方法 这里不再赘述 以上只把分部积分法用于上限为 的积分 现在考虑a和b有限 且 的情形 即 设 而当 时 3渐近方法 3 4积分的渐近展式 当 时 因为 故 采用最徒下降法的思路 首先令 则 3渐近方法 3 5最陡下降法 常遇到这样类型的积分 其中C 是复平面Z 上的路径 在其中假定 缓变 且f和g 均具有适当的正则性 其中u和v是实函数 当S很大时 沿积分路径微小位移所引起的v 的微小变化会引起 注意到 3 5最陡下降法 也就是说 最徒下降法的本质就是尽可能利用这样的积分路径 使被积函数在这个路径上u为最大 而v等于常数 这样可以保证被积函数变化最速下降 也就保证积分值只与u为最大的点 鞍点 附近的邻域有关 从而可以渐近计算 事实上 使v等于常数的路径也就是u变化最快的路径 以下对此证明 3渐近方法 中正弦项的迅速震荡 但如选择积分路径使在其上v为常数 则震荡就会迅速消失 于是被积函数变化最速部分将为 而 显然其主要贡献部分将来自u为最大点的邻域 因而此方法的本 质是尽可能地改变积分路径循着通过u 为最大的点 而v等于常数 的路线进行 3 5最陡下降法 3渐近方法 证明 令 是在 邻近的一点 于是由 得 当 等于常数时 应有 即 这样 注意到柯西 黎曼关系 可得 此式表明 因此 由极值的条件 在 点 v等于常数的方向也正是u 的最大 变化方向 3 5最陡下降法 3渐近方法 为寻找 的最大点 令 因而 故当且仅当在该点 时 取得极值 这样的点称为驻点 曲面 有极大极小值的条件为 现在有 即 故 因而 因为u是解析函数满足拉普拉斯方程 表明驻点不是极值点而是鞍点 它连接曲面的 山谷 和 山脊 沿山脊上升和山谷下降均是u最大变化方向 对我们有意义的是山谷下降路径 即最徒下降路径 因为只有这一路径上在鞍点附近对积分有显著的贡献 所以这种渐近计算的方法称为 最徒下降法 3 5最陡下降法 3渐近方法 鞍点 若 点为鞍点 即 此点的v等于常数的曲线方程为 则 通过 或 其中t是实数 t为正代表下降路径 t为负代表上升路径 由 在 点的Tailor展开式 现在 若 A为正实数 接近 处 则 略去高阶项 3 5最陡下降法 3渐近方法 因为 还可得 由此可以画出 实部 虚部 时的等高线如图所示 如果 则图形将更复杂 可能有三个或更多的山谷在 鞍点相会 3 5最陡下降法 3渐近方法 现在可以假定起止于无限远的积分路径能变形到起点和终点都在山谷的路径 这是积分收敛的要求 积分路径要尽可能地变形到最陡下降路径上沿着山谷的底在鞍点处越过一个山谷进入下一个山谷 一般说来 这种路径由一系列曲线组成 每一个是从鞍点到无穷或到某个奇点 以下假定 来计算一个这种路径对积分的贡献 为此 设 其中 于是最陡下降路径由下式给出 或 t为正实数 其中 取主值 计及 故得 3 5最陡下降法 3渐近方法 上式的不同符号对应于自鞍点出发的两条最陡下降路径 若 号与第三象限的路径有关 与第一象限的路径 有关 考虑负号时所代表的路径如图所示 所得的积分是 负号所对应的路径 其中 是上式中取负号的 z值 另一路径的积分 其中 是上式中取正号的z值 完整的级数太繁 我们将只导出首项 于是 如果把C变形到通过鞍点 其方向如右图 则可以得到 3 5最陡下降法 3渐近方法 由于和 都可用t表示 其中函数f 已假定是缓变的 故 和 均可用 代替 令 引理渐近计算的积分式 则可以得到用 瓦特森 负号所对应的路径 此式即利用最徒下降法得到的积分的渐近展开 如果C通过鞍点的方向与前图示相反的话 结果反号即可 例 求阶乘的斯特林 Stirling 公式 即阶乘的渐近展式 解 已知阶乘的积分表达式 符合前边积分的形式 其中 而积分路径C为实轴 3 5最陡下降法 3渐近方法 它在 时成立 现在我们只考虑s 此积分形式不适合用最陡下降法 但如用sz来代替z就得出 是实数的情形 注意到 有一鞍点 且在该处 在时 因此积分路径应该是和 零点为奇点 两部分 根据公式 3 5最陡下降法 3渐近方法 这就是斯特林公式 3 6驻定相位法 对下列形式的积分 3渐近方法 当参量k 很大时可以用驻定相位法求解 从被积函数的形 式 上看 可当作波的相位 当k很大时 它表示一种迅速的振荡 在积分过程中 这种振荡正负相消 而只有在 的部分 处有平坦 因而对积分的主要贡献来自于 点附近 使 的点称为驻定相位点 所以这种用相位驻定邻近的积分 结果来近似代表整个区间的精确结果的方法称为驻定相位法 为证明对积分的主要贡献来自驻定相位点附近 先看变量z 为实变量x的情形 函数 的驻点 是使 的点 如 而 则称 为 的N级驻点 考察积分 3 6驻定相位法 3渐近方法 如果积分区间 a b 内f x 没有驻点 g x 在 a b 内可微 则可把f作为积分变量 而上面的积分可记为 由f x 反演x可以表示为f的函数 故在积分区间是可微的 经过分部积分 则 3 6驻定相位法 3渐近方法 等号右边第一项在 时趋于零 其量级为 右边第 二项形式上与原积分一样 可微能对它再进行分部积分 积分后的量级也是 但其前面已有系数 故上式等 号右边第二项的量级为 再继续进行分部积分 可见整个 在k很大时量级最多为1 k的小量 如果在积分区间内 有一个一级驻点 则由于 而使 在处失去可微性 因而不能直接进行分部积分 则 后一个积分中把f作为积分变量 则 其中 令 3 6驻定相位法 3渐近方法 它在驻点处 为型 的极限为 从而 也在积分区间内可微 由此 上面后一个积分的量级也是 现在来考虑前一个积分 将 按Taylor级数展开 则 3 6驻定相位法 3渐近方法 略去后面的高阶项 计及g x 在区间内是x的缓变函数 于是上述 积分整个地可写为 令 则 再令 同时考虑到 时积分限 则得 3 6驻定相位法 3渐近方法 由于 得 当 时 第一项的量级为 积分的贡献 所以 与不包含驻点的区间相比较 当时 也就是包含驻点的区间对 含驻点的区间对积分的贡献要重要的多 计及 可能为正或者负 通常把上述结果表示为 如果区间内有多个一级驻点 可分为若干个子区间 使各个驻点都在其中 然后逐个用上式计算 再将结果加起来得到所需结果 对于复变量的情形 积分可写成 于是 即鞍点 3 6驻定相位法 3渐近方法 其中 是给定的积分路径 由 其中C 立即得出在点 从前面的讨论可知对积分的主要贡献来自相位 稳定的区域 即应来自 的极值点附近 所以希望在经过 选出一个特定的方向 沿此方向 相位能最迅速地变化而在 点的所有方向中 点取极值 3 6驻定相位法 3渐近方法 由最陡下降法的分析可知 如取 为常数的路径必为 变化最激烈的路径 可见 在复变函数情形下的驻相法实际 上是上节最陡下降法取共轭路