Обратная матрица

@Savelia · Регистрация: 22.11.2010

Студворк — интернет-сервис помощи студентам

Срочно в СИ нужна программа нахождения обратной матрицы. Что-то я пробовала методом Гаусса у меня на первой же функции ошибка, расскажите какие еще есть альтернативные способы?

silent_1991 · 11.01.2011, 18:31

Вернее всего Гаусс. И просто, и достаточно надёжно. Что за ошибка-то?

@sandye51 · 11.01.2011, 18:48

Не по теме:

silent_1991, деление на 0 наверно))

silent_1991 · 11.01.2011, 18:50

Не по теме:

sandye51, ну так строки-то надо не забывать переставлять при необходимости)))

@sandye51 · 11.01.2011, 18:51

Не по теме:

silent_1991, особенно если матрица вырожденная;)

silent_1991 · 11.01.2011, 18:53

Не по теме:

sandye51, ну а это вообще святое - проверить детерминант перед обращением)))

@Savelia · 11.01.2011, 18:54 **[ТС]**

там код какой-то путанный, я тут по форуму полазила и нашла ваш код, можно его в СИ переделать, там ведь нужно использовать malloc, не поможете мне?

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;
}

silent_1991 · 11.01.2011, 19:05

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
104
105
#include <stdio.h>
#include <stdlib.h>
 
void inversion(double **A, int N)
{
    double temp;
    int i, j, k;
 
    double **E = (double **)malloc(N * sizeof(double *));
 
    for (i = 0; i < N; i++)
        E[i] = (double *)malloc(N * sizeof(double));
 
    for (i = 0; i < N; i++)
        for (j = 0; j < N; j++)
        {
            E[i][j] = 0.0;
 
            if (i == j)
                E[i][j] = 1.0;
        }
 
    for (k = 0; k < N; k++)
    {
        temp = A[k][k];
 
        for (j = 0; j < N; j++)
        {
            A[k][j] /= temp;
            E[k][j] /= temp;
        }
 
        for (i = k + 1; i < N; i++)
        {
            temp = A[i][k];
 
            for (j = 0; j < N; j++)
            {
                A[i][j] -= A[k][j] * temp;
                E[i][j] -= E[k][j] * temp;
            }
        }
    }
 
    for (k = N - 1; k > 0; k--)
    {
        for (i = k - 1; i >= 0; i--)
        {
            temp = A[i][k];
 
            for (j = 0; j < N; j++)
            {
                A[i][j] -= A[k][j] * temp;
                E[i][j] -= E[k][j] * temp;
            }
        }
    }
 
    for (i = 0; i < N; i++)
        for (j = 0; j < N; j++)
            A[i][j] = E[i][j];
 
    for (int i = 0; i < N; i++)
        free(E[i]);
 
    free(E);
}
 
int main()
{
    int N;
    int i, j;
 
    printf("Enter N: ");
    scanf("%d", &N);
 
    double **matrix = (double **)malloc(N * sizeof(double *));
 
    for (i = 0; i < N; i++)
        matrix[i] = (double *)malloc(N * sizeof(double));
 
    for (i = 0; i < N; i++)
        for (j = 0; j < N; j++)
        {
            printf("Enter matrix[%d][%d] = ", i, j);
            scanf("%lf", &matrix[i][j]);
        }
 
    inversion(matrix, N);
 
    for (i = 0; i < N; i++)
    {
        for (j = 0; j < N; j++)
            printf("%lf  ", matrix[i][j]);
 
        printf("\n");
    }
 
    for (i = 0; i < N; i++)
        free(matrix[i]);
 
    free(matrix);
 
    return 0;
}

@sandye51 · 11.01.2011, 19:07

Savelia,

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
#include <stdio.h>
#include <stdlib.h>
 
void inversion(double **A, int N)
{
    double temp;
 
    double **E = (double**)malloc(sizeof(double*) * N);
 
    for (int i = 0; i < N; i++)
        E[i] = (double*)malloc(sizeof(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++)
        free(E[i]);
    free(E);
}
 
int main()
{
    int N;
 
    printf("Enter N: ");
    scanf("%d", &N);
 
    double **matrix = (double **)malloc(N * sizeof(double*));
 
    for (int i = 0; i < N; i++)
        matrix[i] = (double *)malloc(N * sizeof(double));
 
    for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
        {
            printf("Enter matrix[%d][%d] = ", i, j);
            scanf("%lf", &matrix[i][j]);
        }
 
    inversion(matrix, N);
 
    for (int i = 0; i < N; i++)
    {
        for (int j = 0; j < N; j++)
            printf("%lf  ", matrix[i][j]);
 
        printf("\n");
    }
 
    for (int i = 0; i < N; i++)
        free(matrix[i]);
     free(matrix);
 
    return 0;
}

Не по теме:

silent_1991, опередил)

@Savelia · 11.01.2011, 19:09 **[ТС]**

премного благодарна, теперь буду подробнее разбираться)

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .	Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .
И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»	SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 12.02.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .

@Savelia 1 / 1 / 1 Регистрация: 22.11.2010 Сообщений: 42

	Обратная матрица 11.01.2011, 18:07. Показов 2321. Ответов 9 Метки нет (Все метки) Срочно в СИ нужна программа нахождения обратной матрицы. Что-то я пробовала методом Гаусса у меня на первой же функции ошибка, расскажите какие еще есть альтернативные способы? 0

silent_1991 5058 / 3118 / 271 Регистрация: 11.11.2009 Сообщений: 7,044
	11.01.2011, 18:31
	Вернее всего Гаусс. И просто, и достаточно надёжно. Что за ошибка-то? 0

@sandye51
	11.01.2011, 18:48
	Не по теме: silent_1991, деление на 0 наверно)) 0

silent_1991
	11.01.2011, 18:50
	Не по теме: sandye51, ну так строки-то надо не забывать переставлять при необходимости))) 0

@sandye51
	11.01.2011, 18:51
	Не по теме: silent_1991, особенно если матрица вырожденная;) 0

silent_1991
	11.01.2011, 18:53
	Не по теме: sandye51, ну а это вообще святое - проверить детерминант перед обращением))) 0

@Savelia 1 / 1 / 1 Регистрация: 22.11.2010 Сообщений: 42
	11.01.2011, 19:09 [ТС]
	премного благодарна, теперь буду подробнее разбираться) 0