两个矩阵相乘c语言
A. c语言实现矩阵乘法
1,若程序定义为三行三列的矩阵,并在主函数中键盘赋值。请编写函数juzhen(int a[][3],int n),函数功能是:使矩阵左下角元素中的值乘以.例如:若矩阵中的值为:
1 9 7
2 3 8
4 5 6
n=3,则矩阵值变为:
1 9 7
6 9 8
12 15 18
注:二维数组下三角元素每行最后元素的列标和该列的列号一致
#include <stdio.h>
int juzhen(int a[][3],int n)
{
int i,j;
for (i=1;i<3;i++)
for (j=0;j<=i;j++)
a[i][j]=a[i][j]*n;
}
main()
{
int i,j,a[3][3],m;
clrscr();
printf("please enter 9 numbers :\n");
for (i=0;i<3;i++)
for (j=0;j<3;j++)
scanf("%d",&a[i][j]);
printf("the old ju zhen is:\n");
for (i=0;i<3;i++)
{
for (j=0;j<3;j++)
printf("%3d",a[i][j]);
printf("\n");
}
printf("please enter a int number:\n");
scanf("%d",&m);
juzhen (a,m);
printf("the new ju zhen is:\n");
for (i=0;i<3;i++)
{
for (j=0;j<3;j++)
printf("%3d",a[i][j]);
printf("\n");
}
}
2,请编写函数juzhen,该函数的功能是:将三行四列矩阵x乘以四行三列矩阵y,结果放在三行三列矩阵中。矩阵相乘的基本方法是:矩阵Xy中行列下标分别为i,j的元素的值,是矩阵X中第i行上四个元素与矩阵Y中第j列上四个元素对应相乘的积。
X和Y矩阵原型可运行后看到!
#include <conio.h>
#include <stdio.h>
void juzhen (int a[3][4],int b[4][3],int ab[3][3])
{
int i,j,z;
for (i=0;i<3;i++)
for (j=0;j<4;j++)
for (z=0;z<3;z++)
ab[i][z]=ab[i][z]+a[i][j]*b[j][z];/*此得为要填写的内容*/
}
main()
{
int x[3][4]={1,0,1,1,2,1,0,1,1,2,0,3};
int y[4][3]={1,1,1,0,0,0,2,1,1,1,1,3};
int xy[3][3]={0},i,j;
clrscr();
juzhen (x,y,xy);
printf("x ju zhen is:\n");
for (i=0;i<3;i++)
{
for (j=0;j<4;j++)
printf("%3d",x[i][j]);
printf("\n");
}
printf("y ju zhen is:\n");
for (i=0;i<4;i++)
{
for (j=0;j<3;j++)
printf("%3d",y[i][j]);
printf("\n");
}
printf("xy ju zhen is:\n");
for (i=0;i<3;i++)
{
for (j=0;j<3;j++)
printf("%3d",xy[i][j]);
printf("\n");
}
}
由于你没有说要求,我就写了两个,这是两道比较精典的求c语言矩阵乘法程序
B. C语言中两矩阵相乘
#include <stdio.h>
int main()
{
int x,y,z;
int i,j,l;
int sum=0;
printf("请输入前一个矩阵的行数,列数与后一个矩阵的列数:\n");
scanf("%d %d %d",&x,&y,&z);
int a[x][y];
int b[y][z];
int c[x][z];
//读入矩阵
printf("请输入矩阵a:\n");
for(i=0;i<x;i++){
for(j=0;j<y;j++){
scanf("%d",&a[i][j]);
}
}
printf("请输入矩阵b:\n");
for(i=0;i<y;i++){
for(j=0;j<z;j++){
scanf("%d",&b[i][j]);
}
}
printf("矩阵a为:\n");
for(i=0;i<x;i++){
for(j=0;j<y;j++){
printf("%d",a[i][j]);
if(j==y-1)
printf("\n");
else
printf("\t");
}
}
printf("矩阵b为:\n");
for(i=0;i<y;i++){
for(j=0;j<z;j++){
printf("%d",b[i][j]);
if(j==z-1)
printf("\n");
else
printf("\t");
}
}
for(i=0;i<x;i++){
for(l=0;l<z;l++){
for(j=0;j<y;j++){
sum+=a[i][j]*b[j][l];
if(j==y-1){
c[i][l]=sum,sum=0;
}
}
}
}
printf("矩阵a与矩阵b的乘积为:\n");
for(i=0;i<x;i++){
for(j=0;j<z;j++){
printf("%d",c[i][j]);
if(j==z-1)
printf("\n");
else
printf("\t");
}
}
return 0;
}
C. C语言题,求两个矩阵的乘积
简单题,按照矩阵的定义去做即可
源代码如下
//#pragma GCC diagnostic error"-std=c11"
#include<stdlib.h> //有随机数库
#include<malloc.h>
#include<time.h> //用于产生随机数种子
#include<math.h>
#include<string.h>
#include<stdio.h>
#define ELE int
typedef struct Metrix {
ELE *A;
int row, col;
}Metrix;
//初始化矩阵
void init(Metrix *M, int row, int col, int auto_read) {
M->A = (ELE*)calloc(row*col, sizeof(ELE));
M->row = row;
M->col = col;
if (auto_read) {
int i, size = row*col;
for (i = 0; i<size; i++)
scanf("%d", M->A + i);
}
}
#define IDX(M,r,c) (r*(M->col)+c)
ELE* get(Metrix *M, int r, int c) {
return M->A + IDX(M, r, c);
}
int main()
{
int m, l, n, i, j, k;
scanf("%d%d%d", &m, &l, &n);
Metrix L, R;
init(&L, m, l, 1);
init(&R, l, n, 1);
// init(M,m,n,0);
//矩阵乘法走起
for (i = 0; i<m; i++) {
for (j = 0; j<n; j++) {
ELE *lp = get(&L, i, 0), *rp = get(&R, 0, j);
ELE out = 0;
for (k = 0; k<l; k++) {
out += *(rp)*lp[k];
rp += n; //rp换下一行
}
printf("%d ", out);
}
printf(" "); //删去空格再换行。
}
return 0;
}
D. C语言两个矩阵乘积
#include<stdio.h>
int main(){
int a[100][100],b[100][100],c[100][100]={0};
int m1,n1,m2,n2,i,j,k;
printf("请输入矩阵A行数m1 列数n1:\n");
scanf("%d%d",&m1,&n1);
printf("请输入矩阵B行数m2 列数n2:\n");
scanf("%d%d",&m2,&n2);
if(n1!=m2)
printf("\n不可以相乘\n");
else if((m1>100)||(n1>100)||(m2>100)||(n2>100))
printf("越界!!");
else{
for(i=0;i<m1;i++)
for(j=0;j<n1;j++){
printf("a[%d][%d]=:",i,j);
scanf("%d",&a[i][j]);}
printf("\n");
for(i=0;i<m2;i++)
for(j=0;j<n2;j++){
printf("b[%d][%d]=:",i,j);
scanf("%d",&b[i][j]);}
printf("\n矩阵C=A*B:\n");
for(i=0;i<m1;i++)
for(j=0;j<n2;j++){
for(k=0;k<n1;k++)
c[i][j]+=a[i][k]*b[k][j];
printf("%d",c[i][j]);
if(j==n2)
printf("\n");}}
return 0;}
E. C语言求个程序!!!!计算两个矩阵乘积
#include<stdio.h>
#include<stdlib.h>
intmain()
{
unsignedintm,p,n;
inti,j,k,*a,*b,*c;
scanf("%d%d%d",&m,&p,&n);
a=malloc(m*p*sizeof(int));
for(i=0;i<m;i++)
for(j=0;j<p;j++)
scanf("%d",&a[i*p+j]);
b=malloc(p*n*sizeof(int));
for(i=0;i<p;i++)
for(j=0;j<n;j++)
scanf("%d",&b[i*n+j]);
c=calloc(m*n,sizeof(int));
for(i=0;i<m;i++)
for(j=0;j<n;j++)
for(k=0;k<p;k++)
c[i*n+j]+=a[i*p+k]*b[k*n+j];
printf(" ");
for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
printf("%d",c[i*n+j]);
printf(" ");
}
free(a);free(b);free(c);
return0;
}
F. c语言,两矩阵相乘问题
没有细看,是不是维数问题。
给你一个程序,计算:
a[a_row][a_col] * b[b_row][b_col]=c[a_row][b_col]
注意 b_row==a_col
DEBUG 为 1 输出中间结果,方便查错。 可以改为 0。
动态分配数组。
含输入例子。
#include<stdio.h>
#include <stdlib.h>
#define DEBUG 1
void show_array(int *a, int row,int col);
/*---------------------
c[j]][i] = a[j][k] * b[k][i] = c[j][i]
a[c_row][nk]: a[3][2]={1,2,3,4,5,6}
b[nk][c_col]: b[2][3]={1,2,3,4,5,6}
c[c_row][c_col]: c[3][3]={9,12,15,19,26,33,29,40,51}
*---------------------*/
void matrix2(int *a,int *b, int *c, int a_row, int b_col, int b_row)
{
int i,j,k;
int a_col,c_row,c_col;
a_col=b_row; c_row=a_row; c_col=b_col;
if (DEBUG==1){
printf("A[%d][%d]:\n",a_row,a_col);
show_array(a,a_row,a_col);
printf("B[%d][%d]:\n",b_row,c_col);
show_array(b,b_row,b_col);
};
for(j=0;j<c_row;j++){
for(i=0;i<c_col;i++){
c[j*c_col+i]=0;
for(k=0;k<a_col;k++) c[j*c_col+i]+= a[j*a_col+k] * b[k*b_col+i];
};
};
}
main()
{
int i,j,k,tmp;
int a_row,a_col,b_row,b_col,c_row,c_col;
int *a,*b,*c;
printf("please enter a_row a_col of matrix_A: 3 2\n");
scanf("%d %d",&a_row,&a_col);
b_row = a_col;
printf("please enter b_col of matrix_B: 3\n");
scanf("%d",&b_col);
c_row = a_row;
c_col = b_col;
a = (int *) malloc(sizeof(int *) * a_row * a_col);
b = (int *) malloc(sizeof(int *) * b_row * b_col);
c = (int *) malloc(sizeof(int *) * c_row * c_col);
if (!c) { printf("no enought memory for alloc\n");return 0; }
printf("Please input int matrix a[%d][%d] 1 2 3 4 5 6\n",a_row,a_col);
for (j=0;j<a_row;j++)
for (i=0;i<a_col;i++){
scanf("%d",&tmp);
a[j*a_col+i] = tmp;
}
printf("Please input int matrix b[%d][%d] 1 2 3 4 5 6\n",b_row,b_col);
for (j=0;j<b_row;j++)
for (i=0;i<b_col;i++){
scanf("%d",&tmp);
b[j*b_col+i] = tmp;
}
if (DEBUG==1){
printf("A[%d][%d]:\n",a_row,a_col);
show_array(a,a_row,a_col);
printf("B[%d][%d]:\n",b_row,c_col);
show_array(b,b_row,b_col);
};
matrix2( a,b,c, a_row, b_col, b_row);
show_array(c,c_row,c_col);
free(a);free(b);free(c);
return 0;
}
void show_array(int *a, int row,int col)
{
int i,j;
for (j=0;j<row;j++){
for (i=0;i<col;i++) printf("%d ",a[j*col+i]);
printf("\n");
}
printf("--------------------\n");
}
G. 求高手帮忙用c语言求出两个矩阵的乘积
1、首先打开编程软件,新建一个项目,添加一个double.cpp文件,如图所示。
H. C语言求两矩阵乘积
#include"stdio.h"
#defineLANGA4
#defineLANGB3
#defineLANGCLANGB
intmulti(doubleA[][LANGA],doubleB[][LANGB],doubleC[][LANGC],intlanga,intlangb,intlangc){
inti,j,k;
if(langa==LANGB&&langb==LANGA){
for(i=0;i<langa;i++)
for(j=0;j<langc;j++)
for(C[i][j]=k=0;k<langb;k++)
C[i][j]+=A[i][k]*B[k][j];
return0;
}
return-1;
}
intmain(intargv,char*argc[]){
doubleA[LANGB][LANGA],B[LANGA][LANGB],C[LANGB][LANGC];
inti,j;
for(i=0;i<LANGB;i++)
for(j=0;j<LANGA;A[i][j++]=1.0/(i+j+1));
for(i=0;i<LANGA;i++)
for(j=0;j<LANGB;B[i][j]=1.0+j++);
if(multi(A,B,C,LANGB,LANGA,LANGB))
return-1;
for(i=0;i<LANGB;i++){
for(j=0;j<LANGC;printf("%12.4le",C[i][j++]));
printf(" ");
}
return0;
}
运行样例:
I. c语言初学求两个矩阵乘积问题
若然一矩阵的列数与另一矩阵的行数相等,则可定义这两个矩阵的乘积。如A是m×n矩阵和B是n×p矩阵,它们是乘积AB是一个m×p矩阵,其中
(AB)[i,j]=A[i,1]*B[1,j]+A[i,2]*B[2,j]+...+A[i,n]*B[n,j]对所有i及j。
此乘法有如下性质:
(AB)C=A(BC)对所有k×m矩阵A,m×n矩阵B及n×p矩阵C("结合律").
(A+B)C=AC+BC对所有m×n矩阵A及B和n&tiems;k矩阵C("分配律")。
C(A+B)=CA+CB对所有m×n矩阵A及B和k×m矩阵C("分配律")。
要注意的是:可置换性不一定成立,即有矩阵A及B使得AB≠BA。
Q=M*N(Q的行数和M相同,列数和N相同)
其中,M是m1*n1矩阵,N是m2*n2矩阵。当n1=m2时有:
for(i=1;i<=m1;++i)
for(j=1;j<=n2;++j){
Q[i][j]=0;
for(k=1;k<=n1;++k)Q[i][j]+=M[i][k]*N[k][j];
}
J. C语言:编写完整程序:计算两个矩阵的乘积。
#include<iostream.h>
#include<stdlib.h>
#include<iomanip.h>
void caidan();
#define max 12500
typedef struct{
int e;
int i; //行下标
int j; //列下标
}triple;
typedef struct{
triple data[max+1];
int rpos[10];
int m,n,t;
}ts;
//矩阵的转置实现函数
void fasttransposesmatrix(ts M,ts &T)
{
int col,num[10],cpot[10]={0},p,q;
T.m=M.n;
T.n=M.m;
T.t=M.t;
if(T.t){
for(col=1;col<=M.n;++col)
num[col]=0;
for(int t1=1;t1<=M.t;++t1)
++num[M.data[t1].j];
cpot[1]=1;
for(col=2;col<=M.n;++col)
cpot[col]=cpot[col-1]+num[col-1];
for(p=1;p<=M.t;++p){
col=M.data[p].j;
q=cpot[col];
T.data[q].i=M.data[p].j;
T.data[q].j=M.data[p].i;
T.data[q].e=M.data[p].e;
++cpot[col];
}
}
}
//矩阵的乘法Y=M*Q实现函数
void multsmatrix(ts M,ts Q,ts &Y){
int arow,brow=0,p,q,i1,tp,t1,ctemp[10],ccol;
if(M.n!=Q.m)
cout<<"错误!"<<endl;
Y.m=M.m;
Y.n=Q.n;
Y.t=0;
if(M.t*Q.t!=0)
{
for(arow=1;arow<=M.m;++arow)
{
for(i1=1;i1<=M.m;i1++)
ctemp[i1]=0;
Y.rpos[arow]=Y.t+1;
if(arow<M.m)
tp=M.rpos[arow+1];
else
tp=M.t+1;
for(p=M.rpos[arow];p<tp;++p)
{
brow=M.data[p].j;
if(brow<Q.m)
t1=Q.rpos[brow+1];
else
t1=Q.t+1;
for(q=Q.rpos[brow];q<t1;++q)
{
ccol=Q.data[q].j;
ctemp[ccol]+=M.data[p].e*Q.data[q].e;
}
}
for(ccol=1;ccol<=Y.n;++ccol)
if(ctemp[ccol]){
++Y.t;
if(++Y.t>max)
cout<<"错误"<<endl;
Y.data[Y.t].i=arow;
Y.data[Y.t].j=ccol;
Y.data[Y.t].e=ctemp[ccol];
cout<<Y.data[Y.t].i<<Y.data[Y.t].j<<Y.data[Y.t].e<<endl;
}
}
}
}
//将数组转换成三元组
void exchange1(int **a,ts &M,int m,int n)
{
int i1,j1,z=1,pos[10];
for(i1=0;i1<=n;i1++)
{
pos[i1]=0;
M.rpos[i1]=0;
}
pos[0]=1;
for(i1=0;i1<m;i1++)
for(j1=0;j1<n;j1++)
if(a[i1][j1])
{
M.data[z].i=i1+1;
M.data[z].j=j1+1;
M.data[z].e=a[i1][j1];
z++;
pos[i1+1]++;
}
for(i1=0;i1<=m;i1++)
{
M.rpos[i1+1]=M.rpos[i1]+pos[i1];
}
M.t=z-1;
M.m=m;
M.n=n;
}
//将三元组转换成数组形式
void exchange2(int **a,ts &T)
{
int i1,j1,z1;
for(i1=0;i1<T.m;i1++)
for(j1=0;j1<T.n;j1++)
{
a[i1][j1]=0;
}
for(z1=1;z1<=T.t;z1++)
for(i1=0;i1<T.m;i1++)
for(j1=0;j1<T.n;j1++)
if((i1==T.data[z1].i-1)&&(j1==T.data[z1].j-1))
{
a[i1][j1]=T.data[z1].e;
}
}
//输出结果
void output(int **a,int m,int n)
{
int i1,j1;
for(i1=0;i1<m;i1++)
{
for(j1=0;j1<n;j1++)
{
cout<<a[i1][j1]<<',';
}cout<<endl;
}
}
//转置的操作函数
void zhuan(){
int i1,j1,m,n;
cout<<"初始矩阵:"<<endl;
cout<<"请输入你要输入的数组的行数:"<<endl;
cin>>m;
cout<<"请输入你要输入的数组的列数:"<<endl;
cin>>n;
int**pp=new int*[m];
for(i1=0;i1<m;i1++)
pp[i1]=new int[n];
cout<<"请依次输入矩阵的元素:"<<endl;
for(i1=0;i1<m;i1++)
for(j1=0;j1<n;j1++)
{
cin>>pp[i1][j1];
}
cout<<endl;
ts M,T;
exchange1(pp,M,m,n);
fasttransposesmatrix(M,T);
int**pp1=new int*[T.m];
for(i1=0;i1<T.m;i1++)
pp1[i1]=new int[T.n];
exchange2(pp1,T);
cout<<"转置后的矩阵是:"<<endl;
output(pp1,T.m,T.n);
delete pp;
delete pp1;
caidan();
}
//相乘的操作函数
void xiangcheng(){
ts M,Q,Y;
int i1,j1,m,n;
cout<<"请输入你要输入的数组的行数:"<<endl;
cin>>m;
cout<<"请输入你要输入的数组的列数:"<<endl;
cin>>n;
int**p=new int*[m];
for(i1=0;i1<m;i1++)
p[i1]=new int[n];
cout<<"请依次输入矩阵的元素:"<<endl;
for(i1=0;i1<m;i1++)
for(j1=0;j1<n;j1++)
{
cin>>p[i1][j1];
}
cout<<endl;
exchange1(p,M,m,n);
cout<<"请再次输入一个矩阵(即乘数):"<<endl;
cout<<"请输入你要输入的数组的行数:"<<endl;
cin>>m;
cout<<"请输入你要输入的数组的列数:"<<endl;
cin>>n;
int**p2=new int*[m];
for(i1=0;i1<m;i1++)
p2[i1]=new int[n];
cout<<"请依次输入矩阵的元素:"<<endl;
for(i1=0;i1<m;i1++)
for(j1=0;j1<n;j1++)
{
cin>>p2[i1][j1];
}
cout<<endl;
exchange1(p2,Q,m,n);
multsmatrix(M,Q,Y);
int**p1=new int*[Y.m];
for(i1=0;i1<Y.m;i1++)
p1[i1]=new int[Y.n];
exchange2(p1,Y);
cout<<"两矩阵相乘结果为:"<<endl;
output(p1,Y.m,Y.n);
delete p2;
delete p1;
caidan();
}
//菜单
void caidan()
{
int i;
cout<<endl<<endl<<endl<<endl<<endl<<endl<<endl<<endl;
cout<<" 1.矩阵的转置 2.两个矩阵的乘法 3.退出 "<<endl;
cout<<"请选择:"<<endl;
cin>>i;
switch(i){
case 1:zhuan();
break;
case 2:xiangcheng();
break;
case 3: exit(0);
}
}
//主函数
void main()
{
caidan();
}