第四章灰色聚类分析(精)（精编版）

1、第四章 灰色聚类分析在本章中，首先介绍了灰色聚类的概念及其类型。其次对灰色星座聚类、灰色关联聚类、灰色变权聚类和灰色定权聚类的原理和计算方法进行了阐述。最后利用实证分析来分析灰色聚类在渔业科学中的应用。第一节 灰色聚类的概念灰色聚类是根据关联矩阵或灰数的白化权函数将一些观测指标或观测对象聚集成若干个可定义类别的方法。一个聚类可以看作是属于同一类观测对象的集合体。在实际问题中，每个观测对象往往具有许多个特征指标，因而难以进行准确的分类。灰色聚类按聚类方法的不同，可分为灰色星座聚类、灰色关联聚类和灰类白化函数聚类等方法。灰色星座聚类是根据样本自身的属性，利用相似性原理定量地确定样本之间的关系，并按

2、这种关系进行自然聚类。灰色关联聚类主要用于同类因素的归并，以使复杂系统得到简化。通过灰色关联聚类，可以分析出许多因素中是否有若干个因素关系十分密切，以便我们既能够用这些因素的综合平均指标或其中的某一个因素来代表这些因素，同时又使信息不受严重损失，从而使得我们在进行大面积调研之前，通过典型抽样数据的灰色关联聚类，可以减少不必要变量（因素）的收集，以节省成本和经费。灰类白化权函数聚类主要用于检查观测对象是否属于事先设定的不同类别，以便区别对待。从计算工作量来看，灰类白化函数要比灰色关联聚类和星座聚类复杂。第二节 灰色星座聚类一，原理和方法星座聚类在灰色聚类中是一种比较简单易行的聚类方法。其基本原理

3、为：将每个样点按一定的数量关系，点在一个上半圆之中，一个样点用一颗“ 星点” 来表示，同类的样点便组成一个 “ 星座” ，然后勾画出区分不同星座的界线，这样就可以进行分类。实质上，它是将一个样本中的大量信息（或指标值），经过原始数据的变换（极差变换）等手段转化成为无量纲，并成为一个简单的空间坐标比较的问题。一般情况下，星座聚类有如下步骤：（1）对原始指标值进行极差变换，并使变换后的数值均落在0 ，180 的闭区间内。计算公式为：式中 aij变换后的数据，以角度表示；xij 原始数据；xjmax第 j 个变量的最大值；xjmin 第 j 个变量的最小值；（j=1，2，3，p为指标数）（i=1，2

4、，3，n 为样点号）（2）对每个指标，根据其对系统变化的影响程度，分别给一个权数wj，使 0 j 0.5，当r时，则可将 xi 与 xj 视为同类特征。r 可根据实际问题的需要来确定，若r 越接近于 1，则分类越细，每一组中的变量相对地越少；若 r 越小，则分类越粗，这时每一组中的变量相对地越多。二、举例分析我们仍以上一节中的例子进行分析，利用灰色绝对关联度的计算方法（见第三章）进行聚类分析。在本例中，共有7 个样本，每一个样本中有6 个指标，为了节约今后调查和收集资料的成本，我们需要将指标进行归类，以达到精简指标的目的。第一步：进行始点零化利用可求得（见下表）；序号样本 1 样本 2 样本

5、3 样本 4 样本 5 样本 6 样本 7 x00 0.00 0.08 0.09 0.08 0.00 -0.01 -0.10 x01 0.00 -4.07 -16.71 -12.13 -15.01 -24.59 -22.16 x02 0.00 -2.50 13.30 14.90 18.60 20.90 24.10 x03 0.00 -10.34 -18.74 -26.63 -33.10 -40.63 -44.95 x04 0.00 -0.37 -1.06 -1.13 -1.17 -1.31 -1.35 x05 0.00 -0.13 -0.14 -0.22 -0.24 -0.30 -0.29 第

6、二步：求、和=0.19；83.59；77.25；151.92 5.72；1.18 ；第三步：求绝对关联度；则可得绝对关联矩阵a 若我们假定绝对关联度的临界值取0.60，则我们依次可检查出： x5 与 x0 一类，x3 与 x1 一类， x5 与 x4 一类。取标号最小的指标作为各类的代表，可将x5 和x4 并入 x0 中成为一类，这样6个指标的聚类结果为：x5、x4、x0， x3、x1，x2也就是说，我们在以后资料的收集中，只要收集x0、x1、x2 三个指标的数据即可。第四节 灰色变权聚类一，灰色变权聚类的概念及其方法设有 n 个聚类对象， m个聚类指标， s个不同灰类，根据第i（i=1，2，

7、n）对象关于 j（j=1，2，m）指标的样本值 xij （i=1，2，n；j=1，2，m）将第 i 个对象归入第 k（k1，2，s）个灰类之中，称为灰色聚类。现假设 j 指标 k 子类的白化权函数为图 4-2 所示的典型白化权函数，则称，为的转折点。典型白化权函数我们记为：x 0 图 4-2 图 4-3 若白化权函数无第一转折点和第二个转折点，即如图 4-3所示，则称为下限测度白化权函数，记为。若白化权函数的第二转折点和第三个转折点重合，即如图 4-4所示，则称为适中测度白化权函数，记为。若白化权函数无第三转折点和第四个转折点，即如图 4-5所示，则称为上限测度白化权函数，记为。图 4-4 图

8、 4-5 通过上述分析，我们可以得到不同情况下的白化权函数。（1）对于图 4-2 所示的典型白化函数，有（2）对于图 4-3 所示的下限测度白化权函数，有（3）对于图 4-4 所示的适中测度白化权函数，有（4）对于图 4-5 所示的上限测度白化权函数，有对于图 4-2 所示的 j 指标 k 子类白化权函数，令=（+）；对于图 4-3 所示的 j 指标 k 子类白化权函数，令=；对于图 4-4 和图 4-5 所示的 j 指标 k 子类白化权函数，令=；我们称为 j 指标 k子类临界值，称=为 j 指标关于 k 子类的权重。现设 xij 为对象 i 关于指标 j 的标本，为 j 指标 k 子类白化

9、权函数，为 j 指标关于k 子类的权重，则称=为对象 i 属于 k 灰类的灰色变权聚类系数。称 i=（）=为对象 i 的聚类系数向量。现设，则称对象 i 属于灰类 k*。灰色变权聚类适用于指标意义、量纲皆相同的情形。当聚类指标的意义、量纲不同且不同指标的样本值在数量上悬殊较大时，不宜采用灰色变权聚类。二、举例分析设有三个沿海渔区，三个聚类指标分别为海洋捕捞产值、养殖业产值和水产品加工产值，第 i 个渔区第 j 个指标的样本值 xij ( i,j=1,2,3 如矩阵 a 所示：试按高产值、中等产值和低产值作综合聚类。现假定关于指标海洋捕捞产值、养殖业产值和水产品加工产值的白化权函数分别为：，由以

10、上白化权函数分别为：于是，我们就有：，由公式可得：，再由公式=可得，当 i=1 时，则有同理可得：，则由公式可得：由上述表明，第二渔区属于渔业经济中等发达的区域，第一和第三渔区属于高度发达的区域。第五节 灰色定权聚类一，灰色定权聚类的一般方法当聚类指标意义不同、量纲不同，且在数量上悬殊很大时，若我们采用灰色变权聚类可能导致某些指标参与聚类的作用十分微弱。解决这一问题有两种方法：一是采用第二章中有关原始数据处理方法（如初值化或均值化）进行无量纲处理，然后进行聚类。这种方式对所有聚类指标都一视同仁，不能反映不同指标在聚类过程中作用的差异性。另一种方法就是对各聚类指标事先赋权重，赋予权重的方法有很多

11、，一般采用层次分析法。此种方法事先赋予权重的聚类方法，我们称之为灰色定权聚类。灰色定权聚类可按下列步骤进行：第一步：给出 j 指标 k 子类白化权函数（j=1，2，m；k=1，2，s）；第二步：根据定性分析结论确定各指标的聚类权j（j=1，2，m）；第三步：利用第一步和第二步得出的白化权函数（j=1，2，m；k=1，2，s），聚类权 j（j=1，2，m）以及对象 i 关于 j 指标的样本值 xij（i=1，2，n；j=1，2，m）计算出定权聚类系数，i=1，2，s。第四步：若，则断定对象 i 属于灰类 k*。二，举例分析我们同样以上一节中的例子来加以说明。第一步：确定白化权函数与上一节一样。第

12、二步：确定捕捞、养殖和加工各指标的聚类权重。，第三步：由公式进行计算同理可得：，则由公式可得：由上述表明，第二渔区属于渔业经济中等发达的区域，第一属于高度发达的区域，而第三渔区属于渔业不发达区域。第六节 灰色聚类在渔业科学中应用实例分析灰色聚类在渔业科学得到广泛的应用，根据目前所查阅的文献和资料分析，主要在渔业区域经济、鱼类种群划分、鱼类营养价值评价、渔业水域环境评价等方面。现我们结合有关文献的报道，具体来分析灰色聚类在渔业科学中应用的情况。一、灰色聚类在渔业区域经济中的应用在渔业区域经济研究中，依据各区域的自然资源条件和渔业经济的发展水平，来制定今后发展战略是十分重要的。因此，利用灰色聚类方

13、法，对我国沿海地区的渔业经济进行科学分类划区，以便针对性地制定不同类型的发展规划和实行分类指导，对确保我国沿海地区渔业经济的可持续发展具有十分重要的意义。陈新军 2001年在上海水产大学学报上发表了中国沿海省市渔业经济区域类型的初步探讨。文章选取1997年我国渔业发展共11 个省市的统计材料，它们分别是天津、 河北、 辽宁、 上海、 江苏、 浙江、 福建、 山东、 广东、 广西和海南。结合我国水产业的特性和发展状况，建立渔业技术经济综合评价指标，采用了生产总量、生产率等17 个指标，分别是水产品总产量（吨）、海洋捕捞产量（吨）、海洋养殖产量（吨）、淡水捕捞产量（吨）、淡水养殖产量（吨）、远洋渔

14、业产量（吨）、总产值（万元）、水产品加工产量（吨）、人均水产品产量（吨/人）、人均渔业总产值（元/人）、渔民劳均纯收入（元/人）、渔民人均纯收入（元/人）、人均固定资产投入（元/人）、平均海水和淡水养殖单产（千克/公顷）、海水养殖的水面利用率（）、淡水养殖水面的利用率（）（除稻田养殖水面外）。利用灰色星座聚类对中国海省市渔业经济区域类型进行初步划分。1，聚类分析的结果权重取 1/17，经过计算求得各省市（样点）的横坐标（xi 值）和纵坐标（ yi值），如表 4-2。将表 4-2绘制星座聚类图（如图4-6）。表 4-2 聚类分析结果样点天津河北辽宁上海江苏浙江x 值0.52630.58960.2

15、052 0.40640.3354 0.0829y 值0.33960.45680.6318 0.35480.4102 0.6794样点福建山东广东广西海南x 值0.0276-0.39900.0224 0.63670.7256y 值0.57650.58550.5708 0.48120.3004引自陈新军（ 2001）2，计算综合指标值（ z）（见表 43）表 4-3 各省市渔业经济的综合指标值省市天津河北辽宁上海江苏浙江综合值0.52630.58960.20520.40640.33540.0829省市福建山东广东广西海南综合值0.0276-0.39900.02240.63670.7256引自陈新军

16、（ 2001）3，渔区类型划分根据沿海省市 11 个省市样点的分布状况及其综合指标值，其渔业经济发展可以划分为四个区域， i 类为山东省； ii 类为广东、福建、浙江、辽宁四个省，iii 类有江苏、上海、广西、河北、天津5 个省市； iv 类为海南省（表4-4）。-1-0.8-0.6-0.4-0.200.20.40.60.81-1-0.8-0.6-0.4-0.200.20.40.60.81xy? o ? t? o ? -? ? ? ? -? ? ? ? ?1 ? ? ?1 ? ? o ? ?引自陈新军（ 2001）图 4-6 各省市的综合指标值灰色聚类结果表 4-4 渔业经济区域的划分结果类型

17、省市z 值类型省市z 值类型省市z 值i 类山东0.684 ii 类广东0.514 iii 类江苏0.356 ii 类福建0.496 iii 类上海0.308 ii 类浙江0.465 iii 类河北0.253 类型省市z 值ii 类辽宁0.441 iii 类天津0.253 iv 类海南0.163 iii 类广西0.241 引自陈新军（ 2001）4，结论与分析（1）从综合指标值的大小可以看出，我国11个沿海省市的渔业经济发展水平不一致，根据其指标值的不同和聚类分析的结果，发展水平从高到低的依次顺序分别为山东、广东、福建、浙江、江苏、上海、河北、天津和海南。大致可以分为四个发展档次，渔业发展属于

18、发达的仅为山东省，较为发达的有广东、福建、浙江、辽宁，中等发达的有江苏、上海、河北、天津，不发达的仅为海南。（2）山东省是一个渔业大省，海洋渔业发达。1997年水产品总产量达到了610万吨，居全国首位，占全国水产品总产量的16.9%，其中海洋捕捞的产量和海水养殖产量分别为 297和 238万吨，占全国海洋捕捞产量的18.9%和海水养殖产量的30.1%，均居全国首位； 1997年渔业总产值达到3338069万元，渔民人均收入为9500元/人。山东省充分发展养殖业，海水和淡水的养殖水平利用率分别达到51.0%和 51.8%，海水养殖的单产达到1305.1t/km2。远洋渔业也有一定的发展，1997

19、年为 16.6万吨，约占全国远洋渔业总产量的16.1%。（3）东、福建、浙江和辽宁四个省的渔业也较为发达，是重要的渔业省市。其水产品总产量分别居全国的第25位，在 280520万吨之间，其中海洋捕捞产量都在150万吨以上，它们的总产量和海洋捕捞产量约占全国的44.7%和 59.6%。它们的渔业总产值在 130330亿元之间。除浙江省外，广东、福建和辽宁的海水养殖业也有很大的发展，年产量在100万吨以上。海水和淡水的养殖水面也得到了一定利用，特别是福建省的养殖业更为发达，淡水和海水的养殖单产分别达到447.2t/km2和 1798.3 t/km2，水面利用率分别为73.88%和 56.58%。远

20、洋渔业也处在全国的前列，产量在 620万吨之间。（4）江苏、上海、河北、天津和广西五个省市的渔业经济处在全国的中等和中等偏下水平，海洋渔业不作为主要的发展产业。除江苏和广西的水产品产量在200万吨以上外，其余三个省市的产量都在2060万吨之间，但海洋捕捞产量普遍较低，最高产量（广西）为78 万吨，江苏为 71 万吨，最低（天津）不到3 万吨。江苏、天津和广西的淡水养殖业较为发达，其产量占总产量的一半以上，如江苏省淡水的养殖产量达到 150 多万吨。除上海市外（为7.3万吨），其它四省市的远洋渔业都不够发达，产量基本上都在1 万吨以下。（5）海南是我国的第二大岛，发展渔业的自然条件优越。但是由于

21、发展渔业的经济和技术基础较差，其渔业发展处在全国的落后水平。其水产品总产量仅为52万多吨，为全国倒数第三，并且75.8%的产量来自于海洋捕捞，海水养殖的产量只占5.7%。远洋渔业的发展还没有起步，产量只有5吨。海水养殖业的发展还没有很好的展开，海水的水面利用率只有9.47%。（6）当前，我国近海渔业遭到前所未有的压力，中日渔业协定已经生效，中韩、中越的新渔业协定也正在谈判中，大量的渔民和渔船将从传统的渔场撤回。同时，近海渔业资源的衰退，海洋环境污染和赤潮的增多，这些都严重制约着沿海渔业经济的可持续发展。（7）通过聚类结果和上述分析，我们可以清楚地认识到目前我国沿海省市渔业发展所处的地位、存在的

22、问题及潜力。因此，各地必须根据各自的自然资源条件和渔业经济的发展水平，针对性地制定发展规划，为我国水产业的可持续发展提供保障。二、鱼类营养价值评价方面的应用谢骏、潘厚军在 1992年第 19期的珠江水产上发表了几种名优鱼类营养价值的灰色系统综合评价文章。在该文章中，作者采用了能代表鱼类的营养价值的蛋白质含量、脂肪含量、蛋白质的必需氨基酸指数eaai 及呈味氨基酸含量等四参数，利用灰色系统综合评价方法对鳗鲡、红尼罗、斑点叉尾鮰、淡水白鲳、南方大口鲶和草鱼几种名优鱼类营养价值进行评价，以比较其营养价值。作者在聚类分析之前首先根据近30 种淡水鱼类的营养成分统计，初步拟定评价标准如表 4-5。其具体

23、计算步骤如下：表 4-5 我国淡水鱼类营养价值评价标准营养分级评价参数低中高蛋白质含量（ %）1015 20脂肪含量（ %） 1 3 10eaai 2050 80呈味氨基酸含量（%）1020 30引自谢骏、潘厚军（1992）1确定聚类评价因子根据表 4-5 所列项目，聚类分析采用4 个因子，即 x=x1 ，x2，x3，x4= 蛋白质含量，脂肪含量、 eaai 、呈味氨基酸含量 ，营养价值等级分为三个级别，低、中、高分别记为、。表 4-6 六种鱼类营养成分表鱼名项目鳗鲡红尼罗斑点叉尾鮰淡水白鲳南方大口鲶草鱼蛋白质含量（ %）16.51 16.94 16.31 18.75 15.28 17.53

24、脂肪含量（ %）18.31 0.95 2.99 6.68 1.48 1.8 eaai （%）55.43 70.87 83.79 61.54 63.33 42.76 呈味氨基酸含量（%）19.78 29.34 32.62 25.22 19.57 16.77 引自谢骏、潘厚军（1992）2无量纲化计算为了进行比较，需消除表4-5、表 4-6中量纲的影响，故需要进行无量纲化计算，无量纲化采用下列公式计算。其计算结果见表4-7。3确定各类的白化函数根据表 4-8 确定白化函数的阈值jk，fji最大为 1。表 4-8 无量纲化后六种鱼营养成分及标准鱼名鳗鲡红尼罗斑点叉尾鮰淡水白鲳南方大口鲶草鱼i ii

25、iii 项目蛋白质0.826 0.847 0.816 0.938 0.764 0.877 0.5 0.75 1 脂肪1 0.052 0.163 0.365 0.081 0.098 0.055 0.164 0.546 eaai 0.662 0.846 1 0.735 0.756 0.51 0.239 0.597 0.955 呈味氨基酸0.606 0.899 1 0.773 0.599 0.514 0.307 0.613 0.919 引自谢骏、潘厚军（1992）4求各项的聚类权如对于蛋白质含量而言，低营养级权为：11=同样可以求得其他各项各类别的权。5计算聚类系数令第 i 个样品第 k 类的系数为

26、=如对鳗鲡而言，其低营养级的聚类系数是1111=f11(d1111+ f21(d1221+ f31(d1331+ f41(d1441=0+0+0.606 0.118 （1-0.977）+0=0.002 其中 fjk是根据白化函数内插而得。同理可求出鳗鲡关于其它营养级别的聚类系数。12=0.246 13=0.445 6构造聚类向量，按最大原则进行归类：如对鳗鲡而言，其聚类向量11=（0.002，0.246，0.445）13=0.445为最大，故鳗鲡属于高营养级别。各种鱼的聚类系数jk及聚类结果如表 4-8：表 4-8 各样品聚类系数及聚类结果鱼名等级鳗鲡红尼罗斑点叉尾鮰淡水白鲳南方大口鲶草鱼i

27、0.002 0.034 0.088 0 0.044 0.08 ii 0.246 0.113 0.091 0.254 0.253 0.189 iii 0.445 0.374 0.473 0.361 0.079 0.089 聚类结查iii iii iii iii ii ii 引自谢骏、潘厚军（1992）聚类结果表明，除南方大口鲶和鲩鱼属于中营养级外，其余各鱼均属高营养级。由高到低：斑点叉尾鮰（3=0.473 ）鳗鲡（1=0.445 ）红尼罗鱼（2=0.374 ）淡水白鲳（4=0.361 ）南方大口鲶（5=0.253）鲩鱼（6=0.189 ）。聚类结果还表明，五种名优鱼类的营养价值均高过鲩鱼，同属中

28、级的南方大口鲶亦高于鲩鱼。在五种名优鱼类中，斑点叉尾鮰最高（3=0.473 ），鳗鲡较之略低（1=0.445）。红尼罗鱼和淡水白鲳较接近（2=0.374 、4=0.361 ）。另外，根据调查，国内目前水产品市场地销售价格基本如下：鳗鲡斑点叉尾鮰红尼罗鱼淡水白鲳南方大口鲶鲩鱼。说明鱼类市场价值主要取决于营养价值。但鳗鲡的价格高于斑点叉尾鮰可能有以下几个因素日本对鳗鲡需求量大，而且前又无法人工繁殖；不同测定者的测定条件不同也可能造成不一致；每种鱼的大小、不同的喂养方法，不同的测定季节都可能造成一些误差。鱼类的营养价值除以上几成分外，还包括多种维生素和无机盐等，因此，在今后鱼类营养价值评价中，应该选

29、用更多的因素，以便国家符合实际情况。三、在渔业水域环境评价方面（一）湖泊水质富营养化评价杨红在 1995第 6 期的渔业机械仪器发表了灰色聚类法在湖泊水质富营养化评价中的应用，文章选用阳澄东湖有代表性的八个测站（以测站18表示）的六个聚类指标： do，sd，cod，tp，nh3-n，chla的监测数据。并根据阳澄东湖水质的实际情况，将其富营养化程度分为四级，分别（见表4-9）。文章首先对有关指标进行了均值化的无量纲处理（表410、411）。表 4-9 阳澄东胡水质富营养化分级标准营养类型分级污染指标级贫营养（a0-a1）级中营养（ a1-a2）级富营养（a2-a3）级极富营养（a3-a4）sj

30、o do（/l）15-8.5 8.5-5 5-2 2-0 5 sd（m）13.83-2.94 2.94-1.05 1.05-0.38 0.38-0.13 1.05 cod（mg/l）0.96-2.964 2.964-6.285 6.285-13.325 13.325-28.254 6.285 tp（mg/l）0.003-0.055 0.055-0.186 0.186-0.627 0.627-2.106 0.186 nh3-n（mg/l）0.009-0.055 0.055-0.186 0.186-0.627 0.627-2.106 0.186 chla（mg/m3）0.1-1.59 1.59-1

31、0 10-64 64-400 10.00 引自杨红（ 1995）表 4-10 不同污染指标的灰类值营养类型分级污染指标级级中营养级富营养级极富营养s do（/l）11.75 6.75 3.5 1 5 sd（m）8.385 1.995 0.715 0.255 1.05 cod（mg/l）1.962 4.625 9.805 20.700 6.285 tp（mg/l）0.011 0.042 0.065 0.065 0.065 nh3-n（mg/l）0.065 0.855 0.855 0.855 0.855 chla （mg/m3）0.855 0.855 37 232 0.855 引自杨红（ 1995

32、）利用公式 wjk=来确定聚类权重（见表4-12）。表 4-11 无量纲后的灰类营养类型分级污染指标do 2.35 1.35 0.7 0.2 sd 0.855 1.9 0.243 0.243 cod 0.312 0.736 1.56 3.294 tp 0.169 0.646 2.215 7.6 nh3-n 2.215 2.215 2.188 2.349 chla 0.086 0.580 3.7 23.2 引自杨红（ 1995）表 4-12 污染指标对应不同灰类的聚类权营养类型分级污染指标do 0.212 0. 230 0.063 0.005 sd 0.721 0.324 0.062 0.006

33、 cod 0.028 0.126 0.141 0.079 tp 0.015 0.11 0.201 0.18 nh3-n 0.016 0.111 0.198 0.175 chla 0.008 0.099 0.335 0.554 引自杨红（ 1995）其具体的计算步骤如下：1，污染指标聚类白化数：（a0a4指 x 的取值范围），其余各污染指标白化数可如上相应得到。各测点不同污染指标的聚类白化数可表示为fjk。表 4-13 测点 1 各污染指标的相应白化数fjk 污染指标do 1 0.469 0 0 sd 0 0.343 1 0.657 cod 0 0.890 1 0.107 tp 0 1 nh3-

34、n 0 测点 i 1 0.032 chla 0 1 1 0 引自杨红（ 1995）其它各测点可如上相应求得（略）。2，求聚类系数利用公式可得表 4-14 中的聚类系数。表 4-14 阳澄东胡水质富营养化灰色聚类系数营养 jk 分级测点级贫营养级中营养级富营养级极富营养各测点营养等级 jk测点 10.2220.6990.7080.192富营养测点 20.2170.6110.7030.233富营养测点 30.2170.5990.7160.294富营养测点 40.1950.19501950.195中营养测点 50.1950.6950.6240.219中营养测点 60.6600.6600.6600.2

35、17中营养测点 70.2190.6300.726富营养测点 80.1070.7260.1070.107中营养其中 jk是按最大原则确定的。从表 4-14 中我们可以看出，观测点1、2、3 和 7 为富营养，观测点4、5、6 和 8为中营养。（二）水质科学管理上方面的应用刘冬启、熊邦喜在2000年第 3 期的水利渔业灰色聚类分析在水库水质科学管理上的应用一文中，选取tp、tn 及 cod 作为评价因子对石河水库的水质进行科学评价，以及时掌握其水质状况。文章采用了以下水质评价标准如表4-15。具体计算步骤如下：表 4-15 水质评价标准评价因子（极贫营养级）（贫营养级）（中营养级）（富营养级）（极

36、富营养级）tp/gl-1 600 tn/mgl-1 4.60 cod/mg l-1 27.10 引自刘冬启、熊邦喜（2000）1无量纲化处理首先对灰类（评价级别）和原始白化数（各因子实测数据）作无量纲化处理，以便于比较和消除量纲的影响。灰类的无量纲化处理公式：ij=aij/aj ，i=1，2，3；j=1，2，3，4，5。公式中 ij表示无量纲化处理后白化数（实测数据）的标准值；ai表示第 i 个因子参考标准。这里取各因子第iv 级标准值作为相应的ai 值。无量纲化水质标准见表4-16。表 4-16 无量纲化水质评价标准评价因子（极贫营养级）（贫营养级）（中营养级）（富营养级）（极富营养级）tp

37、/gl-1 0.009 0.036 0.209 1 6 tn/mg l-1 0.017 0.05 0.258 1 3.833 cod/mg l-1 0.013 0.051 0.254 1 3.817 引自刘冬启、熊邦喜（2000）2求聚类权用如下公式：由此求得各因子关于各级别的权如表4-17。表 4-17 各因子关于各级别的权评价因子（极贫营养级）（贫营养级）（中营养级）（富营养级）（极富营养级）tp/gl-1 0.231 0.263 0.290 0.333 0.440 tn/mg l-1 0.456 0.365 0.358 0.333 0.281 cod/mg l-1 0.333 0.372

38、 0.352 0.333 0.280 引自刘冬启、熊邦喜（2000）3聚类系数的求取第 k 个聚类对象第 j 个级别的聚类系数为：式中： k=1；j=1，2，3。其中的求取方法略去。结果见表4-18。表 4-18 聚类系数及结果(极贫营养级(贫营养级(中营养级(富营养级(极富营养级0 0 0.770 0.223 0 引自刘冬启、熊邦喜（2000）4，评价结果石河水库的聚类横向量为=（0，0，0.770，0.223，0），其中，13=0.770最大，故石河水库属于iii 级，即为中营养型。四，鱼类种群划分方面的应用陈新军（ 2002）在西北太平洋柔鱼种群结构初步分析一文中，利用灰色变权聚类的方法对西北太平洋柔鱼种群的结构进行了初步划分。文章收集了沪渔公司鱿钓船于 1999 年 10-11 月在西北太平洋 165e 以西海域作业渔场随机采集120 尾柔鱼样本。共测定胴长 (ml，鳍长 (q1，鳍宽 (q2，眼径 (y，右 1 腕长(wn1 ，右 2 腕长(wn2 ，右 3 腕长(wn3 ，右 4 腕长(wn4 ，右触腕穗长 (sl9 个形态特征指标值。利用灰色系统理论中灰色变权数的聚类分析方法，来研究西北太平洋柔鱼的种群结构。文章根据提供的原始数据，分别求得其它度量指标与胴长的比值，即q1/ml 、q2/ml、y/ml 、w