矩阵的运算程序C语言版

1、#include<stdio.h>#defineM20#defineN20floatAMN;floatBMN;floatCMN;inti,j,m,n,p,q;floaty=1.0;voidmain()(printf("#'n”);printf("#欢迎您使用矩阵函数包系统。print*”#系统功能：printf("#计算两个矩阵相加，相减；printf("#数乘矩阵；矩阵转置；求逆#榭#n");矩printf("#阵制用的M”)；printf("#、n”);printf("请选择您需要的运算，按

2、回车键结束n”)；printf*n").printf("1,矩阵相加n”);,printf("2矩阵相减n")；,printf("3数乘矩阵n")；,printf("4矩阵转置n");.printf("5逆矩阵n”)；printf("6,退出n”)；print%'intx;seanf("%d",&x);switch(x)(case 1: printf("请输入矩阵A的行数和列数(用逗号隔开):")；seanf("%d,%d&quo

3、t;,&i,&j);printf("请输入矩阵B的行数和列数(用逗号隔开):")；seanf(”d,%d”,&m,&n);if(i!=m|j!=n)printf("*对不起，您输入的两个矩阵不能相加，请重试.*S”)；elseprintf("请输入矩阵A:fi");for(p=0;p<i;p+)for(q=0;qvj;q+)seanf("%f',&Apq);printf("输出矩阵A：n");for(p=0;p<i;p+)for(q=0;q<j;q+

4、)(printf("%10.2r,Apq);if(q+1)%j=0)printf("n");printf("请输入矩阵B：n");for(p=0;p<i;p+)for(q=0;q<j;q+)scanf(“r,&Bpq);printf("输出矩阵B：n");for(p=0;p<i;p+)for(q=0;q<j;q+)pnntf(H%10.2r,Bpq);if(q+1)%j=0)printf(nnH);>printf(H矩阵A+矩阵B为:n“);计算两个矩阵相加for(p=0;p<i;

5、p+)for(q=0;q<j;q+)Cpq=Apq+Bpq;for(p=0;p<i;p+)for(q=0;q<j;q+)printf("%10.2f',Cpq);if(q+1)%j=0)printf("n");>break;case 2:(printf("请输入矩阵A的行数和列数(用逗号隔开):")；scanf("%d,%d",&i,&j);printf("请输入矩阵B的行数和列数(用逗号隔开):")；scanf("%d,%d",&

6、;m,&n);if(i!=m|j!=n)printf("*对不起，您输入的两个矩阵不能相减，请重试elseprintf("请输入矩阵A:n");for(p=0;p<i;p+)for(q=0;q<j;q+)scanf("%f',&Apq);printf("输出矩阵A:n");for(p=0;p<i;p+)for(q=0;q<j;q+)(printf("%10.2f,Apq);if(q+1)%j=0)printf("n");printf("请输入矩阵B

7、：n");for(p=0;p<i;p+)for(q=0;q<j;q+)scanf("%f',&Bpq);printf("输出第矩阵B:n");for(p=0;p<i;p+)for(q=0;q<j;q+)(printf("%10.2f,Bpq);if(q+1)%j=0)printf("n");>printf("矩阵A-矩阵B为：n");计算两个矩阵相减for(p=0;p<i;p+)for(q=0;q<j;q+)Cpq=Apq-Bpq;for(p=0;

8、p<i;p+)for(q=0;q<j;q+)<printf("%10.2f',Cpq);if(q+1)%j=0)printf("n");break;case 3:floatk;printf("请输入矩阵A的行数和列数(用逗号隔开):")；scanf("%d,%d",&i,&j);printf("请输入矩阵An");for(p=0;p<i;p+)for(q=0;q<j;q+)scanf("%f',&Apq);printf(&qu

9、ot;输出矩阵An");for(p=0;p<i;p+)for(q=0;qvj;q+)(printf("%10.2f',Apq);if(q+1)%j=O)printf("n");)printf("请输入一个实数:n”)；seanf("%f',&k);for(p=0;p<i;p+)/数乘矩阵for(q=0;q<j;q+)Bpq=k*Apq;printf("输出k乘矩阵A的结果山”)；for(p=0;p<i;p+)for(q=0;q<j;q+)(printf("%10

10、.2f',Bpq);if(q+1)%j=O)printf("n"););break;case 4: printfd请输入矩阵A的行数和列数(用逗号隔开)seanf(“d,%d”，&printf("请输入矩阵A:n");for(p=0;p<i;p+)for(q=0;q<j;q+)seanf("%f',&Apq);printf(输出矩阵AnM);for(p=0;p<i;p+)for(q=0;q<j;q+)printf(H%10.2f',Apq);if(q+1)%j=O)printf(H

11、nM);)for(p=0;p<i;p+)/转置for(q=0;q<j;q+)Bqp=Apq;printf("输出矩阵A的转置矩阵:n");for(p=0;p<j;p+)for(q=0;q<i;q+)(printfr%10.2r,Bpq);if(q+1)%i=0)printfCAn'1);/;break;case 5:floataM2*M;floatbN2wN;floatt,x;intk,T;printf("输入方阵的维数：n");请输入方阵，即行和列相等的矩阵。scanf(H%dH,&T);printf("

12、;请输入矩阵:n");for(i=0;i<T;i+)for(j=O;j<T;j+)scanf("%f',&biO)；printf(H原矩阵为：nH);for(i=0;i<T;i+)(for(j=O;j<T;j+)printf(M%1O.3fbi0)；printfCAn'1);/for(i=0;i<T;i+)for(j=0;j<(2*T);j+)ifG<T)ai0=bi0；elseif(j=T+i)ai0=1.O;elseai0=O.O;for(i=0;i<T;i+)(for(k=0;k<T;k+)(if(k!=i)t=aki/aii;for(j=0;jv(2T);j+)x=aij*t;akO=akj-x;for(i=0;i<T;i+)4=aii;ford=0;j<(2*T);j+)aiO=aiO/t;for(i=0;i<T;i+)y=y)皿i；if(y=O)printf("对不起，您输入的矩阵没有逆矩阵，请重新输入。n")；else(for(i=0;i<T;i+)forQ=0;j<T;j+)bij=aiU+T;printf("逆矩阵为:n");f