21.1 Создание и манипулирование диагональными/перестановочными матрицами

Диагональная матрица определяется как матрица, имеющая нулевые элементы вне главной диагонали; то есть, D(i,j) == 0 если i != j. Чаще всего рассматриваются квадратные диагональные матрицы; однако, это определение в равной степени может быть применено к неквадратным матрицам, в этом случае мы обычно говорим о прямоугольной диагональной матрице.

Перестановочная матрица определяется как квадратная матрица, имеющая единственный элемент, равный единице, в каждой строке и каждом столбце; все остальные элементы равны нулю. То есть, существует перестановка (вектор) p такая, что P(i,j) == 1 если j == p(i) и P(i,j) == 0 в противном случае.

Octave предоставляет специальное обращение с действительными и комплексными прямоугольными диагональными матрицами, а также с перестановочными матрицами. Они хранятся как специальные объекты, используя эффективные методы хранения и алгоритмы, что облегчает написание как удобочитаемых, так и эффективных выражений матричной алгебры на языке Octave. Специальное обращение можно отключить, используя функции disable_diagonal_matrix и disable_permutation_matrix.

val = disable_diagonal_matrix ()

old_val = disable_diagonal_matrix (new_val)

disable_diagonal_matrix (new_val, "local")

Запрос или установка значения, определяющего, отключено ли хранение диагональных матриц в специальном формате с экономией места.

Значение по умолчанию — false. Если это значение установлено в true, Octave будет хранить диагональные матрицы как полные матрицы.

При вызове изнутри функции с параметром "local" это значение изменяется локально для функции и всех вызываемых ею подпрограмм. Исходное значение восстанавливается при выходе из функции.

См. также: disable_range, disable_permutation_matrix.

val = disable_permutation_matrix ()

old_val = disable_permutation_matrix (new_val)

disable_permutation_matrix (new_val, "local")

Запрос или установка значения, определяющего, отключено ли хранение перестановочных матриц в специальном формате с экономией места.

Значение по умолчанию — false. Если это значение установлено в true, Octave будет хранить перестановочные матрицы как полные матрицы.

При вызове изнутри функции с параметром "local" это значение изменяется локально для функции и всех вызываемых ею подпрограмм. Исходное значение восстанавливается при выходе из функции.

См. также: disable_range, disable_diagonal_matrix.

Экономия места существенна, как показано в следующем коде.

x = diag (rand (10, 1));
xf = full (x);
sizeof (x)
⇒ 80
sizeof (xf)
⇒ 800