Сравнение двух методов и скорости нахождения решения СЛАУ с комплексными переменными

@ashchukin95 · Регистрация: 17.05.2015

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

Есть 2 метода решения СЛАУ с комплексными переменными. Нужно сравнить скорости нахождения решения двух этих методов между собой. Подскажите, как это можно сделать, если программы написаны в Паскале?

Pascal
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
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
USES
  Crt;
  
TYPE
  Complex=RECORD
            re,im: REAL;
          END;
  MATRIX=ARRAY[1..50, 1..51] OF Complex;
  ROW=ARRAY[1..51] OF Complex;
  
VAR
  i, j, n, m, k, kMax: INTEGER;
  max, r, p: Complex; 
  A: MATRIX;
  AMax: ROW;
 
PROCEDURE Summa(f,g:Complex; var z:complex);
  BEGIN
    z.re:=f.re+g.re;
    z.im:=f.im+g.im;
  END; 
 
PROCEDURE Raznost(f,g:complex; var z:complex);
  BEGIN
    z.re:=f.re-g.re;
    z.im:=f.im-g.im;
  END; 
 
PROCEDURE Proizvedenie(f,g:complex; var z:complex);
  BEGIN
    z.re:=f.re*g.re-f.im*g.im;
    z.im:=g.re*f.im+f.re*g.im;
  END;
 
PROCEDURE Chastnoe(f,g:complex; var z:complex);
  BEGIN
    z.re:=(f.re*g.re+f.im*g.im)/(sqr(g.re)+sqr(g.im));
    z.im:=(f.im*g.re-f.re*g.im)/(sqr(g.re)+sqr(g.im));
  END;
  
FUNCTION Modul(f:complex): REAL;
  BEGIN
    Modul:=sqrt(sqr(f.re)+sqr(f.im));
  END;
 
BEGIN
  ClrScr;
 
// ВВОД МАТРИЦЫ
  WriteLn('Введите размерность матрицы:');
  ReadLn(n);
  WriteLn;
  
  WriteLn('Введите вещественную часть матрицы:');
  FOR i:= 1 TO n DO
    BEGIN
      m:=0;
      FOR j:= 1 TO n+1 DO
        BEGIN
          IF j<n+1 THEN
            BEGIN
              GoToXY(j+m, i+4);
              Read(A[i,j].re);
              m:=m+3;
            END
          ELSE
            BEGIN
              GoToXY(j+m, i+4);
              Write('||  ');
              Read(A[i,j].re);
            END;
        END;
    END; 
  
  WriteLn('Введите мнимую часть матрицы:');
  FOR i:= 1 TO n DO
    BEGIN
      m:=0;
      FOR j:= 1 TO n+1 DO
        BEGIN
          IF j<n+1 THEN
            BEGIN
              GoToXY(j+m, i+5+n);
              Read(A[i,j].im);
              m:=m+3;
            END
          ELSE
            BEGIN
              GoToXY(j+m, i+5+n);
              Write('||  ');
              Read(A[i,j].im);
            END;
        END;
    END;
 
// МЕТОД ГАУССА
  FOR k:=1 TO n DO
    BEGIN
      //поиск максимального
      kMax:=k;
      max.re:=A[k,k].re;
      max.im:=A[k,k].im;
      FOR i:=k TO n DO
        IF Modul(A[i,k])>Modul(max) THEN
          BEGIN
            max:=A[i,k];
            kMax:=i;
          END;
        IF (max.re=0) and (max.im=0) THEN
          BEGIN
            WriteLn('Матрица вырождена!');
            Halt;
          END;
      FOR i:=1 TO n+1 DO
        BEGIN
          AMax[i]:=A[kMax, i];
          A[kMax, i]:=A[k, i];
          A[k, i]:=AMax[i];
        END;
        
      //алгоритм Гаусса
      FOR j:=k+1 TO n DO
        BEGIN
          Chastnoe(A[j,k], A[k,k], r);
          FOR i:=k TO n+1 DO
            BEGIN
              Proizvedenie(A[k,i], r, p);
              Raznost(A[j,i], p, A[j,i]);
            END;
        END;
      FOR j:=1 TO k-1 DO
        BEGIN
          Chastnoe(A[j,k], A[k,k], r);
          FOR i:=k TO n+1 DO
            BEGIN
              Proizvedenie(A[k,i], r, p);
              Raznost(A[j,i], p, A[j,i]);
            END;
        END;
    END;
  FOR i:=1 TO n DO
    BEGIN
      r:=A[i,i];
      FOR j:=i TO n+1 DO
        Chastnoe(A[i,j], r, A[i,j]);
    END;
 
//ВЫВОД МАТРИЦЫ
  WriteLn;
  WriteLn('Вещественная часть диагональной матрицы:');
  FOR i:=1 TO n DO
    BEGIN
      FOR j:=1 TO n DO
        Write(A[i,j].re:5:2, ' ');
      Write('||');
      Write( A[i,n+1].re:5:2);
      WriteLn;
    END;
    
  WriteLn;   
  WriteLn('Мнимая часть диагональной матрицы:');
  FOR i:=1 TO n DO
    BEGIN
      FOR j:=1 TO n DO
        Write(A[i,j].im:5:2, ' ');
      Write('||');
      Write( A[i,n+1].im:5:2);
      WriteLn;
    END;
 
  WriteLn;
  FOR i:=1 TO n DO
    BEGIN
     IF A[i,n+1].im > 0 THEN
       WriteLn('x[',i,']=',A[i,n+1].re:4:2,'+i',A[i,n+1].im:4:2);
     IF A[i,n+1].im < 0 THEN
       WriteLn('x[',i,']=',A[i,n+1].re:4:2,'-i',abs(A[i,n+1].im):4:2);
     IF A[i,n+1].im = 0 THEN
       WriteLn('x[',i,']=',A[i,n+1].re:4:2);
    END;
  WriteLn;
    
END.

Pascal
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
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
USES
 
  Crt;
TYPE
  MATRIX=ARRAY[1..50, 1..51] OF REAL;
  ROW=ARRAY[1..51] OF REAL;
VAR
  i, j, n, m, k, kMax: INTEGER;
  max, r: REAL;
  A, A1, A2: MATRIX;
  AMax: ROW;
BEGIN
  ClrScr;
 
// ВВОД МАТРИЦЫ
  WriteLn('Введите размерность матрицы:');
  ReadLn(n);
  
  WriteLn;
  WriteLn('Введите вещественную часть матрицы:');
  FOR i:= 1 TO n DO
    BEGIN
      m:=0;
      FOR j:= 1 TO n+1 DO
        BEGIN
          IF j<n+1 THEN
            BEGIN
              GoToXY(j+m, i+4);
              Read(A1[i,j]);
              m:=m+3;
            END
          ELSE
            BEGIN
              GoToXY(j+m, i+4);
              Write('||  ');
              Read(A1[i,j]);
            END;
        END;
    END;
    
 WriteLn('Введите мнимую часть матрицы:');
  FOR i:= 1 TO n DO
    BEGIN
      m:=0;
      FOR j:= 1 TO n+1 DO
        BEGIN
          IF j<n+1 THEN
            BEGIN
              GoToXY(j+m, i+5+n);
              Read(A2[i,j]);
              m:=m+3;
            END
          ELSE
            BEGIN
              GoToXY(j+m, i+5+n);
              Write('||  ');
              Read(A2[i,j]);
            END;
        END;
    END;
 
//Создание матрицы A
   FOR i:=1 TO n DO
     FOR j:=1 TO n DO
       BEGIN
         A[i,j]:=A1[i,j];
         A[n+i,n+j]:=A1[i,j];
         A[n+i,j]:=A2[i,j];
         A[i,n+j]:=-A2[i,j];
       END;
   FOR i:=1 TO n DO
     BEGIN
       A[i,n*2+1]:=A1[i,n+1];
       A[n+i,n*2+1]:=A2[i,n+1];
     END;   
 
// МЕТОД ГАУССА
  FOR k:=1 TO n*2 DO
    BEGIN
      //поиск максимального
      kMax:=k;
      max:=A[k,k];
      FOR i:=k TO n*2 DO
        IF abs(A[i,k])>abs(max) THEN
          BEGIN
            max:=A[i,k];
            kMax:=i;
          END;
        IF max=0 THEN
          BEGIN
            WriteLn('Матрица вырождена!');
            Halt;
          END;
      FOR i:=1 TO n*2+1 DO
        BEGIN
          AMax[i]:=A[kMax, i];
          A[kMax, i]:=A[k, i];
          A[k, i]:=AMax[i];
        END;
      //алгоритм Гаусса
      FOR j:=k+1 TO n*2 DO
        BEGIN
          r:=A[j,k]/A[k,k];
          FOR i:=k TO n*2+1 DO
            A[j,i]:=A[j,i]-A[k,i]*r;
        END;
      FOR j:=1 TO k-1 DO
        BEGIN
          r:=A[j,k]/A[k,k];
          FOR i:=k TO n*2+1 DO
            A[j,i]:=A[j,i]-A[k,i]*r;
        END;
    END;
  FOR i:=1 TO n*2 DO
    BEGIN
      r:=A[i,i];
      FOR j:=i TO n*2+1 DO
        A[i,j]:=A[i,j]/r;
    END;
 
//ВЫВОД МАТРИЦЫ
WriteLn;
WriteLn('Диагональная матрица:');
  FOR i:=1 TO n*2 DO
    BEGIN
      FOR j:=1 TO n*2 DO
        Write(A[i,j]:5:2, ' ');
      Write('||');
      Write( A[i,n*2+1]:5:2);
      WriteLn;
    END;
END.

@Puporev · 17.06.2015, 11:40

Сообщение от ashchukin95

если программы написаны в Паскале?

Смотря в каком Паскале, для АВС.net Вам ответили в двух темах.

Новые блоги и статьи Все статьи Все блоги /
Паттерны в Python: Singleton, Factory и Observer py-thonny 26.04.2025 Паттерны проектирования — это проверенные временем решения типовых проблем разработки программного обеспечения. Их история берёт начало с книги "Приёмы объектно-ориентированного проектирования. . . .	Исключения в C#: Stack Overflow, Access Violation и Out of memory stackOverflow 26.04.2025 Исключения в C# — это не только механизм оповещения о проблемах, а целое искусство управления потоком выполнения программы в экстремальных ситуациях. Обычное исключение, например,. . .	Логирование в C# ASP.NET Core с помощью Serilog, ElasticSearch, Kibana stackOverflow 25.04.2025 Помните те времена, когда для анализа проблемы приходилось подключаться к серверу, искать нужный лог-файл среди десятков других и вручную фильтровать тысячи строк в поисках ошибки? К счастью, эти дни. . .	Структура "железный OnKeyUp" вместо антидребезга. Полностью асинхронный счётчик. Hrethgir 25.04.2025 Программа для симуляции схемы - Logisim Evolution В общем какое-то время отвлёкся, так было надо, теперь когда запилю это на verilog и FPGA , досоставлю заявку в ФИПС на полезную модель - не готов. . .	Автоматизация Amazon Web Services (AWS) с Boto3 в Python py-thonny 25.04.2025 Облачные вычисления стали неотъемлемой частью современной ИТ-инфраструктуры, а Amazon Web Services (AWS) занимает лидирующие позиции среди провайдеров облачных услуг. Управление многочисленными. . .
Apache Kafka vs RabbitMQ в микросервисной архитектуре ArchitectMsa 25.04.2025 Современная разработка ПО всё чаще склоняется к микросервисной архитектуре — подходу, при котором приложение разбивается на множество небольших, автономных сервисов. В этой распределённой среде. . .	Параллельное программирование с OpenMP в C++ NullReferenced 24.04.2025 Параллельное программирование — подход к созданию программ, когда одна задача разбивается на несколько подзадач, которые могут выполняться одновременно. Оно стало необходимым навыком для. . .	Цепочки методов в C# с Fluent API UnmanagedCoder 24.04.2025 Современное программирование — это не только решение функциональных задач, но и создание кода, который удобно поддерживать, расширять и читать. Цепочки методов и Fluent-синтаксис в C# стали мощным. . .	Мульти-тенантные БД с PostgreSQL Row Security Codd 23.04.2025 Современные облачные сервисы и бизнес-приложения всё чаще обслуживают множество клиентов в рамках единой программной инфраструктуры. Эта архитектурная модель, известная как мульти-тенантность, стала. . .	Реализация конвейеров машинного обучения с Python и Scikit-learn AI_Generated 23.04.2025 Мир данных вокруг нас растёт с каждым днём, и умение эффективно обрабатывать информацию стало необходимым навыком. Специалисты по машинному обучению ежедневно сталкиваются с задачами предобработки. . .

Опции темы