现在的位置: 首页 > 综合 > 正文

c语言矩阵相乘

2018年05月09日 ⁄ 综合 ⁄ 共 1261字 ⁄ 字号 评论关闭
#include<stdio.h>
int main()
{
float a[100][100],b[100][100],c[100][100];//定义三个数组,分别存储矩阵A,B,C
int m1,n1,m2,n2,i1,j1,i2,j2,i3,j3,i4,j4,k;
float s[100][100]={0};//赋值使数组s元素初值全部为零
printf("请输入矩阵A行数m1,列数n1:");//输入矩阵A行数,列数
scanf("%d,%d",&m1,&n1);
printf("请输入矩阵B行数m2,列数n2:");//输入矩阵B行数,列数
scanf("%d,%d",&m2,&n2);
printf("\n\n");//如果不可以相乘,下面将出现判断,在此换行,便于观看
if(n1!=m2)
printf("不可以相乘!!!");//判断是否可以相乘
printf("\n\n");
if((m1>100)||(n1>100))
printf("数目过多!!!");//控制矩阵A元素数量在数组容纳范围内
else
{
for(i1=1;i1<=m1;i1++)
{
for(j1=1;j1<=n1;j1++)
{
printf("a[%d][%d]=:",i1,j1);
scanf("%f",&a[i1-1][j1-1]);//输入矩阵A元素
}
}
}
printf("\n");//分隔开A,B的元素输入,便于观看
if((m2>100)||(n2>100))
printf("数目过多!!!");
else
{
for(i2=1;i2<=m2;i2++)
{
for(j2=1;j2<=n2;j2++)
{
printf("b[%d][%d]=:",i2,j2);
scanf("%f",&b[i2-1][j2-1]);//输入矩阵B元素
}
}
}
printf("矩阵A:\n");//输出矩阵A,便于观看,检验
for(i3=0;i3<m1;i3++)
{
for(j3=0;j3<n1;j3++)
{
printf("%f",a[i3][j3]);
}
printf("\n");
}
printf("\n");//与矩阵B的输出结果隔开,便于观看
printf("矩阵B:\n");//输出矩阵A,便于观看,检验
for(i4=0;i4<m2;i4++)
{
for(j4=0;j4<n2;j4++)
{
printf("%f",b[i4][j4]);
}
printf("\n");
}
printf("\n");
printf("矩阵C=A*B:\n");
for(i4=0;i4<m1;i4++)
{
for(j4=0;j4<n2;j4++)
{
for(k=0;k<n1;k++)
{
s[i4][j4]=s[i4][j4]+a[i4][k]*b[k][j4];//定义矩阵乘法,相乘时,有一个指标是一样的,都用k
}
c[i4][j4]=s[i4][j4];//定义矩阵乘法
printf("%f",c[i4][j4]);
if(j4==n2)
printf("\n");//控制在列指标到达N时换行
}
}
return0;
}

抱歉!评论已关闭.