Вычисление собственных значений и векторов матриц ---ошибка в коде - не правильно работает программа

@xamelione25 · Регистрация: 04.02.2013

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

Помогите подкорректировать код:

Имеется задание:
1. Локализовать собственные значения матрицы $https://www.cyberforum.ru/cgi-bin/latex.cgi?\begin{vmatrix}5 & -1 & 1\\ 6 & 0 & 2\\ -3 & 1 & -1\end{vmatrix}$ .
2. Определить точные значения собственных чисел и векторов матрицы.
3. Составить блок-схему и программу для вычисления с заданной точностью наибольшего по абсолютной величине собственного числа и соответствующего собственного вектора, используя степенной метод (итерации)

Кликните здесь для просмотра всего текста

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
type vector=array[1..3] of real;
     matrix=array[1..3,1..3] of real;
function sp(const x0,x1:vector; const n:byte):real;
var i:byte; sum:real;
begin 
sum:=0;
for i:=1 to n do 
sum:=sum+x0[i]*x1[i];
sp:=sum;
end;
procedure normir(var x1:vector; const n:byte);
var i:byte; sum:real;
begin 
for i:= 1 to n do
sum:=sum+x1[i]*x1[i];
for i:= 1 to n do
x1[i]:=x1[i]/sqrt(sum);
end;
function norma(const x0,x1:vector; const n:byte):real;
var i:byte; x:vector; sum:real;
begin
for i:=1 to n do
x[i]:=x1[i]-x0[i];
for i:=1 to n do
sum:=sum+sqr(x[i]);
norma:=sqrt(sum);
end;
procedure matproz(const a:matrix; const x0:vector; var x1:vector; const n:byte);
var i,j:byte;
begin
for i:=1 to n do
x1[i]:=0;
for i:=1 to n do
for j:=1 to n do
x1[i]:=x1[i]+x0[j]*a[i,j];
end;
var x0,x1:vector;
a:matrix; l0,l1,e:real;
n,i,j:byte;
begin 
n:=3;
for i:=1 to n do
x0[i]:=1/sqrt(n);
x0[1]:=1;
writeln('Vvedite matricu');
  for i := 1 to n do
    for j := 1 to n do
    begin
      write('a[', i, ',', j, ']=');
      readln(a[i, j]);
    end;
  matproz(a,x0,x1,n);
    l1:=sp(x0,x1,n);
  normir(x1,n);
repeat
  x0:=x1;
  l0:=l1;
  matproz(a,x0,x1,n);
  l1:=sp(x0,x1,n);
 normir(x1,n);
 
  until(abs(l1-l0)<e);
  writeln(x1[1]:5:2,' ',x1[2]:5:2,' ',x1[3]:5:2);
   writeln(l1);
end.

Итог в программе должен быть такой как на скрине в соответствии с пунктами 1 и 2 из задания. Но у меня по какой-то причине не происходит запрос порядка матрицы (возможно он в коде задается уже) и ни каких дальнейших ответов после ввода значений матрицы как на скрине не происходит.

@xamelione25 · 29.12.2018, 18:28 **[ТС]**

Здесь я уже сделал ручной ввод ранга матрицы, но итоговые значения почему-то все также не вылазят

Кликните здесь для просмотра всего текста

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
type vector=array[1..3] of real;
     matrix=array[1..3,1..3] of real;
function sp(const x0,x1:vector; const n:byte):real;
var i:byte; sum:real;
begin 
sum:=0;
for i:=1 to n do 
sum:=sum+x0[i]*x1[i];
sp:=sum;
end;
procedure normir(var x1:vector; const n:byte);
var i:byte; sum:real;
begin 
for i:= 1 to n do
sum:=sum+x1[i]*x1[i];
for i:= 1 to n do
x1[i]:=x1[i]/sqrt(sum);
end;
function norma(const x0,x1:vector; const n:byte):real;
var i:byte; x:vector; sum:real;
begin
for i:=1 to n do
x[i]:=x1[i]-x0[i];
for i:=1 to n do
sum:=sum+sqr(x[i]);
norma:=sqrt(sum);
end;
procedure matproz(const a:matrix; const x0:vector; var x1:vector; const n:byte);
var i,j:byte;
begin
for i:=1 to n do
x1[i]:=0;
for i:=1 to n do
for j:=1 to n do
x1[i]:=x1[i]+x0[j]*a[i,j];
end;
var x0,x1:vector;
a:matrix; l0,l1,e:real;
n,i,j:byte;
begin
writeln('Vvedite rang matrici');
write('n=');
readln(n);
for i:=1 to n do
x0[i]:=1/sqrt(n);
x0[1]:=1;
writeln('Vvedite matricu');
  for i := 1 to n do
    for j := 1 to n do
    begin
      write('a[', i, ',', j, ']=');
      readln(a[i, j]);
    end;
  matproz(a,x0,x1,n);
    l1:=sp(x0,x1,n);
  normir(x1,n);
repeat
  x0:=x1;
  l0:=l1;
  matproz(a,x0,x1,n);
  l1:=sp(x0,x1,n);
 normir(x1,n);
 
  until(abs(l1-l0)<e);
  writeln(x1[1]:5:2,' ',x1[2]:5:2,' ',x1[3]:5:2);
   writeln(l1);
end.

Добавлено через 11 часов 24 минуты
Все ответ найден. Вопрос аннулирован

Новые блоги и статьи Все статьи Все блоги /
Символьное дифференцирование igorrr37 13.02.2026 / * Логарифм записывается как: (x-2)log(x^2+2) - означает логарифм (x^2+2) по основанию (x-2). Унарный минус обозначается как ! */ #include <iostream> #include <stack> #include <cctype>. . .	Камера 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, то после закрытия окошка. . .