矩阵的运算程序(C语言版)

#include#define M 20#define N 20float AMN;float BMN;float CMN;int i,j,m,n,p,q;float y=1.0;void main()printf( #n);printf( # 欢迎您使用矩阵函数包系统。 #n);printf( # 系统功能: #n);printf( # 计算两个矩阵相加,相减; #n);printf( # 数乘矩阵；矩阵转置;求逆 #n);printf( # 矩阵 #n);printf( #n);printf(请选择您需要的运算，按回车键结束n);printf(*n);printf(1，矩阵相加n);printf(2，矩阵相减n);printf(3，数乘矩阵n);printf(4，矩阵转置n);printf(5，逆矩阵 n);printf(6，退出 n);printf(*n);int x;scanf(%d,&x);switch (x)case 1:printf(请输入矩阵A的行数和列数(用逗号隔开):);scanf(%d,%d,&i,&j);printf(请输入矩阵B的行数和列数(用逗号隔开):) ;scanf(%d,%d,&m,&n);if(i!=m|j!=n)printf(*对不起，您输入的两个矩阵不能相加，请重试.*n);else printf(请输入矩阵A:n); for(p=0;pi;p+)for(q=0;qj;q+)scanf(%f,&Apq);printf(输出矩阵A：n);for(p=0;pi;p+)for(q=0;qj;q+) printf(%10.2f,Apq);if(q+1)%j=0)printf(n); printf(请输入矩阵B：n);for(p=0;pi;p+)for(q=0;qj;q+)scanf(%f,&Bpq);printf(输出矩阵B：n);for(p=0;pi;p+)for(q=0;qj;q+) printf(%10.2f,Bpq);if(q+1)%j=0)printf(n); printf(矩阵A+矩阵B为:n); /计算两个矩阵相加for(p=0;pi;p+)for(q=0;qj;q+)Cpq=Apq+Bpq;for(p=0;pi;p+)for(q=0;qj;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,&m,&n);if(i!=m|j!=n)printf(*对不起，您输入的两个矩阵不能相减，请重试.*n);else printf(请输入矩阵A:n); for(p=0;pi;p+)for(q=0;qj;q+)scanf(%f,&Apq);printf(输出矩阵A:n);for(p=0;pi;p+)for(q=0;qj;q+) printf(%10.2f,Apq);if(q+1)%j=0)printf(n); printf(请输入矩阵B：n);for(p=0;pi;p+)for(q=0;qj;q+)scanf(%f,&Bpq);printf(输出第矩阵B:n);for(p=0;pi;p+)for(q=0;qj;q+) printf(%10.2f,Bpq);if(q+1)%j=0)printf(n); printf(矩阵A-矩阵B为：n); /计算两个矩阵相减for(p=0;pi;p+)for(q=0;qj;q+)Cpq=Apq-Bpq;for(p=0;pi;p+)for(q=0;qj;q+)printf(%10.2f,Cpq);if(q+1)%j=0)printf(n);break; case 3:float k; printf(请输入矩阵A的行数和列数(用逗号隔开):);scanf(%d,%d,&i,&j);printf(请输入矩阵An); for(p=0;pi;p+)for(q=0;qj;q+)scanf(%f,&Apq);printf(输出矩阵An);for(p=0;pi;p+)for(q=0;qj;q+) printf(%10.2f,Apq);if(q+1)%j=0)printf(n); printf(请输入一个实数:n);scanf(%f,&k);for(p=0;pi;p+) /数乘矩阵for(q=0;qj;q+)Bpq=k*Apq;printf(输出k乘矩阵A的结果n);for(p=0;pi;p+)for(q=0;qj;q+) printf(%10.2f,Bpq);if(q+1)%j=0)printf(n); ;break; case 4:printf(请输入矩阵A的行数和列数(用逗号隔开):);scanf(%d,%d,&i,&j);printf(请输入矩阵A:n);for(p=0;pi;p+)for(q=0;qj;q+)scanf(%f,&Apq);printf(输出矩阵An);for(p=0;pi;p+)for(q=0;qj;q+) printf(%10.2f,Apq);if(q+1)%j=0)printf(n); for(p=0;pi;p+) /转置for(q=0;qj;q+)Bqp=Apq;printf(输出矩阵A的转置矩阵:n);for(p=0;pj;p+)for(q=0;qi;q+) printf(%10.2f,Bpq);if(q+1)%i=0)printf(n);break; case 5: float aM2*M; float bN2*N; float t,x; int k,T; printf(输入方阵的维数：n); /请输入方阵，即行和列相等的矩阵。scanf(%d,&T);printf(请输入矩阵:n);for(i=0;iT;i+) for (j=0;jT;j+) scanf(%f,&bij); printf(原矩阵为：n); for (i=0;iT;i+) for (j=0;jT;j+) printf(%10.3f,bij); printf(n); for(i=0;iT;i+) for(j=0;j(2*T);j+) if (jT) aij=bij; else if (j=T+i) aij=1.0; else aij=0.0; for(i=0;iT;i+)for(k=0;kT;k+) if(k!=i) t=aki/aii;for(j=0;j(2*T);j+) x=aij*t;akj=akj-x;for(i=0;iT;i+) t=aii;for(j=0;j(2*T);j+)aij=aij/t;for(i=0;iT;i+) y=y*aii;if(y=0)printf(对不起，您输入的矩阵没有逆矩阵，请重新输入。n);else for(i=0;iT;i+) for(j=0;jT;j+)bij=aij+T;printf(逆矩阵为：n); for (i=0;iT;i+)