多功能聚类分析程序.PDFVIP

维普资讯 · 81 · l3暑 l期 l昀-年2H 枷探化探计算枝术 Vo!．13No．1Feb．19gt COMPUTING TECHNIQUES FOR GEOPHYStCAL AND GEOCHEMtCAL EXPLORATION 多 功 能 聚 类 分 析 程 序 赵 玉 珲 (安徽省地矿局332地质 队) 聚类分析是根据标本多种变量测定数据，应用数学方法依其相似性逐步归组成群的分类 手 段。根据处理对象的不同，可分为研究标本之间关系的Q型和研究变量之间关系的R型两 种 。通常为使各变量对分类贡献处于同一量度的相对水平，又需将原始数据进行标准化或规 格化等变换，在成图方面联接二维谱系图 (枝状 图)时，R型与Q型 的原则又完全不 同。由 于这些原因，以往发表的有关 PC一 1500程序都将 两者分开编制 】，以致使用不便。笔者 者开发的多功能聚类分析程序是将R型与Q型分析通过特殊的编程技巧压缩成一个统一的 程 序，井用扩充的新联接表，按一次计算成群法紧邻联结成枝状 图。程序容量为2509个字节， 约 占前述同类单个 (R或Q型)程序的80％，总合程序的45％左右。经测定其计算绘图速度 又较后者提高近一倍 。 一 、 鳙 疆要 点 1． 无论选用R型或Q型分析均用同一数据表，即标本按行、变量按列排歹I】。数据 读 人 后 由程序指令 自动相应转换 (源程序7ooI~)。 2． 对原始数据的各种变换可 自由选择，由专门的子程序完成。如标准化 (400行)和规 格化 (450行)等。 3． 分类统计量亦可 自由选择。如采用距离系数 (500行) 、相似系数 (600行)和相关 系数 (7oo~3)等。 4． 联接水平的隐性转换是本程序特殊设计技巧的核心。众所周知，距离系数是按从小 到大的数值水平联接的，而相似系数和相关系数则恰好相反，且其值按大小波动于1到一12_ 间，但若将各计算值从 1中减去，则可构成。到2_2间的倒序排列，这样就可使用距离系数联 接方式的遁道而不必另行设计 。这种 按倒序联接、按正序打印数值的方式就是隐性转换 (源 程序142行 ，156行等)。它能将R型与Q型变成同一成 图方法 。 5． 新联接表的构成：以往程序多采用逐步形成往 (聚合法)，其框 图逻辑复杂，编程 量大，运算成图较慢。本程序采用 《数学地质引论》 (FA下称 引论)83—85页一次计 算法的基本方法，形成新联接表 (表 1)。解释如下t ①、制图参数中除D列外是顺时针连接各群框形线的距离码 ，其结点坐标依次为 (A，U) ÷ (A(x)，U)一 (Aix)，L)÷(E，L)(源程序270~)。辅助参数主要功能是逻辑判断和 计算距离码 。Aix)，B()，Cix)和S()是源程序 中的一维数组，其意义和计算方法分 别 相 当 《引论联接表的A、 、C、 。距离码乘 以X或y轴方 向上的图形系数再加 修正量 得到其代表的真实绘图距离。 肚墒日期 1990-~8凡 地址：安擞省马鞍 山市向山 (243034) 维普资讯 ·82· 袖探化探计算撞水 l3善 衰1 枝；圈联接表 (以下文实例一计算为侧) 制 凰 参 数 I 辅 助 事 薮 x(Y) D D(x) i u L lA(x) A E Iu， L B(x)c(z) S(x) W Q l 12．5 12 13 11．278 O