矩阵相乘的算法设计.docx

《矩阵相乘的算法设计.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数据结构与算法设计课程实验报告课题矩阵相乘的算法设计专业班级网工专业1405班学号姓名陈晓露指导教师陶跃进目录一、问题描述二、问题分析1、分析最优解的结构2、建立递归关系3、递归实现的复杂性4、算法迭代实现三、结果输出四、实验总结一、问题描述给定n个矩阵｛A1，A2，...，An｝，其中这n个矩阵是可相乘的，i＝1，2，...，n-1。算出这n个矩阵的相乘积A1A2。。。An。补充：如果两个矩阵A和B是可相乘的，那么A的列数要和B的行数是相同的，否则，这两个矩阵是不可相乘的。它们的相乘结果矩阵C的

2、行数是A的行数，而列数是B的列数。设A1,A2,…,An为矩阵序列，Ai为Pi-1×Pi阶矩阵，i=1,2,…,n.确定乘法顺序使得元素相乘的总次数最少.输入：向量P=实例： P=<10,100,5,50> A1:10×100,A2:100×5,A3:5×50乘法次序：(A1A2)A3:10×100×5+10×5×50=7500A1(A2A3):10×100×50+100×5×50=75000搜索空间的规模先将矩阵链加括号分为两部分，即P=A1*A2*...*An=（A1*

3、A2...*Ak）*（Ak+1*...An），则有f(n)=f(1)*f(n-1)+f(2)*f(n-2)+...+f(n-1)*f(1)种方法。动态规划算法输入P=，Ai..j表示乘积AiAi+1…Aj的结果，其最后一次相乘是：m[i,j]表示得到Ai..j的最少的相乘次数。递推方程：为了确定加括号的次序，设计表s[i,j],记录求得最优时最一位置。二、问题分析由于矩阵乘法满足结合律，故连乘积的计算可以有许多不同的计算次序。这种计算次序可以用加括号的方式来确定。若一个矩阵

4、连乘积的计算次序已完全确定，也就是说该连乘积已完全加括号，则我们可以通过反复调用两个矩阵相乘的标准算法计算出矩阵连乘积。1.分析最优解的结构为了方便起见，我们将矩阵连乘AiAi+1。。。Aj记为A[i:j]。经分析，计算A[1:n]的一个最优次序所包含的计算矩阵子链A[1:k]和A[k:n]的次序也是最优的。因此，矩阵连乘计算次序问题的最优解包含着子问题的最优解。2.建立递归关系用矩阵m[n][n]来存放A[i:j]相乘的计算次数，用p[n+1]用来存放矩阵的行数和列数。const int N=5

5、;  int m[N][N]; //m[i][j]存储Ai到Aj的最小乘法次数  int s[N][N];//s[i][j]存储Ai到Aj之间加括号的位置    int RecurMatrixChain(int P[],int i,int j)  {      m[i][j]=;      s[i][j]=i;      if(i==j)          m[i][j]=0;      else{          for(int k=i;k

6、=RecurMatrixChain(P,i,k)+RecurMatrixChain(P,k+1,j)+P[i]*P[k+1]*P[j+1];              if(q

7、,15,5,10,20};      for(int i=0;i

8、.1，所以很多子问题被重复计算了多次。于是，我们想到用自底向上的迭代实现。4.算法迭代实现迭代实现主要思想是子问题由小到大，每个子问题只计算一次，并且把结果保存起来，后来用到这个子问题时，直接代入。void MatrixChain(int P[],int n)  {      int r,i,j,k,t;      for(i=0;i