matlab化工热力学应用,Matlab在化工热力学教学中的应用

Matlab在化工热力学教学中的应用

文 章 编 号 :1007-2853( 2012) 07-0053-04在 《 化 工 热 力 学 》 教 学 中 的 应 用Matlab邱 挺 ，黄 智 贤 ，王 红 星 ，李( 福 州 大 学 化 学 化 工 学 院 ，福 建 福 州 350108)晓摘 要 : 在 化 工 热 力 学 理 论 教 学 过 程 中 ，流 体 热 力 学 性 质 、 相 平 衡 等 计 算 较 困 难 ，经 常 需 要 迭 代 试 差 计 算 ，手 算 费 时 费 力 ． 本 文 探 讨 了 工 程 计 算 软 件 Matlab 用 于 化 工 热 力 学 的 计 算 ，结 果 表 明 : Matlab 易 学 易 用 ，可 有 效 解 决 化 工 热 力 学 计 算 问 题 ，提 高 学 生 的 工 程 计 算 能 力 ．关 键 词 : 化 工 热 力 学 ; Matlab; 状 态 方 程 ; 相 平 衡 ; 迭 代中 图 分 类 号 : G 642 文 献 标 志 码 : A点 ，频 繁 涉 及 到 非 线 性 方 程 和 方 程 组 的 求 解 ［4-5］，需 要 迭 代 试 差 计 算 ． 传 统 手 工 计 算 ，初 值 要 求 高 ， 且 计 算 量 很 大 ，用 普 通 编 程 语 言 如 C、 VB 等 编 程 ［6］，需 要 编 写 数 值 计 算 方 法 ． 而 美 国 MathW orks 公 司 开 发 的 数 学 软 件 Matlab 可 很 好 地 解 决 手 工 计 算 和 普 通 语 言 编 程 的 缺 点 ，很 方 便 地 用 于 非 线 性 方 程 ( 组 ) 计 算 、 数 据 回 归 等 ，广 泛 应 用 于 工 程 计 算 ．Matlab 包 含 了 大 量 的 计 算 算 法 ，拥 有 数 百 个 数 学 运 算 函 数 ，用 户 通 过 调 用 可 方 便 地 进 行 各 种 计 算 ［7］． 函 数 中 所 使 用 的 算 法 都 是 工 程 计 算 中 的 最 经 常 使 用 的 ，且 经 过 各 种 优 化 处 理 ，大 大 提 高 了 计 算 的 准 确 度 ． 一 般 情 况 下 ，用 Matlab 可 替 代 底 层 编 程 语 言 ，如 C 和 C + + ． 在 同 等 计 算 要 求 的 情 况 下 ，使 用 Matlab 计 算 可 极 大 减 少 工 作 量 ． 笔 者 在 数 年 化 工 热 力 学 教 学 过 程 中 ，探 索 了 Matlab 软 件 求 解 化 工 热 力 学 各 类 计 算 问 题 ，主 要 涉 及 流 体 热 力 学 及 汽 液 相 平 衡 等 方 面 的 计 算 ．化 学 工 程 与 工 艺 专 业 是 我 校 首 批 启 动 “ 卓 越工 程 师 教 育 培 养 计 划 ” 的 专 业 ． “ 卓 越 工 程 师 ” 具 备 能 融 会 贯 通 地 利 用 所 学 知 识 ，针 对 不 断 变 化 的 新 情 况 ，通 过 理 论 上 的 探 讨 并 结 合 必 要 的 科 学 验 证 ，提 出 创 新 性 解 决 工 程 问 题 的 能 力 ，同 时 要 具 有 很 强 的 自 学 能 力 ，从 而 不 断 提 高 自 身 的 工 程 素 养 ． 为 了 优 化 “ 卓 越 计 划 ” 人 才 培 养 方 案 ，精 心 设 计 课 程 体 系 ，以 回 归 工 程 应 用 为 指 导 思 想 ，将 “ 专 业 课 教 学 ” 与 “ 工 程 实 践 ” 有 机 融 合 ，本 文 以 化 工 热 力 学 课 程 教 学 为 例 ，从 培 养 学 生 工 程 计 算 能 力 出 发 ， 通 过 布 置 一 些 用 手 工 计 算 根 本 无 法 完 成 、 但 对 理 解 化 工 热 力 学 原 理 和 计 算 方 法 很 重 要 的 作 业 ，让 学 生 把 化 工 数 值 计 算 课 程 中 学 到 的 Matlab 软 件 应 用 到 工 程 实 际 中 ，从 而 使 学 生 具 备 解 决 工 程 复 杂 问 题 计 算 的 能 力 ，为 将 来 成 为 “ 卓 越 工 程 师 ” 打 下 坚 实 的 基 础 ．化 工 热 力 学 是 化 学 工 程 与 技 术 的 一 个 重 要 分 支 ，是 分 离 过 程 的 基 石 ，是 化 学 工 程 与 工 艺 专 业 的 必 修 主 干 课 程 ，是 在 具 备 必 要 的 高 等 数 学 、 物 理 化 学 和 计 算 机 技 术 等 基 础 知 识 之 后 必 修 的 专 业 技 术 基 础 课 ． 该 课 程 将 热 力 学 原 理 应 用 于 化 学 工 程 技 术 领 域 ，是 化 工 过 程 研 究 、 开 发 、 设 计 、 化 工 生 产 、 操 作 的 理 论 基 础 ，涉 及 热 力 学 物 性 参 数 和 质 量 、 能 量 传 递 的 平 衡 问 题 ，在 培 养 学 生 工 程 观 点 ，提 高 学 生 工 程 能 力 方 面 起 着 重 要 的 作 用 ［1-3］． 化 工 热 力 学 主 要 涉 及 到 流 体 热 力 学 性 质 和 流 体 相 平 衡 等 计 算 ，具 有 公 式 多 、 知 识 面 广 、 计 算 量 大 、 手 算 难 等 特Matlab 常 用 求 解 函 数1Matlab 中 单 变 量 非 线 性 方 程 求 解 使 用 fzero函 数 ，非 线 性 方 程 组 求 解 使 用 fsolve 函 数 ，求 无 约 束 多 元 函 数 最 小 值 使 用 fminsearch 函 数 ［7］．fsolve 命 令 是 求 非 线 性 方 程 组 的 解 ，也 可 用 于 求 非 线 性 方 程 的 根 ，最 简 洁 的 调 用 方 式 是 : x = fsolve( @ fun，x0 ) ． 输 入 参 数 : fun 为 用 户 定 义 的 目收 稿 日 期 :2012-05-12基 金 项 目 : 教 育 部 高 等 学 校 特 色 专 业 建 设 点 项 目 ( TS11506) ; 福 州 大 学 第 六 批 高 等 教 育 教 学 改 革 项 目作 者 简 介 :邱 挺 ( 1974-) ，男 ，福 建 福 州 人 ，福 州 大 学 教 授 ，博 士 ，主 要 从 事 传 质 与 分 离 方 面 的 研 究 ．吉 林 化 工 学 院 学 报 2012 年54标 函 数 名 ，即 非 线 性 方 程 组 f ( x) 0 中 的 f ( x)= ( 2) 编 写 求 解 文 件 ( 文 件 名 保 存 为 main． m)的 表 达 式 ，该 函 数 是 一 个 列 向 量 ，x0 是 迭 代 初 值的 向 量 ． 输 出 结 果 : x 为 待 求 解 方 程 组 的 根 向 量 ．fzero 命 令 用 于 求 单 变 量 方 程 的 根 ，所 采 用 的 算 法 主 要 是 二 分 法 、 割 线 法 和 逆 二 次 插 值 法 等 的V0 =［0． 005］; % 初 值V