柴达木盆地盐湖水化学特征的对应聚类分析法.pdf

第2 6 卷第7 期计算机与应用化学V 0 1 2 6，N o 72 0 0 9 年7 月2 8 日C o m p u t e r sa n dA p p l i e dC h e m i s t r yJ u l y，2 0 0 9柴达木盆地盐湖水化学特征的对应聚类分析法李梅兰1，吴启勋2(1 青海大学化工学院，青海，西宁，8 1 0 0 1 6；2 青海民族学院化学系，青海，西宁，8 1 0 0 0 7)摘要：目的：研究柴达木盆地盐湖的水化学特征。方法：对应聚类分析法直接处理柴达木盆地盐湖化学组成分析数据，以深入研究这些盐湖水样本的化学特征与化学组成和样硎来源的关系。结果：将样品点和盐湖化学组成变鼍反映到同一聚类轴的图形上，可实现分类，能够揭示出样品来源与变量之间的内在联系。结论：柴达木盆地盐湖的水化学特征属于硫酸镁亚型和氯化物型。关键词：柴达木盆地；盐湖；水化学特征；对应聚类分析中图分类号：0 6 6 1 1；T Q 0 1 5 9 i0 6-3 9文献标识码：A文章编号：1 0 0 1-4 1 6 0(2 0 0 9)0 7-9 4 3-9 4 7l引言近几年来，一些多元统计数据分析法，如主成分分析、因子分析以及判别分析等相继应用于研究湖水水化学特征等n-l o 。但这些方法的1 个明显局限性就是样本与变量总是分别考虑的。因此，以因子分析为例，2 0 世纪7 0 年代提出了对应分析方法“1 5】。该方法综合了R 型和Q 型因子分析的优点，经过一定的数学变换，把变量和样品的载荷反映在相同的公因子轴上，这样就将变最和样本联系起来便于分析、解释和推断。但变量和样本分类带有随意性。本文尝试将对应分析和聚类分析结合起来应用于研究柴达木盆地盐湖水化学特征，发现柴达木盆地盐湖与其化学组成、分布地之间的关系，为合理解释柴达木盆地盐湖水化学特征提供了科学依据。2 硅应聚类分析法的原理对应聚类分析法的详细数学模型可参见文献 1 6 2 0 。这里简单介绍其原理；设有n 个样本和P 个指标的原始测量数据矩阵x=。且勺 0，其主要计算步骤如下。(1)将原始测量数据矩阵X 按行、列求和，并求出总和L总和r2 蚤荟2；扎2；算。并i 为行和(i=1，2，凡)，茹_ 为列和(J=1，2，p)(2)计算“中心一对称”矩阵Z：乃 r t x p oi=1，2，n；J=l，2，P(3)进行R 型因子分析。A=z 7 Z，取前m 个特征根(贡献率要求大于8 0 以上)，即取m个公因子，计算其相应的单位化特征矢量e。、e：e。，从而得到R 型因子载荷矩阵，记为V。(4)进行Q 型因子分析。可以证明，B=Z Z7 与A 有相同的特征根。因此计算(3)中A 的前m 个特征根计算相应于B矩阵的特征矢量e；=Z e。(i=1，2，m)，并将e+单位化，从而得到Q 型因子载荷矩阵，记为u。r 矿1(5)记分块矩阵F=I，I-瓴)，(i=1，2，P+n；，=LU J1，2，m)，对F 矩阵做Q 型聚类分析。(6)根据具体问题进行分析、解释和推断。3柴达木盆地盐湖水化学特征的对应聚类分析法3 1 数据来源柴达木盆地盐湖水化学特征对应聚类分析研究所需用的原始数据，包括盐湖卤水(含晶间卤水)的化学组成和所属的水化学类型，均引自文献 2 1 。考虑到柴达木盆地盐湖中没有碳酸盐型盐湖，硫酸盐硫酸钠哑型盐湖也只有一例，因此我们随机选取部分硫酸盐硫酸镁亚型盐湖(7 个)和全部氯化物盐型盐湖(4 个)组成对应聚类分析样本(见表1)。3 2 数据预处理化学组成作为对应聚类分析的变虽指标。采用对数变换方法压缩指标后，再采用指标变量除以其标准差的方法，无量纲化处理数据，以保证结论准确。收稿日期：2 0 0 8 1 2-2 5；修回日期：2 0 0 9 0 4-2 8作者简介：李梅兰(1 9 6 9 一)，女，青海人讲师，E m a i l：y o u j i e z i 2 0 0 8 y a h o o c n 等J J勺万方数据计算机与应用化学2 0 0 9，2 6(7)表l 柴达木盆地盐湖的化学组成(m s L)T a b l elC h e m i c a lc o m p o n e n to fs M tl a k e si nQ a r d B a s i n(n【l g L)茶卡盐湖c a k as a l tl a k e柯柯盐湖k e k es a l tl a k e尕海s m a l ls e a 唧e s i m大柴达木湖l a r g eq a r d a ml a k e小柴达木湖B m a q a r d a ml a k e尕斯库勒湖g a s k u hl a k昆特依k u n t i y i钾湖p o t a s s i u ml a k e德宗马海D e z o n g M as e a察尔汁c h a e r h a n东陵湖d o n n i n gl a k e硫酸镁亚型m a g n e s i u ms u l f a t es u b t y p e硫酸镁哑型m a g n e s i u ms u l f a t es u b t y p e硫酸镁亚型s u l f a t es u b t y p e硫酸镁亚型m a g n e s i u ms u l f a t es u b t y p e硫酸镁哑型唧e s i u l ns u l f a t es u b t y p e硫酸镁亚型m a g n e s i u ms u l f a t es u b t y p e氯化物盐型c h l o r i d es a l t氯化物盐型c h l o r i d es a l t硫酸镁亚型m a g n e s i u ms u l f a t es u b t y p e氯化物盐型c h l o r i d es a l t氯化物盐型c h l o r i d es a l t8 02 3 144 7 32 65 0 61 2 4 51 8 77 0 52 36 2 51 7 2 72 6 58 00 4 446 1 82 73 9 62 3 1 01 9 05 1 92 34 3 31 3 1 02 73 3 32 5 08 83 8 732 2 241 7 63 0 4 64 52 6 01 30 8 81 4 1 93 8 O3 3 996 9 74 5 3 01 5 58 9 21 65 6 71 2 2 01 7 5 01 0 6 3 8 635 2 01 34 1 44 1 7 01 83 4 63 l5 8 42 4 7 O0 O7 69 4 845 1 12 95 9 24 0 6 01 7 57 5 74 52 3 54 7 1 01 9 0 05 37 6 51 22 1 93 79 1 772 8 6 02 1 68 5 98 3 540 9 653 6 79 02 1 31 08 2 0 O2 9 3 4 0 64 26 56 8 O2 7 11 3 6 61 0 r 74 3 l3l“1 56 8 91 7 5 11 7 54 2 05 28 4 94 2 6 43 8 6 97 13 6 01 2 儿O2 86 6 7l0 1 7 02 0 l5 5 564 1 42 09 44 2 05 l8 6 0 03 6 5 0 8 02 0O 03 3 2 20 0O O3 3 对应聚类分析计算结果贡献比例达到9 3 4 6)，Q 型因子分析计算结果见表3。对应聚类分析算法采用M A T L A B 编程2 2。拍，已在P c图1 划出1 9 个指标与样品Q 型系统的聚类结果，是以欧氏机(英特尔奔腾双核1 8 6G 处理器)调试通过。距离作为聚类样品的测度，类间采用离差平方和(w a r d)算法运算：R 型因子分析计算结果见表2(取m=3，特征值做出分类的m 功1。表2 主因子的R-型因子载荷T a b l e2C a p a c i t yf o rl o a do fR t y p ef a c t o ro fp r i m a r yf a c t o r a 指标变量特征值诣孺慧尝毡抽,n u o h“,8thedicatorv a r i a b l e si g e嚣：=1 荔第1 主因子载荷f i r s tp r i n c i p a lf a c t o rl o a d i n g第2 主因子载荷s e c o n dp r i n c i p a lf a c t o rl o a d i n g第3 主因子载荷t h i r dp r i n c i p a lf a c t o rl o a d i n g1Na+00263 3 36 1 3 76 1 3 70 0 3 32 6 70 0 l l5 9 20 0 2 32K+00101 4 72 3 6 58 5 0 20 0 3 0 0 70 0 2 46 7 40 0 2 48 8 73M一+00036 2 068 4 49 3 4 60 0 4 69 60 0 1 22 0 20 0 0 39 8 724Ca2+00020 8 564 8 69 8 3 20 0 8 17 0 20 0 0 51 3 010 0 0 87 6 265C100006 0 15 51 4 09 9 7 20 0 3 60 r 7 70 0 0 04 0 45 8一O 0 l l1 4 26S 嘶一0 0 0 01 1 87 50 2 09 9 9 20 0 6 90 2 20 0 1 73 0 60 0 1 97 2 57C O；一0 0 0 00 9 2O 0 0 08l o o 0 00 0 8 76 4 70 0 4 14 6 40 0 3 83 6 78H C O f0 0 0 00 00 0 0 001 0 0 0 00 0 4 1 1 60 0 8 49 0 80 0 1 98 5 7表3 主因子的Q 一型因子载荷T a b l e3C a p a c i t yf o rl o a do fQ-t y p ef a c t o ro fp r i m a r yf a c t o r s 第l 主因子荷载第2 主因子荷载第3 主因子荷载第1 主因子荷载第2 主因子荷载第3 主因子荷载N o f i r s tp r i n c i p a ls e c o n dp r i n c i p a lt h i r dp r i n c i p a lN o f i r s tp r i n c i p a ls e c o n dp r i n c i p a lt h i r dp r i n c i p a lf a c t o rl o a d i n gf a c t o rl o a d i n gf a c t o rl o a d i n gf a c t o rl o a d i n gf a c t o rl o a d i n gf a c t o rl o a d i n g10 0 2 80 7 60 0 1 49 90 0 0 42 2 6770 0 3 20 2 50 0 2 24 7 50 0 1 54 5 220 0 2 33 2 70 0 0 88 7 610 0 0 39 2 4780 0 4 42 7 10 0 2 71 0 6一O 0 2 75 630 0 2 47 8 60 0 1 46 3 30 0 1 21 6 890 0 4 84 2 90 0 1 25 9 50 0 1 56 9 640 0 2 79 5 80 0 1 23 4 10 0 0 52 0 3 71 00 0 1 91 9 3o 0 5 28 2 50 0 3 59 9 250 0 2 28 6 30 0 6 40 7 60 0 2 40 9 91 l一0 1 2 62 60 0 3 42 6 40 0 0 78 2 2l60 0 3 79 7 40 0 0 28 5 22 0 0 1 51 7万方数据2 0 0 9，2 6(7)李梅兰，等：柴达木盆地盐湖水化学特征的对应聚类分析法9 4 5F i 辱lC o r r e p o n d e n c ec l u s t e ro f8v a r i a b l e sa n d11s a m p l e so fs a l tl a k e s 图18 个变量与1 1 个盐湖样品的对应聚类(1 8 为变量，9 1 9 为盐湖样本)4 讨论4 1 关于对应聚类分析中的量纲问题对应聚类分析中的量纲问题，许多专门的多元统计分析著作均未做出明确的表述。本文认为，由原始数据出发分析R 型因子时的特征值有量纲，而Q 型因子无量纲(对应样本，当然无量纲)，2 种分析并不对称。这样，同一物理量用不同的量纲，其数量值有变化，对应分析的结果也可能有变化，将影响最后聚类结论。为此，采用指标变量除以其标准差的方法，处理原始数据无量纲化，以保证结论的准确性。显然，对应聚类分析采用通常的中心化或标准化法处理无量纲化是不行的。因这样处理原始数据后，其列和将等于零，使中心对称矩阵z 变得无意义，也就无法对应聚类分析了。4 2 关于盐湖样本容量随机选取部分硫酸盐硫酸镁亚型盐湖(7 个)和全部氯化物盐型盐湖(4 个)共计1 1 个盐湖组成对应聚类分析样本，是出于统计学上的考虑。有关技术规范指出mJ，研究不同产地样本的某种统计规律，其样本至少1 0 个，才能保证分析结果具有统计学意义。4 3 柴达木盆地盐湖与其化学组成、分布地之间的相关关系对应聚类分析把R 型和Q 型因子分析结合起来，以得出的指标和样本的3 个因子载荷作为聚类分析变量，同时聚类分析。因此，在指标和样本分类的基础上，能够揭示出样本来源与变量之间的内在联系，从而找出柴达木盆地盐湖与其化学组成、分布地(样本来源)之间的关系。图l 中8 个化学组成指标(图l 中1、2、3、4、5、6、7、8 依次对应表1 的8 个化学组成)和1 1 个盐湖样本(图l 中9、1 0、1 1、1 2、1 3、1 4、1 5、1 6、1 7、1 8、1 9 依次对应表1 中1 1 个盐湖样本)可以分成2类。位于图l 左边的指标变量1(N a+)、2(K+)、6(s 暖一)和9(1，表l 中盐湖样本的编号，以下同)1 0(2)、1 1(3)、1 2(4)、1 3(5)、1 4(6)、1 7(9)盐湖样本聚为一类，说明茶卡盐湖、柯柯盐湖、尕海、大柴达木湖、小柴达木湖、尕斯库勒湖、德宗马海等7 个硫酸镁亚型盐湖主要与l(N a+)、2(K+)、6(s o l 一)等3 个化学组成相关。也可以看出指标变量7(c o；一)和8(H C 0 3-)只有在入值较大时才与前些指标变量和样本聚为一类，表明它们与盐湖的水化学特征关系数不大。指标变量3(M 9 2+)、4(c a“)和5(C 1)与1 5(7)、1 6(8)、1 8(1 0)和1 9(1 1)号等盐湖样本聚为一类，表明化学组成3(M 9 2+)、4(c a 2+)和5(c l)对昆特依、钾湖、察尔汗和东陵湖等氯化物盐型盐湖样本水化学特征影响较大。以上揭示的柴达木盆地盐湖与其化学组成、分布地之间的相关关系，与这些盐湖的实际情况相符。4 4 柴达木盆地盐湖水化学特征解释柴达木盆地硫酸镁亚型盐湖主要与1(N a+)、2(K+)和6(s o l。)等3 个化学组成相关，它们正是形成硫酸盐型盐湖的显示和条件；而柴达木盆地氯化物盐型盐湖主要与3(M 9 2+)、4(c a 2+)和5(c l)等3 个化学组成相关，同样也是形成氯化物盐型盐湖的显示和条件，而且显示柴达木盆地氯化物盐型盐湖是以高含量M g 为其特征；化学组成7(C O；)和8(H C O；)与其盐湖水化学特征关系不大，表明柴达木盆地盐湖发育程度较深，没有碳酸盐类型的水化学特征存在。这些结论是与实际情况相符。5结语采用对应聚类分析法研究盐湖水化学特征，可以将指标变量化学组成和盐湖样本联系起来。在找出样本来源与指标变量之间的内在联系的基础上，可以揭示柴达木盆地盐湖与其化学组成、分布地之间的相关关系，为合理解释柴达木盆地盐湖水化学特征提供科学依据。R e f e r e n c e s1L i a oQ i z h ia n dL uX i a o h u a T h ef a c t o ra n a l y s i so fe n v i r o n m e n t a lc h e m i c a lc h a r a c t e r i s t i c si ns a n d a nl a k e JH u a z h o n gU n i vo fS c i T e c h，1 9 9 8，2 6(1 1)：6 5 6 7 2L iL e i W uQ i x u na n dS o n gP i n g M e t h o d so fd i s e r i m i n a n ta n a l y s i sc l a s s i f i c a t i o no fh y d r o e h e m i c a lt y p e so fs a l tl a k e s J o u r n a lo fS a l tL a k eR e s e a r c h，2 0 0 4，1 2(1)：4 6 5 0 3Z h uX i a o j u a na n dP uZ h i x i s o Af a c t o ra n a l y s i so ft h ec h e m i c a lc h a r a c t e r i s t i c so ft h ew a t e ri nx i f l g I l ul a k e T r a n s a c t i o n so fO c e a n o l-o g ya n dL i m n o l o g y，2 0 0 5(1)：6 一1 0 4A nY a h M e t h o d so ff a c t o ra n a l y s i so fh y d r o c h e m i c a lc h a r a c t e r i s t i c so fQ i n g h a iS a l tL a k e s J o u r n a lo fO i n g h 且iN o r m a lU n i v e r s i t y：N a t u r a lS c i e n c e，2 0 0 4(3)：6 1 6 3 5W uQ i x 岫，L iL e ia n dA nY 舳M e t h o d so fm u l t i v a r i a t es t a t i s t i c a lc l a s s i f i c a t i o no fh y d r o c h e m i c a lt y p e so fC h i n ss a l tl a k e s C h i n e s eJ o u r n a lo fA n a l y t i c a lC h e m i s t r y，2 0 0 5，3 3(3)：4 3 3 6C h e nZ h i j u na n d 盔删唱J i n g S t u d yO i lf a c t o ra n a l y s i so fh y d r o c h e m-i s t r yc h a r a c t e r i s t i c so fE b i n u rL a k e J o u r n a lo fS a l tL a k eR e s e a r c h，2 0 0 8，1 6(2)：1 9 2 1 7Z h o uJ i q i n。P a nX i a o l u，Y i nC a i x i aa n dP e n gL i T h ef e a t u r ea n a l y-g i Bo ft r a c ee l e m e n t a Ji nw a t e r-s o l u b l es u b s t a n c eo f t h ec o m p o s i t ef a m i l y C o m p u t e r sa n dA p p l i e dC h e n I i s t r y，2 0 0 0，1 7(6)：5 4 1 5 4 5 8O iJ u n s h e n g，X uH u l b i，Z h o uJ i n g y a n，L aX i a n h u aa n dG u a nJ i n g-525l50吡仉叭n万方数据计算机与应用化学2 0 0 9，2 6(7)h u a S t u d i e so nt h em o u n to ft r a c ee l e m e n t sa n de f t 3 c a c yi nC h i。n e s em e d i c i n a lh e r b sf o rt r e a t i n ge x t e r i o rs y n d r o m e s C o m p u t e r sa n dA p p l i e dC h e m i s t r y，2 0 0 3，2 0(4)：4 4 9 4 5 2 9Y i nC a i x i a，Z h o uJ i q i na n dP e n gL i，n l e 州n c i p a lc o m p o n e n ta n a l y s i so f 仃a c ee l e m e n t a lc o n t e n to np l a n t sf r o mw e s t e r nh i l lo fK u n-n I i l l g T r a c eE l e m e n t sS c i e n c e，1 9 9 9，6(1 1)：2 9 3 2，1 0Y i nC a i x i a，P e n gL ia n dZ h o uJ i q i n，e ta 1 S t u d i e so nm a t h e m a t i c a lm e d i e i n a l c o m p o s i t e sp l a n sf r o mY u n n a n J o u r n a lo fB i o m a t h e m a t i c s，1 9 9 8，1 3(6)：1 0 8 8 1 0 9 3 11D a d d mM C a m p i g l i oCa n dD r l i n gR P r o g r e s si nR-a n dO m o d ea n a l y s i so r r e s p o n d e n c ea n a l y s i sa n di t sa p p l i c a t i o nt ot h es t u d yo fg e o l o g i c a lp r o c e s s e s JE a r t hS e i，1 9 7 4，1 1：1 3 1 1 2G r e e n a c r eM】T h e o r ya n dA p p l i c a t i o no fC o r r e s p o n d e n c eA n a l y s i s N e wY o r k：A c a d e m i cP r e s s。1 9 8 4 1 3X i aL Z h o uXa n dY uL O nt h em a t h e m a t i c a lm o d e lo fc o r r e s p o n-d e n c ea n a l y s i s C h i n a-J a p a ns y m p o s i u mo ns t a t i s t i c s B 嘶i D g，C h i-n a，1 9 8 4：3 4 9 3 5 2 1 4X i aLa n dY a n gY Am e t h o do fv a r i a b l es d e c t i o ni nH a y a s h i st l l i r dm e t h o do fq u a n t i f i c a t i o n O u m a lo ft h eJ a p a n e s eS o c i e t yo fC o m p u r a t i o n a lS t a t i s t i c s，1 9 8 8，l(1)：2 7 一“1 5A n n aTa n dM i c h a e lG D u a ls c a l i n ga n dc o r r e s p o n d e n c ea n a l y s i so fp r e f e r e n c e s，p a i r e dc o m p a r i s o n sa n dr a t i n g s I n t e r n a t i o n a lJ o u n r a lo fR e s e a r c hi nM a r k e t i n g，2 0 0 2，1 9(4)：4 0 1 4 0 5 1 6T a oFa n dX i aLX S t e p w i s ec o r r e s p o n d e n c ea n a l y s i sa n ds t u d yo fp h y s i c a lc h a r a c t e r i s t i c s o fC h i n e s ee t h n i cg r o u p s J o u r n a lo fB i o nA t h a nA r i e s，2 0 0 3，1 8(2)：1 3 9 1 7T a oF C h e nLa n dX i aL C o r r e s p o n d e n c ea n a l y s i so fb o d yf r o mc h a r a c t e r i s t i c so fC h i n e s ee t h n i cg r o u p s D i s e r e ta n dC o n t i n u o u sD y-n a m i c a lS y s t e m sS e r i e sB，2 0 0 4。4(3)：7 6 9 7 7 6 1 8C r e e n a c r eMJ C r r e s p o n d e n c eA n a l y s i si nP r a c t i c e L o n d o n：A c a d e m i cP r e s s，1 9 9 3 1 9B e hEJ Ac o m p a r a t i v es t u d yo f K!o l-e 8f o rc o r r e s p o n d e n c ea n a l y s i sw i t l lo r d e r e dc a t e g o r i e s B i o m e t r i c a lJ o u m a l，1 9 9 8，4 0(4)：4 1 3 4 2 9 2 0W uH a i l o n g，z e l I gB e i w e ia n dY uR u q i n C o r r e s p o n d e n c ef a c t o ra n dc o r r e s p o n d e n c ec l u s t e ra n a l y s i s a p p l i e dt os t u d i e so nt h eq u a l i t yo fC h i n e s et e a C h i n e s eJ o u r n a lo fA n a l y t i c a lC h e m i s t r y，1 9 9 I 1 9(4)：4 5 6 4 5 9 2 1Z h a n gP e n g x i S a l tL a k e si nQ a r d mB a s i n B e i j i n g：S c i e n c eP r e s s，1 9 8 7：8 1 9 1 2 2G a oJ u n b i n M A T L A B5 0L a n g u a g ea n da n dP r o g r a mD e s i g n W u-h a n：H u a z h o n gU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g yP r e s s，1 9 9 9 2 3W a n gM o r a na n dX i a oJ i n g s o n g M A T L A B5 xy uK e x u eJ i s u a n B e i j i n g：T s i n g h u aU n i v e r s i t yP r e s s，2 0 0 0：1 2-6 5 2 4W uXLa n d i nZH D e s i g no fA s s i s t a n tF u z z yS y s t e mo fM A T-L A B X i a n：P u b l i s h i n gC o m p a n yo fX i a nE l e c t r o n i cS c i e n c ea n dT e c h n o l o g yU n i v e r s i t y，2 0 0 2：1 3 1 1 5 8 2 5C l e v eM o l e r N u m e r i c a lC o m p u t i n gw i t hM A T L A B N e wY o r k：t h eS o c i e t yf o r I n d u s t r i a la n dA p p l i e dM a t h e m a t i c s，2 0 0 4 2 6M a t h 眦s S t a t i s t i c sT o o l b o xU s e r 8G u i d e 2 0 0 4 2 7R u s p i n iEH An e wa p p r o a c ht oc l u s t e r i n g I n f o r m a t i o na n dC o n-t r o i。1 9 6 9，1 5：2 2 3 2 2 8J iH o n g i i n，Z h uY o n g z h e n ga n dW uX i s h e n g C o r r e s p o n d e n c eC l u s-t e rA n a l y s i sa n dI t sA p p l i c a t i o ni nE x p l o r a t i o nG e o c h e m s t r y JG e o c h e mE x p l o r，1 9 9 5，5 5：1 3 7 1 4 4 2 9X i n gY a n w e i，W a n gJ i ea n dZ h o n gJ i n b e i S t u d yo nr u l eo fs y n-d r o m ec o m b i n a t i o no f1 0 6 9c a B e 8o fa n g i n ap e e t o r i sb yc l u s t e ra n a l y s i s a n dc o r r e s p o n d i n g-c o r r e l a t i o na n a l y s i s Z h o n g h u aZ h o n gY iY a oz az I I i，2 0 0 7，2 2(1 1)：7 4 7 7 5 0 3 0S t a t eF o o da n dD r u gA d m i n i s t r a t i o n T e c h n i c a lr e q u i r e m e n t so ft h ef i n g e r p r i n ti ni n j e c t i o no fc h i n e s em a t e r i am e d i c s(T e n t a t i v eS t a n d a r d)C h i n e s eT r a d i t i o n a lP a t e n tM e d i c i n e，2 0 0 0，2 2(1 0)：6 7 1 6 7 5 中文参考文献1廖奇志，陆晓华山丹湖环境水化学特征的主因子分析 J 华中理工大学学报1 9 9 8，2 6(1 1)：6 5 6 7 2李磊，吴启勋，宋萍盐湖水化学类型的判别分析研究 J 盐湖研究，2 0 0 4，1 2(1)：4 6 5 0 3朱小娟普智晓星湖水环境化学特征的主因子分析 J 海洋湖沼通报，2 0 0 5(1)：6 1 0 4安燕青海盐湖水化学特征的因子分析研究 J 青海师范大学学报(自然科学版)，2 0 0 4(3)：6 l 一6 3 5吴启勋，张世芝，安燕中国盐湖水化学类型的多元统计分类法 J 分析化学，2 0 0 5，3 3(3)：4 3 3 6陈志军，张晶艾比湖水化学特征的因子分析研究 J 盐湖研究，2 0 0 8。1 6(2)：1 9 2 1 7周纪勤，潘晓露，殷彩霞，彭莉菊科植物水提物微量元素特征分析 J 计算机与应用化学2 0 0 0，1 7(6)：5 4 1-5 4 5 8祁俊生，徐辉碧，周井炎，陆晓华，管竞环解表植物类中药中微量元素与功效关系 J 计算机与应用化学，2 0 0 3，2 0(4)：4 4 9 4 5 2 9殷彩霞，周纪勤，彭莉昆明西山植物微量元素主成分分析 J 微量元素科学，1 9 9 9，6(1 1)：2 9-3 2 1 0 殷彩霞，彭莉，周纪勤，等云南菊科药用植物微量元素数学建模研究 J 生物数学学报，1 9 9 8，1 3(6)：1 0 8 8 1 0 9 3 2 0 吴海龙，曾北危，俞汝勤对应分析用于中国茶叶化学数据的综合分析 J 分析化学，1 9 9 1，1 9(4)：4 5 64 5 9 2 1 张彭熹柴达木盆地盐湖 M 北京：科学出版社，1 9 8 7：8 1 9 1 2 2 高俊斌M A T L A B5 0 语言与程序设计 M 武汉：华中理工大学出版社，1 9 9 9 2 3 王沫然，肖劲松M A T L A B 5 X 与科学计算 M 北京：清华大学出版社，2 0 0 0：1 2 6 5 2 4 吴晓莉，林哲辉M A T L A B 辅助模糊系统设计 M 西安：西安电子科技大学出版社，2 0 0 2：1 3 1 1 5 8 2 9 邢雁伟，王阶，衷敬柏，等采用聚类分析和对应相关方法研究l0 6 9 例冠心病心绞痛证候应证组合规律 J 中华中医药杂志。2 0 0 7，2 2(1 1)：7 4 7 7 5 0 3 0国家食品药品监督管理局中药注射剂色谱指纹图谱研究操作规程指南(试行)J 中成药，2 0 0 0，2 2(1 0)：6 7 1 6 7 5 万方数据2 0 0 9，2 6(7)李梅兰，等：柴达木盆地盐湖水化学特征的对应聚类分析法9 4 7C o r r e s p o n d e n c ec l u s t e ra n a l y s i so fh y d r o c h e m i s t r yc h a r a c t e r i s t i c so fs a l tl a k e si nQ a r d a mb a s