C++矩陣運算的實現(xiàn)簡單

更新時間：2021年09月03日 09:55:30 作者：victor_JZ

本文主要介紹了C++矩陣運算的實現(xiàn)簡單，文中通過示例代碼介紹的非常詳細，具有一定的參考價值，感興趣的小伙伴們可以參考一下

利用C++實現(xiàn)矩陣的構造，通過運算符的重載實現(xiàn)矩陣的乘法、加法等。并且實現(xiàn)矩陣形狀的打印，矩陣的打印。

#include<iostream>
#include<memory>
#include<assert.h>
#include<stdlib.h>
using namespace std;

class Matrix{
    public:
        Matrix(int row, int col);   //構造函數(shù)
        Matrix(int row, int col, int num);//構造函數(shù)重載
        ~Matrix();                  //析構函數(shù)
        Matrix(const Matrix & other); //賦值函數(shù)
        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;
    }
}

矩陣求逆算法及程序實現(xiàn)

　在做課題時，遇到了求多項式問題，利用了求逆方法。矩陣求逆一般使用簡單的算法，還有快速算法如全選主元高斯-約旦消元法，但本文程序主要寫了簡單的矩陣求逆算法定義法之伴隨矩陣求逆公式如下，其中A可逆：A^{-1}=\frac{A^*}{|A|}，其中A^*是A的伴隨矩陣。。

　　1.給定一個方陣，非奇異（不是也可，程序有考慮）；

　　2.由矩陣得到其行列式，求其值如|A|；

　　3.求其伴隨矩陣A^*；

　　4.得到其逆矩陣。

主要函數(shù)如下：

//得到給定矩陣src的逆矩陣保存到des中。
bool GetMatrixInverse(double src[N][N],int n,double des[N][N])
{
    double flag=getA(src,n);
    double t[N][N];
    if(flag==0)
    {
        return false;
    }
    else
    {
        getAStart(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++矩陣運算的實現(xiàn)簡單的文章就介紹到這了,更多相關C++ 矩陣運算內容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關文章希望大家以后多多支持腳本之家！

您可能感興趣的文章:

Qt 添加MSVC2017編譯器的完整教程(保姆級)
本文主要介紹了Qt 添加MSVC2017編譯器的完整教程,文中通過圖文介紹的非常詳細,對大家的學習或者工作具有一定的參考學習價值,需要的朋友們下面隨著小編來一起學習學習吧
2024-06-06
QT5實現(xiàn)UDP通信的示例代碼
本文主要介紹了QT5實現(xiàn)UDP通信的示例代碼，主要使用QUdpSocket類用于實現(xiàn)UDP通信，文中通過示例代碼介紹的非常詳細，具有一定的參考價值，感興趣的小伙伴們可以參考一下
2021-12-12
c語言字符串轉大寫的簡單實例
這篇文章主要介紹了c語言字符串轉大寫的簡單實例，有需要的朋友可以參考一下
2013-12-12
Qt中QSettings配置文件的讀寫和應用場景詳解
這篇文章主要給大家介紹了關于Qt中QSettings配置文件的讀寫和應用場景的相關資料,QSettings能讀寫配置文件,當配置文件不存在時,可生成配置文件,文中通過代碼介紹的非常詳細,需要的朋友可以參考下
2023-10-10
利用Matlab實現(xiàn)陰影柱狀圖的繪制
陰影柱狀圖繪制的代碼MATHWORKS上也有，不過是生成圖片后，識別圖像上不同顏色塊并直接修改像素生成陰影線。這樣的生成方式不可逆且自由度較低，所以本文為大家?guī)淼氖且粋€直接畫線條填充的硬畫版本，希望有所幫助
2022-09-09
C語言入門學習之fgets()函數(shù)和fputs()函數(shù)
fgetc() 和 fputc() 函數(shù)每次只能讀寫一個字符,速度較慢,實際開發(fā)中往往是每次讀寫一個字符串或者一個數(shù)據塊,這樣能明顯提高效率,這篇文章主要給大家介紹了關于C語言入門學習之fgets()函數(shù)和fputs()函數(shù)的相關資料,需要的朋友可以參考下
2021-11-11
用C語言winform編寫滲透測試工具實現(xiàn)SQL注入功能
本篇文章主要介紹使用C#winform編寫滲透測試工具，實現(xiàn)SQL注入的功能。使用python編寫SQL注入腳本，基于get顯錯注入的方式進行數(shù)據庫的識別、獲取表名、獲取字段名，最終獲取用戶名和密碼；使用C#winform編寫windows客戶端軟件調用.py腳本，實現(xiàn)用戶名和密碼的獲取
2021-08-08
C++begin和end運算符的返回迭代器的類型如何判斷?
今天小編就為大家分享一篇關于C++begin和end運算符的返回迭代器的類型如何判斷?，小編覺得內容挺不錯的，現(xiàn)在分享給大家，具有很好的參考價值，需要的朋友一起跟隨小編來看看吧
2019-04-04
常用的C語言編程工具匯總
c語言編程軟件適于編寫系統(tǒng)軟件，是學習編程的同學們的必備軟件。c語言一種非常強大的計算機語言，應用非常廣泛，不僅僅是在軟件開發(fā)上，而且各類科研都會用到c語言。今天小編給大家匯總下C語言的編程工具
2018-01-01
C++中單鏈表操作的示例代碼
這篇文章主要為大家詳細介紹了C++中單鏈表操作的示例代碼，主要包括單鏈表的創(chuàng)建、元素的插入和刪除等。文中的代碼簡潔易懂，需要的可以參考一下
2022-10-10