拆分方程建模分析

1、拆分方程建模分析第1页，共173页，2022年，5月20日，0点33分，星期四1、什么是数学模型？ 数学模型是对于现实世界的一个特定对象，一个特定目的，根据特有的内在规律，做出一些必要的假设，运用适当的数学工具，得到一个数学结构。 简单地说：就是系统的某种特征的本质的数学表达式（或是用数学术语对部分现实世界的描述），即用数学式子（如函数、图形、代数方程、微分方程、积分方程、差分方程等）来描述（表述、模拟）所研究的客观对象或系统在某一方面的存在规律。第2页，共173页，2022年，5月20日，0点33分，星期四2、什么是数学建模? 数学建模是利用数学方法解决实际问题的一种实践。即通过抽象、简化、

2、假设、引进变量等处理过程后，将实际问题用数学方式表达，建立起数学模型，然后运用先进的数学方法及计算机技术进行求解。 观点：“所谓高科技就是一种数学技术”第3页，共173页，2022年，5月20日，0点33分，星期四 数学建模其实并不是什么新东西，可以说有了数学并需要用数学去解决实际问题，就一定要用数学的语言、方法去近似地刻划该实际问题，这种刻划的数学表述的就是一个数学模型，其过程就是数学建模的过程。数学模型一经提出，就要用一定的技术手段（计算、证明等）来求解并验证，其中大量的计算往往是必不可少的，高性能的计算机的出现使数学建模这一方法如虎添翼似的得到了飞速的发展，掀起一个高潮。 数学建模将各种

3、知识综合应用于解决实际题中，是培养和提高同学们应用所学知识分析问题、解决问题的能力的必备手段之一。第4页，共173页，2022年，5月20日，0点33分，星期四 在实际过程中用那一种方法建模主要是根据我们对研究对象的了解程度和建模目的来决定。机理分析法建模的具体步骤大致可见右图。符合实际不符合实际交付使用，从而可产生经济、社会效益实际问题抽象、简化、假设 确定变量、参数建立数学模型并数学、数值地求解、确定参数用实际问题的实测数据等来检验该数学模型建模过程示意图第5页，共173页，2022年，5月20日，0点33分，星期四 模型 数学模型的分类： 按研究方法和对象的数学特征分：初等模型、几何模型

4、、优化模型、微分方程模型、图论模型、逻辑模型、稳定性模型、扩散模型等。 按研究对象的实际领域（或所属学科）分：人口模型、交通模型、环境模型、生态模型、生理模型、城镇规划模型、水资源模型、污染模型、经济模型、社会模型等。三、数学模型及其分类第6页，共173页，2022年，5月20日，0点33分，星期四正是由于认识到培养应用型、研究型科技人才的重要性，而传统的数学竞赛不能担当这个任务，从1983年起，美国就有一些有识之士探讨组织一项应用数学方面的竞赛的可能性。经过论证、争论、争取资金等过程，1985年举行了美国第一届大学生数学建模竞赛。简称MCM竞赛由美国工业与用数学学会和美国运筹学学会联合主办。

5、第7页，共173页，2022年，5月20日，0点33分，星期四从1985年起每年举行一届，时间定为每年的二月下旬或三月初的星期五到星期日举行。 这项竞赛的宗旨是鼓励学生运用所学的知识（数学及其各门科学的知识）去参与解决实际问题的全过程。这些实际问题并不限于某个固定领域，可以涉及非常广泛的、并不固定的范围和领域。第8页，共173页，2022年，5月20日，0点33分，星期四 美国的MCM虽然只是美国的国内 赛，但它欢迎其他国家的大学组队参加，而且越来越多国家的大学参加这一竞赛，因此，在某种意义上它已经是国际比赛。我国最早由北京三所大学组队参加美国的MCM竞赛，继后我国参加此项比赛的大学越来越多。

6、第9页，共173页，2022年，5月20日，0点33分，星期四内容 赛题：工程、管理中经过简化的实际问题 答卷：一篇包含问题分析、模型假设、建立、求解(通常用计算机)、结果分析和检验等的论文形式 3名大学生组队，在3天内完成的通讯比赛 可使用任何“死”材料(图书/互联网/软件等), 但不得与队外任何人讨论（包括上网讨论）宗旨创新意识 团队精神 重在参与 公平竞争标准假设的合理性，建模的创造性，结果的正确性，表述的清晰性。数学建模竞赛内容与形式第10页，共173页，2022年，5月20日，0点33分，星期四我国大学生数学建模竞赛（CUMCM） 1992年中国工业与应用数学学会(CSIAM)开始组

7、织 1994年起教育部高教司和CSIAM共同举办(每年9月) 2008年有31省(市、区）的1022所学校12836队参加 网址： 奖励：全国一等奖（约2%）、全国二等奖（约7%）教育部高教司和CSIAM共同签章 1999年起竞赛分为甲组(本科)、乙组(高职高专组) 优秀论文刊登于次年工程数学学报（ 2000年前为数学的实践与认识）2009 年全国有33个省/市/自治区(包括香港和澳门特区)1137所院校、15046个队（其中甲组12276队、乙组2770队）、4万5千多名来自各个专业的大学生参加竞赛 (其中西藏和澳门是首次参赛)！第11页，共173页，2022年，5月20日，0点33分，星期

8、四 怎样撰写数学建模的论文？1、摘要:问题、模型、方法、结果2、问题重述4、分析与建立模型5、模型求解6、模型检验7、模型推广8、参考文献9、附录实例3、模型假设 返回实例参考解答第12页，共173页，2022年，5月20日，0点33分，星期四差分方程建模处理动态的离散型的问题处理对象虽然涉及的变量(如时间)是连续的，但是从建模的目的考虑，把连续变量离散化更为合适，将连续变量作离散化处理，从而将连续模型(微分方程)化为离散型(差分方程)问题 第13页，共173页，2022年，5月20日，0点33分，星期四第一讲 差分方程基础知识设一阶差分：一阶后向差分一阶前向差分注：差分是导数的一种推广高阶差

9、分：第14页，共173页，2022年，5月20日，0点33分，星期四第一讲 差分方程基础知识性质1(可加性)：微分的性质与运算：性质2(齐次性)：性质3(线性性)：性质3(线性性)：性质4(乘积)：性质4(除法运算)：第15页，共173页，2022年，5月20日，0点33分，星期四第一讲 差分方程基础知识性质1(可加性)：差分的性质与运算：性质2(齐次性)：性质3(线性性)：性质3(线性性)：性质4(乘积)：性质4(除法运算)：注：注意它与微分性质的比较第16页，共173页，2022年，5月20日，0点33分，星期四第一讲 差分方程基础知识差分的性质与运算：第17页，共173页，2022年，5

10、月20日，0点33分，星期四第一讲 差分方程基础知识差分方程及其解：定义2：包含未知函数及差分的方程式称为差分方程注：差分方程的实质就是递推公式线性差分方程：差分方程的解：第18页，共173页，2022年，5月20日，0点33分，星期四常微分方程化为差分方程 用导数近似式替代导数或者说用适当近似式替代含有导数的表达式，可以得到这些近似值满足的代数方程-差分方程 以二阶常微分方程边值问题为例 目的求差分法第19页，共173页，2022年，5月20日，0点33分，星期四一般k阶常系数线性差分方程为差分方程边值问题第20页，共173页，2022年，5月20日，0点33分，星期四二 偏微分方程化为差分

11、方程以二阶椭圆方程的边值问题为例用两族平行坐标轴的直线 正方形网格把区域G剖分 第21页，共173页，2022年，5月20日，0点33分，星期四节点可分三类 1通过该节点的网格线上的相邻四网点都在G内,记 G12在G内部但不属于G1 ，记G23恰在边界上记G3 确定各节点处解的近似值uij，需要建立代数方程，每一节点建立一个代数方程任务第22页，共173页，2022年，5月20日，0点33分，星期四(i，j-1) (i，j+1) (i-1，j)(i，j)(i+1，j)偏导数近似式替代第23页，共173页，2022年，5月20日，0点33分，星期四差分方程 N (i，j) E第24页，共173页

12、，2022年，5月20日，0点33分，星期四 偏导数近似式替代第25页，共173页，2022年，5月20日，0点33分，星期四和分(反差分)周期函数：如果和分：若则称为的和分（反差分），记为或者第26页，共173页，2022年，5月20日，0点33分，星期四积分中的N-L公式和分中的N-L公式定理：若在上有定义，并且，则练习：计算第27页，共173页，2022年，5月20日，0点33分，星期四四 二阶常系数齐次差分方程求法 齐次差分方程 (1)特征方程有两个不相等实根 (2)特征方程有两个相等实根 第28页，共173页，2022年，5月20日，0点33分，星期四例：求解差分方程例：求解差分方程

13、初值问题解：特征方程为从而特征根为所以差分方程的解为第29页，共173页，2022年，5月20日，0点33分，星期四(3)特征方程有一对共轭复根 则为一对复值解.令第30页，共173页，2022年，5月20日，0点33分，星期四例：求解差分方程第31页，共173页，2022年，5月20日，0点33分，星期四非齐次差分方程非齐次方程的通解=非齐次的特解+齐次的通解解法：1.待定系数法2.常数变易法例：兔子问题 在一年的时间里，一对兔子能够生育出多少对兔子？第32页，共173页，2022年，5月20日，0点33分，星期四每对兔子每个月生育出新的一对兔子假设新的一对兔子在一个月之后具有生育能力其次这

14、些兔子都不死亡第n个月开始时兔子对数模型结果Fibonacci数列黄金分割比第33页，共173页，2022年，5月20日，0点33分，星期四 市场经济中的蛛网模型问 题供大于求现象商品数量与价格的振荡在什么条件下趋向稳定当不稳定时政府能采取什么干预手段使之稳定价格下降减少产量增加产量价格上涨供不应求描述商品数量与价格的变化规律数量与价格在振荡第34页，共173页，2022年，5月20日，0点33分，星期四蛛 网 模 型gx0y0P0fxy0 xk第k时段商品数量；yk第k时段商品价格消费者的需求关系生产者的供应关系减函数增函数供应函数需求函数f与g的交点P0(x0,y0) 平衡点一旦xk=x0

15、，则yk=y0, xk+1,xk+2,=x0, yk+1,yk+2, =y0 第35页，共173页，2022年，5月20日，0点33分，星期四xy0fgy0 x0P0设x1偏离x0 x1x2P2y1P1y2P3P4x3y3P0是稳定平衡点P1P2P3P4P0是不稳定平衡点xy0y0 x0P0fg曲线斜率蛛 网 模 型第36页，共173页，2022年，5月20日，0点33分，星期四在P0点附近用直线近似曲线P0稳定P0不稳定方 程 模 型方程模型与蛛网模型的一致第37页，共173页，2022年，5月20日，0点33分，星期四 商品数量减少1单位, 价格上涨幅度 价格上涨1单位, (下时段)供应的

16、增量考察 , 的含义 消费者对需求的敏感程度 生产者对价格的敏感程度小, 有利于经济稳定 小, 有利于经济稳定结果解释xk第k时段商品数量；yk第k时段商品价格经济稳定结果解释第38页，共173页，2022年，5月20日，0点33分，星期四经济不稳定时政府的干预办法1. 使 尽量小，如 =0 以行政手段控制价格不变2. 使 尽量小，如 =0靠经济实力控制数量不变xy0y0gfxy0 x0gf结果解释需求曲线变为水平供应曲线变为竖直第39页，共173页，2022年，5月20日，0点33分，星期四模型的推广 生产者根据当前时段和前一时段的价格决定下一时段的产量。生产者管理水平提高设供应函数为需求函

17、数不变二阶线性常系数差分方程x0为平衡点研究平衡点稳定，即k, xkx0的条件第40页，共173页，2022年，5月20日，0点33分，星期四方程通解(c1, c2由初始条件确定)1, 2特征根，即方程 的根 平衡点稳定，即k, xkx0的条件:平衡点稳定条件比原来的条件 放宽了模型的推广第41页，共173页，2022年，5月20日，0点33分，星期四 银行复利问题 背景所付利息一年内复合n次，即把一年分n个相等的时间段，而所付利息为每一时间段的未尾 .给出一个可以预测在任意给定时间的帐目余额 分析帐目余额与时间直接相关，而时间是离散的本期结束时的总存款等于前一时期余下的本利，及本利得到的利息

18、与第本期内新存入的存款之和 任何时候都可以存款第42页，共173页，2022年，5月20日，0点33分，星期四模型假设1. 储蓄的年利率为 r2. 任何时候都可以存款，但存款利息只从下一时期开始计算，如时间段开始第一天的存款即开始计算利息 t期结束时的总存款 记号第t期内的新存款 第43页，共173页，2022年，5月20日，0点33分，星期四模 型注：上式中n=2时，相应于半年的复利，而n=365则是相应于逐日计算的复利第44页，共173页，2022年，5月20日，0点33分，星期四 抵押贷款买房问题 背景 每户人家都希望有一套属于自己的住房，但又没有足够的资金一次买下。这就产生了贷款买房问

19、题。某新婚夫妇急需一套属于自己的住房。他们看到一则理想的房产广告：“名流花园之高尚住宅公寓，供工薪阶层选择。一次性付款优惠价40.2万元。若不能一次性付款也没关系，只付首期款为15万元，其余每月1977.04元等额偿还，15年还清。(公积金贷款月利息为3.675）。问题公寓原来价多少？每月等额付款如何算出来？第45页，共173页，2022年，5月20日，0点33分，星期四假设贷款期限内利率不变 银行利息按复利计算 记号A（元）：贷款额（本金） n（月）：货款期限r ：月利率B(元) :月均还款额 Ck：第k个月还款后的欠款第46页，共173页，2022年，5月20日，0点33分，星期四模型求解

20、代入n=180、 r=0.003675、 B=1977.04结果： A=260000（元）一次性优惠价9.8折还款总额 利息负担总额第47页，共173页，2022年，5月20日，0点33分，星期四 减肥计划节食与运动背景 多数减肥食品达不到减肥目标，或不能维持 通过控制饮食和适当的运动，在不伤害身体的前提下，达到减轻体重并维持下去的目标分析 体重变化由体内能量守恒破坏引起 饮食（吸收热量）引起体重增加 代谢和运动（消耗热量）引起体重减少 体重指数BMI=w(kg)/l2(m2). 18.5BMI25 超重; BMI30 肥胖.第48页，共173页，2022年，5月20日，0点33分，星期四模型

21、假设1）体重增加正比于吸收的热量每8000千卡增加体重1千克；2）代谢引起的体重减少正比于体重每周每公斤体重消耗200千卡 320千卡(因人而异), 相当于70千克的人每天消耗2000千卡 3200千卡；3）运动引起的体重减少正比于体重，且与运动形式有关； 4）为了安全与健康，每周体重减少不宜超过1.5千克，每周吸收热量不要小于10000千卡。第49页，共173页，2022年，5月20日，0点33分，星期四某甲体重100千克，目前每周吸收20000千卡热量，体重维持不变。现欲减肥至75千克。第一阶段：每周减肥1千克，每周吸收热量逐渐减少，直至达到下限（10000千卡）；第二阶段：每周吸收热量保

22、持下限，减肥达到目标 2）若要加快进程，第二阶段增加运动，试安排计划。1）在不运动的情况下安排一个两阶段计划。减肥计划3）给出达到目标后维持体重的方案。第50页，共173页，2022年，5月20日，0点33分，星期四 确定某甲的代谢消耗系数即每周每千克体重消耗 20000/100=200千卡基本模型w(k) 第k周(末)体重c(k) 第k周吸收热量 代谢消耗系数(因人而异)1）不运动情况的两阶段减肥计划每周吸收20000千卡 w=100千克不变第51页，共173页，2022年，5月20日，0点33分，星期四 第一阶段: w(k)每周减1千克, c(k)减至下限10000千卡第一阶段10周, 每

23、周减1千克，第10周末体重90千克吸收热量为1）不运动情况的两阶段减肥计划第52页，共173页，2022年，5月20日，0点33分，星期四 第二阶段：每周c(k)保持Cm, w(k)减至75千克 1）不运动情况的两阶段减肥计划基本模型第53页，共173页，2022年，5月20日，0点33分，星期四 第二阶段：每周c(k)保持Cm, w(k)减至75千克 第二阶段19周, 每周吸收热量保持10000千卡, 体重按 减少至75千克。第54页，共173页，2022年，5月20日，0点33分，星期四运动 t=24 (每周跳舞8小时或自行车10小时), 14周即可。2）第二阶段增加运动的减肥计划根据资料

24、每小时每千克体重消耗的热量 (千卡): 跑步 跳舞 乒乓 自行车(中速) 游泳(50米/分) 7.0 3.0 4.4 2.5 7.9t每周运动时间(小时)基本模型第55页，共173页，2022年，5月20日，0点33分，星期四3）达到目标体重75千克后维持不变的方案每周吸收热量c(k)保持某常数C，使体重w不变 不运动 运动(内容同前)第56页，共173页，2022年，5月20日，0点33分，星期四2007年全国大学生数学建模A题 中国是一个人口大国，人口问题始终是制约我国发展的关键因素之一。根据已有数据，运用数学建模的方法，对中国人口做出分析和预测是一个重要问题。 近年来中国的人口发展出现了

25、一些新的特点，例如，老龄化进程加速、出生人口性别比持续升高，以及乡村人口城镇化等因素，这些都影响着中国人口的增长。2007年初发布的国家人口发展战略研究报告(附录1) 还做出了进一步的分析。 关于中国人口问题已有多方面的研究，并积累了大量数据资料。附录2就是从中国人口统计年鉴上收集到的部分数据。 试从中国的实际情况和人口增长的上述特点出发，参考附录2中的相关数据（也可以搜索相关文献和补充新的数据），建立中国人口增长的数学模型，并由此对中国人口增长的中短期和长期趋势做出预测；特别要指出你们模型中的优点与不足之处。附录1 国家人口发展战略研究报告附录2人口数据（中国人口统计年鉴中的部分数据）及其说

26、明 第57页，共173页，2022年，5月20日，0点33分，星期四 人口预测是国家工作中的重点，关系着国家的发展方向和命运。我国是一个人口大国，人口问题始终是制约我国发展的关键因素之一。根据已有数据，运用数学建模的方法，对我国人口做出分析和预测是一个重要问题。 近年来我国的人口发展出现了一些新的特点，例如，老龄化进程加速、出生人口性别比持续升高，以及乡村人口城镇化等因素，这些都影响着我国人口的增长。2007年初发布的国家人口发展战略研究报告对此做出了进一步的分析。 从我国的实际情况和人口增长的上述特点出发，参考相关数据，建立我国人口增长的数学模型，并由此对我国人口增长的中短期和长期趋势做出预

27、测；特别注意指出模型中的优点与不足之处。 1. 问题重述 第58页，共173页，2022年，5月20日，0点33分，星期四2. 问题分析 一个社会（国家、省市、地区）人口的变化和随时间的发展过程，是由很多因素决定的，社会制度、自然环境、生活水平、科学文化水平、战争、自然灾害和移民等等，都能严重地影响社会人口的发展过程。然而，婴儿的出生、人口的死亡、居民的迁移却是决定该社会人口变化的直接原因，近年来我国人口发展出现的一些新特点，如老龄化进程加速、出生人口性别比持续升高，以及乡村人口城镇化等因素，都直接或间接地通过这三个现象表现出来。综合考虑这些因素成为构建符合我国国情的人口增长模型关键。 建立模

28、型对人口发展过程进行定量预测，就是根据现有的人口统计资料和原始数据，从当前实际的人口状况出发，并对未来的人口发展过程，提出合理的控制要求和假定，应用科学的方法，预测出未来几年、几十年甚至上百年的人口发展趋势，包括人口总数、人口的性别、年龄和城乡结构，人口出生、死亡和自然增长率的变化以及在未来的人口构成中劳力和抚养水平及老化水平等。 第59页，共173页，2022年，5月20日，0点33分，星期四3. 模型假设针对本问题，建立如下合理的假设： 题中所给数据能反映我国人口变化的基本情况； 一些重大事件，如战争、自然灾害等对人口预测的影响暂不考虑； 所给数据都是年末数据，也即下一年年初数据，如200

29、1年总人口实质上也表示2002年初的总人口； 今年所统计的i岁的人口在下一年年初均为 i+1岁； 生育模式不随时间变化。 第60页，共173页，2022年，5月20日，0点33分，星期四4. 符号说明 第61页，共173页，2022年，5月20日，0点33分，星期四5. 模型建立与求解 5.1. 数据预处理 题中所给 5 年我国人口 1%调查数据是对人口的抽样调查数据，由于数据的不完备性，并不能由它来估计当时的全国总人口数。但基于抽样调查的等概率性，可以认为它所反应的市、镇、乡三个地区的人口比例及男女比例是与实际较为接近的。从中国人口统计年鉴20061, 可以得到20012005 年具体的全国

30、总人口数。进而可以得到各部分人口数。所得数据见表1。 例：第62页，共173页，2022年，5月20日，0点33分，星期四第63页，共173页，2022年，5月20日，0点33分，星期四5.2. 模型一：基于人口迁移的 Logistic阻滞增长模型 第64页，共173页，2022年，5月20日，0点33分，星期四第65页，共173页，2022年，5月20日，0点33分，星期四第66页，共173页，2022年，5月20日，0点33分，星期四第67页，共173页，2022年，5月20日，0点33分，星期四第68页，共173页，2022年，5月20日，0点33分，星期四第69页，共173页，2022

31、年，5月20日，0点33分，星期四第70页，共173页，2022年，5月20日，0点33分，星期四第71页，共173页，2022年，5月20日，0点33分，星期四5.3. 模型二：离散人口发展方程模型 第72页，共173页，2022年，5月20日，0点33分，星期四第73页，共173页，2022年，5月20日，0点33分，星期四第74页，共173页，2022年，5月20日，0点33分，星期四第75页，共173页，2022年，5月20日，0点33分，星期四第76页，共173页，2022年，5月20日，0点33分，星期四第77页，共173页，2022年，5月20日，0点33分，星期四第78页，共1

32、73页，2022年，5月20日，0点33分，星期四第79页，共173页，2022年，5月20日，0点33分，星期四第80页，共173页，2022年，5月20日，0点33分，星期四第81页，共173页，2022年，5月20日，0点33分，星期四第82页，共173页，2022年，5月20日，0点33分，星期四第83页，共173页，2022年，5月20日，0点33分，星期四第84页，共173页，2022年，5月20日，0点33分，星期四第85页，共173页，2022年，5月20日，0点33分，星期四第86页，共173页，2022年，5月20日，0点33分，星期四第87页，共173页，2022年，5月

33、20日，0点33分，星期四第88页，共173页，2022年，5月20日，0点33分，星期四第89页，共173页，2022年，5月20日，0点33分，星期四第90页，共173页，2022年，5月20日，0点33分，星期四第91页，共173页，2022年，5月20日，0点33分，星期四第92页，共173页，2022年，5月20日，0点33分，星期四第93页，共173页，2022年，5月20日，0点33分，星期四第94页，共173页，2022年，5月20日，0点33分，星期四第95页，共173页，2022年，5月20日，0点33分，星期四 Malthus 模型 设xn是某人类群体在第n个时间段(例如

34、年)末时的总数，若在单位时间段内人口相对增长率为r（出生率与死亡率之差）,那么人口增长数与原人口数成正比,从而xn+1 xn r xn即 xn+1 = a xn其中 a=r+1. 差分形式的人口增长模型第96页，共173页，2022年，5月20日，0点33分，星期四这是一个如下线性映射的迭代 f (x) = a x从而 xn= a xn1= a2xn2 = an x0 Malthus的结论：人口增长呈几何级数 约35年增加一倍，与17001961年世界人口统计结果一致 与近年统计结果有误差，由a 1,xn趋向无穷，模型在人口长期预测方面必定是失效的. 差分形式的人口增长模型第97页，共173页

35、，2022年，5月20日，0点33分，星期四 离散形式的阻滞增长模型连续形式的阻滞增长模型 (Logistic模型)t, xN, x=N是稳定平衡点(与r大小无关)离散形式x(t) 某种群 t 时刻的数量(人口)yk 某种群第k代的数量(人口)若yk=N, 则yk+1,yk+2,=N讨论平衡点的稳定性，即k, ykN ？y*=N 是平衡点第98页，共173页，2022年，5月20日，0点33分，星期四离散形式阻滞增长模型的平衡点及其稳定性一阶(非线性)差分方程 (1)的平衡点y*=N讨论 x* 的稳定性变量代换(2)的平衡点第99页，共173页，2022年，5月20日，0点33分，星期四(1)

36、的平衡点 x*代数方程 x=f(x)的根稳定性判断(1)的近似线性方程x*也是(2)的平衡点x*是(2)和(1)的稳定平衡点x*是(2)和(1)的不稳定平衡点补充知识一阶非线性差分方程的平衡点及稳定性第100页，共173页，2022年，5月20日，0点33分，星期四01的平衡点及其稳定性平衡点稳定性x* 稳定x* 不稳定另一平衡点为 x=0不稳定第101页，共173页，2022年，5月20日，0点33分，星期四01/2101的平衡点及其稳定性第102页，共173页，2022年，5月20日，0点33分，星期四初值 x0=0.2数值计算结果b 3.57, 不存在任何收敛子序列混沌现象4倍周期收敛第

37、106页，共173页，2022年，5月20日，0点33分，星期四的收敛、分岔及混沌现象b第107页，共173页，2022年，5月20日，0点33分，星期四2浑沌与遍历性 当c*a4时，Logistic映射进入浑沌区域.反映出的是： 遍历性：点 x0的轨道不趋向任何稳定的周期轨道, 它的轨道在(0,1)(或其中某些区间)内的任何一个子区间(a,b)内都会出现无数次.这是浑沌的 敏感性： 轨道表现出对初始条件的强烈敏感性,即不同初始值，即使它们离得非常近，它们的轨道也终将以某种方式分离. 存在周期小窗口 浑沌区域内某些地方仍有倍周期分叉，例如a3.835附近第108页，共173页，2022年，5月

38、20日，0点33分，星期四 Feigenbaum常数 比值(ck-ck-1)/(ck+1-ck)在k 趋于无穷时，趋于常数 q =4.6692016这常数的意义在于普适性，例如周期3窗口也适用，还适用其他映射任务：验证遍历性、敏感性 周期3窗口的分叉、(结合Feigenbaum常数 )五. 图象方法 蛛网迭代 在以xn为横坐标、xn+1为纵坐标的第一象限作抛物线弧： xn+1a xn(1- xn)第109页，共173页，2022年，5月20日，0点33分，星期四 作图的过程第110页，共173页，2022年，5月20日，0点33分，星期四 任取(0,1)中的点x0,可以通过作图来取得迭代的数值

39、序列xn,从而也通过图象直观地看出由 x0出发的轨道的变化. 这作图的过程颇象蜘蛛织网,故称为蛛网迭代. 第111页，共173页，2022年，5月20日，0点33分，星期四 1a3 从(0,1)中任何初值出发的轨道趋向不动点 (周期1点)第112页，共173页，2022年，5月20日，0点33分，星期四 3a61/2+1 从任何初值出发的轨道趋向周期2点第113页，共173页，2022年，5月20日，0点33分，星期四61/2+1a0，即x(t)单调增加。由x(t0)=0，可以得出 =1，此时， 。当t0，x(t)单调增加，而当tt0时，x(t)0，x(t)单调减小。实际调查表明，销售曲线与L

40、ogistic曲线十分接近，尤其是在销售后期，两者几乎完全吻合。 在销出量小于最大需求量的一半时，销售速度是不断增大的，销出量达到最大需求量的一半时，该产品最为畅销，接着销售速度将开始下降。所以初期应采取小批量生产并加以广告宣传；从有20%用户到有80%用户这段时期，应该大批量生产；后期则应适时转产，这样做可以取得较高的经济效果。 第133页，共173页，2022年，5月20日，0点33分，星期四 观众厅地面设计1 问题的提出在影视厅或报告厅，经常会为前边观众遮挡住自己的视线而苦恼。显然，场内的观众都在朝台上看，如果场内地面不做成前低后高的坡度模式，那么前边观众必然会遮挡后面观众的视线。试建立

41、数学模型设计良好的报告厅地面坡度曲线。第134页，共173页，2022年，5月20日，0点33分，星期四建立坐标系oo处在台上的设计视点bb第一排观众的眼睛到x轴的垂 直距离xyadda第一排观众与设计视点的水平距离d相邻两排的排距视线升高标准x表示任一排与设计视点的水平距离求任一排x与设计视点o的竖直距离函数使此曲线满足视线的无遮挡要求。问题第135页，共173页，2022年，5月20日，0点33分，星期四2 问题的假设观众厅地面的纵剖面图一致，只需求中轴线上地面的起伏曲线即可。同一排的座位在同一等高线上。每个坐在座位上的观众的眼睛与地面的距离相等。每个坐在座位上的观众的头与地面的距离也相等

42、。所求曲线只要使观众的视线从紧邻的前一个座位的人的头顶擦过即可。第136页，共173页，2022年，5月20日，0点33分，星期四3 建模设眼睛升起曲线应满足微分方程初始条件obxyadd1）从第一排起，观众眼睛与o点的连线的斜率随排数的增加而增加，而眼睛升起曲线显然与这些直线皆相交，故此升起曲线是凹的。第137页，共173页，2022年，5月20日，0点33分，星期四2）选择某排和相邻排oyx-dC(x,0)C2(x+d,0)MM2M1xN1ABN相似于D第138页，共173页，2022年，5月20日，0点33分，星期四再计算相似于第139页，共173页，2022年，5月20日，0点33分，

43、星期四4 模型求解 微分不等式（比较定理）设函数定义在某个区域上，且满足1）在D上满足存在唯一性定理的条件；2）在D上有不等式则初值问题与的解在它们共同存在区间上满足第140页，共173页，2022年，5月20日，0点33分，星期四第141页，共173页，2022年，5月20日，0点33分，星期四所求曲线的近似曲线方程（折衷法）折衷法第142页，共173页，2022年，5月20日，0点33分，星期四5 总结与讨论有时只需求近似解。方法利用微分不等式建模；模型讨论obxyadd1）视点移动时升起曲线如何求得？2）怎样减少地面的坡度？调整参数、相邻排错位。3）衡量经济的指标？座位尽量多、升起曲线占

44、据的空间尽量少等。第143页，共173页，2022年，5月20日，0点33分，星期四 碳定年代法考古、地质学等方面的专家常用14C测定法（通常称碳定年代法）来估计文物或化石的年代。 第144页，共173页，2022年，5月20日，0点33分，星期四 14C的蜕变规律14C是一种由宇宙射线不断轰击大气层，使大气层产生中子，中子与氮气作用生成的具有放射性的物质。这种放射性碳可氧化成二氧化碳，二氧化碳被植物所吸收，而植物又作为动物的食物，于是放射性碳被带到各种动植物体内。14C是放射性的，无论在空气中还是在生物体内他都在不断蜕变，这种蜕变规律我们可以求出来。通常假定其蜕变速度与该时刻的存余量成正比。

45、第145页，共173页，2022年，5月20日，0点33分，星期四设在时刻t（年），生物体中14C的存量为x(t),生物体的死亡时间记为t0=0,此时14C含量为x0,由假设，初值问题 （1.1）的解为 （1.2）其中，为常数，k前面的符号表示14C的存量是递减的。（1.2）式表明14C是按指数递减的，而常数k可由半衰期确定， 第146页，共173页，2022年，5月20日，0点33分，星期四若14C的半衰期为T,则有 (1.3)将（1.3）代入（1.2）得 即有 （1.4）第147页，共173页，2022年，5月20日，0点33分，星期四碳定年代法的根据 活着的生物通过新陈代谢不断摄取14C

46、,因而他们体内的14C与空气中的14C含量相同，而生物死亡之后，停止摄取14C，因而尸体内的14C由于不断蜕变而不断减少。碳定年代法就是根据生物体死亡之后体内14C蜕变减少量的变化情况来判断生物的死亡时间的。第148页，共173页，2022年，5月20日，0点33分，星期四碳定年代法的计算由（1.4）解得 (1.5)由于x(0),x(t)不便于测量，我们可把(1.5)作如下修改.对(1.2)式两边求导数，得 (1.6)而 (1.7)第149页，共173页，2022年，5月20日，0点33分，星期四(1.6)和(1.7)两式相除，得 将上式代入(1.5)，得 (1.8)这样由(1.8)可知，只要

47、知道生物体在死亡时体内14C的蜕变速度 和现在时刻t的蜕变速度 ,就可以求得生物体的死亡时间了，在实际计算上，都假定现代生物体中14C的蜕变速度与生物体死亡时代生物体中14C的蜕变速度相同。第150页，共173页，2022年，5月20日，0点33分，星期四马王堆一号墓年代的确定马王堆一号墓于1972年8月出土，其时测得出土的木炭标本的C14平均原子蜕变数为29.78/s,而新砍伐木头烧成的木炭中C14 平均原子蜕变数为38.37/s,又知C14的半衰期为5568年，这样，我们可以把 , ， T=5568 年代入(1.8)，得 这样就估算出马王堆一号墓大约是在2000多年前。第151页，共173

48、页，2022年，5月20日，0点33分，星期四两个注记（1）马王堆中的古代科技之谜素纱蝉衣：两件轻薄的衣服，丝绸，极轻且两千年不腐，南京云锦研究所接受国家科技攻关，用了二十年时间，于1990年成功研制出类似素纱蝉衣的复制品，但该复制品比汉代的还重50克，已不可能再轻了。女尸千年不腐：病理知识：女尸解剖显示患有非常严重的冠心病；肺部有肺结核的钙化，肺部钙化是肺结核痊愈后的表现。2000年后的今天，要想控制肺结核，除自身的第152页，共173页，2022年，5月20日，0点33分，星期四抵抗力要强外，还要有好的营养，要想痊愈是很困难的。两处胆结石，其一在胆总管，有蚕豆大，胆道被堵得水泄不通。三种寄

49、生虫，其中竟有血吸虫，其症状应为腹胀如鼓，骨瘦如柴，但该女子皮下脂肪异常丰满，显然血吸虫被有效的控制住了。该西汉贵妇生前病魔缠身，但从其遗体上未发现长期卧床养病的迹象。一个同时患有这么多疾病的人，能够长期稳定控制病情，在今天也是一个奇迹，说明汉代医术已达到了相当高的水平。第153页，共173页，2022年，5月20日，0点33分，星期四（2）碳定年代法的不足 现在，14C年代测定法已受到怀疑，在2500-10000年前这段时间中与其他断代法的结果有差异。1966年，耶鲁实验室的Minze Stuiver 和加利福尼亚大学圣地亚哥分校的Hans E.Suess在一份报告中指出了这一时期使14C年

50、代测定产生误差的根本原因。在那个年代，宇宙射线的放射强度减弱了，偏差的峰值发生在大约6000年以前。第154页，共173页，2022年，5月20日，0点33分，星期四这两位研究人员的结论出自对Brist/econe松树所作的14C年代测定的结果，因为这种松树同时还提供了精确的年轮断代。他们提出了一个很成功的误差公式，用来校正根据14C断代定出的2300-6000年前这期间的年代：真正的年代=14C年1.4900。第155页，共173页，2022年，5月20日，0点33分，星期四 在研究人口或种群数量的实际增长情况时，有时采用离散化的时间变量更为方便。例如，有些种群具有相对较为固定的繁殖期，按时

51、段统计种群数量更接近种群的实际增长方式。人口增长虽无这种特征，但人口普查不可能连续统计，任何方式的普查都只能得到一些离散时刻的人口总量（指较大范围的普查）。这样，如何建立人口问题的离散模型的问题十分自然地提了出来。第156页，共173页，2022年，5月20日，0点33分，星期四对于k阶差分方程F( n; xn, xn+1, , xn+k ) = 0 (3-6)若有xn = x (n), 满足F(n; x(n), x(n + 1) , , x(n + k ) = 0,则称xn = x (n)是差分方程(3-6)的解, 包含个任意常数的解称为(3-6)的通解, x0, x1, , xk-1为已知时称为(3-6)的初始条件,通解中的任意常数都由初始条件确定后的解称为(3-6)的特解. 若x0, x1, , xk-1已知, 则形如xn+k = g(n; xn, xn+1, , xn+k-1 )的差分方程的解可以在计算机上实现.第157页，共173页，2022年，5月20日，0点33分，星期四 若有常数a是差分方程(3-6)的解, 即F (n; a, a, , a ) = 0,则称 a是差分方程(3-6)的平衡点. 又对差分方程(3-6)的任意由初始条件确定的解 xn= x(n)都有xna (n), 则称这个平衡点a是稳定的. 一阶常系数线性差