第九章马尔可夫概型分析ppt.ppt

* 第九章马尔可夫概型分析 §1 概述 §4遍历定理与极限分布 §5马儿可夫概型检验 §6 应用简例 §2马尔可夫概型 §3马儿可夫链的转移概率 §1 概述 任何一个地质过程都可认为是某些确定性地质因素和随机性地质因素相互作用的产物。也就是说，地质过程都应是确定型和随机型地质过程在时间上和空间上叠加的结果，因此在地质研究工作中，应重视上述两类地质过程的作用。但是，由于产生随机地质过程的机理十分复杂，甚至连随机性地质因素也不容易或不可能观察，这就给用数学校型(一个或几个)来描述随机型地质过程带来了极大的、甚至是在目前条件下不可能克服的困难。这样以来，人们自然就会侧重于确定型地质过程的研究，而使随机型地质过程的研究成为地质工作中的一个薄弱环节。 1949年，维斯捷列马斯在研究复式沉积层形成问题时，首先应用了马尔可夫链。1984年，在莫斯科举行的第27届国际地质会议上，维斯短列乌斯、阿柑特伯格(F．P．Agterbers)等数学地质学家发表了莫斯科宣言，其基本思想是；随机型模型应作为数学地质模型的基础，并且肯定了马尔可夫过程(特别是马尔可夫链)在地质学中应占有特殊的地位。 目前，常用马尔可夫概型分析研究沉积旋回，进行地层对比，查明火山岩系的喷出顺序和侵入杂岩体中各个侵入体形成的先后顺序，划分矿床的成矿期和成矿阶段、揭示各个成矿阶段的空间分布等。 到目前为止，研究随机型地质过程所涉及的数学方法也仅限于马尔可夫概型分析。 随机过程是概率论的基本概念之一，它是依赖于参数t的一族随机变量，记为 这里的参数t一般是时间，T是它的变化范围，随机变量x(ti)也称作随机过程 x(t)在t＝ti ∈T时的状态。 当随机过程在时刻t1所处的状态x(t1)为已知的条件下，若随机过程在时刻t(t＞t1)所处的状态x(t)与随机过程在t1时刻之前发生的状态无关，那么这样的随机过程就称为马尔可夫过程。 用分布函数来描述就是：如果对于参数t的任意n(n＞3)个数值t1＜t2＜…＜tn，在条件x(ti)＝xi，i＝1，2，…，n一1下随机变量x(tn)的分布函数恰好等于在条件x(tn-1)＝xn-1 的分布函数，即 9.2马尔可夫概型 则称x(t)为马尔可夫过程,条件分布函数 是马尔可夫过程的概率模型,上式为马尔可夫过程的转移概率。 马尔可夫过程又称为“无后效随机过程”。所谓无后效性就是在已知随机过程现在状态的情况下，以后它所处的状态与以前它所处的状态无关。可以将其理解为．这个过程的历史对未来的全部影响集中在最后时刻的状态中，即认为过程的任何观测结果只与紧前面的观测结果有关。就是在已经知道过程“现在”的条件下，其“将来”不依赖于“过去”。 状态和参数都是离散的马尔可夫过程称作马尔可夫链，即过程为 马尔可夫链适用于时间离散、状态离散的时间序列。但是，在研究地质过程时，有时能直接确定过程在时间上的先后顺序，有时则只能间接地以空间上的上下、前后、左右关系宋代替；也就是说，在地质过程研究中，有时可以挨到确定的时间序列，有时只能间接地用距离来代替时间参数：但是，只要空间序列有类似于马尔可夫性质的关系存在，则仍然可以应用马尔可夫链对这种序列进行研究。在此，我们将既适用于时间序列又适用于空间序列的马尔可夫概率模型统称为“马尔可夫概型”。 若马尔可夫过程的转移概率随着时间的推移而发生变化，则称其为非齐次或非乎稳马尔可夫过程，转移概率不随时间而变的马尔可夫过程就称为齐次或平稳马尔可夫过程。目前在地质研究中，主要是应用平稳马尔可夫过程。 9.3马儿可夫链的转移概率 设马尔可夫链中可列个发生状态转移的时刻为t1，t2，…tn，…，在已知时刻t＝tn时随机过程xi所处状态为i的条件下，把经过一步转移，即在时刻t＝tn+1(tn+1＞tn)转移到状态j上的概率记为pij， 这个概率称为马尔可夫链的一阶转移概率。 把从状态i到状态j的一阶转移概率pij记为pij(1) pij(k)则是从状态i出发经过k步转移到状态j上的转移概率 对于马尔可夫链来说，转移概率完全括述了它的概率统计特征，因此，如何确定转移概率则成为研究马尔可夫链的一个重要问题。转移概率在理论上是条件概率，而实际应用时则是以转移频率nij/ni．作为条件概率的估计值，即 例：一地层段由五种岩性状态组成，以1表示砾状砂岩和粗砂岩，2表示细砂岩，3表示粉砂岩，4表示粘土岩，5表示灰岩，下面列出了自下而上的岩性状态转移序列，求其转移概率。 1 2 5 1 4 5 4 5 2 4 5 4 3 5 4 1 3 4 3 4 2 3