矩阵特征向量 权重.docx

1、矩阵特征向量 权重矩阵特征向量 权重 项目六 矩阵的特征值与特征向量 实验1 求矩阵的特征值与特征向量 实验目的 学习利用Mathematica(4.0以上版本) 命令求方阵的特征值和特征向量; 能利用软件计算方 阵的特征值和特征向量及求二次型的标准形. 基本命令 1. 求方阵M 的特征值的命令EigenvaluesM 2. 求方阵M 的特征向量的命令EigenvectorsM 3. 求方阵M 的特征值和特征向量的命令EigensystemM 注:在使用后面两个命令时, 如果输出中含有零向量, 则输出中的非零向量才是真正的特 征向量. 4. 对向量组施行正交单位化的命令GramSchmidt

2、使用这个命令, 先要调用“线性代数. 向量组正交化”软件包, 输入 执行后, 才能对向量组施行正交单位化的命令. 命令GramSchmidtA给出与矩阵A 的行向量组等价的且已正交化的单位向量组. 5. 求方阵A 的相似变换矩阵S 和相似变换的约当标准型J 的命令 JordanDecompositionA 注:因为实对称阵的相似变换的标准型必是对角阵. 所以, 如果A 为实对称阵, 则 JordanDecompositionA同时给出A 的相似变换矩阵S 和A 的相似对角矩阵. 实验举例 求方阵的特征值与特征向量. -102 例1. 1 求矩阵A = 12-1. 的特征值与特值向量. 130

3、1. 求矩阵A 的特征值. 输入 A=-1,0,2,1,2,-1,1,3,0 MatrixFormA EigenvaluesA 则输出A 的特征值 -1,1,1 2. 求矩阵A 的特征向量. 输入 A=-1,0,2,1,2,-1,1,3,0 MatrixFormA EigenvectorsA 则输出 -3,1,0,1,0,1,0,0,0 -31 即A 的特征向量为 1, 0. 0 1 3. 利用命令Eigensystem 同时矩阵A 的所有特征值与特征向量. 输入 A=-1,0,2,1,2,-1,1,3,0 MatrixFormA EigensystemA 则输出矩阵A 的特征值及其对应的特征

4、向量. 1/31/3-1/2 例1. 2 求矩阵A = 1/51-1/3的特征值和特征向量的近似值. 61-2 输入 A=1/3,1/3,-1/2,1/5,1,-1/3,6,1,-2; EigensystemA 则屏幕输出的结果很复杂, 原因是矩阵A 的特征值中有复数且其精确解太复杂. 此时, 可采用 近似形式输入矩阵A , 则输出结果也采用近似形式来表达. 输入 A=1/3,1/3,-1/2,1/5,1,-1/3,6.0,1,-2; EigensystemA 则输出 -0.748989+1.27186i,-0.748989-1.27186i,0.831311, 0.179905+0.1921

5、68i,0.116133+0.062477I,0.955675+0.i, 0.179905-0.192168i,0.116133-0.062477i,0.955675+0.i, -0.0872248,-0.866789,-0.490987 从中可以看到A 有两个复特征值与一个实特征值. 属于复特征值的特征向量也是复的; 属于实 特征值的特征向量是实的. 300 例1. 3 已知2是方阵A = 1t 3的特征值, 求t . 123 输入 ClearA,q; A=2-3,0,0,-1,2-t,-3,-1,-2,2-3; q=DetA Solveq=0,t 则输出 t8 即当t =8时,2是方阵A

6、的特征值. 2-12 例1. 4 已知x =(1, 1, -1) 是方阵A = 5a 3的一个特征向量, 求参数a , b 及特征向 -1b -2 量x 所属的特征值. 设所求特征值为t , 输入 ClearA,B,v,a,b,t; A=t-2,1,-2,-5,t-a,-3,1,-b,t+2; v=1,1,-1; B=A.v; SolveB1=0,B2=0,B3=0,a,b,t 则输出 a-3, b0, t-1 即a =-3, b =0时, 向量x =(1, 1, -1) 是方阵A 的属于特征值-1和特征向量. 矩阵的相似变换 411 例1. 5 设矩阵A = 222, 求一可逆矩阵P , 使

7、P -1AP 为对角矩阵. 222 方法1 输入 ClearA,P; A=4,1,1,2,2,2,2,2,2; EigenvaluesA P=EigenvectorsA/Transpose 则输出 0,2,6 0,-1,1,-1,1,1,1,1,1 0 即矩阵A 的特征值为0,2,6. 特征向量为 -1 1 -1可验证P AP 为对角阵, 事实上, 输入 -110-11 , 1与 1, 矩阵P = -111. 1 1 111 InverseP.A.P 则输出 0,0,0,0,2,0,0,0,6 因此, 矩阵A 在相似变换矩阵P 的作用下, 可化作对角阵. 方法2 直接使用JordanDecom

8、position 命令, 输入 jor=JordanDecompositionA 则输出 0,-1,1,-1,1,1,1,1,1,0,0,0,0,2,0,0,0,6 可取出第一个矩阵S 和第二个矩阵, 事实上, 输入 jor1 jor2 则输出 0,-1,1,-1,1,1,1,1,1 0,0,0,0,2,0,0,0,6 输出结果与方法1的得到的结果完全相同. -200-100 例1. 6 已知方阵A = 2x 2与B = 020相似, 求x , y . 311 00y 注意矩阵B 是对角矩阵, 特征值是-1, 2, y . 又矩阵A 是分块下三角矩阵,-2是矩阵A 的特 征值. 矩阵A 与B

9、相似, 则y =-2, 且-1,2也是矩阵A 的特征值. 输入 Clearc,v; v=4,0,0,-2,2-x,-2,-3,-1,1; SolveDetv=0,x 则输出 x0 所以, 在题设条件, x =0, y =-2. 例1. 7 已知二次型 222f (x 1, x 2, x 3) =x 1-2x 2+x 3+2x 1x 2-4x 1x 3+2x 2x 3 (1)求标准形; (2)求正惯性指数; (3)判断二次型是否正定. 输入 A=1,1,-2,1,-2,1,-2,1,1 EigenvaluesA 则输出矩阵A 的特征值为 -3,0,3 22+3y 2所以二次型的标准形为f =3y

10、 1; 正惯性指数为1; 该二次型不是正定的. 例1. 8 求正交变换将二次型 2222f (x 1, x 2, x 3) =x 1+x 2+x 3+x 4+2x 1x 2-2x 1x 4+2x 2x 3-2x 3x 4 化为标准形. 输入 A=1,1,0,-1,1,1,1,0,0,1,1,-1,-1,0,-1,1 MatrixFormA X=x1,x2,x3,x4; ExpandX.A.X P=GramSchmidtEigenvectorsA P.A.InverseP/MatrixForm 则输出所求的正交变换矩阵P 与二次型矩阵A 标准形. 从结果知, 所求二次型的标准型为 2222g =

11、-y 1+y 2+y 3+y 4 实验习题 2-12 1. 求方阵A = 2-1-2的特征值与特征向量. 2-2-1 1111 11-1-12. 求方阵A = 的特征值与特征向量. 1-11-1 1-1-11 101 3. 已知:0是方阵 020的特征值, 求. 10t t 211 4. 设向量x =(1, k , 1) T 是方阵A = 121的特征向量, 求k . 112 0-12 5. 方阵A = 010是否与对角阵相似? 1-11 200200 6. 已知:方阵A = 001与B = 0y 0相似, 01x 00-1 (1)求x 与y ; (2)求一个满足关系P -1AP =B 的方阵

12、P . 124 7. 设方阵A = 2-22, 求正交阵C , 使得B =C T AC 是对角阵. 421 实验2 层次分析法 实验目的 通过应用层次分析法解决一个实际问题, 学习层次分析法的基本原理与方法; 掌握用层次 分析法建立数学模型的基本步骤; 学会用Mathematica 解决层次分析法中的数学问题. 基本原理 层次分析法是系统分析的重要工具之一, 其基本思想是把问题层次化、数量化, 并用数学 方法为分析、决策、预报或控制提供定量依据. 它特别适用于难以完全量化, 又相互关联、 相互制约的众多因素构成的复杂问题. 它把人的思维过程层次化、数量化, 是系统分析的一中 新型的数学方法.

13、运用层次分析法建立数学模型, 一般可按如下四个基本步骤进行. 1. 建立层次结构 首先对所面临的问题要掌握足够的信息, 搞清楚问题的范围、因素、各因素之间的相互 关系, 及所要解决问题的目标. 把问题条理化、层次化, 构造出一个有层次的结构模型. 在这 个模型下, 复杂问题被分解为元素的组成部分. 这些元素又按其属性及关系形成若干层次. 层 次结构一般分三层: 第一层为最高层, 它是分析问题的预定目标和结果, 也称目标层; 第二层为中间层, 它是为了实现目标所涉及的中间环节, 如: 准则、子准则, 也称准则 层; 第三层为最底层, 它包括了为实现目标可供选择的各种措施、决策方案等, 也称方案层

14、. （章栋恩P268图26.1） 图2-1 注:上述层次结构具有以下特点:(1) 从上到下顺序地存在支配关系, 并用直线段表示;(2) 整个层次结构中层次数不受限制. 2. 构造判断矩阵 构造判断矩阵是建立层次分析模型的关键. 假定以上一层的某元素y 为准则, 它所支配 的下一层次的元素为x 1, x 2, , x n , 这n 个元素对上一层次的元素y 有影响, 要确定它们在y 中的比重. 采用成对比较法. 即每次取两个元素x i 和x j , 用a ij 表示x i 与x j 对y 的影响之比, 全部比较的结果可用矩阵A 表示, 即 A =(a ij ) n n , i , j =1, 2

15、, , n . 称矩阵A 为判断矩阵. 根据上述定义, 易见判断矩阵的元素a ij 满足下列性质: a ji =1(i j ), a ij a ii =1, (i =j ) 当a ij 0时, 我们称判断矩阵A 为正互反矩阵. 怎样确定判断矩阵A 的元素a ij 的取值呢? 当某层的元素x 1, x 2, , x n 对于上一层某元素y 的影响可直接定量表示时, x i 与x j 对y 的影响之比可以直接确定, a ij 的值也可直接确定. 但对于大多数社会经济问题, 特别是比较 复杂的问题, 元素x i 与x j 对y 的重要性不容易直接获得, 需要通过适当的量化方法来解决. 通常取数字19

16、及其倒数作为a ij 的取值范围. 这是因为在进行定性的成对比较时, 通常采用 5级制(表1), 在每两个等级之间各有一个中间状态, 共19个尺度, 另外心理学家认为进行成 对比较的因素太多, 将超出人们的判断比较能力, 降低精确. 实践证明, 成对比较的尺度以 72为宜, 故a ij 的取值范围是1, 2, , 9及其倒数. 表1 比较尺度a ij 的取值 x i /x j a ij 相等较强强很强绝对强 13579 3. 计算层次单排序并做一致性检验 层次单排序是指同一层次各个元素对于上一层次中的某个元素的相对重要性进行排序. 具体做法是: 根据同一层n 个元素x 1, x 2, , x

17、n 对上一层某元素y 的判断矩阵A , 求出它们对 于元素y 的相对排序权重, 记为w 1, w 2, , w n , 写成向量形式w =(w 1, w 2, , w n ) T , 称其为A 的层次单排序权重向量, 其中w i 表示第i 个元素对上一层中某元素y 所占的比重, 从而得到 层次单排序. 层次单排序权重向量有几种求解方法, 常用的方法是利用判断矩阵A 的特征值与特征向 量来计算排序权重向量w . 关于正互反矩阵A , 我们不加证明地给出下列结果. (1) 如果一个正互反矩阵A =(a ij ) n n 满足 a ij a jk =a ik (i , j , k =1, 2, ,

18、n ) 则称矩阵A 具有一致性, 称元素x i , x j , x k 的成对比较是一致的; 并且称A 为一致矩阵. (2) n 阶正互反矩阵A 的最大特征根max n , 当=n 时, A 是一致的. (3) n 阶正互反矩阵是一致矩阵的充分必要条件是最大特征值 max =n . 计算排序权重向量的方法和步骤 设w =(1, 2, , n ) T 是n 阶判断矩阵的排序权重向量, 当A 为一致矩阵时, 根据n 阶判断矩阵构成的定义, 有 111 n 12 222 A = 12n (2.1) n n n n 12 因而满足Aw =nw , 这里n 是矩阵A 的最大特征根, w 是相应的特征向量

19、; 当A 为一般的 判断矩阵时Aw =max w , 其中max 是A 的最大特征值(也称主特征根), w 是相应的特征向 量(也称主特征向量). 经归一化(即 i =1n i =1) 后, 可近似作为排序权重向量, 这种方法称为 特征根法. 一致性检验 在构造判断矩阵时, 我们并没有要求判断矩阵具有一致性, 这是由客观事物的复杂性 与人的认识的多样性所决定的. 特别是在规模大、因素多的情况下, 对于判断矩阵的每个元 素来说, 不可能求出精确的i /j , 但要求判断矩阵大体上应该是一致的. 一个经不起推敲 的判断矩阵有可能导致决策的失误. 利用上述方法计算排序权重向量, 当判断矩阵过于偏离

20、一致性时, 其可靠性也有问题. 因此, 需要对判断矩阵的一致性进行检验, 检验可按如下步骤 进行: (1) 计算一致性指标CI -n CI =max (2.2) n -1 当CI =0, 即max =n 时, 判断矩阵A 是一致的. 当CI 的值越大, 判断矩阵A 的不一致的程 度就越严重. (2) 查找相应的平均随机一致性指标RI 表2给出了n (111) 阶正互反矩阵的平均随机一致性指标RI , 其中数据采用了 100150个随机样本矩阵A 计算得到. (3) 计算一致性比例CR CR =CI (2.3) RI 当C R 4. 计算层次总排序权重并做一致性检验 计算出某层元素对其上一层中某

21、元素的排序权重向量后, 还需要得到各层元素, 特别 是最底层中各方案对于目标层的排序权重, 即层次总排序权重向量, 再进行方案选择. 层次 总排序权重通过自上而下地将层次单排序的权重进行合成而得到. 考虑3个层次的决策问题: 第一层只有1个元素, 第二层有n 个元素, 第三层有m 个元 素. 设第二层对第一层的层次单排序的权重向量为 (2) (2) (2) T , 2, , n ) w (2) =(1(3) (3) (3) (3) T =(w k w k 1, w k 2, , w kn ) , k =1, 2, , n 第三层对第二层的层次单排序的权重向量为 (3) 以w k 为列向量构成矩

22、阵: (3) (3) (3) w 11, w 21, , w n 1, (3) (3) w (3) , w 22, , w n (3) (3) (3) 2, (2.4) W (3) =(w 1, w 2, , w n ) = 12 w (3) , w (3) , , w (3) nm m n 2m 1m 则第三层对第一层的层次总排序权重向量为 w (3) =W (3) w (2) (2.5) 一般地, 若层次模型共有s 层, 则第k 层对第一层的总排序权重向量为 w (k ) =W (k ) w (k -1) , k =3, 4, , s (2.6) 其中W (k ) 是以第k 层对第k -1

23、层的排序权向量为列向量组成的矩阵, w (k -1) 是第k -1层对第 一层的总排序权重向量. 按照上述递推公式, 可得到最下层(第s 层) 对第一层的总排序权重 向量为 w (s ) =W (s ) W (s -1) W (3) w (2) (2.7) 对层次总排序权重向量也要进行一致性检验. 具体方法是从最高层到最低层逐层进行 检验. 如果所考虑的层次分析模型共有s 层. 设第l (3l s ) 层的一致性指标与随机一致性 (l ) (l ) (l ) (l ) (l ) (l ) 指标分别为CI 1(n 是第l -1层元素的数目) 与RI 1, 令 , CI 2, , CI n , R

24、I 2, , RI n (l ) (l ) CI (l ) =CI 1, , CI 1w (l -1) (l ) (l ) RI (l ) =RI 1, , RI 1w (l -1) (2.8) (2.9) 则第l 层对第一层的总排序权向量的一致性比率为 CI (l ) (2.10) , l =3, 4, , s RI (l ) 其中CR (2) 为由(2.3)式计算的第二层对第一层的排序权重向量的一致性比率. CR (l ) =CR (l -1) + 当最下层对第一层的总排序权重向量的一致性比率CR (s ) 应用举例 问题 在选购电脑时, 人们希望花最少的钱买到最理想的电脑. 试通过层次分

25、析法建立 数学模型, 并以此确定欲选购的电脑. 1. 建立选购电脑的层次结构模型 （章栋恩P268图26.2 左边加目标层、准则层、方案层字样） 图2-2 该层次结构模型共有三层:目标层(用符号z 表示最终的选择目标); 准则层(分别用符号 y 1, y 2, , y 5表示“性能”、“价格”、“质量”、“外观”、“售后服务”五个判断准则); 方案层(分别用符号x 1, x 2, x 3表示品牌1, 品牌2, 品牌3三种选择方案). 2. 构造成对比较判断矩阵 (1) 建立准则层对目标层的成对比较判断矩阵 根据表1的定量化尺度, 从建模者的个人观点出发, 设准则层对目标层的成对比较判断矩阵为

26、539311/511/221/2 (2.11) A =1/32131 1/91/21/311/31/32131 (2) 建立方案层对准则层的成对比较判断矩阵 11/31/51B 1=311/2, B 2=1/3 1521/51531B 4=1/511/2, B 5=1/3 1/3211/33511/31/5 12, B 3=311/2 1/21152 33 11, 11 3. 计算层次单排序权重向量并做一致性检验 先利用Mathematica 计算矩阵A 的最大特征值及特征值所对应的特征向量. 输入 (*调用只求实数运算的软件包*) A=1.0,5,3,9,3,1/5,1,1/2,2,1/2,

27、1/3,2,1,3,1, 1/9,1/2,1/3,1,1/3,1/3,2,1,3,1; (*以小数形式1.0输入进行近似计算, 可避免精确解太长、太复杂*) T=EigensystemA/Chop (*输入/Chop, 把与零非常接近的数换成零*) 则输出 5.00974,Nonreal,Nonreal,0,0, 0.88126,0.167913,0.304926,0.0960557,0.304926, 0.742882,Nonreal,Nonreal,Nonreal,Nonreal, 0.742882,Nonreal,Nonreal,Nonreal,Nonreal, -0.993398,0,

28、0.0673976,0.0662265,0.0650555, -0.65676,0,0.57431,0.043784,-0.486742 (输出中的Nonreal 表示复数) 从中得到A 的最大特征值max =5. 00974, 及其对应的特征向量 x =(0. 88126, 0. 167913, 0. 304926, 0. 0960557, 0. 304926) T 输入 Clearx; x=T2,1; ww2=x/ApplyPlus,x 则得到归一化后的特征向量 w (2) =(0. 502119, 0. 0956728, 0. 173739, 0. 0547301, 0. 173739)

29、 T -n 计算一致性指标CI =max , 其中n =5, max =5. 00974, 故 n -1 C I =0. 002435. 查表得到相应的随机一致性指标 RI =1. 12 从而得到一致性比率 CI CR (2) =0. 002174 RI 因CR (2) B1=B3=1.0,1/3,1/5,3,1,1/2,5,2,1; 11 B2=TransposeB1; B4=1.0,5,3,1/5,1,1/2,1/3,2,1; B5=1.0,3,3,1/3,1,1,1/3,1,1; T1=EigensystemB1/Chop T2=EigensystemB2/Chop T3=EigensystemB3/Chop T4=EigensystemB4/Chop T5=EigensystemB5/Chop 则输出 3.00369,Nonreal, Nonreal,