Scilab Website | Contribute with GitLab | Mailing list archives | ATOMS toolboxes
Scilab Online Help
2024.1.0 - Русский


bdiag

блоковая диагонализация, обобщённые собственные векторы

Синтаксис

[Ab, X, bs] = bdiag(A)
[Ab, X, bs] = bdiag(A, rmax)

Аргументы

A

вещественная или комплексная квадратная матрица

rmax

вещественное число

Ab

вещественная или комплексная квадратная матрица

X

вещественная или комплексная невырожденная матрица

bs

вектор целых чисел

Описание

[Ab, X, bs] = bdiag(A, rmax) выполняет блоковую диагонализацию матрицы A. bs даёт структуру блоков (соответствующие размеры блоков). X - это изменение базиса, то есть Ab = inv(X)*A*X является блоковой диагональю.

rmax управляет обусловленностью X. Значение по умолчанию - это L1-норма матрицы A.

Чтобы получить диагональную форму (если она существует), выберите большое значение дляrmax (например rmax=1/%eps). В общем случае (для вещественной случайной A), блоки равны (1x1) и (2x2), а X является матрицей собственных векторов.

Примеры

// случай вещественных чисел: блоки 1x1 и 2x2
a = rand(5,5);
[ab, x, bs] = bdiag(a);ab
--> [ab, x, bs] = bdiag(a);ab
 ab  =
   2.4777836   0.          0.          0.         0.
   0.          0.0757087   0.4869251   0.         0.
   0.         -0.5777956  -0.1248605   0.         0.
   0.          0.          0.          0.069654   0.
   0.          0.          0.          0.         0.5341598
// случай комплексных чисел: комплексные блоки 1x1
[ab, x, bs] = bdiag(a+%i*0); ab

Смотрите также

  • schur — [ordered] Schur decomposition of matrix and pencils
  • sylv — Solve the Sylvester equations A*X + X*B = C or A*X*B + X = C
  • spec — собственные значения и собственные вектора матрицы или пучка
  • blockdiag — Создаёт блочную диагональную матрицу из указанных массивов. Соединение диагональной системы блоков.
  • lcmdiag — least common multiple diagonal factorization
Report an issue
<< balanc eigen eigs >>

Copyright (c) 2022-2024 (Dassault Systèmes)
Copyright (c) 2017-2022 (ESI Group)
Copyright (c) 2011-2017 (Scilab Enterprises)
Copyright (c) 1989-2012 (INRIA)
Copyright (c) 1989-2007 (ENPC)
with contributors
Last updated:
Mon Jun 17 17:55:05 CEST 2024