稀疏矩阵的运算(完美版)

专业课程设计I报告（2011 / 2012 学年第二学期）题目稀疏矩阵的转换专业软件工程学生姓名张鹏宇班级学号 09003018指导教师张卫丰指导单位计算机学院软件工程系日期 2012年6月18号指导教师成绩评定表附件：稀疏矩阵的转换一、课题内容和要求1．问题描述设计程序用十字链表实现稀疏矩阵的加、减、乘、转置。

2．需求分析（1）设计函数建立稀疏矩阵，初始化值。

（2）设计函数输出稀疏矩阵的值。

（3）构造函数进行两个稀疏矩阵相加，输出最终的稀疏矩阵。

（4）构造函数进行两个稀疏矩阵相减，输出最终的稀疏矩阵。

（5）构造函数进行两个稀疏矩阵的相乘，输出最终的稀疏矩阵。

（6）构造函数进行稀疏矩阵的转置，并输出结果。

（7）退出系统。

二、设计思路分析（1）设计函数建立稀疏矩阵，初始化值。

（2）设计函数输出稀疏矩阵的值。

（3）构造函数进行两个稀疏矩阵相加，输出最终的稀疏矩阵。

（4）构造函数进行两个稀疏矩阵相减，输出最终的稀疏矩阵。

（5）构造函数进行两个稀疏矩阵的相乘，输出最终的稀疏矩阵。

（6）构造函数进行稀疏矩阵的转置，并输出结果。

（7）退出系统。

三、概要设计为了实现以上功能，可以从3个方面着手设计。

1．主界面设计为了实现对稀疏矩阵的多种算法功能的管理，首先设计一个含有多个菜单项的主控菜单子程序以链接系统的各项子功能，方便用户交互式使用本系统。

本系统主控菜单运行界面如图所示。

2．存储结构设计本系统采用单链表结构存储稀疏矩阵的具体信息。

其中：全部结点的信息用头结点为指针数组的单链表存储。

3．系统功能设计本系统除了要完成稀疏矩阵的初始化功能外还设置了4个子功能菜单。

稀疏矩阵的初始化由函数i typedef int ElemType 实现。

建立稀疏矩阵用void Creat()实现，依据读入的行数和列数以及非零元素的个数，分别设定每个非零元素的信息。

4个子功能的设计描述如下。

（1）稀疏矩阵的加法：此功能由函数void Xiangjia( )实现，当用户选择该功能，系统即提示用户初始化要进行加法的两个矩阵的信息。

然后进行加法，最后输出结果。

（2）稀疏矩阵的乘法：此功能由函数void Xiangcheng( )实现。

当用户选择该功能，系统提示输入要进行相乘的两个矩阵的详细信息。

然后进行相乘，最后得到结果。

（3）稀疏矩阵的转置：此功能由函数void Zhuanzhi( )实现。

当用户选择该功能，系统提示用户初始化一个矩阵，然后进行转置，最终输出结果。

（4）退出：即退出稀疏矩阵的应用系统。

由函数5实现，但用户选择此功能时，系统会提示你是否确实想退出，如果是，则退出，否则继续。

三、模块设计1.模块设计本程序包含1个模块：主程序模块加各功能实现模块。

2．系统子程序及功能设计本系统共设置7个子程序，各子程序的函数名及功能说明如下。

（1）typedef int ElemType ,&[i].j,&[i].e);}for(i=1,k=1;i<=;i++){[i]=k;while[k].i<=i && k<=k++;}}void Xiangjia(TSMatrix A,TSMatrix B,TSMatrix &C,int n){int a,b,temp,l;=;=;a=b=l=1;while(a<= && b<={if[a].i==[b].i){if[a].j<[b].j)[l++]=[a++];else if[a].j>[b].j){[l]=[b]; [l++].e=n*[b++].e;}else{temp=[a].e+n*[b].e;if(temp){[l]=[a];[l].e=temp;l++;}a++;b++;}}else if[a].i<[b].i)[l++]=[a++];else {[l]=[b]; [l++].e=n*[b++].e;} }while(a<=[l++]=[a++];while(b<={[l]=[b]; [l++].e=n*[b++].e;}=l-1;}int Xiangcheng(TSMatrix A,TSMatrix B,TSMatrix &Q) {int arow,brow,ccol,tp,p,q,t;int ctemp[MAXRC+1];if!= return 0;=;=;=0;if*{for(arow=1;arow<=;arow++){for(ccol=1;ccol<=;ccol++)ctemp[ccol]=0;[arow]=+1;if(arow< tp=[arow+1];else tp=+1;for(p=[arow];p<tp;p++){brow=[p].j;if(brow< t=[brow+1];else t=+1;for(q=[brow];q<t;q++){ccol=[q].j;ctemp[ccol]+=[p].e*[q].e;}}for(ccol=1;ccol<=;ccol++){if(ctemp[ccol]){if(++>MAXSIZE) r eturn 0;[].i=arow;[].j=ccol;[].e=ctemp[ccol];}}}}return 1;}void Print_SMatrix(TSMatrix M){int k,l,n;Matrix p;p=&M;for(k=1,n=1;k<=p->hs;k++){for(l=1;l<=p->ls;l++){if(p->data[n].i==k && p->data[n].j==l){printf("%5d",p->data[n].e);n++;}elseprintf("%5d",0);}printf("\n");}printf("\n");}void Zhuanzhi(TSMatrix *a,TSMatrix *b){int q,col,p;b->hs=a->ls;b->ls=a->hs;b->fls=a->fls;if(b->fls){q=1;for(col=1;col<=a->ls;col++)for(p=1;p<=a->fls;p++)if(a->data[p].j==col){b->data[q].i=a->data[p].j;b->data[q].j=a->data[p].i;b->data[q].e=a->data[p].e;++q;}}}void Destory_SMatrix(TSMatrix &M){===0;}void main(){TSMatrix A,B,C;TSMatrix *p=&A,*q=&B;int flag,n;while(1){system("cls");printf("\n\n\n");printf("\t┏━━━━━━━━━━━━━━━━━━━━━━━━━━━┓\n");printf("\t┃ *** 稀疏矩阵的加、减、转、乘 *** ┃\n");printf("\t┣━━━━━━━━━━━━━━━━━━━━━━━━━━━┫\n");printf("\t┃1、稀疏矩阵的加法┃\n");printf("\t┃2、稀疏矩阵的减法┃\n");printf("\t┃3、稀疏矩阵的转置┃\n");printf("\t┃4、稀疏矩阵的乘法┃\n");printf("\t┃5、退出该应用程序┃\n");printf("\t┗━━━━━━━━━━━━━━━━━━━━━━━━━━━┛\n");printf("输入要进行的项目的编号:");scanf("%d",&flag);if(flag==5) break;Creat(A);printf("矩阵A:\n"); Print_SMatrix(A);switch(flag){case 1: C reat(B);n=1;printf("矩阵B:\n");Print_SMatrix(B);if== && =={printf("A+B:\n");Xiangjia(A,B,C,n);Print_SMatrix(C);}else printf("错误!行列不一致\n");break;case 2: Creat(B);n=-1;printf("矩阵B:\n");Print_SMatrix(B);if== && =={printf("A-B:\n");Xiangjia(A,B,C,n);Print_SMatrix(C);}else printf("错误!行列不一致\n");break;case 3: printf("A->B:\n");Zhuanzhi(p,q);Print_SMatrix(B);break;case 4: C reat(B);printf("矩阵B:\n");Print_SMatrix(B);printf("A*B:\n");n=Xiangcheng(A,B,C);if(!n) printf("错误!行列不匹配\n");else Print_SMatrix(C);break;default: printf("输入错误!\n");}Destory_SMatrix(A);Destory_SMatrix(B);Destory_SMatrix(C);getchar();getchar();}printf("\n\t\t\t ***程序已经退出***\n");getchar();}五、测试数据及其结果分析六、调试过程中的问题在进行int Xiangcheng(TSMatrix A,TSMatrix B,TSMatrix &Q) 函数调用的时候进行运算的时候出现了一点小差错。