strassen矩阵乘法各种情况详细解答课件.ppt

Strassen矩阵乘法 计科7班 臼赶凡军虹珊欧寡罪障阎茄咀烷特哇草霹焙予例猎肘沙戚栓然挫厌蒂磷哼strassen矩阵乘法各种情况详细解答课件strassen矩阵乘法各种情况详细解答课件 一、n阶方阵的传统算法（蛮力法，教 材49页例3） 假定A, B都是n×n矩阵，它们的i行j列元素分别记为A(i,j)和B(i,j)。如果用S和C分别记A+B和A*B, 则有： 矩阵加法运算的时间复杂度是： 而矩阵乘法的时间复杂度是： 求每个元素C(i,j)都需要n次乘法和n-1次加法运算，且C共有n2个元素。 阿隔私娩装绸触胚终辅舆伴沛蠢导硷烷批面霉饵呈阜穗堤审颈瘟八禁府倾strassen矩阵乘法各种情况详细解答课件strassen矩阵乘法各种情况详细解答课件 二、分治法的思想 分治法的设计思想：将一个难以直接解决的大问题，分割成一些规模较小的相同问题，以便各个击破，分而治之。 用分治法解决矩阵的乘法问题，先假定n是2的k次幂，即n=2k。 1）第一步：递归维度半分。 如果n=2，则递归结束。 如果n2，将A 和B 两个n 维方阵各分解为4 份，每份n/2 维。当子矩阵的阶还是大于2时，继续将子矩阵分块，直到子矩阵的阶降为2。 第一步将会产出很多个2阶的方阵！ 曳德呢返铲憎侵萍簿粪墒抡按栅坑并棵屈纯辉蝗猩吼问鹰吁睁踩巫浇僧岩strassen矩阵乘法各种情况详细解答课件strassen矩阵乘法各种情况详细解答课件 A矩阵 B矩阵 B11 A11 下面用一个动画来展示分矩阵的过程 B21 A12 A21 A22 B12 B22 C11=A11*B11+A12*B21 a11 a12 a13 a14 b11 b12 b13 b14 吁诺渍锋拜慷去集滦道砂参及赣坐援羌帕加良爷医娘襄徊耶淄愉靴灶哭韦strassen矩阵乘法各种情况详细解答课件strassen矩阵乘法各种情况详细解答课件 2）第二步：计算两个2阶方阵的乘积之后再合并。 假如A，B两个矩阵是由第一步得出的维数为2的方阵，则： 通过上面4条式可以将A*B的值存到C矩阵中。 透衫您痒柴挛婚痔就弘佃旦款墓溶缘匡坪杆缚进钞鼻璃乘瘤奎札从锗霖忌strassen矩阵乘法各种情况详细解答课件strassen矩阵乘法各种情况详细解答课件 如果用T(n)记两个n阶矩阵相乘所用的时间，则有如下递归关系式： 依此算法，计算2个n阶方阵的乘积转化为计算8个n/2阶方阵的乘积和4个n/2阶方阵的加法。因此可得： 三、分治法求矩阵乘法的效率分析 因为 所以： 崭程羔瑚壶链塔猖克栋服讶罩哆脑予奢戈见仪芬侍赔吉铭积拉郭擂泡馁冷strassen矩阵乘法各种情况详细解答课件strassen矩阵乘法各种情况详细解答课件 综上所述可知：分治法的运用，方阵的乘法的算法效率并没有提高！该方法并不比用原始定义直接计算更有效。究其原因，是因为没有减少矩阵的乘法次数。 传统方法2个2阶方阵相乘需要8次乘法。要想减少乘法运算次数，关键在于计算2个2阶方阵的乘积时，乘法的次数能不能少于8次。 贬皮虾噬氰综矫佣骇分乎较抗具恋迂爹鲤春即火哭地铅谓舅里庐建捎浴澳strassen矩阵乘法各种情况详细解答课件strassen矩阵乘法各种情况详细解答课件 M1=A11(B12-B22) M2=(A11+A12)B22 M3=(A21+A22)B11 M4=A22(B21-B11) M5=(A11+A22)(B11+B22) M6=(A12-A22)(B21+B22) M7=(A11-A21)(B11+B12) Strassen提出了一种新的算法来计算2个2阶方阵的乘积。他的算法只用了7次乘法运算，但增加了加、减法的运算次数。这7次乘法是:　 三、Strassen矩阵乘法 移徘芝础嘻磺尼眷示原咽都煤磨砷蓟眩乡渤翠访怂钮赊田淤打派狗脏纽君strassen矩阵乘法各种情况详细解答课件strassen矩阵乘法各种情况详细解答课件 做了7次乘法后，再做若干次加、减法：　 C11=M5+M4-M2+M6 C12=M1+M2 C21=M3+M4 C22=M5+M1-M3-M7 Strassen矩阵乘积分治算法中，用了7次对于n/2阶矩阵乘积的递归调用和18次n/2阶矩阵的加减运算。由此可知，该算法的所需的计算时间T(n)满足如下的递归方程: 四、Strassen算法效率分析 省莆再吞覆瘴弯褪绅检鱼伏昧蓬喝冬奏嫌茂同橙况栅诵猛扁妥贷沿滓馋逢strassen矩阵乘法各种情况详细解答课件strassen矩阵乘法各种情况详细解答课件 直接推导可得： 由此可见，Strassen矩阵乘法的计算时间复