• 沒有找到結果。

美国大学生数学建模 竞赛题解析与研究

N/A
N/A
Protected

Academic year: 2021

Share "美国大学生数学建模 竞赛题解析与研究"

Copied!
45
0
0

加載中.... (立即查看全文)

全文

(1)

美国大学生数学建模 竞赛题解析 与研究

美国大学生数学建模 竞赛题解析与研究

西北工业大学数学建模课题组

定价 00.00元

MCM数学建模丛书

MCM西北工业大学数学建模课题组

冰盖融化问题 数独问题 医疗保健系统评价问题

选区划分问题 飞机就座问题

肾移植问题

第 1 辑

第辑

1

学 科 类 别 : 数 学 / 建 模

http://academic.hep.com.cn

(2)

美国 MCM /ICM 竞赛指导丛书

美国大学生数学建模 竞赛题解析与研究

第 1 辑

Meiguo Dɑxuesheng Shuxue Jiɑnmo Jingsɑiti Jiexi yu Yɑnjiu Di-yi Ji

西北工业大学数学建模课题组

(3)

内容 提 要

本 系 列 丛 书 是 以 美 国 大 学 生 数 学 建 模 竞 赛 ( MCM /ICM ) 赛 题 为 主 要 研究 对 象, 结 合 竞 赛 特 等奖 的 优 秀 论 文 , 对 相 关 的 问 题 做 深 刻 细 致 的 解 析 与研 究。 本 辑 针 对 2007 年 及 2008 年 MCM /ICM 竞 赛 的 6 个 题 目 : 冰 盖 融 化问 题、 数 独 谜 题 生 成 问题 、 医 疗 保 健 系 统 评 估 问 题 、 选 区 划 分 问 题 、 飞 机就 座 问 题 以 及 肾 移 植 问题 进 行 了 解 析 与 研 究

本 书 内 容 新 颖、 实 用性 强 , 目 前 国 内 尚 无 同 类 作 品 。 本 书 可 作 为 指 导 学生 参 加 美 国 大 学 生 数 学建 模 竞 赛 的 主 讲 教 材, 也 可 作 为 本 科 生 、 研 究 生 学习 和 准 备 全 国 大 学 生、 研 究 生 数 学 建 模 竞 赛 的 参 考 书 , 同 时 可 供 研 究 相 关问 题 的 教 师 和 研 究 生 参考 使 用

  图书在版编目 (C I P) 数据

美 国大 学生数 学建 模竞赛 题解 析与研 究.第 1 辑 /西 北工业 大学 数学建 模课 题组.--北京: 高等教育出版社,

2012畅7

( 美国 MCM/ICM 竞赛指导丛书)

ISBN 978-7-04-033845-4

Ⅰ.①美… Ⅱ.①西… Ⅲ.①数学模型-竞赛题-研 究 Ⅳ.①O141畅4-44

中 国版 本图书 馆CIP 数据核字 (2012) 第 156135 号

策划编辑 刘 英    责任编辑 冯 英    封面设计 李卫青    版式设计 范晓红

插图绘制 尹 莉 责任校对 金 辉 责任印制

出版发行 高等教育出版社       咨询电话 400-810-0598 社  址 北京市西城区德外大街4 号 网  址 http://www.hep.edu.cn 邮政编码 100120      http://www.hep.com.cn 印  刷  网上订购 http://www.landraco.com 开  本 787 mm×1092 mm 1/16      http://www.landraco.com.cn

印  张 17 版  次   年 月第1 版

字  数 320 千字 印  次   年 月第 次印刷

购书热线 010-58581118 定  价 45畅00 元

本书如有缺页、倒页、脱页等质量问题,请到所购图书销售部门联系调换 版权所有 侵权必究

物 料 号 33845-00

(4)

“美国 MCM /ICM 竞赛指导丛书”

编审委员会

顾 问

Chris Arney        美国西点军校 主 编

王 杰 美 国麻省 大学 罗威尔 分校

朱 旭 西 安交通 大学

秘 书

王 嘉寅 美 国康涅 狄格 大学/西安交通大学 委 员

Jay Belanger 美 国杜鲁 门州 立大学

陈 秀珍 美 国乔治 华盛 顿大学

冯 国灿 中 山大学

龚 维博 美 国麻省 大学 阿默斯 特分 校

韩 中庚 解 放军信 息工 程大学

李 向阳 美 国伊利 诺伊 理工大 学

杨 新宇 西 安交通 大学

叶 正麟 西 北工业 大学

张 存权 美 国西弗 吉尼 亚大学

(5)

数 学建 模的训 练与 经验能 使同 学们在 解决 问题时 更有创 意, 同时也能帮助 同学们 成为 更为优 秀的 研究生。 “ 美国 MCM/ICM 竞赛指导丛书” 的出版, 将 通过数 学建 模竞赛 题目 和概念 的解 析, 帮助同学们掌握数学建模的技能, 并为 同学们 在今 后的工 作中 获得成 功打 下坚实 的基 础。

数 学建 模是一 种过 程, 也是一种理念, 或者说是一种哲学。 作为过程, 学 生在理 解及 使用建 模过 程或框 架时 需要指 导并 积累经 验。 作为经验, 学生需要 使用不 同的 数 学方 法 (离散、 连续、 线性、 非线性、 随机、 几何及分析) 构 造数学 模型, 从中体验不同的细节及复杂程度。 作为理念, 学生需要发现各种 相关的、 具有挑战性的及有趣的实际问题, 从中培养对数学建模的兴趣, 并认 识到数 学建 模在实 际生 活中的 作用。 数学建模的主要目的是指导学生用建模的 方法解 决实 际问题。 尽管在实际中, 有些问题或许可以使用已有的算法和公式 来求解, 但数学建模的方法比简单使用已有算法和公式能解决更多的问题, 特 别是解 决新 的、 没有固定答案及没有被解决过的问题。

为 了积 累经验, 同学们应尽早地接受数学建模训练, 至少应该在大学低年 级时就 开始, 这样可以在以后的课程学习中进一步强化数学建模能力。 由于数 学建模的综合与交叉 特性, 所以各个专业的学生都能够从数学建模活动中受益。

本套 丛书从将数学 模型 作为 研究 工具 的角 度出 发, 介绍模型的构造, 分析 建模过程, 这些都是帮助学生更好地掌握数学建模技能的重要因素。 数学建模是 充满挑战的高级技能, 更重要的是能够帮助学生更好地成长。 当今世界需要解决 的问题往往很复杂, 所以建立的数学模型也很复杂, 通常需要通过精细的计算和 模拟才能获得解答或 对模型结果的分析与 检验。 由于数据可视化技术的普及, 解 题方法的增加, 所以现在的确是培养更多数学建模高手的最佳时期。

我 希望 同学们 在数 学建模 探索 中取得 进步, 也希望指导教师在使用这套丛 书提供 的例 子及方 法指 导学生 时取 得 很好 的 效 果。 尽管学生的层次可能不同,

但我对 你们 的忠告 是同 样的: 树立你的信心、 发展你的技能, 用你的才能解决 社会中 最具 挑战及 最重 要的问 题。 祝各位建模好运!

Chris Arney, 博士 美国西 点军 校数学 系教 授 ICM 竞赛主席 2011 年 10 月

(6)

Forward

Undergraduate students who receive instruction and experiences in mathematical modeling become better and more creative problem solvers and graduate students.This book series is being published to prepare and educate students on the topics and concepts of mathematical model- ing to help them establish a problem solving foundation for a successful career.

Mathematical modeling is both a process and a mindset or philoso- phy.As a process, students need instruction and experience in under- standing and using the modeling process or framework.As part of their experience, they need to see various levels of sophistication and com- plexity, along with various types of mathematical structures ( discrete,

continuous, liner, nonlinear, deterministic, stochastic, geometric,

and analytic) .As a mindset, students need to see problems that are relevant, challenging, and interesting so they build a passion for the process and its utility in their lives.A major goal in modeling is for students to want to model problems and find their solutions.Recipes for structured or prescribed problem solving ( canned algorithms and formu- las) do exist in the real world, but mathematical modelers can do much more than execute recipes or formulas.Modelers are empowered to solve new, open, unsolved problems.

In order to build sufficient experience in modeling, student expo-

sure must begin as early as possible-definitely by the early undergradu-

ate years. Then the modeling process can be reinforced and used

throughout their undergraduate program.Since modeling is interdisci-

(7)

plinary, students from all areas of undergraduate study benefit from this experience.

The articles and chapters in this series expose the readers to model construction, model analysis, and modeling as a research tool.All these areas are important and build the students摧modeling skills.Mod- eling is a challenging and advanced skill, but one that is empowering and important in student development.In today摧 s world, models are often complex and require sophisticated computation or simulation to provide solutions or insights into model behavior.Now is an exciting time to be a skilled modeler since methodology to provide visualization and find solutions are more prevalent and more powerful than ever before.

I wish the students well in their adventure into modeling and I like- wise wish faculty well as they use the examples and techniques in this book series to teach the modeling process to their students.My advice to all levels of modelers is to build your confidence and skills and use your talents to solve society摧 s most challenging and important problems.

Good luck in modeling!

Chris Arney, PhD

United States Military Academy at West Point Professor of Mathematics

Director of the Interdisciplinary Contest in Modeling October, 2011

Ⅱ Forward

(8)

丛书简介

美 国大 学生 数 学 建 模 竞 赛 ( The Mathematical Contest in Modeling, MCM/

The Interdisciplinary Contest in Modeling, ICM) , 即 “数学建模竞赛” 和 “交叉 学科建 模竞 赛” , 是一 项国际 级的 竞赛 活动, 为现 今各类 数学 建模 竞赛 之 鼻祖。  

1985 年, 在美国教育部的资助下, 在美国针对在校大学生创办了一个名 为 “数学建模竞赛” 的竞赛, 其宗旨是鼓励大学师生对不同领域的各种实际 问题进 行阐 明、 分析并提出解决方案。 它是一种完全公开的竞赛, 参赛形式为 学生三 人组 成一队, 在三天 (72 小时) ( 近年改为四天, 即 96 小时) 内任选 一题, 完成数学建模的全过程, 并就问题的重述、 简化和假设及其合理性的论 述、 数学模型的建立和求解 (及软件)、 检验和改进、 模型的优缺点及其可能 的应用 范围 的自我 评述 等内容 写出 论 文。 MCM/ICM 非常重视解决方案的原创 性、 团队合作与交流以及结果的合理性。 由专家组成的评阅组进行评阅, 评出 优秀论 文。 除了不允许在竞赛期间与团队以外的任何人 ( 包括指导教师) 讨 论赛题 之外, 允许使用图书资料、 互联网上的资料、 任何类型的软件等各种资 料和途 径, 从而为参赛学生提供了广阔的创作空间。 第一届竞赛时, 只有美国 的158 个队参加, 其中只有 90 个队提交了解答论文。 2012 年 MCM/ICM 共有 5 026 个队参加, 其中 MCM 有 3 697 个队, ICM 有 1 329 个队, 遍及五大洲。

MCM/ICM 已经成为最著名的国际大学生竞赛之一, 影响极其广泛。

近 年来, 已有越来越多的中国学生组队参加美国大学生数学建模竞赛, 其 中不乏 被评 为最优 论 文 (Outstanding Winners) 的佼佼者, 这充分显示了中国 大学生 参加MCM /ICM 竞赛的积极性与实力。 学生在准备竞赛的时候, 除了在 指导教 师的 帮助下 阅读 和研究 以往 竞赛的 优胜 论文以 外, 普遍希望能有一些专 门针对 美国 大学生 数学 建模竞 赛的 书籍, 指导和帮助备赛。

“ 美国大学生数学建模竞赛指导丛书” 就是为了满足读者的这一需求而出 版的, 目的是帮助学生学习从全局出发, 不受固定模式的限制, 用建模的手段 解决开 放性 问题的 研究 方法, 并提高写作能力。 丛书的读者对象包括参赛学生 以及对 数学 建模与 算法 感兴趣 的研 究生、 专业人员和业余爱好者。

(9)

我 们邀 请 到 美 国 麻 省 大 学 罗 威 尔 分 校 王 杰 教 授 担 任 丛 书 主 编, 他曾 为 MCM/ICM 命题, 对竞赛具有很多独到的认识。 丛书作者来自各高校, 他们都 是有经 验的 指导 教师 或 参 加 过竞 赛 的 优秀 成 员。 丛书包括一本 枟正确写作美 国大学 生数 学建模 竞赛 论 文枠 和若干辑 枟 美国大学生数学建模竞赛题解析与 研究枠 , 前者为一本指导学生如何正确写作 MCM /ICM 论文的工具书, 后者中 和每一 辑将 讨论若 干赛 题, 包括问题的背景、 分析技巧、 建模与测试方法及算 法设计, 并引导读者列出进一步研究的课题。 目标是培养学生多方面的能力,

如数学、 编程、 写作及课题研究等的训练, 提高学生分析问题、 解决问题的水 平。

丛 书的 出版计 划得 到了 美 国 建 模 专 家 的 广 泛 支 持, ICM 主席 ( 也是 ICM 的发起 人)、 美国西点军校数学系教授 Chris Arney 博士受邀担任丛书顾问并为 丛书作 序。

我 们热 切希望 通 过 本 套 丛 书的 出 版, 进一步活跃我国大学生参加 MCM/

ICM 的积极性, 提高他们的自信心, 并最终取得满意的成绩。 更为重要的, 提 高学生 的研 究解决 实际 问题的 能力。

丛 书 简 介

(10)

前言

美 国大学 生数 学建模 竞赛 (MCM /ICM) 作为一项国际性数学赛事, 尤其 受到中 国大 学生们 的青 睐。 历年赛题的题材广泛, 其应用性、 探究性、 开放性 和挑战 性或 许是吸 引学 生眼球 并积 极参赛 的重 要原因, 同时也引发了国内不少 大学数 学教 师的兴 趣。

2007 年和 2008 年的赛题就有这些特点。 2007 年 MCM 的 A 题 “ Gerryman- dering” 是因探讨美国州众议员选举的公平性而引发的选区划分的社会学研究 问题, 该年 MCM 的 B 题 “ The Airplane Seating Problem”, 则为研究飞机登机 时 间 的 最 优 化 问 题。 2007 年 ICM 的 C 题 “ Organ Transplant: The Kidney Exchange Problem”, 要研究如何建立一个有效机制, 以保证有效利用与开发肾 源。 2008 年 MCM 的 A 题 “ Take A Bath” 要求研究因北极冰雪消融所导致的 一系 列环境、 生态等方面的问题, 该年 MCM 的 B 题 “Creating Sudoku Puzzles”

研究的 是数 独游戏 的反 问题: 数独谜题难度级别的划分与生成算法及其复杂性 分析。 2008 年 ICM 的 C 题 “Finding the Good in Health Care Systems”, 则是研 究在医 疗保 健领域 中, 如何科学合理地评价服务系统方面的问题。

这 些赛 题的前 沿性 和现实 性乃 至趣味 性十 分明显。 对于问题背景的归纳表 述和所 提出 的希望 解决 的问题, 都相当贴近原始问题, 在专业上的简化和在数 学上的 提炼 都比较 少, 问题实际、 明了。 有些问题本身就是一个研究课题, 还 没有成 熟的 研究成 果, 甚至即使有研究结论, 也尚存争议, 未有定论, 有不少 值得进 行探 索和研 究的 方向和 方面, 极具挑战性。 这些赛题往往要求参赛学生 自己查 询、 获取充分的相关信息和数据。 针对赛题的数学建模与算法多种多 样, 开放性强, 参赛学生能够发挥的空间大。

分 析这两 年MCM/ICM 获奖的优秀论文, 它们具有下面的一些特点: 或侧 重于模 型和 方案的 新颖 性, 自成一体; 或建模分析中对于问题的理解和解决有 独特的 视野 和观点; 或表现出有价值的外延与扩展。 论文在建立数学模型解决 问题过 程中, 或者对于问题的全局及各因素的复杂关系有较好的分析与把握;

或者未 必做 得全面, 却往往对其中一两个问题颇有特色的有层次的深入的研 究。 对计算结果似乎比较淡化, 然而强调有说服力的量化结论。 尤其注重新见

(11)

解、 新观点、 新意或创新性。

MCM/ICM 的这些特点, 无疑对于培养大学生优良的科学素质、 独立的研 究能力、 创造性能力都是有益的。 就大学数学教学而言, 或多或少推动着传统 数学教 学模 式的变 化, 使得大学校园里的数学课程的教与学, 变得生动有趣、

富有活 力。 师生走出课堂, 走出书本, 走向实际问题, 走向应用探索, 致力于 开放教 学, 致力于自主学习。

本 辑就2007、 2008 年 MCM /ICM 的 A、 B 与 C 题共 6 道赛题的若干数学建 模方法, 结合当年获奖的优秀论文进行介绍和评析, 希望有益于大学的数学建 模教学 和竞 赛活动。 书中不免有诸多不妥之处, 敬请读者批评指正。 本辑抛砖 引玉, 希冀有更多、 更好的相关书籍问世。

本 辑的第1 章由孙中奎撰写, 肖华勇撰写第 2 章, 徐根玖撰写第 3 章, 雷 佑铭撰 写第4 章, 郭千桥撰写第 5 章, 吕全义撰写第 6 章。 叶正麟负责组稿和 统稿。

这 套丛 书是由 美国 麻省大 学罗 威尔分 校计 算机科 学系王 杰教 授主编, 王杰 教授对 于数 学模型 与算 法有很 深刻 的研究, 对美国大学生数学建模竞赛和学生 的科学 素质 培养饶 有兴 趣。 在本辑的撰写过程中, 王教授给予了富有启示的建 议, 西安交通大学领导提供了诸多便利, 高等教育出版社刘英副编审做了鼎力 策划, 冯英女士做了认真辛勤的编辑工作, 得到了西北工业大学教务处领导及 李辉副 处长 和应用 数学 系孙浩 教授 的有力 支持, 在此一并表示诚挚的感谢。 很 欣赏ICM 主席、 美国西点军校数学系教授 Chris Arney 博士所作的序, 并深表 感谢。

叶正 麟 2012 年 6 月于西安

  言

(12)

目录

第1 章 冰盖融化问题 …… …… ……… …… ……… …… ……… …… ……… 1 1畅1 问题的综述… ……… …… ……… …… ……… …… ……… …… ……… 1 1畅2 问题的背景资料…… …… ……… …… ……… …… ……… …… ……… 2 1畅2畅1 全球气候正在变暖… ……… …… ……… …… ……… …… ……… 2 1畅2畅2 气候变暖的影响 … ……… …… ……… …… ……… …… ……… 4 1畅2畅3 为什么是北极 …… ……… …… ……… …… ……… …… ……… 6 1畅3 问题的数学模型与结果分析…… …… ……… …… ……… …… ……… 7 1畅3畅1 模型一: 多元线性回归模型 …… ……… …… ……… …… ……… 7 1畅3畅2 模型二: 初等数学模型 … …… ……… ……… …… ……… …… 21 1畅4 问题的综合分析与进一步研究的问题 …… ……… …… ……… …… 25 1畅4畅1 问题的综合分析 … ……… …… ……… ……… …… ……… …… 25 1畅4畅2 进一步研究的问题 ……… …… ……… ……… …… ……… …… 28 参 考文 献…… …… ……… …… ……… …… ……… ……… …… ……… …… 29 第2 章 数独谜题生成问题 …… ……… …… ……… ……… …… ……… …… 30 2畅1 数独的概念 ……… …… ……… …… ……… ……… …… ……… …… 30 2畅2 问题的提出 ……… …… ……… …… ……… ……… …… ……… …… 31 2畅3 问题的背景与现状 …… ……… …… ……… ……… …… ……… …… 32 2畅4 人工求解数独的规则 … ……… …… ……… ……… …… ……… …… 34 2畅5 数独难度的划分方法 … ……… …… ……… ……… …… ……… …… 52 2畅5畅1 方法一: 步数法 … ……… …… ……… ……… …… ……… …… 52 2畅5畅2 方法二: 评分法 … ……… …… ……… ……… …… ……… …… 57 2畅5畅3 方法三: 次数加权法 …… …… ……… ……… …… ……… …… 59 2畅5畅4 方法四: 机器求解客观法… …… ……… ……… …… ……… …… 64 2畅6 数独的生成问题 … …… ……… …… ……… ……… …… ……… …… 77 2畅6畅1 终盘生成问题 …… ……… …… ……… ……… …… ……… …… 77

(13)

2畅6畅2 不同级别数独问题的生成… …… ……… ……… …… ……… …… 77 2畅7 问题的综合分析与进一步研究的问题 …… ……… …… ……… …… 78 2畅7畅1 问题的综合分析 … ……… …… ……… ……… …… ……… …… 78 2畅7畅2 进一步研究的问题 ……… …… ……… ……… …… ……… …… 79 2畅7畅3 其他数独…… …… ……… …… ……… ……… …… ……… …… 79 参 考文 献…… …… ……… …… ……… …… ……… ……… …… ……… …… 83 第3 章 医疗保健系统评估问题 ……… …… ……… ……… …… ……… …… 85 3畅1 问题的综述 ……… …… ……… …… ……… ……… …… ……… …… 85 3畅2 问题的背景资料 … …… ……… …… ……… ……… …… ……… …… 88 3畅3 医疗保健系统评价模型及求解 …… ……… ……… …… ……… …… 90 3畅3畅1 模型一: 基于层次分析法的灰色关联评估模型… …… ……… …… 90 3畅3畅2 模型二: 输入输出系统评估模型 …… ……… …… ……… …… 103 3畅3畅3 模型三: 医疗保健系统概率模型 …… ……… …… ……… …… 114 3畅4 问题的综合分析与进一步研究的问题…… ……… …… ……… …… 120 3畅4畅1 三个评估模型的对比…… ……… …… ……… …… ……… …… 120 3畅4畅2 其他评估方法分析… …… ……… …… ……… …… ……… …… 121 3畅4畅3 进一步研究的问题… …… ……… …… ……… …… ……… …… 126 参 考文 献 … ……… …… ……… …… ……… …… ……… …… ……… …… 127 第4 章 选区划分问题…… ……… …… ……… …… ……… …… ……… …… 129 4畅1 问题的综述… …… ……… …… ……… …… ……… …… ……… …… 129 4畅1畅1 问题的提出 ……… …… ……… …… ……… …… ……… …… 129 4畅1畅2 问题的背景资料 … …… ……… …… ……… …… ……… …… 130 4畅2 问题的数学模型与结果分析… ……… …… ……… …… ……… …… 133 4畅2畅1 模型一: 拓扑划分模型 … ……… …… ……… …… ……… …… 133 4畅2畅2 模型二: Voronoi 图划分模型 …… …… ……… …… ……… …… 139 4畅2畅3 模型三: 增长分割模型 … ……… …… ……… …… ……… …… 149 4畅3 问题的综合分析与进一步研究的问题…… ……… …… ……… …… 158 4畅3畅1 问题的综合分析 … …… ……… …… ……… …… ……… …… 158 4畅3畅2 进一步研究的问题… …… ……… …… ……… …… ……… …… 159 参 考文 献 … ……… …… ……… …… ……… …… ……… …… ……… …… 159 第5 章 飞机就座问题…… ……… …… ……… …… ……… …… ……… …… 161 5畅1 问题的综述… …… ……… …… ……… …… ……… …… ……… …… 161 5畅1畅1 问题的提出 ……… …… ……… …… ……… …… ……… …… 161

  录

(14)

5畅1畅2 问题的背景资料 … …… ……… …… ……… …… ……… …… 162 5畅2 问题的数学模型与结果分析… ……… …… ……… …… ……… …… 164 5畅2畅1 模型一: 随机模型… …… ……… …… ……… …… ……… …… 164 5畅2畅2 模型二和模型三: 数组模型和图论模型 ……… …… ……… …… 177 5畅2畅3 模型四: 元胞自动机模型 ……… …… ……… …… ……… …… 184 5畅2畅4 模型五: 模拟 (仿真) 模型 …… …… ……… …… ……… …… 193 5畅2畅5 下机模型 … ……… …… ……… …… ……… …… ……… …… 204 5畅3 问题的综合分析与进一步研究的问题…… ……… …… ……… …… 206 5畅3畅1 问题的综合分析 … …… ……… …… ……… …… ……… …… 206 5畅3畅2 进一步研究的问题… …… ……… …… ……… …… ……… …… 209 参 考文 献 … ……… …… ……… …… ……… …… ……… …… ……… …… 209 第6 章 肾移植问题… …… ……… …… ……… …… ……… …… ……… …… 212 6畅1 问题的综述… …… ……… …… ……… …… ……… …… ……… …… 212 6畅2 问题的背景资料… ……… …… ……… …… ……… …… ……… …… 217 6畅2畅1 肾移植概况 ……… …… ……… …… ……… …… ……… …… 217 6畅2畅2 影响移植的因素 … …… ……… …… ……… …… ……… …… 218 6畅2畅3 扩大肾源的方法 … …… ……… …… ……… …… ……… …… 219 6畅2畅4 肾移植的伦理与政策问题 ……… …… ……… …… ……… …… 220 6畅3 问题一: 美国移植网络的模拟……… …… ……… …… ……… …… 221 6畅3畅1 方法一: 随机排队服务模型 …… …… ……… …… ……… …… 222 6畅3畅2 方法二: 多排队多服务模型 …… …… ……… …… ……… …… 225 6畅3畅3 方法三: 带约束微分方程模型 … …… ……… …… ……… …… 229 6畅3畅4 计算结果与分析 … …… ……… …… ……… …… ……… …… 230 6畅3畅5 政策的改进 ……… …… ……… …… ……… …… ……… …… 230 6畅3畅6 几种方法的评价 … …… ……… …… ……… …… ……… …… 232 6畅4 问题二: 其他国家的移植政策及影响…… ……… …… ……… …… 232 6畅4畅1 欧洲各国的移植政策…… ……… …… ……… …… ……… …… 233 6畅4畅2 关于当前肾移植政策的建议 …… …… ……… …… ……… …… 233 6畅5 问题三: 肾循环交换…… …… ……… …… ……… …… ……… …… 235 6畅5畅1 环与链 …… ……… …… ……… …… ……… …… ……… …… 235 6畅5畅2 优化模型 … ……… …… ……… …… ……… …… ……… …… 236 6畅5畅3 算法… …… ……… …… ……… …… ……… …… ……… …… 237 6畅5畅4 计算结果与分析 … …… ……… …… ……… …… ……… …… 238 6畅5畅5 肾源的进一步利用… …… ……… …… ……… …… ……… …… 238 6畅5畅6 优化方法的进一步讨论 … ……… …… ……… …… ……… …… 239

Ⅲ 目  录

(15)

6畅6 问题四: 患者的选择…… …… ……… …… ……… …… ……… …… 240 6畅6畅1 方法一: 期望存活时间目标 …… …… ……… …… ……… …… 240 6畅6畅2 方法二: 期望存活时间动态处理 …… ……… …… ……… …… 242 6畅6畅3 方法三: 考虑存活 5 年的概率 … …… ……… …… ……… …… 243 6畅6畅4 三种方法的评价 … …… ……… …… ……… …… ……… …… 244 6畅7 问题五: 伦理与政策问题…… ……… …… ……… …… ……… …… 244 6畅7畅1 伦理问题 … ……… …… ……… …… ……… …… ……… …… 244 6畅7畅2 政策问题 … ……… …… ……… …… ……… …… ……… …… 245 6畅7畅3 结论… …… ……… …… ……… …… ……… …… ……… …… 247 6畅8 问题六: 捐赠者的决定… …… ……… …… ……… …… ……… …… 248 6畅8畅1 方法一: 消费决定模型 … ……… …… ……… …… ……… …… 248 6畅8畅2 方法二: 考虑风险模型 … ……… …… ……… …… ……… …… 251 6畅8畅3 两种方法的评价 … …… ……… …… ……… …… ……… …… 253 6畅9 问题的综合分析与进一步研究的问题…… ……… …… ……… …… 253 6畅9畅1 问题的综合分析 … …… ……… …… ……… …… ……… …… 253 6畅9畅2 进一步研究的问题… …… ……… …… ……… …… ……… …… 253 参 考文 献 … ……… …… ……… …… ……… …… ……… …… ……… …… 254

  录

(16)

第 1 章 冰盖融化问题

1畅1 问题的综述

“ Take a Bath” 是 2008 年美 国大学 生数学 建模 竞赛 (The Mathematical Contest in Modeling, MCM2008) 的 A 题, 研究的是北极冰雪消融所导致的一 系列环 境、 生态等方面的问题。 该题目如下。

大 陆遭受 侵蚀

研 究因 全球气 温升 高所造 成的 北极冰 盖融 化对大 陆的影 响。 具体而言, 为 未来50 年中每 10 年因北极冰盖融化对佛罗里达州海岸所造成的影响建立模 型, 尤其要关注对大型都市所造成的影响。 就此问题给出适宜的回答, 对所用 数据作 详细 讨论应 是重 要的组 成部 分。

问 题的 原文如 下:

Take a Bath

Consider the effects on land from the melting of the north polar ice cap due to the predicted increase in global temperatures.Specifically, model the effects on the coast of Florida every ten years for the next 50 years due to the melting, with particular attention given to large metropolitan areas.Propose appropriate responses to deal with this.A careful discussion of the data used is an important part of the an- swer.

( 题 目 来 源 : ht tp: //www.comap.com /undergraduat e/contests /mcm /contests /2008 /problems /)

显 然, 出题人希望参赛者能够收集与北极冰雪变化相关的数据, 讨论北极 冰雪融 化对 佛罗里 达海 岸今后50 年的影响。 看到题目, 可能每个人都会产生 这样的 两个 疑问: ①为什么研究的是北极冰雪融化的影响? 毕竟南极冰盖无论 从规模 上还 是面积 上都 远远超 过北 极。 ②为什么是对美国佛罗里达州海岸的影 响? 众所周知, 受到海平面上升影响最大的是荷兰以及太平洋上的一些岛国。

第二个 问题 很容易 回答, 因为出题人是美国人, 而且佛罗里达州海岸是美国人

(17)

享受海 滩度 假的最 大的 和最佳 的地 点之一, 同时也是美国发生飓风最频繁的地 区之一。 第一个问题的答案则隐藏在下面的数据和资料中。

1畅2 问题的背景资料

越 来越 多的证 据表 明, 全球气候正在变暖, “我们正居住在一个急速变暖 的星球 上!” 实际上, 早在 1992 年, 科学家关怀社会联盟的 1 700 位世界顶尖 科学家 包括 大多数 诺贝 尔科学 奖得 主, 就已经在极力呼吁世界各国联手遏制全 球暖化, “地球上的二氧化碳浓度已逼近警戒线!” “ 我们可以避免灾难发生的 时间已 经剩 下不到 十年 抑或 几十 年” “ 我们所面临的威胁一旦无法挽救, 人类 的希望 将非 常渺茫。”[ 1]

全 球升 温会带 来什 么影响? 英国作家林纳斯所著的 枟6 度变化, 一个愈来 愈热星 球变 化的未 来枠 一书中呈现了该景象: 如果全球

升 温1 ℃, 无冰海域吸收更多的热气, 加速全球暖化效应; 地球表面 1 /3 的水资 源流 失; 低海岸地区遭海水淹没。

升 温2 ℃, 欧洲居民中暑而亡; 森林被大火吞噬; 处于逆境的植物开始释 放碳, 不再具备吸碳功能; 有 1 /3的物种濒临灭绝。

升 温3 ℃, 从植物和土壤中排出的碳物质, 加速全球暖化效应; 亚马逊热 带雨林 荡然 无存; 超级飓风袭击沿海城市; 非洲闹饥荒。

升 温4℃, 永冻土无止境地溶解, 造成全球暖化效应一发不可收拾; 英国 大部分 地方 也因严 重的 水患而 不适 合居住; 地中海区域成为废墟。

升 温5 ℃, 甲烷从海床窜出, 加速全球暖化效应; 两极冰层融化; 人类逐 食物而 居, 但徒劳无功, 形同野生动物在这片土地上苟延残喘。

升 温6 ℃, 地球上的生物会在狂风、 山洪暴发、 硫化氢毒气及带着原子弹 般威力 的甲 烷火球 流窜 地表的 灾难 中完全 灭迹; 唯一存活下来的只有霉菌。

这 是危 言耸听 吗?

现 在有 三个问 题: 全球气候是否在变暖? 如果全球变暖, 会有什么影响?

因变暖 所造 成的冰 雪融 化最严 重的 是哪些 地方?

1畅2畅1 全球气候正在变暖

时 至今 日, 仍然有很多人在质疑一个问题: 全球气候真的在变暖吗? 联合 国政 府 间 气 候 变 化 专 门 委 员 会 ( Intergovemmental Panel on Climate Change,

IPCC) 的回答是: “气候系统变暖是毋庸置疑的, 目前从全球平均气温和海温 升高, 大范围积雪和冰融化, 全球平均海平面上升的观测中可以看出气候系统 变暖是 明显 的” (IPCC 2007 年评估报告[2])。 根据 IPCC 的结论和观测数据,

自1850 年以来, 在 1995—2006 年的 12 年中, 有 11 个年份位列最暖的 12 个年份

2 第1 章 冰盖融化问题

(18)

之中。 而在 1906—2005 年的 100 年中, 温度以 0畅74 ℃的趋势线性增长。

在 图1-1 中, 已观测到全球平均地表温度、 分别来自验潮仪 ( 深色) 和 卫星 (浅色) 的全球平均海平面以及 3—4 月北半球积雪的变化。 图中所有变 化差异 均相 对于1961—1990 年的相应平均值。 各平滑曲线表示 10 年平均值,

各圆点 表示 年平均 值。 阴影区为不确定性区间, 根据已知的不确定性和时间序 列综合 分析 估算得 出 ( 图和数据均来源于文献 [2])。

图1-1 温度、 海平面和北半球积雪变化

( 图 片来 源 : IPCC 2007 年 评 估 报 告 )

1畅2 问题的背景资料 3

(19)

图1-2 是观测到的全球海洋和大陆温度变化与使用自然强迫和人为强迫 的气候 模式 模拟结 果的 比较。 相对于 1901—1950 年相应的平均值, 图 1 -2 给 出了1906—2005 年观测到的 10 年平均值 ( 黑线), 对应于该 10 年的中心绘 制。 虚线部分表示空间覆盖率低于 50%。 深色阴影表示仅使用太阳活动和火 山自然 强 迫的5 个气候模式 19 个模拟试验结果的 5% ~95%可信度范围。 浅 色阴影 表示 同时使 用自 然强迫 和人 为强迫 的14 个气候模式 58 个模拟试验结果 的5% ~95%可信度范围。

图1-2 全球海洋和大陆的温度变化

( 图 片来 源 : IPCC 2007 年 评 估 报 告 )

显 然, 从图 1-2 中可以看到, 不管是全球平均温度, 还是全球陆地、 全 球海洋 的平 均温度 都呈 现波动 上升 趋势, 尤其是最近几十年, 这种上升趋势更 加明显, 平均上升幅度更大。 而且, 北极的冰雪覆盖面积和范围不断缩小, 全 球平均 海平 面高度 不断 攀升。

1畅2畅2 气候变暖的影响

全 球气 候变暖 会 带 来什 么 影 响? 温度的升高带来的影响是全面的、 深入 的, 将涉及自然、 生态、 环境、 社会、 经济等各领域。 下面的具体影响摘自

4 第1 章 冰盖融化问题

(20)

IPCC 在 2007 年的评估报告。

1畅区域尺度的变化

瞯变暖区域扩大。 陆地上和北半球大部分高纬度地区变暖幅度最大, 南半 球海洋 地区 和北大 西 洋 部分 地 区 变暖 幅 度 最 小。 近期观测到的各种趋势仍在 持续。

瞯冰雪面积缩小。 大部分多年冻土区域的融化深度增加, 海冰面积退缩。

在利用SRES 情景所作的某些预估中, 到 21 世纪后半叶, 北冰洋夏季后期的 冰几乎 全部 消失。

瞯热极端事件、 热浪以及强降水的频率很可能增加。

瞯热带气旋强度可能增加, 全球热带气旋数量减少具有较小的可信度。

瞯温带风暴路径向极地推移, 造成风、 降水和温度形态的变化。

瞯高纬度地区降水很可能增加, 大部分亚热带陆地区域降水可能减少, 已 观测到 的趋 势仍在 持续。

2畅受影响的系统和行业 瞯特定的生态系统

① 陆地: 对变暖敏感的苔原、 北方森林和山区, 受降雨减少影响的地中 海类型 的生 态系统 和受 降雨减 少影 响的热 带雨 林。

② 海岸带: 受到多重压力的红树林和盐沼。

③ 海洋: 受到多重压力的珊瑚礁, 对变暖敏感的海冰生物群落。

瞯由于降雨和蒸散的变化, 中纬度一些干旱地区和干旱的热带地区以及那 些依靠 冰雪 融化地 区的 水资源。

瞯面临可用水减少的低纬度地区的农业。

瞯地势低洼的沿海系统, 受到海平面上升的威胁和极端天气事件加大的 风险。

瞯适应能力低的人群的身体健康。

IPCC 的表述可能比较专业, 或许下面的一组数据更能触动人类 “坚强”

的神经。

瞯到 2050 年, 世界上将有 2 /3 的北极熊死亡, 因为它们赖以生存的北极 海冰受 到全 球变暖 的影 响而逐 渐消 融。 此外, 格陵兰冰盖消融加快还将在当地 引发地 震。 (来源于 http: //www.chinadaily.com.cn /hqbl/2007-09/10 /content_

6093810.htm? anchor =1)

瞯几亿非洲居民和几千万拉美居民在不到 20 年内将出现饮用水短缺。 到 2050 年, 超过 10 亿的亚洲居民将面临水源短缺。 到 2080 年, 水源短缺将威 胁11 亿 ~32 亿人。

瞯到 2030 年, 全球变暖所造成的穷人死亡率将升高。 疟疾、 登革热以及 食用污 染海 产品所 造成 疾病的 发生 率上升。

1畅2 问题的背景资料 5

(21)

瞯到 2050 年, 欧洲的小冰 川将会消失, 陆地大冰川将 会显著缩小。 到 2100 年, 欧洲超过一半的植物物种将会濒临灭绝。

瞯到 2080 年, 由于全球变暖而遭受饥荒的人数将在 2 亿 ~6 亿。

瞯到 2080 年, 随着海平面的上升, 每年遭受洪灾的人数约有 1 亿。

瞯美国城市的烟霾现象将会严重。 21 世纪中期由于气候变化导致的臭氧 污染将 比20 世纪 90 年代增加 4畅5%, 导致健康问题。 (以上数据来源于 2007 年IPCC 报告和 http://blog.66wz.com /user3 /cutycai/208185.html)

1畅2畅3 为什么是北极

毫 无疑问, 全球暖化将直接导致地球上冰雪 (如北极浮冰、 格陵兰冰盖、

南极冰 盖等) 的消融———这将促使海平面进一步升高。

那 么为 什么在 “Take a Bath” 中要求参赛队研究北极冰盖而不是南极冰盖 融化的 影响 呢? 毕竟南极冰盖无论从面积还是质量来说都远比北极的冰盖大。

原因如 下。

瞯全球温度普遍升高, 但是北半球较高纬度地区温度升幅较大。 在过去的 100 年中, 北极温度升高的速率几乎是全球平均速率的 2 倍。

瞯1978 年以 来的卫 星资料 显示, 北极 年平 均海冰 面积 已经 以每 10 年 2畅7% (2畅1% ~3畅3%) 的速率退缩, 夏季的海冰退缩率较大, 为每 10 年退 缩7畅4% (5% ~9畅8%)。

瞯自 1900 年以来, 北半球季节性冻土的面积减少了大约 7%, 春季冻土 面积的 减幅 高达15%。 自 20 世纪 80 年代以来, 北极多年冻土层上层温度普 遍升高, 达到 3 ℃。

瞯预估在 2100 年之后, 格陵兰冰盖的退缩将会继续导致海平面上升。 如 果这种 地表 冰物质 的 负 平衡 ( 即冰物质损失的速率比由于降水增加而使冰物 质增加 的速 率更 快) 持续数千年, 最终将会导致格陵兰冰盖完全消融, 并导 致海平 面上 升约7 m。

瞯美国宇航局 2007 年 12 月卫星资料显示, 北极冰层的厚度比以前减少了 23%; 航海资料则显示北极冰层比 20 世纪 50 年代减少了 50%; 格陵兰表面 冰层融 化的 速度是15 年前的 4 倍。

如 图1-3 所示为从 1979—2007 年北极冰层的演化。

与 此同 时, IPCC 在 2007 年的报告中指出: “ 根据当前的全球模式研究,

预估南 极冰 盖将会 维持 在非常 寒冷 的状态, 不至于会出现表层大范围融化。 而 且由于 降雪 增 加, 冰质量会增加。” 因此, 北极冰雪融化所导致的一 系列环 境、 生态等问题将是今后一段时间关注的焦点之一———这可能也是数学及其应 用联合 会 (The Consortium for Mathematics and its Applications, COMAP) 将北 极冰雪 融化 所造成 的影 响作为2008 年国际数学建模竞赛题的原因吧。

6 第1 章 冰盖融化问题

(22)

图1-3 1979—2007 年北极冰层的演化: 北极的冰雪覆盖明显在缩小

(图 片 来 源 自 http: //www畅5 joys畅com /better/sos2010 /)

1畅3 问题的数学模型与结果分析

“ Take a Bath” 是个非常开放的问题。 北极冰雪融化所带来的影响有很多,

可以研 究其 对生态 生物 系统的 影响、 对气候系统的影响、 对人类社会生活的影 响、 对社会经济活动的影响, 等等。 而且数学模型的选取更是多种多样, 所获 得的结 论也 可以是 各不 相同的。

本 书选 取了两 个模 型: 多元线性回归模型和初等数学模型, 下面分别进行 介绍。 这 两 个 模 型 来 自 于 两 份 非 常 优 秀 的 答 卷, 一 份 是 当 年 的 特 等 奖

(Outstanding Winner) 的论 文, 另一份是当 年的一等奖 ( Meritorious Winner)

的论文。 希望通过对这两个模型的介绍, 能够引发读者对该问题的研究兴趣和 对环境 问题 的深度 思考。

1畅3畅1 模型一: 多元线性回归模型

本 小节 将介绍 如何 运用多 元 线 性回 归 模型[3]来 预 测北 极 冰 雪 融化 所 带 来 的可能 影响。 首先, 收集了大量数据, 如全球温度数据、 温室气体排放数据、

佛罗里 达州 附近 海平 面 数 据 (佛罗里达东海岸的圣彼得斯堡和彭萨科拉, 佛 罗里达 西海 岸的 费南 迪 纳 比 奇和 弗 吉 尼亚 岛) 等, 并对数据进行了必要的整 理和分 析。 其次, 分别建立模型对格陵兰岛冰盖消融体积和北极海冰消融体积 进行了 计算。 最后, 建立了多元线性回归模型对未来 50 年的全球平均气温、

全球CO排 放量以 及平 均海平 面上 升幅度 进行 了预测。 通过求解模型, 得到了 未来50 年北极冰雪融化 (主要是格陵兰冰盖和北极海冰的融化) 对佛罗里达 海岸区 域的 影响, 并特别讨论了由于海平面上升所导致的海岸线侵蚀问题。

一、 数据的收集和处理 1畅全球温度数据

由 于建 立模型 的需 要, 收集了 1950—2001 年的全球温度数据, 所有数据 都来 源于 地 球政 策 研究 所 网站 ( http://www畅earth-policy.org /)。 从图 1 -4 可

1畅3 问题的数学模型与结果分析 7

(23)

以看出, 从 1950—1973 年, 全球平均气温总是在 14 ℃上下波动, 变化非常平 缓。 但是从 1973—2007 年, 全球平均气温上升了大约 1 ℃, 接近 15 ℃。 因 此, 全球变暖的趋势非常明显。

图1-4 从 1950 年到2001 年的全球温度演化图

( 图 片 来 源 文献 [ 3] )

2畅温室气体排放数据

温室 气体的排放被公认 为是 导致 全球 变暖 的主 要因 素之 一, 图 1-5 是从美 国能 源 部 二 氧 化 碳 信 息 分 析 中 心 (Carbon Dioxide Information Analysis Center,

CDIAC) 收集的温室气体数据。

图1-5 不同温室气体对热量的吸收率以及在大气中的成分分布

( 图 片 来 源 文献 [ 3] )

从 图1-6 可以看出, 不管是对热量的吸收率还是在大气中的含量, CO都 是温室 气体 的主要 成分。 图 1-6 给出了 1950—2001 年 CO的 排放 量。

3畅佛罗里达州附近海平面数据

选 取佛 罗里达 州四 个沿海 城市 作为观 测点, 它们分别是佛罗里达东海岸的 圣彼得 斯堡 和彭萨 科 拉, 佛罗里达西海岸的费南迪纳比奇和弗吉尼亚岛。 图

8 第1 章 冰盖融化问题

(24)

1-7显示了所选取的四个城市。 图 1-8 则给出了所选择的四个城市附近海平面 的年平 均数 据。

图1-6 从1950 年到 2001 年 CO的排放量

( 图 片 来 源 文献 [ 3] )

图1-7 所选取的四个城市

( 图 片 来 源 文献 [ 3] )

1畅3 问题的数学模型与结果分析 9

(25)

图1-8 所选取的四个城市附近海平面的数据

( 图 片 来 源 文献 [ 3] )

二、 冰雪消融体积的计算

随 着温 度的升 高, 北极冰雪将从低纬度到高纬度不断融化。 北极冰雪的构 成主要 包括 格陵兰 冰盖、 北极海冰以及其他一些较小的冰河和冻土。 作为北极 冰雪的 主体, 格陵兰冰盖和北极海冰的融化将对全球海平面的上升起到关键作 用。 如何合理计算格陵兰冰盖和北极海冰的融化量将直接关系到整个解决方案 的成败。 在原答卷中, 作者们根据阿基米德定律 ( Archimedes law, 又称为阿 基米德 原理) 和冰水体积关系分别建立模型估计了格陵兰冰盖和北极海冰的 融化量。

( 一) 格陵兰冰盖消融体积的计算 1畅基本原理

格 陵兰 冰盖和 北极 海冰的 主要 成分是 不同 的, 格陵兰冰盖的主要成分是淡 水, 密度是 900 kg /m; 北极海冰的主要成分是海水, 密度是 910 ~915 kg /m。 全球的 不断 暖化将 导致 这两类 冰雪 的显著 消融, 并都将导致海平面的升高。 根 据阿基 米德 定律, 有

F=ρgV (1畅1)

0

1 第1 章 冰盖融化问题

(26)

式中, F表示 液体 的 浮 力, ρ表 示 液 体的 密 度, g 是重力加速度, V 是 溢 出 的液体 的体 积。

因 为格 陵兰冰 盖覆 盖在大 陆上, 所以当冰雪消融时, 融化的淡水将直接注 入海水 中, 由于冰雪和淡水的密度差异, 1 个体积单位的冰雪融化将得到 0畅9 个体积 单位 的淡水, 即

ΔVw a te r=0畅9ΔVic e (1畅2)

忽略淡 水和 海水在 密度 上的微 小差 别, 则

ΔVse a=ΔVw ate r (1畅3)

北 极海 冰的计 算将 复杂许 多, 因为海冰是浮在海面上的, 90%的体积沉浸 在水面 以下, 因此当海冰融化的时候, 并不是简单地增加海水总量进而抬升海 平面。

海 冰是 漂浮在 海面 上的, 根据阿基米德定律, 重力和浮力是相等的, 故有 F=Gsi (1畅4)

式中, Gsi是 海冰 的重量, 它可以进一步表示为

Gsi=ρsigV (1畅5)

这样, 海冰融化对海平面上升的贡献为

ΔVs i=ΔVsiwa ter-ΔV (1畅6)

式中, ΔVsiw ate r是海冰 融化 带来的 海水 增量, 可以根据式 (1畅2) 计算出来。

代 入ρsi=915 kg /m及ρse a=1 030 kg /m, 并应用式 (1畅1)、 式 (1畅2) 可 以得到

ΔVs ea≈0畅026 7ΔVs i (1畅7)

式 (1畅7) 表明一个体积单位的海 冰融化将带来 0畅026 7 个体积单位 的海水 增量。

显 然, 根据式 (1畅2) 和式 (1畅7), 只要知道了格陵兰冰盖的冰雪融化量 和北极 海冰 的冰雪 融化 量就可 以大 致计 算 出 北极 冰 雪融 化 所带 来 的海 水 增量, 进而计 算出 平均海 平面 的上升 幅度。

2畅格陵兰冰盖的年融化量计算

全 球气 温的不 断升 高, 造成格陵兰冰盖的不断融化。 在原答卷中, 作者们 收集了1992—2003 年共 11 个年份由 ERS-1 /ERS-2 卫星发布的数据 ( 卫星图 像)。 如图 1-9 所示为格陵兰冰盖在 1992 年和 1993 年夏季的卫星图像, 不同 的深浅 代表 了不同 的冰 雪厚度。

根 据格 陵兰冰 盖 的 卫 星图 像, 可以估算出每个单元 ( 色块) 上的冰盖变 化量 (冰雪融化量), 因此, 整个格陵兰岛冰盖的年变化量就可以由下式计算 出来

ΔVic eg

i =0n Si· Δhi (1畅8)

1 1畅3 问题的数学模型与结果分析 1

(27)

图1-9 格陵兰冰盖在 1992 年和1993 年夏季的卫星图像

( 图 片 来 源 文献 [ 3] )

式中, Si 表示第i 个单元所表示的实际地表面积, Δhi 表 示第 i 个单元区域的 冰雪厚 度变 化量, n 是整个格陵兰冰盖上总的单元数目。

显 然每 个Si对 应唯 一的一 个四 元数组 (u, v, s, t), 如图 1-10 所示。

图1-10 单元与四元数组 ( u, v, s, t) 的对应关系

( 图 片 来 源 文献 [ 3] )

利 用曲 面积分 可以 比较精 确 地 计算 出Si 的 面 积, 进而计算出该单元冰雪 的融化 量, 但是计算过程过于复杂。 简单的做法是, 根据问题的实际背景, 将 这样 一 个 不 太 规 则 的 曲 面 近 似 为 一 个 长 方 形, 这 样 其 面 积 在 t-s → 0,

u-v →0时就可以近似地用下式计算 Si≈ t-s360 · u-v

180 · cos u+cos v

2 · L· L (1畅9)

式中, L 和L 分别表 示Si 所对应 的纬 度和经 度的 跨度。

根 据从 英 国气 象 局 哈 德 利 中 心 ( Met Office Hadley Center ) 获得的数据,

可以得 到在每 个1×1 区域 ( 这里的 “1” 分别表示 1 度跨度的纬度和 1 度跨度 2

1 第1 章 冰盖融化问题

(28)

的经度) 的冰雪覆盖率 (见图 1-11) 。 假设每个 1 ×1 区域上的冰雪覆盖率为 ηi, 则格陵兰冰盖的年变化量可以精确表示为

ΔVic e g

i =0n Si· Δhi· ηi (1畅10)

图1-11 1×1 区域上的冰雪覆盖率

( 图 片 来 源 文献 [ 3] )

但 是, 英国气象局哈德利中心 (Met Office Hadley Center) 的数据并不完 善, 在某些区域上的数据存在缺失, 这就要求用户能够合理地估计、 修复这些 缺失的 数据。 通过比较格陵兰岛的等高线图和上述卫星图像, 可以发现格陵兰 岛的冰 雪覆盖 规律 (范围、 厚度等) 与该岛的地形分布大致相符, 因此根据 该岛的 地形 可以设 计一 种简单 有 效的 方 法来 估 计、 修复缺失的数据 (见图 1-

12) 。 具体步骤如下。

第 一步, 根据抛物线拟合 -42°经度线上 1×1 区域的数据, 该抛物线可以

根据 (lno rth, 0), (0, hpea k), ( lso uth, m) 三个点确定。

y =-hpe ak

lno rth· x+hpe a k (1畅11)

式中, lno rth是 图1-12 (b) 中从 O 到 A 的距离; hpe ak是格 陵 兰 岛 的 海 拔 高度;

lso uth是 从O 到 B 的距离, O、 A、 B 的位置如图 1-12 (b) 所示。

第 二步, 采用与第一步相同的步骤拟合出不同纬度的抛物线。

经 过对 缺失数 据 的修 复, 再应用式 (1畅10) 就可以计算出格陵兰岛每年 大约有5畅171×1011 的冰 雪融化 成淡 水并注 入太 平洋。

3 1畅3 问题的数学模型与结果分析 1

(29)

图1-12 格陵兰岛的地形图和解决方法

( 图 片 来 源 文献 [ 3] )

( 二) 北极海冰消融体积的计算 1畅北极海冰数据

北 极海 冰数据 可以 从 美 国国 家 冰 雪 数 据 中心 ( National Snow and Ice Data Center, NSIDC) 获得。 图 1-13 显示了从 1979—2007 年北极海冰体积的演化,

显然海 冰的 区域和 体积 在不断 缩减。

图1-13 北极海冰变化曲线

( 图 片 来 源 文献 [ 3] )

4

1 第1 章 冰盖融化问题

(30)

2畅面积-厚度关系模型

如 图1-14 所示为北极海冰的面积和厚度分布。 从图中不难发现, 海冰的 面积在 迅速 缩小, 越靠近加拿大和格陵兰岛的海冰厚度越大。

图1-14 北极海冰的分布

( 图 片 来 源 文献 [ 3] )

通 过观 察北极 海冰 的厚度 和面 积, 可以发现北极海冰的分布类似于一个半 圆, 当温度下降时北极海冰的面积和厚度将同步增加。 这里建立了一个几何模 型来估 算融 化的北 极海 冰体积。 图 1-15 是海冰体积变化示意。 海冰的融化体 积等于 上下 两个抛 物面 所界定 部分 的体积, 可以利用积分计算出来。

图1-15 北极海冰体积变化示意

( 图 片 来 源 文献 [ 3] )

假 设在 较低温 度下, 新扩展的海冰面积是 S, 在较高温度下, 缩减的海 冰面积 是 S。 根据前面的假设, 海冰的范围是一个半圆, 那么该半径可以计 算如下

R= 2S /π (1畅12)

冰体侧 面轮 廓形状 由 (0, h) , (±R) 确定

5 1畅3 问题的数学模型与结果分析 1

(31)

y =-hR· x+h (1畅13)

融化的 海冰 体积可 以通 过积分 计算

ΔVice s=π·

hh-yh · Rdy (1畅14)

根据式 (1畅14), 北极海冰总的融化体积是

ΔVic es=πh2 · (R-R) (1畅15)

式中, h 是海冰最厚处的厚度, 约等于 5畅5 m。

最 后, 应用式 (1畅15) 可以计算出北极海冰融化将增加 2畅6 ×10 的 水 到海洋 中, 与格陵兰岛相比, 这个数值无疑是非常小的。

三、 多元线性回归模型

基 于前 面收集 的数 据和已 有结 果, 这里将建立多元线性回归模型来预测由 于北极 冰雪 融化 (主要是格陵兰冰盖和北极海冰的融化) 所导致的全球海平 面上升 幅度。 如图 1-16 所示为使用多元线性回归模型对海平面上升幅度进行 预测的 流程。

图1-16 使用多元线性回归模型对海平面上升幅度进行预测的流程

( 图 片 来 源 文献 [ 3] )

众 所周 知, CO等 温室气 体的 排放以 及冰 雪消融 对全 球气温 的升 高 具有 重 要影响。 同时, 全球气温的升高和冰雪消融对海平面的升高也具有重要的影 响。 因此, 首先需要根据因素之间的相关性选取对海平面上升影响最显著的 因素。

6

1 第1 章 冰盖融化问题

(32)

瞯温室气体和全球气温变化之间的相关性。

瞯北极冰雪融化和全球气温之间的相关性。

瞯北极冰雪融化和海平面变化之间的相关性。

瞯全球气温变化和海平面变化之间的相关性。

相 关性 分析的 结果 见表1-1。

表1-1 相关性分析

因素 温室气体和温度 冰雪融化量和温度 冰雪融化量和海平面 温度和海平面

相关系数 0畅762 0畅481 0畅793 0畅553

显 然, 全球气温直接受到温室气体排放和北极冰雪融化的影响。 因此, 这 里首先 建立 对全球 温度 变化进 行预 测的多 元线 性回归 模型

ΔT =α+α· ΔVic e+α· ΔMC O2+ε (1畅16)

式中, ΔVice表示 北极冰 雪融 化量, ΔMC O2表 示CO 的排放 量。

全 球海 平面的 上升 将受到 北极 冰雪融 化和 全球气 温变化 的直 接影响。 类似 地, 可以建立多元线性回归模型

Δhse a=β+β· ΔVice+β· ΔT+ε (1畅17)

为 了确 定上述 多元 线性回 归模 型的系 数, 这里选取了佛罗里达州四个城市

(佛罗里达东海岸的圣彼得斯堡和彭萨科拉, 佛罗里达西海岸的费南迪纳比奇 和弗吉 尼亚 岛) 附近的海面的观测数据。 回归方程 (1畅16) 和 (1畅17) 的系 数可以 使用SPSS 软件确定

α=14畅269, α=0畅436, α=-0畅254 以及

β=0畅993, β=0畅021, β=0畅008

显 著性检 验的 结果表 明, α和α、 β 和β 都是显 著的, 而且线性关系也 比较显 著———这也说明选择多元线性回归模型对全球平均气温和全球海平面上 升幅度 进行 预测是 比较 合理的。

四、 模型的求解和结果分析

根 据收 集的数 据, 求解建立的模型, 可以预测未来 50 年中全球海平面的 上升幅 度、 佛罗里达附近海平面的上升幅度以及对海岸线的侵蚀。

( 一) 全球海平面的上升幅度

求 解预 测 模 型 (1畅16 ), 得到全球气温 在未来 50 年的变化预 测值 ( 图 1-17)。  

7 1畅3 问题的数学模型与结果分析 1

(33)

图1-17 全球气温预测

( 图 片 来 源 文献 [ 3] )

  IPCC 报告中全球未来 50 年的温度预测见表 1-2。

表1-2 IPCC 报告中全球未来 50 年的温度预测

年份 2010 2020 2030 2040 2050

温度℃ 15畅43 15畅56 15畅68 15畅99 16畅28 使 用预 测模型 (1畅17) 对未来 50 年由于北极冰雪融化所导致的海平面上 升幅度 进行 了预测, 具体结果如图 1-18 所示。

图1-18 北极冰雪融化所导致的海平面上升幅度

( 图 片 来 源 文献 [ 3] )

根 据上 面的计 算结 果不 难发 现, 在未来 50 年中, 北极冰雪融化将导致海 平面显 著上 升。

8

1 第1 章 冰盖融化问题

參考文獻

相關文件

推 荐 期 刊 : Journal of Cross-cultural Psychology, Journal of International Business Studies, Management and Organization Review

[r]

李胜宏 第44届IMO中国队领队、浙江大学教授、博士生导师 李伟固 国家集训队教练、北京大学教授、博士生导师 刘鸿坤 第31、32届IMO中国队副领队、华东师范大学教授 刘诗雄

微积分的创立是数学发展中的里程碑, 它的发展 和广泛应用开启了向近代数学过渡的新时期, 为研究 变量和函数提供了重要的方法和手段. 运动物体的瞬

为此, 我们需要建立函 数的差商与函数的导数间的基本关系式, 这些关系式称为“微分学中值定理”...

在科学史上 , 第一个发表能量转化与守恒原理的是德国青年 医生迈尔。迈尔的研究是从生命运动开始的。 1840 年,迈尔随船从

2.

与大学生一起做研究不是一个受关注的话题,因 为不少资深教育家和德高望重的老师不认为在校