kmeans图像聚类操作

经典资源分配网络RAN程序

24-5-3 轴系5(3)-深沟球轴承寿命计算.m

【实例简介】SNOPT使用说明

【实例简介】根据标准生成随机路谱

双馈风力发电机（DFIG_kq_chuantong.slx）

【实例简介】MATLAB 层次分析法

第一章  误差与范数1.1  误差的来源1.2  误差和有效数字1.3  误差估计的基本方法1.4  数值计算中应注意的问题第二章  非线性方程（组）的数值解法2.1  方程（组）的根及其MATLAB命令2.1.2  求解方程（组）的solve命令 2.1.4  求解方程（组）的fsolve命令2.2  搜索根的方法及其MATLAB程序2.2.1  作图法及其MATLAB程序2.2.2  逐步搜索法及其MATLAB程序2.3  二分法及其MATLAB程序2.3.1  二分法求方程根的迭代次数的MATLAB命令    2.3.2  二分法的MATLAB程序2.4  迭代法及其MATLAB程序2.4.2  迭代法的MATLAB程序1  2.4.5  迭代法的MATLAB程序22.5  迭代过程的加速方法及其MATLAB程序2.5.2  加权迭代的MATLAB程序2.5.4  艾特肯(Aitken)加速方法的MATLAB程序2.6  牛顿(Newton)切线法及其MATLAB程序2.6.2  牛顿切线法的收敛性及其MATLAB程序2.6.3  牛顿切线法的MATLAB程序2.6.4  求2.6.6  牛顿切线法的加速及其两种MATLAB程序2.7  割线法及其MATLAB程序2.7.2  割线法的MATLAB程序2.8  抛物线法及其MATLAB程序2.8.2  抛物线法的MATLAB程序2.9  求解非线性方程组的牛顿法及其MATLAB程序2.9.1  求解二阶非线性方程组的牛顿法及其MATLAB程序2.9.2  求解阶非线性方程组的牛顿法及其MATLAB程序2.9.3  求解阶非线性方程组的拟牛顿法及其MATLAB程序第三章  解线性方程组的直接方法3.1  方程组的逆矩阵解法及其MATLAB程序3.1.3  线性方程组有解的判定条件及其MATLAB程序3.2  三角形方程组的解法及其MATLAB程序3.2.2  解三角形方程组的MATLAB程序3.3  高斯消元法和列主元素消元法及其MATLAB程序3.3.1  高斯消元法及其MATLAB程序3.3.2  列主元素消元法及其MATLAB程序3.4 LU分解法及其MATLAB程序3.4.1  判断矩阵LU分解的充要条件及其MATLAB程序3.4.2  直接LU分解法及其MATLAB程序3.4.3  含交换阵的LU分解法及其MATLAB程序3.4.4  正定对称矩阵的Cholesky分解及其MATLAB程序3.4.5  正定对称矩阵的Cholesky分解及其MATLAB程序3.5  求解线性方程组的方法及其MATLAB实现3.5.1  求解线性方程组的Cholesky分解及其MATLAB实现3.5.2  求解线性方程组的直接分解法及其MATLAB实现3.5.3  解线性方程组的选主元的方法及其MATLAB实现3.6  误差分析及其两种MATLAB程序3.6.1  用MATLAB软件作误差分析3.6.2  求P 条件数和讨论解的性态的MATLAB程序3.6.3  用P范数讨论解和的性态的MATLAB程序第四章  解线性方程组的迭代法4.1  迭代法和敛散性及其MATLAB程序4.1.2  用谱半径判别迭代法产生的迭代序列的敛散性的MATLAB程序4.2  雅可比（Jacobi）迭代及其MATLAB程序4.2.2  判别雅可比迭代收敛性的MATLAB程序4.2.3  雅可比迭代的两种MATLAB程序4.3  高斯—塞德尔（Gauss-Seidel）迭代及其MATLAB程序    4.3.3  高斯—塞德尔迭代两种MATLAB程序4.4  解方程组的超松弛迭代法及其MATLAB程序4.4.2  超松弛迭代法收敛性及其MATLAB程序4.4.3  超松弛迭代法的MATLAB程序第五章  矩阵的特征值与特征向量的计算5.2  幂法及其MATLAB程序5.2.2  幂法的MATLAB程序5.3  反幂法和位移反幂法及其MATLAB程序5.3.3  原点位移反幂法的两种MATLAB程序5.4  雅可比(Jacobi)方法及其MATLAB程序5.4.3  雅可比方法的MATLAB程序5.5  豪斯荷尔德方法及其MATLAB程序5.5.1  豪斯荷尔德方法及其MATLAB程序5.5.2  矩阵约化为上Householder阵及其MATLAB程序5.5.3  实对称矩阵的三对角化及其MATLAB程序5.6  QR方法及其MATLAB程序5.6.2  QR方法的MATLANB程序5.6.5  最末元素位移QR法计算实对称矩阵特征值及其MATLAB程序5.6.6  求根位移QR方法计算实对称矩阵A的特征值及其MATLAB程序5.7  广义特征值问题及其MATLAB程序5.7.2  用MATLAB计算型的广义特征值和特征向量第六章  函数的插值方法6.1  插值问题及其误差6.1.2  与插值有关的MATLAB 函数6.2  拉格朗日（Lagrange）插值及其MATLAB程序6.2.1  线性插值及其MATLAB程序6.2.2  抛物线插值及其MATLAB程序6.2.3  次拉格朗日（Lagrange）插值及其MATLAB程序6.2.5  拉格朗日多项式和基函数的MATLAB程序6.2.6  拉格朗日插值及其误差估计的MATLAB程序6.3  牛顿（Newton）插值及其MATLAB程序6.3.3  牛顿插值多项式、差商和误差公式的MATLAB程序6.3.4  牛顿插值及其误差估计的MATLAB程序6.3.5  牛顿插值法的MATLAB综合程序6.4  埃尔米特（Hermite）插值及其MATLAB程序6.4.3  埃尔米特插值多项式的MATLAB程序6.5  分段插值及其MATLAB程序6.5.1  高次插值的振荡和MATLAB程序6.5.3  分段线性插值的MATLAB程序6.5.4  作有关分段线性插值图形的MATLAB程序6.5.5  用MATLAB计算有关分段线性插值的误差6.6  分段埃尔米特（Hermite）插值及其MATLAB程序6.6.2 Hermite插值的MATLAB程序6.6.3  作有关分段Hermite插值图形的MATLAB程序6.6.4  用MATLAB计算有关分段Hermite插值的误差6.7  三次样条插值及其MATLAB程序6.7.4  用一阶导数计算的几种样条函数和MATLAB程序6.7.6  用MATLAB计算三次样条6.7.7  几种作三次样条有关图像的MATLAB程序6.7.8  用MATLAB计算有关分段三次样条的误差6.8  高元插值及其MATLAB程序6.8.1  MESHGRID命令的功能和调用格式6.8.2  单调数据点上的二元插值及其MATLAB程序6.8.3  三元插值及其MATLAB程序第七章  函数逼近与曲线（面）拟合7.1  曲线拟合、误差及其MATLAB实现7.2  曲线拟合的线性最小二乘法及其MATLAB程序7.3  函数的选取及其MATLAB程序7.4  多项式拟合及其MATLAB程序7.5  拟合曲线的线性变换及其MATLAB程序7.6  函数逼近及其MATLAB程序7.7  三角多项式逼近及其MATLAB程序7.8  随机数据点上的二元拟合及其MATLAB程序7.9  随机数据点上的元拟合及其MATLAB程 第八章  数值微分8.2  一阶导数的数值计算及其MATLAB程序……………………………………1178.2.1  差商求导及其MATLAB程序8.2.2  中心差商公式求导四种MATLAB程序8.2.3 Richardson外推法求导及其MATLAB程序8.2.4  牛顿多项式求导及其MATLAB程序8.2.5 diff和gradient函数在数值求导中的应用8.3  高阶导数的数值计算及其MATLAB程序8.3.1  插值或拟合高阶数值导数及其MATLAB程序8.3.2  高阶泰勒数值导数及其MATLAB程序8.4  数值梯度和数值偏导数的计算及其MATLAB程序8.4.1  梯度和偏导数的数值计算及其MATLAB程序8.4.2  计算雅克比矩阵及其行列式的MATLAB方法第九章  数值积分9.1  积分的符号计算9.1.1  定积分的MATLAB符号计算9.1.2  变限积分的MATLAB符号计算9.2  数值积分的思想及其MATLAB程序9.2.3  矩形公式的MATLAB运算9.3  插值型数值积分及其MATLAB 程序9.3.2  梯形公式的四种MATLAB程序9.3.4  辛普生（Simpson）数值积分的MATLAB程序9.3.6  牛顿-科茨（Newton-Cotes）数值积分和误差分析的三种MATLAB程序9.3.7  利用三次样条求表格型数值积分的MATLAB方法9.3.8  利用拉格朗日插值等方法求表格型数值积分的MATLAB方法9.4  龙贝格公式及其MATLAB程序9.4.2  龙贝格积分的MATLAB程序9.5  自适应积分及其MATLAB程序9.6  高斯（Gauss）型积分公式及其MATLAB程序9.6.2  在上的高斯-- 勒让德积分公式及其MATLAB程序9.6.3  在上的高斯-- 勒让德积分公式及其MATLAB程序9.6.4 Radau和Lobatto积分公式及其MATLAB程序9.7  反常积分的计算及其MATLAB程序9.7.1  无穷积分的符号计算及其MATLAB程序9.7.2  无穷积分的近似计算及其MATLAB程序9.7.3  无界函数反常积分的符号计算的五种MATLAB程序9.7.4  无界函数反常积分的近似计算的两种MATLAB程序9.8  多重积分的计算及其MATLAB程序9.8.1  二重积分的符号计算及其MATLAB程序9.8.2  二重积分的梯形公式及其MATLAB程序9.8.3  矩形域上的辛普生公式及其MATLAB程序9.8.4  一般域上二重积分的数值计算及其MATLAB程序9.8.5  三重积分的计算及其MATLAB程序第十章  常微分方程（组）求解10.1 常微分方程（组）的MATLAB符号求解10.1.1 MATLAB求常微分方程（组）的通解10.1.2  MATLAB求常微分方程（组）的特解10.1.3  线性常微分方程组的MATLAB解法10.3  欧拉方法及其MATLAB程序10.3.2  向前欧拉方法的三种MATLAB程序10.3.4  向后欧拉方法的MATLAB程序10.4  改进的欧拉方法及其MATLAB程序10.4.2  梯形公式的两个MATLAB程序10.4.4  改进的欧拉公式的MATLAB程序10.5  龙格—库塔方法10.5.2  二阶龙格—库塔方法及其MATLAB程序10.5.3  三阶龙格—库塔方法及其MATLAB程序10.5.4  四阶龙格—库塔方法及其MATLAB程序10.5.5  自适应龙格—库塔方法及其MATLAB库函数10.6  线性多步法及其MATLAB程序10.6.2  亚当斯显式公式及其MATLAB程序10.6.3  亚当斯隐式公式及其MATLAB程序10.6.4  米尔恩公式及其MATLAB程序10.6.5  汉明公式及其MATLAB程序10.6.6  预测-校正系统及其MATLAB程序10.7  一阶（高）阶微分方程（组）的数值解及其MATLAB程序10.7.1  一阶微分方程组的数值解及其MATLAB程序10.7.2  高阶微分方程（组）的数值解及其MATLAB实现10.8  边值问题的数值解及其MATLAB程序10.8.1  打靶法及其MATLAB程序10.8.2  有限差分方法及其MATLAB程序10.8.3  求解常微分方程（组）边值问题数值解的MATLAB库函数

模糊自适应PID控制是在PID算法的基础上，以误差e和误差变化率ec作为输入，利用模糊规则进行模糊推理，查询模糊矩阵表进行参数调整，来满足不同时刻的e和ec对PID参数自整定的要求。 【源码目录】 自适应模糊PID ├── Fuzzypid2.mdl ├── __Previews │   └── work222.mdlPreview ├── fight.mdl ├── fuzzyPID.fis ├── mohupid1.fis └── work222.mdl 1 directory, 6 files

【实例简介】

权重确定方法归纳多指标综合评价是指人们根据不同的评价目的，选择相应的评价形式 据此选择多个因素或指标，并通过一定的评价方法将多个评价因素或指标转化为能反映评价对象总体特征的信息，其中评价指标与权重系数确定将直接影响综合评价的结果。按照权数产生方法的不同多指标综合评价方法可分为主观赋权评价法和客观赋权评价法两大类，其中主观赋权评价法采取定性的方法由专家根据经验进行主观判断而得到权数，然后再对指标进行综合评价，如层次分析法、综合评分法、模糊评价法、指数加权法和功效系数法等。客观赋权评价法则根据指标之间的相关关系或各项指标的变异系数来确定权数进行综合评价，如熵值法、神经网络分析法、TOPSIS法、灰色关联分析法、主成分分析法、变异系数法等。两种赋权方法特点不同，其中主观赋权评价法依据专家经验衡量各指标的相对重要性，有一定的主观随意性，受人为因素的干扰较大，在评价指标较多时难以得到准确的评价。客观赋权评价法综合考虑各指标间的相互关系，根据各指标所提供的初始信息量来确定权数，能够达到评价结果的精确 但是当指标较多时，计算量非常大。下面就对当前应用较多的评价方法进行阐述。一、变异系数法（一）变异系数法简介变异系数法是直接利用各项指标所包含的信息，通过计算得到指标的权重。是一种客观赋权的方法。此方法的基本做法是：在评价指标体系中，指标取值差异越大的指标，也就是越难以实现的指标，这样的指标更能反映被评价单位的差距。例如，在评价各个国家的经济发展状况时，选择人均国民生产总值(人均GNP)作为评价的标准指标之一，是因为人均GNP不仅能反映各个国家的经济发展水平，还能反映一个国家的现代化程度。如果各个国家的人均GNP没有多大的差别，则这个指标用来衡量现代化程度、经济发展水平就失去了意义。由于评价指标体系中的各项指标的量纲不同，不宜直接比较其差别程度。为了消除各项评价指标的量纲不同的影响，需要用各项指标的变异系数来衡量各项指标取值的差异程度。各项指标的变异系数公式如下： 式中：是第项指标的变异系数、也称为标准差系数；是第项指标的标准差；是第项指标的平均数。各项指标的权重为： （二）案例说明例如，英国社会学家英克尔斯提出了在综合评价一个国家或地区的现代化程度时，其各项指标的权重的确定方法就是采用的变异系数法。案例：利用变异系数法综合评价一个国家现代化程度时的指标体系中的各项指标的权重。数据资料是选取某一年的数据，包括中国在内的中等收入水平以上的近40个国家的10项指标作为评价现代化程度的指标体系，计算这些国家的变异系数，反映出各个国家在这些指标上的差距，并作为确定各项指标权重的依据。其标准差、平均数数据及其计算出的变异系数等见表1-1。   表1-1      现代化水平评价指标的权重指标人均GNP农业占GDP的比重第三产业占GDP比重非农业劳动力比重城市人口比重人口自然增长率平均预期寿命成人识字率大学生占适龄人口比重每千人拥有医生总(美元)(%)(%)(%)(%)(%)(岁)(%)（%）(人)和平均数11938.49.35254.860.82669.7920.721472.63293.3436.5562.446—标准差7966.277.31612.940.1719.3390.83195.3759.0520.4771.314—变异系数0.6670.7820.2360.2060.2771.1530.0740.0970.560.5374.59权重0.1450.170.0510.0450.060.2510.0160.0210.1220.1171 计算过程如下：（1）先根据各个国家的指标数据，分别计算这些国家每个指标的平均数和标准差；（2）根据均值和标准差计算变异系数。即：这些国家人均GNP的变异系数为： 农业占GDP比重的变异系数： 其他类推。（3）将各项指标的变异系数加总： （4）计算构成评价指标体系的这10个指标的权重：人均GNP的权重： 农业占GDP比重的权重： 其他指标的权重都以此类推。（三）变异系数法的优点和缺点当由于评价指标对于评价目标而言比较模糊时，采用变异系数法评价进行评定是比较合适的，适用各个构成要素内部指标权数的确定，在很多实证研究中也多数采用这一方法。缺点在于对指标的具体经济意义重视不够，也会存在一定的误差。二、层次分析法（一）层次分析法概述人们在对社会、经济以及管理领域的问题进行系统分析时，面临的经常是一个由相互关联、相互制约的众多因素构成的复杂系统。层次分析法则为研究这类复杂的系统，提供了一种新的、简洁的、实用的决策方法。层次分析法(AHP法) 是一种解决多目标的复杂问题的定性与定量相结合的决策分析方法。该方法将定量分析与定性分析结合起来，用决策者的经验判断各衡量目标能否实现的标准之间的相对重要程度，并合理地给出每个决策方案的每个标准的权数，利用权数求出各方案的优劣次序，比较有效地应用于那些难以用定量方法解决的课题。（二）层次分析法原理层次分析法根据问题的性质和要达到的总目标，将问题分解为不同的组成因素，并按照因素间的相互关联影响以及隶属关系将因素按不同层次聚集组合，形成一个多层次的分析结构模型，从而最终使问题归结为最低层(供决策的方案、措施等)相对于最高层(总目标)的相对重要权值的确定或相对优劣次序的排定。层次分析法的特点是在对复杂的决策问题的本质、影响因素及其内在关系等进行深入分析的基础上，利用较少的定量信息使决策的思维过程数学化，从而为多目标、多准则或无结构特性的复杂决策问题提供简便的决策方法。尤其适合于对决策结果难于直接准确计量的场合。（三）层次分析法的步骤和方法 1. 建立层次结构模型利用层次分析法研究问题时，首先要把与问题有关的各种因素层次化，然后构造出一个树状结构的层次结构模型，称为层次结构图。一般问题的层次结构图分为三层，如图所示。最高层为目标层（O）：问题决策的目标或理想结果，只有一个元素。中间层为准则层（C）：包括为实现目标所涉及的中间环节各因素，每一因素为一准则，当准则多于9个时可分为若干个子层。最低层为方案层（P）：方案层是为实现目标而供选择的各种措施，即为决策方案。一般说来，各层次之间的各因素，有的相关联，有的不一定相关联；各层次的因素个数也未必一定相同．实际中，主要是根据问题的性质和各相关因素的类别来确定。         层次分析法所要解决的问题是关于最低层对最高层的相对权重问题，按此相对权重可以对最低层中的各种方案、措施进行排序，从而在不同的方案中作出选择或形成选择方案的原则。2. 构造判断（成对比较）矩阵构造比较矩阵主要是通过比较同一层次上的各因素对上一层相关因素的影响作用．而不是把所有因素放在一起比较，即将同一层的各因素进行两两对比。比较时采用相对尺度标准度量，尽可能地避免不同性质的因素之间相互比较的困难。同时，要尽量依据实际问题具体情况，减少由于决策人主观因素对结果造成的影响。设要比较个因素对上一层（如目标层）的影响程度，即要确定它在中所占的比重。对任意两个因素和，用表示和对的影响程度之比，按1～9的比例标度来度量．于是，可得到两两成对比较矩阵，又称为判断矩阵，显然,因此，又称判断矩阵为正互反矩阵．比例标度的确定：取1-9的9个等级，取的倒数,1-9标度确定如下：= 1，元素与元素对上一层次因素的重要性相同；= 3，元素比元素略重要；= 5，元素比元素重要；= 7， 元素比元素重要得多；= 9，元素比元素的极其重要；，元素与的重要性介于与之间；，当且仅当。由正互反矩阵的性质可知，只要确定的上（或下）三角的个元素即可。在特殊情况下，如果判断矩阵的元素具有传递性，即满足则称为一致性矩阵，简称为一致阵．3. 层次单排序及一致性检验3.1相对权重向量确定（1）和积法取判断矩阵个列向量归一化后的算术平均值，近似作为权重，即类似地，也可以对按行求和所得向量作归一化，得到相应的权重向量。（2）求根法（几何平均法）将的各列（或行）向量求几何平均后归一化，可以近似作为权重，即（3）特征根法设想把一大石头分成个小块，其重量分别为，则将块小石头作两两比较，记的相对重量为，于是可得到比较矩阵显然，为一致性正互反矩阵，记，即为权重向量．且则这表明为矩阵的特征向量，且为特征根．事实上：对于一般的判断矩阵有，这里是的最大特征根，为对应的特征向量．将作归一化后可近似地作为的权重向量，这种方法称为特征根法。注：现有软件求得最大特征根与特征向量。3.2一致性检验通常情况下，由实际得到的判断矩阵不一定是一致的，即不一定满足传递性和一致性．实际中，也不必要求一致性绝对成立，但要求大体上是一致的，即不一致的程度应在容许的范围内．主要考查以下指标：（1）一致性指标：．（2）随机一致性指标：，通常由实际经验给定的，如表2-1。表2-1  随机一致性指标 （3）一致性比率指标：，当时，认为判断矩阵的一致性是可以接受的，则对应的特征向量可以作为排序的权重向量。此时其中表示的第个分量。4.计算组合权重和组合一致性检验（1）组合权重向量设第层上个元素对总目标（最高层）的排序权重向量为第层上个元素对上一层(层)上第个元素的权重向量为则矩阵是阶矩阵，表示第层上的元素对第层各元素的排序权向量．那么第层上的元素对目标层（最高层）总排序权重向量为或    对任意的有一般公式其中是第二层上各元素对目标层的总排序向量．（2）组合一致性指标设层的一致性指标为，随机一致性指标为则第层对目标层的（最高层）的组合一致性指标为组合随机一致性指标为组合一致性比率指标为当时，则认为整个层次的比较判断矩阵通过一致性检验．（四）案例说明实例：人们在日常生活中经常会碰到多目标决策问题，例如假期某人想要出去旅游，现有三个目的地（方案）：风光绮丽的杭州（ ）、迷人的北戴河（）和山水甲天下的桂林（）。假如选择的标准和依据（行动方案准则）有5个景色，费用，饮食，居住和旅途。1.建立层次结构模型 目标层　  准则层   2.构造判断矩阵      构造所有相对于不同准则的方案层判断矩阵（1）相对于景色   （2）相对于费用    （3）相对于居住    （4）相对于饮食      （5）相对于旅途   3. 层次单排序及一致性检验3.1用matlab求得判断矩阵的最大特征根与特征向量：，对应于的正规化的特征向量为：判断矩阵的最大特征值与特征向量判断矩阵的最大特征值与特征向量判断矩阵的最大特征值与特征向量判断矩阵的最大特征值与特征向量判断矩阵的最大特征值与特征向量4.一致性检验对于判断矩阵进行一致性检验：查表知平均随机一致性指标RI，从而可检验矩阵一致性：同理，对于第二层次的景色、费用、居住、饮食、旅途五个判断矩阵的一致性检验均通过。利用层次结构图绘出从目标层到方案层的计算结果： 5.层次总排序各个方案优先程度的排序向量为：决策结果是首选旅游地为 其次为，最后为。（五）优点与缺点人们在进行社会的、经济的以及科学管理领域问题的系统分析中，面临的常常是一个由相互关联、相互制约的众多因素构成的复杂而往往缺少定量数据的系统。层次分析法为这类问题的决策和排序提供了一种新的、简洁而实用的建模方法。在应用层次分析法研究问题时，遇到的主要困难有两个：（i）如何根据实际情况抽象出较为贴切的层次结构；（ii）如何将某些定性的量作比较接近实际定量化处理。层次分析法对人们的思维过程进行了加工整理，提出了一套系统分析问题的方法，为科学管理和决策提供了较有说服力的依据。但层次分析法也有其局限性，主要表现在：（i）它在很大程度上依赖于人们的经验，主观因素的影响很大，它至多只能排除思维过程中的严重非一致性，却无法排除决策者个人可能存在的严重片面性。（ii）当指标量过多时，对于数据的统计量过大，此时的权重难以确定。AHP 至多只能算是一种半定量（或定性与定量结合）的方法。三、熵值法（一）熵值法的原理在信息论中，熵是对不确定性的一种度量。信息量越大，不确定性就越小，熵也就越小；信息量越小，不确定性越大，熵也越大。根据熵的特性，我们可以通过计算熵值来判断一个事件的随机性及无序程度，也可以用熵值来判断某个指标的离散程度，指标的离散程度越大，该指标对综合评价的影响越大。（二）算法实现过程1.数据矩阵   其中为第个方案第个指标的数值。2. 数据的非负数化处理  由于熵值法计算采用的是各个方案某一指标占同一指标值总和的比值，因此不存在量纲的影响，不需要进行标准化处理，若数据中有负数，就需要对数据进行非负化处理。此外，为了避免求熵值时对数的无意义，需要进行数据平移：对于越大越好的指标：   对于越小越好的指标：为了方便起见，仍记非负化处理后的数据为3.计算第项指标下第个方案占该指标的比重4.计算第项指标的熵值   5.计算第项指标的差异系数。对于第项指标，指标值的差异越大，对方案评价的作用越大，熵值就越小。6.求权数7.计算各方案的综合得分（三）熵值法的优缺点熵值法是根据各项指标值的变异程度来确定指标权数的，这是一种客观赋权法，避免了人为因素带来的偏差，但由于忽略了指标本身重要程度，有时确定的指标权数会与预期的结果相差甚远，同时熵值法不能减少评价指标的维数。四、主成分分析法（一）主成分分析法简介  主成分分析是将多个变量通过线性变换以选出较少个数重要变量的一种多元统计分析方法，又称主分量分析。在实际问题中，为了全面分析问题，往往提出很多与此有关的变量（或因素），因为每个变量都在不同程度上反映这个课题的某些信息。但是，在用统计分析方法研究这个多变量的课题时，变量个数太多就会增加课题的复杂性。人们自然希望变量个数较少而得到的信息较多。在很多情形，变量之间是有一定的相关关系的，当两个变量之间有一定相关关系时，可以解释为这两个变量反映此课题的信息有一定的重叠。  主成分分析是对于原先提出的所有变量，建立尽可能少的新变量，使得这些新变量是两两不相关的，而且这些新变量在反映问题的信息方面尽可能保持原有的信息。信息的大小通常用离差平方和或方差来衡量。（二）主成分分析原理主成分分析是数学上对数据降维的一种方法。其基本思想是设法将原来众多的具有一定相关性的指标，，…，（比如个指标），重新组合成一组较少个数的互不相关的综合指标来代替原来指标。那么综合指标应该如何去提取，使其既能最大程度的反映原变量所代表的信息，又能保证新指标之间保持相互无关（信息不重叠）。设表示原变量的第一个线性组合所形成的主成分指标，即,由数学知识可知，每一个主成分所提取的信息量可用其方差来度量，其方差越大，表示包含的信息越多。常常希望第一主成分 所含的信息量最大，因此在所有的线性组合中选取的应该是，，…，的所有线性组合中方差最大的，故称为第一主成分。如果第一主成分不足以代表原来个指标的信息，再考虑选取第二个主成分指标，为有效地反映原信息，已有的信息就不需要再出现在中，即与要保持独立、不相关，用数学语言表达就是其协方差，所以是与不相关的，，…，的所有线性组合中方差最大的，故称为第二主成分，依此类推构造出的、、……、为原变量指标，，…，第一、第二、……、第m个主成分。（三）利用主成分确定权重现举例说明： 假设我们对反映某卖场表现的4项指标（实体店、信誉、企业形象、服务）进行消费者满意度调研。调研采取4级量表，分值越大，满意度越高。现回收有效问卷2000份，并用SPSS录入了问卷数据。部分数据见下图。图4-1 主成分确定权重示例数据（部分）  1、操作步骤： Step1：选择菜单：分析——降维——因子分析 Step2：将4项评价指标选入到变量框中 Step3：设置选项，具体设置如下：  2、 输出结果分析 按照以上操作步骤，得到的主要输出结果为表1——表3，具体结果与分析如下：表4-1  KMO 和 Bartlett 的检验      表3是对本例是否适合于主成分分析的检验。KMO的检验标准见图2。图4-2   KMO检验标准   从图3可知，本例适合主成分分析的程度为‘一般’，基本可以用主成分分析求权重。表4-2  解释的总方差      从表4可知，前2个主成分对应的特征根>1，提取前2个主成分的累计方差贡献率达到94.513% ，超过80%。因此前2个主成分基本可以反映全部指标的信息，可以代替原来的4个指标（实体店、信誉、企业形象、服务）。表4-3   成份矩阵      从表3可知第一主成分与第二主成分对原来指标的载荷数。例如，第一主成分对实体店的载荷数为0.957。 3、确定权重    指标权重等于以主成分的方差贡献率为权重，对该指标在各主成分线性组合中的系数的加权平均的归一化， 因此，确定指标权重需要知道三点：①指标在各主成分线性组合中的系数②主成分的方差贡献率③指标权重的归一化。 （1）指标在不同主成分线性组合中的系数 用表4-3中的载荷数除以表4-2中第1列对应的特征根的开方。 例如，在第一主成分 的线性组合中，实体店的=0.957/(2.775)1/2 ≈0.574。 按此方法，基于表4-3和表4-2的数据，在excel中可分别计算出各指标在两个主成分线性组合中的系数（见下图）。图4-3 各指标在两个主成分线性组合中的系数     由此得到的两个主成分线性组合如下： （2）主成分的方差贡献率 表4中“初始特征值”的“方差%”表示各主成分方差贡献率，方差贡献率越大则该主成分的重要性越强。  因此，方差贡献率可以看成是不同主成分的权重。 由于原有指标基本可以用前两个主成分代替，因此，指标系数可以看成是以这两个主成分方差贡献率为权重，对指标在这两个主成分线性组合中的系数做加权平均。 按上述思路，实体店这个指标的系数为：  这样，我们可以用excel计算出所有指标的系数（见下图）图4-4 所有指标在综合得分模型中的系数 由此得到综合得分模型为：（3）指标权重的归一化    由于所有指标的权重之和为1，因此指标权重需要在综合模型中指标系数的基础上归一化 。    图4-5 指标权重的确定 上图显示了我们基于主成分分析，最终所得到的指标权重。（四）优点和缺点主成分分析法具有处理多个具有一定相关性变量的能力，因此，主成分分析法使用与任何领域的多变量分析。主成分分析法对于各评价指标排序，可以直观的分析出起决定性作用和对综合评价结果影响较大的评价指标。但是，它对于主要指标的依赖性过大，对研究所选取指标体系是一个考验。