课程设计矩阵运算系统

wen滨江学院windows 程序设计综合实验课程设计题目矩阵综合运算系统学生姓名晏文涛学号***********院系电子工程系专业信息工程指导教师方忠进二Ｏ一二年12 月16 日设计了一个矩阵运算系统，该矩阵运算系统具有普通矩阵相加、相减、相乘及稀疏矩阵转置等功能。

本运算系统以Microsoft Visual C++ 6.0 作为系统开发工具，采用算数表达式处理算法来实现了矩阵的加、减、乘等混合运算和稀疏矩阵的转置矩阵运算。

系统操作简单，界面清晰，便于用户使用。

关键词：普通矩阵; 运算; VC6.01 课题描述 (1)2 设计过程 (1)3 程序编码 (3)4 测试 (10)总结 (12)参考文献 (13)1 课题描述矩阵运算系统是一个非常重要的运算,很多软件开发公司都开发了这个运算系统。

现在我们用C 语言编出这个运算系统。

它的原理是对于输入的矩阵，进行相加、相乘以及相减。

另外一个是稀疏矩阵的转置运算系统，按提示输入数值即可得到所要求的稀疏矩阵的转置矩阵。

运行环境：Visual C++ 6.02 设计过程经过对程序设计题目的分析可知，整个程序的设计实现大致分为四个模块，其中每一个模块对应一个函数，他们的功能分别是：1）矩阵相加运算函数(ADD)，主要实现将两矩阵相加的功能；2）矩阵相乘运算函数(MUL)，主要实现将两矩阵相乘的功能；3)矩阵相减函数（SNB）;实现的功能是矩阵之间的减法4）稀疏矩阵矩阵转置函数(TRANPOSE) 实现的功能是将稀疏矩阵进行转置。

在这些函数当中，第1、2、4个函数的实现严格按照题目的要求，而第3个函数为自行设计的函数。

程序的一次运行当中可以循环执行所有的功能，并根据需要终止程序的执行。

在这个程序中，将各个功能以子程序模块的形式编写。

这样使所编写的程序简单明了，逻辑性思维表达明确，具有很强的可读性。

流程图如下：1）矩阵相乘流程图如图2.1所示：图2.12）矩阵相加流程图如图2.2所示图2.23）矩阵相减流程图如图2.3所示图2.34）矩阵相减、相加、相乘MAIN函数如图2.4所示图2.45）稀疏矩阵的转置流程图如图2.5所示图2.53 程序编码/*矩阵相加、相减以及相乘*/#include "stdio.h"#include "stdlib.h"#define M 3#define N 3//指针数组int a[M][N];int b[M][N];int c[M][N];void matrixMul(int b[][M],int c[][N]);void matrixAdd(int b[][N],int c[][N]);void matrixSub(int b[][N],int c[][N]);int main(){int i,j,temp=0;printf("Please input int matrix b[%d][%d]\n",M,N);for(i=0;i<N;i++)for(j=0;j<N;j++){scanf("%d",&temp);b[i][j]=temp;}printf("Please input int matrix c[%d][%d]\n",M,N);for(i=0;i<M;i++){for(j=0;j<N;j++){scanf("%d",&temp);c[i][j]=temp;}}//输出原矩阵printf("Now print resource matrix b[%d][%d]=",M,N);for(i=0;i<M;i++){printf("\n");for(j=0;j<N;j++)printf("%d ",b[i][j]);}printf("\n");printf("Now print resource matrix c[%d][%d]=",M,N);for(i=0;i<M;i++){printf("\n");for(j=0;j<N;j++)printf("%d ",c[i][j]);}if (M == N){matrixMul(b,c); //调用矩阵相乘//输出矩阵相乘结果printf("\n");printf("Now printm matrixMul results matrix a[%d][%d]=B*C:",M,N);for(i=0;i<M;i++){printf("\n");for(j=0;j<N;j++)printf("%d ",a[i][j]);}}else{printf("matrix Error,check and try again ");}//输出矩阵相加结果matrixAdd(b,c);printf("\n");printf("Now printm matrixAdd results matrix a[%d][%d]=B+C:",M,N);for(i=0;i<M;i++){printf("\n");for(j=0;j<N;j++)printf("%d ",a[i][j]);}//输出矩阵相减结果matrixSub(b,c);printf("\n");printf("Now printm matrixSub results matrix a[%d][%d]=B-C:",M,N);for(i=0;i<M;i++){printf("\n");for(j=0;j<N;j++)printf("%d ",a[i][j]);}return 0;}void matrixMul(int b[][M],int c[][N]){int i,j,k;for(i=0;i<M;i++)for(j=0;j<N;j++){for(k=0;k<N;k++)a[i][j]+=b[i][k]*c[k][j];}}void matrixAdd(int b[][M],int c[][N]){int i,j;for(i=0;i<M;i++){for(j=0;j<N;j++){a[i][j]=b[i][j]+c[i][j];}}}void matrixSub(int b[][M],int c[][N]){int i,j;for(i=0;i<M;i++){for(j=0;j<N;j++){a[i][j]=b[i][j]-c[i][j];}}}/*矩阵的快速转置*/#include <stdio.h>#include <stdlib.h>#include <process.h>#define MAXSIZE 200 /*矩阵中最大非零元的个数*/typedef struct triple{int i; /*行标,本程序中从1开始的*/int j; /*列标，本程序中从1开始的*/int e; /*非零元*/}Triple; /*三元组定义*/typedef struct tabletype{int mu; /*矩阵的行数*/int nu; /*列数*/int tu; /*非零元个数*/Triple data[MAXSIZE]; /*非零元的三元组表,本程序中是从data[1]开始使用的*/}Tabletype; /*三元组线性表*//*以下为函数声明，注意和书本上的参数类型不同，我用的形参全为指针*/void CreatSMatrix(Tabletype *); /*生成矩阵*/void DestroySMatrix(Tabletype *); /*销毁矩阵*/void out_matrix(Tabletype *); /*打印矩阵*/int FastTransposeSMatrix(Tabletype *,Tabletype *); /*快速转置算法*/int main( void ) /*主函数*/{char ch;Tabletype a = {0,0,0,{0,0,0}}; /*初始化为0，便于输入数据时的无效检测*/Tabletype b; /*声明矩阵b*/while(1){printf(" @@@@@@@@@@@@@@本程序的功能是实现稀疏矩阵的快速转置@@@@@@@@@@@@@@@\n");printf(" @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@CreatSMatrix(&a);printf("The source Matrix:\n");out_matrix(&a);if(FastTransposeSMatrix(&a,&b)) /*若a不为零矩阵则转置a，存入b中*/{ printf("After TransposeSMatrix: \n");out_matrix(&b);}else{printf("The matrix is zeros:\n");out_matrix(&a);}/*以下为程序控制*/printf("Input 'q' to quit and 'c' run again:");do{if((ch = getchar()) == 'q' || ch == 'Q'){DestroySMatrix(&a);DestroySMatrix(&b);exit(0);}}while((ch!='C') && (ch!='c'));system("cls");}return 1;}void CreatSMatrix(Tabletype *a){int i;printf("请输入矩阵的行数、列数和非零元个数，用空格间隔:");scanf("%d%d%d",&(a->mu),&(a->nu),&(a->tu));for(i=1;i<= a->tu;){printf("请输入矩阵中第%d个非零元（按行标、列标、值的顺序，空格间隔）:",i);scanf("%d%d%d",&(a->data[i].i),&(a->data[i].j),&(a->data[i].e));if(a->data[i].i < 1 || a->data[i].i > a->mu || a->data[i].j < 1 || a->data[i].j > a->nu) /*下标越界*/{printf("注意：下标越界输入数据无效!\n请重新输入:行标范围：1--%d,列标范围1--%d\n",a->mu,a->nu);continue;if( ((a->data[i].i) < (a->data[i-1].i)) ||(((a->data[i].i) == (a->data[i-1].i)) && ((a->data[i].j) <= (a->data[i-1].j)))) /*非按行顺序输入*/{printf("注意：输入数据无效!\n请按照按行存储的顺序输入数据\n");continue;}i++;}}void DestroySMatrix(Tabletype *a){ /* 销毁稀疏矩阵a*/(*a).mu=0;(*a).nu=0;(*a).tu=0;}void out_matrix(Tabletype *a) /* 打印矩阵*/{int i,j,k = 1;for(i = 1 ;i <= a->mu; i++){for(j = 1; j<= a->nu; j++){ /*判断是否为非零元*/if((a->data[k].i == i)&&(a->data[k].j == j)){printf("%4d",a->data[k].e);k++;}elseprintf("%4d",0);}printf("\n");}}i nt FastTransposeSMatrix(Tabletype *a,Tabletype *b){int p,q,col;int *num;int *cpot;b->mu = a->nu; /*原矩阵的行数为新矩阵的列数，原列数为新行数，非零元个数不变*/ b->nu = a->mu;b->tu = a->tu;num=(int *)malloc((b->nu+1)*sizeof(int)); /* 生成两个辅助数组*/ cpot=(int *)malloc((b->nu+1)*sizeof(int));if(b->tu) /*若a不为零矩阵*/{for(col = 0;col <a->nu;col++) /*初始化矩阵a的每列中非零元的个数均为0*/num[col] = 0;for(col = 0; col <=a->tu ; col++)/*统计每列中非零元的个数*/num[a->data[col].j]++;cpot[1] = 1; /*确定每列中第一个非零元的位置*/for(col = 2;col <= a->nu; col++)cpot[col] = num[col-1]+cpot[col-1];for(p = 1; p <= a->tu; p++) /*p为a-data的下标*/{col = a->data[p].j; /*交换元素*/q = cpot[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++;cpot[col]++;}free(num); /*释放两个辅助数组*/free(cpot);return 1; /*转置成功*/}else /*a为零矩阵*/return 0;}4 测试1）矩阵相乘的运行结果如图4.1所示图4.1 矩阵乘法的运行结果截图2）矩阵相加的运行结果如图4.2所示：图4.2 矩阵加法的运行结果截图3）矩阵减法的运行结果如图4.3所示图4.3 矩阵减法的运行结果截图4）稀疏矩阵转置的运行结果如图4.4所示:图4.4图4.4 稀疏矩阵转置的运行结果截图5.总结矩阵的相加相乘、以及稀疏矩阵的转置的开发与实践是我在大学的第一次课程设计，总体上来说在这次的课程设计当中我本人受益匪浅我了解到了如何利用C语言来开发与实现一些我们常遇到的问题，总之这次课程设计让我知道了我们学习C语言要同平时常遇的问题结合起来这样才能更好的掌握所学的知识，学以致用！参考文献[1] 谭浩强.C程序设计教程[M].北京:清华大学出版社,2007[2]《数据结构》教材[M].高等教育出版社[3] 严蔚敏，吴伟民.数据结构（C语言版）[M].北京:清华大学出版社,2008。