(******************************************************************************)
(**)                        ОТДЕЛ ПроLL+;
(******************************************************************************
 * НАЗНАЧЕНИЕ: 
 *   разложение симметрической положительно определенной матрицы A на две
 *   треугольные матрицы L и L', где L - это нижняя треугольная матрица,
 *   а L' - матрица, транспонированная к L
 *
 * ПРИМЕЧАНИЯ: 
 *     Исходная матрица записана в файле "ПроМатр4.дан".
 *   В этом файле 1-я строка содержит размер матрицы (число строк и столбцов),
 *   а последующие строки содержат значения элементов матрицы.
 ******************************************************************************)
ИСПОЛЬЗУЕТ
  Матем  ИЗ "...\Отделы\Числа\",
  Матр   ИЗ "...\Отделы\Числа\",
  МатрВВ ИЗ "...\Отделы\Числа\",
  Читать ИЗ "...\Отделы\Обмен\",
  Вывод  ИЗ "...\Отделы\Обмен\";

(******************************************************************************)
ЗАДАЧА Перемножить(L-,A+:Матр.Вид); (* A=L*L' *)
ПЕР
  i,j,k,n:ЦЕЛ;
  m:Матем.Вещ;
УКАЗ
  n:=РАЗМЕР(L);
  ОТ i:=0 ДО n-1 ВЫП
    ОТ j:=0 ДО i ВЫП
      m:=0;
      ОТ k:=0 ДО j ВЫП
        m:=m+L[i,k]*L[j,k]
      КОН;
      A[i,j]:=m;
      ЕСЛИ i # j ТО
        A[j,i]:=m 
      КОН
    КОН 
  КОН
КОН Перемножить;

(******************************************************************************)
ЗАДАЧА Голова;
ПЕР
  n:ЦЕЛ;         (* размер матрицы *)
  A:Матр.Доступ; (* исходная матрица *)
  L:Матр.Доступ; (* матрица L *)
  п:Читать.Поток;
УКАЗ
  п:=МатрВВ.ОткрытьПоток("ПроМатр4.дан");
  A:=МатрВВ.ЧитатьМатр(п);
  Читать.Закрыть(п);
  n:=РАЗМЕР(A^);
  СОЗДАТЬ(L,n,n);
  Вывод.Цепь("^Исходная матрица:^");
  МатрВВ.ВыводМатр("%10.5f",A^);
  ЕСЛИ Матр.РазложитьНаLL(A^,L^) = 0 ТО
    Вывод.Цепь("^L матрица:^");
    МатрВВ.ВыводМатр("%10.5f",L^);
    Перемножить(L^,A^);
    Вывод.Цепь("^L*L' матрица:^");
    МатрВВ.ВыводМатр("%10.5f",A^)
  ИНАЧЕ
    Вывод.Цепь("^Матрица не является положительно определённой!^")
  КОН
КОН Голова;

(******************************************************************************)
УКАЗ
  Голова
КОН ПроLL.

▲ Вопросы, замечания и предложения высылайте на atimopheyev@yahoo.com или Издателю Глагола: