矩阵有理标准形的性质研究及其应用
- 格式:pdf
- 大小:273.61 KB
- 文档页数:8
第31卷第2期 2015年4月 大 学 数 学
COLLEGE MATHEMATICS Vo1.31,№.2
Apr.2015
矩阵有理标准形的性质研究及其应用 李绍刚 (1.桂林电子科技大学数学与计算科学学院; 2.广西密码学与信息安全重点实验室,广西桂林541004)
[摘要]研究了矩阵的有理标准形理论,给出了友矩阵的主要性质和判断其相似对角化的条件,证明了 线性定常系统的能控性,研究了伪随机数的生成过程和周期性,并对其在自动控制和计算机科学领域的应用 做了一定的解读,最后用实例验证系统的有效性. [关键词]友矩阵;有理标准形;状态方程;能控系统;伪随机数 [中图分类号]O151.21 [文献标识码]C [文章编号]1672—1454(2015)02—0076—08
1 引 言 矩阵的相似标准形理论应用非常广泛,复数域上任何一个矩阵均相似于一个若当标准形,任何数域 上的矩阵则相似于一个有理标准形,可以证明一个复矩阵的有理标准形和其若当标准形是相似的,见文 献[1],由于没有数域上的限制,有理标准形在很多学科得到了广泛的应用,比如在微分方程组的计算 E2],群的构造E3],对于有理标准形的Matlab软件计算,文献[4]给出了相应的算法,在计算机上可以轻 松实现,以上均可以说明其应用是比较广泛的. 本文首先研究有理标准形的友矩阵的性质,其次研究了有理标准形在自动控制领域中状态方程的 计算方面、计算机科学中的伪随机数生成器方面等的应用,最后进行了总结和归纳.
2 基本理论 2.1 基本概念 定义1嘲 对数域P上的一个多项式d( )==: +a +…+n 称矩阵
A:= 为多项式 ( )的伴侣阵,又称作友矩阵. 定义2[ 下列准对角矩阵
B== A1 A2 (2)
[收稿日期]2014一】2—20; [修改日期]2015—03—28 [基金项目]广西区精品课程线性代数项目;线性代数课程与工程应用型人才培养相结合的混合教学模式研究项目
0 O0 11 Oo O第2期 李绍刚:矩阵有理标准形的性质研究及其应用 77 称为P上的一个有理标准形矩阵,其中A 分别是数域P上某些多项式d ( )( 一1,2,…,s)的伴侣阵, 且满足d ( )『dz( )『..·I d ( ),d ( ),d ( ),…,d ( )的次数之和为n.
友矩阵在有理标准形中起着至关重要的作用,其主要性质总结归纳如下:
A= 性质2有理标准形矩阵B的不变因子为1,1,…,1,d ( ),d ( ),…,d ( ),其中1的个数等于 性质3多项式 ( )的友矩阵A的行列式det(A)一(一1) 口 .利用性质3易知
A 一 性质5 多项式 (A)的友矩阵A的特征多项式与最小多项式为 ( )一-厂( )= +al 一 +…+口 . 性质1—5易证,证明从略. 2.3基本定理 定理1 数域P上 ×n方阵A在P上相似于唯一的一个有理标准形B,称为A的有理标准形. 有理标准形的具体计算步骤如下: 步骤1 求出A的不变因子为d ( ),d ( ),…,d ( ); 步骤2对每个次数大于0的不变因子d ( )写出其对应的友矩阵A ( =1,2,… );
步骤3按顺序写出A的有理标准形 A1 A2
利用以上3个步骤可以方便快捷的计算方阵的有理标准形,举例如下: 2.4基本例题 例1 设6阶矩阵A的不变因子为1,1,1, 一1,( 一1)( +2),( 一1) ( +2),求A的有理 标准形. 解 由不变因子可知三个伴侣阵为 r o 1 O] A 一 ], A =[; 一二 ], A。一I。 。 I,
L一2 3 O_I 故可得A的有理标准形为
O 1 2 —1 O 1 0 O O 1 2 3 0
0 o ~ o 0 o ~ ~ ~ ~ ~ ~ o。;o 一 o .。 一 0~0 78 大 学 数 学 第31卷 例2求矩阵 厂0 A—I 3 L一1
1 —1_ 2 0 l 1 —1 l 的有理标准形. 解A的不变因子为 dl( )一1,d2( )一1,d3( )一( 一1)( 十42+2)一 。+32 一22—2. 故A的有理标准形为
3重要应用 3.1在自动控制领域中的应用 3.1.1基本概念与引理 定义3[ 以状态变量 (£),z。( ),…, (f)构成的 维空间称为状态空间. 系统在任意时刻的状态向量 ( )在状态空间中是一个点,系统随时间的变化过程,使 ( )在状态 空间中描绘出一条轨迹. 定义4_6 将反映系统动态过程的微分方程或传递函数,转换成一阶微分方程组的形式,并利用矩 阵和向量的数学工具,将一阶微分方程组用一个式子来表示,这就是状态方程. JC(t)===A ( )+Bu(t), (£)一 ( ), 其中 ( )和主(£)表示状态向量及其一阶导数,U(£)和Y(t)分别表示系统的输人变量和输出变量, A∈ ,口∈ ,C∈ . 系统的能控性指的是控制作用对被控系统状态进行控制的可能性,对于给定的系统称其完全能控 是指对任意初始时刻的任意一个初始状态总可以找到一个容许控制使系统在有限时间内达到目标 状态. 引理1 L6](代数判据) 线性定常系统单输入系统 x( )一A (z)+bu( ), A∈ ,b∈ ”. 其状态完全能控的充条件是由A,b构成的能控性判别矩阵 Q ===(b,Ab,A b…A一 ) 是行满秩矩阵,即有rank(Q )= . 引理2[。 线性定常系统经线性变换后状态能控性保持不变. 引理3[6](模态判据) 线性定常系统单输入系统 主(£)一A ( )+bu(t), AE ,bE吨 . 其状态完全能控的充要条件是系统矩阵A的所有特征值J=【 (i一1,2… )满足 rank(2 E—A, )一72, i=1,2… . 代数判据的优点是方便计算,缺点是无法确定状态空问中哪些变量(特征值、极点)能控;模态判据 的优点是易于确定哪些变量能控,缺点是需要求系统的特征值.这两类方法优势互补,可根据实际情况 选择合适的方法. 3.1.2状态方程的建立 下面以控制系统状态方程为例给出线性定常连续与离散系统状态空间方程的构造过程. 为方便起见,记号如下:
解 求 的 子 因 变 不 阵 矩 于 赖 依 ]●,副算 O 一计 的 1 O 2形 O 0 2准
标
理 有 的 阵 钳 现 发 难 不 题 例 个 两 -Lj 以 从 第2期 李绍刚:矩阵有理标;隹形的性质研究及其应用 79 (£)一 , A一 0 0 O —a (是+1)一 1 0
0 口n 1( +1)
2( +1)
z 一1( +1) (是-I-1)
(忌)一
情形1 线性定常连续系统状态方程 (i) +n1 ‘ 一 +以2 选取 和Y“’为系统状态向量,记
(t)一 . 西一 C—E1,0,0…o l,d—b。. +…+n 一1 +n Y—U. z1一 , z2一 , z3=:: 姑 ,…, = 一”
则上述方程可化为 个一阶常微分方程 Z1 Z2, 2 Z3,…, Z 一a 1一a 一1Z2一a 一2X3…一al ”十 ’ 写成矩阵形式为 f ( )===Ax(t)+6“(£), lY(£)=== (£). (ii)作用函数含导数项
Y +口1 +nz +…+口 一l +口 一6。 +61 十…+6 一1 十6 选取 和Y“’为系统状态向量,记 zl— , 322一 , 3一 ,…,z 一 ’则上述方程可化为 个一阶常微分方程 lz1 2+h1 U, 2一lz3+h2 ,…, zn一一a"z1一a”1z2一a"2z3…一a1 n+h U, 其中 h1一b1一a1bo,h2一(62一a2b0)一nl h1…h 一(6 一a bo)一口 一1h1一…一alh 一1. 矩阵形式为 f (f)一A ( )+bu, l Y(£)===Cx( )+du. 情形2线性定常离散系统状态方程 (i)作用函数不含未来项 (愚+ )+a1 ( + 一1)+…+a 选取状态变量为 1(忌)一.y(是), 2(是)一y(k+1), 根据离散系统差分方程可得状态方程 1 (忌+1)+a Y(矗)一“(是).
,z (志)一y(k+ 一1), z1(忌+1)一zz(忌), z2(是+1)一z。(愚+1),
z (志+1)一一a z1(忌)一a 1z2(忌)一…一a2z 一1( )一a1z (忌)+u(k). 矩阵形式为 f (走+1)一Ax(k+1)+bu(是), Iy(k)===Cx(k). (ii)作用函数含未来项
l h2 : h一1 h
l 2 一 ” z z; ●)) 忌) O O ; 一 1 2 一 " z . O 1;O 一
z;