矩阵的常见操作
矩阵的转置
对于一个mxn的矩阵A[m][n],其转置矩阵是一个nxm的矩阵B[n][m],且B[i][j]=A[j][i],0<=i<n,0<=j<m。实现代码如下:
/* 矩阵的转置 */
/* A[][MAXSIZE]指的是原矩阵A;B[][MAXSIZE]指的是转置后的矩阵;m指的是是原矩阵A的行数;n指的是矩阵A的列数 */
void trsmat(int A[][MAXSIZE],int B[][MAXSIZE],int m,int n){
for(int i=0;i<m;i++){
for(int j=0;j<n;j++){// 循环遍历A矩阵的所有元素
B[j][i]=A[i][j];// 转置
}
}
}
完整代码:
#include<stdio.h>
#include<stdlib.h>
#define MAXSIZE 20
/* 矩阵的转置 */
/* A[][MAXSIZE]指的是原矩阵A;B[][MAXSIZE]指的是转置后的矩阵;m指的是是原矩阵A的行数;n指的是矩阵A的列数 */
void trsmat(int A[][MAXSIZE],int B[][MAXSIZE],int m,int n){
for(int i=0;i<m;i++){
for(int j=0;j<n;j++){// 循环遍历A矩阵的所有元素
B[j][i]=A[i][j];// 转置
}
}
}
/* 打印矩阵 */
/* A[][MAXSIZE]指的是矩阵;m指的是矩阵的行数;n指的是矩阵的列数 */
void printMat(int A[][MAXSIZE],int m,int n){
for(int i=0;i<m;i++){
for(int j=0;j<n;j++){
printf("%d\t",A[i][j]);
}
printf("\n");
}
}
int main(){
int A[][MAXSIZE]= {
{0,0,0,1},
{0,0,3,2},
{1,0,0,0},
{0,2,0,0},
{1,0,0,0}
};
int m=5,n=4;// 指的是矩阵的行数和列数
printf("原矩阵A:\n");
printMat(A,m,n);// 打印原矩阵
int B[MAXSIZE][MAXSIZE];
trsmat(A,B,m,n);// 矩阵的转置
printf("转置后的矩阵B:\n");
printMat(B,n,m);
return 0;
}
运行结果:
矩阵的相加
两个尺寸都是mxn的矩阵相加后得到依然是一个mxn的矩阵,相加规则为C[i][j]=A[i][j]+B[i][j],其中C[i][j]为结果矩阵C中的元素,A[i][j]和B[i][j]都是操作数矩阵A和B中的元素,可见矩阵相加操作为两矩阵对应位置上的元素逐一相加。
核心代码:
/* 矩阵的相加 */
/* C[][MAXSIZE]指的是相加的结果矩阵;A[][MAXSIZE]指要参与相加运算的矩阵;B[][MAXSIZE]指要参与相加运算的矩阵;m指的是矩阵的行数;n指的是矩阵的列数 */
void addmat(int C[][MAXSIZE],int A[][MAXSIZE],int B[][MAXSIZE],int m,int n){
for(int i=0;i<m;i++){
for(int j=0;j<n;j++){
C[i][j]=A[i][j]+B[i][j];// 矩阵元素的相加
}
}
}
完整代码如下:
#include<stdio.h>
#include<stdlib.h>
#define MAXSIZE 20
/* 矩阵的相加 */
/* C[][MAXSIZE]指的是相加的结果矩阵;A[][MAXSIZE]指要参与相加运算的矩阵;B[][MAXSIZE]指要参与相加运算的矩阵;m指的是矩阵的行数;n指的是矩阵的列数 */
void addmat(int C[][MAXSIZE],int A[][MAXSIZE],int B[][MAXSIZE],int m,int n){
for(int i=0;i<m;i++){
for(int j=0;j<n;j++){
C[i][j]=A[i][j]+B[i][j];// 矩阵元素的相加
}
}
}
/* 打印矩阵 */
/* A[][MAXSIZE]指的是矩阵;m指的是矩阵的行数;n指的是矩阵的列数 */
void printMat(int A[][MAXSIZE],int m,int n){
for(int i=0;i<m;i++){
for(int j=0;j<n;j++){
printf("%d\t",A[i][j]);
}
printf("\n");
}
}
int main(){
int A[][MAXSIZE]= {
{0,0,0,1},
{0,0,3,2},
{1,0,0,0},
{0,2,0,0},
{1,0,0,0}
};
int m=5,n=4;// 指的是矩阵的行数和列数
printf("原矩阵A:\n");
printMat(A,m,n);// 打印原矩阵
int B[][MAXSIZE]= {
{1,1,1,0},
{1,1,0,0},
{1,0,1,0},
{0,2,0,1},
{1,0,1,1}
};
printf("原矩阵B:\n");
printMat(B,m,n);// 打印原矩阵
int C[MAXSIZE][MAXSIZE];
addmat(C,A,B,m,n);// 矩阵的相加
printf("相加后的矩阵C:\n");
printMat(C,m,n);// 打印相加后的矩阵
return 0;
}
运行结果如下:
矩阵的相乘
假设两个矩阵A与B相乘,结果为C,C中第i行与第j列上的元素为A中第i行的元素与B中第j列元素对应相乘并且求和的结果。A和B两矩阵可以相乘的条件是A的列数必须等于B的行数。其中A的尺寸是mxn,B的尺寸为nxk。
核心代码:
/* 矩阵的相乘 */
/* C[][MAXSIZE]指的是相乘的结果矩阵;A[][MAXSIZE]指要参与相乘运算的矩阵;B[][MAXSIZE]指要参与相乘运算的矩阵;m指的是A矩阵的行数;n指的是A矩阵的列数也是B矩阵的行数;k指的是B矩阵的列数 */
void mutmat(int C[][MAXSIZE],int A[][MAXSIZE],int B[][MAXSIZE],int m,int n,int k){
for(int i=0;i<m;i++){
for(int j=0;j<k;j++){
C[i][j]=0;
for(int h=0;h<n;h++){
C[i][j]+=A[i][h]*B[h][j];// 相乘
}
}
}
}
完整代码如下:
#include<stdio.h>
#include<stdlib.h>
#define MAXSIZE 20
/* 矩阵的相乘 */
/* C[][MAXSIZE]指的是相乘的结果矩阵;A[][MAXSIZE]指要参与相乘运算的矩阵;B[][MAXSIZE]指要参与相乘运算的矩阵;m指的是A矩阵的行数;n指的是A矩阵的列数也是B矩阵的行数;k指的是B矩阵的列数 */
void mutmat(int C[][MAXSIZE],int A[][MAXSIZE],int B[][MAXSIZE],int m,int n,int k){
for(int i=0;i<m;i++){
for(int j=0;j<k;j++){
C[i][j]=0;
for(int h=0;h<n;h++){
C[i][j]+=A[i][h]*B[h][j];// 相乘
}
}
}
}
/* 打印矩阵 */
/* A[][MAXSIZE]指的是矩阵;m指的是矩阵的行数;n指的是矩阵的列数 */
void printMat(int A[][MAXSIZE],int m,int n){
for(int i=0;i<m;i++){
for(int j=0;j<n;j++){
printf("%d\t",A[i][j]);
}
printf("\n");
}
}
int main(){
int A[][MAXSIZE]= {
{0,0,0,1},
{0,0,3,2},
{1,0,0,0},
{0,2,0,0},
{1,0,0,0}
};
int m=5,n=4;// 指的是矩阵的行数和列数
printf("原矩阵A:\n");
printMat(A,m,n);// 打印原矩阵
int B[][MAXSIZE]= {
{1,1,1,0,1},
{1,1,0,0,0},
{1,0,1,0,1},
{1,0,1,1,1}
};
int k=5;
printf("原矩阵B:\n");
printMat(B,n,k);// 打印原矩阵
int C[MAXSIZE][MAXSIZE];
mutmat(C,A,B,m,n,k);// 矩阵的相乘
printf("相乘后的矩阵C:\n");
printMat(C,m,n);// 打印相乘后的矩阵
return 0;
}
运行结果如下:
矩阵的相减
原理和相加一样,只是换成了减号而已。
核心代码如下:
/* 矩阵的相减 */
/* C[][MAXSIZE]指的是相减的结果矩阵;A[][MAXSIZE]指要参与相减运算的矩阵;B[][MAXSIZE]指要参与相减运算的矩阵;m指的是矩阵的行数;n指的是矩阵的列数 */
void submat(int C[][MAXSIZE],int A[][MAXSIZE],int B[][MAXSIZE],int m,int n) {
for(int i=0; i<m; i++) {
for(int j=0; j<n; j++) {
C[i][j]=A[i][j]-B[i][j];// 矩阵元素的相减
}
}
}
完整代码如下:
#include<stdio.h>
#include<stdlib.h>
#define MAXSIZE 20
/* 矩阵的相减 */
/* C[][MAXSIZE]指的是相减的结果矩阵;A[][MAXSIZE]指要参与相减运算的矩阵;B[][MAXSIZE]指要参与相减运算的矩阵;m指的是矩阵的行数;n指的是矩阵的列数 */
void submat(int C[][MAXSIZE],int A[][MAXSIZE],int B[][MAXSIZE],int m,int n) {
for(int i=0; i<m; i++) {
for(int j=0; j<n; j++) {
C[i][j]=A[i][j]-B[i][j];// 矩阵元素的相减
}
}
}
/* 打印矩阵 */
/* A[][MAXSIZE]指的是矩阵;m指的是矩阵的行数;n指的是矩阵的列数 */
void printMat(int A[][MAXSIZE],int m,int n) {
for(int i=0; i<m; i++) {
for(int j=0; j<n; j++) {
printf("%d\t",A[i][j]);
}
printf("\n");
}
}
int main() {
int A[][MAXSIZE]= {
{0,0,0,1},
{0,0,3,2},
{1,0,0,0},
{0,2,0,0},
{1,0,0,0}
};
int m=5,n=4;// 指的是矩阵的行数和列数
printf("原矩阵A:\n");
printMat(A,m,n);// 打印原矩阵
int B[][MAXSIZE]= {
{0,0,0,1},
{0,0,3,2},
{1,0,0,0},
{0,2,0,0},
{1,0,0,0}
};
printf("原矩阵B:\n");
printMat(B,m,n);// 打印原矩阵
int C[MAXSIZE][MAXSIZE];
submat(C,A,B,m,n);// 矩阵的相减
printf("相减后的矩阵C:\n");
printMat(C,m,n);// 打印相减后的矩阵
return 0;
}
运行结果如下: