现在来做一些现代里面的矩阵运算:
矩阵的加减运算,矩阵的数乘运算
/**
MatrixAddition.c
实现两个矩阵相加
*/
#include<stdio.h>
const int COLS=3;
void printMatirx(int *pArray,int rows,int cols);
void printMatirx2(int (*pArray)[COLS],int rows);
int main()
{
int A[3][4]={{15,10,9,12},
{18,14,8,7},
{16,13,6,11}};
printf("矩阵A=\n");
printMatirx(A,3,4);
int B[3][4]={{4,3,5,2},
{0,9,6,1},
{5,7,2,6}};
printf("矩阵B=\n");
printMatirx(B,3,4);
int C[3][4];
int i,j;
//矩阵相加:两个矩阵必须行数和列数一样才能相加,
//和矩阵的每个元素分别是两个矩阵对应的元素的和
printf("矩阵A+矩阵B=\n");
for(i=0;i<3;i++)//控制行
{
for(j=0;j<4;j++)
{
C[i][j]=A[i][j]+B[i][j];
// printf("%3d",C[i][j]);//输出结果
}
// printf("\n");
}
int (*p)[COLS]=C;
printMatirx2(p,3);
/*结果应为:
19 13 14 14
18 23 14 8
21 20 8 17
*/
printf("矩阵C-矩阵A=\n");
for(i=0;i<3;i++)//控制行
{
for(j=0;j<4;j++)
{
printf("%3d",C[i][j]-A[i][j]);//输出结果
}
printf("\n");
}
/*结果应为:
4 3 5 2
0 9 6 1
5 7 2 6
*/
//矩阵数乘:
int D[3][4];
printf("矩阵D:\n");
for(i=0;i<3;i++)//控制行
{
for(j=0;j<4;j++)
{
D[i][j]=1;
printf("%d\t", D[i][j]);//输出结果
}
printf("\n");
}
int mul;
printf("矩阵D数乘以:");
scanf("%d",&mul);
printf("矩阵D数乘以%d=\n",mul);
for(i=0;i<3;i++)//控制行
{
for(j=0;j<4;j++)
{
D[i][j]=D[i][j]*mul;
printf("%d\t", D[i][j]);//输出结果
}
printf("\n");
}
return 0;
}
//使用一级指针来访问二维数组
void printMatirx(int *pArray,int rows,int cols)
{
int i;
int j;
for(i=0;i<rows;i++)
{
for(j=0;j< cols;j++)
{
printf("%d\t",*(pArray+i*cols+j));
}
printf("\n");
}
}
//使用指向一维数组的指针,访问二维数组
void printMatirx2(int (*pArray)[COLS],int rows)
{
int i,j;
for(i=0;i<rows;i++)
{
for(j=0;j<COLS;j++)
{
printf("%d\t",pArray[i][j]);
/*
[]运算符的结合方向是由左向右,pArray[1][2]就等价于(* (pArray + 1))[2],
而由于pArray是数组指针,而且数组的长度为3,
所以* (pArray + 1):表示跨过3个元素,而也就是到array第二行
即表示array[1]这个数组,则pArray[1][2]则就完全等价于array[1][2]。
*/
}
printf("\n");
}
}
结果:
矩阵A=
15 10 9 12
18 14 8 7
16 13 6 11
矩阵B=
4 3 5 2
0 9 6 1
5 7 2 6
矩阵A+矩阵B=
19 13 14
14 18 23
14 8 21
矩阵C-矩阵A=
4 3 5 2
0 9 6 1
5 7 2 6
矩阵D:
1 1 1 1
1 1 1 1
1 1 1 1
矩阵D数乘以:666
矩阵D数乘以666=
666 666 666 666
666 666 666 666
666 666 666 666
Process returned 0 (0x0) execution time : 9.059 s
Press any key to continue.