十四矩阵的特征值与特征向量.doc

《十四矩阵的特征值与特征向量.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、天水师范学院数学与统计学院实验报告实验项目名称矩阵的特征值与特征向量所属课程名称高等数学实验实验类型线性代数实验日期2011.12.14班级08级数应七班学号姓名吴向阳成绩一、实验概述：【实验目的】学习掌握利用Mathematica(4．0以上版本)命令求方阵的特征值和特征向量；利用特征值求二次型的标准形【实验原理】(1)命令Eigenvalues，给出方阵M的特征值．(2)命令Eigenvectors[M]，给出方阵M的特征向量．但是时输出中含有零向量，此时输出中的非零向量才是真正的特征向量．(3)命

2、令eigensystem[M]，给出方阵m的特征值和特征向量．同样，有时输出的向量中含有零向量．(4)调用“线性代数．向量组正交化”软件包命令是<

3、、分析）1．求方阵的特征值与特征向量Clear[M];M={{1,2,3},{2,1,3},{3,3,6}};Eigenvalues[M]Eigenvectors[M]Eigensystem[M]用命令Eigenvalues[M]立即求得方阵M的特征值，命令Eigenvectors[M]立即求得方阵M的特征向量，命令Eigensystem[M]立即求得方阵的特征值和特征向量例14.2求方阵M的特征值和特征向量(*Example14.2*)G={{1/3,1/3,-1/2},{1/5,1,-1/3},{6

4、,1,-2}};Eigensystem[G]G={{1/3,1/3,-1/2},{1/5,1,-1/3},{6.0,1,-2}};Eigensystem[G]例14．，已知：是方阵的特征值，求t．(*Example14.3*)Clear[A,q];A={{2-3,0,0},{-1,2-t,-3},{-1,-2,2-3}};q=Det[A];Solve[q0,t]例14．4已知z=(1，1，一I)是方阵A=J5(*Example14.4*)Clear[A,B,v,a,b,t];A={{t-2,1,-2},

5、{-5,t-a,-3},{1,-b,t+2}};v={1,1,-1};B=A.v;Solve[{B[[1]]0,B[[2]]0,B[[3]]0},{a,b,t}]2-矩阵的相似变换若n阶方阵A有n个线性无关的特征向量，则A与对角阵相似．实对称阵总与对角阵相似，且存在正交阵P，使P’AP为对角阵．命令EigenVectors[A]与Eigensystem[A]给出还未经过正交化和单位化的特征向量因此要对特征向量进行正交化和单位化所用的命令是GramSchmidt[]不过首先要输人调用软件包<

6、Algebra＼Orthogonalization．m的命令例145设方阵A=，求一可逆阵P，使P-1AP为对角阵．Clear[A,p];A={{4,1,1},{2,2,2},{2,2,2}};Eigenvalues[A];p=Eigenvectors[A]//TransposeInverse[p].A.pjor=JordanDecomposition[A]jor[[1]]jor[[2]]例146方阵A是否与对角阵相似?Clear[A];A={{1,0},{2,1}};Eigensystem[A]例14

7、．7Clear[x,v];v={{4,0,0},{-2,2-x,-2},{-3,-1,1}};Solve[Det[v]0,x]例14．8对实对称矩阵A=<

8、标准二次型的矩阵为f=Table[x[j],{j,4}].A.Table[x[j],{j,4}]//Simplifyf/.Table[x[j]®(p.Table[y[j],{j,4}])[[j]],{j,4}]//Simplify【实验结论】（结果）用Mathematica(4．0以上版本)命令求方阵的特征值和特征向量很方便，从而便于利用特征值求二次型的标准形【实验小结】（收获体会）学会了掌握利用Mathematica(4．0以上版本)命令