矩阵理论 第八讲 矩阵函数的求法

摘要:
第八课是关于矩阵函数的解。首先,使用Jordan标准形式求解矩阵函数。对于矩阵的多项式,我们推导出:多项式事实上,上述结果不仅对矩阵的多项式有效,而且对矩阵的幂级数也有效。这引出了矩阵函数的另一种定义和计算方法。1.定义:设n阶矩阵A的Jordan正则形式为J,并且存在非奇异矩阵P,使得:对于函数f(z),如果以下函数是有意义的,则矩阵函数f(A)被称为有意义的;2.矩阵函数的计算(步骤):找到Jordan正则形式和A的变化

第八讲 矩阵函数的求法

一、利用Jordan标准形求矩阵函数。

对于矩阵的多项式,我们曾导出矩阵理论 第八讲 矩阵函数的求法第1张矩阵理论 第八讲 矩阵函数的求法第2张:多项式

矩阵理论 第八讲 矩阵函数的求法第3张

矩阵理论 第八讲 矩阵函数的求法第4张

实际上,以上结果不仅对矩阵的多项式成立,对矩阵的幂级数也成立。由此引出矩阵函数的另一种定义及计算方法。

1. 定义:设n阶矩阵A的Jordan标准形为J

矩阵理论 第八讲 矩阵函数的求法第5张矩阵理论 第八讲 矩阵函数的求法第6张

有非奇异矩阵P使得:矩阵理论 第八讲 矩阵函数的求法第7张

对于函数f(z),若下列函数

矩阵理论 第八讲 矩阵函数的求法第8张矩阵理论 第八讲 矩阵函数的求法第9张

均有意义,则称矩阵函数f(A)有意义,且

矩阵理论 第八讲 矩阵函数的求法第10张

矩阵理论 第八讲 矩阵函数的求法第11张2. 矩阵函数的求法(步骤):

矩阵理论 第八讲 矩阵函数的求法第12张    求出A的Jordan标准形及变换矩阵P,矩阵理论 第八讲 矩阵函数的求法第13张

矩阵理论 第八讲 矩阵函数的求法第14张    对于J的各Jordan块矩阵理论 第八讲 矩阵函数的求法第15张求出矩阵理论 第八讲 矩阵函数的求法第16张,即计算出

矩阵理论 第八讲 矩阵函数的求法第17张

并按照顺序构成矩阵理论 第八讲 矩阵函数的求法第18张

矩阵理论 第八讲 矩阵函数的求法第19张

矩阵理论 第八讲 矩阵函数的求法第20张 合成矩阵理论 第八讲 矩阵函数的求法第21张

矩阵理论 第八讲 矩阵函数的求法第22张 矩阵乘积给出矩阵理论 第八讲 矩阵函数的求法第23张

需要说明的是,计算结果与Jordan标准形中Jordan块的顺序无关。

例1 (教材P176例3-8). 矩阵理论 第八讲 矩阵函数的求法第24张,求矩阵理论 第八讲 矩阵函数的求法第25张

[解] 1矩阵理论 第八讲 矩阵函数的求法第26张求出JP

矩阵理论 第八讲 矩阵函数的求法第27张

2矩阵理论 第八讲 矩阵函数的求法第28张 求出矩阵理论 第八讲 矩阵函数的求法第29张并构成矩阵理论 第八讲 矩阵函数的求法第30张

矩阵理论 第八讲 矩阵函数的求法第31张

f(1)=1,矩阵理论 第八讲 矩阵函数的求法第32张

矩阵理论 第八讲 矩阵函数的求法第33张

3矩阵理论 第八讲 矩阵函数的求法第34张 合成矩阵理论 第八讲 矩阵函数的求法第35张

4矩阵理论 第八讲 矩阵函数的求法第36张矩阵理论 第八讲 矩阵函数的求法第37张矩阵理论 第八讲 矩阵函数的求法第38张

说明:

(1)矩阵理论 第八讲 矩阵函数的求法第39张

可见这样的矩阵理论 第八讲 矩阵函数的求法第40张确与矩阵理论 第八讲 矩阵函数的求法第41张构成反函数;

(2)矩阵函数的种类不仅是我们介绍的这种,如辛矩阵。以

矩阵理论 第八讲 矩阵函数的求法第42张为例,以我们这里的定义,矩阵理论 第八讲 矩阵函数的求法第43张,但矩阵理论 第八讲 矩阵函数的求法第44张

亦满足矩阵理论 第八讲 矩阵函数的求法第45张,即B也可以看作某种矩阵理论 第八讲 矩阵函数的求法第46张

二、利用零化多项式求解矩阵函数.

利用Jordan标准型求解矩阵函数的方法比较复杂,它需要求J

P。下面我们介绍根据零化多项式求解矩阵函数的一种方法。

定律:n阶方阵A的最小多项式等于它的特征矩阵的第n个(也就是最后一个)不变因子矩阵理论 第八讲 矩阵函数的求法第47张。(可参见张远达《线性代数原理》P215)

设n阶方阵A的不变因子反向依次为矩阵理论 第八讲 矩阵函数的求法第48张矩阵理论 第八讲 矩阵函数的求法第49张,由它们给出的初等因子分别为

矩阵理论 第八讲 矩阵函数的求法第50张矩阵理论 第八讲 矩阵函数的求法第51张矩阵理论 第八讲 矩阵函数的求法第52张

由于矩阵理论 第八讲 矩阵函数的求法第53张,故

1矩阵理论 第八讲 矩阵函数的求法第54张矩阵理论 第八讲 矩阵函数的求法第55张必定出现在矩阵理论 第八讲 矩阵函数的求法第56张中;

2矩阵理论 第八讲 矩阵函数的求法第57张矩阵理论 第八讲 矩阵函数的求法第58张矩阵理论 第八讲 矩阵函数的求法第59张

根据上述定理,A的最小多项式

矩阵理论 第八讲 矩阵函数的求法第60张

A的最小多项式为其零化多项式,

矩阵理论 第八讲 矩阵函数的求法第61张

矩阵理论 第八讲 矩阵函数的求法第62张,则可见矩阵理论 第八讲 矩阵函数的求法第63张可以由矩阵理论 第八讲 矩阵函数的求法第64张线性表示,从而矩阵理论 第八讲 矩阵函数的求法第65张亦可由矩阵理论 第八讲 矩阵函数的求法第66张线性表示。所以,矩阵函数f(A)若存在,也必定可由矩阵理论 第八讲 矩阵函数的求法第67张线性表示。

因此,我们定义一个系数待定的(m-1)次多项式矩阵理论 第八讲 矩阵函数的求法第68张,根据以上论述,适当选择系数矩阵理论 第八讲 矩阵函数的求法第69张,就可以使f(A)=g(A).

又,假设J,P分别为A的Jordan标准形及相应变换矩阵:矩阵理论 第八讲 矩阵函数的求法第70张

矩阵理论 第八讲 矩阵函数的求法第71张矩阵理论 第八讲 矩阵函数的求法第72张矩阵理论 第八讲 矩阵函数的求法第73张 f(J)=g(J)矩阵理论 第八讲 矩阵函数的求法第74张矩阵理论 第八讲 矩阵函数的求法第75张

矩阵理论 第八讲 矩阵函数的求法第76张矩阵理论 第八讲 矩阵函数的求法第77张矩阵理论 第八讲 矩阵函数的求法第78张

由于g(A)为待定系数的多项式,上面就成为关于矩阵理论 第八讲 矩阵函数的求法第79张的线性方程组。且方程的个数为矩阵理论 第八讲 矩阵函数的求法第80张等于未知数个数,正好可以确定矩阵理论 第八讲 矩阵函数的求法第81张

由此给出根据最小多项式求矩阵函数的一般方法。

1矩阵理论 第八讲 矩阵函数的求法第82张 求出最小多项式

矩阵理论 第八讲 矩阵函数的求法第83张

或者特征多项式矩阵理论 第八讲 矩阵函数的求法第84张

2矩阵理论 第八讲 矩阵函数的求法第85张 形式上写出待定多项式

矩阵理论 第八讲 矩阵函数的求法第86张

或者矩阵理论 第八讲 矩阵函数的求法第87张

3矩阵理论 第八讲 矩阵函数的求法第88张求解关于矩阵理论 第八讲 矩阵函数的求法第89张的线性方程组

矩阵理论 第八讲 矩阵函数的求法第90张矩阵理论 第八讲 矩阵函数的求法第91张

或者矩阵理论 第八讲 矩阵函数的求法第92张

4矩阵理论 第八讲 矩阵函数的求法第93张求出g(A),即可得f(A)=g(A).

从推导的过程看,似乎不仅最小多项式可用于矩阵函数的计算,一般的零化多项式也可以,其中以特征多项式最为方便。

例2. 采用新方法计算矩阵理论 第八讲 矩阵函数的求法第94张的函数 矩阵理论 第八讲 矩阵函数的求法第95张。(矩阵理论 第八讲 矩阵函数的求法第96张

[解] 1矩阵理论 第八讲 矩阵函数的求法第97张矩阵理论 第八讲 矩阵函数的求法第98张矩阵理论 第八讲 矩阵函数的求法第99张

2矩阵理论 第八讲 矩阵函数的求法第100张矩阵理论 第八讲 矩阵函数的求法第101张

3矩阵理论 第八讲 矩阵函数的求法第102张方程组为

矩阵理论 第八讲 矩阵函数的求法第103张矩阵理论 第八讲 矩阵函数的求法第104张

矩阵理论 第八讲 矩阵函数的求法第105张矩阵理论 第八讲 矩阵函数的求法第106张

矩阵理论 第八讲 矩阵函数的求法第107张矩阵理论 第八讲 矩阵函数的求法第108张

4矩阵理论 第八讲 矩阵函数的求法第109张矩阵理论 第八讲 矩阵函数的求法第110张

矩阵理论 第八讲 矩阵函数的求法第111张矩阵理论 第八讲 矩阵函数的求法第112张

 

 

矩阵理论 第八讲 矩阵函数的求法第113张

矩阵理论 第八讲 矩阵函数的求法第114张矩阵理论 第八讲 矩阵函数的求法第115张

与Jordan标准形方法完全一致。

 

 

作业: P163 6

免责声明:文章转载自《矩阵理论 第八讲 矩阵函数的求法》仅用于学习参考。如对内容有疑问,请及时联系本站处理。

上篇NGUI系列教程一[Go] 在golang中使用正则表达式捕获子表达式下篇

宿迁高防,2C2G15M,22元/月;香港BGP,2C5G5M,25元/月 雨云优惠码:MjYwNzM=

相关文章

Laplace算子和Laplacian矩阵

1 Laplace算子的物理意义 Laplace算子的定义为梯度的散度。 在Cartesian坐标系下也可表示为: 或者,它是Hessian矩阵的迹: 以热传导方程为例,因为热流与温度的梯度成正比,那么温度的梯度的散度就是热量的损失率。 由此可见,Laplace算子可用于表现由于物质分布不均引起的物质输送。  2 Laplace算子的数学意义 现在...

矩阵求逆

LuoguP4783 思路: 求A的逆矩阵,把A和单位矩阵I放在一个矩阵里 对A进行加减消元使A化成单位矩阵 此时原来单位矩阵转化成逆矩阵 原理大概就是 A(逆) * [A I] = [I A(逆)] Code: 1 #include <bits/stdc++.h> 2 #define ll long long 3 using namespa...

R数据挖掘 第二篇:基于距离评估数据的相似性和相异性

聚类分析根据对象之间的相异程度,把对象分成多个簇,簇是数据对象的集合,聚类分析使得同一个簇中的对象相似,而与其他簇中的对象相异。相似性和相异性(dissimilarity)是根据数据对象的属性值评估的,通常涉及到距离度量。相似性(similarity)和相异性(dissimilarity)是负相关的,统称为临近性(proximity)。 在聚类分析中,聚类...

采用DCT进行图像压缩

一、实验目的:熟悉变换编码的图像压缩方法 二、实验内容:以所给图像为例,采用dct进行图像压缩编码,通过改变模板矩阵中非零元素的个数,得到不同缩编码图像, 根据公式 ,编写程序计算原图像和dct变换后得到的图像之间的均方误差。用到的matlab函数为im2double,dctmtx,blkproc。 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~...

61 相机投影原理、相机模型中的坐标系统以及标定方法

0 引言   世界坐标系下的点如何投影到CCD镜头上,通过成像的方式得到点在二维图像上的像素坐标值,这是摄影测量中的一个基础而核心的问题。这个问题中核心的东西有两个:1、坐标系的定义及其空间转换矩阵  2、成像中的误差 下面将从这两个角度对这个问题进行详细探讨。文章参考了 # 陈建平: 《相机成像原理》PPThttps://blog.csdn.net/...

OpenCV-C++ Mat对象的使用

这一篇文章中主要介绍有关Mat类的使用,主要包括: 通过构造函数定义Mat对象; Mat对象常用属性以及方法 创建一个空白图像并进行赋值 创建一个小数组 通过构造函数定义Mat对象 Mat中提供了很多的构造函数,这里使用两种方式: Mat(Size size, int type); Mat(int rows, int cols, int type, c...