C++矩阵运算的实现简单
短信预约 -IT技能 免费直播动态提醒
利用C++实现矩阵的构造，通过运算符的重载实现矩阵的乘法、加法等。并且实现矩阵形状的打印，矩阵的打印。
#include<iostream>
#include<memory>
#include<assert.h>
#include<stdlib.h>
using namespace std;
class Matrix{
public:
Matrix(int row, int col); //构造函数
Matrix(int row, int col, int num);//构造函数重载
~Matrix(); //析构函数
Matrix(const Matrix & other); //赋值函数
Matrix operator*(const Matrix& other); //矩阵相乘
Matrix operator+(const Matrix& other); //矩阵相加
Matrix operator-(const Matrix& other); //矩阵相减
int **a = nullptr; //初始化一共空指针
int row, col;
void shape(); //打印矩阵形状
void Ma_pri(); //打印矩阵
};
int main(){
Matrix a(2,1); //构造一个（2,1）矩阵
Matrix b(1,2); //构造一个（1，2）矩阵
a.a[0][0] = 4; //初始化矩阵
a.a[1][0] = 2;
b.a[0][0] = 3;
b.a[0][1] = 5;
a.shape(); //矩阵形状打印
b.shape();
Matrix c = a*b; //矩阵相乘
c.shape();
c.Ma_pri(); //矩阵打印
Matrix d(3,3,1);
d.Ma_pri();
system("pause");
return 0;
}
Matrix::Matrix(int row, int col){
this->row = row;
this->col = col;
this->a = new int*[row];
for(int i=0;i<this->row;i++){
a[i] = new int[this->col];
}
}
Matrix::Matrix(int row, int col, int num){
this->row = row;
this->col = col;
this->a = new int*[row];
for(int i=0;i<this->row;i++){
a[i] = new int[this->col];
}
for(int i = 0; i < this->row; i++){
for(int j =0; j <this->row; j++){
this->a[i][j] = num;
}
}
}
Matrix::~Matrix(){
for(int i=0;i<this->row;i++){
if(a[i] != nullptr){
delete[] a[i];
a[i] = nullptr;
}
}
if(a != nullptr){
delete[] a;
a = nullptr;
}
}
Matrix::Matrix(const Matrix& other){
row = other.row;
col = other.col;
a = new int*[row];
for(int i=0;i<row;i++){
a[i] = new int[col];
memcpy(a[i], other.a[i],sizeof(int)*col);
}
}
Matrix Matrix::operator*(const Matrix& other){
if(this->col != other.row){
cout<<"shape error"<<endl;
exit(0);
}
Matrix m(this->row,other.col);
for(int i=0; i<this->row; i++){
for(int j=0;j<other.col;j++){
int sum = 0;
for(int k=0;k<this->col;k++){
sum += this->a[i][k] * other.a[k][j];
}
m.a[i][j] = sum;
}
}
return m;
}
Matrix Matrix::operator+(const Matrix& other){
if(this->col != other.col or this->row != other.row){
cout<<"shape error"<<endl;
exit(0);
}
Matrix m(this->row,this->col);
for(int i = 0;i < this->row; i++){
for(int j = 0; j < this-> col; j++){
m.a[i][j] = this->a[i][j] + other.a[i][j];
}
}
return m;
}
Matrix Matrix::operator-(const Matrix& other){
if(this->col != other.col or this->row != other.row){
cout<<"shape error"<<endl;
exit(0);
}
Matrix m(this->row,this->col);
for(int i = 0;i < this->row; i++){
for(int j = 0; j < this-> col; j++){
m.a[i][j] = this->a[i][j] - other.a[i][j];
}
}
return m;
}
void Matrix::shape(){
cout<<"("<<this->row<<","<<this->col<<")"<<endl;
}
void Matrix::Ma_pri(){
for(int i = 0; i < this->row; i++){
for(int j =0; j <this->row; j++){
cout<<this->a[i][j]<<" ";
}
cout<<endl;
}
}
矩阵求逆算法及程序实现
在做课题时，遇到了求多项式问题，利用了求逆方法。矩阵求逆一般使用简单的算法，还有快速算法 如全选主元高斯-约旦消元法，但本文程序主要写了简单的矩阵求逆算法定义法之伴随矩阵求逆公式如下，其中A可逆：A^{-1}=\frac{A^*}{|A|}，其中A^*是A的伴随矩阵。。
1.给定一个方阵，非奇异（不是也可，程序有考虑）；
2.由矩阵得到其行列式，求其值如|A|；
3.求其伴随矩阵A^*；
4.得到其逆矩阵。
主要函数如下：
//得到给定矩阵class="lazy" data-src的逆矩阵保存到des中。
bool GetMatrixInverse(double class="lazy" data-src[N][N],int n,double des[N][N])
{
double flag=getA(class="lazy" data-src,n);
double t[N][N];
if(flag==0)
{
return false;
}
else
{
getAStart(class="lazy" data-src,n,t);
for(int i=0;i<n;i++)
{
for(int j=0;j<n;j++)
{
des[i][j]=t[i][j]/flag;
}
}
}
return true;
}计算|A|:
//按第一行展开计算|A|
double getA(double arcs[N][N],int n)
{
if(n==1)
{
return arcs[0][0];
}
double ans = 0;
double temp[N][N]={0.0};
int i,j,k;
for(i=0;i<n;i++)
{
for(j=0;j<n-1;j++)
{
for(k=0;k<n-1;k++)
{
temp[j][k] = arcs[j+1][(k>=i)?k+1:k];
}
}
double t = getA(temp,n-1);
if(i%2==0)
{
ans += arcs[0][i]*t;
}
else
{
ans -= arcs[0][i]*t;
}
}
return ans;
}计算伴随矩阵：
//计算每一行每一列的每个元素所对应的余子式，组成A*
void getAStart(double arcs[N][N],int n,double ans[N][N])
{
if(n==1)
{
ans[0][0] = 1;
return;
}
int i,j,k,t;
double temp[N][N];
for(i=0;i<n;i++)
{
for(j=0;j<n;j++)
{
for(k=0;k<n-1;k++)
{
for(t=0;t<n-1;t++)
{
temp[k][t] = arcs[k>=i?k+1:k][t>=j?t+1:t];
}
}
ans[j][i] = getA(temp,n-1);
if((i+j)%2 == 1)
{
ans[j][i] = - ans[j][i];
}
}
}
}到此这篇关于C++矩阵运算的实现简单的文章就介绍到这了,更多相关C++ 矩阵运算内容请搜索编程网以前的文章或继续浏览下面的相关文章希望大家以后多多支持编程网！
免责声明：
① 本站未注明“稿件来源”的信息均来自网络整理。其文字、图片和音视频稿件的所属权归原作者所有。本站收集整理出于非商业性的教育和科研之目的，并不意味着本站赞同其观点或证实其内容的真实性。仅作为临时的测试数据，供内部测试之用。本站并未授权任何人以任何方式主动获取本站任何信息。
② 本站未注明“稿件来源”的临时测试数据将在测试完成后最终做删除处理。有问题或投稿请发送至: 邮箱/279061341@qq.com QQ/279061341
