中文字幕av专区_日韩电影在线播放_精品国产精品久久一区免费式_av在线免费观看网站

溫馨提示×

C++矩陣類中如何實現矩陣的行列式計算

c++
小樊
103
2024-08-29 22:38:39
欄目: 編程語言

在C++中,可以通過創建一個名為Matrix的類來實現矩陣的行列式計算。以下是一個簡單的示例:

#include<iostream>
#include<vector>

class Matrix {
public:
    // 構造函數
    Matrix(int rows, int cols) : rows_(rows), cols_(cols) {
        data_.resize(rows);
        for (int i = 0; i< rows; ++i) {
            data_[i].resize(cols);
        }
    }

    // 設置矩陣元素
    void setElement(int row, int col, double value) {
        if (row >= 0 && row< rows_ && col >= 0 && col< cols_) {
            data_[row][col] = value;
        } else {
            std::cerr << "Error: Invalid index."<< std::endl;
        }
    }

    // 獲取矩陣元素
    double getElement(int row, int col) const {
        if (row >= 0 && row< rows_ && col >= 0 && col< cols_) {
            return data_[row][col];
        } else {
            std::cerr << "Error: Invalid index."<< std::endl;
            return 0;
        }
    }

    // 計算行列式
    double determinant() const {
        if (rows_ != cols_) {
            std::cerr << "Error: The matrix is not square."<< std::endl;
            return 0;
        }

        if (rows_ == 1) {
            return data_[0][0];
        }

        if (rows_ == 2) {
            return data_[0][0] * data_[1][1] - data_[0][1] * data_[1][0];
        }

        double det = 0;
        for (int i = 0; i< cols_; ++i) {
            Matrix submatrix(rows_ - 1, cols_ - 1);
            for (int j = 1; j< rows_; ++j) {
                for (int k = 0; k< cols_; ++k) {
                    if (k < i) {
                        submatrix.setElement(j - 1, k, data_[j][k]);
                    } else if (k > i) {
                        submatrix.setElement(j - 1, k - 1, data_[j][k]);
                    }
                }
            }
            det += ((i % 2 == 0) ? 1 : -1) * data_[0][i] * submatrix.determinant();
        }
        return det;
    }

private:
    int rows_;
    int cols_;
    std::vector<std::vector<double>> data_;
};

int main() {
    Matrix matrix(3, 3);
    matrix.setElement(0, 0, 1);
    matrix.setElement(0, 1, 2);
    matrix.setElement(0, 2, 3);
    matrix.setElement(1, 0, 4);
    matrix.setElement(1, 1, 5);
    matrix.setElement(1, 2, 6);
    matrix.setElement(2, 0, 7);
    matrix.setElement(2, 1, 8);
    matrix.setElement(2, 2, 9);

    std::cout << "The determinant of the matrix is: "<< matrix.determinant()<< std::endl;

    return 0;
}

這個示例中的Matrix類包含了一個二維向量data_來存儲矩陣的元素。構造函數用于初始化矩陣的行數和列數。setElement和getElement方法用于設置和獲取矩陣元素。determinant方法用于計算矩陣的行列式。

注意:這個示例僅適用于方陣(行數和列數相等的矩陣)。對于非方陣,需要使用其他方法來計算行列式。

0
大荔县| 金昌市| 天津市| 山东省| 天气| 会泽县| 曲靖市| 观塘区| 班玛县| 罗山县| 新田县| 平乐县| 虞城县| 金秀| 涟水县| 米脂县| 德惠市| 东兰县| 类乌齐县| 柳江县| 商丘市| 喜德县| 图木舒克市| 沾益县| 贵州省| 吉首市| 东莞市| 南充市| 图片| 承德市| 神池县| 安顺市| 天长市| 霍山县| 绩溪县| 甘肃省| 桑植县| 苗栗市| 永春县| 车致| 枣强县|