矩阵的加法和乘法

编写一个矩阵的加法函数、一个矩阵的乘法函数，然后调用函数计算矩阵的和与积，并输出。
为了使函数更加通用，本次课程设计使用一维数组来存放矩阵。问题的关键是如何将一维数组与二维矩阵的元素之间建立对应关系。
问题拓展：
1、增加求矩阵的转置、求矩阵的秩、求正方阵的行列式，求矩阵的鞍点元素的位置与值。
2、将矩阵进行初等变换，行初等变换、列初等变换，变换成上三角矩阵或下三角矩阵。

#include<iostream>
using namespace std;
#include<cstdlib>
#include<ctime>
#include<iomanip>
void MatAdd(double a[],double b[],double c[],int n) //矩阵的加法
{
for(int i=0;i<n;i++)
c[i]=a[i]+b[i];
}
void MatMull(double a[],double b[],double c[],int m,int n,int p)//矩阵的乘法
{
int i,j,k;
for(i=0;i<m;i++)
for(j=0;j<p;j++)
{
c[i*p+j]=0;
for(k=0;k<n;k++)
c[i*p+j]+=a[i*n+k]*b[k*p+j];
}
}
void ShowMat(double a[],int m,int n)
//输出矩阵
{
for(int i=0;i<m;i++)
{
for(int j=0;j<n;j++)
cout<<setw(8)<<fixed<<setprecision(2)<<a[i*n+j];
cout<<endl;
}
cout<<endl;
}
void MatRand(double a[],int n)
//使用随机数产生一个矩阵
{
for(int i=0;i<n;i++)
a[i]=(rand()%100)/10.0;
}
void fun1()
//测试矩阵的加法
{
double a[4][5],b[4][5],c[4][5];
MatRand(a[0],20);
MatRand(b[0],20);
MatAdd(a[0],b[0],c[0],20);
ShowMat(a[0],4,5);
ShowMat(b[0],4,5);
ShowMat(c[0],4,5);
}
void fun2() //测试矩阵的乘法
{
double x[4][5],y[5][8],z[4][8];
MatRand(x[0],4*5);
MatRand(y[0],5*8);
MatMull(x[0],y[0],z[0],4,5,8);
ShowMat(x[0],4,5);
ShowMat(y[0],5,8);
ShowMat(z[0],4,8);
}
void main()
{
srand(time(0));
fun1();
fun2();
}

最后

以上就是忧虑宝马为你收集整理的矩阵的加法和乘法的全部内容，希望文章能够帮你解决矩阵的加法和乘法所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。