数学地质第七章 有序地质量最优分割法.ppt

第三节 多元有序数据的最优分割 得正规化数据矩阵 第三节 多元有序数据的最优分割 根据正规化数据，将样品段{ i,…，j}的段直径定义为 （7-9） 式中 （1≤ i≤j≤n， β=1，2，…，p） （7-10） 第三节 多元有序数据的最优分割 若n个有序样品分为K段，每段内有nk个样品，则多元有序数据最优分割的原理与单元有序数据最优分割一样，使组内离差平方和 应当指出，样品的段直径除了用式（7-9）定义外，还可用其它方法定义。 第四节 最优分割法的计算步骤 1、数据正规化 设原始数据矩阵为 第四节 最优分割法的计算步骤 将X中的元素xil变换为 （i=1,2,…，n；l=1，2，…，p） 得正规化数据矩阵 第四节 最优分割法的计算步骤 2、计算段直径矩阵D 其中 因为 第四节 最优分割法的计算步骤 故只须计算 个 ，得 第四节 最优分割法的计算步骤 3、计算全部分割的组内离差平方和（或段直径和）及各种分段的最优分割 1）最优二段分割 由D矩阵对每一个m=n，n-1，…，2计算相应的组内离差平方和，为 找出最小值，确定响应的最优二段分割点，即 第四节 最优分割法的计算步骤 2）最优三段分割 根据D矩阵及最优二段分割结果，对每一个m=n，n-1，…，3计算相应的三段分割的组内离差平方和，为 然后求出最小值，并确定相应的最优三段分割点，为 第四节 最优分割法的计算步骤 3）最优K段分割 根据矩阵D及最优K-1段分割计算结果，对于每一个m=n，n-1，…，K分别计算相应的K段分割的组内离差平方和，为 找出最小值，并确定相应的最优K段分割点，即 第四节 最优分割法的计算步骤 4 绘制W-K曲线 在W-K曲线上，选择曲线拐点对应的K值（取整）作为最终分段数。 第四节 最优分割法的计算步骤 例7-1 某煤矿所采煤层的煤质牌号为主焦煤。在煤巷中见一火成岩墙侵入煤层，致使煤质发生变化，为弄清楚煤质变化情况，从火成岩附近每隔0.5m依次取一煤样，获得6个有序煤样的镜煤最大反射率数据为 试进行最优分割。 第四节 最优分割法的计算步骤 此样本的可能分割法共有25-1=31种，今要在这31种分割中找出一种最优的分割（类内差别小，类间差别大）。其作法如下： （1）对原始数据进行正规化变换后得正规化数据，为 （2）计算直径矩阵D，即 第四节 最优分割法的计算步骤 第四节 最优分割法的计算步骤 （3）最优二段分割。由D对于m=6，5，4，3，2时，计算 第四节 最优分割法的计算步骤 第四节 最优分割法的计算步骤 第四节 最优分割法的计算步骤 第四节 最优分割法的计算步骤 从而得到6个样品的最优二段分割为： 其中， 为分割点。 （4）最优三段分割。即对于m=6，5，4，3时，计算 当m=6时，则 分割点为： 第四节 最优分割法的计算步骤 第四节 最优分割法的计算步骤 第四节 最优分割法的计算步骤 第四节 最优分割法的计算步骤 （5）最优四段分割。即对于m=6，5，4时，计算 当m=6时，则 第四节 最优分割法的计算步骤 第四节 最优分割法的计算步骤 第四节 最优分割法的计算步骤 （6）最优五段分割。即对于m=6，5时，计算 第四节 最优分割法的计算步骤 当m=5时，则 其中 则其分割点为： 第四节 最优分割法的计算步骤 由组内离差平方和（W）与分段段数（K）的关系图（图7-2）中看出，三段最优分割是最终分段结果。 图7-2 W-K曲线图 第四节 最优分割法的计算步骤 最优分割法可以用来对地质量进行分段（或分类）。尤其在地层划分、水系中化学成分与水系沉积物化学分带、变质岩相带、蚀变带、矿化带等划分中用得最多。有时对一组没有明显顺序的变量也可给予一定的排列顺序进行最优分割。例如把各种煤的镜煤最大平均反射率排序进行最优分割，研究烟煤变质阶段的合理划分等。 本章内容结束 谢谢！ * * 第七章 有序地质量最优分割法 杨永国中国矿业大学资源与地球科学学院 内容提要: 第一节 概述 第二节 单元有序数据的最优分割 第三节 多元有序数据的最优分割 第四节 最优分割法的计算步骤 第一节 概述 地质数据中有相当多是有序的。这些按一定顺序排列的地质变量，叫做有序地质量。例如，沿地层露头剖面采集的岩石标本；钻孔取出的岩芯样品； 这类数据的特点是样品的前后次序不能变更。所以，一些不