稀疏矩阵链表实现——原创.doc

大连理工大学实验报告 学院（系）： 专业： 班级： 姓 名： 学号： 组： ___ 实验时间： 实验室： 实验台： 实验名称 实验目的和要求 Given two sparse matrices, design an algorithm to accomplish multiplication of sparse matrices, and implement it in C language by using circularly linked lists. Requirements: 1) Use stream (file) to input sparse matrices instead of console I/O 2) Analyze the time complexity of your algorithm 3) Submit the document explaining your algorithm as well as the source code. 实验原理和内容 存储结构： 在循环链表中，矩阵的每个非零元素对应着一个含有五个域的结点，这五个域分别为该非零元素在稀疏矩阵中的行号、列号、元素的数值以及两个指针，其中一个指针指向同一列的下一个非零元素的结点，另一个指针指向同一行的右边一个非零元素的结点。 循环结构：将每行的非零元素结点链接成循环链表，又将每列的非零元素结点链接成循环链表，就构成了十字形的链接结构。 头结点：对于每行、每列的循环链表都有一个空表头结点hdnode[i]，以利于元素的插入和删除运算。这些空表头结点的行号、列号都标成零，以便和其它结点相区别。 整个十字链表有一个总的空表头结点node，在一般结点标以行号、列号之处，此结点是标出矩阵的行数m和列数n。有一个指针指向这个总的空表头结点。 因各列、各行的空表头结点中的行号和列号都是零，且每列空表头结点只用到向下指针，每行空表头结点只用到向右指针，故可将这组空表头结点合用。 由于数值域也没有用，可将空表头结点的数值域改为一个指针域，将各个空表头结点也链接成一个循环链表。 e.g.稀疏矩阵M及其对应的十字链表如图所示。 基本思想： 矩阵A与矩阵B相乘时，首先要保证A的列数与B的行数相等，矩阵的第i行与矩阵B的第j列相乘求和结果将存储在矩阵C的[i][j]位置上，而且第i行的第k个元素与对应的第j列的第k个元素才可以相乘。这样A中的每个非零元素将与B中每个非零元素比较上述条件是否成立，进行矩阵乘法操作。 主要仪器设备 Visual studio编程环境 实验步骤与操作方法 矩阵乘法流程图 头结点设置 //构建循环链表框架，即头结点的个数，和next指针的初始化 x=node_1-right for(i=0;inode_1-u.entry.row;i++)//矩阵A每一行遍历 last=hdnode[i]; y=node_2-right; for(j=0;jnode_2-u.entry.col;j++) //矩阵B每一列遍历，每进行一次循环 产生一个新的结点 产生C的一个元素 for(p=x-right;p!=x;p=p-right) for(q=y-down;q!=y;q=q-down) if(p-u.entry.col==q-u.entry.row)//判断是否可以对应相乘 temp-u.entry.value+=p-u.entry.value*q-u.entry.value; temp-u.entry.value!=0 设置结点指针//头结点的next free(temp) //相乘求和为0 temp结点的right y=y-u.next;//B下一列 x=x-u.next;//A下一行 其他指针设置 输入输出调用教材中mread、mwrite函数 实验数据记录和处理 MATLAB中的结果： A= (1,4) -3 (1,5) -20 (2,5) 33 B= (1,1) 3