В функции обратной матрицы, не могу найти ошибку

@2garin_95 · Регистрация: 19.05.2013

У меня есть программа для нахождения обратной матрицы написанная на Dev C++.

C++

#include <iostream>
#include <conio.h>
 
void inversion(double **A, int N)
{
    double temp;
 
    double **E = new double *[N];
 
    for (int i = 0; i < N; i++)
        E[i] = new double [N];
 
    for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
        {
            E[i][j] = 0.0;
 
            if (i == j)
                E[i][j] = 1.0;
        }
 
    for (int k = 0; k < N; k++)
    {
        temp = A[k][k];
 
        for (int j = 0; j < N; j++)
        {
            A[k][j] /= temp;
            E[k][j] /= temp;
        }
 
        for (int i = k + 1; i < N; i++)
        {
            temp = A[i][k];
 
            for (int j = 0; j < N; j++)
            {
                A[i][j] -= A[k][j] * temp;
                E[i][j] -= E[k][j] * temp;
            }
        }
    }
 
    for (int k = N - 1; k > 0; k--)
    {
        for (int i = k - 1; i >= 0; i--)
        {
            temp = A[i][k];
 
            for (int j = 0; j < N; j++)
            {
                A[i][j] -= A[k][j] * temp;
                E[i][j] -= E[k][j] * temp;
            }
        }
    }
 
    for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
            A[i][j] = E[i][j];
 
    for (int i = 0; i < N; i++)
        delete [] E[i];
 
    delete [] E;
}
 
int main()
{
    int N;
 
    std::cout << "N: ";
    std::cin >> N;
 
    double **matrix = new double *[N];
 
    for (int i = 0; i < N; i++)
        matrix[i] = new double [N];
 
    for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
        {
            std::cout << "[" << i << "][" << j << "] = ";
            std::cin >> matrix[i][j];
        }
 
    inversion(matrix, N);
 
    for (int i = 0; i < N; i++)
    {
        for (int j = 0; j < N; j++)
            std::cout << matrix[i][j] << "  ";
 
        std::cout << std::endl;
    }
 
    for (int i = 0; i < N; i++)
        delete [] matrix[i];
 
    delete [] matrix;
 
    getch();
    return 0;
}

Я взял эту функцию и вставил ее в C++ Builder 6 с небольшим редактированием, поменял тип данных с double на float
Вот код:

C++

void inversion(float **A, int N)
{
float temp;
float **E = new float *[N];
for (int i = 0; i < N; i++)
E[i] = new float [N];
for (int i = 0; i < N; i++)
for (int j = 0; j < N; j++)
{
E[i][j] = 0.0;
if (i == j)
E[i][j] = 1.0;
}
for (int k = 0; k < N; k++)
{
temp = A[k][k];
for (int j = 0; j < N; j++)
{
A[k][j] /= temp;
E[k][j] /= temp;
}
for (int i = k + 1; i < N; i++)
{
temp = A[i][k];
for (int j = 0; j < N; j++)
{
A[i][j] -= A[k][j] * temp;
E[i][j] -= E[k][j] * temp;
}
}
}
for (int k = N - 1; k > 0; k--)
{
for (int i = k - 1; i >= 0; i--)
{
temp = A[i][k];
for (int j = 0; j < N; j++)
{
A[i][j] -= A[k][j] * temp;
E[i][j] -= E[k][j] * temp;
}
}
}
for (int i = 0; i < N; i++)
for (int j = 0; j < N; j++)
A[i][j] = E[i][j];
for (int i = 0; i < N; i++)
delete [] E[i];
delete [] E;
}

Далее создаю динамический массив, и заполняю его с обычного. После того мне надо динамический массив обратить
я это делаю так:

C++

float **xxt1 = new float*[kt];
for(int i = 0; i < kt; i++)
{
xxt1[i] = new float[kt];
}
 
for(i=0;i<kt;i++)
{for(j=0;j<kt;j++)
{xxt1[i][j]=xxt[i][j];}}
 
inversion(xxt1,kt);
 
for(i=0;i<kt;i++)
{for(j=0;j<kt;j++)
{StringGrid6->Cells[j+1][i+1]=FloatToStrF(xxt1[i][j],ffFixed,10,2);}}

но когда я вывожу оно мне неправильные значения выводит, точнее оно выводит такое, что просто жуть.

Помогите найти ошибку, а то послезавтра сдавать а программа не работает.

@volvo · 21.05.2013, 13:10

2garin_95, единственное, чего я в твоем коде не вижу - это чем заполняется созданный динамический массив. Ведь если там была "жуть", то после вызова inversion() там будет "инвертированная жуть"

Кстати, оригинальная функция (тестировалось на GCC) выдает при инвертировании вот такого массива:

C++

const int N = 5;
double data[N][N] =
{{6, 7, 23, 84, 15},
{21, 22, 23, 24, 25},
{31, 42, 53, 64, 75},
{41, 52, 63, 74, 85},
{10, 20, 30, 40, 50}
};

вот такие результаты:

Код

0.0648539  -2.25359e+014  -2.83143e+014  3.39001e+014  -3.8908e+013
-0.103614  2.71579e+014  3.82757e+014  -4.36224e+014  3.16556e+013
0.0214895  4.49058e+013  -1.2285e+013  -2.17473e+013  3.29449e+013
0.00844721  -3.11367e+012  8.8703e+012  -3.83775e+012  -5.22444e+012
0.00882332  -8.80124e+013  -9.61996e+013  1.22808e+014  -2.0468e+013

, так что... Делай выводы.

@2garin_95 · 21.05.2013, 13:43 **[ТС]**

Вот заполнения этого массива

C++

int i, j;
for(int i=0; i<kt; i++)
{for(int j=0; j<kt+1; j++){
if(j==0)
x[i][j]=1;
else
x[i][j]=pow(z[i],j);}}
 
for(i=0;i<kt;i++)
{for(j=0;j<kt;j++)
{StringGrid3->Cells[j+1][i+1]=FloatToStrF(x[i][j],ffFixed,5,2);}}
for(int i=0; i<kt+1; i++)
{for(int j=0; j<kt+1; j++){
xt[j][i]=x[i][j];}}
for(i=0;i<kt;i++)
{for(j=0;j<kt;j++)
{StringGrid4->Cells[j+1][i+1]=FloatToStrF(xt[i][j],ffFixed,5,2);}}
 
for(i=0;i<kt;i++)
{for(j=0;j<kt;j++)
xxt[i][j]=0;}
 
for(i=0;i<kt;i++)
{for(j=0;j<kt;j++)
{for(int k=0;k<kt-1;k++)
{xxt[i][j]=xxt[i][j]+(xt[i][k]*x[k][j]);}}}
 
for(i=0;i<kt;i++)
{for(j=0;j<kt;j++)
{StringGrid5->Cells[j+1][i+1]=FloatToStrF(xxt[i][j],ffFixed,10,2);}}

Массив Z задается рандомно от 0.01 до 5.99

@volvo · 21.05.2013, 14:05

С матрицей, определитель которой не нулевой, все нормально выводит:

C++

{{2, 9, 9, 4},
 {2, -3, 12, 8},
 {4, 8, 3, -5},
 {1, 2, 6, 4}};

, результат:

В функции обратной матрицы, не могу найти ошибку

Так что проблема - с твоей матрицей, которую ты хочешь инвертировать.

@2garin_95 · 22.05.2013, 10:04 **[ТС]**

UI, возможно, но в консоли обратную матрицу находит.

@volvo · 22.05.2013, 14:09

Сообщение от 2garin_95

в консоли обратную матрицу находит.

Приведи матрицу, от которой ищется обратная - посмотрим, что она там находит. Тот код, что ты привел 3-им постом, меня мало интересует, значения массива Z мне неизвестны.

@2garin_95 0 / 0 / 0 Регистрация: 19.05.2013 Сообщений: 7
	22.05.2013, 10:04 [ТС]	5
	UI, возможно, но в консоли обратную матрицу находит. 0

@volvo Супер-модератор 30599 / 19972 / 7801 Регистрация: 22.10.2011 Сообщений: 34,773 Записей в блоге: 6
	22.05.2013, 14:09	6
	Сообщение от 2garin_95 в консоли обратную матрицу находит. Приведи матрицу, от которой ищется обратная - посмотрим, что она там находит. Тот код, что ты привел 3-им постом, меня мало интересует, значения массива Z мне неизвестны. 0

Опции темы