数学地质实验指导书（教材）(11)

2） 找出最优二段分割点

Wm(2;?1(m))?min{Wm(2;j)} （1?j?m?1） （8－6）

3） N个样品的最优二段分割为：

{x1,x2,…..xa1(N)} { xa1(N)+1…..xN} （8－7）

最优二段分割点为a1(N)。

（2）最优三段分割

1） 由D和最优二段分割结果，对每个m = N,N-1,…..3 计算相应的组内离差平方和。

Wm(3;?1(j),j)?Wj(2;?1(j))?d(j?1,m) （j=2,3,…,m-1） （8－8）

2） 确定最优三段分割及分割点

Wm(3;?1(m),?2(m))?min{Wm(3;?1(j),j)} （2?j?m?1） （8－9）

N个点的最优三段分割为：

{x1,x2,x?1(N)?1}{x?1(N)?1...x?2(N)}{x?(N)?1...xN} （8－10）

最优三段分割点为：?1(N),?2(N)

（3）最优K段分割

1） 由D和K－1段分割计算的结果，对每个m = N,N-1,….k分别计算相应的组内离差平方和。

Wm(k;?1(j),?2(j),...,?k?2(j),j)?Wj(k?1;?1(j),?2(j),...,?k?2(j))?d(j?1,N)

(j=k-1,k,…,m-1） （8－11）

2） 确定最优K段分割及分割点

Wm(k;?1(m),?2(m),...,?k?1(m))?min{Wm(k;?1(j),?2(j),...,?k?1(j),j)}

（k?1?j?m?1） （8－12）

N个样品的最优K段分割为：

{x1...x?1(N)}{x?1(N)?1...x?2(N)}...{x?(K?2)(N)...x?(K?1)(N)}{x?(K?1)...xN} （8－13）

最优K段分割点为：?1(N),?2(N),...?(k?1)(N)

5、给定K值或绘制W-K曲线，确定最终分段段数。

33

三、程序设计框图

开始挑选最小变差（8－5式）打开标准格式的原始数据文件（包括N和最高段数K）确定最优L段分割（8－6式）数据正规化（实验一）L＝L+1计算L段分割的组内离差平方和（8－11式）Y计算段径矩阵D数据正规化（实验一）挑选最小变差（8－12式）计算段径矩阵D（8－4式）确定最优L段分割（8－13式）L＝2计算L段分割的组内离差平方和（8－5式）L≤KN输出2至K段的最优分割点和最优分割结束