11.4矩阵的其它运算
第11章 矩阵及其运算
其中aij位于矩阵的第i 行第j 列,称为矩阵 的第i 行第j 列元素,或(i,j)元。
通常我们用大写字母A,B,C 表示矩阵。 为了明确矩阵的行数m 和列数n,我们可用带 下标的大写字母表示矩阵,如Am×n 表示一个 m 行n 列的矩阵。
另外,还可将矩阵表示成A = (aij)或A= (aij)m×n 的形式,其中aij为矩阵的(i,j)元。
本书中矩阵中各元素均为实数。
11.1.2 几种特殊矩阵
1. 方阵
行数与列数相等的矩阵称为方阵,若 行数m 等于列数n,通常称为n 阶方阵或n 阶矩阵,其中行列数n 称为此方阵的阶数。 如
2. 上(下)三角形矩阵
如果n 阶方阵主对角线下方元素都等 于零,则称此矩阵为上三角形矩阵。
如果n 阶方阵主对角线上方元素都等 于零,则称此矩阵为下三角形矩阵。
定义5 设A = (aij)m×n,称A 的最高阶非零 子式的阶数r 为矩阵A 的秩,记为秩(A)=r 或r(A)=r。
若A =0,则其任何阶子式均为零,规定 r(0)=0。
例4 求下列矩阵的秩
解 (1)按定义,找A 的最高阶非零子式,从 低向高找。
定理4 矩阵的初等变换不改变矩阵的秩。
则称这样的阶梯形矩阵为简化的阶 梯形矩阵。
11.4.2 初等矩阵
定义3 n 阶单位矩阵E,经过一次初等变 换所得到的矩阵,称为初等矩阵。
3种初等变换对应3种初等矩阵。
定理1 设矩阵Am×n ,则(1)对A 进行一次 初等行变换,相当于用一个m 阶的同类型 的初等矩阵左乘A 。
(2)对A 进行一次初等列变换,相当于用一 个n 阶的同类型的初等矩阵右乘A 。
例如
则A 为n 阶上三角形矩阵,B 为m 阶 下三角形矩阵。
矩阵的运算
§2 矩阵的运算一、矩阵的相等、加、减、数乘、乘法、转置与共轭(A +B )=A +B (kA )=kA (k 为任意复数) (AB )τ=BA (反序定律)(A 1A 2...A s )=τττ12...A A A s(A k )=(A )k (k 为整数)二、 矩阵的初等变换与初等矩阵设I =⎥⎥⎥⎥⎤⎢⎢⎢⎢⎡10101,称为单位矩阵.用数k(0)乘矩阵的第i 列(或行)初等变换具有性质:1° 任何矩阵(a ij )都可经过有限次初等变换化为对角矩阵(a ij )⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡0001012° 初等变换不改变矩阵的秩.三、 矩阵的微积分假设矩阵A 的元素a ij 都是参数t 的函数,那末1° 矩阵A 的导数定义为⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡==t a t a ta t a t a tat a t a t a A tA mn m m n n d d ...d d d d ............d d ...d d d d d d ...d d d d d d 212222111211同样可定义矩阵的高阶导数. 2° 矩阵A 的积分定义为⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡=⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰t a t a ta t at at a t a t a ta t A mn m m n nd ...d d ............d ...d d d ...d d d 212222111211同样可定义矩阵的多重积分.四、 特殊矩阵[零矩阵与零因子] 元素a ij 全为零的矩阵称为零矩阵,记作O =(0)=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡0...00............0 (00)0 (00)零矩阵具有性质:O +A =A +O =A OA =AO =OA +(-A )=O ,-A 称为A 的负矩阵若A ,B 为非零矩阵,即A ≠O ,B ≠O ,而AB =O ,则称矩阵A 为矩阵B 的左零因子,矩阵B 为矩阵A 的右零因子,例如A =⎥⎦⎤⎢⎣⎡--1111,B =⎥⎦⎤⎢⎣⎡--1111 AB =⎥⎦⎤⎢⎣⎡--1111⎥⎦⎤⎢⎣⎡--1111=⎥⎦⎤⎢⎣⎡0000=O[对角矩阵] 主对角线以外的元素都是零(d ij =0,i ≠j )的方阵称为对角矩阵,记作D =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡n d d d 0...021=diag(d 1,d 2,...,d n )=[ d 1 d 2 ... d n ] 对角矩阵具有性质: 1° 左乘BDB =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡n d d d 0021⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n n n b b b b b b b b b .....................212222111211=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n n n n n n b d b d b d b d b d b d b d b d b d ............... (2)12222221211121111 =)(ij i b d 2° 右乘BBD =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n n n b b b b b b b b b (2)12112111211⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡n d d d 0021=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n n n n n n b d bd b d b d b d bd b d b d b d (2211222)22111122111 3° 两个对角矩阵的和、差、积仍为对角矩阵.[数量矩阵] d i =d (i =1,2,...,n )的对角矩阵称为数量矩阵,记作D =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡d d d00 =[d d... d ]显然DB =BD =dB .[单位矩阵] d =1的数量矩阵称为单位矩阵,记作 I =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡10101 =「1 1 ... 1」显然IB =BI =B .[对称矩阵] 满足条件a ij =a ji (i ,j =1,2,...,n )的方阵A =(a ij )称为对称矩阵.例如A =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--423261315 是对称矩阵.对称矩阵具有性质: 若A ,B 都是对称矩阵,则A A=τ,且A -1(使A -1=A -1A =I 的矩阵.详见本节,六),A m (m 为正整数),A +B 仍是对称矩阵.[实对称矩阵]实对称矩阵按其特征值(本节,七)可分为正定矩阵,半正定矩阵、负定矩阵、半负定矩阵和不定矩阵,它们的定义与充分必要条件如下[反对称矩阵] 满足条件⎩⎨⎧-=jiij a a 0 )()(j i j i ≠= (i ,j =1,2,...,n )的方阵A =(a ij )称为反对称矩阵.例如A =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---023201310 是反对称矩阵.反对称矩阵具有性质:1° 若A ,B 都是反对称矩阵,则A τ=-A ,且A -1, A +B 仍是反对称矩阵,A m 为⎩⎨⎧反对称矩阵对称矩阵)()(为奇数为偶数m m2° 任意方阵A 都可分解为一个对称矩阵B =(b ij )与一个反对称矩阵C =(c ij )之和,即A =B +C只需取b ij =21 (a ij +a ji ),c ij =21(a ij -a ji )(i ,j =1,2,...n )[埃尔米特矩阵] 满足条件A τ=A的方阵A 称为埃尔米特矩阵.例如A =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+--++-4232231212215i i i i i i 是埃尔米特矩阵.埃尔米特矩阵具有性质:若A ,B 都是埃尔米特矩阵,则1-A ,A +B 仍是埃尔米特矩阵.若A 又是实方阵(即a ij 全为实数),则A 就是对称矩阵.[反埃尔米特矩阵] 满足条件A τ=A -的方阵A 称为反埃尔米特矩阵.例如A =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-+--+-05250212210i i i i i i 是反埃尔米特矩阵.反埃尔米特矩阵具有性质: 若A ,B 都是反埃尔米特矩阵,则1-A , A +B 仍是反埃尔米特矩阵.若A 又是实方阵,则A 就是反对称矩阵.[正交矩阵] 满足条件A τ=1-A的方阵A 称为正交矩阵.例如 A =⎥⎦⎤⎢⎣⎡-θθθθcos sin sin cos 是正交矩阵.正交矩阵具有性质:若A =(a ij )和B 都是正交矩阵,则 1° 1-A , AB 仍是正交矩阵. 2° det A =±1.3° ⎩⎨⎧=∑=011n k jk ik a a )()(j i j i ≠=⎩⎨⎧=∑=011n k kj ki a a )()(j i j i ≠=[酉(U )矩阵] 满足条件1-=A A τ的方阵A 称为酉(U )矩阵.例如:A =⎥⎦⎤⎢⎣⎡00i i 是酉矩阵.酉矩阵具有性质:若A =(a ij )和B 都是酉矩阵,则 1° A -1,AB 仍是酉矩阵. 2° det A ∙det A =1.3° 若A 又是实方阵,则A 是正交矩阵.[带型矩阵] 满足条件a ij =0 )(m j i >-的方阵A =(a ij )称为带型矩阵.2m +1称为带宽.一般形式为A =⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--++++nn mn n n m n n n n m a a a a a a a,,1,11,11,11100[三角矩阵] 满足条件a ij =0 (i >j )的方阵A =(a ij )称为上三角形矩阵,一般形式为A =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n a a a a a a 022211211 满足条件()j i b ij <=0的方阵()ij b B =称为下三角形矩阵,一般形式为B =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n b b b b b b 212221110 三角形矩阵具有性质:1° 任何秩为r 的方阵C 的前r 个顺序的主子式不为0时,C 可表为一个上三角形矩阵A与一个下三角形矩阵B 的乘积,即C =AB2° 上(或下)三角形矩阵的和、差、积及数乘仍是上(或下)三角形矩阵.[分块矩阵] 用水平和垂直虚线将矩阵A 中的元素的阵列分成小块(称为子阵),A 就成为分块矩阵.例如A =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡333231232221131211a a a a a a a a a =⎥⎦⎤⎢⎣⎡22211211B B B B 式中B 11=⎥⎦⎤⎢⎣⎡22211211a a a a,B 12=⎥⎦⎤⎢⎣⎡2313a a B 21=[]3231a a , B 22=[]33a 它们都是A 的子阵. 进行分块矩阵的运算时,可将子阵当作通常矩阵的元素看待.这些运算指加、减、乘法、数乘、转置与共轭等.[分块对角矩阵] 主对角线上的子阵都是方阵,其余子阵都是零矩阵的分块矩阵称为分块对角矩阵.一般形式为A =⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡kkB O B O O O B 2211 分块对角矩阵A 的逆矩阵A -1和A 的行列式可以用下面简单公式求出A -1=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡---1122111KK B OB O Bdet A =det B 11·det B 22·...·det B kk注意,一般分块矩阵的行列式不能用把子阵当作通常矩阵的元素的方法来计算,例如把四阶方阵化为分块矩阵A =⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡44434241343332312423222114131211...........................a a a a a a a a a a a a a a a a =⎥⎦⎤⎢⎣⎡22211211B B B B 一般det A =det B 11·det B 22-det B 21·det B 12不成立(参见§1,二,3中的四阶行列式).五、 相似变换[相似变换] 如果有一非奇异矩阵X (即det X ≠0)使得B =1-X AX那末称矩阵A 与矩阵B 相似,也称A 经相似变换化为B ,记作A ~B .它具有下列性质: 1° A ~A ,AA .2° 若A ~B ,则BA .3° 若A ~C ,B ~C ,则A ~B .4° 1-X (A 1+ A 2+...+ A m )X =1-X A 1X + 1-X A 2X + ...+ 1-X A m X 5° 1-X (A 1 A 2 ...A m )X =1-X A 1 X ·1-X A 2 X ·... ·1-X A m X 6° 1-X A m X =( 1-X AX )m7° 若)(A f 为矩阵A 的多项式,则1-X )(A f X =)(1AX X f -8° 若A ~B ,则A 与B 的秩相同,即rank A =rank B . A 与B 的行列式相同,即det A =det B .A 与B 的迹(定义见本节,七)相同,即tr A =tr B . A 与B 具有相同的特征多项式和特征值(本节,七).[正交变换] 若Q 为正交矩阵(即1-Q =Q τ),则称Q τAQ 为矩阵A 的正交变换,其性质与相似变换类似.特别还有性质: 对称矩阵A 经正交变换后仍是对称矩阵.[旋转变换] 取正交矩阵U 为)(p)(qU pq =(u ij )=)()(11cos sin 11sin cos 11q p ⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡θθ-θθ 即u pp =u qq =θcosu pq =-u qp =θsin u ii =1 (i ≠p,q )u ij =0 (i,j ≠p,q;i ≠j ) 这时称B =pq pq AU U τ为A 的旋转变换,称为旋转角,如果A 是对称矩阵,那末B 的元素b ij 与A 的元素a ij 有 如下对应关系:⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=θ+θ=θ-θ=θ-θ+θθ-==θ+θθ+θ=θ+θθ-θ=ijijqj pj qj qj pj pj pq qq pp qp pqqq pq pp qq qq pq pp pp a b a a b a a b a a a b b a a a b a a a b cos sin sin cos )sin (cos cos sin )(cos cos sin 2sin sin cos sin 2cos 222222)其他元素(),(),(q p j q p j ≠≠同时有性质:∑=nj i ija1,2=∑=nj i ij b 1,2∑=ni iia 12∑=≤ni ii b 12 若取旋转角pqpp qq a a a 2cot arc 21-=θ则旋转变换使0==qp pq b b六、 逆矩阵[逆矩阵及其性质] 若方阵A ,B 满足等式AB=BA=I (I 为单位矩阵)则称A 为B 的逆矩阵,或称B 为A 的逆矩阵,记作A=1-B 或B=1-A这时A,B 都称为可逆矩阵(或非奇异矩阵,或满秩矩阵).否则称为不可逆矩阵(或奇异矩阵,或降秩矩阵).可逆矩阵具有性质:1° 若A,B 为可逆矩阵,则AB 仍为可逆矩阵,且111)(---=A B AB (反序定律)一般地,若A 1 ,A 2 ,…,A s 为可逆矩阵,则=-121)(s A A A 11121---A A A s2° 矩阵A 可逆的充分必要条件是:det A ≠0.3° 若矩阵A 可逆,则det 1-A ≠0 且 det 1-A =(det 1)-A11)(--A =A , 111)(---=A a aA (a ≠0)1)(-τA =(1-A )τ,()()11--=A A4° 矩阵A 可逆的充分必要条件是:矩阵A 的特征值全不为零.[伴随矩阵与逆矩阵表达式] 设A ij 为矩阵A =(a ij )的第i 行第j 列元素a ij 的代数余子式,则矩阵A *=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn nnn n A A A A A A A A A (2122212)12111称为矩阵A 的伴随矩阵.若A 为非奇异矩阵,即det A ≠0,则A 的逆矩阵表达式为AA A det *1=-注意,A *的第i 行第j 列元素是A 的第j 行第i 列元素的代数余子式.[对角矩阵的逆矩阵] 对角矩阵D =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡n d d d 0...021, d i ≠0 (i =1,2,...,n )的逆矩阵为D -1=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡---112110...0n d d d 显然对角矩阵的逆矩阵仍是对角矩阵.[三角形矩阵的逆矩阵] 三角形矩阵L =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n l l l l l l ...............0...0...21222111, 00=≠ij ii l l )(),...,2,1(i j n i >= 的逆矩阵为1-L =P =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n p p p p pp ...............0...0 (02)1222111 式中iiii l p 1=(i =1,2,...,n )∑-=-=11i jk kj ikiiij p ll p⎪⎪⎭⎫ ⎝⎛+=-=n j i n j ,...,11,...,2,1 0=ij p)(i j >显然非奇异下(上)三角形矩阵的逆矩阵仍是下(上)三角形矩阵.[正定矩阵的逆矩阵] 1° 高斯—若当法正定矩阵A =(a ij )的逆A -1=(b ij )可由下列递推公式求出:)1(11)(1-=k k nnaa, )1(11)1(1)(1,----=k k jk j n aa a, )1(11)1(1)(,1---=k k i k ni a a a)1(11)1(1)1(1)1()(1,1-------=k k jk i k ij k j i aa a a a )2,...,1,,(-=n n j i ij n ij a a =)((k=1,2,...,n )最后得到)(n ijij a b = 式中n 为该正定矩阵A 的阶. 2° 三角阵法 其步骤如下:(1) 把正定矩阵A =(a ij )表示为A =ΛD Λτ式中D 为实的非奇异对角矩阵D =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡n d d d 0021为实的非奇异下三角矩阵.Λ=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡λλλλ-1111,2121n n n n是的转置矩阵.d i (i =1,2,...,n )与λij (i =2,...,n;j=1,…,n )由下面递推公式算出:0=ij λ)(i j > 1=λii ),...,2,1(n i =∑-=-=11j k jk ik ij ij x a x λ)1,...,2,1;,...,2(-==i j n ijij ij d x =λ)1,...,2,1;,...,2(-==i j n i∑-=-=11i k ik ik ii i x a d λ),...,2,1(n i =(2)求出D 的逆矩阵1-D =⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡n d d d 11121(3)求出Λ的逆矩阵1-Λ=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡1112121 n n ρρρ 式中⎪⎩⎪⎨⎧=-=∑-=11ii i jk kjik ij ρρλρ ),...,2,1(),...,2,1;1,...,2,1(n i n j j i n j =++=-=(4)求出A 的逆矩阵1-A =(ΛD 1)-τΛ=(1-Λ)τ1-D 1-Λ =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n n n βββββββββ212222112111式中∑==nik kkjki ij d ρρβ ),,2,1;,,2,1(n i i j ==注意,这种方法的好处是避免了求平方根的运算.[分块矩阵的逆矩阵] 设非奇异矩阵A 的分块矩阵为A =⎥⎦⎤⎢⎣⎡22211211B B B B 式中B 11,B 22为方子阵,那末A 的逆矩阵A -1=⎥⎦⎤⎢⎣⎡22211211C C C C由下面公式求出111211211111111212221221211112112111212222)(-------=-=-=-=B B C B C B B C C C B B C B B B B C[初等变换法求逆矩阵] 设1-A =1212222111211...........................-⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n n n a a a a a a a a a =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n n n b b b b b b b b b 212222111211=B 对矩阵⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡100010001212222111211 nn n n n n a a a a a a a a a 作一系列行的初等变换,使虚线左边一块矩阵化为单位矩阵,而右边一块单位矩阵就变为A 的逆矩阵B =A -1,即⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n n n b b b b b b b b b212222111211100010001[逆矩阵的近似求法] 设10-A 为矩阵A 的初始近似逆矩阵,可由下列迭代公式求出更精确的逆矩阵:)2(1111---+-=n n n AA I A A (n=0,1,2,...)式中I 为与A 同阶的单位矩阵.[计算机求逆程序的检验矩阵] 用下列n 阶非奇异矩阵及其逆矩阵,来检验大矩阵求逆的计算程序.A =⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡++-+------+-++222210221211210002112100002112122100021222n n n n n n1-A =⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--------n n n n n n n n n n n n n13211432341223111221七、 特征值与特征矢量[特征值与特征矢量] 对n 阶方阵A =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n n n a a a a a a a a a 212222111211 和n 维非零列矢量α=(a 1,a 2,...,a n )τ如果有一个数λ,使得A α=λα则称λ为矩阵A 的特征值(特征根),α为矩阵A 的特征值λ所对应的特征矢量. 矩阵A 的所有特征值中绝对值最大的一个称为A 的第一特征值.[特征矩阵特征多项式特征方程] n 阶方阵A =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n n n a a a a a a a a a 212222111211 的特征矩阵定义为=-I A λ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---λλλnn n n n n a a a a a aa a a212222111211 式中I 为n 阶单位矩阵.行列式|A -λI |称为矩阵A 的特征多项式,记作()=|-A λI |方程()=0称为矩阵A 的特征方程.[矩阵的迹与谱] n 阶方阵A 的主对角线上各元素之和称为A 的迹,记作∑==ni ii a A 1tr特征方程()=0的n 个根1,2,...,n 就是矩阵A 的n 个特征值.集合{1,2,...,n }称为矩阵A 的谱,记作ch A .线性齐次方程组0)(=-αλI A i的非零解便是矩阵A 的特征值i 所对应的特征矢量.[特征值与特征矢量的性质]1° 设1,2,...,n 为n 阶方阵A 的n 个特征值,则A k 的特征值为k n k k λλλ,,,21 (k 为正整数). A 的逆矩阵A -1的特征值为11211,,,---n λλλ .A 的伴随矩阵A *的特征值为A A A n 11211,,,---λλλ .2° n 阶方阵A 的n 个特征值之和等于A 的迹,矩阵A 的n 个特征值之积等于A 的行列式,即1+2+...+n =a 11+a 22+...+a nn12...n =A由此可以推出矩阵可逆的另一充分必要条件是:A 的所有特征值都不为零. 3° 若i 是特征方程的k 重根,则对应于i 的线性无关的特征矢量的个数不大于k .当i 为单根时,对应于i 的线性无关特征矢量只有一个.4° 矩阵A 的不同特征值所对应的特征矢量线性无关.若n 阶方阵A 对应于特征值1,2,...,s 的线性无关的特征矢量分别有k 1,k 2,...,k s个,则这∑=s i i k 1个特征矢量线性无关,且n k si i ≤∑=1.5° 实对称矩阵的特征值都是实数,并且有 n 个线性无关(而且是正交)的特征矢量. 6° 矩阵的特征值在相似变换下保持不变,特别,A τ与A 具有相同的特征值.[求第一特征值的迭代法] 在实际问题中,往往不要求算出矩阵A 的全部特征值,只需算出第一特征值,用迭代法计算如下:⎩⎨⎧=λ=α++b αα)0()1()1(1)(k k k A )2,1,0( =k 假定当ε<-+)1()(m m αα时,可以认为(k ) ≈(m +1),那末迭代到m k =即可.这时)1(1+m λ为矩阵A 的第一特征值的近似值,(m +1)为所对应的特征矢量.[求实对称矩阵的雅可比法] 设n 阶实对称矩阵A =(a ij )的特征值是1,2,...,n ,则必存在一正交矩阵Q ,使得Q τAQ =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡λλλn 0021为对角矩阵.正交矩阵Q 可用一系列旋转矩阵的积来逼近:Q =∏pq U式中)()(11cos sin 11sin cos 11)()()(q p u U q p ij pq⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-==θθθθ取pqpp qq a a a 2cot arc 21-=θ因为在这种旋转变换下,消去了矩阵中位于第p 行第q 列(p ≠q )交点上的元素(见本节,五),而矩阵所有元素的平方和保持不变,而且对角线上的元素的平方和增大,因而非对角线元素的平方和随之减小,因此,当旋转次数足够大时,可使非对角线元素的绝对值足够小.对于预先给定的精度>0,如果|a ij |<(i ≠j ),则可认为a ij ≈0.于是得到求矩阵A 的特征值与特征矢量的具体迭代方法.1° 按以下递推公式求特征值1,2,...,n :⎪⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎪⎨⎧+=θ=⎪⎪⎩⎪⎪⎨⎧<+->-+=θ=⎪⎩⎪⎨⎧<ςς++ς-≥ςς++ς=θ=-=θ=ς--2221212)()()(1sin )0(11)0(112tan )0()1()0()1(tan 22cot k k k k k k k k k kk k k k k k k k pq k pp k qq k t t s t t t t t t v t a a a⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧===≠≠=≠-+=≠+-=+=-=+++++),2,1(),,2,1,(),,,()()()()()1()1()()()()1()()()()1()()()1()()()1( k n j i a a q p j q p i a a q j a a s a a p j a a s a a a t a a a t a a ij ij kijk ijk qj k k pj k k qj k qj k pj k k qj k k pj k pj k pqk k qq k qq k pqk k pp k pp υυ假定当)()(j i a m ij ≠<ε时,可以认为0)(≈m ij a ,则迭代到1-=m k 即可.而取)(m iia 作为i的近似值:),,2,1(n i a miii =≈λ2° 求特征矢量 从1°有m m m m U U AU U U U 1111-- τττ=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡n λλλ0021记P m =U 1…U m-1U m则AP m = P m ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡n λλλ0021所以P m 为特征矢量矩阵.P m 由下列递推公式算出:)1,,2,1(),,2,1,(),,2,1(),()()()1()()1()()()()1()()()()1(-=⎪⎪⎪⎩⎪⎪⎪⎨⎧===≠=-+=+-=+++m k n j i u u n i q p j u u u u s u u u u s u u ijij k ijk ij k iq k k ip k k iq k iq k ip k k iq k k ip k ip υυ最后得到 )()(m ij m u P =即 τ),,,()()(2)(1)(m ni m i m i m i u u u u =为对应于特征值i 的特征矢量的近似值.[求对称三对角矩阵特征值的方法]1° 相似变换法 设A 为n 阶对称三对角矩阵:A =⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--n n n d e e d e e d e e d 113222111(1)经过相似变换1211211)(U U U I t A U U U A n k k n k --+-=τττ式中I 为单位矩阵,t k 为适当选定的常数,U i 为雅可比旋转矩阵:)1()(1111)1()(+⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-=+i i c s s c U i i ii i iiτi U 为U i 的转置矩阵.又A 1=A ,A k +1与k k t A -I 相似,且A m 与∑-=-111m j j I t A 相似.因此,若A m 的特征值为),,2,1()(n i m i =λ,则A 1的特征值i (i=1,2,...,n )为∑-=+=11)(m j j m ii t λλ(i =1,2,…,n )假定当),,2,1()(n i e m i =<ε时,可认为0)(≈m i e ,那末可适当选择s i ,c i ,使得当m 充分大时,A m 在该精度下化为对角线矩阵;其特征值),,2,1()()(n i d m i m i =≈λ.)(m i d (i=1,2,...,n )可由下列递推公式算出:()())1,,2,1;1,2,,2,1(,)]([)(//g ])()[(0,,)(1)(1)1(1)(1)(1)1(1)(1)(1)1(1)()()(1)()()(1)1(1)(1)()()()()(1)()()(1)(1)()(1)(1(k)1)()(1(k)1212)(2)(1)(1)()(-=--=⎪⎩⎪⎨⎧===-++=--=====+==-=+++++++++++++++++++++m k n n i q s e q c d r s e t d s g c s h d g s t d c q r e s r q c q c h e c c q rs c t d q k k k k k k k i k i k i k k i k i k i k i k i k i k i k i k i k k i k i k i k i k i k i k i k i k i k i k i i k i k i i k ik i k i k nk n k k n k nt k 的选择对收敛速度影响较大,取t k 为二阶矩阵⎥⎦⎤⎢⎣⎡)(2)(1)(1)(1k k k k d e e d 的接近于)(1k d 的那个特征值,即t k =⎪⎩⎪⎨⎧≥ββ++β-<ββ+-β-)0()1/()0()1/(2)(1)(12)(1)(1k k k k e d e d式中 )(1)(1)(22k k k e d d -=β 2° 二分法 设A 为n 阶对称三对角矩阵(如(1)式),对任意,设序列q 1()=d 1-q i ()=),,2()()(121n i q e d i i i =----λλ中q i ()<0的个数为N ()(在这些关系式中,对于某些i ,如果q i -1()=0,则只需用适当小的数代替即可),则N ()等于矩阵A 的小于的特征值的个数.假定矩阵A 的第k 个特征值k (1≤2≤… ≤k ≤…≤n )在区间[u ,υ]中,令21υ+=u r ,当N (r 1)≥k 时,则k ∈[u , r 1];当N (r 1)<k 时,则k ∈[ r 1,v ];…依此类推,m步之后,k 包含在宽度为mu2-υ的区间中.m 充分大时,便可得到所求的特征值.八、 矩阵多项式与最小多项式[矩阵多项式] 设i a (i=1,2,...,n )为某一数域(实数域或复数域)中的数,A 为这个数域上的n 阶方阵,则表示式f (A )=a 0I+a 1A+...+a n A n称为矩阵A 的多项式,式中I 为n 阶单位矩阵.如果矩阵A 使得f (A )=O那末称A为多项式f(λ)=a0λ+ a1λ+ ...+a nλn的根.[哈密顿-凯莱定理] 任一方阵都是它的特征多项式的根.[最小多项式及其性质] 以矩阵A为根的非零多项式f(λ)中,存在首项系数为1次数最低的多项式(λ),它就称为矩阵A的最小多项式.最小多项式具有性质:1°任一方阵仅有一个最小多项式;2°任一以A为根的多项式f(λ)都可被A的最小多项式(λ)所整除.特别,任一方阵的最小多项式可整除其特征多项式;3°方阵A的特征多项式的根都是A的最小多项式的根:4°相似矩阵具有相同的特征多项式和最小多项式.。
《矩阵及其运算 》课件
幂法
通过迭代计算矩阵A的幂 ,最终得到特征值和特征 向量。
反迭代法
利用已知的特征向量x, 通过反迭代计算得到对应 的特征值λ。
06
应用实例
在物理中的应用
线性变换
矩阵可以表示线性变换,如平移、旋转、缩放等,在物理中广泛应 用于描述物体运动和力的作用。
振动分析
矩阵可以用于分析多自由度系统的振动,通过矩阵表示系统的运动 方程,简化计算过程。
详细描述
矩阵乘法要求第一个矩阵的列数等于第二个矩阵的行数,并 且结果矩阵的行数等于第一个矩阵的行数,列数等于第二个 矩阵的列数。在计算过程中,对应元素相乘并求和,得到新 矩阵的一个元素。
矩阵的转置
总结词
矩阵的转置是将原矩阵的行变为列,列变为行的一种运算。
详细描述
矩阵的转置可以通过交换原矩阵的行和列得到,也可以通过计算元素的代数余 子式得到。转置后的矩阵与原矩阵的行列式值相等,但元素的位置发生了变化 。
《矩阵及其运算》PPT课件
目 录
• 矩阵的定义与性质 • 矩阵的运算 • 矩阵的逆与行列式 • 矩阵的秩与线性方程组 • 特征值与特征向量 • 应用实例
01
矩阵的定义与性质
矩阵的基本概念
矩阵的定义
矩阵是一个由数字组成的矩 形阵列,通常表示为二维数 组。
矩阵的元素
矩阵中的每个元素都有行标 和列标,表示其在矩阵中的 位置。
回带法
在消元过程中,每一步都需要回带, 以确保解的正确性。
解的判定
当系数矩阵的秩等于增广矩阵的秩时 ,线性方程组有唯一解;否则,无解 或有无数多解。
线性方程组的解的结构
解的表示
线性方程组的解可以表示为一个向量与自由变量 的线性组合。
矩阵的运算及其运算规则
矩阵基本运算及应用201700060牛晨晖在数学中,矩阵是一个按照长方阵列排列的复数或实数集合。
矩阵是高等代数学中的常见工具,也常见于统计分析等应用数学学科中。
在物理学中,矩阵于电路学、力学、光学和量子物理中都有应用;计算机科学中,三维动画制作也需要用到矩阵。
矩阵的运算是数值分析领域的重要问题。
将矩阵分解为简单矩阵的组合可以在理论和实际应用上简化矩阵的运算。
在电力系统方面,矩阵知识已有广泛深入的应用,本文将在介绍矩阵基本运算和运算规则的基础上,简要介绍其在电力系统新能源领域建模方面的应用情况,并展望随机矩阵理论等相关知识与人工智能电力系统的紧密结合。
1矩阵的运算及其运算规则1.1矩阵的加法与减法1.1.1运算规则设矩阵,,则简言之,两个矩阵相加减,即它们相同位置的元素相加减!注意:只有对于两个行数、列数分别相等的矩阵(即同型矩阵),加减法运算才有意义,即加减运算是可行的.1.1.2运算性质满足交换律和结合律交换律;结合律.1.2矩阵与数的乘法1.2.1运算规则数乘矩阵A,就是将数乘矩阵A中的每一个元素,记为或.特别地,称称为的负矩阵.1.2.2运算性质满足结合律和分配律结合律:(λμ)A=λ(μA);(λ+μ)A =λA+μA.分配律:λ(A+B)=λA+λB.1.2.3典型举例已知两个矩阵满足矩阵方程,求未知矩阵.解由已知条件知1.3矩阵与矩阵的乘法1.3.1运算规则设,,则A与B的乘积是这样一个矩阵:(1) 行数与(左矩阵)A相同,列数与(右矩阵)B相同,即.(2) C的第行第列的元素由A的第行元素与B的第列元素对应相乘,再取乘积之和.1.3.2典型例题设矩阵计算解是的矩阵.设它为可得结论1:只有在下列情况下,两个矩阵的乘法才有意义,或说乘法运算是可行的:左矩阵的列数=右矩阵的行数;结论2在矩阵的乘法中,必须注意相乘的顺序.即使在与均有意义时,也未必有=成立.可见矩阵乘法不满足交换律;结论3方阵A和它同阶的单位阵作乘积,结果仍为A,即.1.3.3运算性质(假设运算都是可行的)(1) 结合律.(2) 分配律(左分配律);(右分配律).(3) .1.3.4方阵的幂定义:设A是方阵,是一个正整数,规定,显然,记号表示个A的连乘积.1.4矩阵的转置1.4.1定义定义:将矩阵A的行换成同序号的列所得到的新矩阵称为矩阵A的转置矩阵,记作或.例如,矩阵的转置矩阵为.1.4.2运算性质(假设运算都是可行的)(1)(2)(3)(4) ,是常数.1.4.3典型例题利用矩阵验证运算性质:解;而所以.定义:如果方阵满足,即,则称A为对称矩阵.对称矩阵的特点是:它的元素以主对角线为对称轴对应相等.1.5方阵的行列式1.5.1定义定义:由方阵A的元素所构成的行列式(各元素的位置不变),称为方阵A的行列式,记作或.1.5.2运算性质(1) (行列式的性质)(2) ,特别地:(3) (是常数,A的阶数为n)思考:设A为阶方阵,那么的行列式与A的行列式之间的关系为什么不是,而是?不妨自行设计一个二阶方阵,计算一下和.例如,则.于是,而2光伏逆变器的建模光伏并网逆变器是将光伏组件输出的直流电转化为符合电网要求的交流点再输入电网的关键设备,是光伏系统并网环节中能量转换与控制的核心。
矩阵的简单运算公式
矩阵的简单运算公式矩阵是线性代数中的重要概念,广泛应用于数学、物理、计算机等各个领域。
矩阵的运算涉及到加法、减法、数乘和乘法等操作,下面将介绍一些简单的矩阵运算公式。
1. 矩阵加法矩阵加法是指两个矩阵按照相同位置的元素进行相加的运算。
设矩阵A和矩阵B分别为m行n列的矩阵,其加法公式为:C = A + B其中C为相加后的结果矩阵,C的每个元素等于A和B对应位置元素的和。
2. 矩阵减法矩阵减法是指两个矩阵按照相同位置的元素进行相减的运算。
设矩阵A和矩阵B分别为m行n列的矩阵,其减法公式为:C = A - B其中C为相减后的结果矩阵,C的每个元素等于A和B对应位置元素的差。
3. 数乘数乘是指将矩阵的每个元素乘以一个常数。
设矩阵A为m行n列的矩阵,k为常数,其数乘公式为:C = kA其中C为数乘后的结果矩阵,C的每个元素等于k乘以A相应位置的元素。
4. 矩阵乘法矩阵乘法是指两个矩阵按照一定规律进行的乘法运算。
设矩阵A为m行p列的矩阵,矩阵B为p行n列的矩阵,其乘法公式为:C = AB其中C为乘法的结果矩阵,C的第i行第j列的元素等于矩阵A的第i行与矩阵B的第j列的对应元素的乘积之和。
以上是矩阵的几种简单运算公式,在实际运用中可以通过这些公式进行各种复杂的矩阵运算。
矩阵运算在线性代数、图像处理、数据分析等领域具有广泛的应用,依靠这些运算公式可以很方便地对矩阵进行操作和计算。
需要注意的是,在进行矩阵运算时,要确保参与运算的矩阵具有相同的行列数,否则运算无法进行。
此外,矩阵运算具有交换律、结合律和分配律等基本性质,可以根据需要灵活运用。
总之,矩阵的简单运算公式包括加法、减法、数乘和乘法等操作,这些公式可以帮助我们对矩阵进行各种运算和计算。
掌握这些运算公式,并善于应用,将会对求解复杂问题起到很大的帮助作用。
结构力学11.4 结构的原始刚度矩阵
§11-4结构的原始刚度矩阵整体分析:即建立求解基本未知量的结构刚度方程。
而位移法中求解的基本未知量是结点位移,包括线位移与角位移。
考虑如上图所示刚架,有4个结点,3个单元,受结点荷载的作用。
至于非结点荷载作用的情况,需要将其转化为等效的结点荷载,这将在后面的内容中进行专门介绍。
这里,暂只考虑结点荷载作用的情况。
各单元的局部坐标系与整体坐标系如下图所示。
各单元的单元刚度矩阵,进行坐标变换后,得到整体坐标系下各单元的单元刚度矩阵如下1221][22211211⎥⎦⎤⎢⎣⎡=①①①①①k k k k k ,2332][33322322⎥⎦⎤⎢⎣⎡=②②②②②k k k k k ,43][44433433⎥⎦⎤⎢⎣⎡=③③③③③k k k k k 34每个结点,有x 方向线位移、y 方向线位移与角位移3个位移分量。
结构的结点位移列向量为⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=1111}{ϕv u Δ,⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=2222}{ϕv u Δ,⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=3333}{ϕv u Δ,⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=4444}{ϕv u Δ⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧=4321}{ΔΔΔΔΔ与结点位移列向量对应的结点外力(包括荷载和反力)列向量为结构刚度方的建立前面学习位移法时,已经知道,位移法方程,即结构的刚度方程,就是结点的平衡方程。
所以,通过结点的平衡条件,可建立结构的刚度方程。
下面,以结点2为例,如图示。
结点2的3个平衡方程为②①222x x x F F F +=,②①222y y y F F F +=,②①222M M M +=写成矩阵形式,有⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧+⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧②②②①①①222222222M F F M F F M F F y x y x y x ,即,}{}{}{222②①F F F +=单元杆端力,可用杆端位移表示为}]{[}]{[}{2221212①①①①①δδk k F +=,}]{[}]{[}{3232222②②②②②δδk k F +=得到,}]{[}]){[]([}]{[}{323222221212Δk Δk k Δk F ②②①①+++=此即结点2的平衡方程。
矩阵常见运算
矩阵的基本运算公式加法,减法,数乘,转置,共轭和共轭转置。
1、矩阵的加法满足A+B=B+A;(A+B)+C=A+(B+C)。
在两个数的加法运算中,在从左往右计算的顺序,两个加数相加,交换加数的位置,和不变。
A+B+C=A+C+B。
加法定理一个是指概率的加法定理,讲的是互不相容事件或对立事件甚至任意事件的概率计算方面的公式;另一个是指三角函数的加法定理。
2、把矩阵A的行和列互相交换所产生的矩阵称为A的转置矩阵,这一过程称为矩阵的转置。
设A为m×n阶矩阵(即m行n列),第i 行j 列的元素是a(i,j),即:A=a(i,j)定义A的转置为这样一个n×m阶矩阵B,满足B=b(j,i),即a(i,j)=b (j,i)(B的第i行第j列元素是A的第j 行第i列元素),记A'=B。
3、矩阵乘法是一种根据两个矩阵得到第三个矩阵的二元运算。
二元运算属于数学运算的一种。
二元运算需要三个元素:二元运算符以及该运算符作用的两个变量。
如四则运算的加、减、乘、除均属于二元运算。
如在运算1 + 2之中,二元运算符为“+”,而该运算符作用的操作数分别为1与2。
二元运算只是二元函数的一种,由于它被广泛应用于各个领域,因此受到比其它函数更高的重视。
矩阵运算公式大全
矩阵运算公式大全一、矩阵基本概念和性质1.矩阵的定义:一个m×n的矩阵A是由m行n列的数排成的一个矩形阵列,其中每个数称为矩阵的一个元素。
2. 矩阵元素的表示:A=[a_ij]_{m×n},其中a_ij表示矩阵A的第i行第j列的元素。
3. 矩阵的加法和减法:给定两个相同阶的矩阵A=[a_ij]_{m×n}和B=[b_ij]_{m×n},则它们的和A+B=[a_ij+b_ij]_{m×n}和差A-B=[a_ij-b_ij]_{m×n}定义为对应元素相加或相减得到的结果。
4. 矩阵的数乘:给定一个矩阵A=[a_ij]_{m×n}和一个实数k,则kA=[ka_ij]_{m×n}定义为矩阵A的每个元素乘以实数k得到的结果。
5. 矩阵的乘法:给定一个m×n的矩阵A和一个n×p的矩阵B,它们的乘积AB=[c_ij]_{m×p}定义为矩阵A的第i行与矩阵B的第j列的对应元素乘积之和。
二、矩阵的转置和逆1. 矩阵的转置:给定一个m×n的矩阵A=[a_ij]_{m×n},它的转置记作A^T,其中A^T=[a_ji]_{n×m},即将矩阵A的行变为列,列变为行。
2.矩阵的逆:给定一个n×n的矩阵A,如果存在一个n×n的矩阵B,使得AB=BA=I,其中I是n阶单位矩阵,则称矩阵A是可逆的,矩阵B称为矩阵A的逆矩阵,记作A^{-1}。
三、矩阵的特殊类型1.零矩阵:所有元素都为0的矩阵,记作0。
2.单位矩阵:对角线上的元素都为1,其余元素都为0的矩阵,记作I。
3.对角矩阵:非对角线上的元素都为0的矩阵。
4.上三角矩阵:下三角元素都为0的矩阵。
5.下三角矩阵:上三角元素都为0的矩阵。
6. 对称矩阵:对于任意元素a_ij,有a_ij=a_ji的矩阵,记作A^T=A。
7. 反对称矩阵:对于任意元素a_ij,有a_ij=-a_ji的矩阵,记作A^T=-A。
矩阵的运算的所有公式
矩阵的运算的所有公式矩阵是数学中一个重要的概念,研究矩阵的运算公式对于理解线性代数和计算机图形学等领域都至关重要。
以下是矩阵的运算公式的详细介绍:1.矩阵的加法:对于两个相同大小的矩阵A和B,它们的加法定义为:C=A+B,其中C的元素等于A和B对应元素的和。
2.矩阵的减法:对于两个相同大小的矩阵A和B,它们的减法定义为:C=A-B,其中C的元素等于A和B对应元素的差。
3.矩阵的数乘:对于一个矩阵A和一个标量k,它们的数乘定义为:B=k*A,其中B的元素等于A的对应元素乘以k。
4.矩阵的乘法:对于两个矩阵A和B,它们的乘法定义为:C=A*B,其中C的元素等于A的行向量与B的列向量的内积。
5.矩阵的转置:对于一个矩阵A,它的转置定义为:B=A^T,其中B的行等于A的列,B的列等于A的行,且B的元素和A的对应元素相同。
6.矩阵的逆:对于一个可逆矩阵A,它的逆定义为:A^{-1},使得A*A^{-1}=I,其中I是单位矩阵。
7.矩阵的行列式:对于一个方阵A,它的行列式定义为:,A,是A的元素的代数余子式之和。
8.矩阵的迹:对于一个方阵A,它的迹定义为:tr(A),是A的主对角线上元素之和。
9.矩阵的转置乘法:对于两个矩阵A和B,它们的转置乘法定义为:C=A^T*B,其中C的元素等于A的列向量与B的列向量的内积。
10.矩阵的伴随矩阵:对于一个方阵A,它的伴随矩阵定义为:adj(A),是A的代数余子式构成的矩阵的转置。
11.矩阵的秩:对于一个矩阵A,它的秩定义为:rank(A),是A的线性无关的行或列的最大数量。
12.矩阵的特征值和特征向量:对于一个方阵A,它的特征值是满足方程det(A - λI) = 0的λ值,特征向量是对应于特征值的非零向量。
13.矩阵的奇异值分解(SVD):对于一个矩阵A,它的奇异值分解定义为:A=U*Σ*V^T,其中U和V 是正交矩阵,Σ是一个对角线上元素非负的矩阵。
14.矩阵的广义逆矩阵:对于一个矩阵A,它的广义逆矩阵定义为:A^+,使得A*A^+*A=A,其中A*A^+和A^+*A均为投影矩阵。
11.4 连续梁的整体刚度矩阵
1 i1
(4i1+4i2)∆2 ∆2 2i2∆3
2 i2
2i2∆2
0
1 i1
2 i2
4i2∆3 ∆3
F1 F2 F3
4i1
2i1
0
∆1 ∆2 ∆3
2
=
2i1 4i1+4i2 2i2 0 2i2 4i2
{F}=[K]{∆}
一、单元集成法的力学模型和基本概念
分别考虑每单元对{ 的单独贡献, 分别考虑每单元对{F}的单独贡献,整体刚度矩阵由单元直接集成
{} {}
8
e (3)单刚 [k]e和单元贡献 [K] 中元素的对应关系 单元贡献矩阵是单元刚度矩阵,利用“单元定位向量” 单元贡献矩阵是单元刚度矩阵,利用“单元定位向量”进行 “换码重排位”
1 (1) 1 2 (2)
(1)
单元
[k] = (2)
2 2
1
4i1 2i1 2i1 4i1
2 (1) 3 (2)
9
三、单元集成法的实施 (定位并累加) 定位并累加)
置零, ]=[0]; (1)将[K]置零,得[K]=[0]; 的元素在[ 中按{ 定位并进行累加得[ (2)将[k] 的元素在[K ]中按{λ} 定位并进行累加得[K ]=[K] 的元素在[ 中按{ 定位并进行累加, (3)将[k] 的元素在[K ]中按{λ} 定位并进行累加,得 [K ]=[K] +[K] ; 对所有单元循环一遍, 按此作法对所有单元循环一遍 最后即得整体刚度矩阵[ 按此作法对所有单元循环一遍,最后即得整体刚度矩阵[K]。 3 2 1 3 1 2 1 2 3 2 1 4i 1 0 2i1 1 1 4i 2i 0 1 0 0 0 0 0
1)单元 ①对结点力 的单独贡献 单元 对结点力{F}的单独贡献 1 F1 1 i1 {F}1 = [ F11 F21 F31 ]T
矩阵的基本运算
矩阵的基本运算矩阵是数学中非常重要的一个概念,它在各个领域都有着广泛的应用。
矩阵的基本运算包括矩阵的加法、减法、数乘和矩阵的乘法等。
本文将围绕这些基本运算展开讨论。
首先,我们来讲解矩阵的加法。
如果两个矩阵A和B的维数相同,即都是m行n列的矩阵,那么它们可以相加。
矩阵的加法运算是将对应位置的元素相加得到新的矩阵。
即若A=(a_{ij}),B=(b_{ij}),则A+B=(a_{ij}+b_{ij})。
例如,给定两个矩阵A和B如下:A = [1 2 3][4 5 6]B = [7 8 9][10 11 12]则A与B的和C为:C = [1+7 2+8 3+9][4+10 5+11 6+12]简化运算后,C的结果为:C = [8 10 12][14 16 18]接下来我们讨论矩阵的减法。
矩阵的减法运算与加法类似,也是将对应位置的元素相减得到新的矩阵,即若A=(a_{ij}),B=(b_{ij}),则A-B=(a_{ij}-b_{ij})。
例如,给定两个矩阵A和B如下:A = [1 2 3][4 5 6]B = [7 8 9][10 11 12]则A与B的差D为:D = [1-7 2-8 3-9][4-10 5-11 6-12]简化运算后,D的结果为:D = [-6 -6 -6][-6 -6 -6]矩阵的数乘是指将一个矩阵的每个元素都乘以一个实数。
即若A=(a_{ij})是一个m行n列的矩阵,k是一个实数,那么kA=(ka_{ij})。
例如,给定一个矩阵A和一个实数k如下:A = [1 2 3][4 5 6]k = 2则kA的结果为:kA = [2*1 2*2 2*3][2*4 2*5 2*6]简化运算后,kA的结果为:kA = [2 4 6][8 10 12]最后我们来讨论矩阵的乘法。
矩阵的乘法运算是指矩阵与矩阵之间进行乘法运算,得到一个新的矩阵。
矩阵的乘法有一定的规则,即若A是一个m行n列的矩阵,B是一个n行p列的矩阵,那么它们可以相乘,得到一个m行p列的矩阵C。
矩阵及其运算
矩阵及其运算矩阵的线性运算(1)加法:两同型矩阵()ij m nA a ´=与()ij m nB b ´=的和矩阵为()ij ij m nA B a b ´+=+.(2)数乘法:数k 与矩阵()ij m nA a ´=的数量乘积矩阵()ij m nkA ka ´=.(3)m s ´矩阵()ij m s C c ´=称为矩阵()ik m n A a ´=与()kj n s B b ´=的乘积. 其中1122ij i j i j in nj c a b a b a b =+++ (1,2,,;1,2,,i m j s ==L L ).(4)k kA A A A =6447448L 个为n 阶方阵A 的k 次幂,特别规定0A E =.(5)1110()m m m m f A a A a Aa A a E --=++++L (i a 为数)为方阵A 的多项式.矩阵的转置以()ij m nA a ´=的行为列,列为行构成的n m ´矩阵()T ji n m A a ´=为A 的转置矩阵.A 是n 阶方阵,如果TA A =,称A 为对称矩阵;如果TA A =-,称A 为反对称矩阵.方阵的行列式以n 阶方阵A 的元素构成的行列式ij na 称为方阵A 的行列式.记为A 或det A .矩阵运算注意事项:利用运算定义和运算律进行运算.注意(ⅰ)第一个矩阵的列数与第二个矩阵的行数相等时两矩阵乘积才有意义.(ⅱ)由于乘法没有交换律,在进行两个矩阵乘积时,矩阵因子的顺序不能变. (ⅲ)矩阵的乘法不满足消去律. 即AB AC =且A O ¹,不一定有B C =;BA CA =且A O ¹,不一定有B C =.特别地,,A B O = 且A O ≠,不一定有B O =.(ⅳ)我们在做多个矩阵乘积时经常使用乘法结合律()()A BC AB C =. (ⅴ),A B 分别是,m n ns 创矩阵,则()TTTAB B A =.(ⅵ)只有方阵才定义行列式;矩阵是数表,行列式是数值,这是它们之间的本章区别.可逆矩阵1.设A 是n 阶方阵,如果存在n 阶方阵B ,使得AB BA E ==,称A 为可逆矩阵,称B 为A 的一个逆矩阵.2.可逆矩阵的逆矩阵唯一.3. 设()ij n nA a ⨯=,称由以A 的第(1,2,,i i n = 行元素在A 中的代数余子式(1,2,,)ij A j n = 为第i 列元素构成的矩阵*A ()jin nA ⨯=为A 的伴随矩阵.4.设A 是n 阶方阵,*A 是A 的伴随矩阵,则**AA A A A E ==. 5. n 阶方阵A 是可逆的充分必要条件为0A ≠.而且1*1A A A-=(*1A A A-=).6.可逆矩阵具有如下运算性质:(ⅰ)A 是n 阶可逆矩阵,A 的逆矩阵1A -也可逆,且()11A A --=; (ⅱ)A 是n 阶可逆矩阵,k 是非零数,则kA 可逆,且()111kA kA---=;(ⅲ),A B 都是n 阶可逆矩阵,那么A B 也可逆,且111()AB B A ---=;(ⅳ)A 是n 阶可逆矩阵,T A 也可逆,且()()11TTA A--=;(ⅴ)A 是m 阶可逆矩阵,,B C 都是m n ⨯矩阵,且AB AC =,则B C =,A 是n 阶可逆矩阵,,BC 都是m n ⨯矩阵,且BA CA =,则B C =.常用解题方法:(设A 是n 阶方阵)1.求可逆矩阵的逆矩阵:1*1A A A-=(用于阶数较低的具体给定的数字矩阵求逆或理论证明)2.利用定义证明矩阵可逆,或求满足给定方程的矩阵A 的逆矩阵:找到n 阶方阵B ,使得AB BA E ==,则A 可逆,且1A B -=.注意 *A 的第(1,2,,)i i n = 列元素是A 的第(1,2,,)i i n = 行元素在A 的代数余子式;*A 的第(1,2,,)i i n = 行元素是A 的第(1,2,,)i i n = 列元素在A 的代数余子式.分块矩阵及其运算用若干条横线和纵线将矩阵分成若干小矩阵,以小矩阵为元素的矩阵表示形式称为分块矩阵.我们将这些小块称为矩阵的子块. 1 加法 对两个m n ´矩阵()ij m nA a ´=,()ij m nB b ´=进行同样分块,则A B +为对应块相加得到的分块矩阵; 2 数乘法 设()ij m nA a ´=是一个m n ´矩阵,k 是一个数,将kA 为由k 数乘每个子块矩阵得到的分块矩阵;3乘法 设()ik m n A a ´=, ()kj n t B b ´=,将,A B 分块为111212122212.....................p p s s sp A A A A A A A A A A 骣÷ç÷ç÷ç÷ç÷ç÷=ç÷ç÷÷ç÷ç÷ç÷ç桫,111212122212.....................q q p p pq B B B B B B B B B B 骣÷ç÷ç÷ç÷ç÷ç÷=ç÷ç÷÷ç÷ç÷ç÷ç桫,则1121212212.....................qq s s sq C C CC C C A B C C C 骣÷ç÷ç÷ç÷ç÷ç÷=ç÷ç÷÷ç÷ç÷ç÷ç桫.其中ik A 为i k m n ´矩阵,kj B 为k j n t ´矩阵,1122i j i j ip p j j i B A B C A A B =+++L . 4 转置 设()ij m nA a ´=是一个m n ´矩阵,将A 分块为111212122212.....................t t s s st A A A A A A A A A A 骣÷ç÷ç÷ç÷ç÷ç÷=ç÷ç÷÷ç÷ç÷ç÷ç桫,则112111222212.....................T T Ts T TT T s T TT t tst A A A A A A A A A A 骣÷ç÷ç÷ç÷ç÷ç÷=ç÷ç÷÷ç÷ç÷ç÷ç桫. 注意事项:1.利用分块矩阵表示矩阵或进行矩阵运算只是为了表达简便.分块矩阵的运算与普通数字元素的运算法则和运算律是类似的;2.第一个矩阵列的分块方式与第二个矩阵行的分块方式必须相同,即ik A 列数必须等于kj B 的行数,这时两分块矩阵的乘积才有意义;3.由于矩阵乘法没有交换律,作分块矩阵乘法时,一定要注意子块的前后顺序不能换.即上面的ik kj A B 绝对不能写成kj ik B A .4.分块矩阵的转置不仅要将子块为元素构成的矩阵看成普通矩阵进行转置,还要将每块转置.初等变换与初等矩阵初等变换1.对调两行(或列);2.以数(0)k k ≠乘矩阵某一行(或列)的所有元素; 2.把矩阵的某行(或列)所有元素乘一个数加到另一行(或列)对应位置的元素上. 矩阵的等价标准形1.称具有如下特点m n ⨯矩阵A 为行阶梯形矩阵:(ⅰ)A 的前()r r m ≤行,每行元素均不全为0,后m r -行元素都为零; (ⅱ)第(1)k k r ≤≤行的第一个非零元素为kkj a ,且满足12r j j j <<< .如果行阶梯形矩阵还满足:(ⅲ)第(1)k k r ≤≤行的第一个非零元素1kkj a =,且(1,2,,)k kj a k r = 所在的k j 列的其它元素都为0,就称A 为行最简形矩阵.2.任何矩阵都可以经过初等行变换化为行阶梯形矩阵,进而可以化为行最简形矩阵.3.任何一个m n ⨯矩阵A 都可以通过初等变换化为()()()()rr n r m r rm r n r E O OO ⨯--⨯-⨯-⎛⎫⎪ ⎪⎝⎭型矩阵. 初等矩阵1.由单位矩阵E 经过一次初等变换得到的矩阵称为初等矩阵.2.设A 一个m n ⨯矩阵,A 左乘一个m 阶初等矩阵相当于对A 作一次相应的初等行变换,A 右乘一个n 阶初等矩阵相当于对A 作一次相应的初等列变换.3.A 与B 等价当且仅当存在可逆矩阵P 与可逆矩阵Q ,使得PBQ A =.4.n 阶方阵A 可逆当且仅当A 可以写成一些初等矩阵的乘积.5.n 阶方阵A 可逆当且仅当A 可以只用初等行变换化为单位矩阵.常用解题方法1.用初等行变换求可逆矩阵A 的逆矩阵的求法:()()1AE EA-−−−−→初等行变换.2.用初等行变换求矩阵方程A X B =(A 可逆)的求法:()()1A B EA B -−−−−→初等行变换则1X AB -=即可求得.矩阵基础部分习题1. A 是行等和矩阵(各行元素之和都相等),且A 可逆,证明:1A-也是行等和矩阵.证明:设A 的行和为a ,则11,(0)11A a a 骣骣鼢珑鼢珑鼢珑鼢= 珑鼢珑鼢珑鼢鼢珑桫桫M M ,所以111111A a --骣骣鼢珑鼢珑鼢珑鼢=珑鼢珑鼢珑鼢鼢珑桫桫M M , 即1A-是行等和矩阵.2.设实矩阵111213212223313233a a a A a a a a a a ⎛⎫⎪= ⎪ ⎪⎝⎭,且TA A O =,证明A O =.3.已知()33ij A a ´=可逆,将A 的第2列加上第3列的5倍,然后第1列减去第2列的2倍得到B , 求1B A -解:11121511B A 骣骣鼢珑鼢珑鼢珑鼢=-珑鼢珑鼢珑鼢鼢珑桫桫, 111111211151BA ----骣骣鼢珑鼢珑鼢珑鼢=-珑鼢珑鼢珑鼢鼢珑桫桫,11111211151B A ---骣骣鼢珑鼢珑鼢珑鼢=-珑鼢珑鼢珑鼢鼢珑桫桫1112112115151骣骣骣鼢 珑 鼢 珑 鼢 珑 鼢 ==珑 鼢 珑 鼢 珑 鼢鼢珑 --桫桫桫. 4.设A 为3阶矩阵,将A 的第2行加到第1行得B ,再将B 的第1列的-1倍加到第2列得C ,记110010001P ⎛⎫⎪= ⎪ ⎪⎝⎭,则 (A )1.C P AP -= (B )1.C P A P -=(C ).TC P AP =(D ).TC PAP = 【 】5. ,A B 都为n 阶方阵,且A B AB += 1) 证明:A E -可逆2)证明:AB BA =3)如果130210002B 骣-÷ç÷ç÷ç÷=ç÷ç÷ç÷÷ç桫,求A 1)证明:由A B AB +=有()A A E B O --=,所以()A E A E B E ---=-,即()()A E B E E --=,所以A E -可逆,且1()A E B E --=-2)由()()A E B E E --=,有()()B E A E E --=, 所以()()()()A E B E B E A E --=--,即有AB BA =3)由1()A E B E --=-有()11100030010200001001A EB E --骣骣-鼢珑鼢珑鼢珑鼢=+-=+珑鼢珑鼢珑鼢鼢珑桫桫1100102210011010001033001001002骣骣鼢珑鼢珑鼢珑鼢珑骣鼢珑÷鼢ç珑÷鼢ç÷珑鼢ç÷鼢=+-=-珑ç÷鼢珑ç÷鼢珑ç÷鼢珑÷ç鼢桫珑鼢珑鼢珑鼢鼢珑鼢珑桫桫6.已知121210121,()(),00210003A B A E A E -骣--÷ç÷ç÷ç÷-ç÷ç÷==-+ç÷ç÷÷ç÷ç÷ç÷ç桫求1()B E -- 解:由1()(),B A E A E -=-+有()(),A E B A E -=+,AB B A E --=()()2A B E B E E ---=,所以()()11()2,()2A EB E E B E A E ---=-=-7. 方阵A 满足223A A E O +-= 1)求证:4A E +可逆,并求其逆 2)讨论:()A nE n Z ++ 的可逆性1)证明:由于223A A E O +-=,有()()425A E A E E +-=-,所以4A E +可逆,其逆为()()11425A E A E -+=--.2)由于223A A E O +-=,有()()()()222323(3)(1)A nE A n E n n E n n E n n E 轾轾+--=--+=---=--+臌臌,所以3,1n n 构-时,()A nE n Z ++ 可逆,其逆为()()()11223A nE A n E n n -轾+=--臌--+当3n =时,如果A E O -=,则34A E E +=可逆;如果A E O - ,由()()3A E A E O +-=,所以()3A E X O +=有非零解, 所以()3r A E n +<, 3A E +不可逆.当1n =-时,如果3A E O +=,则4A E E -=-可逆; 如果3A E O + ,与上面相同有()r A E n -<,A E -不可逆. 8.,,A B A B +均为n 阶可逆矩阵,求证11A B--+也可逆,并求其逆证明:()11,A A BB B A --+=+所以()()1111ABAB A B ----+=+,所以11AB--+可逆,且()()1111A B B B A A ----+=+.9.设A 为n 阶非零矩阵,E 为n 阶单位矩阵,30A =,则( ) A .E A -不可逆,E A +不可逆;B. E A -不可逆,E A +可逆; C. E A -可逆,E A +可逆; D. E A -可逆,E A +不可逆10.设16,A X A A X A -=+ 其中100310041007A ⎛⎫ ⎪⎪⎪= ⎪ ⎪ ⎪ ⎪⎝⎭,求X . 11.设,A B 均为2阶方阵,**,A B 分别为,A B 的伴随矩阵,若3,2A B ==则O A BO ⎛⎫⎪⎝⎭的伴随矩阵为( ) A .**32OA B O ⎛⎫⎪ ⎪⎝⎭;B. **23O A B O ⎛⎫⎪ ⎪⎝⎭;C. **23OB A O ⎛⎫⎪ ⎪⎝⎭;D. **32OB A O ⎛⎫⎪ ⎪⎝⎭12.A 是3阶方阵,B 是2阶方阵,且2A =-,1B =,则23A O OB=- ;*2A= .13.33A R ´Î,且()**16,det 0,AA =>求()det 2A -解:()()2*22**16,det 2,AA A A ====()3det 2(2)det 16A A -=-=-14.设A 是n 阶方阵,3A =,*A 是A 的伴随矩阵,则1*2AA--=15.已知A ,B 都是3阶方阵,且9A =-,3AB E O +=,求B 及12AO OB -⎛⎫⎪⎝⎭.16.设矩阵2112A ⎛⎫=⎪-⎝⎭,E 为2阶单位矩阵,矩阵B 满足2B A B E =+,则B =. 17. ()ij n A a =可逆,且1,A A-所以元素都是整数,证明:11A =-或解:由于行列式是取自不同行、不同列元素乘积的代数和,且1,A A -所以元素都是整数,所以1,A A-都是整数,又111,AAA A--==所以11A =-或18.如果A B B A =,则称矩阵B 与A 可交换,求与矩阵A 可交换的矩阵具有的形式.100002000010011A ⎛⎫ ⎪⎪= ⎪ ⎪⎝⎭; 矩阵的秩1.设A 是一个m n ⨯矩阵,如果A 中存在r 阶子式不为零,而所有1r +阶子式(如果有的话)全为零,我们称r 为矩阵A 的秩,记为()R A 或秩()A .2. 矩阵的秩具有如下性质: (ⅰ)()0R A =当且仅当A O =; (ⅱ)()()T R A R A =;(ⅲ)()()R kA R A =,其中k 为非零数;(ⅳ)n 阶方阵A 的秩()R A n =的充分必要条件0A ≠; (ⅴ)n 阶方阵A 可逆的充分必要条件为()R A n =. 3. 行阶梯形矩阵A 的非零行的行数等于A 的秩. 4.初等变换不改变矩阵的秩.5.矩阵P ,Q 可逆,则()()R PAQ R A =.6. 设A 是秩为r 的m n ⨯矩阵,则存在m 阶可逆矩阵P 和n 阶可逆矩阵Q ,使得rE O PAQ OO ⎛⎫=⎪⎝⎭. 常用方法1.求矩阵A 的秩:利用矩阵的初等变换将矩阵A 化为阶梯形矩阵,阶梯数即为矩阵A 的秩.2.如果A 是n 阶方阵,0A ≠时()R A n =.求元素含有参数的方阵A 的秩时,先求出0A ≠时的参数取值,此时()R A n =; 对于使0A =的参数再特别讨论.1.已知111111A λλλ⎛⎫⎪= ⎪ ⎪⎝⎭,讨论λ为何值时(1)()1R A =;(2) ()2R A =;(3)()3R A =. 2.设A 是m n ⨯矩阵,证明: ()1R A =当且仅当存在m 维列向量矩阵α和n 维行向量矩阵Tβ,使得TA αβ=.3.设A 是秩为1的3阶方阵,证明:存在不全为零的数1,23,a a a 和不全为零的数1,23,b b b ,使得111213212223313233a b a b a b A a b a b a b a b a b a b ⎛⎫⎪= ⎪ ⎪⎝⎭;并求100A 4.(2010)2. 设,AB 分别为m n ⨯和n m ⨯矩阵,E 为单位矩阵,且AB E =,则 A .秩A m =,秩B m =;B. 秩A m =,秩B n = ; C. 秩A n =,秩B m = ;D. 秩A n =,秩B n =5. ,A B 都为n 阶实方阵,G 为n m ´矩阵,()r G n =,证明:如果AG BG = 则A B =证明:方法一:()r G n =,则G 有n 个线性无关的列向量,可以通过列交换使的G 的前n 列线性无关,即存在可逆矩阵Q (是换位矩阵的乘积)使得12(,)GQ G G =, 其中1G 为n 阶可逆方阵,由AG BG =,有AG Q BG Q =,所以1212(,)(,)A G G B G G = 因此1212(,)(,)AG AG BG BG =,所以11AG BG =,即有A B =方法二:由AG BG =,有()0A B G -=,所以G 的列向量都是齐次线性方程组()0A B x -=的解,所以()0A B x -=有n 个线性无关的解向量,所以()0r A B -=,即A B =.方法三:利用()()Tr A R AA =由于AG BG =,则TTA G GB G G =,而()()Tr G G r G n ==,所以TG G 是n 阶可逆方阵,所以A B =.6.A 是m n ´实矩阵,证明:()()()T T r A r A A r AA == 证明:先证明齐次线性方程组0Ax =与0TA Ax =同解如果0x 是0Ax =的解,则00Ax =,所以00T A A x =,所以0x 是0TA Ax =的解反之,0x 是0TA Ax =,则00T A A x =,有000T Tx A A x =,即()000TA x A x =,所以00Ax =即0x 是0Ax =的解,所以0Ax =与0TA Ax =同解.所以()(),T n r A n r A A -=- 即()()T r A r A A =线性方程组的解部分知识概要1.线性方程组Ax b =(x 是n 维列向量)的系数矩阵为A ,增广矩阵为()A A b =,则:(ⅰ)线性方程组无解的充分必要条件为()()R A R A <; (ⅱ)线性方程组有唯一解的充分必要条件为()()R A R A n ==; (ⅲ)线性方程组有无穷多解的充分必要条件为()()R A R A n =<. 2.齐次线性方程组0Ax =(x 是n 维列向量)永远有零解. (ⅰ)0Ax =只有零解的充分必要条件为()()R A R A n ==; (ⅱ)0Ax =有非零解的充分必要条件为()()R A R A n =<.3.矩阵方程A X B =(X 是n s ⨯矩阵)的系数矩阵为A ,增广矩阵为()A A B =,则关于方程A X B =的解有与1中相同的结论.1.讨论,,,a b t k 取何值时下列方程组有解?并在有解时求解1)123423423412340221(3)2321x x x x x x x x a x x b x x x ax ì+++=ïïïï++=ïíï-+--=ïïï+++=-ïïî;解:11110111101011101221012210122101320132001013211012310010a b a b a b aa a 骣骣骣---鼢珑 鼢 珑 鼢 珑 鼢 珑 鼢 珑 鼢 珑 鼢 珑 鼢 -------+鼢 珑 鼢 珑 鼢 珑 鼢珑 ------桫桫桫1,1a b =?时,方程无解1,1a b ==-时,方程有无穷多解,此时,10111101110122101221001010000000100a b a 骣骣------鼢珑鼢珑鼢珑鼢珑鼢珑鼢®珑鼢珑鼢-+鼢珑鼢珑鼢珑鼢珑-桫桫,1342341122x x x x x x ì=-++ïïíï=--ïî 方程组的通解为121234111122(010001x x c c x x 骣骣骣骣-÷鼢 ç珑 ÷鼢 ç珑 ÷鼢 ç珑 ÷鼢 --ç珑 ÷鼢 ç珑 ÷鼢 =++ç珑 ÷鼢 ç÷珑 鼢 ÷鼢 ç珑 ÷鼢 ç珑 ÷鼢 ç珑 鼢÷珑 ç桫桫桫桫12,c c 为任意数) 1a ¹时,方程有唯一解,此时110001101111101110122110122101001211001010010110010001010010001b a b a b a b b a a a 骣+÷ç-+÷ç骣 ---ç-÷÷çç骣---÷÷çç÷÷ç÷çç÷÷ç÷+ç÷ç÷÷ç÷琪-ç÷ç÷÷çç÷ç÷÷ -ç÷ç+ç÷÷÷ççç÷÷-+÷çç÷÷ç÷ç?÷÷ç÷-çç÷÷÷ç÷çç÷ç-÷桫ç÷ç-桫ççç桫÷÷÷÷÷ 所以12341111121110b x a b x a b x a x ì+ïï=-+ïï-ïïï+ï=-ï-íïï+ïï=ï-ïïï=ïïî;2)1231231234324ax x x x bx x x bx x ì++=ïïïï++=íïï++=ïïî解:1141101110211311310121214001001a a b a b b bbb骣骣骣---鼢珑 鼢 珑 鼢 珑 鼢 珑 鼢 珑 鼢 珑 鼢鼢珑 桫桫桫()()()00112210114210121012001101b a a a a bb¹骣÷ç÷ç骣÷-------ç÷÷çç÷÷çç÷÷çç÷÷ ç÷ç÷ç÷÷çç÷÷ç÷÷çç桫÷ç÷ç÷ç桫()1001421012101a ab b 骣÷ç----÷ç÷ç÷ç÷ç÷ç÷®ç÷÷ç÷ç÷ç÷ç÷ç÷ç桫()142100*********a b ab a b b¹骣--÷ç÷ç÷ç--÷ç÷ç÷ç÷ç÷®÷ç÷ç÷ç÷ç÷ç÷ç÷÷ç桫()()42100111210011010b ab a b b a b b骣--÷ç÷ç÷ç--÷ç÷ç÷ç÷ç÷-÷ç ç÷ç÷--ç÷ç÷ç÷ç÷÷ç÷ç÷ç÷桫,所以1,0a b 构时方程组有唯一解()()12342111211b ab a b x b x a b x b骣--÷ç÷ç÷ç--÷ç÷骣ç÷ç÷÷ç÷ç÷ç-÷÷çç÷= çç÷÷çç÷÷--çç÷÷ç÷ç÷桫ç÷ç÷÷ç÷ç÷ç÷桫 当0b =时,方程组无解,当0,1b a?时有()11001420002114113101210121214110101a a ab b b bbb 骣骣鼢珑-----鼢珑骣鼢珑鼢÷ç珑鼢÷ç珑鼢÷ç珑鼢÷ ç珑鼢÷琪珑÷鼢ç珑÷鼢÷ç珑桫鼢珑鼢珑鼢珑桫桫所以11,0,2a b b=构时无解;11,2a b ==时,114000011310121214012ab b骣骣鼢珑鼢珑鼢珑鼢®珑鼢珑鼢珑鼢鼢珑桫桫,无穷多解.解为123212001x x c x 骣骣骣-鼢珑 鼢 珑 鼢 珑 鼢 =+珑 鼢 珑 鼢 珑 鼢 鼢 珑 桫桫桫(c 为任意数)3).设2221212aaa A aa ⎛⎫ ⎪ ⎪=⎪ ⎪ ⎪⎝⎭,且A 满足A X B =,123(,,),(1,0,0)T T x x x x B == (I)求证(1)nA n a =+;(II )a 为何值时,方程组有唯一解,并求1x (III )a 为何值时,方程组有无穷多解,并求通解4).1101011A λλλ⎛⎫⎪=- ⎪ ⎪⎝⎭,11a b ⎛⎫ ⎪= ⎪ ⎪⎝⎭,AX b =存在两个不同解, (1)求,a λ;2)求AX b =的通解。
矩阵及其运算PPT课件
设
1 A 1
1 1
1 1 2 3 1, B 1 2 4
1 1 1 0 5 1
求3AB 2A及 AT B
2 13 22
0 5 8
答案:3AB 2A 2 17 20 , AT B 0 5 6.
第22页/共24页
六、方阵的行列式
2010年期末考题(I)
二、选择(每题4分,共16分)
1、设A与B均为n阶方阵,则下列结论中成立的是( B )
A. |AB|=0,则A=0或B=0; B. |AB|=0,则|A|=0或|B|=0; C. AB=0,则A=0或B=0; D. AB≠0,则|A|≠0或|B|≠0;
T ,
则An ____1___12.
1 3
23
3
n
1
2
1
2 3
矩阵拆分相乘
3
3
1
2
第13页/共24页
2012年期末考试题
二项式法
1
4、设A
0
0 0
2012年期末考试题
0
1
,
则A
n
n nn1
_0____n .
0
0
n(n 1) n2
2
nn1
n
五.(10分)(线性代数I,36学时专业学时做 )设
转置矩阵的运算性质 (1) (AT)T = A; (2) (A+B)T = AT + BT;
(3) (A)T = AT;
(4) (AB)T = BTAT;
第2页/共24页
由n 阶方阵A 的元素所构成的行列式叫做方阵A 的行列式, 记作 | A | 或 detA .
矩阵的简单运算公式
矩阵的简单运算公式矩阵是线性代数中一个非常重要的概念,广泛应用于各个领域,如数学、物理、工程等。
矩阵的运算是对矩阵进行各种操作的过程,包括加法、减法、乘法等。
本文将介绍矩阵的简单运算公式,并给出相应的例子,以帮助读者更好地理解矩阵运算的基本原理。
一、矩阵的加法矩阵的加法是指将两个矩阵的对应元素相加,依次得到一个新的矩阵。
加法的具体操作如下:设A和B为两个相同阶数的矩阵,即A和B的行数和列数相等。
则它们的和记作C=A+B,C的每个元素ci,j等于A和B相应元素的和,即ci,j = ai,j + bi,j。
举个例子,假设有两个矩阵:A = [1 2 3][4 5 6]B = [7 8 9][10 11 12]则A和B的和矩阵C为:C = A + B = [1+7 2+8 3+9][4+10 5+11 6+12]二、矩阵的减法矩阵的减法是指将两个矩阵的对应元素相减,得到一个新的矩阵。
减法的操作与加法类似,不同之处在于相减而非相加。
设A和B为两个相同阶数的矩阵,即A和B的行数和列数相等。
则它们的差记作D=A-B,D的每个元素di,j等于A和B相应元素的差,即di,j = ai,j - bi,j。
继续以上面的矩阵A和B为例,它们的差矩阵D为:D = A - B = [1-7 2-8 3-9][4-10 5-11 6-12]三、矩阵的数乘矩阵的数乘是指将一个矩阵的每个元素都乘以一个常数,得到一个新的矩阵。
数乘的具体操作如下:设A为一个矩阵,k为一个常数。
则A乘以k的结果记作E=kA,E 的每个元素ei,j等于A相应元素的k倍,即ei,j = k * ai,j。
继续以上面的矩阵A为例,假设k=2,则矩阵A乘以2的结果E为:E = 2A = [2*1 2*2 2*3][2*4 2*5 2*6]四、矩阵的乘法矩阵的乘法是指将一个矩阵与另一个矩阵相乘,得到一个新的矩阵。
乘法的操作稍微复杂一些,需要满足一定的条件。
设A是一个m×n的矩阵,B是一个n×p的矩阵,则AB的结果是一个m×p的矩阵。
矩阵的运算及其运算规则
矩阵的运算及其运算规则在数学和许多相关领域中,矩阵是一种非常重要的工具。
它不仅在理论研究中发挥着关键作用,还在实际应用中有着广泛的用途,比如图像处理、数据分析、物理学等。
要深入理解和运用矩阵,就必须掌握矩阵的运算及其运算规则。
首先,让我们来了解一下矩阵的加法。
矩阵加法是指两个具有相同行数和列数的矩阵,将它们对应位置的元素相加所得到的新矩阵。
比如说,有矩阵 A = 1 2; 3 4 和矩阵 B = 5 6; 7 8,那么 A + B = 1 + 52 + 6;3 + 74 + 8 = 6 8; 10 12。
这里要注意的是,只有当两个矩阵的行数和列数都完全相同,才能进行加法运算。
接下来是矩阵的减法,它与加法类似,只是将对应位置的元素相减。
例如,矩阵 A = 1 2; 3 4,矩阵 B = 05 1; 15 2,那么 A B = 1 05 2 1; 3 15 4 2 = 05 1; 15 2 。
然后是矩阵的数乘运算。
这是指用一个数乘以矩阵中的每一个元素。
假设矩阵 A = 1 2; 3 4,如果用 2 去乘以矩阵 A,那么得到的新矩阵就是 2A = 2×1 2×2; 2×3 2×4 = 2 4; 6 8 。
再来说说矩阵的乘法。
矩阵乘法相对来说要复杂一些。
不是任意两个矩阵都能相乘的,只有当第一个矩阵的列数等于第二个矩阵的行数时,它们才能相乘。
假设矩阵 A 是 m×n 的矩阵,矩阵 B 是 n×p 的矩阵,那么它们相乘得到的矩阵 C 是 m×p 的矩阵。
具体的计算方法是,矩阵C 中第 i 行第 j 列的元素等于矩阵 A 的第 i 行元素与矩阵 B 的第 j 列对应元素相乘的和。
例如,矩阵 A = 1 2; 3 4,矩阵 B = 5 6; 7 8,那么AB = 1×5 + 2×7 1×6 + 2×8; 3×5 + 4×7 3×6 + 4×8 = 19 22; 43 50 。
矩阵运算法则PPT课件
是 A 的逆矩阵,
利用待定系数法
则
AB 2 1 a b 1 0
1 0 c d 0 1
2a c 2b d 1 0 a b 0 1
第33页/共78页
2a c 1,
2b
d a
0, 0,
b 1,
又因为 AB
a 0,
b 1,
c
1,
d 2.
内容提要
• 矩阵的下列运算的性质与应用 • 乘法 • 转置 • 初等变换 •逆
第1页/共78页
乘法
定义
设矩阵
A
aij
,B
mn
bij
,那么
sn
矩阵A与矩阵B的乘积是一个m n矩阵 C s
cij mn ,其中cij ai1b1 j ai2b2 j aisbsj= aikbkj k1
1设
A=
aaa123,,,
1 1 1
a1, 2 a2, 2 a3, 2
a1, a2, a3,
333
计算并总结规律。
(1)
1 0 0
0 1 0
001 A
(2)
A
1 0 0
0 1 0
001
第13页/共78页
(3)
1 0 0
0 0 1
010
A
(4)
A
1 0 0
0 0 1
010
(5)
1 0 0
0 k 0
a1, a2, a3,
222
第18页/共78页
初等矩阵的概念
定义 由单位 E矩阵经过一次初等变换得到 的方阵称为初等矩阵.
三种初等变换对应着三种初等方阵. 1. 对调两行或两列; 2.以数 k 0 乘某行或某列; 3.以数 k 乘某行(列)加到另一行(列)上去.
矩阵的各种运算详解
一、矩阵的线性运算定义1 设有两个矩阵和,矩阵与的和记作, 规定为注:只有两个矩阵是同型矩阵时,才能进行矩阵的加法运算. 两个同型矩阵的和,即为两个矩阵对应位置元素相加得到的矩阵.设矩阵记,称为矩阵的负矩阵, 显然有.由此规定矩阵的减法为.定义2 数与矩阵A的乘积记作或, 规定为数与矩阵的乘积运算称为数乘运算.矩阵的加法与矩阵的数乘两种运算统称为矩阵的线性运算. 它满足下列运算规律:设都是同型矩阵,是常数,则(1)(2) ;(3)(4)(5)(6)(7)(8)注:在数学中,把满足上述八条规律的运算称为线性运算.二、矩阵的相乘定义3设矩阵与矩阵的乘积记作, 规定为其中,(记号常读作左乘或右乘.注: 只有当左边矩阵的列数等于右边矩阵的行数时, 两个矩阵才能进行乘法运算.若,则矩阵的元素即为矩阵的第行元素与矩阵的第列对应元素乘积的和. 即.矩阵的乘法满足下列运算规律(假定运算都是可行的):(1)(2)(3)(4)注: 矩阵的乘法一般不满足交换律, 即例如, 设则而于是且从上例还可看出: 两个非零矩阵相乘, 可能是零矩阵, 故不能从必然推出或此外, 矩阵乘法一般也不满足消去律,即不能从必然推出例如, 设则但定义4如果两矩阵相乘, 有则称矩阵A与矩阵B可交换.简称A与B可换.注:对于单位矩阵, 容易证明或简写成可见单位矩阵在矩阵的乘法中的作用类似于数1.更进一步我们有命题1设是一个n阶矩阵,则是一个数量矩阵的充分必要条件是与任何n阶矩阵可换。
命题2设均为n阶矩阵,则下列命题等价:(1)(2)(3)(4)三、线性方程组的矩阵表示设有线性方程组若记则利用矩阵的乘法, 线性方程组(1)可表示为矩阵形式:(2)其中矩阵称为线性方程组(1)的系数矩阵. 方程(2)又称为矩阵方程.如果是方程组(1)的解, 记列矩阵则,这时也称是矩阵方程(2)的解; 反之, 如果列矩阵是矩阵方程(2)的解, 即有矩阵等式成立, 则即也是线性方程组(1)的解. 这样, 对线性方程组(1)的讨论便等价于对矩阵方程(2)的讨论. 特别地, 齐次线性方程组可以表示为将线性方程组写成矩阵方程的形式,不仅书写方便,而且可以把线性方程组的理论与矩阵理论联系起来,这给线性方程组的讨论带来很大的便利.四、矩阵的转置定义6把矩阵的行换成同序数的列得到的新矩阵, 称为的转置矩阵, 记作(或). 即若则.矩阵的转置满足以下运算规律(假设运算都是可行的):(1)(2)(3)(4)五、方阵的幂定义5设方阵, 规定称为的次幂.方阵的幂满足以下运算规律(假设运算都是可行的):(1)(2)注: 一般地,为自然数命题3 设均为n阶矩阵,则有为自然数,反之不成立。
矩阵的基本运算和应用
矩阵的乘法
两个矩阵相乘,需要满足第一个矩阵的列数等于第二个矩阵的行数,结
果矩阵的行数等于第一个矩阵的行数,列数等于第二个矩阵的列数。乘
法运算遵循特定的运算法则。
特殊类型矩阵
方阵
01 行数和列数相等的矩阵称为方
阵。
零矩阵
02 所有元素都为零的矩阵称为零
矩阵。
对角矩阵
03 除主对角线外,其他元素都为
零的方阵称为对角矩阵。
矩阵乘法运算
乘法定义
设A为m×n矩阵,B为n×s矩阵,那么称m×s矩阵C为矩阵A 与B的乘积,记作C=AB。
运算步骤
矩阵乘法运算时,先将第一个矩阵的每一行分别与第二个 矩阵的每一列相乘,再将得到的积相加,得到结果矩阵的 对应元素。
运算性质
矩阵乘法一般不满足交换律,但满足结合律和分配律,且 单位矩阵作为乘法的单位元。
特征选择
基于矩阵分解等方法,选取对模型训练有重要贡 献的特征。
主成分分析(PCA)原理及实现
主成分分析(PCA)原理及实现
计算协方差矩阵。
对原始数据进行标准化处 理。
实现步骤
01
03 02
主成分分析(PCA)原理及实现
01
对协方差矩阵进行特征值分解,得到特征值和特征 向量。
02
选择前k个最大特征值对应的特征向量组成矩阵W。
• 克拉默法则:如果线性方程组的系数矩阵A的行列式|A|不等 于零,则该线性方程组有唯一解,且解可以通过系数矩阵A 和常数项向量b的行列式计算得到。
克拉默法则求解线性方程组
具体步骤 构造系数矩阵A和常数项向量b。 计算系数矩阵A的行列式|A|。
克拉默法则求解线性方程组
对于每一个未知数,将系数矩阵A中对应列替换为常数项向量b,得到新的矩阵B,并计算其行列式|B| 。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
7
在后一种情况下, 故正交矩阵的定义可改为: 在后一种情况下, * = AT 故正交矩阵的定义可改为: A 矩阵A 矩阵A称为正交矩阵 ⇔ AT = A−1
⇔ AAT = E
其中最后的定义是数学中最常用的定义. 其中最后的定义是数学中最常用的定义. cos α − sin α 例如 R = sin α cos α cos α sin α −1 满足 R' = =R − sin α cosα 故R是正交矩阵. 是正交矩阵. 对应于U矩阵和正交矩阵的线性变换是所谓的酉 对应于U 变换和正交变换. 变换和正交变换.
5
对于矩阵的转置运算,下列性质成立: 对于矩阵的转置运算,下列性质成立:
(1)( AT )T = A; ( 3)( kA)T = kAT ;
( 2)( A + B )T = AT + BT ;
特别注意顺序! (4)( AB) = B A 特别注意顺序!
T T T
共轭转置矩阵的概念与性质 定义 A是复矩阵,若把其所有元素都改成共轭 是复矩阵, 复数,得到矩阵A 复数,得到矩阵A的共轭矩阵 A ; 的共轭转置矩阵. A 的转置矩阵 A 叫A的共轭转置矩阵. A* 即 A* = AT 注意区分伴随矩阵! 记作 注意区分伴随矩阵! 2 2 − 3i 例如 3i − 3i 2i 5 A = − 2i 0 , A = 2i 0 , A* = 5 4−i 5 4+ i 2 0 4+i 6
§11.4 矩阵的其它运算 11.
1.矩阵的线性运算 2.转置矩阵与共轭转置矩阵 3.U矩阵与H矩阵 .U矩阵与H 矩阵与 4.对角矩阵与准对角矩阵
1
一、矩阵的线性运算
定义 设A和B是两个同型矩阵,则把两个矩阵 是两个同型矩阵, 的对应元素相加而得到的矩阵称为矩阵A 的对应元素相加而得到的矩阵称为矩阵A 的和,记作A+B. 与B的和,记作A+B. 3 2 6 − 3 A = − 2 0 B = 3 7 例如: 例如: 5 4 4 − 2 −2 3 + 6 2 − 3 9 −1 则 A+ B = − 2 + 3 0 + 7 = 1 7 5 + 4 4 − 2 9 2 必须注意:只有两个同型的矩阵才能相加! 必须注意:只有两个同型的矩阵才能相加!
∑ | aij |
j =1
n
2
= ∑ aij ⋅ aij = 1
j =1
n
每一行的元素与另一行的对应元素的共轭复数 的乘积之和为0 的乘积之和为0
∑ aij ⋅ akj = 0
j =1
n
i≠k
9
而第三式说明酉矩阵的列是满足正交条件的. 而第三式说明酉矩阵的列是满足正交条件的. (2)酉矩阵的行列式的模是1. 酉矩阵的行列式的模是1 只须对等价定义中的第二式或第三式取行列式 即可知. 即可知.
其中 aij = a ji 这里的矩阵A 这里的矩阵A就 是一个对称矩阵 x1 而 x = x2 x 3
+ a31 x3 x1 + a32 x3 x2 + a33 x3 x3 a11 x1 + a12 x2 + a13 x3 = ( x1 x2 x3 ) a21 x1 + a22 x2 + a23 x3 a x + a x + a x 31 1 32 2 33 3 a11 a12 a13 x1 = ( x1 x2 x3 ) a21 a22 a23 x2 a a32 a33 x3 31 12 = xT Ax
T
对矩阵的共轭转置运算有下列性质: 对矩阵的共轭转置运算有下列性质: (1)( A*)* = A; ( 2)( A + B )* = A * + B*;
(3)(kA)* = kA*; (4)( AB)* = B* A*
特别注意(3)和 ! 特别注意 和(4)!
三、U矩阵与H矩阵 、U矩阵与H 矩阵与
10
2.H矩阵与对称矩阵 .H矩阵与对称矩阵
厄米特Hermite矩阵 矩阵 厄米特
如果方阵A 则称为H矩阵. 定义 如果方阵A满足 A = A, 则称为H矩阵. * 特别地, 是实矩阵,则是对称矩阵. 特别地,如A是实矩阵,则是对称矩阵. 即对称矩阵是特殊的厄米特矩阵. 即对称矩阵是特殊的厄米特矩阵. 由定义可知:H矩阵中, 由定义可知:H矩阵中,关于主对角线对称的元素 :H矩阵中 是互相共轭的复数, 是互相共轭的复数,即 aij = aji , 而对角线上元素全 是实数. 是实数. 而对称矩阵的条件可改写为 AT = A, 其关于主对角线对称的元素必相等, 其关于主对角线对称的元素必相等,即 aij = aji . 对称阵在讨论二次型时起着关键的作用. 对称阵在讨论二次型时起着关键的作用.
11
2 2 f ( x1 , x2 , x3 ) = a11 x1 + 2a12 x1 x2 + a22 x2 如: 2 + 2a13 x1 x3 + 2a23 x2 x3 + a33 x3 可以重新 f ( x1 , x2 , x3 ) = a11 x1 x1 + a12 x1 x2 + a13 x1 x3 排列为 + a21 x2 x1 + a22 x2 x2 + a23 x2 x3
2
矩阵加法的运算规律
(1) A + B = B + A;
− a11 − a 21 (3 ) − A = ⋯ −a m1 − a12 − a 22 ⋯ − am1
交换律
(2) ( A + B ) + C = A + (B + C ). 结合律
⋯ − a1 n ⋯ − a2 n = (− aij ), ⋯ ⋯ 称作矩阵A 称作矩阵A ⋯ − a mn 的负矩阵. 的负矩阵.
下面我们谈谈准对角矩阵的概念及最简单的性质. 下面我们谈谈准对角矩阵的概念及最简单的性质. 阶方阵, 设A为n阶方阵,则形如 阶方阵 A 1 O A2 A= , ⋱ O As n 的矩阵, 阶方阵, 的矩阵,其中 A 为 ni × ni 阶方阵, 1 + n2 +⋯+ ns = n i 表示适当类型的零矩阵. 而O表示适当类型的零矩阵. 这样的矩阵称为准对角矩阵或分块对角矩 阵.n ,⋯, n ) 型 (n1, 2 s 由名称即可知道, 由名称即可知道,准对角矩阵不过是分块矩阵的特殊 情形.关于一般的分块矩阵的内容请参看其他教材. 情形.关于一般的分块矩阵的内容请参看其他教材.
14
矩阵分块是为了简化对较高阶数的矩阵的计算. 矩阵分块是为了简化对较高阶数的矩阵的计算.因为 分块矩阵的运算非常简单. 分块矩阵的运算非常简单. 对同类型的分块对角矩阵 0 0 A B1 1 A2 B2 A= B = ⋱ ⋱ 0 0 Ak Bk 如果下面的运算都有意义,则运算性质有: 如果下面的运算都有意义,则运算性质有: 0 A + B1 1 A + B2 2 (1) A + B = ⋱ 0 Ak + Bk
1.U矩阵与正交矩阵 .U矩阵与正交矩阵
Unitary, Hermite matrix
下面要介绍的是另外几种特殊的矩阵. 下面要介绍的是另外几种特殊的矩阵. 若一个矩阵A 称为酉矩阵. 定义 若一个矩阵A满足条件 A* = A−1, 称为酉矩阵. 特别地,如果A是实矩阵,则称为正交矩阵. 特别地,如果A是实矩阵,则称为正交矩阵. 由此可知:正交矩阵不过是实的U矩阵. 由此可知:正交矩阵不过是实的U矩阵.
二、矩阵的转置与共轭转置 Transpose, conjugate
定义 把一个矩阵的行变成列,而列变成行所得到 把一个矩阵的行变成列, 的新的矩阵叫原矩阵的转置矩阵.矩阵A 的新的矩阵叫原矩阵的转置矩阵.矩阵A的 转置矩阵记作 A' 或 AT 3 2 3 − 2 5 T 例如: 例如: A = − 2 0 则 A'= A = 5 4 2 0 4 cos α − sin α cos α sin α R= R' = = R −1 sin α cos α − sin α cosα 旋转变换
、 为矩阵, 为数) (设 A、B为矩阵,λ , µ为数)
(1) (λµ ) A = λ (µA); (2) (λ + µ ) A = λA + µA; (3 ) λ ( A + B ) = λA + λB .
(4)λ ( AB ) = (λA) B = A(λB )
当然AB要有意义 当然 要有意义 矩阵相加与数乘矩阵合起来,统称为矩阵的线 矩阵相加与数乘矩阵合起来,统称为矩阵的线 性运算. 性运算. 4
U = a ⇒ U * = a,
UU * = E ⇒ U U * = a ⋅ a = a = 1
2
(3)两个同阶U矩阵的乘积是U矩阵. 两个同阶U矩阵的乘积是U矩阵.
− (U1U 2 )* = U 2 * U 1* = U 2 1U 1−1 = (U 1U 2 )−1