Найти обратную матрицу

@KillerMak · Регистрация: 23.09.2015

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

Ребят немогу в конда сделал провереку, ошибку выдаёт, не могу найти ошибку помогите!!

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
#include <iostream>
#include <cmath>
#include <algorithm>
#include <vector>
#include <map>
#include <stack>
#include <deque>
#include <set>
#include <string>
#include <limits>
#include <fstream>
 
using namespace std;
      // Алгорит обращения 
vector < vector <double> > inverse(vector < vector <double> > a) {
    int  n = a.size();
    vector < vector <double> > ans(n, vector <double> (n, 0));  
    for (int i = 0; i < n; i++){
        ans[i][i] = 1.0;
    }  
    //Вычисление матрицы
    for (int i = 0; i < n; i++){
        int row = i;
        double mx = a[i][i];
        for(int k = i + 1; k < n; k++){
            if (abs(a[k][i]) > mx){
                row = k;
                mx = abs(a[k][i]);
            }
        }
        if (row != i) {
            swap(a[row], a[i]);
            swap(ans[row], ans[i]);
        }
        for (int j = i+1; j < n; j++){
            double e = a[j][i]/a[i][i];
            for (int k = 0; k < n; k++){
                a[j][k] -= e*a[i][k];
                ans[j][k] -= e*ans[i][k];
            }
        }
    }
    // Обратная матрица 
    for (int i = n - 1; i >= 0; i--) {
        for (int j = i - 1; j >= 0; j--){
            double e = a[j][i]/a[i][i];
            for (int k = 0; k < n; k++){
                a[j][k] -= e*a[i][k];
                ans[j][k] -= e*ans[i][k];
            }
        }
        for (int j = 0; j < n; j++) {
            ans[i][j] /= a[i][i];
        }
    }
    return ans;
}
 
 
int main(){
    // Ввод исходной матрици 
    setlocale (LC_ALL, "RUS");
    int matrix, k;
    int c;
    cout << "Введите размер матрицы:\n" << endl;
    cin >> matrix;                
    vector < vector <double> > a(matrix, vector <double> (matrix, 0)), ans;
 cout << "Введите матрицу:\n" << endl;
    for (int i = 0; i < matrix; i++){
        for (int j = 0; j < matrix; j++){
            scanf("%lf", &a[i][j]);
        }
    } 
    ans = inverse(a);
    // Ввывод результата
 cout << "Обратная матрица:\n" << endl;
    for (int i = 0; i < matrix; i++) {
        for (int j = 0; j < matrix; j++) {
            printf("%.15lf ", ans[i][j]);
        }
        puts("");
  //Проверка
        int H;
        for (int i = 0; i < matrix; i++){   //умножение матриц начальную на полученую должны получить еденичную
        for (int j = 0; j < matrix; j++){
            H[i][j] = 0;
            for (int t = 0; t < n; t++){
                H[i][j] += a[i][j] * ans[i][j]);
                printf("%1.2f", Ans[i][j]); 
             cout << " " ;
            }
        }
    }
 
    return 0;
}

@zss · 30.11.2015, 16:32

Ошибки компиляции я исправил:

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
106
107
108
109
110
111
112
113
114
115
#include <iostream>
#include <cmath>
//#include <algorithm>
#include <vector>
//#include <map>
//#include <stack>
//#include <deque>
//#include <set>
//#include <string>
//#include <limits>
//#include <fstream>
 
using namespace std;
// Алгорит обращения 
vector < vector <double> > inverse(vector < vector <double> > a) 
{
    int  n = a.size();
    vector < vector <double> > ans(n, vector <double> (n, 0));  
    for (int i = 0; i < n; i++)
    {
        ans[i][i] = 1.0;
    }  
    //Вычисление матрицы
    for (int i = 0; i < n; i++)
    {
        int row = i;
        double mx = a[i][i];
        for(int k = i + 1; k < n; k++)
        {
            if (abs(a[k][i]) > mx)
            {
                row = k;
                mx = abs(a[k][i]);
            }
        }
        if (row != i) 
        {
            swap(a[row], a[i]);
            swap(ans[row], ans[i]);
        }
        for (int j = i+1; j < n; j++)
        {
            double e = a[j][i]/a[i][i];
            for (int k = 0; k < n; k++)
            {
                a[j][k] -= e*a[i][k];
                ans[j][k] -= e*ans[i][k];
            }
        }
    }
    // Обратная матрица 
    for (int i = n - 1; i >= 0; i--) 
    {
        for (int j = i - 1; j >= 0; j--)
        {
            double e = a[j][i]/a[i][i];
            for (int k = 0; k < n; k++)
            {
                a[j][k] -= e*a[i][k];
                ans[j][k] -= e*ans[i][k];
            }
        }
        for (int j = 0; j < n; j++) {
            ans[i][j] /= a[i][i];
        }
    }
    return ans;
}
 
 
int main()
{
    // Ввод исходной матрици 
    setlocale (LC_ALL, "RUS");
    int matrix;
    //int c;
    cout << "Введите размер матрицы:\n" << endl;
    cin >> matrix;                
    vector < vector <double> > a(matrix, vector <double> (matrix, 0)), ans,H;
    cout << "Введите матрицу:\n" << endl;
    for (int i = 0; i < matrix; i++)
    {
        for (int j = 0; j < matrix; j++)
        {
            scanf("%lf", &a[i][j]);
        }
    } 
    ans = inverse(a);
    // Ввывод результата
    cout << "Обратная матрица:\n" << endl;
    for (int i = 0; i < matrix; i++) 
    {
        for (int j = 0; j < matrix; j++) 
        {
            printf("%.15lf ", ans[i][j]);
        }
        puts("");
        //Проверка
        //int H;
        for (int i = 0; i < matrix; i++)
        {   //умножение матриц начальную на полученую должны получить еденичную
            for (int j = 0; j < matrix; j++)
            {
                H[i][j] = 0;
                for (int t = 0; t < 3/*n*/; t++) // что такое n? 
                {
                    H[i][j] += a[i][j] * ans[i][j];
                    printf("%1.2f",ans[i][j]); 
                    cout << " " ;
                }
            }
        }
    }
    return 0;
}

Дальше работайте с помощью отладчика
https://www.cyberforum.ru/cpp-... 62479.html

Новые блоги и статьи Все статьи Все блоги /
Символьное дифференцирование igorrr37 13.02.2026 / * Программа принимает математическое выражение в виде строки и выдаёт его производную в виде строки и вычисляет значение производной при заданном х Логарифм записывается как: (x-2)log(x^2+2) -. . .	Камера 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. Пошагово создадим проект для загрузки изображения. . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL3_image 8Observer8 10.02.2026 Содержание блога Библиотека SDL3_image содержит инструменты для расширенной работы с изображениями. Пошагово создадим проект для загрузки изображения формата PNG с альфа-каналом (с прозрачным. . .	Установка Qt-версии Lazarus IDE в Debian Trixie Xfce volvo 10.02.2026 В общем, достали меня глюки IDE Лазаруса, собранной с использованием набора виджетов Gtk2 (конкретно: если набирать текст в редакторе и вызвать подсказку через Ctrl+Space, то после закрытия окошка. . .