项目六 矩阵的特征值与特征向量.docVIP

  1. 1、本文档共15页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
项目六 矩阵的特征值与特征向量 实验1 求矩阵的特征值与特征向量 实验目的 学习利用Mathematica(4.0以上版本)命令求方阵的特征值和特征向量;能利用软件计算方 阵的特征值和特征向量及求二次型的标准形. 基本命令 1.求方阵的特征值的命令Eigenvalues[M] 2.求方阵的特征向量的命令Eigenvectors[M] 3.求方阵的特征值和特征向量的命令Eigensystem[M] 注:在使用后面两个命令时,如果输出中含有零向量,则输出中的非零向量才是真正的特 征向量. 4.对向量组施行正交单位化的命令GramSchmidt 使用这个命令,先要调用“线性代数.向量组正交化”软件包,输入 <<LinearAlgebra\Orthogonalization.m 执行后,才能对向量组施行正交单位化的命令. 命令GramSchmidt[A]给出与矩阵的行向量组等价的且已正交化的单位向量组. 5.求方阵的相似变换矩阵和相似变换的约当标准型的命令 JordanDecomposition[A] 注:因为实对称阵的相似变换的标准型必是对角阵. 所以,如果为实对称阵,则 JordanDecomposition[A]同时给出的相似变换矩阵和的相似对角矩阵. 实验举例 求方阵的特征值与特征向量. 例1.1 求矩阵的特征值与特值向量. 1.求矩阵A的特征值. 输入 A={{-1,0,2},{1,2,-1},{1,3,0}} MatrixForm[A] Eigenvalues[A] 则输出A的特征值 {-1,1,1} 2.求矩阵A的特征向量. 输入 A={{-1,0,2},{1,2,-1},{1,3,0}} MatrixForm[A] Eigenvectors[A] 则输出 {{-3,1,0},{1,0,1},{0,0,0}} 即A的特征向量为 3.利用命令Eigensystem同时矩阵A的所有特征值与特征向量. 输入 A={{-1,0,2},{1,2,-1},{1,3,0}} MatrixForm[A] Eigensystem[A] 则输出矩阵A的特征值及其对应的特征向量. 例1.2 求矩阵的特征值和特征向量的近似值. 输入 A={{1/3,1/3,-1/2},{1/5,1,-1/3},{6,1,-2}}; Eigensystem[A] 则屏幕输出的结果很复杂,原因是矩阵的特征值中有复数且其精确解太复杂.此时,可采用 近似形式输入矩阵,则输出结果也采用近似形式来表达. 输入 A={{1/3,1/3,-1/2},{1/5,1,-1/3},{6.0,1,-2}}; Eigensystem[A] 则输出 {{-0.748989+1.27186i,-0.748989-1.27186i,0.831311}, {{0.179905+0.192168i,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}}} 从中可以看到有两个复特征值与一个实特征值.属于复特征值的特征向量也是复的;属于实 特征值的特征向量是实的. 例1.3 已知2是方阵的特征值,求. 输入 Clear[A,q]; A={{2-3,0,0},{-1,2-t,-3},{-1,-2,2-3}}; q=Det[A] Solve[q==0,t] 则输出 {{t8}} 即当时,2是方阵的特征值. 例1.4 已知是方阵的一个特征向量,求参数及特征向 量所属的特征值. 设所求特征值为,输入 Clear[A,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; Solve[{B[[1]]==0,B[[2]]==0,B[[3]]==0},{a,b,t}] 则输出 {{a-3, b0, t-1}} 即时,向量是方阵的属于特征值-1和特征向量. 矩阵的相似变换 例1.5 设矩阵,求一可逆矩阵,使为对角矩阵. 方法1 输入 Clear[A,P]; A={{4,1,1},{2,2,2},{2,2,2}}; Eigenvalues[A] P=Eigenvectors[A]//Transpose 则输出 {0,2,6} {{0,-1,1},{-1,1,1},{1,1,1}} 即矩阵A的特征值为0,2,6.特征向量为,与,矩阵. 可验证为对角阵, 事实上,输入 Inverse[P].A.P 则输出 {{0,0,0},{0,2,0},{0,0,6}} 因此,矩阵在相似变换矩

文档评论(0)

***** + 关注
实名认证
文档贡献者

擅长各种类型文档 找不到的文档可私信

1亿VIP精品文档

相关文档