有序样品的最优分割的算法及其在MATLAB中的实现

《有序样品的最优分割的算法及其在MATLAB中的实现》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、有序样品的最优分割算法及其在Matlab中的实现一、有序样品聚类——最优分割的概念地质数据中，有些样品有一定的排列顺序，如沿地层剖面采集的岩石标本，由钻孔取得的岩芯样品，由测井曲线所得的数据，由岩体中心到围岩的蚀变剖面的样品等，它们是有序地质变量，在对这些有序样品进行分类时，不能打乱样品的前后次序。所以，一些不考虑样品排列顺序的数学处理方法，对此并不适用。有序样品的聚类分析就是对有序样品进行分段的统计方法。对n个有序样品进行n-1分割，就可能有2种划分方法，这每一种分法成为一种分割。在所有的这些分割中，有一种分割使得各段内部之间差异性最小，而短语段之间差

2、异性最大。这种对n个样品分段并使组内离差平方和最小的分割方法，成为最优分割法。这类问题的提法如下：设有一批（N个）按一定顺序排列的样品，每个样品测得p项指标，其原始资料矩阵:x11x12⋯x1Nx21x22⋯x1NX（p×N）=⋮⋮⋱⋮xp1xp2⋯xpN其中元素xij表示第j个样品的第i个指标的观测值。现在要把此N个样品按顺序（不破坏序列的连续性）进行分割（分段或者分类）。其所有可能的分割法共有123N-1N-1CN-1+CN-1+CN-1+…+CN-1=2-1种。现在要求在所有分割中找出一种分割法，这种分割法使得各段内样品之间的差异最小，而各段之间的

3、差异最大。各段内部差异最小，即各段内数值变化最小，段内数值变化可用变差或者极差来表示，比如样品段｛xi、xi+1、xi+2、…、xj｝:变差：j2dij=α=i[xα−x(i,j)]j1xi,j=xαj−i+1α=1dij表示样品段｛xi、xi+1、xi+2、…、xj｝内样品间的差异情况，dij小表示段内各样品之间数值比较接近，反之，dij大表示段内各样品数值之间的差异大。极差：pdij=α=1(maxi≤β≤jxαβ−mini≤β≤jxαβ)对于单指标情况，则1/7dij=（maxi≤β≤jxβ−mini≤β≤jxβ）要各段内部的差异最小，即所分成各段

4、变差的总和（即段内离差平方和，称为总变差）为最小。总变差分解公式：S总=S段间+S段内Smnl(x2总=l=1j=1jl−x)=mnl[(x−x)+(x−x)]2l=1j=1jlll=mnl(x−x)2+mnl(x−x)2+2mnl(x−x)(x−x)l=1j=1jlll=1j=1ll=1j=1ijll=mnl(x−x)2+mnl(x−x)2l=1j=1jlll=1j=1l=S段内+S段间其中：D=2mnl(x−x)(x−x)l=1j=1ijll=2m(x−x)nl(x−x)l=1ijlj=1l=2mx−x(nx−nx)l=1lllll=0Smnl(x2

5、段内=l=1j=1jl−xl)为段内离差平方和Smnl(x2段间=l=1j=1l−x)为段间离差平方和所以S段间=S总-S段内对给定的N个样品，S总是个固定的量。若使段内离差平方和为最小，则段间离差平方和必为最大。所以，使总变差（段内离差平方和）为最小的分割法就是最优的分割法。二、有序样品聚类的最优分割意义最优分割在地质研究中是一个非常有用的手段，只要地质体的某些地球化学特征存在规律性的差异，采用最优分割的数学树立方法，就能按顺序在最理想的地方进行分段。通过对地层中采集的具某些地球化学特征样品的最优分割，能在地层的划分对比中发挥重要的辅助作用；在找矿过程

6、中，该方法更显得天独厚的优势，它能进行蚀变、矿化及矿体的准确分带；对岩浆岩相带划分及演化序列的研究也十分有效。第四纪地层岩相和厚度变化很大，给第四纪地层划分和对比地层带来了困难，但第四系地层是不同阶段和不同的环境条件下形成的，如重矿物、微量元素等地球化学特征的分布规律与一定的沉积阶段和沉积环境相对应，因此，2/7最优分割法的数据处理，可能是第四纪地层划分和对比的有效方法。三、最优分割的计算步骤及其计算公式1.数据正规化设原始资料矩阵为x11x12⋯x1Nx21x22⋯x1NX（p×N）=⋮⋮⋱⋮xp1xp2⋯xpN将矩阵X中的元素xij变换为：xij−m

7、in1≤j≤n{xij}zij=(i=1、2、…、p；j=1、2、…、N)max1≤j≤n{xij}−min1≤j≤n{xij}而得矩阵Z(p×N)=[zij]Matlab代码：function[std]=std1(vector)%对矩阵进行标准化%vector为待分割矩阵max1=max(vector);%对列求最大值min1=min(vector);%对列求最小值[a,b]=size(vector);%矩阵大小,a为行数,b为列数forj=1:bstd(j)=(vector(j)-min1)/(max1-min1);end2.计算极差（或变差）矩阵由

8、上述极差（或变差）计算公式得到矩阵：d12d13⋯d1Nd23…d2ND=⋱⋮d