- 1、本文档共15页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 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}}
因此,矩阵在相似变换矩
您可能关注的文档
- 项目编号:TPCG2011-33.doc
- 项目编号:XZC2015-035.doc
- 项目编号:受理单位填写).doc
- 项目编号:皖C32015GC018.doc
- 项目管理相关文件汇编.doc
- 项目教学法在Flash课程中的应用.doc
- 项目类别:关键共性技术研发.doc
- 项目类型:指南代码:.doc
- 项目名称: 大同利美特(上海)有限公司马鞍山分公司.doc
- 项目名称: 新建机动车辆检测站项目.doc
- 地理-内蒙古锡林郭勒盟2023~2024学年第一学期高三年级全盟统考试题和答案.docx
- 历史|湖北元月调考暨湖北省部分市州元月高三期末联考历史试卷及答案.docx
- 历史-内蒙古锡林郭勒盟2023~2024学年第一学期高三年级全盟统考试题和答案.docx
- 数学(理科)-内蒙古锡林郭勒盟2023~2024学年第一学期高三年级全盟统考试题和答案.docx
- 历史丨金太阳24~239C湖北省十堰市2024届高三上学期1月调研考试历史试卷及答案.docx
- 青岛某高层办公楼地下防水施工方案(sbs卷材).docx
- 土木工程实习报告表格.doc
- 工程流体力学教学课件ppt作者闻建龙工程流体力学习题答案部分.doc
- 围手术期管理制度(1)_0.doc
- 2023年公务员多省联考《申论》题(内蒙古旗县卷)及参考答案.doc
文档评论(0)