- 1、本文档共58页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
矩阵及特征值和特征向量二次型
实验三 一、特征值与特征向量 二、 二次型化标准型 三、 正定二次型的判定 A = [4 –2 2;-2 1 –1/2;2 –1/2 1]; [ P , T ] = schur (A) P = 0.5458 -0.0000 0.8379 0.5925 0.7071 -0.3859 -0.5925 0.7071 0.3859 T = -0.3423 0 0 0 0.5000 0 0 0 5.8423 答:所作的正交变换为: 二次型的标准型为: 1. 顺序主子式判断法 ⑴ 求二次型 F=X’AX 的矩阵 A 的各阶顺序 主子式 Di (i=1,2,3…..); ⑵ 判断 Di 是否大于0 . 程序:建立函数文件 shxu.m function [C,M] =shxu(A) % C为A的各阶顺序主子式组成的向量 % M为判定向量: if C(i)0, then M(i)=1; % others M(i)=0 n=size(A); C=[ ]; M=[ ]; for i=1:n(1) A1=A([1:i],[1:i]); D=det(A1); C=[C D]; if D0 m=1; else m=0; end M=[M,m]; end 2、特征值判别法 ⑴ 求二次型 f =X’AX 的矩阵 A 的全部特征 值 (i=1,2,……); ⑵ 判断 是否大于 0 . 程序:建立函数文件 tezh.m function [ T , M ] = tezh (A) n=size(A); T=(eig(A))’ ; M=[ ]; for i =1:n(1) if T(i)0 m=1; else m=0; end M=[M,m]; end 例8 判定下列二次型是否正定 解 二次型矩阵 方法一 顺序主子式 A = [1 –1 2 1;-1 3 0 –3;2 0 9 –6;1 –3 –6 19] ; [C,M] = shxu (A) 答:此二次型是正定的。 C = 1 2 6 24 M = 1 1 1 1 方法二 特征值法 T = 0.0643 2.2421 7.4945 22.1991 M = 1 1 1 1 A = [ 1 –1 2 1;-1 3 0 –3;2 0 9 –6;1 –3 –6 19] [ T , M ] = tezh (A) 答:此二次型是正定的。 * * 矩阵的特征值和特征向量 二次型 实验目的 1、学会用MATLAB软件求矩阵的特征值和特征向量 2、学会用MATLAB软件将二次型化为标准型 3、通过用MATLAB软件编程来判断二次型的正定性 其中:D为由特征值构成的对角阵,V为由特征向量作为列向量构成的矩阵。且使 AV=VD 成立 用Matlab计算特征值和特征向量的命令如下: d=eig(A) 仅计算A的特征值(以向量形式d存放) [V,D]=eig(A) trace(A) 计算矩阵A的迹 例1:求方阵 的特征值、特征向量和迹 解:
您可能关注的文档
最近下载
- 《无人机艺术摄影》课件——项目5 无人机拍摄视频.ppt VIP
- HD21.2 S3-1997+A1-2002 额定电压450/750V及以下热塑性绝缘电缆 第2部分:试验方法(中文).pdf VIP
- 基于体卫融合的儿童青少年运动指南(2024) .pdf VIP
- GB_T 18920-2020 城市污水再生利用城市杂用水水质.pdf VIP
- 陀螺仪减震装置及无人机.pdf VIP
- 12J003室外工程(标准图集).docx
- 项目管理办公室PMO实践指南.docx VIP
- 必威体育精装版医疗器械分类标准目录.pptx VIP
- 风电场项目风机及附属设备吊装工程(混塔风机吊装施工方案).pptx VIP
- DLT 5265_水电水利工程混凝土搅拌楼安全操作规程.pdf VIP
文档评论(0)