分块矩阵教案.pptVIP

1、矩阵分块的方法 1、分块加法 分块矩阵 一. 分块矩阵的运算规则 二. 分块矩阵的一些例子 矩阵分块，是矩阵运算的一个重要方法，可将大规模矩阵的运算化为若干小矩阵进行计算。 在矩阵某些行之间插入横线，某些列之间插入纵线，将矩阵分割成若干个小矩阵，每个小矩阵称为矩阵的子块；以子块为元素的矩阵，称为分块矩阵。 例如 即 说明 (1). 矩阵分块时，同一个矩阵可以有不同的分块方法，应根据需要进行选择。 2、矩阵分块一般形式 矩阵A = ( aij )m×n，在行方向分s块，列方向分t块，称A为s×t分块矩阵，第k行l列子块Akl是mk×nl阶矩阵。 各子块行数 各子块列数 说明 (2). 矩阵分块三原则：体现原矩阵特点，依据问题需要，子块可以作元素运算。 一、分块矩阵的运算规则 设A、B是m×n阶矩阵，采用相同的分块法分块将A、B分块如下： 注. 分块矩阵运算中，每个子块具有二重性：一是分块矩阵的元素；二是本身是矩阵。 2、分块数乘 设A是m×n阶矩阵，任意分块，k是常数，则定义 3、分块乘法 设A是m×l阶矩阵， B是l×n阶矩阵， 即A的列数 = B 的行数 即A的列分块法 = B 的行分块法 分块A = ( Auv )s×r B = ( Bvw )r×t 则A与B的乘积C = ( Cuw ) 是s×t阶分块矩阵，满足 注. 分块矩阵乘积AB中，每个子块： （1）作为分块阵元素参与运算 （2）作为矩阵也要满足乘法条件 例1. 用分块矩阵法求AB 解： 则 又 于是 说明 (3). 矩阵分块的目的，是让矩阵的计算过程更简单，计算量更少。 4、分块转置 设矩阵A = ( Aij ) 是s×r 阶分块矩阵 例1的计算量比较： 直接进行矩阵乘积需要的四则运算次数 用分块矩阵进行矩阵乘积需要的四则运算次数 合计32次 说明：分块转置中，每个子块一方面作为分块阵元素要转置；另一方面作为矩阵本身也要转置。 分外层、内层双重转置 特别地，对于列分块矩阵： 二、一些特殊的分块矩阵 1. 2阶分块上（下）三角形矩阵求逆 例2. 求下列2阶分块逆矩阵 其中A11， A22可逆矩阵 其中B12， B21可逆矩阵 解(1) ：设A的分块逆矩阵为 得到4个矩阵方程组 求解该方程组，得 (2) （解略，请仿(1)方法自行求解） 设A1, A2, … , As均为方阵（不一定同阶），则称下面的A为分块对角矩阵 2. 分块对角矩阵 如果矩阵A1, A2, … , As均可逆，则分块对角矩阵A可逆，且其逆矩阵为 说明：分块对角阵的逆矩阵，与对角矩阵的逆矩阵形式类似。 3. 矩阵乘积AB，A不分块，B按列分块 设矩阵A、B分别是s×n 和n×t 阶矩阵，A不分块，B按列分块，即 则 例3. 求解下列矩阵方程 说明：矩阵方程AX = B 可看成 t 个线性方程组 Ax1 = b1, Ax2 = b2, …, Axt = bt 其中B = ( b1, b2, …, bt )， X = ( x1, x2, …, xt ) 解：对增广矩阵( A, B )进行初等行变换 r2+r1 r3-2r1 -r2 r1-2r2 r3+r2 于是方程组Ax1 = b1有解 当且仅当λ= 0 时，Ax2 = b2有解 所以矩阵方程AX = B 在参数λ= 0 时，有解： 说明：利用增广矩阵的初等行变换，可以对矩阵方程AX = B 的 t 个线性方程组同时进行求解。 4. 矩阵乘积AB，A按列分块，B每个元素为块 （1）设矩阵A是s×n 矩阵，X 是n×1矩阵： 将A按列分块，即 则 我们将表达式 称为向量