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

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

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Найти индекс максимального элемента матрицы (Borland C++) http://www.cyberforum.ru/cpp-beginners/thread736202.html
1. Дано квадратную матрицу действительных чисел размерности 3. Вывести индексы максимального элемента.
C++ Работа с ASCII Что покажется на экране? ASCII* l 108 a 97 b 98 s 115 c 99 i 105 e 101 L 76 http://www.cyberforum.ru/cpp-beginners/thread736198.html
Отсортировать по возрастанию произведения отрицательных элементов строк матрицы C++
Помогите пожалуйста с задачей по 2-мерному массиву. Условия задачи: отсортировать по возрастанию произведения отрицательных элементов строк матрицы. Саму матрицу можно рандом делать... Добавлено через 1 час 47 минут апппппп
Найти номер последней из строк матрицы, содержащей максимальное кол-во одинаковых элементов C++
Дана матрица MxN, найти номер последней из ее строк содержащей максимальное кол-во одинаковых элементов! Помогите пожалуйста, ввести матрицу получается а вот что дальше не пойму #include <iostream> #include <stdlib.h> using namespace std; void main () { setlocale(LC_ALL, "Russian"); int a,n,m,i,j,S;
C++ Ввод-Вывод в файл http://www.cyberforum.ru/cpp-beginners/thread736181.html
Помогите написать 2 приложения, одно записывает введенный в консоле текст в файл .txt (объекты задаются в форме массива), второе выводит из .txt в консоль текст в двоичном коде (объекты задаются в форме списка)
C++ Ошибка кода Здравствуйте. Компилятор (Dev-C++) пишет ошибку, незнаю почему в коде (отрывок, есле нужен весь код выставлю весь) struct clientData_m = {clientData_m.accNum = "1", clientData_m.accNum = "12", clientData_m.accNum = "2", clientData_m.accNum = "23", clientData_m.accNum = "3", clientData_m.Name = "f", clientData_m.Name = "g", clientData_m.Name = "d", clientData_m.Name... подробнее

Показать сообщение отдельно
Mishaka
0 / 0 / 0
Регистрация: 25.09.2012
Сообщений: 32

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

18.12.2012, 01:18. Просмотров 1655. Ответов 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
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";
    }
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 13:19. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru