线性代数矩阵乘法用C++代码实现,再也不用担心算错啦!

c++实现矩阵相乘

    • 矩阵乘法的定义
    • 矩阵相乘的代码
    • 完整的程序
    • 实现线代其它操作的参考链接

线性代数矩阵相乘一旦阶数多了算的可真是让人CPU疼,但可爱的计算机是不怕累又超听话,所以用一个c++程序帮助你轻松解决矩阵相乘!!

矩阵乘法的定义

设A为m * n的矩阵,B为n * p的矩阵,那么称m * p的矩阵C为矩阵A与B的乘积,记作C=AB ,其中矩阵C中的第i行第j列元素可以表示为:
在这里插入图片描述
线性代数矩阵乘法用C++代码实现,再也不用担心算错啦!_第1张图片

矩阵相乘的代码

这个算法很简单,直接利用矩阵相乘的定义及运算法则即可:
话不多说了,直接上代码:

vector<vector<double> > multiply(vector<vector<double> > matrix1,vector<vector<double> > matrix2)
{
        
    if(matrix1[0].size()!=matrix2.size())
    {
     
    	cout<<" 这两个矩阵不能相乘";
    	exit(1);
		 
	}
	int n = matrix1.size(),m=matrix2[0].size();
	int l =matrix2.size();
	vector<vector<double> > matrix3(n);
	for(int i=0;i<n;i++)
	for(int j=0;j<m;j++)
	{
        double t =0;
	   //执行一次相乘
		for(int k=0;k<l;k++)
		{
     
			t+=matrix1[i][k]*matrix2[k][j];
		}
		matrix3[i].push_back(t);
	}
	return matrix3;
}

完整的程序

给出完整的程序去供直接使用

#include
using namespace std;
vector<vector<double> > multiply(vector<vector<double> > matrix1,vector<vector<double> > matrix2)
{
        
    if(matrix1[0].size()!=matrix2.size())
    {
     
    	cout<<" 这两个矩阵不能相乘";
    	exit(1);
		 
	}
	int n = matrix1.size(),m=matrix2[0].size();
	int l =matrix2.size();
	vector<vector<double> > matrix3(n);
	for(int i=0;i<n;i++)
	for(int j=0;j<m;j++)
	{
        double t =0;
		for(int k=0;k<l;k++)
		{
     
			t+=matrix1[i][k]*matrix2[k][j];
		}
		matrix3[i].push_back(t);
	}
	return matrix3;
}
void output(vector<vector<double> > det)
{
     
	int m=det.size(),n=det[0].size();
	for(int i=0;i<m;i++)
	{
     
		for(int j=0;j<n;j++)
	      cout<<det[i][j]<<"  ";
	      cout<<endl;
	}
}
int main()
{
        int n1,m1;
	cout<<"Enter matrix1(n1*m1) n1,m1:";
	cin>>n1;cin>>m1;
	vector<vector<double> > matrix1(n1);//初始化det1为n1行 
	
	for(int i=0;i<n1;i++)
	for(int j=0;j<m1;j++)
	{
     
		double t;
		cin>>t;
		matrix1[i].push_back(t); 
	}
	
	int n2,m2;
	cout<<"Enter matrix2(n2*m2) n2,m2:";
	cin>>n2;cin>>m2;
	vector<vector<double> >matrix2(n2);//初始化det1为n1行 
	
	for(int i=0;i<n2;i++)
	for(int j=0;j<m2;j++)
	{
     
		double t;
		cin>>t;
		matrix2[i].push_back(t); 
	}
	output(multiply(matrix1,matrix2));
 } 

实现线代其它操作的参考链接

  • 线性代数行列式求值/矩阵相乘/求矩阵的逆,一个c++程序全部解决
  • 线性代数求行列式的值(用C++程序轻松解决)
  • 让c++程序助你轻松求矩阵的逆

你可能感兴趣的:(线性代数,线性代数,c++,算法)