Форум программистов, компьютерный форум CyberForum.ru

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 33, средняя оценка - 4.79
Helo
4 / 4 / 0
Регистрация: 04.07.2013
Сообщений: 52
#1

Обратная матрица. Метод Гаусса—Жордана - C++

18.07.2013, 12:45. Просмотров 4108. Ответов 2
Метки нет (Все метки)

Здравствуйте. Пишу программку вычисления матриц и вот эта функция ( см. код ) компилируется и работает без ошибок, но считает не правильно. Кто хорошо знает математику и программирование подскажите, что я делаю не так?
Ссылка на вики обратная матрица http://ru.wikipedia.org/wiki/%D0%9E%...B8%D1%86%D0%B0
Код проблемной функции. Язык C++
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
void Invers_Array ( double ** & Matrix, const short int & size )
{
    double a, b, sum = 0;
    double ** Invers_Matrix = new double * [ size ];     /* Единичная матрица */
    for ( int i = 0; i < size; i ++ )
    {
        Invers_Matrix [i] = new double [ size ];
        for ( int j = 0; j < size; j ++ )
            Invers_Matrix [i][j] = 0;
            Invers_Matrix [i][i] = 1;
    }
    for ( int i = 0; i < size; i ++ )     /* Приведём матрицу к треугольной при помощи алгоритма Гаусса */
    {
        a = Matrix [i][i];
        for ( int j = i + 1; j < size; j ++ )
        {
            b = Matrix [j][i];
            for ( int k = 0; k < size; k ++ )
            {
                Matrix [j][k] = Matrix [i][k] * b - Matrix [j][k] * a;
                Invers_Matrix [j][k]= Invers_Matrix [i][k] * b - Invers_Matrix [j][k] * a;
            }
        }
    }
    for ( int i = 0; i < size; i ++ )     /* Обратный ход вычисления элементов обратной матрицы */
    {
        for ( int j = size - 1; j >= 0; j -- )
        {
            for ( int k = size - 1; k > j; k -- )
                sum += Matrix [j][k] * Invers_Matrix [k][i];
                if ( Matrix [j][j] == 0 )
                {
                    for ( int i = 0; i < size; i ++ )
                        delete [] Invers_Matrix [i];
                        delete [] Invers_Matrix;
                }
                Invers_Matrix [j][i] = ( Invers_Matrix [j][i] - sum ) / Matrix [j][j];
        }
    }
    printf ( "\n%s", " Обратная матрица \n\n" );
    Output_Array ( Invers_Matrix, size );
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
18.07.2013, 12:45     Обратная матрица. Метод Гаусса—Жордана
Посмотрите здесь:

C++ Обратная матрица на С
C++ Обратная матрица
C++ Обратная матрица
C++ Обратная матрица
C++ СЛАУ методом Жордана-Гаусса
C++ Метод гаусса жордана в с++
Вылет программы для вычисления СЛАУ методом Жордана-Гаусса при исполнении C++
C++ Метод Йордана-Гаусса. Как изменяется обратная матрица в силу таких преобразований?
Решение СЛУ методом Жордана-Гаусса, заменить массив на динамический C++
C++ Обратная матрица
C++ Метод Жордана - прокомментировать код программы
C++ обратная подстановка в методе Гаусса

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Aleator
 Аватар для Aleator
65 / 66 / 3
Регистрация: 04.05.2011
Сообщений: 1,163
18.07.2013, 12:58     Обратная матрица. Метод Гаусса—Жордана #2
Как именно не правильно? Выложи скрины.
Helo
4 / 4 / 0
Регистрация: 04.07.2013
Сообщений: 52
18.07.2013, 13:02  [ТС]     Обратная матрица. Метод Гаусса—Жордана #3
Цитата Сообщение от Aleator Посмотреть сообщение
Как именно не правильно? Выложи скрины.
Не сходится с ответом
ответ такой должен быть
-2 1
1.5 -0.5
Чем выше ранг матрицы тем больше расхождений с ответом
Миниатюры
Обратная матрица. Метод Гаусса—Жордана  
Yandex
Объявления
18.07.2013, 13:02     Обратная матрица. Метод Гаусса—Жордана
Ответ Создать тему
Опции темы

Текущее время: 03:13. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru