基于MATLAB主成分分析的攀枝花市地表水水质评价.pdfVIP

PANZHIHUA SCI—TECHINFORMATION · 环境与科技 · 基于MATLAB主成分分析的攀枝花市 地表水水质评价 林 武 攀枝花市环境监测站 摘 要 根据攀枝花市2008—2010年的实测地表水质监测资料，选取 l5个具有代表性的 污染因子，采用基于MATLAB的主成分分析法对所有监测数据进行了分析。结果表明：从原始 数据 中提取 占总方差 86．33％的4个因子来反映水体的污染程度 ，主要关注指标为氨氮、高锰 酸盐指数、粪大肠茵群、溶解氧、铜和pH，主成分分析方法能真实、准确分析 出河流各水质指标 的内在关系。 关键词 攀枝花市 水质现状 MATLAB 主成分分析法 目前，常用的水质评价方法主要有简单 法 。设有n个环境样本，每一个样本有 m 指数法、分级加权平均法、综合污染指数法、 个指标，一般nm，于是形成样本矩阵： 模糊数学法、普通概率统计法等，上述方法虽 u l2 1 然都有一定的数学基础和理论依据，但水质 X2l X22 X2 X = 系统是一个 由多因子构成的复杂系统，水质 评价受诸多指标因子的影响，从而使上述方 瓦 瓦 法在进行水质评价时表现出一定的局限性。 x表示第 i个样本中第J个指标的值，主 主成分分析方法 (PCA法)是一种将多维因 成分分析的实质是通过多维空间中坐标轴的 子纳入同一系统中进行定量化研究、理论 比 旋转变换，找出P个(一般Pm)主成分。为 较完善的多元统计分析方法，在解决很多实 使样本降维后引起的方差最小，需找出坐标 际问题时取得了较好的效果。本文采用 PCA 轴的最佳旋转。因此需要计算矩阵x的协方 法，运用 MATLAB软件统计工具箱，对攀枝 差矩阵C=(c) 。C是一个实对称矩阵， 花市金沙江河流3个监测断面，雅砻江河流2 由c的特征值 (A)的特征向量 ()可以组成 个监测断面，2008～2010年水质监测指标统 一 个正交方阵L=(2) 。 计值 进行了分析评价。 根据特征值和特征向量的概念有：CL= 一 、 研究方法介绍 A。两边分别右乘矩阵 ，并令 LX=Y，有 主成分分析是一种数学变换方法，从应 CLX=A上 Cl，=AY。这样得到由变量X= 用的角度看 ，主成分分析实际上是把多个指 (。， ，…， )线性组合成新变量 Y=(Y， 标化为少数几个综合指标的一种统计分析方 )，2，…，)，并且 l，也是特征值A对应的特征 34 第37卷 第3期 攀枝花科技与信息 (总第 110期) 向量 ，有 ： 性剂、硫化物等9项指标监测值均在检测限 YI ，l1 ，12 1 附近，所有断面监测结果几乎无变化，而且均 I21 I22 小于地表水环境质量三类标准。因此，在数 Y2 Y= (2) ： 据分析中不列入这8项指标，选择剩余的l4 ● ● ● ● ● 项污染因子作为评价指标，pH (x1)、溶解氧 y 一 般Y^。，Y2，…，Y 和其特征值是按方差 ()(