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

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

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 13, средняя оценка - 4.62
Mishaka
0 / 0 / 0
Регистрация: 25.09.2012
Сообщений: 32
18.12.2012, 01:18     Обратная матрица методом Гаусса-Жордана #1
Хотел получить обратную матрицу проводимостей которая является частью программы, которая считает токи в ветвях методом узловых напряжений. За основу взял код который вычисляет обратную матрицу методом Гаусса-Жордана.
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
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
#include <iostream>
 
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 << "Enter 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 << "Enter matrix[" << 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;
 
    std::cin.get();
    return 0;
}
Переделал под комплексные числа и матрица стала вычислять что зря, ничего не сходится!
Пожалуйста проверьте, где я не так написал и где ошибка...
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
43
44
45
46
47
48
49
50
51
52
for(int i=0;i<n;i++){                     //заполняется единичная матрица
        for(int j=0;j<n;j++){
            if (i==j){Ed[i][j]=std::complex<double>(1,0);}
            else {Ed[i][j]=std::complex<double>(0,0);}
        }
    }       
        for (int k = 0; k < n; k++)  //обратная матрица
    {
        temp = std::complex<double>(AYT[k][k]);
        for (int j=0;j<n;j++)
        {
            AYT[k][j]=std::complex<double>(AYT[k][j]/temp);
            Ed[k][j]=std::complex<double>(Ed[k][j]/temp);
        }
 
        for (int i=k+1;i<n;i++)
        {
            temp = std::complex<double>(AYT[i][k]);
            for (int j=0;j<n;j++)
            {
                AYT[i][j]=std::complex<double>(AYT[i][j]-(AYT[k][j]*temp));
                Ed[i][j]=std::complex<double>(Ed[i][j]-(Ed[k][j]*temp));
            }
        }
    }
 
    for (int k=n-1;k>0;k--)
    {
        for (int i=k-1;i>=0;i--)
        {
            temp = std::complex<double>(AYT[i][k]);
            for (int j=0;j<n;j++)
            {
                AYT[i][j]=std::complex<double>(AYT[i][j]-(AYT[k][j]*temp));
                Ed[i][j]=std::complex<double>(Ed[i][j]-(Ed[k][j]*temp));
            }
        }
    }
 
    for (int i=0;i<n;i++){
        for (int j=0;j<n;j++){
            AYT[i][j] = std::complex<double>(Ed[i][j]);}
    }
 
       
    cout << "OBRATNAYA!!!!!!!!!!!!!!!!!!!!!!\n\n\n\n";
    for (int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            cout<<AYT[i][j]<<"\t";
        }
        cout<<"\n";
    }
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
18.12.2012, 01:18     Обратная матрица методом Гаусса-Жордана
Посмотрите здесь:

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

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

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