矩阵n次方地几种求法地归纳
矩阵n次方的几种求法的归纳
矩阵n 次方的几种求法1.利用定义法()(),,ij kj s nn mA aB b ⨯⨯==则(),ij s mC c ⨯=其1122...ij i j i j in nj c a b a b a b =+++1nik kj k a b ==∑称为A 与B 的乘积,记为C=AB ,则由定义可以看出矩阵A与B 的乘积C 的第i 行第j 列的元素等于第一个矩阵A 的第i 行与第二个矩阵B 的第j 列的对应元素乘积之和,且由定义知:第一个矩阵的列数与第二个矩阵的行数要相[]1同。
例1:已知矩阵34125310210134A ⨯⎛⎫ ⎪=- ⎪ ⎪⎝⎭,445130621034510200B ⨯⎛⎫⎪⎪= ⎪⎪⎝⎭,求AB解:设C AB ==()34ij c ⨯,其中1,2,3i =;1,2,3,4j =由矩阵乘积的定义知:111526533032c =⨯+⨯+⨯+⨯=121122543231c =⨯+⨯+⨯+⨯= 131321553030c =⨯+⨯+⨯+⨯=14102051305c =⨯+⨯+⨯+⨯=21150623101c =-⨯+⨯+⨯+⨯= 22110224129c =-⨯+⨯+⨯+⨯= 23130125107c =-⨯+⨯+⨯+⨯= 24100021102c =-⨯+⨯+⨯+⨯= 310516334015c =⨯+⨯+⨯+⨯= 320112344222c =⨯+⨯+⨯+⨯= 330311354016c =⨯+⨯+⨯+⨯= 34001031403c =⨯+⨯+⨯+⨯=将这些值代入矩阵C 中得:C AB ==34323130519721522163⨯⎛⎫ ⎪⎪ ⎪⎝⎭则矩阵A 的n 次方也可利用定义的方法来求解。
2.利用矩阵的分块来求解这类方法主要是把一个大矩阵看成是由一些小矩阵组成,就如矩阵由数组成的一样在运算中将这些小矩阵当做数一样来处理,再由矩阵乘法的定义来求解这些小矩阵的乘积所构成的矩阵。
求矩阵的n次幂有如下几个常用方法
求矩阵的n次幂有如下几个常用方法1、求矩阵的n次幂的矩阵乘法法:求矩阵的n次幂的矩阵乘法法是用矩阵的乘法来求n次幂的一种方法,假设n>1。
令A为一个n阶矩阵,将A^n表示为A•A•…•A(n个A表示n次乘积),这样就可以用矩阵的乘法运算,把矩阵的n次幂表示出来。
这种方法适合任意阶数的矩阵,但是运算量大,一般在n大于4时会给计算机造成较大压力。
快速乘法法是将连乘拆成若干小段,用平方法计算这些小段,最后把平方结果合成出原来的积,这样就可以利用矩阵的平方法降低运算的复杂度,近似时间复杂度仅为O(logn)。
遗传算法(GA)是一种模拟自然辅助搜索算法,其可利用遗传运算(Genetic Operation)求解难以用传统算法求解的复杂问题,也可用来求矩阵的n次幂。
此方法通过使用遗传运算对n次幂矩阵A求解,其中有“选择(selection)”、“交叉(crossover)”、“变异(mutation)”等随机算法组成,在一定时间内,做出一定代数运算就能求出矩阵的n次幂,这种方法的效率取决于遗传算子的设计,但是因为这种方法涉及较少的运算,所以可能运算效率会很高。
线性矩阵分解法是把矩阵A事先分解成正交矩阵和对角矩阵的向量形式,将n次幂矩阵A^n分解成m分,从而减少计算量,缩短计算时间。
这种方法可以有效减少计算过程的数量,但对于大矩阵来说,可能由于分解矩阵的复杂度过高而无法令效率上升。
树结构法是一种求解n次方矩阵A的技术,它是建立树,由树的叶节点求出矩阵A的n次方。
由于每一层都有一个乘积,树结构法可以有效减少计算次数,较为高效。
通常来说,这种方法的复杂度降低到O(logn)。
总之,上面提到的几种方法都可以用来求矩阵的n次幂,根据矩阵的阶数和n的大小,可以合理选择合适的算法,从而提高求解效率。
矩阵n次方的几种求法的归纳
矩阵n 次方的几种求法1.利用定义法()(),,ij kj s nn mA aB b ⨯⨯==则(),ij s mC c ⨯=其1122...ij i j i j in nj c a b a b a b =+++1nik kj k a b ==∑称为A 与B 的乘积,记为C=AB ,则由定义可以看出矩阵A 与B 的乘积C 的第i 行第j 列的元素等于第一个矩阵A 的第i 行与第二个矩阵B 的第j 列的对应元素乘积之和,且由定义知:第一个矩阵的列数与第二个矩阵的行数要相[]1同。
例1:已知矩阵34125310210134A ⨯⎛⎫ ⎪=- ⎪ ⎪⎝⎭,445130621034510200B ⨯⎛⎫⎪⎪= ⎪⎪⎝⎭,求AB解:设C AB ==()34ij c ⨯,其中1,2,3i =;1,2,3,4j =由矩阵乘积的定义知:111526533032c =⨯+⨯+⨯+⨯=121122543231c =⨯+⨯+⨯+⨯= 131321553030c =⨯+⨯+⨯+⨯=14102051305c =⨯+⨯+⨯+⨯=21150623101c =-⨯+⨯+⨯+⨯= 22110224129c =-⨯+⨯+⨯+⨯= 23130125107c =-⨯+⨯+⨯+⨯= 24100021102c =-⨯+⨯+⨯+⨯= 310516334015c =⨯+⨯+⨯+⨯= 320112344222c =⨯+⨯+⨯+⨯= 330311354016c =⨯+⨯+⨯+⨯= 34001031403c =⨯+⨯+⨯+⨯=将这些值代入矩阵C 中得:C AB ==34323130519721522163⨯⎛⎫ ⎪⎪ ⎪⎝⎭则矩阵A 的n 次方也可利用定义的方法来求解。
2.利用矩阵的分块来求解这类方法主要是把一个大矩阵看成是由一些小矩阵组成,就如矩阵由数组成的一样在运算中将这些小矩阵当做数一样来处理,再由矩阵乘法的定义来求解这些小矩阵的乘积所构成的矩阵。
求矩阵的n次方的方法
求矩阵的n次方的方法简介矩阵的n次方是指将一个矩阵连乘n次的结果。
求矩阵的n次方是在很多数学和工程问题中都会遇到的核心计算任务之一。
本文将介绍几种常见的求矩阵的n次方的方法,包括矩阵乘法运算的定义、直接求解法、分治法以及特征分解法等。
不同的方法有不同的适用场景和时间复杂度,我们将对每种方法进行详细的探讨。
1. 矩阵乘法运算的定义在开始讨论求矩阵的n次方之前,我们首先需要了解矩阵乘法运算的定义。
给定两个矩阵A和B,它们的乘积AB定义为:这里的AB是一个n行p列的矩阵,其中第i行第j列的元素可以通过矩阵A的第i行和矩阵B的第j列的对应元素相乘并求和得到。
2. 直接求解法直接求解法是最直观也最容易理解的一种方法。
我们可以通过连乘n次矩阵A自身来求得矩阵的n次方,即。
具体的求解步骤如下: 1. 初始化一个单位矩阵I,它的大小与矩阵A相同。
2. 循环进行n次矩阵乘法运算,每次将结果保存在I中。
3. 当循环结束后,I即为矩阵A的n次方。
以下是使用直接求解法求解矩阵的n次方的示例代码:def matrix_power(A, n):I = [[1 if i == j else 0 for j in range(len(A))] for i in range(len(A))]for _ in range(n):I = matrix_multiply(I, A)return Idef matrix_multiply(A, B):n, m, p = len(A), len(A[0]), len(B[0])result = [[0 for _ in range(p)] for _ in range(n)]for i in range(n):for j in range(p):for k in range(m):result[i][j] += A[i][k] * B[k][j]return result直接求解法的时间复杂度为O(n^3)。
矩阵n次方的几种求法的归纳复习过程
矩阵n次方的几种求法的归纳收集于网络,如有侵权请联系管理员删除矩阵n 次方的几种求法1.利用定义法()(),,ij kj s nn mA aB b ⨯⨯==则(),ij s mC c ⨯=其1122...ij i j i j in nj c a b a b a b =+++1nik kj k a b ==∑称为A 与B 的乘积,记为,则由定义可以看出矩阵A 与B 的乘积C 的第i 行第j 列的元素等于第一个矩阵A 的第i 行与第二个矩阵B 的第j 列的对应元素乘积之和,且由定义知:第一个矩阵的列数与第二个矩阵的行数要相[]1同。
例1:已知矩阵34125310210134A ⨯⎛⎫ ⎪=- ⎪ ⎪⎝⎭,445130621034510200B ⨯⎛⎫⎪⎪= ⎪⎪⎝⎭,求解:设C AB ==()34ij c ⨯,其中1,2,3i =;1,2,3,4j = 由矩阵乘积的定义知:111526533032c =⨯+⨯+⨯+⨯=121122543231c =⨯+⨯+⨯+⨯= 131321553030c =⨯+⨯+⨯+⨯=14102051305c =⨯+⨯+⨯+⨯= 21150623101c =-⨯+⨯+⨯+⨯= 22110224129c =-⨯+⨯+⨯+⨯=收集于网络,如有侵权请联系管理员删除23130125107c =-⨯+⨯+⨯+⨯= 24100021102c =-⨯+⨯+⨯+⨯= 310516334015c =⨯+⨯+⨯+⨯= 320112344222c =⨯+⨯+⨯+⨯=330311354016c =⨯+⨯+⨯+⨯= 34001031403c =⨯+⨯+⨯+⨯=将这些值代入矩阵C 中得:C AB ==34323130519721522163⨯⎛⎫⎪⎪ ⎪⎝⎭ 则矩阵A 的n 次方也可利用定义的方法来求解。
2.利用矩阵的分块来求解这类方法主要是把一个大矩阵看成是由一些小矩阵组成,就如矩阵由数组成的一样在运算中将这些小矩阵当做数一样来处理,再由矩阵乘法的定义来求解这些小矩阵的乘积所构成的矩阵。
矩阵n次方通用解法
矩阵n次方通用解法介绍矩阵的n次方运算是矩阵乘法的重要应用之一,它在数学、计算机科学和工程领域都有广泛的应用。
本文将深入探讨矩阵n次方的通用解法,包括计算过程、优化方法以及一些应用案例。
矩阵乘法回顾在进一步探讨矩阵n次方之前,我们先回顾一下矩阵乘法。
对于两个矩阵A和B,它们的乘积C可以通过以下公式计算:C = A * B其中,A是一个m行n列的矩阵,B是一个n行p列的矩阵,C是一个m行p列的矩阵。
矩阵乘法的计算规则是,C的第i行第j列元素等于A的第i行与B的第j 列对应元素的乘积之和。
矩阵的1次方和0次方矩阵的1次方就是矩阵本身,即:A^1 = A。
矩阵的0次方定义为单位矩阵,即:A^0 = I。
矩阵的n次方对于一个矩阵A,它的n次方可以通过连续进行n次矩阵乘法来计算,即:A^n = A * A * A * … * A然而,直接按照这种方法计算矩阵的n次方在效率上并不高。
接下来,我们将介绍一个通用解法,可以更高效地计算矩阵的n次方。
矩阵的n次方通用解法为了高效计算矩阵的n次方,我们可以利用矩阵乘法的性质。
假设我们要计算矩阵A的2n次方,即A(2^n)。
我们可以通过以下步骤来逐步计算:1.计算 A2、A4、A^8、…,直到 A(2n)。
–这可以通过每次将矩阵平方来实现,即 A(2i) = (A(2(i-1)))^2,其中i从1递增到n。
2.根据 A(2n) 的定义,将其展开为累积乘积的形式,即:–A(2n) = A(2(n-1)) * A(2(n-1)) * … * A(2(n-1)),总共有 2^(n-1) 个 A(2(n-1))。
通过以上步骤,我们可以高效地计算矩阵的n次方。
下面是一个具体的计算演示:以计算矩阵A的8次方为例,即 A^8。
根据通用解法,我们先计算出 A2、A4 和 A^8,然后根据 A^8 的定义展开累积乘积。
具体计算过程如下:1.计算 A^2:–A^2 = A * A2.计算 A^4:–A^4 = (A^2) * (A^2)3.计算 A^8:–A^8 = (A^4) * (A^4)4.展开 A^8 的累积乘积:–A^8 = A^4 * A^4–A^8 = (A^2 * A^2) * (A^2 * A^2)–A^8 = (A * A) * (A * A) * (A * A) * (A * A)通过以上计算,我们可以得到矩阵A的8次方。
矩阵n次方的几种求法的归纳
1
32 31 30 5
C
AB
=
1
9
7
2
15 22 16 3 34
则矩阵 A 的 n 次方也可利用定义的方法来求解。
2.利用矩阵的分块来求解
这类方法主要是把一个大矩阵看成是由一些小矩阵组成,就如矩阵
由数组成的一样在运算中将这些小矩阵当做数一样来处理,再由矩阵
乘法的定义来求解这些小矩阵的乘积所构成的矩阵。即设
sin
sin cos
3
cos
sin
sin 2 cos
cos
sin
sin
cos
cos cos
2 2
cos sin
sin 2 sin 2
sin cos
cos 2 sin sin 2 cos
cos 2 cos sin 2 sin
cos 3 sin 3
sin 3
B11 B21
,其中
E11
1 0 0
0 1 0
0
0 1
0 6 0 6
2 5
1 2
A12
1
2
3
,
A22
0
8
6
,
B11
4
1
5 0
,
B21
4 0
2
6
由矩阵乘积法则知:
AB=
B11 A12 B21 A21B11 A22 B21
42
由矩阵加法和乘积法则 知1 :
4 2 0 6 52
1 0 0 2 5 1 0 0 2 5
解:将
A
0 0
1 0
0 1
0
0
0
1
3
求矩阵a的n次方的方法
求矩阵a的n次方的方法求矩阵a的n次方是一个非常常见的线性代数问题,在很多数学和工程领域都有广泛的应用。
矩阵的n次方表示将矩阵自乘n次,即a * a * a * ... * a。
本文将一步一步地回答如何求解矩阵的n次方,以及该方法的应用。
一、矩阵的乘法在介绍求解矩阵的n次方之前,让我们先回顾一下矩阵的乘法。
对于两个相容的矩阵A和B,它们的乘积AB的定义如下:[AB]ij = Σ(Aik * Bkj),其中k为1到n的取值范围。
简而言之,AB的第i行第j列的元素等于A的第i行的元素与B的第j列的元素逐一相乘后再求和。
二、矩阵的平方在求解矩阵的n次方之前,我们先看一个特殊情况:矩阵的平方。
矩阵的平方是将矩阵自身乘以自己,即A的平方为A * A。
算法如下:1. 建立一个与A维度相同的零矩阵C;2. 对于C的第i行第j列的元素,将A的第i行与A的第j列进行乘法运算,得到C的第i行第j列的元素;3. 返回矩阵C,即A的平方。
三、任意正整数次方的求解接下来是核心部分:求解任意正整数次方的矩阵。
我们可以利用矩阵的平方来递归求解。
算法如下:1. 如果n等于1,直接返回矩阵A;2. 如果n是偶数,将n除以2,得到商为m,然后计算A的平方的m次方,即A的2^m次方;3. 如果n是奇数,将n减1变成偶数,并继续计算A的n-1次方;4. 将A的2^m次方乘以A的n-1次方,得到A的n次方。
这个算法利用了二分法,将n的复杂度从n降低到了对数级别。
通过递归地求解A的平方的平方,可以大幅度减小计算量,并且提高了效率。
四、应用举例矩阵的n次方的求解方法在数值计算、线性回归分析、图像处理等领域有广泛的应用。
下面举例说明其中两个应用。
1. 图像处理在图像处理中,矩阵通常被用来表示像素值。
通过将像素矩阵与特定的变换矩阵相乘,可以实现图像的旋转、缩放、平移等操作。
这里的变换矩阵就是原始像素矩阵的n次方。
通过求解原始像素矩阵的n次方,可以实现对图像的高效处理。
矩阵n次方通用解法
矩阵n次方通用解法矩阵n次方通用解法矩阵是线性代数中的重要概念,它在数学、物理、计算机科学等领域都有广泛的应用。
矩阵的n次方也是一个重要的问题,因为它涉及到很多实际问题中的计算。
本文将介绍矩阵n次方通用解法。
一、矩阵乘法在介绍矩阵n次方通用解法之前,我们需要先了解矩阵乘法。
对于两个矩阵A和B,它们的乘积C为:C(i,j) = ∑(k=1 -> n)A(i,k)*B(k,j)其中C(i,j)表示C矩阵第i行第j列元素,n表示A和B的列数相同。
二、暴力求解最简单的方法是通过暴力求解来计算矩阵n次方。
例如,对于一个2x2的矩阵A和一个正整数n,我们可以通过以下方式计算A^n:result = Afor i in range(n-1):result = result * A这种方法可以得到正确的结果,但是时间复杂度为O(n^3),当n较大时会非常耗时。
三、分治法分治法是一种常见的优化算法,在计算矩阵n次方时也可以使用。
假设我们要计算A^n,我们可以将其分解为两个子问题:计算A^(n/2)和(A^(n/2))^2。
然后再通过矩阵乘法将两个子问题的结果合并起来即可得到A^n。
该算法的时间复杂度为O(n^3logn),比暴力求解要快很多。
四、矩阵快速幂矩阵快速幂是一种更加高效的算法,它可以将时间复杂度降低到O(n^3logn)。
具体来说,我们可以先将指数n转换为二进制形式,例如:n = 13 -> 1101然后根据二进制形式中1的位置来计算矩阵的乘积。
以计算A^13为例,我们可以这样做:result = Ibase = Afor i in range(k):if n & (1 << i):result = result * basebase = base * base其中I表示单位矩阵,k表示二进制位数。
该算法的时间复杂度为O(n^3logn),比分治法还要快一些。
五、应用举例矩阵n次方通用解法在实际问题中有广泛应用。
求矩阵的n次方的方法
求矩阵的n次方的方法一、前言矩阵是数学中的一个重要概念,它在线性代数、微积分等领域中都有广泛的应用。
其中,求矩阵的n次方是一个常见的问题。
本文将介绍几种常用的方法来求矩阵的n次方。
二、矩阵乘法在介绍求矩阵的n次方之前,我们先来回顾一下矩阵乘法。
假设有两个矩阵A和B,它们分别是m×k和k×n的矩阵,则它们的乘积C=A×B是一个m×n的矩阵,其中C[i][j]表示A[i][1]×B[1][j]+A[i][2]×B[2][j]+...+A[i][k]×B[k][j]。
三、暴力法最简单直接的方法就是暴力法,即将原始矩阵连乘n次。
假设原始矩阵为A,则其n次方为An=A×A×...×A(共n个A相乘)。
这种方法虽然简单易懂,但时间复杂度为O(nm3),当n或m较大时会非常耗时。
四、幂运算为了提高效率,我们可以使用幂运算来计算矩阵的n次方。
假设原始矩阵为A,则其n次方可以表示为An=A^(log2(n))×A^(log2(n))×...×A^(log2(n))×A^(n-2^k)(其中k=log2(n)),即将n表示成二进制数,将每一位对应的幂运算结果相乘。
例如,当n=13时,13的二进制为1101,那么An=A^8×A^4×A^1。
这种方法的时间复杂度为O(m3log2(n)),相比暴力法有了很大的提升。
五、分治法分治法也是一种常用的方法来求矩阵的n次方。
假设原始矩阵为A,则我们可以将其划分成四个子矩阵:A11、A12、A21和A22,每个子矩阵都是原始矩阵的一部分。
则原始矩阵的n次方可以表示为:An = ( A11^n A12^n )( A21^n A22^n )其中,每个子矩阵的n次方可以通过递归调用求解。
具体地,我们可以按以下步骤来计算An:1. 将原始矩阵划分成四个子矩阵;2. 递归计算每个子矩阵的n/2次方;3. 根据公式计算An。
待定系数法求矩阵n次方
待定系数法求矩阵n次方
待定系数法是一种求解矩阵n次方的方法,主要应用于求解矩阵的高次幂。
以下是待定系数法求矩阵n次方的具体步骤:
1. 假设矩阵A的n次方可以表示为:An = A + c1A^2 + c2A^3 + ... + cn-1A^(n-1) + cnA^n,其中ci(i=1,2,...,n)为待定系数。
2. 求解方程组。
根据矩阵乘法的性质,An与A的每一列有关,因此可以列出n个方程,分别表示An的每一列与A的每一列的关系。
3. 解方程组,得到待定系数ci。
可以使用高斯消元法、矩阵分解法等求解方程组。
4. 将求得的待定系数ci代入假设式中,得到矩阵A的n次方:An =
A + c1A^2 + c2A^3 + ... + cn-1A^(n-1) + cnA^n。
5. 验证结果。
将求得的An代入原矩阵方程,检查是否满足矩阵乘法的性质。
需要注意的是,待定系数法适用于求解方阵矩阵的高次方,对于非方阵矩阵,需要先将其化为方阵形式,然后再应用待定系数法。
此外,待定系数法在求解过程中涉及到矩阵的多次乘法运算,可能较为繁琐,但对于一些特殊情况(如矩阵特征值、特征向量已知)来说,待定系数法是一种有效的求解方法。
在MATLAB中,可以使用命令求解矩阵的n次方。
例如,对于一个2阶矩阵A,可以按照以下方式计算其平方:
MATLAB
A = [1 2; 3 4];
An = A^2
待定系数法也可以应用于求解高次幂,只需将待定系数ci代入矩阵乘法公式即可。
矩阵a的n次方
矩阵a的n次方
一般有以下几种方法:
1、计算A^2,A^3找规律,然后用归纳法证明。
2、若r(A)=1,则A=αβ^T,A^n=(β^Tα)^(n-1)A
注:β^Tα =α^Tβ = tr(αβ^T)
3、分拆法:A=B+C,BC=CB,用二项式公式展开。
适用于B^n易计算,C的低次幂为零:C^2或C^3=0
4、用对角化 A=P^-1diagP
A^n = P^-1diag^nP
扩展资料:
将一个矩阵分解为比较简单的或具有某种特性的若干矩阵的和或乘积,矩阵的分解法一般有三角分解、谱分解、奇异值分解、满秩分解等。
在线性代数中,相似矩阵是指存在相似关系的矩阵。
相似关系是两个
矩阵之间的一种等价关系。
两个n×n矩阵A与B为相似矩阵当且仅当存
在一个n×n的可逆矩阵P。
一个矩阵A的列秩是A的线性独立的纵列的极大数目。
类似地,行秩
是A的线性无关的横行的极大数目。
通俗一点说,如果把矩阵看成一个个
行向量或者列向量,秩就是这些行向量或者列向量的秩,也就是极大无关
组中所含向量的个数。
求矩阵n次幂的方法
求矩阵n次幂的方法
矩阵的幂可以通过递归方式计算,也可以使用矩阵乘法实现。
下面分别介绍两种方法:
1. 递归方式
设矩阵 A 的 n 次幂为 A^n,则 A^1 等于矩阵 A 本身,A^2 等于 A 乘以 A^1,即 A^2 = A * A^1,以此类推。
因此,求矩阵 A 的n 次幂可以递归地进行以下计算:
a1 = a * a^(n-1)
a2 = a^2 * a^(n-2)
a3 = a^3 * a^(n-3)
...
an = an-1 * a^(n-1)
其中,a 为矩阵 A 的某一行或列向量,a^(n-1) 为矩阵 A 的
n-1 次幂。
当递归计算到 n=1 时,返回矩阵 A 本身;当 n>1 时,返回矩阵 A^n。
2. 矩阵乘法实现
另一种方法是使用矩阵乘法来计算矩阵的幂。
设矩阵 A 的 n 次幂为 A^n,则 A^n 可以表示为:
A^n = (A*A*...*A)^(n-1)
其中,*表示矩阵乘法,^(n-1) 表示 n-1 次方。
因此,求矩阵 A 的 n 次幂可以依次执行以下计算:
a1 = a * a^(n-1)
a2 = a^2 * a^(n-2)
a3 = a^3 * a^(n-3)
...
an = an-1 * a^(n-1)
其中,a 为矩阵 A 的某一行或列向量,a^(n-1) 为矩阵 A 的n-1 次幂。
这种方法的优点是在大规模矩阵计算中效率较高,但需要使用额外的内存来存储乘数矩阵。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
矩阵n 次方的几种求法1.利用定义法()(),,ij kj s nn mA aB b ⨯⨯==则(),ij s mC c ⨯=其1122...ij i j i j in nj c a b a b a b =+++1nik kj k a b ==∑称为A 与B 的乘积,记为C=AB ,则由定义可以看出矩阵A 与B 的乘积C 的第i 行第j 列的元素等于第一个矩阵A 的第i 行与第二个矩阵B 的第j 列的对应元素乘积之和,且由定义知:第一个矩阵的列数与第二个矩阵的行数要相[]1同。
例1:已知矩阵34125310210134A ⨯⎛⎫ ⎪=- ⎪ ⎪⎝⎭,445130621034510200B ⨯⎛⎫⎪⎪= ⎪⎪⎝⎭,求AB解:设C AB ==()34ij c ⨯,其中1,2,3i =;1,2,3,4j =由矩阵乘积的定义知:111526533032c =⨯+⨯+⨯+⨯=121122543231c =⨯+⨯+⨯+⨯= 131321553030c =⨯+⨯+⨯+⨯=14102051305c =⨯+⨯+⨯+⨯=21150623101c =-⨯+⨯+⨯+⨯= 22110224129c =-⨯+⨯+⨯+⨯= 23130125107c =-⨯+⨯+⨯+⨯= 24100021102c =-⨯+⨯+⨯+⨯= 310516334015c =⨯+⨯+⨯+⨯= 320112344222c =⨯+⨯+⨯+⨯= 330311354016c =⨯+⨯+⨯+⨯= 34001031403c =⨯+⨯+⨯+⨯=将这些值代入矩阵C 中得:C AB ==34323130519721522163⨯⎛⎫ ⎪⎪ ⎪⎝⎭则矩阵A 的n 次方也可利用定义的方法来求解。
2.利用矩阵的分块来求解这类方法主要是把一个大矩阵看成是由一些小矩阵组成,就如矩阵由数组成的一样在运算中将这些小矩阵当做数一样来处理,再由矩阵乘法的定义来求解这些小矩阵的乘积所构成的矩阵。
即设()(),,ij kj s nn mA aB b ⨯⨯==把A ,B 分解成一些小矩阵:1111l t tl A A A A A ⎛⎫ ⎪=⎪ ⎪⎝⎭,1111r l lr B B B B B ⎛⎫⎪= ⎪ ⎪⎝⎭,其中ij A 是i j s n ⨯小矩阵且1,2...i t =,1,2...j l =,且12...t s s s s +++= ,12...l n n n n +++=;ij B 是j k n m ⨯小矩阵且1,2...j l =,1,2...k r =;且12...l n n n n +++=,12...r m m m m +++=;令C AB ==1111r t tr C C C C ⎛⎫⎪⎪ ⎪⎝⎭,其中ij C 是i j s m ⨯小矩阵且1,2...i t =,1,2,...,j r =,且12...t s s s s +++=,12...r m m m m +++=;其中1122...ij i j i j il lj C A B A B A B =+++。
这里我们应注意:矩阵A 列的分法必须与矩阵B 行的分法一[]1致。
例2:已知矩阵4510025010130012800006A ⨯⎛⎫ ⎪⎪= ⎪ ⎪⎝⎭,521245104206B ⨯⎛⎫⎪ ⎪⎪=⎪ ⎪ ⎪⎝⎭,求AB 解:将45451002510025010130101300128001280000600006A ⨯⨯⎛⎫⎛⎫⎪⎪ ⎪⎪== ⎪ ⎪ ⎪⎪⎪⎝⎭⎝⎭11122122EA A A ⎛⎫⎪⎝⎭写成 12124545101042420606B ⎛⎫⎛⎫⎪ ⎪ ⎪ ⎪ ⎪ ⎪== ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭1121B B ⎛⎫ ⎪⎝⎭写成,其中11100010001E ⎛⎫⎪= ⎪ ⎪⎝⎭ 12251328A ⎛⎫ ⎪= ⎪ ⎪⎝⎭,()2206A =,11124510B ⎛⎫⎪= ⎪ ⎪⎝⎭,214206B ⎛⎫= ⎪⎝⎭由矩阵乘积法则知:AB=1112212111222142B A B A B A B ⨯+⎛⎫⎪+⎝⎭由矩阵加法和乘积法则[]1知:42936825AB 952036⨯⎛⎫ ⎪⎪= ⎪ ⎪⎝⎭则矩阵A 的n 次方的求解也可利用以上方法来求解。
3.利用数学归纳法求解这种方法与矩阵定[]1义和数学归纳[]3法相结合,从而找出规律再求解,但是这种方法比较适合低阶且有规律的方阵n 次方的运[]2算。
例3:已知A=cos sin sin cos θθθθ-⎛⎫⎪⎝⎭,求nA 解:当2n =时2cos sin cos sin cos sin sin cos sin cos sin cos θθθθθθθθθθθθ---⎛⎫⎛⎫⎛⎫= ⎪ ⎪⎪⎝⎭⎝⎭⎝⎭2222cos 2sin 2cos sin 2cos sin sin 2cos 22cos sin cos sin θθθθθθθθθθθθ-⎛⎫--⎛⎫==⎪ ⎪-⎝⎭⎝⎭ 当3n =时32cos sin cos sin cos sin sin cos sin cos sin cos θθθθθθθθθθθθ---⎛⎫⎛⎫⎛⎫=⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭cos 2cos sin 2sin cos 2sin sin 2cos cos 2sin sin 2cos cos 2cos sin 2sin θθθθθθθθθθθθθθθθ---⎛⎫=⎪+-⎝⎭cos3sin 3sin 3cos3θθθθ-⎛⎫=⎪⎝⎭所以假设n A =cos sin sin cos n n n n θθθθ-⎛⎫⎪⎝⎭当1k =时成立,假设当1k n =-时成立;则当k n =时1cos sin cos sin sin cos sin cos n n A θθθθθθθθ---⎛⎫⎛⎫=⎪ ⎪⎝⎭⎝⎭()()()()cos 1sin 1cos sin sin 1cos 1sin cos n n n n θθθθθθθθ---⎛⎫-⎛⎫=⎪⎪--⎝⎭⎝⎭由矩阵乘法定及三角函数知:n A =cos sin sin cos n n n n θθθθ-⎛⎫⎪⎝⎭则假设成立。
所以n A =cos sin sin cos n n n n θθθθ-⎛⎫⎪⎝⎭4.利用分拆法求解这类方法主要是将一个矩阵分解成一个单位矩阵和另外一个矩阵之和再求[]1解,且另外这个矩阵的n 次方计算起来比较简[]2单。
例4:已知A=110011001⎛⎫ ⎪⎪ ⎪⎝⎭,求n A解:A E B =+,其中010001000B ⎛⎫⎪= ⎪ ⎪⎝⎭,矩阵E 为单位阵且2E E =EB BE B ==;故 n A =()122+C C C nn nn n n E B E B B B +=+++由2010010001001001000000000000B ⎛⎫⎛⎫⎛⎫⎪⎪ ⎪== ⎪⎪ ⎪ ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭23010010001010001001000001000000000000B ⎛⎫⎛⎫⎛⎫⎛⎫⎪ ⎪ ⎪⎪== ⎪ ⎪ ⎪⎪ ⎪ ⎪ ⎪⎪⎝⎭⎝⎭⎝⎭⎝⎭000000000⎛⎫⎪= ⎪ ⎪⎝⎭则3n ≥时,n B =0。
故122n nn A E C B C B =++由矩阵加法运算法则[]1知:n A =211011001n n C n ⎛⎫+ ⎪+ ⎪ ⎪⎝⎭5.利用相似矩阵求解(利用对角矩阵来求)定义:设矩阵A ,B 为数域P 上两个n 级矩阵,如果可以找到数域P 上的n 级可逆阵X ,使得矩阵1B X AX -=,就说A 与B 相[]1似。
如果矩阵A 或B 有一个可以化成对角矩阵则计算比较简便。
而判断矩阵A 可对角化的条件[]1有:1)矩阵A 可对角化的必要条件是矩阵A 有n 个不同的特征值2)矩阵A 可对角化的充要条件是矩阵A 有个n 线性无关的特征向量 3)在复数域上矩阵A 没有重根而求矩阵A 的特征值和特征向量的方法[]1有:1)求矩阵A 特征多项式E A λ-在数域P 中的全部根,这些根是矩阵A 的全部特征值。
把这些所求的特征值逐个的代入方程组()0E A X λ-=中,对于每一个特征值,解方程组()0E A X λ-=,求出一组基础解系,那么这个基础解系就是属于这个特征值的特征向量。
再利用判别法判断矩阵A 是否可对角化。
例5:已知矩阵33122212221A ⨯⎛⎫ ⎪=- ⎪ ⎪--⎝⎭,求n A解:易知矩阵的A 特征多项式E A λ-=122212221λλλ------由行列式计算方法知:E A λ-=()()()()()213113λλλλλ--=-+-所以矩阵A 的特征值为1,1,3-。
当特征值为1时,解方程()0E A X -=,由齐次线性方程组的计算方法知:()0E A X -=的基础解系为1a =()111'-;所以矩阵A 属于特征值1的全部特征向量为()1111k '-,其中1k ≠0。
当特征值为1-时,解方程()0E A X --=,由齐次线性方程组的计算方法知:()0E A X --=的基础解系为2a =()110'-;所以矩阵A 属于特征值1-的全部特征向量为()2110k '-,其中2k ≠0。
当特征值为3时,解方程()30E A X -=,由齐次线性方程组的计算方法知:()30E A X -=的基础解系为3a =()011'-,所以矩阵A 属于特征值3的全部特征向量为()3011k '-,其中3k ≠0。
则由矩阵A 可对角化的条件知:矩阵A 可对角化且对角阵为B =100010003⎛⎫⎪- ⎪ ⎪⎝⎭令123C a a a →→→'⎛⎫= ⎪⎝⎭=33110111101⨯⎛⎫⎪-- ⎪ ⎪-⎝⎭,由求逆矩阵的方法知:1111011110C -⎛⎫ ⎪=-- ⎪ ⎪⎝⎭因为线性变换在不同基下所对应的矩阵是相似的知:1C AC B -= 所以()11nn n C AC C A C B --==,则()333310010001001000303nnn n B ⨯⨯⎛⎫⎛⎫⎪ ⎪=-=- ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭ 由1n n A CB C -=,由矩阵的乘法运算法则知:()()()()3311111311311113131nnn n nn nn n A ⨯⎛⎫---- ⎪ ⎪=--+--- ⎪-- ⎪⎝⎭2)对方阵A ,设()()1F E A λλ'=-,对()()1n F E λ做初等变换,化成()()()D P λλ其中()D λ为上三角阵,则矩阵()D λ主对角线上元素乘积的λ的多项式的根即为A 的特征根i λ。
对矩阵A 的任一特征根i λ,代入()()()D P λλ中,若()i D λ中非零向量构成一满秩矩阵,则()i D λ行向量所对应的()i P λ中的行向量i ξ即为i λ的特征向量;否则,继续施行初等行变换,使得()i D λ中非零向量构成一满秩矩阵,则()i D λ中零向量所对应的()i P λ中的行向量i ξ即为i λ的特征向[]8量。