ЛИНЕЙНАЯ АЛГЕБРА 12 DET Определитель матрицы. Аналитическое нахождение определителя смотрите пункт 6а. Функция r = det (A) вычисляет определитель символьной и числовой матрицы с результатом в символьном или числовом виде соответственно.
Пример:
A = sym([9 1 3; -2 -2 1; 7 2 2])
A =
[ 9, 1, 3]
[ -2, -2, 1]
[ 7, 2, 2]
» r = det(A)
r =
-13
13 . POLY Характеристический полином матрицы. СПРАВКА
Алгебраические элементарные функции, определенные формулами, содержащими конечное число алгебраических или тригонометрических операций, производимых над аргументом, функцией, постоянными путем сложения, вычитания, умножения, деления, возведением в степень, извлечением корня, логарифмированием, взятием прямых или обратных тригонометрических функций, называются целыми функциями или полиномами.
Функция p = poly (A) вычисляет коэффициенты характеристического полинома матрицы А. Если матрица представлена в символьном виде, то и полином выдается в символьном виде, используя по умолчанию переменную х.
Функция p = poly(A,v) позволяет указать независимую переменную.
Пример: A = gallery(3) ,где gallery – тестовая матрица в программе MATLAB
A =
-149 -50 -154
537 180 546
-27 -9 -25
» p = poly(A)
p =
1 -6 11 -6
» g = poly(sym(A))
g =
x^3-6*x^2+11*x-6
» s = poly(sym(A),sym('z'))
s =
z^3-6*z^2+11*z-6
14 .DIAG Формирование или извлечение диагоналей матрицы. СПРАВКА
Вычисление определителя упрощается, если в матрице элементы, расположенные по одну сторону от главной диагонали, равны нулю. Такой определитель называется определителем треугольного вида. Из свойства матриц имеем: среди всех матриц размера n n существует единственная, что ее произведение на произвольную квадратную матрицу А слева или справа не изменит матрицу А АЕ = ЕА. Непосредственным умножением матрицы А на матрицу Е можно проверить, что матрицей Е, удовлетворяющей указанному равенству, является матрица, у которой элементы, расположенные на главной диагонали, равны единице, а остальные равны нулю (единичная матрица). Для использования этих свойств необходима проверка элементов матрицы по диагоналям.
Функция X =diag(v) формирует квадратную матрицу Х с вектором v на главной диагонали.
Функция Х = diag(v,k) формирует квадратную матрицу Х порядка length(v)+abs(k), где length – длина вектора, abs – значение комплексного числа, c вектором v на k – й диагонали.
Функция v = diag(X) извлекает из матрицы Х главную диагональ.
Функция v = diag(X,k) извлекает из матрицы Х диагональ с номером k; при k 0, это номер k-й верхней диагонали, при k 0 это номер k-й нижней диагонали.
Пример:
syms a b c
» v = [a b c]
v =
[ a, b, c]
» diag(v)
ans =
[ a, 0, 0]
[ 0, b, 0]
[ 0, 0, c]
» diag(v,-2)
ans =
[ 0, 0, 0, 0, 0]
[ 0, 0, 0, 0, 0]
[ a, 0, 0, 0, 0]
[ 0, b, 0, 0, 0]
[ 0, 0, c, 0, 0]
» syms x y z
» A = [a,b,c; 1,2,3; x,y,z]
A =
[ a, b, c]
[ 1, 2, 3]
[ x, y, z]
» diag(A)
ans =
[ a]
[ 2]
[ z] syms a b c
» A = [a,b,c; 1,2,3; a+1,b+2, c+3]
A =
[ a, b, c]
[ 1, 2, 3]
[ a+1, b+2, c+3]
» tril(A)
ans =
[ a, 0, 0]
[ 1, 2, 0]
[ a+1, b+2, c+3]
» tril(A,1)
ans =
[ a, b, 0]
[ 1, 2, 3]
[ a+1, b+2, c+3]
» tril(A,-1)
ans =
[ 0, 0, 0]
[ 1, 0, 0]
[ a+1, b+2, 0]
» diag(A,1)
ans =
[ b]
[ 3]
|