数据结构课程设计41号李家福.doc

数据结构课程设计 稀疏矩阵运算器 班级：计算机应用10-1班 学号： 10602101041 姓名： 李 家 福 2011-12-27 1、设计综述 用三元组实现稀疏矩阵的转置、相加、相乘 （1）用creat函数创建三元组； （2）用print函数打印计算结果； （3）用add函数实现稀疏矩阵的相加运算； （4）用mult 函数实现稀疏矩阵的相乘运算； （5）用?Transm函数实现稀疏矩阵的转置运算； （6）用menu函数创建一个菜单 以顺序存储结构来表示三元组表，则可得到稀疏矩阵的一种压缩存储方式——三元组顺序表。 //-----------稀疏矩阵的三元组顺序存储表示------------- #define MAXSIZE 100 /*假设非零元个数的最大值为20*/ typedef struct {int i,j; /*该非零元的行下标和列下标*/ int v; }Triple; typedef struct {Triple data[MAXSIZE+1]; /*非零元三元组表,data[0]未用*/ int mu,nu,tu; /*矩阵的行数,列数和非零元个数*/ }TSMastrix; 在此，data域中表示非零元得三元组是以行序为主序顺序排列的，这样有利于进行某些矩阵运算。 2、算法实现描述 算法的设计思想? a) 相加运算 对于两个稀疏矩阵相加，即行与行，列与列相加 b）相乘运算 若设Q=M*N 其中M是m1*n1矩阵，N是m2*n2矩阵，只有当n1=m2时才可以相乘。 乘积矩阵Q中元素 Q（i，j）=∑M（i，k）*N（k，j） 1≤i≤m1，1≤j≤n2 在算法中，不论M（i，k）和 N（k，j）的值是否为零，都要进行一次乘法运算，而实际上，这两者有一个值为零时，其乘积也为零。因此，在对稀疏矩阵进行运算时，应免去这种无效操作，只需在M.data和N.data中找到相应的各对元素(即M.data中的j值和 N.data中的i值相等的各对元素)相乘即可。 ?c) 转置运算 对于一个m*n的矩阵M，它的转置矩阵T是一个n*m的矩阵，且T（i，j）=M（j，i）， 1≤i≤n，1≤j≤m。完成一个稀疏矩阵的转置分为三步： （1）将矩阵的行列值相互交换； （2）将每个三元组中的i和j相互调换； （3）重排三元组之间的次序便可实现矩阵的转置； 稀疏矩阵各种运算的性质变换 a)加法运算 两个稀疏矩阵的加和矩阵仍然是稀疏矩阵 b）乘法运算 两个稀疏矩阵的乘积矩阵不是稀疏矩阵 c）转置运算 一个稀疏矩阵的转置矩阵仍然是稀疏矩阵 ????时间复杂度的分析 a)加法运算 由于两矩阵行列相等，故时间复杂度为O（行*列） b）乘法运算 设M是m1*n1矩阵，N是m2*n2矩阵，则时间复杂度是O（m1*n1*n2） c）转置运算 时间复杂度和原矩阵的列数及非零元的个数的乘积成正比，即O（mu*nu） 3、总体设计（10%） 4、程序源代码 #includestdio.h #includestdlib.h #includeconio.h #includemalloc.h #define MAXSIZE 100 /*假设非零元个数的最大值为20*/ typedef struct {int i,j; /*该非零元的行下标和列下标*/ int v; }Triple; typedef struct {Triple data[MAXSIZE+1]; /*非零元三元组表,data[0]未用*/ int mu,nu,tu; /*矩阵的行数,列数和非零元个数*/ }juzheng; void creat(juzheng T)/*创建三元组*/ { int k,i,j; printf(\n 请输入矩阵\n); do{ printf(\n请输入矩阵的行数，列数，非零元素个数\n); scanf(%d%d%d,T.mu,T.nu,T.tu); if((T.mu0||T.mu20)||(T.nu0||T.nu20)||((T.tu/(T.mu*T.nu))0.05)||((T.tuT.mu*T.nu)||T.tuMAXSIZE)) printf(\n 不满足，请重新输入!!\n); }while((T.mu0||T.mu20)||(T.nu0||T.nu20)||((T.tu/(T.mu*T.nu))