# Вторая пачка часть 186

ВЫЧИСЛЕННЫЕ СОБСТВЕННЫЕ ЗНА-ЧЕНИЯ РАЗМЕЩАЮТСЯ В ПОРЯДКЕ ВОЗРАСТАНИЯ В МАССИВЕ D НА МЕСТЕ ДИАГОНАЛЬНЫХ ЭЛЕМЕНТОВ. СОБСТВЕННЫЕ ВЕКТОРЫ ФОРМИРУЮТСЯ В МАССИВЕ Z [1:N, 1:N] НА МЕСТЕ МАТРИЦЫ ОРТОГОНАЛЬНЫХ ПРЕОБРАЗОВАНИЙ Z. ПРИ РАБОТЕ ПРОЦЕДУРЫ ПРОИСХОДИТ ВЫХОД ПО КОМАНДЕ EXIT, ЕСЛИ ДЛЯ ОПРЕДЕЛЕНИЯ СОБСТВЕННОГО ЗНАЧЕНИЯ ПО¬ТРЕ¬БО¬ВАЛОСЬ БОЛЕЕ 30 ИТЕРАЦИЙ ;) LABEL CONT, ROOT,NEXTIT; VAR I, J, K, L, M : INTEGER; B, C, F, G, H, P, R, S : REAL; BEGIN FOR I := 2 TO N DO E[I-1] := E[I]; FOR I := 1 TO N DO FOR J := 1 TO N DO IF I=J THEN ZA[I,J] := 1.0 ELSE ZA[I,J] := 0.0; E[N] := 0.0; B:= 0.0; F := 0.0; (; ОПРЕДЕЛЕНИЕ САМОГО МАЛЕНЬКОГО ЧИСЛА. ПОЛУЧЕННОЕ ЧИСЛО ЗАВИСИТ ОТ ТИПА ИСПОЛЬЗУЕМОЙ ЭВМ ;) FOR L := 1 TO N DO BEGIN J := 0; H := MASHEPS * (ABS(D[L])+ABS(E[L])); IF B<H THEN B:= H; (;ОПРЕДЕЛЕНИЕ ПРЕНЕБРЕЖИМО МАЛОГО ПОДДИАГОНАЛЬНОГО ЭЛЕМЕНТА. ЗАВИСИТ ОТ ТИПА ИСПОЛЬЗУЕМОЙ ЭВМ;) FOR M := L TO N DO IF ABS(E[M])<=B THEN GOTO CONT; CONT: IF M=L THEN GOTO ROOT; NEXTIT: IF J=30 THEN EXIT; INC(J); (; ФОРМИРОВАНИЕ СДВИГА ;) G := D[L]; P:= (D[L+1]-G) / (2*E[L]); R:= SQRT (P*P+1.0); IF P<0.0 THEN D[L] := E[L]/(P-R) ELSE D[L] := E[L]/(P+R); H := G - D[L]; FOR I := L+1 TO N DO D[I] := D[I] - H; F := F + H; (; QL - ПРЕОБРАЗОВАНИЕ ;) P:= D[M]; C := 1; S := 0; FOR I := M-1 DOWNTO L DO BEGIN G := C*E[I]; H := C*P; IF ABS(P)>=ABS(E[I]) THEN BEGIN C := E[I]/P; R := SQRT(C*C+1); E[I+1] := S*P*R; S := C / R; C := 1.0 / R; END ELSE BEGIN C:= P/E[I]; R := SQRT(C*C+1); E[I+1] := S*E[I]*R; S := 1.0/R; C := C / R; END; P := C*D[I] - S*G; D[I+1] := H + S*(C*G + S*D[I]); (; ФОРМИРОВАНИЕ СОБСТВЕННОГО ВЕКТОРА ;) FOR K := 1 TO N DO BEGIN H := ZA[K,I+1]; ZA[K,I+1] := S*ZA[K,I] + C*H; ZA[K,I] := C*ZA[K,I] - S*H; END; END; E[L] := S*P; D[L] := C*P; IF ABS(E[L])>B THEN GOTO NEXTIT; ROOT: D[L] := D[L] + F; END; (; РАЗМЕЩЕНИЕ СОБСТВЕННЫХ ЗНАЧЕНИЙ И СОБСТВЕННЫХ ВЕКТОРОВ;) FOR I := 1 TO N DO BEGIN K := I; P := D[I]; FOR J := I+1 TO N DO IF D[J] < P THEN BEGIN K:=J; P:=D[J]; END; IF K<>I THEN BEGIN D[K] := D[I]; D[I] := P; FOR J := 1 TO N DO BEGIN P := ZA[J,I]; ZA[J,I] := ZA[J,K]; ZA[J,K] := P; END; END; END; END; {$F+} PROCEDURE TRBAK (N,M1,M2:INTEGER;AZ:MAS1; E:MAS; VAR ZA:MAS1); (; ВЫЧИСЛЕНИЕ СОБСТВЕННЫХ ВЕКТОРОВ ИСХОДНОЙ СИМ¬МЕТ¬РИ¬ЧЕС¬КОЙ МАТРИЦЫ, КОГДА ИЗВЕСТНЫ СОБСТВЕННЫЕ ВЕКТОРЫ ТРЕХ¬ДИАГОНАЛЬНОЙ МАТРИЦЫ, РАЗМЕЩЕННЫЕ В МАССИВЕ Z[1:N,M1:M2]

Индекс
Элементарные функции Линейные уравнения Нелинейные уравнения Случайные числа