如何使用C代码对一维矩阵求逆?

问题描述 投票:-1回答:1

你好,我用高斯约旦进行了1天,但我没有使用

我想找到一维矩阵逆。我找到行列式,但是我不知道这个矩阵的逆数

你好,亲爱的朋友

我们的矩阵:

double A[] = {6, 6 ,2, 4, 9 ,7, 4, 3 ,3};
double B[] = {6, 6 ,2, 4, 9 ,7, 4, 3 ,3};
double Final[9];
Function to calculate determinant:

int Inverse(double *A, double *C, int N){

    int n = N;
        int  i, j, k;
        float a[10][10] = { 0.0 };
        double C[9] = { 0.0 };
    float pivot = 0.0;
    float factor = 0.0;
    double sum = 0.0;  ```  variables



    for (k = 1; k <= n - 1; k++)
    {

        if (a[k][k] == 0.0)
        {

            printf("error");

        }
        else
        {
            pivot = a[k][k];
            for (j = k; j <= n + 1; j++)
                a[k][j] = a[k][j] / pivot;

            for (i = k + 1; i <= n; i++)
            {
                factor = a[i][k];

                for (j = k; j <= n + 1; j++)
                {
                    a[i][j] = a[i][j] - factor * a[k][j];
                }
            }

        }

        if (a[n][n] == 0)

            printf("error");

        else
        {

            C[n] = a[n][n + 1] / a[n][n];

            for (i = n - 1; i >= 1; i--)
            {

                sum = 0.0;

                for (j = i + 1; j <= n; j++)

                    sum = sum + a[i][j] * C[j];

                C[i] = (a[i][n + 1] - sum) / a[i][i];

            }

        }

    }

    for (i = 1; i <= n; i++)
    {

        printf("\n\tx[%1d]=%10.4f", i, C[i]);

    }

    system("PAUSE");
    return 0;

}

尽管我非常努力,但在C编程中找不到1x1维矩阵的反面。输出始终生成0。在我可能会出错的地方,您可以帮我吗?谢谢。

texture1d
1个回答
1
投票

似乎您正在使用C作为输出参数(以存储反函数);但是,您还要在函数中声明一个具有相同名称的局部变量。这将导致局部变量遮蔽(即:隐藏)输出参数;因此,您对函数中的C所做的更改不会影响调用函数看到的C

要解决此问题,您需要从函数中删除行double C[9] = {0};

© www.soinside.com 2019 - 2024. All rights reserved.