网站大量收购闲置独家精品文档,联系QQ:2885784924

《矩阵连乘实验报告.docxVIP

  1. 1、本文档共8页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
华北电力大学科技学院实 验 报 告实验名称 矩阵连乘问题 课程名称 计算机算法设计与分析 专业班级:软件12K1 学生姓名:吴旭学 号:121909020124 成 绩:指导老师:刘老师 实验日期:2014.11.14实验内容矩阵连乘问题,给定n个矩阵{A1,A2,…,An},其中Ai与Ai+1是可乘的,i=1,2,3…,n-1。考察这n个矩阵的连乘A1,A2,…,An。主要思想由于矩阵乘法满足结合律,故计算矩阵的连乘积可以有许多不同的计算次序。这种计算次序可以用加括号的方式来确定。若一个矩阵连乘积的计算次序完全确定,也就是说该连乘积已经完全加括号,则可依此次序反复调用2个矩阵相乘的标准算法计算出矩阵连乘积。完全加括号的矩阵连乘积可递归的定义为:单个矩阵是完全加括号的;矩阵连乘积A是完全加括号的,则A可表示为2个完全加括号的矩阵连乘积B和C的乘积并加括号,即A=(BC)。运用动态规划法解矩阵连乘积的最优计算次序问题。按以下几个步骤进行分析最优解的结构设计求解具体问题的动态规划算法的第1步是刻画该问题的最优解的结构特征。为方便起见,将矩阵连乘积简记为A[i:j]。考察计算A[1:n]的最优计算次序。设这个计算次序矩阵在Ak和Ak+1之间将矩阵链断开,,则其相应的完全加括号方式为((A1…Ak)(Ak+1…An))。依此次序,先计算A[1:k]和A[k+1:n],然后将计算结果相乘得到A[1:n]。建立递归关系设计动态规划算法的第二步是递归定义最优值。对于矩阵连乘积的最优计算次序问题,设计算A[i:j],,所需的最少数乘次数为m[i][j],原问题的最优值为m[1][n]。当i=j时,A[i:j]=Ai为单一矩阵,无需计算,因此m[i][i]=0,i=1,2,…n。当ij时,可利用最优子结构性质来计算m[i][j]。m[i][j]=m[i][k]+m[k+1][j]+pi-1pkpj。由于在计算时并不知道断开点k的位置,所以k还未定。计算最优值根据计算m[i][j]的递归式,容易写一个递归算法计算m[1][n]。动态规划法解决此问题,可依据递归式以自底向上的方式进行计算,在计算过程中保存已解决的子问题答案。每个子问题只计算一次,而在后面需要时只要简单查一下,从而避免大量的重复计算,最终得到多项式时间的算法matrixChain。(见实验代码部分)构造最优解算法matrixChain只计算出最优值,并没有给出最优解。但是matrixChain已经记录了构造最优解所需的全部信息。S[i][j]中的数表明,计算矩阵链A[i:j]的最佳方式应在矩阵Ak和Ak+1之间断开,最优加括号方式为(A[i:k])(A[k+1:j])。依次构造最优解。(算法见实验代码部分)实验结果结果验证对实验结果进行验证,4个矩阵分别是A1[35*15],A2[15*5],A3[5*10],A4[10*20]。依递归式有:M[1][4]=min=7125 且k=3。计算结果正确,证明所编写的程序可正确算出最优解。实验代码#includestdio.h#define N 100//定义最大连乘的矩阵个数是100void matrixChain(int p[],int m[N+1][N+1],int s[N+1][N+1])/*用m[i][j]二维数组来存储Ai*.....Aj的最少数乘次数,用s[i][j]来存储使Ai.....Aj获得最少数乘次数对应的断开位置k,需要注意的是此处的N+1非常关键,虽然只用到的行列下标只从1到N,但是下标0对应的元素默认也属于该数组,所以数组的长度就应该为N+1*/{int n=N;//定义m,s数组的都是n*n的,不用行列下标为0的元素,但包括在该数组中for(int i=1;i=n;i++)m[i][i]=0;/*将矩阵m的对角线位置上元素全部置0,此时应是r=1的情况,表示先计算第一层对角线上个元素的值*/for(int r=2;r=n;r++)//r表示斜对角线的层数,从2取到n{for(int i=1;i=n-r+1;i++)//i表示计算第r层斜对角线上第i行元素的值{int j=i+r-1;//j表示当斜对角线层数为r,行下标为i时的列下标m[i][j]=m[i+1][j]+p[i-1]*p[i]*p[j];//计算当断开位置为i时对应的数乘次数s[i][j]=i;//断开位置为ifor (int k=i+1;kj;k++){int t=m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j];/*计算断开位置k为从i到j(不包括i和j)的所有取值对应的(Ai*.....*Ak)*(Ak+1*.....Aj)的数乘次数*/if(tm[i][j]){m[i][j]=t;//将Ai*....Aj的最少

文档评论(0)

WHpwbonDHw + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档