No6 矩阵函数33.ppt

  1. 1、本文档共47页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
No6 矩阵函数33

高等代数与矩阵分析;第六章 矩阵函数;§1、矩阵的多项式与极小多项式;于是有;我们称上面的表达式为矩阵多项式 的Jordan表示。其中;;;二、 零化多项式、最小多项式;定义6.3:已知 ,在 的零化多项式中,次数最低且首项系数为1的零化多项式称为 的最小多项式,通常记为 。;例 2 :已知一个Jordan块;因此有;例 2 :求下列矩阵的最小多项式;所以其最小多项式为 。;(3)该矩阵的Jordan标准形为;定义4:设 为 的 个互不相同的特征值, 为其最 小多项式且有;容易求得矩阵 的最小多项式为;容易求得矩阵 的最小多项式为;(1)设 ,如果 有定义,那么 是否也有定义?;§2、矩阵函数及其计算;函数 在矩阵 的谱上有定义,如果存在多项式 且满足;如果函数 在矩阵 的谱上有定义,那么;求 的Jordan表示并计算 ;当 时,可得 从而有;当 时,可得 于是有;求 的Jordan表示并计算;当 时,可得;当 时,可得 故;定理3:设函数 与函数 在矩阵 的影谱上都有定义,那么 的充分必要条件是 与 在 的影谱上的值完全相同。;二、矩阵函数的计算(归纳);Jordan分解法;其中;特征值分解法;例 3 求矩阵函数 、 和 ,其中;当 时 ,则;当 时;%ex601.m A=[-1 1 0 ;-4 3 0; 1 0 2]; expm(A) %调用expm函数 %expm uses the Padé approximation with scaling % and squaring.;%ex601.m(续) A=[-1 1 0 ;-4 3 0; 1 0 2]; syms t % 声明符号变量t S=expm(A*t); S=simple(S) % 简化矩阵函数的结果;当 时;%ex601.m(续) A=[-1 1 0 ;-4 3 0; 1 0 2]; syms t % 声明符号变量t A1=sin(A*t) % 内置函数sin(A)给出错误结果;%ex601.m(续) A=[-1 1 0 ;-4 3 0; 1 0 2]; j=sqrt(-1); %虚数单位i A2=(expm(j*A*t) - expm((-1)*j*A*t))/(2*j) %利用Euler公式,调用函数expm A2=simple(A2) % 简化矩阵函数的结果 ;?? 4 求矩阵函数 和 ,其中; 因此相似矩阵为;%ex602.m A=[4 6 0;-3 -5 0;-3 -6 1]; syms t % 声明符号变量t S=expm(A*t); S=simple(S) % 简化矩阵函数的结果;%ex602.m (续) A=[4 6 0;-3 -5 0;-3 -6 1]; j=sqrt(-1); %虚数单位i B=(expm(j*A)+expm((-1)*j*A))/2 %利用Euler公式,调用函数expm;%ex602.m(续) A=[4 6 0;-3 -5 0;-3 -6 1]; [P,D]=eig(A); j=sqrt(-1); C= P*((expm1(j*D)+expm1((-1)*j*D))/2)*inv(P) +eye(size(A)) %函数expm1返回e(

文档评论(0)

pfenejiarz + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档