数学建模ch离散模型.pptx





《数学建模ch离散模型.pptx》由会员分享,可在线阅读,更多相关《数学建模ch离散模型.pptx(203页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、一、过河问题 问题 有三名商人各带一名随从要乘一条小船过河,这条船每次最多只能容纳两个人,并且由于某种原因,商人们总是提防着随从们,预感到一旦在任何地方只要随从人数多于商人数,就会对商人构成危害。但是由于商人们控制着如何乘船的指挥权,所以商人们就可以制定一个过河方案,以确保商人们的安全。试求出这个方案。第1页/共203页 建模 设在渡河过程中,此岸的商人个数为 随从个数为 以 表示此岸的状态向量,即在 中有一部分对商人是安全的,称为容许状态集合,记为 即有第2页/共203页 在上图中,实点即表示为容许状态的集合.乘船的方案称为决策,仍然用向量 来表示,即 名商人和 名随从同坐一条船.在这些决策
2、中,有第3页/共203页是符合条件的,称为容许决策。容许决策的全体组成集合构成容许决策的集合,记为 在这个问题中,容许决策的集合为 小船从此岸到彼岸的一次航行,会使两岸的状态发生一次变化,此称为状态的转移。用第4页/共203页表示状态的转移。其中 用表示在状态 下的决策。当 为奇数时,表示从此岸到彼岸,当 为偶数时,表示从彼岸到此岸。所以公式称为状态转移公式。所以,该问题转变成寻找一系列的决策 使状态 按由初始状态经过有限次的转移达到第5页/共203页 建立坐标系统,并在坐标平面上建立的刻度单位。做网格线,网格线上的每一个交点代表一个状态(用实点表示)。黄色曲线弧表示向彼岸渡人,绿色曲线弧表示
3、从彼岸返回。容许决策 表现为从一个实点向另一个实点的转移。当 为奇数时,容许决策表现的是向下及向左的移动,当 为偶数时容许决策表现的是向上及向右的移 解模第6页/共203页动。整个状态的转移用下面的表格来表示。第7页/共203页序号序号状态状态决策决策序号序号状态状态决策决策172839410511612第8页/共203页 分析 从上表中可以看到,该方案是可行的。第9页/共203页二、差分方程 本节介绍离散模型中的一种重要类型差分方程及相应的解法.第10页/共203页 1.差分方程 离散模型的基本形式是:下一变量为由当前值、先前值及 构成的函数,即 方程称为差分方程.例1 某产品当年的产量与前
4、一年的产量关系为:当年产量在去年产量增加 则相应的差分方程为第11页/共203页 一类比较简单的差分方程为具体形式为:第12页/共203页 在差分方程中,若仅有 的值,即方程具有形式则方程称为一阶差分方程.若只有 的值,则方程称为二阶差分方程.一般差分方程的阶定义为出现在方程中最高阶与最低阶的差值.例如一个二阶线性方程为所谓同类线性差分方程指的是方程具有形式第13页/共203页而方程则不是同类方程.第14页/共203页例2 设在一场战斗中,交战双方为 军和 军,军的一个单位一次可摧毁 军的 个单位,军的一个单位可摧毁 的 个单位.设经过 次战斗后,所剩下的人数,则有此是由两个变量关连的差分方程
5、,经过转化,方程可变为第15页/共203页由此得到的是一个二阶的差分方程.对与给定的初始值 经过若干次的迭代,最终能求出 应用 拥有10000人的 军和拥有5000人的 展开一场战斗,军的杀伤力为 军的杀伤力为 用上式预测战斗结果.分析 为使问题简化,以1000人为一个单位,即第16页/共203页代入上式,有结果-0.160.611.397.557.657.8686542.213.08458.198.659.25103210第17页/共203页当战斗进行到第6次时,已经没有意义了,此时 军还有7550人.问题:如果有俘虏产生,该问题又该如何?第18页/共203页 2.矩阵表示 由若干个变量构成
6、的线性差分方程,在很多情况下可通过矩阵西形式表现出来.例如在上面的两军交战的模型中,相应的矩阵形式为而相应的进程关系可表达为第19页/共203页其中,矩阵 又称为状态转移矩阵.第20页/共203页应用 设有某种生物,一年后成熟并有繁殖能力.此种生物在状态 时,有状态 年龄在1岁之前的幼虫总数;年龄在12岁之间的壮年虫总数;年龄在2岁以上的老年虫总数;不同年龄段的此种生物的繁殖和死亡率有下表所示:第21页/共203页组别组别繁殖率繁殖率死亡率死亡率00.10.30.20.10.3由此得到状态转移矩阵及关系表达式第22页/共203页是否进入稳定状态,看矩阵的最大特征值 若 则总数将无限制增长,总数
7、将趋于稳定,总数将逐渐减少.在Matlab下,求出特征值,值为第23页/共203页此说明基本进入稳定状态.造成这一现象的主要原因是繁殖率偏低.事实上,若将繁殖率改为 则特征值为第24页/共203页 3.模型应用问题1 定期存款的计算 某人在银行存款,设期初为 之后在每年的年底再存入定额为 的钱.年息为 若 是存款 年后的总额,则有记 则有第25页/共203页此方程为一阶线性差分方程,其数值解为 在Matlab下,编写如下程序,运行后得:15年后存款总额为 元.例如:某人期初存入2万元,计划每年年底再存入5000,年利息为 求15年后在银行的存钱总额.如果要让存款达到20万元,求出年数.第26页
8、/共203页存款曲线图示第27页/共203页 问题2 贷款还款计划 问题的提出 当今人们消费,经常会遇到分期付款消费等现象.对实际贷款额和银行利率,该如何制定相应的还款计划.模型分析 设 为 年后所欠的钱数,为每月偿还的钱数,为还清所需的年数,是与欠款有关的年利息率.则 表示欠款总额,由条件,知下一年的欠款钱数为第28页/共203页其中:模型建立 由上式,得当 时则有,即由此得第29页/共203页 应用:设某人贷款20万元,年利率为 计划15年还清,求每月应还的钱额.解 将数据代入上式,即得 元.第30页/共203页还贷曲线图第31页/共203页 讨论:当利率改变时,还贷计划应该做怎么样的调整
9、.例如当利率下降 还款是否也下降 当利率上升 而你又不想改变每月支付的金额,则还款年数将 增加到多少?解 在求解公式中,将年利率改为 得 元,降幅为 若利率上升 则每月的还款额为 涨幅为 若不想增加还款额的话,则还款期限为18年.第32页/共203页 进一步的讨论:是否存在一辈子都还不清的可能?在上面的情况下,假设贷款者仍然每月还款 元,则当年利率上升到多少,还款者要永远还下去?在关系式取 则有第33页/共203页 问题3 生物种群的生长问题 某种生物群种的增长情况遵从著名的Logistic方程:其中 是该种第 代的个体总量与该群体所能达到的最大个体容量之比.为比例系数.注意到该方程是一个一阶
10、的非线性差分方程.第34页/共203页在Matlab下,对初始值 的不同取值,得到如下的结果:第35页/共203页第36页/共203页第37页/共203页第38页/共203页 问题4 减肥与饮食问题 一个女子每天摄入2500卡的食物,其中1200卡用于基本的新陈代谢,并用每公斤体重16卡为每日锻炼的消耗,其它剩余的转换成脂肪.设10000卡等效为1公斤脂肪.在星期六上午,她的准确体重为57.1526公斤.在周三她饱餐了一顿,摄入了3500卡的食物.试建立一个数学模型求第 天的体重 并用它来作以下预测:到星期六时她的体重第39页/共203页为保持体重不变,每日应摄入的热量;周后她可减轻到的最小体
11、重;若她想在8周后减轻到50.802344公斤,每天摄入量应如何限制.第40页/共203页 模型分析 我们不考虑由于吸入氧气而产生的热量.由于新陈代谢及锻炼所需的热量为 故每天能转化成脂肪的热量为 因而相应的脂肪增加值为从而第 的体重为 第41页/共203页 设周六的体重为 则在下周二的体重为 但当天多摄入了1000卡的热量,故周三的实际体重为以它为初始值,可计算出周六的体重为设每天的摄入量为常数 则模型为,欲使体重不变,即有第42页/共203页若要维持体重不变,则 卡.若完全拒绝饮食,但基本的新陈代谢需要维持及锻炼照常,因而有从而有关系代入 周所需要的天数,则可得到那时的体重.第43页/共2
12、03页 使用中的模型,或改写成另一个形式欲使体重8周后成为 代入上式,得第44页/共203页代入 即可得到每天应摄入的热量数.下表给出了当 时体重在第 周时的情况.第三列表示希望在第 周时体重为某值的热量摄入量变化值表.第45页/共203页 周周 (卡)(卡)549.952501043.1411562030.6116093019.411759409.41833500.451877第46页/共203页 问题5 数列与黄金分割 问题的提出 提出了这样一个问题:一对小兔子二个月后可以生兔子,而成熟兔子每月可生一对小兔子.假如去年12月底养一对小兔子,问到今年年底共有多少对兔子.13个月份中的兔子的对
13、数如下表所示:第47页/共203页12123456112358137891011122133558914423313个月份中的兔子的对数第48页/共203页 以 表示各月份中兔子的对数,则有关系 定义 称数列为 为 数列,若数列有关系:法国数学家奇拉特在 死后400年时证明了第49页/共203页 从而证明了第50页/共203页 黄金分割与斐氏数列 黄金分割据说是由达芬奇首先提出的.所谓黄金分割指的是按中外比割的.即在下图中满足关系记则有第51页/共203页即有 解此方程得第52页/共203页第53页/共203页第54页/共203页三、马氏链及其应用 1.一个简单的例子 我们知道,人寿保险公司最
14、为关心的是投保人的健康与疾病以及相应的风险。通过下面的例子我们来看保险公司是如何处理这类问题的。第55页/共203页 问题的提出 设 表示年龄的时段,假定在一年中,今年健康而明年患病的概率是 而今年患病明年转为健康的概率为 假设一个人在投保时处于健康状态,我们来研究若干年之后他分别处于这两种状态的概率。第56页/共203页 建模 用随机变量 表示第 年的状态,表示健康,表示疾病。以 表示第 年状态为 的概率。即以 表示今年状态处于 明年状态处于 的概率,即第57页/共203页由全概率公式得到:即由假设,第58页/共203页再由于投保人处于健康状态,即由此得到若投保人在开始时处于疾病状态,即则有
15、第59页/共203页 从两张表中可以看到,无论投保人在初始时处于什么状态,当时间趋于无穷大时,该时刻的状态趋于稳定,且与初始值无关。即第60页/共203页两种状态的转移概率 意义 若将众多投保人处于两种状态的比例,视为投保人处于两种状态的概率,例如健康人占3/4,病人占1/4,即 则同样可计算出第61页/共203页 由上面的分析可以看出,对于给定的状态转移概率,时的状态概率,趋向于稳定值,该值与初始值无关,这是马氏链的重要性质。第62页/共203页 把人的死亡看作第三种状态,用 来表示,相应的转移概率如下图表示。三种状态的转移概率 仍以 表示状态为 时的概率,表示状态转移概率,即有第63页/共
16、203页平行于式,有 设投保人在期初处于健康状态,则由可计算出若干年后他处于各个状态的概率。第64页/共203页 表中最后一列数据是通过预测得到的。从表中的数据又可以看到,无论投保人在期初处于什么状态,当 时,总有第65页/共203页 2.马尔可夫链 假设 1.系统是随时间的发展而离散为 2.在任何时刻,系统的状态为有限多个。在时间 时,系统的状态的 的取值为 3.在时刻 时系统处于各状态的概率只与时刻 时系统所处的概率与转移概率有关。满足以上三个假设的系统的随机发展过程称为马尔可夫过程或马氏链。第66页/共203页 设在时刻 时系统处于状态 的概率为 行向量称为状态概率向量,由概率的意义,向
17、量应该满足及第67页/共203页 设在时刻 处于状态 的系统转移到 时刻处于 的概率为 它应该满足 1.2.引如概率转移矩阵第68页/共203页由假设3,再由全概率公式得第69页/共203页 用矩阵的方法来表示的话,可以写成 简单地可以写成由此可得系统在时刻 时的状态向量为其中 为时刻 时系统的状态概率向量,又称为状态初始向量。第70页/共203页例 在前两例中,初始向量与概率转移矩阵分别为我们通过下面的例子具体说明:第71页/共203页上式表明在时刻 时投保人处于患病状态的概率为:从上面的例子中可以看出,对于马氏链模型,最重要的是构造状态 及概率转移矩阵 由此对于给定的初始状态 由可计算出任
18、意时刻 的状态第72页/共203页 正则链定义 一个有 个状态的马氏链,如果存在正整数 使从任意状态 经 次的转移,能以大于零的概率到达状态 则称这样的链为正则链.定理1 设马氏链的转移矩阵为 则该链为正则链的充分必要条件是存在 使得 第73页/共203页定理2 正则链存在唯一的极限状态概率满足 与初始状态概率 无关,且 及第74页/共203页例1 设 则由此确定的马氏链为正则链。令 满足式,即有由此得到方程组第75页/共203页联系则得到故方程组的解为第76页/共203页这和前面的结果是相吻合的。第77页/共203页例2 设因第78页/共203页故由此确定的马氏链是正则链。令由方程,确定方程
19、组第79页/共203页从方程中解出 即第80页/共203页 吸收链 定义 如果存在某个状态转移概率 则称状态 是吸收的.如果马氏链中含有吸收状态,并且从每一个非吸收状态出发都可以达到某个吸收状态,则称这个马氏链为吸收链。例如在前面三个状态的转移概率中,转移概率矩阵为第81页/共203页并且从每个状态最终都转移到第三种状态,因而这样的链是吸收链。注 吸收链的特征是:任一状态一旦进入该状态就将停留在该状态。第82页/共203页 含有 个吸收状态和 非吸收状态的吸收链的状态转移概率矩阵的标准形式是其中 是单位矩阵。第83页/共203页定理3 对于具有标准形式的状态转移概率矩阵,有如下的性质:矩阵 具
20、有零极限,即矩阵 可逆且记 则矩阵的第 行元素之和值是从非吸收状态出发被某个吸收状态吸收之前的平均转移次数。第84页/共203页记 则矩阵 的元素 是从非吸收状态 出发而被状态 吸收的概率。第85页/共203页 在前面的例2中,将 改写成则第86页/共203页则第87页/共203页 应用 基因遗传问题 生物的外部特征是由生物体内的基因决定的。基因分优势与劣势基因两种。分别表示为 对于生物的某个外部特征,体内有两个基因与之对应。由于体内的每个基因都可以是两种基因之一,因此体内的基因对类型可能有三种:分别被称为优种、混种和劣种。按基因理论:含优种和混种的基因个体类型,其外部特征呈优势;而含劣势基因
21、类型的个体,其外部特征呈劣势。第88页/共203页 生物在繁殖时,后代随机地继承父亲和母亲的两个基因中的各一个而形成自己的基因对。因此后代成为优种、劣种、混种基因类型的概率是不同的。下面讨论两种基因繁殖后代的情况第89页/共203页 一、永远与混种繁殖后代的情况 假设一个个体是优种,而另一个个体是混种,则它们的直接后代成为优种 、混种 、劣种 的概率分别为 假设一个个体是混种,而另一个个体是混种,则它们的直接后代成为优种 、混种 、劣种 的概率分别为第90页/共203页 假设一个个体是劣种,而另一个个体是混种,则它们的直接后代成为优种 、混种 、劣种 的概率分别为由此得到概率转移矩阵第91页/
22、共203页由前面的例2知该链为正则链,极限状态概率向量为 上式表明,经过长时间的繁殖过程,后代的外部特征呈优势的概率是优种和混种概率的和,这个量与初始的个体所含基因的种类无关。第92页/共203页 2.近亲繁殖的结果 假设最初的父母可以是优种、混种或劣种,它们有大量的后代,这些后代又随机地雌雄交配后代,今来分析它们后代的演变情况。由于每次繁殖都是随机地配对父亲和母亲,而父亲和母亲可以是 中的一种,组合后就有 六种状态,分别记为当父母都是优种 时,后代必然是优种 因此有第93页/共203页 同理,当父母都是劣种时,后代只能是劣种,由此得 当父母一方为 而另一方为 时,当前状态可能是 因而再次配对
23、产生的可能结果有因此,有第94页/共203页 当父母方为 对时,其后代只可能是 因而再次配对之后之可能产生 所以 当父母方为 对时,其后代可能是 甲 乙 因而相应的概率为第95页/共203页所以概率转移矩阵为第96页/共203页从上面中可以看到状态1和状态2是吸收状态。所以该链为吸收链。由前面的计算公式得到的行和第97页/共203页根据矩阵 和 的性质,上式表明从状态3出发经过第98页/共203页代后它们的后代都会变成优种或劣种,从状态3出发其后代全变为优种的概率为 上面表明:近亲繁殖的后代变成劣种的可能性很大。第99页/共203页三、钢琴销售的存储策略 问题的提出 一家商店根据以往经验,平均
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学 建模 ch 离散 模型

限制150内