C22-矩阵的乘法-乘方运算的C程序

《C22-矩阵的乘法-乘方运算的C程序》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、//矩阵的乘法、乘方运算的C语言程序//求两个矩阵的乘积及方阵的乘方#include//getch#include//printf,scanf#include//time#include//srand,rand//============抽象数据类型矩阵的表示============//typedefstruct{//矩阵的定义double*p;//矩阵数据存放的首地址introw;//行数intcol;//列数}Matrix;//============抽象

2、数据类型矩阵的实现============//voidInitMat(Matrix&m)//初始化为空{m.p=NULL;m.row=0;m.col=0;}voidErrExit(char*s){//遇错退出并提示printf("%s!",s);getch();exit(0);}voidInitMat(Matrix&m,intr,intc){//初始化指针及行列数if(r<=0

3、

4、c<=0)ErrExit("矩阵型行列数错");m.p=newdouble[r*c];m.row=r;m.col=c;}voidRandMat(

5、Matrix&m){//给矩阵赋随机值for(inti=0;i

6、.col;j++)if(i==j)*(m.p+i*m.col+j)=1;else*(m.p+i*m.col+j)=0;}voidPrintMat(Matrixm,char*str){//输出矩阵的所有值printf("矩阵%s为：",str);for(inti=0;i

7、p;m.p=NULL;}voidMulMat(Matrix&m,Matrixa,Matrixb){//求两矩阵的积if(a.col!=b.row)ErrExit("相乘失配");//矩阵失配不能相乘InitMat(m,a.row,b.col);//申请空间ClearMat(m);//清空矩阵为全0for(inti=0;i

8、+j);}voidCopyMat(Matrixd,Matrixs){//拷贝矩阵s至dif(d.row!=s.row

9、

10、d.col!=s.col)ErrExit("拷贝扮配");for(inti=0;i

11、(x,a.row,a.col);UnitMat(x);//构造单位阵xif(n==0){InitMat(m,a.row,a.col);CopyMat(m,x);}for(unsignedi=1;i<=n;i++){if(m.p)Destroy(m);//MulMat后必须销毁mMulMat(m,x,a);CopyMat(x,m);}Destroy(x);}//下面的程序似乎更简单，但是没有清理内存voidPowMat2(Matrix&m,Matrixa,unsignedn){//求矩阵的乘方InitMat(m,a.row,a.co

12、l);UnitMat(m);for(unsignedi=1;i<=n;i++)MulMat(m,m,a);}//=====================程序测试=================//voidmain(){Matrixm,a,b,c