Вычислить для обратной матрицы сумму эллементов, расположенных над главной диагональю и под ней. С++

@Abysss · Регистрация: 09.01.2012

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

Программный код для формирования матрицы, а так же нахождения обратной матрицы

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
#include <conio.h>
#include <math.h>
#include <stdio.h>
#include <iostream>
 
using namespace std;
    int nn,i,j;
double A[8][8], L[8][8];
int main()
{
 
cout<<"Dimension of the matrix (<=8):"; 
cin>>nn;
cout<<"Enter the matrix of the element: "<<endl;
 
 
 
 for (i=0; i<nn; i++) 
   for (j=0; j<nn; j++)    {
       cout<<"\nElement:"<<i+1<<":"<<j+1<<endl;                                            
     float k; 
     cin>>k; 
     A[i][j] = k;
   }
 
 
 
 
cout<<"The formation of the identity matrix"<<endl;
 for (i=0; i<nn; i++) 
   for (j=0; j<nn; j++)
       if (i!=j) L[i][j]=0;      
   else L[i][j]=1; 
 j = 0;
 int m=0;
 int n;
 do    { 
 
          m = j; 
    for (n=j+1; n<nn; n++) 
      if (fabs(A[m][j])<fabs(A[n][j])) m = n;
        double prm=0;
 
 
 
    for (n=0; n<nn; n++)
   
 
    prm = A[j][j];
    for (n=0; n<nn; n++)
    { A[j][n] = double (A[j][n]/prm);
      L[j][n] = double (L[j][n]/prm);
    }
    for (i=0; i<=nn; i++)
      if (i!=j)
      { prm = A[i][j];
    for (n=j; n<=3; n++) A[i][n] = A[i][n]- double (prm*A[j][n]);
    for (n=0; n<=3; n++) L[i][n] = L[i][n]- double (prm*L[j][n]);
      }
    j++;
    }
  while (j<nn);
 
 
 
 cout<<"\nInverse matrix:\n";
 for (i=0; i<nn; i++)
   for (j=0; j<nn; j++)
   { if (j==0) cout<<"\n";  
            else cout<<"\t";
     printf("%f", L[i][j]);
   }

Комментарий модератора

Используйте теги форматирования кода!

@BumerangSP · 09.01.2012, 16:58

Сделал немного по-своему:

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
#include <conio.h>
#include <math.h>
#include <stdio.h>
#include <iostream>
#include <time.h>
#include <conio.h> 
using namespace std;
        int nn,i,j;
double A[8][8], L[8][8];
int main()
{
 srand(time(NULL));
cout<<"Dimension of the matrix (<=8):"; 
cin>>nn;
 
 
 
 
 for (i=0; i<nn; i++) 
   for (j=0; j<nn; j++)    
   {
           cout<<"\nElement:"<<i+1<<":"<<j+1<<endl;                                                                                        
 
   
     A[i][j] = rand()%10;
     printf("%f",A[i][j]);
   }
 
 
 
 
cout<<"The formation of the identity matrix"<<endl;
 for (i=0; i<nn; i++) 
   for (j=0; j<nn; j++)
           if (i!=j) L[i][j]=0;          
   else L[i][j]=1; 
 j = 0;
 int m=0;
 int n;
 do    { 
 
              m = j; 
    for (n=j+1; n<nn; n++) 
      if (fabs(A[m][j])<fabs(A[n][j])) m = n;
            double prm=0;
 
 
 
    for (n=0; n<nn; n++)
   
 
    prm = A[j][j];
    for (n=0; n<nn; n++)
    { A[j][n] = double (A[j][n]/prm);
      L[j][n] = double (L[j][n]/prm);
    }
    for (i=0; i<=nn; i++)
      if (i!=j)
      { prm = A[i][j];
        for (n=j; n<=3; n++) A[i][n] = A[i][n]- double (prm*A[j][n]);
        for (n=0; n<=3; n++) L[i][n] = L[i][n]- double (prm*L[j][n]);
      }
    j++;
    }
  while (j<nn);
 
 
 
 cout<<"\nInverse matrix:\n";
 for (i=0; i<nn; i++)
   for (j=0; j<nn; j++)
   { 
    if (j==0) 
     cout<<"\n";  
      else cout<<"\t";
     printf("%5.2f", L[i][j]);
   }
  
  
   printf("\n");
   printf("\n");
   float sum1=0;      
    for (i=0; i<nn; i++)
   {     
    for (j=0; j<nn; j++) 
     {
                               if (j<i) 
                               {
                               sum1=sum1+L[i][j];
      printf("%5.2f ",L[i][j]);   //Вывод элементов под диагональю
      }
      }
    printf("\n");  
   }
   printf("\n");
   printf("\n");
   float sum2=0;        
    for (i=0; i<nn; i++)
   {     
    for (j=0; j<nn; j++) 
     {
                               if (j>i) 
                               {
                               sum2=sum2+L[i][j];
      printf("%5.2f ",L[i][j]);  //Вывод элементов над диагональю
      }
      }
    printf("\n");  
   }
printf("\nPod diagonal`y: %5.2f",sum1);
printf("\nNad diagonal`y: %5.2f",sum2);     
getch();
return 0;
}

Элементы вводятся рандомно, и округление стоит до сотых.

@Abysss · 09.01.2012, 17:09 **[ТС]**

ММм Спасибо огромное))))) Осталось только переписать на С++.
А не подскажете ещё как вывести на печть треугольную матрицу, имеющую большую сумму?? то есть необходимо сравнить полученные результаты над и под диагональю, и вывести на печать треугольную матрицу, сумма эллементов которой больше.

Новые блоги и статьи Все статьи Все блоги /
Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка SDL3, Box2D, FreeType и SDL3_ttf из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .	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 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.

@Abysss 0 / 0 / 0 Регистрация: 09.01.2012 Сообщений: 6
	09.01.2012, 17:09 [ТС]
	ММм Спасибо огромное))))) Осталось только переписать на С++. А не подскажете ещё как вывести на печть треугольную матрицу, имеющую большую сумму?? то есть необходимо сравнить полученные результаты над и под диагональю, и вывести на печать треугольную матрицу, сумма эллементов которой больше. 0