2 目 录 Euler 方 法 龙 格 库 塔 方 法 两 点 边 值 问 题

Size: px
Start display at page:

Download "2 目 录 2.4.1 Euler 方 法................................ 39 2.4.2 龙 格 库 塔 方 法............................ 40 2.4.3 两 点 边 值 问 题..........................."

Transcription

1 目 录 第 零 章 绪 论 5 第 一 章 计 算 机 数 和 误 差 计 算 机 数 及 其 表 示 舍 入 误 差 对 计 算 的 影 响 减 法 的 计 算 第 二 章 物 理 学 中 的 常 用 数 值 方 法 数 值 积 分 方 法 梯 形 法 和 辛 普 生 方 法 反 常 积 分 的 计 算 高 斯 积 分 方 法 哥 西 主 值 的 计 算 急 速 振 荡 函 数 的 积 分 高 维 积 分 的 计 算 固 体 热 容 量 的 计 算 方 程 的 求 根 和 最 优 化 方 法 二 分 法 牛 顿 法 和 割 线 法 一 维 最 优 化, 黄 金 分 割 法 高 维 最 优 化, 共 扼 梯 度 方 法 约 束 最 优 化 最 小 二 乘 法 及 曲 线 拟 合 函 数 的 求 值 多 项 式 的 求 值 和 级 数 求 和 连 分 式 的 求 值 常 微 分 方 程 的 数 值 解 法

2 2 目 录 Euler 方 法 龙 格 库 塔 方 法 两 点 边 值 问 题 插 值 和 逼 近 多 项 式 插 值 分 式 有 理 函 数 插 值 和 外 推 三 次 样 条 插 值 列 表 函 数 的 积 分 Padé 插 值 与 外 推 发 散 级 数 的 Cesáro 求 和 及 其 推 广 Borel 求 和 Bézier 逼 近 多 元 函 数 的 插 值 及 逼 近 快 速 付 里 叶 变 换 付 里 叶 变 换 离 散 付 里 叶 变 换 快 速 付 里 叶 变 换 附 录 : 特 殊 函 数 的 计 算 误 差 函 数 的 计 算 指 数 积 分 的 计 算 整 数 阶 贝 塞 尔 函 数 及 虚 宗 量 贝 塞 尔 函 数 的 计 算 球 谐 函 数 的 计 算 椭 圆 积 分 的 计 算 附 录 : 高 斯 积 分 的 节 点 和 权 重 第 三 章 数 值 线 性 代 数 主 元 消 去 法 解 线 性 代 数 方 程 组 LU 分 解 法 三 对 角 方 程 组 的 求 解 实 对 称 矩 阵 的 本 征 值 和 本 征 向 量 Householder 方 法 QL 算 法 第 四 章 电 磁 场 的 计 算 Maxwell 方 程, 边 值 问 题 差 分 和 差 商

3 目 录 二 维 Possion 方 程 的 五 点 差 分 格 式 差 分 方 程 的 求 解 简 单 迭 代 法 Jaobi 方 法 逐 次 迭 代 法 Gauss Seidal 迭 代 方 法 逐 次 超 松 弛 迭 代 法 (Suessive Over-Relaxation) 变 分 方 法 复 习 有 限 元 方 法 的 理 论 基 础 用 有 限 元 方 法 求 解 二 维 Laplae 方 程 第 五 章 Monte Carlo 方 法 随 机 变 量 及 其 分 布 赝 随 机 数 的 产 生 用 Monte Carlo 方 法 计 算 积 分 自 相 似 性 与 分 形 扩 散 限 制 粘 结 的 计 算 机 模 拟 高 分 子 的 模 拟 和 无 规 行 走 第 六 章 逾 渗 和 统 计 物 理 问 题 逾 渗 简 介 逾 渗 的 计 算 机 模 拟 和 分 析 正 则 系 综 和 统 计 模 型 Ising 模 型 Lenard-Jones 模 型 统 计 物 理 的 Monte Carlo 模 拟, Metroplis 方 法 Ising 模 型 的 Monte Carlo 模 拟 经 典 统 计 物 理 复 习 The miro anonial ensemble(nve) The anonial ensemble(nvt) The NPT ensemble(npt) The grand anonial ensemble( Thermodynamial relations Correlation funtions Time orrelation funtion and transport oeffiients 液 体 模 型 的 Monte Carlo 模 拟 Monte Carlo Simulation of The Lenard -Jones Model 自 由 能 计 算

4 4 目 录 王 -Landau 方 法 及 其 它 分 子 动 力 学 方 法 General proedure of MD (NVE ensemble) Simulation of Lennard-Jones liquids Simulation of Hard-sphere systems Simulation of Langevin dynamis Long range fore and Ewald summation 第 七 章 原 子 结 构 的 计 算 原 子 结 构 问 题 变 分 法 Viral 定 理 和 Hellmann-Feynman 定 理 轨 道 近 似 和 Hartree Fok 方 程 统 计 近 似 和 X α 方 程 径 向 方 程 的 求 解 方 法 计 算 程 序

5 第 零 章 绪 论 不 论 是 理 论 物 理, 还 是 实 验 物 理 都 离 不 开 数 值 计 算, 例 如, 量 子 力 学 中 简 谐 振 子 的 波 函 数 可 以 用 厄 米 多 项 式 来 表 示, 但 为 了 得 到 波 函 数 的 数 值, 仍 然 需 要 作 数 值 计 算 ; 实 验 上 测 量 到 的 数 据 要 进 行 数 据 处 理, 也 涉 及 到 数 据 拟 合, 分 析 并 计 算 误 差 等 一 系 列 数 值 计 算. 大 家 在 过 去 的 学 习 中 已 经 知 道, 牛 顿 力 学 方 程 只 有 二 体 问 题 是 可 解 的 ( 自 由 粒 子 或 小 振 动 等 特 殊 问 题 除 外 ), 三 体 以 上 的 问 题 曾 折 磨 了 全 世 界 的 许 多 优 秀 的 数 学 家 和 理 论 物 理 学 家, 但 最 终 也 未 得 到 解 析 解 ; 麦 克 斯 韦 方 程 组 只 有 在 一 些 非 常 特 殊 的 条 件 下 才 能 求 得 解 析 解 ; 量 子 力 学 的 薛 定 格 方 程, 除 了 氢 原 子 和 简 谐 振 子 外, 还 没 有 一 个 真 实 的 物 理 问 题 可 以 找 到 解 析 解 ; 统 计 物 理 告 诉 我 们, 只 要 求 得 了 系 统 的 配 分 函 数, 则 一 切 物 理 量 都 可 通 过 求 微 分 得 到, 然 而, 除 了 可 数 的 几 个 二 维 模 型 外, 没 有 一 个 实 际 的 物 理 系 统 的 配 分 函 数 被 求 出 来 过. 一 些 大 规 模 的 物 理 实 验, 耗 资 在 数 亿 美 元, 如 果 没 有 充 分 的 论 证 就 去 做 的 话, 其 浪 费 是 十 分 巨 大 的. 对 于 这 样 一 类 问 题, 大 规 模 的 数 值 计 算 往 往 可 以 发 挥 重 要 作 用. 大 规 模 数 值 计 算 也 提 供 了 理 论 物 理 和 实 验 物 理 之 间 的 桥 梁. 几 乎 所 有 的 实 际 物 理 问 题 都 无 法 得 到 精 确 解, 发 展 各 种 理 论 模 型 和 近 似 方 法 就 成 为 理 论 物 理 学 的 重 要 方 面. 理 论 模 型 的 正 确 性 最 终 需 要 经 受 实 践 的 检 验, 但 对 于 一 个 确 定 的 物 理 模 型, 其 对 应 的 实 验 体 系 往 往 很 难 获 得, 为 了 检 验 针 对 一 个 确 定 模 型 的 近 似 方 法, 最 好 能 够 得 到 这 个 模 型 的 精 确 解, 而 数 值 模 拟 方 法 正 好 能 够 提 供 这 样 的 结 果. 大 量 的 实 验 体 系 往 往 包 含 多 种 因 素, 这 些 因 素 有 些 可 以 排 除, 而 有 些 非 常 难 以 排 除. 如 果 一 个 理 论 的 出 发 点 包 含 了 所 有 的 实 验 因 素, 则 这 样 的 理 论 几 乎 肯 定 得 不 到 解 析 结 果. 但 是, 数 值 模 拟 方 法 一 般 总 能 够 把 各 种 因 素 包 含 进 去, 从 而 能 够 更 好 的 模 拟 实 验 环 境. 需 要 指 出 的 是, 长 期 以 来, 我 们 总 有 一 些 误 解, 一 讲 到 计 算, 似 乎 就 需 要 超 级 计 算 机. 事 实 上, 现 在 使 用 的 笔 记 本 计 算 机, 各 种 微 型 计 算 机 的 计 算 能 力 往 往 比 10 年 前 的 超 级 计 算 机 的 计 算 能 力 更 强. 用 一 台 目 前 似 乎 已 经 被 人 看 不 起 的 IBM-PC/XT 微 机 (IBM-PC 是 基 于 intel 芯 片 的 第 一 代 个 人 计 算 机, 大 约 1980 年 推 出, 其 CPU 是 intel-8086 芯 片, 有 64k 内 存, 两 个 5 吋 软 驱, 机 器 带 有 MS-DOS 1.0 和 CP/M 两 个 操 作 系 统,IBM-PC/XT 是 IBM-PC 的 扩 展, 增 加 了 一 个 10M 的 硬 盘, 并 且 把 内 存 增 加 到 256K), 就 可 以 很 容 易 地 解 决 三 体 问 题, 四 体 问 题,. 所 有 元 素 周 期 表 上 的 原 子 的 能 级 和 波 函 数 也 可 以 相 当 精 确

6 6 第 零 章 绪 论 地 计 算 出 来. 如 果 你 拥 有 一 台 目 前 市 场 上 一 般 配 置 的 微 机, 你 就 可 以 模 拟 一 些 二 维 和 三 维 模 型 系 统 的 配 分 函 数, 如 Lenard-Jones 流 体, 三 维 Ising 模 型 等, 可 以 计 算 元 素 固 体 和 简 单 化 合 物 晶 体 的 电 子 结 构, 声 子 谱 等. 进 一 步, 你 如 果 拥 有 或 能 够 使 用 到 大 型 计 算 机, 或 利 用 微 机 构 成 计 算 集 群, 你 就 可 以 作 真 实 流 体 的 流 体 力 学 计 算, 可 以 计 算 复 杂 边 界 条 件 下 的 三 维 电 磁 场, 可 以 计 算 三 维 真 实 系 统 的 配 分 函 数,. 同 样 你 也 能 够 在 大 型 物 理 实 验 开 始 前 作 很 多 模 拟, 对 实 验 的 各 种 情 况 有 一 个 清 楚 的 概 念. 本 课 程 的 目 的 在 于 对 物 理 中 的 数 值 计 算 进 行 一 些 入 门 介 绍, 使 大 家 在 学 完 本 课 程 后, 在 组 织 一 些 较 大 规 模 的 计 算 时 心 中 有 数, 少 走 弯 路. 课 程 分 两 部 分, 一 部 分 是 基 本 算 法, 主 要 介 绍 计 算 机 数 的 特 点, 舍 入 误 差, 计 算 稳 定 性 等 基 本 概 念 和 一 些 常 用 的 数 值 计 算 方 法 ; 第 二 部 分 则 结 合 物 理 问 题, 讲 一, 二 个 专 题, 使 大 家 学 到 一 些 组 织 较 大 规 模 计 算 的 步 骤 和 技 巧. 学 习 本 课 程 要 求 学 过 四 大 力 学 ( 经 典 力 学, 电 动 力 学, 热 力 学 和 统 计 物 理 学 及 量 子 力 学 ) 和 数 学 物 理 方 法, 掌 握 至 少 一 门 计 算 机 高 级 语 言 ( 如 C, PASCAL, FORTRAN, Java 等 ). 计 算 是 一 门 实 践 性 很 强 的 课 程, 上 机 实 习 是 本 课 程 的 一 个 有 机 组 成 部 分. 关 于 本 讲 义 的 一 点 说 明 : 这 是 笔 者 在 1993 年 为 开 设 计 算 物 理 课 程 准 备 的 讲 义 年, 上 海 交 通 大 学 应 用 物 理 系 代 理 系 主 任 张 仲 渊 教 授 找 到 我, 告 诉 我 决 定 在 本 科 生 高 年 级 开 设 计 算 物 理 课 程, 并 安 排 由 笔 者 开 课. 在 接 到 这 一 任 务 后, 笔 者 翻 阅 了 当 时 国 内 出 版 的 几 本 计 算 物 理 教 材, 发 现 都 过 于 专 门, 不 适 于 做 为 本 科 高 年 级 的 教 材, 而 更 适 合 于 某 一 专 业 的 研 究 生 教 材. 为 此, 我 只 好 试 图 自 己 写 一 本 讲 义, 以 应 急 需. 基 本 算 法 部 分 主 要 参 考 了 W. H. Press, S. A. Teukolsky, W. T. Vetterling and B. P. Flannery 所 著 Numerial Reipes, The Art of Sientifi Computing 一 书, 当 时, 这 本 书 出 版 时 间 还 不 算 太 长, 但 已 经 成 为 大 多 数 经 常 从 事 数 值 计 算 的 物 理 学 工 作 者 的 手 头 必 备 书, 我 自 己 当 时 也 买 了 一 本, 书 价 是 我 的 两 个 月 的 工 资 这 本 书 一 方 面 讲 解 算 法, 另 一 方 面 对 每 一 个 算 法 给 出 了 透 明, 简 洁 的 源 程 序 尽 管 这 些 程 序 的 效 率, 安 全 性 方 面 比 一 些 成 熟 的 软 件 包 要 差, 但 由 于 其 简 洁 清 楚, 特 别 适 合 于 修 改 后 放 到 自 己 的 程 序 中 去 讲 义 的 内 容 曾 在 上 海 交 通 大 学 应 用 物 理 系 讲 授 多 次, 各 个 年 级 的 同 学 积 极 参 与 了 本 课 程 的 学 习 并 对 课 程 的 讲 法 提 出 了 很 多 非 常 宝 贵 的 意 见 和 建 议, 这 些 建 议 在 讲 义 的 成 文 中 对 笔 者 有 很 大 帮 助, 在 此 表 示 感 谢 年 后, 物 理 系 决 定 把 计 算 物 理 课 程 改 为 双 语 课 程, 这 本 讲 义 就 算 完 成 了 其 历 史 使 命. 在 教 学 过 程 中, 我 也 发 现 这 一 讲 义 对 一 些 问 题 的 讲 解 过 于 简 洁, 而 内 容 的 深 度 也 不 太 适 合 本 科 生 的 学 习. 这 一 段 时 间, 国 际 上 出 现 了 几 本 比 较 合 适 的 教 材, 所 以,

7 7 当 时 认 为 这 一 讲 义 不 需 要 继 续 修 订, 也 就 放 在 一 边 了. 随 后, 研 究 生 的 课 程 做 了 一 些 改 革, 增 加 了 36 学 时 的 数 值 计 算, 任 课 老 师 感 觉 这 个 讲 义 的 内 容 作 为 研 究 生 数 值 计 算 的 课 程 比 较 合 适, 遂 采 用 这 本 讲 义 的 电 子 版 教 学 年 10 月, 致 远 学 院 的 叶 曦 副 院 长 和 负 责 物 理 班 教 学 的 郑 杭 教 授 找 我, 希 望 我 为 致 远 学 院 的 同 学 开 设 一 门 物 理 中 数 值 计 算 的 课 程 注 意 到 致 远 学 院 已 经 有 至 少 两 门 相 关 课 程, 这 门 课 程 应 该 更 多 强 调 物 理 方 面 为 了 这 个 课 程, 利 用 寒 假, 我 修 订 了 讲 义, 主 要 是 对 原 来 的 一 些 表 述 做 了 修 改 和 补 充, 同 时 尽 可 能 强 调 物 理 方 面 的 内 容. 鉴 于 目 前 已 经 存 在 大 量 优 秀 的 数 学 软 件, 如 Matlab, Maple 等. 大 量 的 小 规 模 的 计 算 往 往 可 以 通 过 这 些 软 件 方 便 进 行, 因 此, 此 次 修 订 时 时 也 适 当 提 到 这 方 面 的 内 容 讲 义 中 提 到 的 一 些 程 序 将 放 在 课 程 主 页 上, 这 些 程 序 中, 有 些 来 自 前 面 提 到 的 Numerial Reipes, 有 些 是 我 曾 经 多 次 使 用 过 的, 也 有 些 是 专 为 此 讲 义 而 准 备 的, 由 于 时 间 关 系, 没 有 经 过 仔 细 调 试 和 考 验, 因 此 在 用 于 实 际 问 题 时, 一 定 要 经 过 仔 细 试 算, 以 免 造 成 不 必 要 的 损 失, 同 时, 我 也 愿 在 此 声 明, 笔 者 对 因 使 用 此 讲 义 中 的 程 序 而 导 致 的 任 何 直 接 的 和 间 接 的 损 失 将 不 负 担 任 何 责 任. 所 有 程 序 都 以 FORTRAN 77 写 出, 鉴 于 C 语 言 在 实 际 工 作 中 使 用 的 越 来 越 广 泛, 同 时 也 是 很 多 同 学 喜 欢 的 语 言, 部 分 程 序 已 经 翻 译 为 C 语 言. 由 于 程 序 是 用 FORTRAN 写 的, 翻 译 中 自 然 留 有 FORTRAN 的 痕 迹. 在 上 次 讲 义 完 成 到 这 次 修 订 之 间, 面 向 对 象 的 程 序 设 计 逐 步 走 向 主 流, 其 代 表 性 语 言 是 ++ 和 Java, 建 议 有 志 于 从 事 计 算 物 理 的 同 学 注 意 这 一 点.

8 8 第 零 章 绪 论

9 第 一 章 计 算 机 数 和 误 差 在 这 一 章 中, 我 们 将 给 出 电 子 数 字 计 算 机 中 数 的 表 示, 计 算 规 则, 舍 入 误 差 及 算 法 稳 定 性 等 概 念, 在 处 理 上, 并 不 追 求 数 学 上 的 严 密, 而 是 用 一 些 启 发 式 的 推 导 和 例 子 来 说 明 主 要 概 念, 需 要 进 一 步 钻 研 的 同 学, 可 以 参 看 有 关 计 算 数 学 的 书 籍. 1.1 计 算 机 数 及 其 表 示 我 们 在 初 中 的 数 学 上 就 已 经 知 道, 数 有 整 数, 实 数 和 复 数, 整 数 可 以 取 0, ±1, ±2,..., ±n,..., ±, 实 数 则 是 整 个 数 轴 上 的 连 续 统. 数 学 分 析 就 是 建 立 在 实 数 系 上 的. 实 数 的 运 算 满 足 加 法 和 乘 法 的 交 换 律, 接 合 律 等 一 般 的 运 算 律. 而 计 算 机 所 能 表 达 的 数 的 全 体, 则 是 一 个 离 散 的, 有 限 的 集 合. 通 常 计 算 机 整 数 的 取 值 范 围 为 N 到 N 1, N 的 数 值 随 计 算 机 不 同 而 异, 如 IBM-PC 机 的 N = 实 数 也 有 一 个 取 值 范 围, 在 IBM-PC 上, 实 数 的 绝 对 值 最 大 约 为 10 38, 最 小 约 为 10 38, 大 于 最 大 值 的 数 将 造 成 上 溢, 小 于 最 小 值 的 数 将 造 成 下 溢.( 在 很 多 计 算 机 系 统 中, 下 溢 作 为 0 来 处 理 ). 计 算 机 实 数 在 其 范 围 内 的 分 布 是 非 常 不 均 匀 的, 在 靠 近 0 的 地 方 较 密, 离 开 0 越 远 越 稀 疏. 计 算 机 实 数 一 般 不 满 足 实 数 的 运 算 律. 计 算 机 中 的 实 数 通 常 用 规 格 化 的 二 进 制 浮 点 数 来 表 示, 例 如 , 和 分 别 表 示 为 , 和 由 于 计 算 机 只 有 有 限 的 字 长, 所 以 只 能 表 示 有 限 个 实 数, 例 如, 对 于 一 个 字 长 为 五 的 计 算 机, 将 无 法 区 分 , 既 然 计 算 机 的 实 数 是 一 个 有 限 的 集 合, 那 么 初 始 数 据 和 中 间 结 果 都 有 可 能 不 在 这 个 集 合 中, 于 是 必 须 用 计 算 机 实 数 去 近 似 表 示 相 应 的 实 数, 从 而 引 进 舍 入 误 差. 如 果 按 照 四 舍 五 入 的 原 则 去 近 似 表 示, 则 由 计 算 机 表 示 一 个 实 数 时 引 进 的 舍 入 误 差 不 大 于 最 后 一 位 数 的 一 半, 但 很 多 计 算 机 采 用 了 只 舍 不 入 的 方 法, 则 由 计 算 机 数 表 示 一 个 实 数 时 引 进 的 舍 入 误 差 最 大 可 为 最 后 一 个 二 进 制 位.

10 10 第 一 章 计 算 机 数 和 误 差 1.2 舍 入 误 差 对 计 算 的 影 响 在 数 值 计 算 过 程 中, 一 般 每 一 步 都 有 舍 入 误 差, 舍 入 误 差 的 积 累 有 可 能 使 整 个 计 算 失 败. 例 如, 考 虑 下 述 问 题, 计 算 积 分 利 用 分 部 积 分 法 得 到 或 1 E n = 1 0 x n e x 1 dx, n = 1, 2,..., 20. x n e x 1 dx = x n e x n x n 1 e x 1 dx, E n = 1 n E n 1 容 易 算 出 E 1 = 1/e. 利 用 上 述 递 推 公 式, 在 IBM-PC 微 机 上 用 双 精 度 计 算 结 果 如 下 : E 1 = E 2 = E 3 = E 4 = E 5 = E 6 = E 7 = E 8 = E 9 = E 10 = E 11 = E 12 = E 13 = E 14 = E 15 = E 16 = E 17 = E 18 = E 19 = E 20 = 被 积 函 数 在 积 分 区 间 [0, 1] 的 内 部 总 是 正 的, 因 此 积 分 值 不 可 能 为 负, 更 进 一 步, 由 于 E n = 1 0 x n e x 1 dx 1 0 x n dx = 1 n + 1 所 以, 必 有 E 20 1/21. 显 然, 上 述 计 算 结 果 是 不 正 确 的, 而 问 题 就 出 在 舍 入 误 差 上 面. 我 们 对 此 作 一 个 简 单 的 分 析. 双 精 度 数 可 表 示 14 位 10 进 制 数, 所 以 E 1 的 误 差 约 为 ϵ = ( 也 可 能 大 于 此 数, 这 里 仅 作 一 数 量 级 的 估 计 ). 在 递 推 过 程 中, 这 个 误 差 将 被 传 播, 放 大, 最 后 导 致 计 算 结 果 的 错 误. 误 差 的 传 播 规 律 为 E 2 = 1 2(E 1 + ϵ) = 1 2E 1 2ϵ E 3 = 1 3(1 2E 1 2ϵ) = 1 3(1 2E 1 ) + 3!ϵ

11 1.2 舍 入 误 差 对 计 算 的 影 响 11 计 算 到 E 15 时 的 误 差 为 15!ϵ 10 3, 而 计 算 到 E 18 时 的 误 差 已 为 18!ϵ 6.0, 已 经 远 大 于 E 18 的 上 限 1/19. 现 在 我 们 考 虑 把 递 推 关 系 改 写 成 E n 1 = 1 E n, n =, 4, 3, 2. n 利 用 关 系 E n 1 n + 1 取 E 30 = 0.0 开 始 计 算, 此 时 误 差 为 ϵ 1/ , 作 一 次 递 推 后 得 到 E 29, 其 误 差 减 小 到 在 计 算 E 20 时, 误 差 已 减 小 到 20! 31! 准 确 到 14 位 有 效 数 字, 初 始 误 差 的 影 响 已 完 全 消 除 了. 而 每 一 步 的 舍 入 误 差, 都 在 其 后 的 迭 代 中 减 小 和 消 除 下 表 是 计 算 结 果. E 30 = 0.0 E 29 = E 28 = E 27 = E 26 = E 25 = E 24 = E 23 = E 22 = E 21 = E 20 = E 19 = E 18 = E 17 = E 16 = E 15 = E 14 = E 13 = E 12 = E 11 = E 10 = E 9 = E 8 = E 7 = E 6 = E 5 = E 4 = E 3 = E 2 = E 1 = 由 上 例 可 以 看 出 舍 入 误 差 对 计 算 的 影 响, 如 果 我 们 在 计 算 中 不 去 分 析 算 法, 即 使 编 出 了 正 确 的 程 序, 也 可 能 得 出 错 误 的, 有 时 甚 至 是 荒 谬 的 结 果.

12 12 第 一 章 计 算 机 数 和 误 差 1.3 减 法 的 计 算 也 许, 大 家 认 为 减 法 是 一 种 非 常 简 单 的 计 算, 但 是 在 用 计 算 机 做 减 法 时, 如 不 小 心, 就 很 容 易 出 错, 问 题 就 出 在 舍 入 误 差 上 面. 我 们 考 虑 下 面 的 例 子 : 对 于 大 的 x, 计 算 x + 1 x (1.1) 取 x = , 在 一 个 具 有 6 位 十 进 制 有 效 数 字 的 计 算 机 上, 计 算 结 果 为 0.001, 而 准 确 到 6 位 十 进 制 有 效 数 字 的 精 确 结 果 为 , 计 算 中 损 失 了 5 位 有 效 数 字. 对 于 更 大 的 x, 则 可 能 根 本 得 不 到 正 确 的 结 果. 但 是, 如 果 把 计 算 公 式 改 为 ( x + 1 x x 1 x) = (1.2) x x x x 就 可 以 得 到 较 好 的 结 果. 又 如 对 于 大 的 N, 计 算 N+1 N dx = artan(n + 1) artan(n) 1 + x2 时, 可 以 使 用 等 价 的 公 式 : ( ) 1 artan(n + 1) artan(n) = artan 1 + N(N + 1) 同 理, 对 于 小 的 ϵ ( sin (x + ϵ) sin(x) = 2 os x + ϵ ) ( ϵ sin 2 2) 上 面 分 析 的 是 一 些 简 单 的 例 子, 在 实 际 问 题 的 计 算 中, 常 常 会 碰 到 一 些 十 分 复 杂 的 问 题, 即 使 有 经 验 的 计 算 研 究 人 员 也 经 常 在 减 法 计 算 中 犯 错 误. 练 习 : 分 析 下 面 的 表 达 式, 给 出 对 任 何 x( 只 要 在 所 给 函 数 定 义 域 内 ) 都 能 得 出 正 确 结 果 的 计 算 公 式. os 2 x sin 2 x 1 + x2 1 x 2 ln(x) 1 e 2x e x 1 v2 1 x2 sin(x) x tg(x)

13 第 二 章 物 理 学 中 的 常 用 数 值 方 法 在 这 一 章 和 随 后 的 几 章 中, 我 们 将 介 绍 几 种 常 用 算 法, 这 些 内 容 应 该 在 数 学 课 中 讨 论, 但 目 前 的 数 学 课 程 并 不 包 含 这 些 属 于 计 算 数 学 课 程 的 内 容, 所 以 我 们 先 作 一 些 介 绍. 2.1 数 值 积 分 方 法 积 分 分 为 定 积 分 和 不 定 积 分 两 种, 我 们 在 这 里 只 讨 论 定 积 分, 不 定 积 分 的 计 算 可 作 为 变 上 限 的 定 积 分 来 计 算. 考 虑 定 积 分 I = 它 代 表 f(x) 曲 线 下 的 面 积. b a f(x)dx (2.1) 梯 形 法 和 辛 普 生 方 法 在 微 积 分 中, 我 们 知 道 积 分 是 由 求 和 的 极 限 定 义 的, 这 就 是 第 一 种 计 算 积 分 的 方 法 梯 形 法. 把 区 间 [a, b] 分 成 一 些 子 区 间 [a = x 0 < x 1 < x 2 < < x n = b], 则 (2.1) 的 积 分 I 可 近 似 为 : I = 1 n 1 (x i+1 x i )(f(x i ) + f(x i+1 )) (2.2) 2 i=0 在 实 际 计 算 中, 常 采 用 等 间 距 分 法, 即 x i+1 x i = h = b a n 为 一 常 数, 上 面 公 式 可 简 化 为 : I = h n 1 (f(x i ) + f(x i+1 )) 2 i=0 为 了 节 省 工 作 量, 实 际 使 用 的 计 算 公 式 为 : n 1 I = h f(x i ) + h 2 (f(x 0) + f(x n )) (2.3) i=1 当 子 区 间 越 来 越 小, h 0 时, 求 和 的 结 果 将 趋 于 积 分 值 由 于 我 们 事 先 并 不 知 道 积 分 的 结 果, 所 以 我 们 也 就 不 知 道 求 和 结 果 和 积 分 结 果 之 间 的 差 别 为 此, 我 们 需 要 估 计 计

14 14 第 二 章 物 理 学 中 的 常 用 数 值 方 法 算 误 差, 一 个 直 观 的 估 计 是 比 较 不 同 大 小 的 子 区 间 的 结 果 之 差 为 了 在 每 次 子 区 间 变 小 时 能 够 用 到 前 次 的 计 算 结 果, 我 们 可 以 从 n = 1 开 始 计 算, 每 次 把 区 间 数 增 加 一 倍, 则 前 一 次 计 算 的 函 数 结 果 仍 然 有 用, 通 过 比 较 连 续 两 次 的 计 算 结 果 之 差 是 否 小 于 给 定 的 误 差 限 来 结 束 计 算. 梯 形 方 法 实 际 上 把 每 一 个 子 区 间 上 的 函 数 近 似 为 线 性 函 数 αx + β, 然 后 计 算 其 积 分 值, 再 把 结 果 加 起 来.( 作 为 一 个 练 习, 请 自 行 证 明 ), 因 此, 我 们 说 梯 形 方 法 是 具 有 线 性 代 数 精 度 的 方 法. 下 面 考 虑 具 有 二 次 代 数 精 度 的 方 法, 辛 普 生 方 法. 先 考 虑 把 积 分 区 间 [a, b] 分 为 二 个 子 区 间, 分 点 为 (a, a + h, a + 2h = b), h = b a 2, 在 区 间 上 把 被 积 函 数 用 二 次 函 数 来 近 似, 即 取 f(x) α(x a) 2 + β(x a) + γ, 其 中 α, β, γ 可 利 用 分 点 上 的 函 数 值 f(a), f(a + h), f(b) 表 示 为 : 简 单 地 计 算 得 到 : b a α = β = f(b) 2f(a + h) + f(a) 2h 2 f(b) + 4f(a + h) 3f(a) 2h γ = f(a) (α(x a) 2 + β(x a) + γ) dx = h (f(b) + 4f(a + h) + f(a)) (2.4) 3 这 就 是 著 名 的 辛 普 生 求 积 公 式. 显 然, 它 具 有 二 次 代 数 精 度. 实 际 计 算 中, 把 被 积 区 间 [a, b] 先 分 为 N 个 子 区 间, 然 后 在 每 个 子 区 间 中 用 辛 普 生 求 积 公 式 计 算, 并 把 最 后 结 果 加 起 来. 其 计 算 公 式 为 : b a f(x)dx h (f(a) + 4f(a + h) + 2f(a + 2h) + 4f(a + 3h) + + f(b)) (2.5) 3 这 是 一 个 在 实 际 计 算 中 常 用 的 方 法, 它 对 于 一 般 的 有 限 区 间 的 积 分 都 能 给 出 较 好 的 正 确 结 果. 辛 普 生 方 法 有 各 种 不 同 的 实 现, 我 们 在 这 里 给 出 简 单 的 一 种, 取 每 个 子 区 间 均 相 等, 从 二 个 子 区 间 开 始, 每 次 子 区 间 数 目 加 倍, 这 样 可 以 在 每 次 计 算 时 使 用 前 次 的 计 算 值. 这 种 实 现 的 缺 点 是 在 整 个 积 分 区 间 均 匀 取 点, 对 于 在 部 分 区 间 变 化 剧 烈 而 在 其 余 地 方 变 化 平 缓 的 被 积 函 数 计 算 效 率 不 高. 这 一 实 现 可 以 充 分 利 用 它 与 梯 形 法 的 关 系. 假 定 n = 2 k, 用 S k 表 示 用 2 k 个 子 区 间 梯 形 公 式 求 得 的 结 果. 则 式 (2.5) 可 写 为 : b a f(x)dx 1 3 [4h(f(a + h) + f(a + 2h) + + f(b h) f(a) f(b)) 2h(f(a + 2h) + f(a + 4h) + + f(b 2h) f(a) f(b)) + h(f(a) + f(b)) 2h(f(a) + f(b)) + h(f(a) + f(b))] = 1 3 (4S k S k 1 ) (2.6)

15 2.1 数 值 积 分 方 法 15 计 算 定 积 分 还 有 很 多 方 法, 如 龙 贝 格 方 法, Newton-Cotes 方 法 等, 我 们 将 不 在 此 讨 论, 有 兴 趣 的 同 学 可 参 看 有 关 计 算 数 学 的 书 籍. Quadpak(netlib.org) 中 包 含 了 若 干 非 常 稳 定 的 求 积 分 的 程 序, 建 议 在 实 际 工 作 中 使 用 在 Matlab 中, 梯 形 法, 辛 普 森 均 作 为 标 准 函 数 提 供, 实 际 应 用 时, 只 要 写 一 个 被 积 函 数 的 m 程 序, 调 用 积 分 程 序 即 可 常 用 的 几 个 求 积 分 的 函 数 是 quad, quadl, quadgk 等 其 中 quad 使 用 自 适 应 的 辛 普 森 算 法,quadl 应 用 龙 贝 格 加 速 来 改 善 计 算 效 率 和 提 高 精 度 例 如, 计 算 积 分 1 先 定 义 函 数 fun 并 存 入 fun.m funtion y=fun(x) end y=sin(x.^2)./(1.+x.^4); 然 后 在 Matlab 的 命 令 窗 口 运 行 0 sin(x 2 ) 1 + x 4 dx >> format long >> quad(@fun,0,1) 结 果 为 : ans = 反 常 积 分 的 计 算 反 常 积 分 分 为 两 类, 一 类 是 积 分 区 间 有 限, 在 积 分 区 间 内 被 积 函 数 有 奇 点, 另 一 类 是 积 分 区 间 为 无 限. 对 于 二 者 兼 而 有 之 的 积 分, 可 以 分 为 两 个 或 多 个 积 分 来 处 理, 使 每 个 积 分 为 上 述 两 类 之 一. 积 分 区 间 内 含 有 奇 点 的 积 分 对 于 积 分 区 间 内 含 有 奇 点 的 积 分, 根 据 奇 点 的 性 质, 可 采 用 下 面 的 方 法 处 理. 可 去 奇 点, 设 x 0 为 一 可 去 奇 点, 且 知 道 在 x 0 点 被 积 函 数 的 极 限, 则 只 要 在 计 算 函 数 时, 在 x 0 的 一 个 邻 域 内 用 极 限 取 代 函 数 值 即 可. 例 如, sin(x) 在 x = 0 有 可 去 奇 点, 但 当 x x 0 时, sin(x) x 1 x2 6 +, 在 计 算 函 数 时, 可 使 用 下 面 的 语 句 :

16 16 第 二 章 物 理 学 中 的 常 用 数 值 方 法 IF(ABS(X).LT.1.0E-4) THEN F=1.0-X*X/6.0 ELSE F=SIN(X)/X ENDIF 极 限 方 法, 有 时 我 们 并 不 知 道 奇 点 处 函 数 的 极 限 表 达 式, 有 时 甚 至 不 知 道 奇 点 是 否 可 去, 此 时 可 用 极 限 方 法 来 逼 近. 例 如, 若 已 知 x 0 为 奇 点, 对 于 积 分 b x 0 f(x)dx = lim r x 0 b r f(x)dx 定 义 一 个 收 敛 于 x 0 的 序 列 b > r 1 > r 2 >, 例 如 可 取 r n = x n, 则 上 述 积 分 可 写 为 一 个 求 和 : b b r1 f(x)dx = f(x)dx + f(x)dx + f(x)dx + f(x)dx x 0 r 1 r 2 r 3 r 4 等 式 中 每 一 项 都 是 正 常 积 分, 当 r n r n+1 f(x)dx ϵ 时, 终 止 计 算. 如 果 上 式 不 收 敛, 则 很 可 能 原 积 分 不 存 在. 消 除 奇 点, 有 些 奇 点 可 以 通 过 变 量 替 换 或 对 被 积 函 数 进 行 变 换 而 消 去. 我 们 看 几 个 r2 r3 例 子 : 1 作 变 量 替 换 x = t 2, dx = 2t dt, 则 有 : 或 把 积 分 写 为 : 1 0 os(x) x dx = 1 1 第 二 项 积 分 中 的 奇 点 已 成 为 可 去 奇 点 os(x) x dx os(x) x dx = os(t 2 )dt 1 1 os(x) 1 1 dx + dx = 2 + x 0 x 0 I(x) = x 0 e t 1 t 在 t = 1 的 邻 域, 被 积 函 数 与 e 1 /(1 t) 很 相 像, 因 此 变 为 : I(x) = x 0 e 1 x 1 t dt + 0 所 有 的 奇 异 性 都 在 第 一 项, 而 这 一 项 可 解 析 处 理. os(x) 1 dx x e t e 1 dt = 1 x 1 t e ln 1 x + e t e 1 dt 1 t 除 上 面 例 子 演 示 的 方 法 外, 还 可 用 分 部 积 分 法 消 除 奇 点, 这 里 不 再 一 一 举 例 了. 0

17 2.1 数 值 积 分 方 法 17 积 分 区 间 为 无 限 的 积 分 下 面 我 们 讨 论 积 分 区 间 为 无 限 的 情 形. 为 具 体 起 见, 只 讨 论 形 为 0 f(x)dx 的 积 分, 其 它 形 式 的 无 穷 区 间 的 积 分 均 可 化 为 上 面 的 形 式. 变 量 替 换 法 : 通 过 变 量 替 换, 可 把 无 穷 区 间 的 积 分 化 为 有 限 区 间 上 的 积 分, 如 令 x = ln(t), dx = dt t, 则 有 : 0 f(x)dx = 1 0 f( ln(t) 1 dt = t 0 g(t) dt (2.7) t g(t) 这 里 g(t) = f( ln(t)), 若 在 t = 0 的 邻 域 内 有 界, 则 上 式 的 积 分 成 为 一 个 正 常 积 分, t 否 则 就 是 前 面 讨 论 过 的 有 限 区 间 上 的 反 常 积 分, 可 以 用 前 面 的 方 法 来 处 理. 上 述 变 换 对 于 当 x 时, f(x) 以 e kx 的 形 式 衰 减 的 函 数, 计 算 结 果 十 分 令 人 满 意. 其 它 常 用 的 变 换 还 有 x = 以 及 t = tanh x 等. t, dx = 1 dt 1 t (1 t) 2 极 限 法 : 由 定 义 0 f(x)dx = f ( ) t 1 dt (2.8) 1 t (1 t) 2 r f(x)dx = lim f(x)dx r 0 定 义 一 个 趋 于 的 序 列 0 < r 1 < r 2 <, 例 如 可 取 r n = 2 n, 则 上 述 积 分 可 写 为 一 个 求 和 : 0 f(x)dx = r1 0 f(x)dx + r2 r 1 f(x)dx + r3 r 2 f(x)dx + (2.9) 式 中 每 一 项 都 是 正 常 积 分, 当 r n+1 f(x)dx ϵ 时, 终 止 计 算. 如 果 上 式 不 收 敛, 则 很 可 能 原 积 分 不 存 在. 如 果 能 找 到 尾 项 r n R 需 进 行 到 某 一 足 够 大 的 r n 处, 再 加 上 尾 项. f(x)dx 的 一 个 较 好 的 渐 近 表 达 式, 则 上 述 积 分 只 还 有 一 些 其 它 的 有 效 方 法, 这 里 不 再 一 一 列 举, 我 们 下 面 就 要 讲 到 的 高 斯 积 分 方 法 提 供 了 另 一 种 处 理 反 常 积 分 的 工 具 高 斯 积 分 方 法 前 面 的 积 分 公 式 ( 梯 形 公 式 及 辛 普 生 公 式 ) 的 代 数 精 度 都 较 低, 一 个 自 然 的 问 题 是, 能 否 构 造 出 高 代 数 精 度 的 公 式, 回 答 是 肯 定 的, 这 就 是 高 斯 积 分 方 法. 我 们 将 不 讨 论 高 斯 积 分 方 法 的 理 论 基 础, 而 只 是 给 出 有 关 的 计 算 公 式, 对 理 论 基 础 有 兴 趣 的 同 学 可 以 参 看 有 关 的 计 算 数 学 的 书 籍.

18 18 第 二 章 物 理 学 中 的 常 用 数 值 方 法 有 限 区 间 的 高 斯 型 积 分 公 式 : 高 斯 勒 让 德 积 分 公 式 : 1 1 f(x)dx = n w i f(x i ) (2.10) 此 处, x i 为 n 阶 勒 让 德 多 项 式 的 根, 通 常 称 为 积 分 节 点, w i 称 为 积 分 权 重. 等 权 重 切 贝 雪 夫 积 分 公 式 : 这 里 x i 为 积 分 节 点, 权 重 为 2/n. 1 第 一 类 高 斯 切 贝 雪 夫 积 分 公 式 : f(x)dx = 2 n f(x) 1 x 2 dx = π n n f i=1 i=1 n f(x i ) (2.11) i=1 ( os ( )) (2i + 1)π 2n (2.12) 第 二 类 高 斯 切 贝 雪 夫 积 分 公 式 : 1 1 f(x) 1 x 2 dx = π n + 1 n ( ) ( ( )) iπ iπ sin 2 f os n + 1 n + 1 i=1 如 果 积 分 区 间 不 是 [ 1, 1], 而 是 [a, b], 可 利 用 变 换 ( ) ( ) b a b + a x = t (2.13) (2.14) 将 积 分 区 间 变 为 [ 1, 1], 然 后 利 用 上 述 公 式 进 行 计 算. 在 实 际 计 算 时, 常 把 积 分 区 间 分 为 几 个 小 区 间, 在 每 一 个 小 区 间 内 使 用 高 斯 积 分 公 式, 并 把 结 果 加 起 来. 练 习 : 试 编 写 上 述 四 种 高 斯 积 分 公 式 的 程 序, 并 利 用 所 编 写 的 程 序 计 算 积 分 1 1 x + 1.5dx, x dx, 1 0 dx x 2 0 e os2 (x) dx, ( 积 分 所 需 的 节 点 及 权 重 在 后 面 的 表 中 给 出 ). 1 0 sin(x) x dx 无 限 区 间 的 高 斯 型 积 分 公 式 : 高 斯 拉 盖 尔 积 分 公 式 : 0 e x f(x)dx = n w i f(x i ) (2.15) i=1

19 2.1 数 值 积 分 方 法 19 高 斯 厄 米 积 分 公 式 : e x2 f(x)dx = n w i f(x i ) (2.16) i=1 这 里 x i 分 别 为 n 阶 拉 盖 尔 多 项 式 和 厄 米 多 项 式 的 根. 其 数 值 和 对 应 的 权 重 的 数 值 在 附 录 的 表 中 给 出. 练 习 : 试 编 写 上 述 二 种 高 斯 积 分 公 式 的 程 序, 并 利 用 所 编 程 序 计 算 积 分 0 e x 1 + x 4 dx, 0 e x2 os(x)dx 各 阶 高 斯 积 分 公 式 的 取 点 并 不 重 和, 当 增 加 分 点 时, 所 有 的 函 数 值 都 要 重 算 1965 年 Kronrod 提 出 了 结 合 高 斯 方 法 和 加 速 收 敛 技 巧 的 算 法 这 一 算 法 也 可 以 处 理 有 一 定 奇 异 性 的 被 积 函 数, 但 需 要 把 奇 点 放 在 积 分 限 上 在 Matlab 中 已 实 现 这 一 算 法, 其 调 用 命 令 是 quadgk. 例 如, 计 算 积 分 定 义 函 数 funsin, 并 存 入 funsin.m 0 sin x x e x dx funtion y=funsin(x) end y=sin(x).*exp(-x)./x; 在 Matlab 命 令 窗 口 键 入 如 下 命 令 format long quadgk(@funsin,0,inf) ans = 这 一 积 分 的 精 确 结 果 是 π 4 =

20 20 第 二 章 物 理 学 中 的 常 用 数 值 方 法 哥 西 主 值 的 计 算 在 物 理 计 算 中, 哥 西 主 值 的 计 算 占 了 很 重 要 的 位 置, 如 物 质 的 介 电 函 数 的 实 部 和 虚 部 之 间 满 足 著 名 的 Kramers-Kronig 关 系 : ϵ(ω) = ϵ 1 (ω) + iϵ 2 (ω) ϵ 1 (ω) = π P sϵ 2 (s) s 2 ω ds 2 ϵ 2 (ω) = 2ω π P 0 0 ϵ 1 (s) ds (2.17) s 2 ω2 这 个 关 系 是 一 个 严 格 的 关 系, 它 来 源 于 因 果 关 系, 是 非 常 普 遍 的. 在 物 理 上, ϵ 2 (ω) 的 测 量 和 计 算 都 比 较 容 易, 因 此, 一 般 不 直 接 测 量 或 计 算 ϵ 1 (ω), 而 是 利 用 上 述 关 系 由 ϵ 2 (ω) 的 数 据 进 行 计 算. 上 述 Kramers-Kronig 关 系 中 出 现 的 P 就 代 表 哥 西 主 值 积 分, 设 a < < b, f(x) 在 x = 的 邻 域 内 无 界, 则 哥 西 主 值 积 分 定 义 为 : P b a f(x)dx = lim r 0 [ r a f(x)dx + a 不 失 一 般 性, 我 们 取 = 0 并 将 积 分 取 成 f(x)dx 的 形 式. 令 a b +r ] f(x)dx (2.18) g(x) = 1 [f(x) f( x)], 2 h(x) = 1 [f(x) + f( x)] 2 则 f(x) = g(x) + h(x) 由 于 g(x) 是 奇 函 数 而 h(x) 是 偶 函 数, 于 是 : 因 此, P b a = = 2 r a r a a r f(x)dx = 2 lim f(x)dx + g(x)dx + h(x)dx r 0 a r a r a r f(x)dx g(x)dx + h(x)dx = 2 a 0 r a h(x)dx = h(x)dx + a 0 a r h(x)dx [f(x) + f( x)]dx 这 样 我 们 就 把 哥 西 主 值 积 分 的 计 算 变 为 在 x = 0 处 有 奇 点 的 一 个 反 常 积 分 的 计 算. h(x) 也 可 能 在 x = 0 处 无 奇 点, 此 时 只 要 计 算 一 个 常 规 积 分. 例 如 前 面 的 Kramers-Kronig 关 系 中 积 分 的 计 算 可 化 为 : 2 π P 0 sϵ 2 (s) s 2 ω ds = 1 2 2π P ϵ 2 (s) s ω ds = 1 π 0 ϵ 2 (ω + x) ϵ 2 (ω x) dx x

21 2.1 数 值 积 分 方 法 21 2ω π P ϵ 1 (s) 0 s 2 ω ds = 1 2 2π P ϵ 1 (s) s ω ds = 1 ϵ 1 (ω + x) ϵ 1 (ω x) dx π 0 x 在 得 到 上 述 关 系 时, 我 们 利 用 了 ϵ 1 为 偶 函 数 而 ϵ 2 为 奇 函 数 的 事 实. 当 哥 西 积 分 区 间 内 的 无 界 邻 域 不 止 一 个 时, 可 以 把 积 分 区 间 分 为 几 段, 使 得 每 一 段 中 只 包 含 一 个 无 界 邻 域, 然 后 在 每 一 段 用 上 述 方 法 计 算. 分 另 一 种 计 算 主 值 积 分 的 方 法 是 直 接 从 定 义 出 发, 设 f(x) 在 积 分 区 间 内 无 奇 点, 则 积 b P a f(x) x x 0 dx = x0 f(x) b f(x) x0 + f(x) dx + dx + P dx a x x 0 x 0 + x x 0 x 0 x x 0 = I + lim I ϵ, ϵ 0 (2.19) I 为 第 一 行 右 方 第 一 和 第 二 项 之 和, 是 正 常 积 分, 可 用 前 述 任 一 方 法 计 算. I ϵ, = x0 ϵ x 0 把 f(x) 在 x 0 点 展 开, 并 逐 项 积 分, 得 到 取 极 限 ϵ 0, 得 到 f(x) x x 0 dx + x0 + x 0 +ϵ f(x) x x 0 dx I ϵ, = 2f (x 0 )( ϵ) + f (3) (x 0 )( 3 ϵ 3 ) + I 0,δ = 2 这 样, 原 主 值 积 分 为 I 0, + I. 练 习 : f (2n+1) 2n+1 (x 0 ) (2n + 1)!(2n + 1). (2.20) n=0 试 用 上 述 方 法 计 算 主 值 积 分 1 e x P 0 x 0.5, 对 于 不 同 的, 在 (2.20) 中 取 两 项, 比 较 计 算 结 果 与 精 确 结 果. ( 到 16 位 的 精 确 结 果 为 ) 另 外 一 种 常 用 的 计 算 主 值 积 分 的 方 法 是 基 于 如 下 公 式 1 lim η 0 + x iη = P 1 + iπδ(x) (2.21) x 式 中 η 0 + 表 示 η 从 大 于 0 的 方 向 趋 于 0. 这 一 公 式 当 然 是 在 积 分 的 意 义 下 理 解 利 用 上 述 结 果, 就 可 以 得 到 计 算 主 值 积 分 的 一 个 方 法 : 例 如 要 计 算 P b a F (x) dx b 取 一 个 小 量 η, 计 算 F (x iη) dx, 其 实 部 给 出 P b F (x) dx 的 依 赖 于 η 的 近 似 值, 通 a a 过 减 小 η, 可 以 得 到 精 度 要 求 下 收 敛 的 结 果 例 如, 对 于 前 面 的 练 习 题, 分 别 取 不 同 的 η, 计 算 结 果 如 下

22 22 第 二 章 物 理 学 中 的 常 用 数 值 方 法 η 积 分 值 i i i i 急 速 振 荡 函 数 的 积 分 换 在 各 种 变 换 特 别 是 付 里 叶 变 换 中, 经 常 会 遇 到 急 速 振 荡 函 数 的 积 分, 例 如, 付 里 叶 变 b 又 如 付 里 叶 贝 塞 尔 变 换, a f(x) os(nx)dx, 1 这 里 0 < γ 1 < γ 2 < 是 贝 塞 尔 函 数 的 根. 一 般 我 们 可 把 这 类 函 数 写 为 : I(t) = b a 0 b a f(x)xj n (γ m x)dx f(x) sin(nx)dx f(x)k(x, t)dx, < a < b < (2.22) 式 中 K(x, t) 为 一 振 荡 核 而 f(x) 为 非 振 荡 部 分. 由 于 振 荡 型 函 数 在 积 分 区 间 内 多 次 取 几 乎 相 等 的 正 值 和 负 值, 如 果 用 通 常 的 方 法 计 算, 则 由 于 正 负 数 相 加, 造 成 有 效 数 字 的 损 失, 几 乎 得 不 到 正 确 的 结 果. 下 面 先 给 出 最 常 遇 到 的 振 荡 函 数 的 积 分 付 里 叶 系 数 的 一 种 计 算 方 法, 然 后 再 讨 论 较 为 一 般 的 方 法. 付 里 叶 系 数 的 计 算 : 付 里 叶 系 数 的 复 数 形 式 为 : b 对 上 式 分 部 积 分, 假 定 f(x) 存 在 直 到 n 阶 的 导 数, 可 以 得 到 : b a n 1 f(x)e isx dx = e isa k=0 a f(x)e isx dx (2.23) ( ) k+1 i n 1 f (k) (a) e isb s 当 n 趋 于 无 穷 大 时, 余 项 R n 趋 于 0, 而 求 和 项 就 给 出 了 积 分 的 近 似 值. 练 习 : 写 出 (2.24) 的 实 部 和 虚 部, 即 k=0 ( ) k+1 i f (k) (b) + R n (2.24) s b a f(x) os sxdx 和 b a f(x) sin sxdx 的 积 分 公 式.

23 2.1 数 值 积 分 方 法 23 上 述 公 式 需 要 计 算 f(x) 的 直 到 n 1 阶 的 导 数, 这 对 于 较 为 复 杂 的 函 数 来 说, 曾 是 一 个 很 重 的 负 担. 不 过 目 前 计 算 机 代 数 的 发 展, 计 算 导 数 已 经 没 有 多 少 困 难. 为 了 看 出 这 一 方 法 的 效 力, 我 们 来 看 一 个 例 子. 计 算 I = 2π 0 x os(x) sin 30xdx 用 公 式 (2.24) 的 虚 部, 取 n = 3, 由 f(x) = x os(x), a = 0, b = 2π, f(0) = 0, f(2π) = 2π, f (0) = 0, f (2π) = 2π, 代 入 得 到 : 2π I = x os(x) sin 30xdx = 1 ( ) π + ( 2π) = 而 该 积 分 的 精 确 值 为 基 于 Euler-Malaurin 公 式 的 梯 形 法 : 可 以 证 明 ( 见 王 竹 溪, 郭 敦 仁 特 殊 函 数 概 论 第 一 章 ), 函 数 f(x) 的 积 分 可 以 写 成 a+mh [ ] f(a) f(a + mh) f(x)dx = h + f(a + h) + + f(a + (m 1)h) + a 2 2 n ( 1) k B k h 2k [ + f 2k 1 (a + mh) F 2k 1 (a) ] + R n (2k)! 其 中 k=1 R n = θ ( 1)n+1 B n+1 h [ (2n+1) f (2n+1) (a + mh) f (2n+1) (a) ] (2n + 2)! (2.25) 这 里 0 θ 1,B n 是 Bernoulli 数 如 果 f(x) 是 周 期 函 数,f(a + mh) = f(a), 且 f(x) 存 在 直 到 (2n + 1) 阶 的 导 数, 则 显 然 有 f (2k+1) (a) = f (2k+1) (a + mh), 于 是, 利 用 梯 形 法 可 以 得 到 精 确 的 结 果 作 为 例 子, 取 不 同 的 点, 试 用 梯 形 法 计 算 如 下 积 分, 并 于 精 确 结 果 比 较 1 Bessel 函 数 的 积 分 表 示 之 一 是 sin 10πxdx, 精 确 值 = J n (t) = 1 π π 计 算 n = 2, t = 1, 3, 7, 10 时 Bessel 函 数 的 值 0 os (t sin x nx)dx 由 Euler-Malaurin 公 式 可 知, 如 果 一 个 函 数 满 足 f (a) = f (a + mh), f (3) (a) = f (3) (a+mh),, 则 梯 形 公 式 也 可 以 得 到 很 精 确 的 结 果 例 如 函 数 f(x) = exp(x 2 (1 x 2 )), 在 0 和 1 满 足 上 述 关 系, 从 而, 可 以 用 梯 形 公 式 计 算 1 0 exp(x 2 (1 x 2 ))dx 这 个 积 分 的 精 确 结 果 是 , 请 分 别 用 梯 形 法 和 辛 普 森 法 计 算 并 比 较 其 结 果

24 24 第 二 章 物 理 学 中 的 常 用 数 值 方 法 子 区 间 法, 在 一 般 情 况 下, 如 果 K(x, t) = 0 的 根 容 易 求 出, 且 记 为 a x 1 < x 2 < < x n b, 则 可 用 普 通 的 积 分 规 则 计 算 每 一 个 子 积 分 xi+1 x i, 并 把 结 果 加 起 来. 用 这 种 方 法 得 到 的 一 般 是 一 个 交 错 级 数, 且 相 邻 项 的 绝 对 值 比 较 接 近, 其 求 和 可 用 Euler 变 换 计 算. (Euler 变 换 将 在 后 面 讨 论 ) 在 讲 了 样 条 插 值 之 后, 我 们 再 介 绍 一 种 基 于 样 条 插 值 的 计 算 积 分 (2.23) 的 方 法. 更 多 的 方 法 可 以 参 看 (P. J. Davis and P. Rabinowitz,Methods of Numerial Integration, Aademi Press, 1984) 高 维 积 分 的 计 算 高 维 积 分 的 计 算 原 则 上 可 以 通 过 化 成 累 次 积 分 来 进 行, 例 如, 对 于 积 分 I = g(x, y, z) dxdydz (2.26) 其 中 Ω 为 积 分 区 域, 总 可 以 化 为 I = b dx Ω d(x) dy f(x,y) a (x) e(x,y) g(x, y, z) (2.27) 这 相 当 于 依 次 计 算 三 个 一 重 积 分, 每 一 积 分 以 其 前 一 积 分 作 为 被 积 函 数. 稍 加 考 虑 就 可 发 现, 这 一 过 程 的 计 算 量 随 积 分 重 数 的 增 加 而 变 得 无 法 完 成. 如 每 重 积 分 取 10 个 点 ( 这 是 一 个 非 常 小 的 数 字 ), 则 n 重 积 分 就 要 计 算 10 n 个 函 数 值. 因 此, 这 一 方 法 只 能 用 于 积 分 重 数 较 低 的 积 分. 高 维 积 分 还 可 以 用 高 斯 方 法 计 算, 这 只 要 把 每 重 积 分 用 高 斯 积 分 公 式 代 入 即 可. 对 于 n > 3 的 高 维 积 分, 计 算 量 总 是 相 当 大 的, 如 果 你 对 计 算 精 度 要 求 不 高, 那 么 可 以 利 用 后 面 将 要 讲 到 的 Monte Carlo 方 法. 另 外, 如 果 可 能, 请 尽 可 能 地 解 析 积 出 内 层 积 分, 以 减 小 计 算 量 固 体 热 容 量 的 计 算 在 德 拜 近 似 下, 固 体 的 热 容 量 由 下 式 给 出 ( ) 3 T Θd /T C V = 9Nk Θ D 0 ξ 4 e ξ (e ξ 1) 2 dξ 式 中 Θ D 为 德 拜 温 度. 为 了 求 得 热 容 量, 需 要 计 算 上 式 中 的 积 分, 这 一 积 分 只 有 当 Θ d /T 1 或 Θ d /T 1 时 才 能 解 析 积 出, 一 般 情 况 下 只 能 进 行 数 字 计 算. 考 虑 f(x) = x 0 ξ 4 e ξ (e ξ 1) 2 dξ

25 2.2 方 程 的 求 根 和 最 优 化 方 法 25 利 用 我 们 学 到 的 方 法, 可 以 对 不 同 的 x 求 出 其 值, 并 做 出 f(x) 的 表 或 图 形. 为 了 计 算 积 分, 注 意 到 ξ = 0 是 被 积 函 数 的 可 去 奇 点, 因 此, 应 在 该 点 把 被 积 函 数 展 开 为 ξ 的 级 数 ; 当 ξ, 被 积 函 数 指 数 减 小. 已 知 f( ) = 0 ξ 4 e ξ (e ξ 1) dξ = π4 我 们 可 以 根 据 x 的 值 的 大 小 分 别 用 两 种 方 法 计 算, 取 一 值 x, 当 x < x 时, 用 原 式 计 算, 当 x > x 时, 可 以 计 算 f( ) f 1 (x), 其 中 f 1 (x) = x ξ 4 e ξ (e ξ 1) 2 dξ x 可 以 选 择 的 大 一 点, 这 样 f 1 (x) 的 数 值 比 较 小, 计 算 精 度 要 求 可 以 低 一 些. 我 们 选 择 x = 10, x < x 时, 用 辛 普 生 方 法 计 算 f(x), 被 积 函 数 为 { ξ 2 ξ ξ < 0.05 ξ 4 e ξ (e x 1) ξ x 当 x > x 时, 计 算 f 1 (x), 对 积 分 作 变 量 变 换 ξ = η + x 则 f 1 (x) = e x 0 (η + x) 4 e η (1 e (η+x) ) 2 dη 这 一 积 分 可 用 高 斯 拉 盖 尔 方 法 计 算. 作 为 练 习, 请 同 学 完 成 余 下 的 计 算 并 作 出 f(x) x 的 图 形 2.2 方 程 的 求 根 和 最 优 化 方 法 方 程 的 求 根 问 题 是 物 理 计 算 中 最 常 遇 到 的 问 题 之 一, 对 于 低 次 代 数 方 程, 已 经 找 到 了 解 的 公 式, 但 除 了 二 次 方 程 之 外, 这 些 公 式 并 无 多 少 实 用 价 值, 而 对 于 即 使 如 os(x) x = 0 这 样 简 单 的 方 程, 也 无 法 找 到 解 析 形 式 的 解. 数 值 计 算 几 乎 是 求 解 任 何 有 实 际 意 义 的 问 题 的 唯 一 方 法. 在 这 一 节 中, 我 们 将 讨 论 几 种 求 解 f(x) = 0 的 根 的 数 字 方 法. 最 优 化 方 法 除 了 本 身 的 重 要 性 之 外, 也 是 一 种 求 解 多 元 方 程 的 有 效 方 法, 因 此, 最 优 化 方 法 的 介 绍 也 是 本 节 的 一 个 重 要 内 容 二 分 法 二 分 法 是 求 解 方 程 的 最 安 全, 也 是 概 念 上 最 简 单 的 方 法, 它 可 以 做 为 后 面 将 要 讲 到 的 方 法 失 败 后 的 最 后 保 险. 假 定 我 们 知 道 在 区 间 a < x < b 之 内 只 有 f(x) = 0 的 一 个 根, 则 显 然 f(a) 与 f(b) 异 号, 把 区 间 分 半, 即 令 = (b + a)/2, 计 算 f() 的 值, 若 f(a) 与 f() 同 号, 则 根 在 [, b] 之 间, 否 则, 根 在 [a, ] 之 间, 对 根 所 在 的 区 间 重 复 上 述 过 程 直 到 达 到 所 需 精 度.

26 26 第 二 章 物 理 学 中 的 常 用 数 值 方 法 牛 顿 法 和 割 线 法 在 所 有 一 元 函 数 求 根 方 法 中, 牛 顿 法 也 许 是 最 值 得 推 荐 的 一 个. 其 想 法 是 这 样 的, 假 定 f(x) = 0 的 根 为 x, 而 我 们 有 一 个 对 x 的 猜 测 值 x n, f(x) 可 在 x n 附 近 展 开 为 泰 勒 级 数 : f(x) = f(x n ) + f (x n )(x x n ) + 略 去 高 次 项, 令 上 式 为 0, 可 解 出 x, 将 此 值 作 为 x 的 一 个 更 好 的 近 似, 记 为 x n+1, 则 x n+1 = x n f(x n) f (x n ). (2.28) 这 样, 给 定 一 个 x 的 初 始 值 x 0, 利 用 下 面 的 迭 代 公 式 反 复 迭 代 : x n+1 = x n f(x n) f (x n ) (2.29) 如 果 迭 代 收 敛, 则 必 收 敛 于 f(x) = 0 的 一 个 根 x. 这 一 迭 代 公 式 (2.29) 就 是 著 名 的 牛 顿 迭 代 公 式. 可 以 证 明, 如 果 迭 代 到 第 n 步 时 的 误 差 为 ϵ n = x n x, 则 第 n + 1 步 的 误 差 为 f (x ) ϵ n+1 = ϵ 2 n, 也 就 是 说, 牛 顿 法 是 平 方 收 敛 的. 证 明 如 下 : 由 迭 代 公 式 (2.29), 2f (x ) x n+1 = x n f(x n) f (x n ) x + ϵ n ϵ nf (x ) ϵ2 nf (x ) f (x ) + ϵ n f (x ) ( x + ϵ n ϵ n + 1 f (x ) 2 f (x ) ϵ2 n f ) (x ) f (x ) ϵ2 n = x + 1 f (x ) 2 f (x ) ϵ2 n 二 分 法 是 线 性 收 敛 的, 请 证 明 之. 但 是, 另 一 方 面, 对 于 任 一 初 始 值, 牛 顿 法 并 不 总 是 收 敛, 事 实 上, 牛 顿 法 的 收 敛 范 围 是 较 小 的, 因 此, 在 实 际 使 用 时, 必 须 找 一 个 较 好 的 x 0, 使 x 0 尽 可 能 地 靠 近 x, 这 可 以 把 二 分 法 与 牛 顿 法 结 合 起 来 以 完 成. 试 计 算 10. 令 f(x) = x 2 10, f(x) = 0 的 根 就 是 所 求 结 果 取 初 始 值 x 0 = 3, f (x) = 2x. 迭 代 结 果 如 下 :

27 2.2 方 程 的 求 根 和 最 优 化 方 法 27 牛 顿 法 的 一 个 重 要 的 缺 点 是 需 要 计 算 函 数 的 导 数, 对 于 一 些 复 杂 的 函 数 来 说, 这 是 一 件 十 分 麻 烦 的 工 作. 为 了 即 保 持 较 高 的 收 敛 速 度, 又 避 免 麻 烦 的 导 数 计 算, 可 以 对 牛 顿 法 稍 作 变 形, 用 差 分 代 替 导 数, 从 而 得 到 下 面 的 割 线 法. x n+1 = x n x n x n 1 f(x n ) f(x n 1 ) f(x n) (2.30) 这 一 迭 代 方 式 需 要 二 个 初 始 值 x 0, x 1, 可 以 取 为 所 猜 测 的 解 附 近 的 二 个 点. 割 线 法 的 收 敛 速 度 比 牛 顿 法 稍 慢, 但 省 去 了 导 数 的 计 算, 对 于 求 导 数 比 较 困 难 的 函 数 特 别 有 用. 作 为 一 个 编 程 练 习, 请 写 出 这 一 算 法 的 程 序 并 作 为 你 自 己 的 库 程 序. 练 习 : 证 明 割 线 法 的 收 敛 阶 为 ϵ n+1 ϵ n 另 一 个 不 需 要 计 算 导 数 的 方 法 是 所 谓 的 Stewenson 方 法, 算 法 如 下 : x n+1 = x n f(x n ) 2 f(x n ) f(x n f(x n )) (2.31) 对 于 单 零 点, 这 一 算 法 的 收 敛 阶 为 2,( 请 证 明 ), 而 且 不 需 要 计 算 导 数, 是 一 个 较 好 的 算 法. 试 计 算 10. 令 f(x) = x 2 10, f(x) = 0 的 根 就 是 所 求 结 果 取 初 始 值 x 0 Stewenson 法 迭 代 结 果 如 下 : = 3. 用 一 维 最 优 化, 黄 金 分 割 法 这 一 节 开 始 我 们 讨 论 最 优 化 方 法, 或 计 算 一 目 标 函 数 的 极 值 的 方 法. 最 优 化 方 法 在 物 理 学 中 有 十 分 广 泛 的 应 用, 函 数 的 求 根 问 题 也 可 化 为 一 个 最 优 化 问 题, 例 如, 求 函 数 f(x) = 0 的 问 题 就 与 求 目 标 函 数 为 F (x) = f(x) 2 的 极 小 值 的 问 题 等 价, 而 方 程 组 f i (x 1, x 2, x 3,, x n ) = 0, i = 1, 2, 3,, n 的 求 解 问 题 则 与 F (x 1, x 2,, x n ) = n w i f i (x 1, x 2,, x n ) 2 i=1

28 28 第 二 章 物 理 学 中 的 常 用 数 值 方 法 的 最 优 化 问 题 等 价. 这 里 w i, i = 1, 2,, n 为 一 组 权 重 参 数, 适 当 的 选 则 可 使 计 算 更 为 容 易 进 行, 但 对 根 的 位 置 没 有 影 响. 这 一 节 我 们 先 考 虑 一 个 自 变 量 的 目 标 函 数 的 最 优 化 问 题. 假 定 我 们 有 一 目 标 函 数 F (x), 我 们 的 任 务 是 寻 找 一 个 x = x, 使 得 目 标 函 数 F (x) 最 小. 这 一 工 作 可 以 分 两 步 来 完 成, 首 先, 我 们 要 找 出 目 标 函 数 的 极 小 点 所 在 的 大 致 位 置, 或 者 更 确 切 一 些, 我 们 要 找 三 个 点 a < b <, 使 得 目 标 函 数 F (b) < F (a), F (b) < F (). 第 二 步 则 在 区 间 [a, ] 内 寻 找 F (x) 的 最 小 值. 第 一 步 原 则 上 是 十 分 简 单 的, 假 定 有 两 点 a, b, 我 们 可 以 比 较 F (a) 和 F (b) 以 确 定 一 个 函 数 减 小 的 方 向, 不 失 一 般 性, 假 定 F (b) < F (a), 然 后 沿 减 小 方 向 前 进 一 步, 到 达 新 的 一 点, 一 般 取 = b (b a), 如 果 F () < F (b), 则 令 a = b, b =, 继 续 前 进, 寻 找 新 的 ; 如 果 F () > F (b), 则 我 们 的 目 的 已 经 达 到. 下 面 考 虑 第 二 步 的 计 算, 即 从 三 个 点 a < b <, F (b) < F (a), F (b) < F () 出 发, 寻 找 F (x) 的 极 小 点 x min. 这 里 介 绍 黄 金 分 割 法. 寻 找 极 小 点 的 过 程, 实 际 上 就 是 用 一 种 方 法 不 断 缩 小 区 间 [a, ], 并 保 证 极 小 点 一 直 位 于 区 间 之 内. 这 一 过 程 可 以 这 样 来 完 成, 假 定 已 经 作 了 n 步, 得 到 三 个 点 a, b, 和 对 应 的 目 标 函 数 值 F (a) > F (b), F (b) < F (), 不 失 一 般 性, 假 定 b > b a, 在 区 间 [b, ] 内 选 一 点 x, 如 果 F (x) < F (b), 则 作 代 换 a = b, b = x, 新 的 缩 小 了 的 区 间 为 [a, ]( 即 原 [b, ], 丢 掉 的 区 间 为 原 [a, b]); 如 果 F (x) > F (b), 则 作 代 换 = x, 新 的 缩 小 了 的 区 间 为 [a, ]( 即 原 [a, x], 丢 掉 的 区 间 为 原 [x, ]). 问 题 在 于 如 何 选 择 x, 一 个 自 然 的 想 法 是, 不 论 F (x) < F (b) 或 F (x) > F (b), 丢 掉 的 区 间 都 相 同, 即 选 择 x, 使 得 b a = x w( a) (2.32) 这 里 w 是 一 个 比 例 因 子, 我 们 现 在 来 确 定 它. 由 于 要 求 每 一 步 都 满 足 上 述 关 系, 显 然 有 x b = w( b) = w(( a) (b a)) = w( a) w 2 ( a) = (w w 2 )( a) (2.33) 另 一 方 面, 由 (2.32) 可 得 (b a) + ( x) = ( a) (x b) = 2w( a) (2.34) 即 由 (2.33) 和 (2.35), 我 们 得 到 : x b = (1 2w)( a) (2.35) w 2 3w + 1 = 0 (2.36) 由 此 解 得, w = , 1 w w 通 常 被 称 为 黄 金 分 割 数 或 黄 金 数. 一 般 在 开 始 计 算 时, a, b, 三 者 并 不 满 足 黄 金 分 割 关 系, 但 几 次 迭 代 后, 这 一 关 系 便 可 满 足, 并 在 其 后 一 直 保 持.

29 2.2 方 程 的 求 根 和 最 优 化 方 法 高 维 最 优 化, 共 扼 梯 度 方 法 这 一 节 给 出 多 维 问 题 的 最 优 化 方 法. 对 于 一 个 含 有 n 个 自 变 量 的 目 标 函 数 F (x 1, x 2, x n ), 我 们 的 目 的 是 要 求 出 一 组 (x 1, x 2,, x n) 使 目 标 函 数 达 到 极 小. 为 了 讨 论 问 题 的 方 便, 下 面 用 x 代 表 n 维 空 间 的 一 个 点, 最 优 化 问 题 也 可 以 表 述 为, 寻 找 n 维 空 间 的 一 点 x 使 得 目 标 函 数 F (x ) 最 小. 高 维 最 优 化 问 题 有 很 多 求 解 方 法, 其 基 本 思 路 是 相 同 的, 这 就 是 从 一 点 x 0 出 发, 按 照 某 种 规 定 的 方 向 p 0, 用 一 维 最 优 化 的 方 法 寻 找 函 数 F (x) 的 极 小 点 x 1, 继 续 这 一 过 程 直 到 达 到 最 优 点. 若 第 i 步 到 达 x i, 设 选 定 方 向 p i, 则 新 的 x i+1 可 如 下 求 得. 实 际 上, 我 们 只 要 找 一 参 数 t i, 使 得 F (x i + t i p i ) = min t F (x i + tp i ) (2.37) 则 x i+1 = x i + t i p i 便 为 所 求. 因 此, 高 维 最 优 化 的 的 不 同 方 法 实 际 上 是 选 择 不 同 的 一 维 寻 找 方 向. 一 个 直 观 的 推 断 就 是 在 每 一 点 选 择 该 点 目 标 函 数 下 降 最 快 的 方 向 p i = F (x i ) (2.38) 作 为 寻 找 方 向, 这 一 方 法 称 为 最 陡 下 降 法. 直 观 的 想 象 可 能 会 认 为 最 陡 下 降 方 向 是 一 种 理 想 的 寻 找 方 向, 但 事 实 并 非 如 此, 最 陡 下 降 方 向 只 表 示 在 x i 附 近 的 下 降 性 质, 而 对 整 个 最 优 化 过 程 来 说, 我 们 寻 找 的 是 全 局 的 最 优 方 向. 作 为 一 个 例 子, 考 虑 如 下 函 数 的 最 优 化 : F = x y 2 (2.39) 显 然, 这 一 函 数 的 极 小 点 是 (0, 0). 函 数 (2.39) 的 负 梯 度 为 g = (2x, 20y) T (2.40) 如 果 从 点 (2, 1) 出 发, 利 用 最 陡 下 降 法 计 算 所 得 逐 x 如 下 : (2, 1), ( , ), ( , ), ( , ), ( , ), ( , ), ( , ), ( , ), ( , ), ( , ), ( , ), ( , ), ( , ), ( , ), ( , ), ( ). 显 然, 收 敛 是 很 慢 的. 最 陡 下 降 方 向 并 不 是 理 想 的 下 降 方 向, 为 了 改 进 最 优 化 的 收 敛 速 度, 需 要 寻 求 其 他 方 向. 为 了 判 断 一 个 寻 找 方 向 的 好 坏, 应 该 有 一 个 标 准, 注 意 到 一 般 函 数 在 最 小 点 附 近 近 似 于 二 次 函 数, 因 此 如 果 一 个 方 法 对 二 次 函 数 比 较 有 效, 则 可 望 其 对 一 般 函 数 也 比 较 有 效, 如 果 一 个 方 法 对 二 次 函 数 的 效 果 都 不 好 的 话, 很 难 期 望 它 对 一 般 函 数 会 有 好 的 效 果.

30 30 第 二 章 物 理 学 中 的 常 用 数 值 方 法 计 算 实 践 也 证 实 了 这 一 论 断. 下 面, 我 们 就 从 二 次 函 数 寻 求 一 种 比 较 有 效 的 下 降 方 向. 为 此, 我 们 先 给 出 几 个 数 学 定 理 ; 定 理 一 : 若 N 维 欧 几 里 德 空 间 中 的 向 量 q 与 N 个 线 性 独 立 向 量 p 1, p 2,, p N 都 正 交, 则 向 量 q 必 为 0. 这 一 定 理 的 结 论 是 显 然 的. 定 义 : A 共 轭 向 量. 设 A 是 一 个 N N 的 对 称 正 定 矩 阵, p, q 是 两 个 N 维 向 量, 若, p T Aq = 0 (2.41) 则 称 向 量 p 和 向 量 q 互 为 A 共 轭 或 互 为 A 正 交. 定 理 二 : 若 A 为 N N 对 称 正 定 矩 阵, p 1, p 2,, p N 为 A 共 轭 的 N 维 非 零 向 量, 则 此 向 量 组 必 为 线 性 独 立. 证 明 : 设 向 量 组 p 1, p 2,, p N 之 间 存 在 线 性 关 系 α 1 p 1 + α 2 p α N P N = 0 对 i = 1, 2,, N, 用 p T i A 左 乘 上 式 得 α i p T i Ap i = 0 因 为 p i 0, A 正 定, 所 以 p T i Ap i > 0 从 而 必 有 α i = 0, i = 1, 2,, N 因 此, 向 量 组 p 1, p 2,, p N 线 性 独 立. 有 了 以 上 数 学 准 备, 我 们 来 考 虑 二 次 函 数 的 最 优 化 问 题, 我 们 的 目 的 是 找 一 种 下 降 方 向, 能 够 在 有 限 步 达 到 二 次 函 数 的 极 小 点. 考 虑 二 次 N 元 函 数 F (x) = a + b T x xt Ax (2.42) 其 中 A 为 一 正 定 矩 阵. 函 数 (2.44) 的 梯 度 为 G(x) = b + Ax (2.43) 记 x i 点 的 梯 度 为 g i = G(x i ), 从 x 0 点 出 发, 对 于 第 一 个 寻 找 方 向, 除 梯 度 外 没 有 别 的 信 息, 我 们 就 取 p 0 = g 0, 沿 此 方 向 找 到 极 小 点 x 1 并 可 求 得 这 一 点 的 梯 度 g 1, 因 为 x 1 是 沿 p 0 方 向 的 极 小 点, 这 一 点 的 梯 度 方 向 一 定 与 p 0 垂 直,g 1 p 0, 从 而 g 1 g 0. 现 在 有 了 两 个

31 2.2 方 程 的 求 根 和 最 优 化 方 法 31 方 向,g 1 和 g 0, 最 陡 下 降 法 选 取 新 的 方 向 为 g 1. 我 们 保 留 原 来 方 向 的 一 点 信 息, 选 新 的 寻 找 方 向 为 这 两 个 方 向 的 组 合 p 1 = g 1 α 0 g 0 (2.44) 并 要 求 p 1 与 p 0 为 A 共 轭, 则 应 有 ( g 1 α 0 g 0 ) T Ap 0 = 0 (2.45) 由 于 x 1 = x 0 + t 0 p 0 则 g 1 g 0 = A(x 1 x 0 ) = t 0 Ap 0 一 般, 若 第 i 步 的 寻 找 方 向 为 p i, 则 同 理 可 证 g i+1 g i = A(x i+1 x i ) = t i Ap i (2.46) 因 此, 方 程 (2.45) 成 为 ( g 1 α 0 g 0 ) T (g 1 g 0 ) = 0 (2.47) 由 此 解 得 这 里 假 定 g 0 0, 否 则, 若 g 0 = 0, 则 x 0 就 是 所 求 的 极 值 点. α 0 = gt 1 g 1 g T 0 g 0 (2.48) 把 α 0 代 入 (2.44), 便 得 到 p 1. 沿 p 1 求 出 极 小 点 x 2 并 算 出 g 2. 由 于 p 0 和 p 1 为 A 共 轭, 因 此 由 (2.46) g T 0 (g 2 g 1 ) = 0 考 虑 到 g 0 与 g 1 正 交, g T 0 g 1 = 0, 由 上 式 可 知 g 0 与 g 2 也 正 交. 因 x 2 是 沿 p 1 方 向 的 极 小 点, 故 g 2 p 1, 也 就 是 g T 2 ( g 1 α 0 g 0 ) = 0 即 g 2 g 1. g 0, g 1 和 g 2 构 成 一 个 正 交 矢 量 组, 我 们 可 以 在 这 个 正 交 组 构 成 的 三 维 空 间 寻 求 与 p 0, p 1 均 为 A 共 轭 的 寻 找 方 向 p 2. 令 p 2 = g 2 α 1 g 1 β 0 g 0 (2.49) 要 求 p 2 分 别 与 p 0 和 p 1 为 A 共 轭, 也 就 是 要 求 如 下 方 程 成 立 ( g 2 α 1 g 1 β 0 g 0 ) T (g 1 g 0 ) = 0 (2.50) ( g 2 α 1 g 1 β 0 g 0 ) T (g 2 g 1 ) = 0 (2.51) (2.52)

32 32 第 二 章 物 理 学 中 的 常 用 数 值 方 法 由 此 解 出 于 是 α 1 = gt 2 g 2 g T 1 g 1, β 0 = α 0 α 1. (2.53) p 2 = g 2 + α 1 p 1 (2.54) 现 假 定 已 求 得 p 0, p 1,, p j, 其 中 p 0 = g 0 p i = g i + α i 1 p i 1 α i 1 = g T i g i /g T i 1g i 1 i = 1, 2,, j 这 里 假 定 g 0, g 1,, g j 都 不 为 零 且 构 成 一 个 正 交 系. 沿 p j 求 出 x j+1 并 算 出 g j+1 = G(x j+1 ), 设 g j+1 0, 则 仿 前 可 证, g 0, g 1,, g j+1 也 构 成 正 交 系. 令 p j+1 = g j+1 β 0 g 0 β 1 g 1 β j 1 g j 1 α j g j (2.55) 并 要 求 p j+1 与 p 0, p 1,, p j 为 A 共 轭, 即 ( g j+1 β 0 g 0 β 1 g 1 β j 1 g j 1 α j g j ) T (g j+1 g j ) = 0 (2.56) 得 到 α j = g T j+1g j+1 /g T j g j, β j 1 = α j α j 1, β 1 = α j α j 1 α 1, β 0 = α j α j 1 α 1 α 0. 且 β 0 g 0 β 1 g 1 β j 1 g j 1 α j g j = α j α 1 (α 0 g 0 + g 1 ) β 2 g 2 α j g j = α j (α 1 p 1 g 2 ) α j g j = = α j p j

33 2.2 方 程 的 求 根 和 最 优 化 方 法 33 即 p j+1 = g j+1 + α j p j, α j = g T j+1g j+1 /g T j g j (2.57) 根 据 归 纳 法, 在 (2.57) 中 令 j = 0, 1, 2,, 所 得 方 向 为 A 共 轭 方 向. 上 述 方 法 最 多 只 要 N 步 就 可 找 到 函 数 (2.44) 的 极 小 点, 这 是 因 为 所 有 的 g j 互 相 正 交, 即 g T N g j = 0 对 所 有 j = 0, 1, 2,, N 1 都 成 立, 而 我 们 的 问 题 在 N 维 空 间, 故 必 有 g N = 0, 亦 即 x N 为 函 数 (2.44) 的 极 小 点. 这 就 是 所 谓 共 轭 梯 度 方 法, 理 论 上, 对 于 二 次 函 数, 只 要 N 次 迭 代 就 可 以 达 到 极 小 点, 但 实 际 计 算 时, 由 于 舍 入 误 差 的 影 响, 一 般 N 次 迭 代 并 不 能 达 到 极 小 点. 对 于 一 般 的 目 标 函 数, 有 限 次 (N 次 ) 的 迭 代 也 一 般 不 能 达 到 极 小 点, 因 此 实 际 计 算 时, 从 一 初 始 点 出 发, 迭 代 N 次 后, 以 所 得 结 果 为 出 发 点, 重 新 开 始 计 算, 直 到 收 敛 为 止. 共 轭 梯 度 方 法 由 于 占 用 内 存 小, 方 法 简 单, 编 程 方 便, 目 前 在 物 理 问 题 中 应 用 比 较 多. 作 为 例 子, 我 们 再 考 虑 前 面 的 例 题, 对 函 数 x y 2, 取 x 0 = (2, 1) T, 则 g 0 = (4, 20) T, 取 p 0 = g 0, 求 得 x 1 = ( , ) T, g 1 = ( , ) T, α 0 = , 由 此 可 构 造 出 p 1 = ( , ) T, 沿 p 1 方 向 求 得 x 2 = ( , ) T, 即 已 经 达 到 极 小 点. 习 题 : 编 制 共 轭 梯 度 法 的 程 序 并 求 下 面 函 数 的 极 小 值 点 及 极 小 值, F (x 1,, x n ) = 1 + n [ 100(xi x 2 i 1) 2 + (1 x i 1 ) 2] (2.58) i=2 取 n = 10, 选 择 初 始 点 为 x i = i/10, i = 1, 2,, 约 束 最 优 化 前 面 几 节 所 讲 的 是 无 约 束 最 优 化 问 题, 在 实 际 物 理 问 题 中, 还 经 常 会 遇 到 有 约 束 最 优 化 问 题, 其 数 学 提 法 是, 给 定 一 目 标 函 数 F (x 1, x 2,, x n ), 在 一 定 的 约 束 下 求 解 其 极 小 值, 约 束 一 般 写 成 一 组 函 数. 即 { F (x1, x 2,, x n ) = Min G i (x 1, x 2,, x n ) = 0 (i = 1, 2,, m) (2.59) 式 中 m < n, m 个 方 程 G i (x 1, x 2,, x n ) 代 表 m 个 约 束. 解 决 约 束 优 化 问 题 的 通 常 做 法 是 构 造 一 个 新 的 函 数, 把 约 束 优 化 问 题 改 为 一 无 约 束 最 优 化 问 题, 然 后 利 用 前 面 的 方 法 求 解. 通 常 的 做 法 是, 定 义 F = F + m λ i G 2 i (2.60) i=1

34 34 第 二 章 物 理 学 中 的 常 用 数 值 方 法 对 于 合 适 选 择 的 λ i, (i = 1, 2,, m), F 的 极 小 值 就 代 表 了 由 方 程 (2.59) 定 义 的 约 束 最 优 化 问 题. 例 题 在 参 数 x 1, x 2,, x n 的 一 个 限 制 区 域 D 内 求 解 F 的 极 小 值. 为 此, 我 们 把 求 解 区 域 扩 展 到 整 个 区 域, 如 果 F 在 D 外 无 定 义, 则 应 扩 充 其 定 义. 同 时 定 义 函 数 G 为 { 0 {x} D G = i λ ix 2 i {x} D (2.61) 则 通 过 计 算 F + G 的 极 小 值 就 可 得 到 问 题 的 解. 在 使 用 这 种 方 法 时 应 注 意 如 果 在 D 外 F 会 变 小 时, G 的 选 择 应 使 得 F + G 在 D 外 总 是 较 快 增 加 的. 式 (2.61) 给 出 的 仅 仅 是 一 种 可 能 的 选 择, 实 际 上 G 的 选 择 可 以 有 无 穷 多 种 最 小 二 乘 法 及 曲 线 拟 合 在 物 理 实 验 中, 常 常 要 测 量 一 对 物 理 量 之 间 的 关 系, 在 有 些 情 况 下, 一 对 物 理 量 之 间 的 函 数 关 系 是 已 知 的, 但 包 含 若 干 个 未 知 参 数, 若 用 x 和 y 代 表 这 一 对 物 理 量, 则 有 y = f(x; a 1, a 2,, a m ) (2.62) 式 中 a 1, a 2,, a m 为 一 组 未 知 参 量, f 的 函 数 形 式 是 已 知 的, 在 这 种 情 况 下, 曲 线 拟 合 的 任 务, 就 是 根 据 (x, y) 的 N 个 观 测 点 (x 1, y 1 ), (x 2, y 2 ),, (x N, y N ) 来 确 定 参 数 a i, (i = 1, 2,, m), 如 果 测 量 中 没 有 误 差, 则 每 一 组 (x i, y i ) 都 应 准 确 地 落 在 理 论 曲 线 上, 即 应 有 y i = f(x i ; a 1, a 2,, a m ) (i = 1, 2,, N) (2.63) 为 了 确 定 诸 a i 的 值, 只 要 选 出 (2.63) 中 的 m 个 方 程 求 解 即 可, 但 实 际 测 量 总 是 有 误 差 的, (x i, y i ) 并 不 准 确 落 在 理 论 曲 线 上, 在 N > m 的 情 况 下, 方 程 组 (2.63) 成 为 矛 盾 方 程 组, 不 能 用 解 方 程 的 方 法 求 参 数 a 1, a 2,, a m, 只 能 用 曲 线 拟 合 的 方 法 根 据 带 有 误 差 的 观 测 值 点 求 得 理 论 曲 线 参 数 的 估 计 值. 在 很 多 实 际 问 题 中, x 和 y 这 两 个 物 理 量 中 总 有 一 个 量 的 测 量 精 度 比 另 一 个 高 得 多, 其 测 量 误 差 可 以 忽 略. 把 可 以 忽 略 测 量 误 差 的 量 选 作 自 变 量 x, 其 观 测 值 可 看 作 是 准 确 的, 对 于 每 个 x i, 对 应 的 y 的 观 测 值 y i 是 一 个 随 机 变 量 ( 见 第 五 章 ), 其 误 差 可 由 其 方 差 来 估 计. 物 理 中 还 经 常 遇 到 另 外 一 种 曲 线 拟 合 的 任 务 : 物 理 量 y 和 x 间 的 函 数 形 式 未 知, 需 要 从 观 测 点 求 出 y 和 x 的 函 数 关 系 的 一 个 经 验 公 式. 一 般 的 作 法 是 选 取 某 一 函 数 族 {

35 2.2 方 程 的 求 根 和 最 优 化 方 法 35 ϕ i (x) } 的 一 个 线 性 组 合, y = i a i ϕ i (x) (2.64) 然 后 用 曲 线 拟 合 方 法 求 出 系 数 a i. 函 数 族 { ϕ i (x) } 通 常 从 物 理 考 虑 来 选 取, 一 种 常 用 的 方 法 是 选 取 { 1, x, x 2,, x m 1 } 用 于 拟 合 m 个 参 数 ; 也 可 选 取 三 角 函 数 { 1, sin x, os x, sin 2x, os 2x,, sin 2mx, os 2mx} 用 于 拟 合 2m + 1 个 参 数. 这 种 方 法 称 为 线 性 拟 合 方 法, 即 y 线 性 地 依 赖 于 参 数 a i, 对 于 有 些 物 理 问 题, y 和 x 的 函 数 关 系 中 非 线 性 地 依 赖 于 参 数 a i, 例 如 y = a 1 e a2x + a 3 e a4x sin(a 5 x + a 6 ) 就 是 一 个 具 有 6 个 参 数 的 非 线 性 拟 合 问 题. 曲 线 拟 合 通 常 采 用 最 小 二 乘 法, 如 果 x 的 测 量 误 差 可 以 忽 略, 观 测 点 i 的 残 差 δ i 定 义 为 在 该 点 y 的 观 测 值 与 理 论 值 之 差, 即 δ i = y i f(x i ; a 1, a 2,, a m ) (2.65) 若 观 测 值 y i 的 方 差 为 σi 2, 则 定 义 其 权 重 因 子 为 记 χ 2 为 观 测 点 残 差 的 加 权 平 方 和, χ 2 = N w i δi 2 = i=1 w i = 1 σ 2 i N w i (y i f(x i ; a 1, a 2,, a m )) 2 (2.66) i=1 曲 线 拟 合 的 最 小 二 乘 法 就 是 以 χ 2 为 目 标 函 数 的 优 化 问 题, 即 寻 找 一 组 参 数 a i, (i = 1, 2,, m) 使 χ 2 为 最 小. 在 观 测 点 的 方 差 未 知 的 情 况 下, 可 以 取 w i = 1, 也 可 以 在 每 一 观 测 点 多 测 量 几 次, 以 其 平 均 值 y i 作 为 该 点 的 观 测 值, 而 以 y 2 i y i 2 作 为 其 方 差. 在 曲 线 拟 合 问 题 中, 最 困 难 的 是 选 取 合 适 的 经 验 公 式, 这 需 要 结 合 物 理 考 虑 和 对 测 量 数 据 的 观 察, 也 依 赖 于 对 各 种 函 数 及 其 图 像 的 明 确 概 念. 例 题 在 某 一 能 谱 测 量 中, 得 到 下 表 所 示 数 据, 试 给 出 经 验 公 式 并 作 出 曲 线 拟 合. x y x y x y 把 表 中 数 据 点 到 坐 标 纸 上 ( 见 图 2.1 中 圈 点 ), 可 以 看 出 数 据 有 二 个 峰, 为 此, 我 们 试 用 Lorentz 形 曲 线 去 拟 合, 即 取 经 验 公 式 为 y = a 1 (x a 2 ) 2 + a 2 3 a 4 + (x a 5 ) 2 + a 2 6

36 36 第 二 章 物 理 学 中 的 常 用 数 值 方 法 图 2.1: 用 上 式 和 上 表 中 的 数 据 计 算 χ 2, 并 用 共 轭 梯 度 方 法 优 化 χ 2, 求 得 各 参 数 值 为 a 1 = a 2 = a 3 = a 4 = a 5 = a 6 = 在 图 2.1 中, 我 们 用 实 线 画 出 了 拟 合 的 结 果, 虚 线 给 出 了 二 个 峰 的 图 像. 2.3 函 数 的 求 值 在 数 值 计 算 中, 经 常 要 求 一 些 函 数 的 值, 一 些 常 用 的 初 等 函 数 如 三 角 函 数, e 指 数 函 数, 双 曲 函 数, 对 数 函 数 等, 均 已 写 入 计 算 机 语 言 的 定 义 做 为 内 部 函 数, 在 使 用 时 只 需 调 用 即 可, 而 其 它 的 函 数 和 表 达 式 则 要 自 己 计 算, 本 章 将 给 出 一 些 在 物 理 计 算 中 经 常 遇 到 的 函 数 及 表 达 式 的 计 算 方 法 和 部 分 程 序. 一 般 说 来, 计 算 机 的 内 部 函 数 都 是 由 专 家 设 计 的, 其 效 率 和 精 度 都 很 好, 因 此, 如 果 计 算 机 语 言 已 经 提 供 了 所 需 的 函 数, 则 应 尽 量 使 用 这 些 函 数, 下 面 提 供 的 方 法 仅 仅 作 为 计 算 机 系 统 不 提 供 这 些 函 数 时 的 补 充 多 项 式 的 求 值 和 级 数 求 和 所 谓 多 项 式 是 指 由 下 式 给 出 的 表 达 式 : p(x) = a 0 + a 1 x + a 2 x 2 + a 3 x a n x n 为 了 计 算 p(x) 在 给 定 x 时 的 数 值, 我 们 用 (n), n = 0, 1, 2,, n 来 代 表 多 项 式 的 系 数, 则 可 以 写 出 计 算 语 句. 以 n = 4 为 例, 直 接 用 FORTRAN 语 言 写 出 的 算 式 为 :

37 2.3 函 数 的 求 值 37 p=(0)+(1)*x+(2)*x**2+(3)*x**3+(4)*x**4 这 显 然 不 是 一 个 好 的 算 法 ( 为 什 么?). 较 好 的 算 法 是 下 面 二 者 之 一 : p=(0)+x*((1)+x*((2)+ x*((3)+x*(4)))) p=((((4)*x+(3))*x+(2))*x+(1))*x+(0) 当 n 较 大 时, 可 采 用 下 面 的 语 句 段 : P=C(N) DO 11 J=N-1,1,-1 P=P*X+C(J) 11 CONTINUE 在 科 研 中, 还 常 常 碰 到 级 数 求 和 的 问 题, 例 如 计 算 S = i=0 u i 我 们 当 然 不 可 能 计 算 到 无 穷 多 项, 所 以 上 述 级 数 必 需 在 某 处 截 断, 如 果 级 数 收 敛 很 慢, 就 要 计 算 很 多 项, 一 方 面 占 用 大 量 的 CPU 时 间, 更 重 要 的 是 当 计 算 的 项 数 非 常 多 时, 舍 入 误 差 的 积 累 有 可 能 完 全 破 坏 整 个 计 算. 为 此 必 须 使 用 加 速 收 敛 的 技 巧. Aitken 算 法 是 常 用 的 加 速 技 巧 之 一. 它 是 对 相 邻 的 三 个 部 分 和 进 行 处 理, 得 到 一 个 更 好 的 近 似. 记 部 分 和 为 : S n = n u i (2.67) i=0 则 由 S n 1, S n, S n+1 可 计 算 出 S n = S n+1 (S n+1 S n ) 2 S n+1 2S n + S n 1 (2.68) 对 n = 1, 2, 3,, 先 求 出 S n, 同 时 算 出 S n, 一 般 S n 更 快 地 收 敛 于 S. 在 实 际 使 用 时, (2.68) 最 好 用 写 出 的 式 子 计 算, 一 些 代 数 上 等 价 的 式 子 可 能 导 致 较 大 的 舍 入 误 差. Euler 变 换 对 于 由 下 式 定 义 的 交 错 级 数 ( 1) s u s (2.69) 这 里 u s > 0. Euler 变 换 是 一 个 强 有 力 的 算 法, 它 由 下 式 给 出 ( 取 n 为 偶 数 ) 其 中 : s=0 ( 1) s u s = u 0 u 1 + u 2 u n 1 + s=0 s=0 ( 1) s 2 s+1 ( s u n ) (2.70) u n u n+1 u n 2 u n ( u n ) = u n+2 2u n+1 + u n (2.71)

38 38 第 二 章 物 理 学 中 的 常 用 数 值 方 法 为 向 前 差 分. 式 (2.70) 可 推 导 如 下 : 定 义 位 移 算 符 E 为 : u n+1 Eu n 由 于 (1 + )u n = u n + u n+1 u n = u n+1 = Eu n, 因 此 有 算 符 等 式 E = 1 + 而 u 0 u 1 + u 2 u 3 + = (1 E + E 2 E 3 + )u 0 = E u 0 = 1 1 u = ( 1) s 2 s+1 s u 0 2 s= 连 分 式 的 求 值 所 谓 连 分 式, 是 指 由 下 式 给 出 的 表 达 式, f = b 0 + b 1 + b 2 + a 1 b 3 + a 2 a 3 a 4 b 4 + a 5 b 5 + 数 学 家 们 对 它 做 过 很 多 研 究, 这 是 一 个 非 常 有 趣 的 函 数 表 示 方 法. 除 了 上 面 的 表 达 式 外, 人 们 还 经 常 使 用 另 一 种 记 法. 这 里, 系 数 a n 和 b n 也 可 以 是 x 的 函 数. a 2 a 3 a 4 a 5 f = b 0 + a 1 b 1 + b 2 + b 3 + b 4 + b 5 + 由 于 连 分 式 只 能 从 右 向 左 求 值, 因 此 在 具 体 计 算 之 前, 很 难 判 断 一 个 无 限 的 连 分 式 应 该 在 何 处 截 断, 一 种 直 接 的 方 法 就 是 猜 一 个 截 断 点, 求 出 一 个 值, 再 猜 一 个 更 大 一 点 的 截 断 点 计 算, 通 过 比 较 不 同 截 断 处 的 结 果 来 判 断 是 否 收 敛. 这 显 然 不 是 一 个 好 的 方 法. 另 一 种 在 实 践 中 常 用 的 方 法 是 用 一 个 分 式 有 理 函 数 来 近 似 连 分 式, 用 f n 表 示 计 算 到 a n 和 b n 的 连 分 式 的 值, 则 : 这 里 A n, B n 由 下 面 的 递 推 公 式 给 出 : A 1 1 B 1 0 f n = A n B n (2.72) A 0 b 0 B 0 1 (2.73) A j = b j A j 1 + a j A j 2 B j = b j B j 1 + a j B j 2 j = 1, 2,, n

39 2.4 常 微 分 方 程 的 数 值 解 法 39 练 习 试 证 明 递 推 关 系 (2.73) 并 写 出 计 算 机 程 序 2.4 常 微 分 方 程 的 数 值 解 法 常 微 分 方 程 在 物 理 学 中 具 有 重 要 意 义, 牛 顿 第 二 定 律 就 表 示 为 一 组 二 阶 常 微 分 方 程, 第 一 章 中 提 到 的 著 名 的 三 体 问 题 就 是 由 九 个 二 阶 常 微 分 方 程 组 成 的 方 程 组 ; 在 量 子 力 学 中 和 电 动 力 学 中, 由 薛 定 谔 方 程 及 Maxwell 方 程 分 离 变 量 得 到 的 所 谓 Sturm Liouville 本 征 值 问 题 也 是 二 阶 常 微 分 方 程 的 求 解 问 题. 在 这 一 章 中, 我 们 将 给 出 几 种 数 值 求 解 常 微 分 方 程 的 常 用 方 法 Euler 方 法 对 于 常 微 分 方 程 的 初 值 问 题 dy dx = f(x, y) y(0) = η (2.74) 数 值 求 解 的 第 一 步 是 把 方 程 (2.74) 变 为 一 个 差 分 方 程, 为 此, 让 x 取 一 组 分 立 值 x n, n = 0, 1, 2,, x 0 = 0, x n+1 = x n + h, 这 里 h 称 为 积 分 的 步 长, 把 y(x) 在 x n 处 展 开, 有 y(x n+1 ) = y(x n ) + h dy dx + O(h 2 ) (2.75) x=xn 利 用 方 程 (2.74), 便 可 得 到 y(x n+1 ) = y(x n ) + hf(x n, y n ) + O(h 2 ) (2.76) 这 就 是 Euler 差 分 格 式, 从 x = 0 出 发, 利 用 这 一 方 程 一 步 一 步 地 向 前 计 算, 就 可 以 得 到 任 一 x n 处 的 函 数 值 y n y(x n ). Euler 方 法 是 一 个 显 式 方 法, 即 第 n + 1 步 的 y n+1 只 依 赖 于 x n, y n 的 值, 而 且 每 一 步 只 需 计 算 一 次 f(x n, y n ) 的 值, 这 是 Euler 方 法 的 优 点, 但 另 一 方 面, Euler 方 法 也 存 在 严 重 的 不 足, 因 为 这 一 方 法 是 一 次 的 ( 每 一 步 的 误 差 项 为 h 2 的 数 量 级, y(x N ) Nh 2 h), 所 以 h 必 须 取 得 足 够 小, 以 保 证 有 满 意 的 精 度, 但 步 长 太 小 时, 积 分 到 一 预 定 的 x 则 要 做 很 多 步 的 计 算, 舍 入 误 差 的 积 累 将 会 使 计 算 结 果 严 重 失 真. 由 于 这 一 原 因, 在 实 际 计 算 中 几 乎 不 使 用 Euler 方 法, 而 常 常 利 用 其 推 导 上 的 简 单 和 概 念 上 非 常 清 楚 而 用 于 教 学.

第 9 卷 江 南 大 学 学 报 人 文 社 会 科 学 版 Z 第 2 期 掌握 是指在 表 层 知 识 教 学 过 程 中 学 生 对 表 层 知 识 的 掌 想 方法有所悟 有所体会 5 数学思想 方法教学是循环往 握 学生掌握 了 一 定 量 的 数 学 表 层 知 识 是 学 生 能 够

第 9 卷 江 南 大 学 学 报 人 文 社 会 科 学 版 Z 第 2 期 掌握 是指在 表 层 知 识 教 学 过 程 中 学 生 对 表 层 知 识 的 掌 想 方法有所悟 有所体会 5 数学思想 方法教学是循环往 握 学生掌握 了 一 定 量 的 数 学 表 层 知 识 是 学 生 能 够 江南大学学报 人文社会科学版 第9卷 Z第2期 2010 年 12 月 犑 狅狌 狉狀犪 犾狅 犳犑 犻 犪狀犵狀犪狀犝狀 犻 狏 犲 狉 狊 犻 狋 犎狌犿犪狀 犻 狋 犻 犲 狊牔 犛狅 犮 犻 犪 犾犛 犮 犻 犲狀犮 犲 狊 狔 Vo l 9 Z2 Dec 2010 数学教学 初中数学教学研究 雍兴智 南部县四龙乡小学 四川 南部 637300 不论我们选教什么 学 科 务 必 使 学 生

More information

#!!

#!! 浅 谈 日 本 的 文 化 外 交 吴 咏 梅 近 年 来 日 本 实 行 借 助 动 漫 影 视 等 软 实 力 促 进 与 海 外 的 相 互 理 解 和 友 好 输 出 其 价 值 观 和 提 高 国 家 形 象 的 新 型 文 化 外 交 战 略 为 扩 大 其 对 国 际 事 务 的 影 响 力 提 高 国 际 地 位 实 现 政 治 大 国 的 目 标 而 做 出 了 种 种 努 力 日

More information

證 明 : 令 φ(x f(x, ydy, 則 φ(x + x φ(x x f x (ξ, ydy f x (ξ, y f x (x, y dy f x (x, ydy f(x + x, y f(x, y d dy f x (x, ydy x f x (x, ydy, ξ ξ(y 介 於 x, x

證 明 : 令 φ(x f(x, ydy, 則 φ(x + x φ(x x f x (ξ, ydy f x (ξ, y f x (x, y dy f x (x, ydy f(x + x, y f(x, y d dy f x (x, ydy x f x (x, ydy, ξ ξ(y 介 於 x, x 微 分 與 積 分 的 交 換 積 分 設 f 在 [a, b] [, d] 上 連 續, 問 d dx f(x, y? f(x, ydy x 首 先 (1 式 兩 邊 必 須 有 意 義 f(x, ydy 必 須 對 x 可 導 若 f 及 x f(x, ydy 積 分 必 須 存 在 x f 在 [a, b] [, d] 上 連 續, 則 ( 及 (3 式 成 立, 下 面 的 定 理 告 訴

More information

( β ) () () R () R ) β ( ) ( ( β ) () R [ ] C( ) f ( ) g( ) ( f ( ) g( )) f ( ) g( ) d () () C( ) R [ ] R[ ] () 4 H ξ ( ) < H (Hlert) ) ( k β ) ( kβ )

( β ) () () R () R ) β ( ) ( ( β ) () R [ ] C( ) f ( ) g( ) ( f ( ) g( )) f ( ) g( ) d () () C( ) R [ ] R[ ] () 4 H ξ ( ) < H (Hlert) ) ( k β ) ( kβ ) 4 5 : V R V ( β ) ) ( β ) ( β ) ; ) ( k β ) k( β ) ; : ) ( β γ ) ( γ ) ( β γ ) ; 4) ( ) ( ) β γ V k V R ) β ( ) ( ( β ) () () R () R ) β ( ) ( ( β ) () R [ ] C( ) f ( ) g( ) ( f ( ) g( )) f ( ) g( ) d

More information

總目186-運輸署

總目186-運輸署 管 制 人 員 : 運 輸 署 署 長 會 交 代 本 總 目 下 的 開 支 二 零 一 六 至 一 七 年 度 預 算... 28.585 億 元 二 零 一 六 至 一 七 年 度 的 編 制 上 限 ( 按 薪 級 中 點 估 計 的 年 薪 值 ) 相 等 於 由 二 零 一 六 年 三 月 三 十 一 日 預 算 設 有 的 1 5 3 6 個 非 首 長 級 職 位, 增 至 二 零

More information

Microsoft Word - BD-QH2004.doc

Microsoft Word - BD-QH2004.doc 从 线 性 方 程 谈 起 ( 清 华 大 学 数 学 科 学 系 章 梅 荣 教 授 ) 烟 台 大 学 与 北 京 大 学 和 清 华 大 学 有 着 密 切 的 联 系, 我 在 北 大 读 书 时 的 不 少 同 学 朋 友 还 在 烟 台 大 学 工 作 我 应 邀 来 烟 台 大 学 作 这 个 演 讲, 受 到 了 朋 友 们 的 热 情 接 待, 我 感 到 非 常 高 兴, 也 深

More information

標準 BIG 中文字型碼表 A 0 9 B C D E F 一 乙 丁 七 乃 九 了 二 人 儿 入 八 几 刀 刁 力 匕 十 卜 又 三 下 丈 上 丫 丸 凡 久 么 也 乞 于 亡 兀 刃 勺 千 叉 口 土 士 夕 大 女 子 孑 孓 寸 小 尢 尸 山 川 工 己 已 巳 巾 干 廾

標準 BIG 中文字型碼表 A 0 9 B C D E F 一 乙 丁 七 乃 九 了 二 人 儿 入 八 几 刀 刁 力 匕 十 卜 又 三 下 丈 上 丫 丸 凡 久 么 也 乞 于 亡 兀 刃 勺 千 叉 口 土 士 夕 大 女 子 孑 孓 寸 小 尢 尸 山 川 工 己 已 巳 巾 干 廾 標準 BIG 中文字型碼表 A 0 9 B C D E F B C D ± E F A 0 9 B C D E F 兙 兛 兞 兝 兡 兣 嗧 瓩 糎 0 B 9 Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ Ⅵ Ⅶ C Ⅷ Ⅸ Ⅹ 〡 〢 〣 〤 〥 〦 〧 〨 〩 十 卄 卅 D B C D E F G H I J K L M N O P Q E R S T U V W X Y Z a b c d e f g F h i

More information

! # % &# % %( ) #&#! % ( ) # +.! / 0 1 2 3 4 % )! & + 5 6 % #!& # & 7! 8 #&% 9! : ; +>?1 & 9! & % ( 2 2 ; & #! #! % & % & 2! %!+ &/ Α Β # % & # & #% % ( = #& #! #! & #

More information

ii

ii i 概 率 统 计 讲 义 原 著 : 何 书 元 课 件 制 作 : 李 东 风 2015 年 秋 季 学 期 ii 目 录 第 一 章 古 典 概 型 和 概 率 空 间 3 1.1 试 验 与 事 件............................ 3 1.2 古 典 概 型 与 几 何 概 型....................... 7 1.2.1 古 典 概 型.........................

More information

素 4 在上述学 者 观 点 的 基 础 上本 文 认 为 员 工 需 要 一 理论与假设 同时具备创新意 识 创 新 能 力 创 新 动 机 和 创 新 机 会 才 能产生创新行为 创新 意 识 是 指 员 工 能 通 过 对 组 织 环 境 一 组织的创新战略与员工的创新行为 的解读认识到创新

素 4 在上述学 者 观 点 的 基 础 上本 文 认 为 员 工 需 要 一 理论与假设 同时具备创新意 识 创 新 能 力 创 新 动 机 和 创 新 机 会 才 能产生创新行为 创新 意 识 是 指 员 工 能 通 过 对 组 织 环 境 一 组织的创新战略与员工的创新行为 的解读认识到创新 中国人民大学学报 5 年第 5 期 N 55 OURNAL OF RENMN UNVER TY OF CHNA 组织创新战略如何转化为 员工创新行为 以中关村T 行业为例 苏中兴 张雨婷 曾湘泉 摘要 随 着 市 场 竞 争 加 剧 和 劳 动 力 成 本 快 速 上 升 越 来 越 多 的 中 国 企 业 开 始 采 用 创 新 战 略 来 应 对挑战 组织特定的战略需要匹配的员工行为作为支撑因此激发员工的创新

More information

5-2微积分基本定理

5-2微积分基本定理 s v s v. T, ]. s v ; n lim v ξi i λ i F ', [ T s T s T T T T T v s T s T v F F T v v T T T s v ξ i i i [ T, T] . b Φ [, b] [, b] [, b] . [, b], y y Φ b o Φ ' [, b] o, bφ Φ ξ b + h + h Φ + h + + h Φ Φ +

More information

Microsoft Word - chead095.doc

Microsoft Word - chead095.doc 管 制 人 员 : 康 乐 及 文 化 事 务 署 署 长 会 交 代 本 总 目 下 的 开 支 二 O 一 O 至 一 一 年 度 预 算... 54.547 亿 元 二 O 一 O 至 一 一 年 度 的 编 制 上 限 ( 按 薪 级 中 点 估 计 的 年 薪 值 ) 相 等 于 由 二 O 一 O 年 三 月 三 十 一 日 预 算 设 有 的 8 049 个 非 首 长 级 职 位,

More information

d y dy P x Q x y 0. dx dx d d P x Q x C C 1y1 y dx dx d d P x Q x C 1y 1 dx dx d d P x Q x C y 0. dx dx d x 1dx F. ox1 dt dt d x1 1dx1 x 0 1 F 1 dt dt d x 1dx x 0 F dt dt d y 1dy y F 0 1 F1 y x1 x. dt

More information

untitled

untitled 5 55-% 8-8 8-5% - 7 7 U- n n lim n n n n lim n n n n un n n k k n k u n k k k n k ` n k u n k k n n n lim n n.7. 8 f f e f f f e. f f e f f f f f f f [ e ] f e e www.tsinghuatutor.om 5 79 755 f f f f [

More information

<4D F736F F D20CAFDD1A7D1A7D4BABACFB2A2B3F6C6ACCEC4BCFE2E646F63>

<4D F736F F D20CAFDD1A7D1A7D4BABACFB2A2B3F6C6ACCEC4BCFE2E646F63> 课 程 质 量 标 准 汇 编 ( 数 学 学 院 2010) 教 务 处 编 印 PDF 文 件 使 用 "pdffactory Pro" 试 用 版 本 创 建 www.fineprint.cn 目 录 高 等 代 数 方 法 选 讲 课 程 简 介 教 学 大 纲 考 核 大 纲 1 数 学 分 析 方 法 选 讲 课 程 简 介 教 学 大 纲 考 核 大 纲 7 数 学 分 析 课 程 简

More information

ABC Amber Text Converter

ABC Amber Text Converter 释 迦 牟 尼 佛 广 传 白 莲 花 论 全 知 麦 彭 仁 波 切 著 下 一 页 Page 1 一 功 德 品 上 一 页 首 页 下 一 页 一 功 德 品 总 体 说 来, 以 业 力 感 召 而 轮 回 于 三 界 中 之 众 生, 无 有 能 摆 脱 三 大 痛 苦 抛 开 束 缚 的 自 由, 因 而 也 得 不 到 任 何 究 竟 的 安 乐 自 在 正 因 为 他 们 从 无 始

More information

为 了 进 一 步 研 究 这 种 新 射 线 的 性 质, 搞 清 楚 这 个 不 速 之 客 的 真 实 身 份, 伦 琴 在 玻 璃 管 与 屏 幕 之 间 放 了 一 本 比 较 厚 的 书, 结 果 照 样 可 以 看 到 荧 光 随 后, 他 又 把 一 块 薄 木 板 放 在 了 书

为 了 进 一 步 研 究 这 种 新 射 线 的 性 质, 搞 清 楚 这 个 不 速 之 客 的 真 实 身 份, 伦 琴 在 玻 璃 管 与 屏 幕 之 间 放 了 一 本 比 较 厚 的 书, 结 果 照 样 可 以 看 到 荧 光 随 后, 他 又 把 一 块 薄 木 板 放 在 了 书 伦 琴 的 发 现 19 世 纪, 是 人 类 历 史 上 一 个 非 常 重 要 的 发 展 时 期 工 业 革 命 的 兴 起, 带 动 了 基 础 理 论 的 研 究 ; 反 过 来, 基 础 理 论 的 研 究 成 果, 又 进 一 步 推 动 了 工 业 的 发 展 特 别 是 19 世 纪 中 期, 随 着 电 力 的 广 泛 应 用, 人 们 对 于 生 产 过 程 中 出 现 的 放

More information

22 8529 專 用 或 主 要 用 於 第 8525 至 8528 節 所 屬 器 具 之 零 件 23 8308 用 於 衣 服 靴 鞋 帳 蓬 手 提 包 旅 行 用 品 或 其 他 已 製 作 品 之 卑 金 屬 搭 鈕 帶 搭 鈕 之 框 架 帶 扣 帶 扣 搭 鈕 眼 環 眼 及 其

22 8529 專 用 或 主 要 用 於 第 8525 至 8528 節 所 屬 器 具 之 零 件 23 8308 用 於 衣 服 靴 鞋 帳 蓬 手 提 包 旅 行 用 品 或 其 他 已 製 作 品 之 卑 金 屬 搭 鈕 帶 搭 鈕 之 框 架 帶 扣 帶 扣 搭 鈕 眼 環 眼 及 其 附 表 一 項 次 4 位 碼 稅 號 我 國 對 哥 倫 比 亞 出 口 符 合 哥 國 零 關 稅 項 目 清 單 ( 以 2014 年 海 關 稅 則 前 4 碼 排 名 前 300 項 產 品 為 比 對 基 礎 ) 中 文 名 稱 2013/01-2013/12 2014/01-2014/12 單 位 : 金 額 ( 美 元 ) 名 次 增 減 比 (%) 全 部 貨 品 出 口 總 額

More information

.. (,, ),(): ( (,,, (, (,,, (), ( ): (, ), ( ): (,,, (,,, (,,, ),( ): (,,, (, (,,, (,, (, ),(): ~ (, ~ (,, ~ (, ~ (,,, ), ( ), ( ): ( (,, ~ (,, ), ~ (

.. (,, ),(): ( (,,, (, (,,, (), ( ): (, ), ( ): (,,, (,,, (,,, ),( ): (,,, (, (,,, (,, (, ),(): ~ (, ~ (,, ~ (, ~ (,,, ), ( ), ( ): ( (,, ~ (,, ), ~ ( .. 中 国 热 带 医 学 杂 志 年 卷 主 题 索 引 (,,, ),( ): (, (,,, ( (,,, ~ (, ), (,,, (,, (,,, (, ),( ): (, (,, (,, (, (,,, ),( ): ~ (,,, (,, (,,, (,,, ),( ): ( ), ( ): - (, (, (, (,,, ),(): (,,, ), α- (,,, ), (, (,,,

More information

3.2 導 函 數 其 切 線 (tangent line) 為 通 過 P, 且 其 斜 率 為 m 的 直 線, 即 y = f(a) + m(x a) (3) 其 法 線 (normal line) 為 通 過 P 且 與 切 線 垂 直 的 直 線, 即 y = f(a) 1 (x a) m

3.2 導 函 數 其 切 線 (tangent line) 為 通 過 P, 且 其 斜 率 為 m 的 直 線, 即 y = f(a) + m(x a) (3) 其 法 線 (normal line) 為 通 過 P 且 與 切 線 垂 直 的 直 線, 即 y = f(a) 1 (x a) m 第 3 章 微 分 (Differentiation) 目 錄 3.1 切 線................................... 25 3.2 導 函 數.................................. 26 3.3 微 分 公 式................................. 28 3.4 連 鎖 律..................................

More information

产 科 麻 醉 指 南 推 荐 剖 宫 产 常 规 采 用 椎 管 内 麻 [ ] 醉, 全 身 麻 醉 (, ) 用 于 产 科 易 发 生 胎 儿 窒 息, 所 以 目 前 主 要 用 于 有 椎 管 内 麻 醉 禁 忌 证 的 产 妇 全 麻 产 科 中 新 生 儿 的 安 全 问 题 一 直

产 科 麻 醉 指 南 推 荐 剖 宫 产 常 规 采 用 椎 管 内 麻 [ ] 醉, 全 身 麻 醉 (, ) 用 于 产 科 易 发 生 胎 儿 窒 息, 所 以 目 前 主 要 用 于 有 椎 管 内 麻 醉 禁 忌 证 的 产 妇 全 麻 产 科 中 新 生 儿 的 安 全 问 题 一 直 中 国 病 理 生 理 杂 志, ( ): - [ 文 章 编 号 ] - ( ) - - 全 身 麻 醉 对 脐 血 乳 酸 以 及 胎 盘 体 视 学 的 影 响 段 善 娥, 张 清 德 ( 暨 南 大 学 附 属 第 一 医 院 麻 醉 科,, 彭 雪 梅, 高 雪 松 产 科, 广 东 广 州 ) [ 摘 要 ] 目 的 : 通 过 研 究 全 身 麻 醉 ( ) 与 腰 硬 联 合 麻

More information

行政院及各所屬機出國報告

行政院及各所屬機出國報告 臺 北 市 政 府 所 屬 各 機 關 因 公 出 國 人 員 出 國 報 告 書 出 國 類 別 :( 教 育 考 察 ) 臺 北 市 95 年 度 高 職 國 文 科 教 學 輔 導 團 大 陸 參 訪 考 察 報 告 書 服 務 機 關 : 臺 北 市 政 府 教 育 局 臺 北 市 政 府 研 考 會 編 號 欄 出 國 人 員 : 林 騰 蛟 副 局 長 ( 其 他 如 團 員 名 冊 )

More information

untitled

untitled 9 考 研 数 学 试 题 详 解 与 评 析 水 木 艾 迪 考 研 辅 导 班 教 务 电 话 :6755\87885 9 年 全 国 硕 士 研 究 生 入 学 统 一 考 试 数 学 二 试 题 一 选 择 题 :~8 小 题 每 小 题 分 共 分 下 列 每 小 题 给 出 的 四 个 选 项 中 只 有 一 项 符 合 题 目 要 求 把 所 选 项 前 的 字 母 填 在 题 后 的

More information

616 622 605 618. 644 757 j.. 795^ (. 806 820. 878. 880. 894 906 Li 909 958. 977. 10M ( 1036 ( )..

616 622 605 618. 644 757 j.. 795^ (. 806 820. 878. 880. 894 906 Li 909 958. 977. 10M ( 1036 ( ).. 221 P0). 208. 207 ) 196 199 208. 210 Zhi) 226 Γ 236 ( 257. 328. 535. 550. m (. 607 616 622 605 618. 644 757 j.. 795^ (. 806 820. 878. 880. 894 906 Li 909 958. 977. 10M ( 1036 ( ).. 1065 1068 1129 1131

More information

數學教育學習領域

數學教育學習領域 高 中 数 学 课 程 补 充 资 料 013/14 学 年 就 读 中 四 学 生 适 用 013 ( 空 白 页 ) 目 录 页 数 1. 概 论 1 1.1 背 景 1 1. 关 注 事 项 及 考 虑 因 素 1 1.3 短 期 方 案 摘 要 1 1.4 评 核 设 计 概 要. 修 订 后 的 高 中 数 学 课 程 学 习 内 容 3.1 修 订 后 的 必 修 部 分 学 习 内 容

More information

G(z 0 + "z) = G(z 0 ) + "z dg(z) dz z! # d" λ "G = G(z 0 ) + #cos dg(z) ( & dz ) * nv #., - d+ - - r 2 sin cosds e / r # ddr 4.r 2 #cos! "G = G(z 0 )

G(z 0 + z) = G(z 0 ) + z dg(z) dz z! # d λ G = G(z 0 ) + #cos dg(z) ( & dz ) * nv #., - d+ - - r 2 sin cosds e / r # ddr 4.r 2 #cos! G = G(z 0 ) 2005.7.21 KEK G(z 0 + "z) = G(z 0 ) + "z dg(z) dz z! # d" λ "G = G(z 0 ) + #cos dg(z) ( & dz ) * nv #., - d+ - - r 2 sin cosds e / r # ddr 4.r 2 #cos! "G = G(z 0 ) + #cos dg(z) ( & dz ) * nv 2+ + ds -

More information

$$() * * ) ) + + +, ) - ),,, ) ). /, ) ) ). /01(). /,,. / ) ), ) ), + + ) ), ) ) ) ) ), $ ( ) $ $ $ ( ) * $ $ * * (, -. -/01/. (, -. * $ ) ( + $ $ ( ) $ ** $ $ $ $ ** ** + $ ), $ $ ( )) * ( * + $ $ (

More information

山 东 大 学 信 号 与 系 统 和 数 字 信 号 处 理 (833) 考 研 内 部 精 华 资 料...27 山 东 大 学 信 号 与 系 统 和 数 字 信 号 处 理 (833)(70% 信 号 与 系 统,30% 数 字 信 号 处 理 不 含 滤 波 器 设 计 )/ 考 研 内

山 东 大 学 信 号 与 系 统 和 数 字 信 号 处 理 (833) 考 研 内 部 精 华 资 料...27 山 东 大 学 信 号 与 系 统 和 数 字 信 号 处 理 (833)(70% 信 号 与 系 统,30% 数 字 信 号 处 理 不 含 滤 波 器 设 计 )/ 考 研 内 山 东 大 学 839 生 物 化 学 ( 生 ) 全 套 考 研 资 料 ......7 山 东 大 学 生 物 化 学 (839)( 生 )/ 考 研 内 部 精 华 资 料...7 山 东 大 学 生 物 化 学 ( 生 )(839) 考 研 内 部 精 华 资 料...8 山 东 大 学 661 考 古 综 合 全 套 考 研 资 料

More information

Microsoft Word - whfq fm_new_.doc

Microsoft Word - whfq fm_new_.doc 图 灵 数 学 统 计 学 丛 书 The Calculus Lifesaver:All the Tools You Need to Excel at Calculus 普 林 斯 顿 微 积 分 读 本 [ 美 ] Adrian Banher 著 杨 爽 赵 晓 婷 高 璞 译 北 京 图 书 在 版 编 目 (CIP) 数 据 普 林 斯 顿 微 积 分 读 本 / ( 美 ) 班 纳 (Banner,

More information

untitled

untitled , 5 55-% 8-8 8-45% - 7 7 U- 4 4 4in 5 co 5 4in lim lim co 5 5 in a, d, a in lim lim 4.4..4. d www.inghaor.com 5 79 755 d d 5 7 7 7..7 4 c 7 5, d d 9.9, 4- A E B BA B E B BA B E B A E E B A E E 4 A E B

More information

R C + C = C R = 1 N (R-1) (N-1) R (N-1) N 10 + 7 4913 3 2 3 10 300 1000 2 a Λ = 3913 3agΛ 3 10 7 = 210 2 2 b Λ 7 = 49 3a2 + 3ab + b 2 2 2 Λ 559 3913Λ ( 3a + 3ab + b ) b 3 4913 = 17 223 22 71 7 3 10

More information

總目160-香港電台

總目160-香港電台 管 制 人 員 : 廣 播 處 長 會 交 代 本 總 目 下 的 開 支 二 零 一 五 至 一 六 年 度 預 算... 8.231 億 元 二 零 一 五 至 一 六 年 度 的 編 制 上 限 ( 按 薪 級 中 點 估 計 的 年 薪 值 ) 相 等 於 由 二 零 一 五 年 三 月 三 十 一 日 預 算 設 有 的 691 個 非 首 長 級 職 位, 增 至 二 零 一 六 年 三

More information

996,,,,,,, 997 7, 40 ; 998 4,,, 6, 8, 3, 5, ( ),, 3,,, ;, ;,,,,,,,,,

996,,,,,,, 997 7, 40 ; 998 4,,, 6, 8, 3, 5, ( ),, 3,,, ;, ;,,,,,,,,, ,, AB,, ( CIP) /, 000 ( /, ) ISBN 704009448 F47 CIP ( 000) 86786 55 00009 0064054588 ht tp www hep edu cn ht tp www hep com cn 006404048 787960/ 6 05 370 000 730,, 996,,,,,,, 997 7, 40 ; 998 4,,, 6, 8,

More information

土木工程施工

土木工程施工 8 9 10 11 12 V K S = V K S = V V 2 1 3 1 13 N H H H H H N i i i i i 4 ) ( 1 4 3 2 1 0 = + + + = ( ) + + + = 4 3 2 1 0 4 3 2 4 1 H H H H N H 2 3 0.002 h n = = H n H ' n ) 4 4 3 2 1 2 h h h h a V +

More information

A t 1 2 L A t 2 1.44 1 2 (2) N n + + n N q 1.2 n N - + = - 100% = 100% = 100 = PER = 1 mg = mg + mg 60 % = ( kg) 100 % = ( kg) 100 24 ( mg)

More information

dn = kn ( 5 1) dt t = 2 303 1 k tg N 0. ( 5 2) N S m + M v = V (5 3) K S dx = µ X dt 5 4 S µ = µ m K + S ( 5 5) S ds 1 dx 1 = = µ X ( 5 6) dt Y dt Y x/ s x/ s ds 1 = + + ( ) dt Y X mx 1 dp

More information

99710b44zw.PDF

99710b44zw.PDF 10 1 a 1 aa bb 4 + b ± b 4ac x a 1 1 CBED DC(BC ED) (a b) DAE CBA DAE 1 ab ABE c 1 1 (ab) c ab 3 4 5 5 1 13 7 4 5 9 40 41 11 60 61 13 84 85 m 1 m + 1 m m ( m 1 ) ( m +1 = ) () m AB (m n ) n

More information

.0 引 言

.0 引 言 第 二 章 Bayes 决 策 理 论.0 引 言. 基 于 最 小 错 误 率 的 Bayes 决 策. 基 于 最 小 风 险 的 Bayes 决 策.3 正 态 分 布 的 最 小 错 误 率 Bayes 决 策.4 说 明 .0 引 言 .0 引 言 统 计 决 策 理 论 根 据 每 一 类 总 体 的 概 率 分 布 决 定 决 策 边 界 Bayes 决 策 理 论 是 统 计 决 策

More information

koji-13.dvi

koji-13.dvi 26 13 1, 2, 3, 4, 5, 6, 7 1 18 1. xy D D = {(x, y) y 2 x 4 y 2,y } x + y2 dxdy D 2 y O 4 x 2. xyz D D = {(x, y, z) x 1, y x 2, z 1, y+ z x} D 3. [, 1] [, 1] (, ) 2 f (1)

More information

總目100-海事處

總目100-海事處 管 制 人 员 : 海 事 处 处 长 会 交 代 本 总 目 下 的 开 支 二 零 一 六 至 一 七 年 度 预 算... 1 3. 1 0 7 亿 元 二 零 一 六 至 一 七 年 度 的 编 制 上 限 ( 按 薪 级 中 点 估 计 的 年 薪 值 ) 相 等 于 由 二 零 一 六 年 三 月 三 十 一 日 预 算 设 有 的 1 3 7 6 个 非 首 长 级 职 位, 增 至

More information

. () ; () ; (3) ; (4).. () : P.4 3.4; P. A (3). () : P. A (5)(6); B. (3) : P.33 A (9),. (4) : P. B 5, 7(). (5) : P.8 3.3; P ; P.89 A 7. (6) : P.

. () ; () ; (3) ; (4).. () : P.4 3.4; P. A (3). () : P. A (5)(6); B. (3) : P.33 A (9),. (4) : P. B 5, 7(). (5) : P.8 3.3; P ; P.89 A 7. (6) : P. () * 3 6 6 3 9 4 3 5 8 6 : 3. () ; () ; (3) (); (4) ; ; (5) ; ; (6) ; (7) (); (8) (, ); (9) ; () ; * Email: huangzh@whu.edu.cn . () ; () ; (3) ; (4).. () : P.4 3.4; P. A (3). () : P. A (5)(6); B. (3) :

More information

13WuYW_4questions

13WuYW_4questions EuMath (/008 來 自 身 邊 的 四 個 小 問 題 胡 奕 偉 麗 水 學 院 數 學 系 這 是 一 組 來 自 筆 者 身 邊 的 問 題 問 題 平 凡, 問 題 簡 單, 問 題 3 略 見 抽 象, 問 題 4 則 源 遠 流 長, 被 稱 為 亞 里 斯 多 德 旋 輪 悖 論 平 凡 的 問 題 呼 喚 靈 活 的 思 維, 處 理 方 法 要 創 新 ; 貌 似 簡 單

More information

参 2

参 2 1 参 2 閲 3 4 5 { } i τ 6 7 参 参 8 9 10 11 12 13 14 15 [ ) 0, ( ], st τ { } i τ = t t ( t = 0) i i i 1 0 { } τ i { } t τ t i { } { } τ i i { } i t 0 t 1 t 2 t 3 t 4 t n τ 1 τ 2 τ3 τ 4 () λ t 16 ( ( ) ) P

More information

基础班讲义

基础班讲义 年 考 研 数 学 基 础 班 讲 义 高 等 数 学 第 一 章 函 数 极 限 连 续 一 函 数 函 数 的 概 念 : 函 数 的 性 态 : 单 调 性 奇 偶 性 周 期 性 有 界 性 有 界 性 : 定 义 : M >, I, M ; 复 合 函 数 与 反 函 数 函 数 的 复 合, 求 反 函 数 4 基 本 的 初 等 函 数 与 初 等 函 数 基 本 初 等 函 数 :

More information

① ⑰ ⒀ ⒐ ② ⑱ ⒁ ⒑ 〡 〱 ぁ ③ ⑲ ⒂ ⒒ 〢 〲 あ ④ ⑳ ⒃ ⒓ 〣 〳 ぃ ⑤ ⑴ ⒄ ⒔ 〤 〴 い ⑥ ⑵ ⒅ ⒕ 々 〥 〵 ぅ ⑦ ⑶ ⒆ ⒖ 〆 〦 う ⑧ ⑷ ⒇ ⒗ 〇 〧 ぇ ⑨ ⑸ ⒈ ⒘ 〨 ⑩ ⑹ ⒉ ⒙ 〩 ⑪ ⑺ ⒊ ⒚ ⓪ ⑯ ⑿ ⒏ え ぉ お

① ⑰ ⒀ ⒐ ② ⑱ ⒁ ⒑ 〡 〱 ぁ ③ ⑲ ⒂ ⒒ 〢 〲 あ ④ ⑳ ⒃ ⒓ 〣 〳 ぃ ⑤ ⑴ ⒄ ⒔ 〤 〴 い ⑥ ⑵ ⒅ ⒕ 々 〥 〵 ぅ ⑦ ⑶ ⒆ ⒖ 〆 〦 う ⑧ ⑷ ⒇ ⒗ 〇 〧 ぇ ⑨ ⑸ ⒈ ⒘ 〨 ⑩ ⑹ ⒉ ⒙ 〩 ⑪ ⑺ ⒊ ⒚ ⓪ ⑯ ⑿ ⒏ え ぉ お ! " @ P Q R ` a b p q r # S c s $ % & ' ( G H T U V W X d e f g h t u v w x ) I Y i y * : J Z j z + ; K [ k {, -. / < = >? L M N O \ ] ^ _ l m n o } ~ À Ð à ð ΐ Π ΰ π ϐ Ϡ ϰ Ѐ А Р а р Ѐ Ѡ Ѱ Ҁ Ґ Ҡ Ұ Ӏ Ӑ

More information

untitled

untitled 6 6 6 5 55-6% 8-8 8-5% 6-6 7 6 6 7 U- l lim co l ~ co ~ 6 6 - - ce [] 6 6 6 7 Σ z Σ ddz dzd z dd π : z Gau www.ighuauor.com 5 6796 6755 6 Σ Z z Z 6 z Σ 6dddz ΩΣΣ Ω 6V V π 6 π ddz dzd z dd Σ Σ z dz 6 678

More information

绪论

绪论 微 分 方 程 数 值 解 法 信 息 与 计 算 科 学 专 业 核 心 课 程 数 学 科 学 学 院 课 程 基 本 情 况 课 程 名 称 : 微 分 方 程 数 值 解 法 法 Nmerical Soltio of Eqatio 课 程 类 型 : 专 业 核 心 课 课 程 设 立 :995: 年 开 课 范 围 : 信 息 与 计 算 科 学 专 业 学 分 4 学 时 64 必 修

More information

在 上 述 物 理 模 型 中 ( 三 隻 猴 子 的 重 量 都 一 樣 ), 考 慮 底 下 四 個 問 題 : () 當 三 股 力 量 處 於 平 衡 狀 態, 而 且 F 點 處 於 ABC 的 內 部 時, 利 用 力 的 向 量 和 為 零 的 觀 念, 求 角 度 AFB, BFC,

在 上 述 物 理 模 型 中 ( 三 隻 猴 子 的 重 量 都 一 樣 ), 考 慮 底 下 四 個 問 題 : () 當 三 股 力 量 處 於 平 衡 狀 態, 而 且 F 點 處 於 ABC 的 內 部 時, 利 用 力 的 向 量 和 為 零 的 觀 念, 求 角 度 AFB, BFC, 許 教 授 講 故 事 許 志 農 / 國 立 台 灣 師 範 大 學 數 學 系 在 數 學 教 學 中, 有 這 樣 一 道 數 學 應 用 問 題 : 在 哪 裡 建 學 校, 可 使 附 近 的 三 個 村 子 A, 與 C 的 三 位 學 生 到 學 校 所 走 路 程 的 和 最 小? 此 問 題 實 質 為 : 給 平 面 上 A, B, C 三 點, 試 尋 求 一 點 F, 使 距

More information

lim f(x) lim g(x) 0, lim f(x) g(x),

lim f(x) lim g(x) 0, lim f(x) g(x), 2016 11 14 1 15 lim f(x) lim g(x) 0, lim f(x) g(x), 0 0. 2 15 1 f(x) g(x) (1). lim x a f(x) = lim x a g(x) = 0; (2). a g (x) f (x) (3). lim ( ). x a g (x) f(x) lim x a g(x) = lim f (x) x a g (x). 3 15

More information

虎克定律實驗 楊勝斐

虎克定律實驗           楊勝斐 虎 克 定 律 實 驗 楊 勝 斐 1. 目 地 : 測 試 彈 簧 的 虎 克 定 律, 並 從 彈 簧 作 簡 諧 運 動 的 特 性, 以 求 其 彈 性 係 數. 儀 器 : 彈 簧 一 條, 砝 碼 一 組, 虎 克 定 律 儀 一 組, 碼 錶 一 個, 米 尺 一 支 圖 1. 虎 克 定 律 儀 器 全 圖 與 零 件 圖 1 圖. 彈 簧 受 外 力 產 生 反 彈 力 與 變 形

More information

範 例 1.1 試 解 出 下 列 微 分 方 程 dx = y. 不 嚴 謹 做 法 : 把 微 分 方 程 改 寫 為 y = dx. 兩 邊 同 時 積 分 y = 之 後 可 以 推 得 : ln y = X + C, 兩 邊 同 時 取 exp 之 後 可 以 得 到 y = Ce x.

範 例 1.1 試 解 出 下 列 微 分 方 程 dx = y. 不 嚴 謹 做 法 : 把 微 分 方 程 改 寫 為 y = dx. 兩 邊 同 時 積 分 y = 之 後 可 以 推 得 : ln y = X + C, 兩 邊 同 時 取 exp 之 後 可 以 得 到 y = Ce x. 微 分 方 程 法 蘭 克 老 師 1 微 分 方 程 1.1 可 分 離 微 分 方 程 假 設 M(x), N(y) 都 是 定 義 在 某 個 區 間 上 的 連 續 函 數 我 們 希 望 解 以 下 類 型 的 常 微 分 方 程 以 不 嚴 謹 的 方 法 我 們 可 以 把 (1.1) 改 寫 成 M(x) N(y) = 0. (1.1) dx N(y) = M(x)dx. (1.2)

More information

2 a i, f, b j, Z 6. {x 1, x 2,, x n }Du = ϕx 1, x 2,, x n 7.1 F x 1, x 2,, x n ; ϕ,,,, 0 u = ϕx 1, x 2,, x n 7.1D 7.1n = 2 F x, y, z,, = 0 7.4 z = ϕx,

2 a i, f, b j, Z 6. {x 1, x 2,, x n }Du = ϕx 1, x 2,, x n 7.1 F x 1, x 2,, x n ; ϕ,,,, 0 u = ϕx 1, x 2,, x n 7.1D 7.1n = 2 F x, y, z,, = 0 7.4 z = ϕx, 1. 7.1 ux 1, x 2,, x n n 2 u, u,, u F x 1, x 2,, x n ; u, u, u,, u = 0 7.1 2. F u, u, u,, u a 0 x 1, x 2,, x n u + 3. fx 1, x 2,, x n 0 a i x 1, x 2,, x n u i = fx 1, x 2,, x n 7.2 a 0 x 1, x 2,, x n u

More information

MATHEMATICAL MODELING SPONSORED BY: SHUMO.COM COMPILED BY: Mathematical Modeling Editors Group HOMEPAGE:

MATHEMATICAL MODELING SPONSORED BY: SHUMO.COM COMPILED BY: Mathematical Modeling Editors Group  HOMEPAGE: MATHEMATICAL MODELING S H U M O SHUMO.COM 2004. 1 Vol. 1 No. 1 MATHEMATICAL MODELING SPONSORED BY: SHUMO.COM COMPILED BY: Mathematical Modeling Editors Group www.shumo.com HOMEPAGE: www.shumo.com mmjournal@yeah.net

More information

微 分 方 程 是 经 典 数 学 的 一 个 重 要 分 支, 常 用 来 描 述 随 时 间 变 化 的 动 态 系 统, 被 广 泛 应 用 于 物 理 学 工 程 数 学 和 经 济 学 等 领 域. 实 际 上, 系 统 在 随 时 间 的 变 化 过 程 中, 经 常 会 受 到 一 些

微 分 方 程 是 经 典 数 学 的 一 个 重 要 分 支, 常 用 来 描 述 随 时 间 变 化 的 动 态 系 统, 被 广 泛 应 用 于 物 理 学 工 程 数 学 和 经 济 学 等 领 域. 实 际 上, 系 统 在 随 时 间 的 变 化 过 程 中, 经 常 会 受 到 一 些 不 确 定 微 分 方 程 研 究 综 述 李 圣 国, 彭 锦 华 中 师 范 大 学 数 统 学 院, 湖 北 4379 黄 冈 师 范 学 院 不 确 定 系 统 研 究 所, 湖 北 438 pengjin1@tsinghua.org.cn 摘 要 : 不 确 定 微 分 方 程 是 关 于 不 确 定 过 程 的 一 类 微 分 方 程, 其 解 也 是 不 确 定 过 程. 本 文 主

More information

K526-ML

K526-ML 目 录 C O N T E N T S 21 世 纪 的 养 生 大 趋 势 寻 找 最 自 然 的 自 我 保 健 法 1. 新 世 纪, 新 期 望, 需 要 新 的 养 生 观 点 / 2 2. 生 物 节 律 养 生 应 运 而 生 / 3 3. 注 重 生 物 节 律 养 生, 为 健 康 加 油 / 5 第 课 探 索 充 满 节 律 的 世 界 课 前 热 身 : 判 断 2 分 钟

More information

4. 计 算 积 分 : ż ż βi fdl = f(x(t), y(t), z(t)) a x 1 (t) 2 + y 1 (t) 2 + z 1 (t) 2 dt L i α i ż ż βi 或 者 在 二 维 情 形 中 fdl = f(x(t), y(t)) a x 1 (t) 2 +

4. 计 算 积 分 : ż ż βi fdl = f(x(t), y(t), z(t)) a x 1 (t) 2 + y 1 (t) 2 + z 1 (t) 2 dt L i α i ż ż βi 或 者 在 二 维 情 形 中 fdl = f(x(t), y(t)) a x 1 (t) 2 + 微 积 分 B2 曲 面 曲 线 积 分 小 结 马 晓 光 2014 年 5 月 15 日 1 第 一 型 曲 线 曲 面 积 分 这 一 部 分 的 积 分 区 域 是 没 有 定 向 的 解 题 的 关 键 是 计 算 长 度 微 元 dl 和 面 积 微 元 ds 1.1 第 一 型 曲 线 积 分 积 分 区 域 是 一 条 曲 线 L, 可 以 在 二 维 平 面 内, 也 可 以 在

More information

号 决 议, 邀 请 收 到 大 会 长 期 邀 请 参 加 其 主 持 下 召 开 的 一 切 国 际 会 议 的 会 议 和 工 作 的 各 组 织 代 表, 以 观 察 员 身 分 参 加 ; 按 照 大 会 一 九 七 四 年 十 二 月 十 曰 第 3280(Х Х 1Х ) 号 决 议,

号 决 议, 邀 请 收 到 大 会 长 期 邀 请 参 加 其 主 持 下 召 开 的 一 切 国 际 会 议 的 会 议 和 工 作 的 各 组 织 代 表, 以 观 察 员 身 分 参 加 ; 按 照 大 会 一 九 七 四 年 十 二 月 十 曰 第 3280(Х Х 1Х ) 号 决 议, 联 合 国 海 上 货 物 运 送 会 议 最 后 文 件 1 联 合 国 大 会 在 审 议 了 联 合 Й 国 际 贸 今 法 委 (п 奂 一 九 七 六 年 # 九 М Л 代 工 作 报 告 中 载 有 海 上 货 物 运 送 公 约 草 案 的 第 四 章 之 后, 以 一 九 七 六 年 十 二 月 十 五 日 第 3 1 / 1 0 0 ^ 决 议, 决 定 于 一 九 七 八 年

More information

微积分 授课讲义

微积分 授课讲义 2018 10 aiwanjun@sjtu.edu.cn 1201 / 18:00-20:20 213 14:00-17:00 I II Taylor : , n R n : x = (x 1, x 2,..., x n ) R; x, x y ; δ( ) ; ; ; ; ; ( ) ; ( / ) ; ; Ů(P 1,δ) P 1 U(P 0,δ) P 0 Ω P 1: 1.1 ( ). Ω

More information

Microsoft Word - 实验习题N.doc

Microsoft Word - 实验习题N.doc 数 学 实 验 实 验 练 习 题 汇 总 实 验 数 学 建 模 初 步 实 验 目 的 通 过 解 决 简 化 的 实 际 问 题 学 习 初 步 的 数 学 建 模 方 法, 培 养 建 模 意 识 实 验 内 容. 怎 样 解 决 下 面 的 实 际 问 题? 包 括 需 要 哪 些 数 据 资 料, 要 做 些 什 么 观 察 试 验 以 及 建 立 什 么 样 的 数 学 模 型 等 :

More information

C35N32.dvi

C35N32.dvi 數 學 傳 播 35 卷 3 期, pp. 11-21 數 學 的 詩 篇 一 一 Fourier 分 析 林 琦 焜 深 入 研 究 大 自 然 是 所 有 數 學 發 現 最 富 饒 的 來 源, 不 僅 對 於 決 定 良 好 的 目 標 有 好 處, 也 有 助 於 排 除 含 糊 的 問 題 無 用 的 計 算 這 是 建 立 分 析 學 本 身 的 手 段, 也 協 助 我 們 發 現

More information

第6章

第6章 第 6 章 控 制 系 统 的 设 计 6. 引 言 我 们 已 经 在 前 面 几 章 学 习 了 分 析 控 制 系 统 的 方 法, 其 最 终 目 的 是 为 了 能 够 设 计 控 制 系 统 控 制 系 统 的 分 析 是 研 究 给 定 系 统 的 动 态 和 稳 态 特 性 控 制 系 统 的 设 计 则 是 一 个 相 反 的 过 程, 即 根 据 生 产 工 艺 的 要 求 设

More information

<4D6963726F736F667420576F7264202D20B9FAB7C0BFC6D1A7BCBCCAF5B4F3D1A7B9D8D3DA32303135C4EAD5D0C9FAD7A8D2B5B8FCB8C4B5C4B9ABB8E6>

<4D6963726F736F667420576F7264202D20B9FAB7C0BFC6D1A7BCBCCAF5B4F3D1A7B9D8D3DA32303135C4EAD5D0C9FAD7A8D2B5B8FCB8C4B5C4B9ABB8E6> 2015 年 硕 士 研 究 生 入 学 考 试 自 命 题 科 目 考 试 大 纲 科 目 代 码 :830 科 目 名 称 : 辐 射 防 护 与 安 全 一. 考 试 要 求 主 要 考 查 学 生 对 辐 射 防 护 与 安 全 基 本 概 念 的 理 解 与 掌 握 ; 对 辐 射 剂 量 辐 射 探 测 辐 射 防 护 辐 射 安 全 的 理 解 与 掌 握 ; 以 及 运 用 基 本

More information

<4D6963726F736F667420576F7264202D20CAA5B2C5D1A7CFB0CDF8CBAED3A120C4A3B0E52E646F63>

<4D6963726F736F667420576F7264202D20CAA5B2C5D1A7CFB0CDF8CBAED3A120C4A3B0E52E646F63> 官 方 网 站 : 圣 才 学 习 网 www.xuexi.o 免 费 咨 询 热 线 :46-3-9 物 理 化 学 主 要 公 式 及 使 用 条 件 第 一 章 气 体 的 VT 关 系 主 要 公 式 及 使 用 条 件. 理 想 气 体 状 态 方 程 式 V ( / M ) RT nrt 或 V ( V / n) RT 式 中,V,T 及 n 单 位 分 别 为 Pa, 3,K 及 ol

More information

一 耀 州 青 瓷 的 裝 飾 手 法 與 紋 飾 種 類 耀 州 窯 的 裝 飾 紋 樣, 豐 富 多 變, 而 且 題 材 內 容 廣 泛, 組 合 形 式 多 樣, 圖 案 形 象 優 美, 令 人 賞 心 悅 目, 並 且 反 映 了 當 時 社 會 的 審 美 趣 味 和 理 想 裝 飾

一 耀 州 青 瓷 的 裝 飾 手 法 與 紋 飾 種 類 耀 州 窯 的 裝 飾 紋 樣, 豐 富 多 變, 而 且 題 材 內 容 廣 泛, 組 合 形 式 多 樣, 圖 案 形 象 優 美, 令 人 賞 心 悅 目, 並 且 反 映 了 當 時 社 會 的 審 美 趣 味 和 理 想 裝 飾 宋 代 耀 州 青 瓷 的 紋 飾 風 格 與 意 義 曾 肅 良 英 國 萊 斯 特 大 學 博 物 館 學 博 士 國 立 台 灣 師 範 大 學 美 術 研 究 所 助 理 教 授 摘 要 中 國 的 飲 茶 之 風, 興 於 唐 而 盛 於 宋, 特 別 是 宋 代 宮 廷 禁 苑 和 地 方 官 吏 文 人 學 士 的 尚 茶 崇 茶, 以 品 茶 為 雅 尚 的 觀 念 與 作 法, 使

More information

Microsoft Word - 3-11.专业主干课程和主要专业课程的教学大纲.doc

Microsoft Word - 3-11.专业主干课程和主要专业课程的教学大纲.doc ... 3... 7... 10... 31... 36... 41... 49... 52... 57... 62... 65... 71... 73... 98... 105... 108... 131... 137... 148... 159... 174... 176... 192... 195... 228... 235... 242... 243 ... 247... 255... 265...

More information

Microsoft Word - 第一篇第三章_3.doc

Microsoft Word - 第一篇第三章_3.doc 第 三 章 生 产 环 境 与 健 康 生 产 劳 动 是 人 类 获 得 健 康 的 必 要 条 件 之 一, 良 好 的 生 产 环 境 和 劳 动 条 件 不 仅 对 劳 动 者 的 健 康 有 利, 也 是 促 进 生 产 顺 利 发 展 的 有 利 因 素 但 在 不 良 的 劳 动 条 件 下, 由 于 生 产 环 境 存 在 着 各 种 职 业 危 害 因 素, 则 可 使 职 业 人

More information

1 2 m v e 2 ö e m e m e m e m e e m m 1 1840 e m e m 2 v r Å Å Å 9999 10000 2 n λ = b( 2 2 ) n 2 Å 1 1854 1919 λ 1 1 1 2 2 λ = R ( Z H n ) 1 1 1 2 2 λ R H ( ) n f ni Z Z E ν = h mvr = n h 2π mvr = nh

More information

数值分析ch3.ppt

数值分析ch3.ppt 数 值 积 分 和 数 值 微 分 I = f(x)dx 由 微 积 分 学 基 本 定 理, 当 f(x) 在 [,] 上 连 续 时, 存 在 原 函 数 F(x), 由 Newto-Leiits 公 式 I= f(x)dx = F() F() 有 时, 用 上 面 的 方 法 计 算 定 积 分 有 困 难.. 不 易 求 f(x) 的 原 函 数 F(x).f(x) 的 原 函 数 表 达

More information

07HW1.mps

07HW1.mps 第 3 章 混 凝 土 结 构 的 设 计 方 法 学 习 目 标 (1) 对 结 构 的 功 能 要 求 极 限 状 态 及 结 构 可 靠 度 的 基 本 概 念 有 一 般 的 了 解 ; (2) 掌 握 结 构 设 计 中 基 本 术 语 的 定 义, 例 如 结 构 的 作 用 作 用 效 应 荷 载 的 标 准 值 和 设 计 值 混 凝 土 和 钢 筋 的 标 准 强 度 和 设 计

More information

前 言 指压按摩这种以 手 为治病工具的医疗方法 可以说是人类最早的医学形式之一 这在世界各 古老民族的医学发展史中都可以得到佐证 因为 人们以手摩擦 按揉和捶打自己或同伴的肢体 以 御寒取暖及解除疲劳 食胀与各种伤痛所造成的不适感等行为 是与人类生来就有的自卫防御本能 有关 中国的按摩疗法历史悠久 可以追溯到上古的黄帝时代 那时的按摩称案杬 明朝以前称为按 摩 明朝以后一般称为推拿 发展至今日 推拿按摩成为一门独立的学科

More information

#. #. # #. /0* # # # # # /0* /12345416 04787948: # # )*+,- *:87712 # # # # */0* # # # # # ) # * /0* # )*+,- # )*+,- * ) ) * ) )*+,- # # # /0* # # # /0

#. #. # #. /0* # # # # # /0* /12345416 04787948: # # )*+,- *:87712 # # # # */0* # # # # # ) # * /0* # )*+,- # )*+,- * ) ) * ) )*+,- # # # /0* # # # /0 - 5 # ) # *+,-./0123 ) # 011* -1 /## # 112 2 2 # # 3)4+1-*# # # /..,252)2 # -))+ - # -))2 1 - # # # /.. # # /.. # # # /.. ## # # ))* +, # # 112 2 2 -.# # # # /.. # # # # /.. # /.. # #,252)2 -))2 1 # #

More information

得 到 了 補 償. 對 於 武 姜 而 言, 莊 公 與 自 己 的 關 係 並 不 親 密, 而 共 叔 段 又 是 自 己 向 來 疼 愛 有 加 的 兒 子, 所 以, 對 莊 公 提 出 再 怎 麼 無 理 的 要 求, 武 姜 也 不 會 覺 得 有 什 麼 不 妥 之 處, 而 對 共

得 到 了 補 償. 對 於 武 姜 而 言, 莊 公 與 自 己 的 關 係 並 不 親 密, 而 共 叔 段 又 是 自 己 向 來 疼 愛 有 加 的 兒 子, 所 以, 對 莊 公 提 出 再 怎 麼 無 理 的 要 求, 武 姜 也 不 會 覺 得 有 什 麼 不 妥 之 處, 而 對 共 左 傳 - 鄭 伯 克 段 於 鄢 人 物 心 理 1021141 林 詩 倩 一. 緒 論 鄭 伯 克 段 於 鄢, 及 共 叔 段 之 亂, 是 魯 隱 公 元 年, 即 公 元 前 722 年, 春 秋 初 年 在 鄭 國 國 內 發 生 的 一 場 內 亂. 武 姜 成 為 武 公 夫 人 並 先 後 為 武 公 生 下 了 兩 個 兒 子, 長 子 莊 公 由 於 腳 先 出 來 造 成

More information

标题

标题 第八章 微积分的核心 导数与微分 内容提要:17 世纪初期,笛 卡 儿 提 出 变 量 和 函 数 的 概 念 由 此,客 观 世 界 的 运 动 变 化过程可以用数学来描述 稍后,牛顿和莱布尼兹基于直观的无穷小量,分别独立地建立 了微积分学 到了 19 世纪,柯西和维尔斯特 拉 斯 建 立 了 极 限 理 论,康 托 尔 等 建 立 了 严 格 的实数理论,使微积分学得以严密化 微积分是人类智慧的伟大结晶,极大地推动了数学

More information

杭 州 安 鸿 科 技 股 份 有 限 公 司 财 务 报 表 附 注 2013 年 1 月 1 日 至 2015 年 10 月 31 日 一 公 司 基 本 情 况 杭 州 安 鸿 科 技 股 份 有 限 公 司 ( 以 下 简 称 公 司 或 本 公 司 ) 前 身 系 杭 州 安 鸿 科 技 有 限 公 司 ( 以 下 简 称 安 鸿 有 限 公 司 ), 安 鸿 有 限 公 司 系 由 谢

More information

㲁 㱾 㲂 㱿 㲃 㲀 ኇห ᓆ㠬 敗 敜高雄市佛教會於九十八年六月六日假元亨寺舉辦 HINI新流 感消災祈安法會 由理事長淨明法師率領全體理 監事暨高雄市各寺 院法師 信眾等雲集一堂 虔誠持誦 藥師琉璃光如來本願功德經 大悲咒 一0八遍 以及佛前大供 願祈風調雨順 國泰民安 吉星 高照 遠離疫區 如意保安寧 攝影 王清江 ᐙ ᙘ 3 ኇห ᓆ㠬 妙林 第二十一卷六月號 雙(月刊 )西元二OO九年六月三十日

More information

中華民國第51屆中小學科學展覽會

中華民國第51屆中小學科學展覽會 中 華 民 國 第 53 屆 中 小 學 科 學 展 覽 會 作 品 說 明 書 國 中 組 化 學 科 030217 精 油 如 解 液, 飛 蚊 莫 摧 殘 - 真 薰 傳 學 校 名 稱 : 臺 中 市 立 居 仁 國 民 中 學 作 者 : 指 導 老 師 : 國 二 許 哲 睿 國 二 黃 冠 綸 王 翠 妃 陳 惠 如 國 二 黃 靖 惟 關 鍵 詞 : 精 油 白 線 斑 蚊 壹 摘

More information

1911 年 武 汉 起 义, 广 东 独 立 胡 汉 民 任 总 督, 陈 任 广 东 军 政 府 外 交 部 副 部 长 陈 不 愿 做 官, 几 个 月 后 即 辞 职 1915 年 与 李 煜 堂 设 立 上 海 保 险 公 司, 陈 任 主 席 1921 年 孙 中 山 就 任 非 常 大

1911 年 武 汉 起 义, 广 东 独 立 胡 汉 民 任 总 督, 陈 任 广 东 军 政 府 外 交 部 副 部 长 陈 不 愿 做 官, 几 个 月 后 即 辞 职 1915 年 与 李 煜 堂 设 立 上 海 保 险 公 司, 陈 任 主 席 1921 年 孙 中 山 就 任 非 常 大 近 代 新 会 名 人 事 迹 张 云 田 : 新 会 县 双 水 区 人, 中 国 同 盟 会 员 华 侨 镇 南 关 起 义 烈 士 张 云 田 少 年 受 其 父 教 育, 精 通 文 翰, 其 时 深 受 外 国 嘲 笑 中 华 民 族 为 东 亚 病 夫 之 辱, 因 而 弃 文 就 武, 中 武 秀 才 中 年 时 结 交 三 合 会 兄 弟, 立 志 革 清 兴 华, 参 加 孙 中

More information

BB.S92

BB.S92 目 录 怎 样 护 理 失 眠 的 患 者? (1) 对 长 期 卧 床 病 人 失 眠 者 怎 样 进 行 家 庭 护 理? (3) 安 眠 药 是 一 位 危 险 的 朋 友 吗? (5) 心 理 疗 法 为 什 么 可 以 治 疗 失 眠? (8) 失 眠 治 疗 的 误 区 有 哪 些? (9) 解 梦 为 什 么 可 以 治 疗 失 眠 症? (11) 为 什 么 暗 示 疗 法 与 催

More information

Microsoft PowerPoint - 华罗庚直接法(2).ppt

Microsoft PowerPoint - 华罗庚直接法(2).ppt 历 史 上 的 今 天 1. 1887 年 11 月 17 日, 二 战 著 名 将 领 蒙 哥 马 利 诞 辰 伯 纳 德 劳 蒙 哥 马 利 是 英 国 著 名 的 军 事 家 元 帅, 第 二 次 世 界 大 战 期 间 英 国 武 装 部 队 杰 出 的 领 导 人 之 一 1887 年 11 月 17 日, 蒙 哥 马 利 出 生 在 伦 敦 肯 宁 登 区 圣 马 克 教 区 的 一 个

More information

巫月娥网络生态视阈下感知网络商业伦理对顾客忠诚的影响 在网络商业环境中遵循的伦理准则 以往研究对 一 引言 网络商业伦理的维度没有达成一致的结论被学 者和消费者认为最重要的网络商业伦理问题是安 根据中国电子商务研究中心发布 年中 全性和隐私 此外有学者还研究网络销售中 国电子商务市场数据监测报告 截

巫月娥网络生态视阈下感知网络商业伦理对顾客忠诚的影响 在网络商业环境中遵循的伦理准则 以往研究对 一 引言 网络商业伦理的维度没有达成一致的结论被学 者和消费者认为最重要的网络商业伦理问题是安 根据中国电子商务研究中心发布 年中 全性和隐私 此外有学者还研究网络销售中 国电子商务市场数据监测报告 截 经济 管理 重庆理工大学学报 社会科学 年第 卷第 期 C U T V N j 网络生态视阈下感知网络商业伦理 对顾客忠诚的影响 巫月娥 闽南师范大学 管理学院福建 漳州 摘要 借鉴顾客忠诚的认知 情感 行为模式 针对网络生态系统中出现的伦理问题 探讨顾客 感知的网络商业伦理对顾客满意和顾客忠诚的影响 以网络商业伦理研究文献为基础 提出顾 客感知网络商业伦理包括 隐私保护 安全保障 销售诚信 交易履行和服务公正

More information

. 142 135 145 48 1 5 ξ λ ϕ1 ϕ2 ξ ε η ε ξ ε ε 1 4 499 940 43 81 385 1539 14 4279634 21007820 142 143 1979 1978

More information

<4D6963726F736F667420576F7264202D20CFEEC4BFB5B3C8BACDC5B9A4D7F7B2DFBBAED6B8B5BCCAD6B2E15F315F2E646F63>

<4D6963726F736F667420576F7264202D20CFEEC4BFB5B3C8BACDC5B9A4D7F7B2DFBBAED6B8B5BCCAD6B2E15F315F2E646F63> 目 录 党 支 部 概 述 一 党 支 部 的 性 质 01 二 党 支 部 的 作 用 01 三 党 支 部 的 基 本 任 务 01 四 党 支 部 的 组 织 形 式 02 1. 党 员 大 会 02 2. 党 支 部 委 员 会 02 3. 党 小 组 03 五 相 关 制 度 03 党 支 部 书 记 一 党 支 部 书 记 的 产 生 04 二 党 支 部 书 记 的 地 位 04 三

More information

= 6000000 6000000 6000000 x = x = 8000000 x = 8000000 0 00 x = 00 x = 00 y = 0 00 y = 0 00 6000000 9000000 6 80 60 000 8000 = 8 000 6000 = 6 000 0000 00 A. B. C. 0000 000000 A. B. C. 0 00 000 A. B. C.

More information

吴 燕 / 国 民 政 府 时 期 四 川 县 级 司 法 人 事 制 度 改 革 研 究 司 法 公 署, 或 称 法 院 司 法 处 ; 在 人 事 任 命 上, 司 法 厅 省 政 府 地 方 军 事 当 局 高 等 法 院 均 可 委 派 县 级 司 法 人 员 ; 审 判 权 的 大 小

吴 燕 / 国 民 政 府 时 期 四 川 县 级 司 法 人 事 制 度 改 革 研 究 司 法 公 署, 或 称 法 院 司 法 处 ; 在 人 事 任 命 上, 司 法 厅 省 政 府 地 方 军 事 当 局 高 等 法 院 均 可 委 派 县 级 司 法 人 员 ; 审 判 权 的 大 小 专 题 论 文 国 民 政 府 时 期 四 川 县 级 司 法 人 事 制 度 改 革 研 究 * 吴 燕 内 容 提 要 国 民 政 府 时 期, 最 高 司 法 行 政 当 局 力 图 对 县 级 司 法 人 事 制 度 进 行 改 革, 并 在 此 基 础 上 推 动 县 法 院 的 普 及 四 川 各 县 司 法 机 构 在 省 高 院 的 部 署 下 实 施 了 改 革 在 人 事 任 命

More information

Ⅰ Ⅱ Ⅲ Ⅳ

Ⅰ Ⅱ Ⅲ Ⅳ Ⅰ Ⅱ Ⅲ Ⅳ Ⅲ Ⅳ Ⅴ Ⅵ { { { { { ~ ~ α { α α α ~ ~ ~ ~ ~ ~ ~ ~ ⅠⅡⅢ Ⅰ Ⅱ Ⅲ Ⅰ Ⅱ Ⅲ { { { { { { { { { { { 珔 珔 珔 珔 ~ ~ ~ ~ ~ ~ { [ ] [ ] { 1 1 1 1 σ σ σ σ { 1 2 σ { { ⅠⅡⅢ Ⅰ Ⅱ v v v { Δ Δ Δ Δ Δ Δ { { δ δ δ {

More information

高等数学A

高等数学A 高等数学 A March 3, 2019 () 高等数学 A March 3, 2019 1 / 55 目录 1 函数 三要素 图像 2 导数 导数的定义 基本导数表 求导公式 Taylor 展开 3 积分 Newton-Leibniz 公式 () 高等数学 A March 3, 2019 2 / 55 函数 y = f(x) 函数三要素 1 定义域 2 值域 3 对应关系 () 高等数学 A March

More information

国防科学技术大学

国防科学技术大学 国 防 科 学 技 术 大 学 2016 年 硕 士 研 究 生 入 学 考 试 自 命 题 及 复 试 科 目 考 试 大 纲 国 防 科 学 技 术 大 学 研 究 生 院 2015 年 6 月 目 录 航 天 科 学 与 工 程 学 院 ( 联 系 人 : 李 人 杰,0731 84573119)... 1 科 目 代 码 :812 科 目 名 称 : 空 气 动 力 学... 1 科 目 代

More information

飞机结构设计

飞机结构设计 第 2 章 飞 机 的 外 载 荷 南 京 航 空 航 天 大 学 飞 机 设 计 技 术 研 究 所 2.1 飞 机 结 构 上 的 主 要 载 荷 2.2 不 同 飞 行 条 件 下 的 过 载 2.3 其 它 载 荷 情 况 2.4 疲 劳 载 荷 2.5 飞 机 设 计 规 范 简 介 2.1 飞 机 结 构 上 的 主 要 载 荷 飞 机 在 飞 行 起 飞 着 陆 地 面 维 护 等 使

More information

12天 本 會103年 模 範 郵 工 董 麗 珍 趙 美 珍 2人 參 加 梁 周昆法 歐陽陪興 林青豊 林秀蓮 曾文俊 甯鎮美 鄭麗娟 周肖梅 陳宏 103 年 11 月 23 日 板 橋 分 會 假 西 湖 渡 假 益 周 錦 燕 等12人 奉 准 升 遷 申 請 中 華 郵 政 村 舉 辦2

12天 本 會103年 模 範 郵 工 董 麗 珍 趙 美 珍 2人 參 加 梁 周昆法 歐陽陪興 林青豊 林秀蓮 曾文俊 甯鎮美 鄭麗娟 周肖梅 陳宏 103 年 11 月 23 日 板 橋 分 會 假 西 湖 渡 假 益 周 錦 燕 等12人 奉 准 升 遷 申 請 中 華 郵 政 村 舉 辦2 總公司分會 王瑞春 103 年 9 月 13 14 日總會假新竹市立高商 辦 理 103 年 中 華 郵 政 工 會 球類錦標賽 本會組桌 球隊 羽球隊 保齡球 隊 慢速壘球隊參加 桌球隊 羽球隊獲得全 區 甲 組 第4名 慢 速 壘 球 隊 獲 得 全 區 乙 組 第 2名 103 年 9 月 27 日 本 會 辦 理 王 功 搭 鐵 牛 車 103年11月22日婦女工作委員會暨郵工運動委員會 聯合假臺北市建

More information

WX16d22q

WX16d22q 习 近 平 对 开 展 两 学 一 做 学 习 教 育 作 出 重 要 指 示 强 调 突 出 问 题 导 向 确 保 取 得 实 际 成 效 把 全 面 从 严 治 党 落 实 到 每 一 个 支 部 两 学 一 做 学 习 教 育 工 作 座 谈 会 在 京 召 开 刘 云 山 出 席 并 讲 话 近 日, 中 共 中 央 总 书 记 国 家 主 席 中 央 军 委 主 席 习 近 平 对 在

More information

课 程 结 构 : 一 规 章 制 度 撰 写 二 劳 动 合 同 订 立 变 更 三 工 作 内 容 绩 效 管 理 四 违 纪 违 规 问 题 员 工 处 理 2

课 程 结 构 : 一 规 章 制 度 撰 写 二 劳 动 合 同 订 立 变 更 三 工 作 内 容 绩 效 管 理 四 违 纪 违 规 问 题 员 工 处 理 2 HR 如 何 精 准 预 防 用 工 风 险 和 化 解 劳 动 争 议 主 讲 : 钟 永 棣 2013 年 5 月 12 日 浙 江 杭 州 1 课 程 结 构 : 一 规 章 制 度 撰 写 二 劳 动 合 同 订 立 变 更 三 工 作 内 容 绩 效 管 理 四 违 纪 违 规 问 题 员 工 处 理 2 新 劳 动 法 下, 人 力 资 源 管 理 之 道 与 术 1 打 造 和 谐 高

More information

3.1 ( ) (Expectation) (Conditional Mean) (Median) Previous Next

3.1 ( ) (Expectation) (Conditional Mean) (Median) Previous Next 3-1: 3.1 ( )........... 2 3.1.1 (Expectation)........ 2 3.1.2............. 12 3.1.3 (Conditional Mean)..... 17 3.1.4 (Median)............ 22 Previous Next First Last Back Forward 1 1.. 2. ( ): ( ), 3.

More information

whitepaper.dvi

whitepaper.dvi π + π ϕ ϕ ϕ ϕ = ) cos( ) cos( cos cos sin sin cos 3 3 0 1 1 1 3 θ θ θ 3 3 V V q d ϕ ϕ ϕ ϕ ŵϕ ST 1+ ST n n s s + ω ςω + T m p / V K m p T V K (z - 1) ) - z(z α V 1 1 X X C LC L di d dt di q dt 1 1 R sl

More information

Microsoft Word - 吴教普〔2016〕19号.doc

Microsoft Word - 吴教普〔2016〕19号.doc 吴 教 普 2016 19 号 关 于 2016 年 初 中 毕 业 生 体 育 考 试 的 实 施 方 案 各 初 中 学 校 : 根 据 中 共 中 央 国 务 院 关 于 加 强 青 少 年 体 育 增 强 青 少 年 体 质 的 意 见 ( 中 发 2007 7 号 ) 精 神 以 及 苏 州 市 教 育 局 关 于 组 织 2016 年 苏 州 市 初 中 毕 业 生 体 育 考 试 的

More information

1965 20 4 37 19 1985 28 19 4 45 1986 38 4 61 1 1 4 114 2 1 1 4 141 4 104 2 3 171 2 3 4 3 4 3 146 3 1956 133 1974 4 1960 1974 4 1965 5 1977 1 1 1984 6 3 1965 5 1959 1961 2 1950 5 1951 1974 8 4 105 19

More information

江 苏 科 技 大 学 809 机 械 设 计 全 套 考 研 资 料 <2016 年 最 新 考 研 资 料 >......18 江 苏 科 技 大 学 810 机 械 原 理 全 套 考 研 资 料 <2016 年 最 新 考 研 资 料 >......19 江 苏 科 技 大 学 机 械 原

江 苏 科 技 大 学 809 机 械 设 计 全 套 考 研 资 料 <2016 年 最 新 考 研 资 料 >......18 江 苏 科 技 大 学 810 机 械 原 理 全 套 考 研 资 料 <2016 年 最 新 考 研 资 料 >......19 江 苏 科 技 大 学 机 械 原 江 苏 科 技 大 学 825 西 方 经 济 学 全 套 考 研 资 料 ......3 江 苏 科 技 大 学 西 方 经 济 学 (825) 全 套 考 研 资 料...3 江 苏 科 技 大 学 理 论 力 学 (801) 全 套 考 研 资 料...3 江 苏 科 技 大 学 801 理 论 力 学 全 套 考 研 资 料

More information

025-

025- < 篇 名 > 本 草 求 真 書 名 : 本 草 求 真 作 者 : 黃 宮 繡 朝 代 : 清 年 份 : 西 元 1644-1911 年 < 目 錄 > < 篇 名 > 凡 例 內 容 : 一 本 草 一 書 首 宜 分 其 形 質 氣 味 次 宜 辨 其 經 絡 臟 腑 終 宜 表 其 證 治 功 能 曆 觀 諸 書 無 不 備 載 然 理 道 不 明 意 義 不 疏 徒 將 治 效 彰 著

More information

042-

042- < 篇 名 > 雷 公 炮 製 藥 性 解 書 名 : 雷 公 炮 製 藥 性 解 作 者 : 李 士 材 朝 代 : 明 年 份 : 西 元 1588~1655 年 < 目 錄 > < 篇 名 > 序 內 容 : 余 讀 仲 景 之 敘 醫, 輒 為 之 掩 卷, 蓋 其 感 生 死 之 芒 忽, 篤 君 父 之 危 殆 賤 名 利 之 浮 榮, 冀 年 壽 以 沒 世 傷 哉 其 言, 焉 得

More information