稀疏矩阵的运算完美版Word下载.docx
《稀疏矩阵的运算完美版Word下载.docx》由会员分享,可在线阅读,更多相关《稀疏矩阵的运算完美版Word下载.docx(15页珍藏版)》请在冰豆网上搜索。
指导教师签名
日期
备注
评分等级有五种:
优秀、良好、中等、及格、不及格
附件:
稀疏矩阵的转换
一、课题内容和要求
1.问题描述
设计程序用十字链表实现稀疏矩阵的加、减、乘、转置。
2.需求分析
(1)设计函数建立稀疏矩阵,初始化值。
(2)设计函数输出稀疏矩阵的值。
(3)构造函数进行两个稀疏矩阵相加,输出最终的稀疏矩阵。
(4)构造函数进行两个稀疏矩阵相减,输出最终的稀疏矩阵。
(5)构造函数进行两个稀疏矩阵的相乘,输出最终的稀疏矩阵。
(6)构造函数进行稀疏矩阵的转置,并输出结果。
(7)退出系统。
二、设计思路分析
三、概要设计
为了实现以上功能,可以从3个方面着手设计。
1.主界面设计
为了实现对稀疏矩阵的多种算法功能的管理,首先设计一个含有多个菜单项的主控菜单子程序以链接系统的各项子功能,方便用户交互式使用本系统。
本系统主控菜单运行界面如图所示。
2.存储结构设计
本系统采用单链表结构存储稀疏矩阵的具体信息。
其中:
全部结点的信息用头结点为指针数组的单链表存储。
3.系统功能设计
本系统除了要完成稀疏矩阵的初始化功能外还设置了4个子功能菜单。
稀疏矩阵的初始化由函数itypedefintElemType
实现。
建立稀疏矩阵用voidCreat()实现,依据读入的行数和列数以及非零元素的个数,分别设定每个非零元素的信息。
4个子功能的设计描述如下。
(1)稀疏矩阵的加法:
此功能由函数voidXiangjia()实现,当用户选择该功能,系统即提示用户初始化要进行加法的两个矩阵的信息。
然后进行加法,最后输出结果。
(2)稀疏矩阵的乘法:
此功能由函数voidXiangcheng()实现。
当用户选择该功能,系统提示输入要进行相乘的两个矩阵的详细信息。
然后进行相乘,最后得到结果。
(3)稀疏矩阵的转置:
此功能由函数voidZhuanzhi()实现。
当用户选择该功能,系统提示用户初始化一个矩阵,然后进行转置,最终输出结果。
(4)退出:
即退出稀疏矩阵的应用系统。
由函数5实现,但用户选择此功能时,系统会提示你是否确实想退出,如果是,则退出,否则继续。
三、
模块设计
1.模块设计
本程序包含1个模块:
主程序模块加各功能实现模块。
2.系统子程序及功能设计
本系统共设置7个子程序,各子程序的函数名及功能说明如下。
(1)typedefintElemType
//初始化矩阵
(2)voidCreat(TSMatrix&
M)
//建立矩阵
(3)voidPrint_SMatrix(TSMatrixM)
//输出矩阵的信息
以下编号(4)-(6)是稀疏矩阵的基本操作。
依次是:
相加,相乘,转置等。
(4)voidXiangjia(TSMatrixA,TSMatrixB,TSMatrix&
C,intn)
//把A和B两个矩阵相加,结果是C
(5)voidXiangcheng(TSMatrixA,TSMatrixB,TSMatrix&
Q)
//把A和B两个矩阵相乘,结果是Q
(6)voidZhuanzhi(TSMatrix*a,TSMatrix*b)
//把A转置
(7)voidmain()
//主函数。
设定界面的颜色,大小和窗口的标题,调用工作区模块函数
四、详细设计
#include<
stdio.h>
malloc.h>
stdlib.h>
#defineMAXSIZE40//假设非零元素个数的最大值为40
#defineMAXRC20//假设矩阵的最大行数为20
typedefintElemType;
typedefstruct
{
inti,j;
//非零元的行下标和列下标
ElemTypee;
//非零元的值
}Triple;
Tripledata[MAXSIZE+1];
intrpos[MAXRC+1];
//各行第一个非零元在三元组的位置表
inths,ls,fls;
}TSMatrix,*Matrix;
voidCreat(TSMatrix&
M)
inti,k;
for(i=1;
i<
=MAXRC+1;
i++)
M.rpos[i]=0;
printf("
请输入矩阵的行数、列数和非零元个数(以空格隔开):
"
);
scanf("
%d%d%d"
&
M.hs,&
M.ls,&
M.fls);
=M.fls;
i++)
{
请用三元组形式输入矩阵的元素(行列非零元素):
M.data[i].i,&
M.data[i].j,&
M.data[i].e);
}
for(i=1,k=1;
=M.hs;
{
M.rpos[i]=k;
while(M.data[k].i<
=i&
&
k<
=M.fls)
k++;
}
voidXiangjia(TSMatrixA,TSMatrixB,TSMatrix&
inta,b,temp,l;
C.hs=A.hs;
C.ls=A.ls;
a=b=l=1;
while(a<
=A.fls&
b<
=B.fls)
if(A.data[a].i==B.data[b].i)
if(A.data[a].j<
B.data[b].j)
C.data[l++]=A.data[a++];
elseif(A.data[a].j>
{C.data[l]=B.data[b];
C.data[l++].e=n*B.data[b++].e;
}
else{
temp=A.data[a].e+n*B.data[b].e;
if(temp)
C.data[l]=A.data[a];
C.data[l].e=temp;
l++;
a++;
b++;
elseif(A.data[a].i<
B.data[b].i)
else{C.data[l]=B.data[b];
=A.fls)
while(b<
C.fls=l-1;
intXiangcheng(TSMatrixA,TSMatrixB,TSMatrix&
intarow,brow,ccol,tp,p,q,t;
intctemp[MAXRC+1];
if(A.ls!
=B.hs)return0;
Q.hs=A.hs;
Q.ls=B.ls;
Q.fls=0;
if(A.fls*B.fls)
for(arow=1;
arow<
=A.hs;
arow++)
for(ccol=1;
ccol<
=Q.ls;
ccol++)
ctemp[ccol]=0;
Q.rpos[arow]=Q.fls+1;
if(arow<
A.hs)tp=A.rpos[arow+1];
elsetp=A.fls+1;
for(p=A.rpos[arow];
p<
tp;
p++)
brow=A.data[p].j;
if(brow<
B.hs)t=B.rpos[brow+1];
elset=B.fls+1;
for(q=B.rpos[brow];
q<
t;
q++)
ccol=B.data[q].j;
ctemp[ccol]+=A.data[p].e*B.data[q].e;
}
ccol++)
if(ctemp[ccol]){
if(++Q.fls>
MAXSIZE)return0;
Q.data[Q.fls].i=arow;
Q.data[Q.fls].j=ccol;
Q.data[Q.fls].e=ctemp[ccol];
return1;
voidPrint_SMatrix(TSMatrixM)
intk,l,n;
Matrixp;
p=&
M;
for(k=1,n=1;
k<
=p->
hs;
k++)
for(l=1;
l<
ls;
l++)
if(p->
data[n].i==k&
p->
data[n].j==l)
%5d"
p->
data[n].e);
n++;
else
0);
\n"
voidZhuanzhi(TSMatrix*a,TSMatrix*b)
{
intq,col,p;
b->
hs=a->
ls=a->
fls=a->
fls;
if(b->
fls)
q=1;
for(col=1;
col<
=a->
col++)
for(p=1;
if(a->
data[p].j==col)
data[q].i=a->
data[p].j;
data[q].j=a->
data[p].i;
data[q].e=a->
data[p].e;
++q;
voidDestory_SMatrix(TSMatrix&
M)
M.hs=M.ls=M.fls=0;
voidmain()
TSMatrixA,B,C;
TSMatrix*p=&
A,*q=&
B;
intflag,n;
while
(1)
system("
cls"
\n\n\n"
\t┏━━━━━━━━━━━━━━━━━━━━━━━━━━━┓\n"
\t┃***稀疏矩阵的加、减、转、乘***┃\n"
\t┣━━━━━━━━━━━━━━━━━━━━━━━━━━━┫\n"
\t┃1、稀疏矩阵的加法┃\n"
\t┃2、稀疏矩阵的减法┃\n"
\t┃3、稀疏矩阵的转置┃\n"
\t┃4、稀疏矩阵的乘法┃\n"
\t┃5、退出该应用程序┃\n"
\t┗━━━━━━━━━━━━━━━━━━━━━━━━━━━┛\n"
输入要进行的项目的编号:
%d"
flag);
if(flag==5)break;
Creat(A);
矩阵A:
Print_SMatrix(A);
switch(flag)
case1:
Creat(B);
n=1;
矩阵B:
Print_SMatrix(B);
if(A.hs==B.hs&
A.ls==B.ls)
A+B:
Xiangjia(A,B,C,n);
Print_SMatrix(C);
elseprintf("
错误!
行列不一致\n"
break;
case2:
n=-1;
A-B:
case3:
A->
B:
Zhuanzhi(p,q);
case4:
A*B:
n=Xiangcheng(A,B,C);
if(!
n)printf("
行列不匹配\n"
elsePrint_SMatrix(C);
default:
输入错误!
Destory_SMatrix(A);
Destory_SMatrix(B);
Destory_SMatrix(C);
getchar();
getchar();
\n\t\t\t***程序已经退出***\n"
五、测试数据及其结果分析
六、调试过程中的问题
在进行intXiangcheng(TSMatrixA,TSMatrixB,TSMatrix&
Q)函数调用的时候进行运算的时候出现了一点小差错。
七、专业课程设计总结
总通过本次实验,在我的学习理解上有很大的促进作用。
单单对稀疏矩阵,我了解了由于矩阵在程序中常使用二维阵列表示,二维阵列的大小稀疏矩阵与使用的存储器空间成正比,如果多数的元素没有数据,则会造成存储器空间的浪费,为此,必须设计稀疏矩阵的阵列储存方式,利用较少的存储器空间储存完整的矩阵数据。
二维数组Amn中有N个非零元素,若N<
<
m*n,则称A为稀疏矩阵。
由于稀疏矩阵中含有很多的0元素,在计算机中存储会浪费很多的空间,因此我们通常采用压缩存储的方法。
其他要求:
1、最后阶段请认真完成课程设计报告的电子文档和封面(见专业课程设计I报告封面.doc)。
2、每个学生最后提交的材料:
①专业课程设计I报告的电子文档(起名为:
专业课程设计I报告_学号.doc),单独作为一个文件,验收时拷到老师U盘。
②软件压缩文件一个,专业课程设计I_源代码_学号.rar,验收时拷到老师U盘。