И другие программы этой серии
Пример: для многочлена
x^3+x^2-6x-8 вектор полиномиальных коэффициентов r имеет следующий
вид:
>> r=[1,1,-6,-8]
r=
1
1 -6 -8
>> A=compan(r) % nii?iai?aa?uay iao?eoa
A=
-1
6
8
1
0
0
0
1
0
>> eig(compan(r)) % ei?ie iiiai?eaia
ans =
-2.0000
2.5616
-1.5616
• triu(X) – возвращает неизменной верхнюю треугольную часть матрицы O,
включая элементы главной диагонали, и заменяет нулями остальные элементы;
• triu(X,k) – возвращает неизменной верхнюю треугольную часть матри
цы O, начиная с k ой диагонали. При k=0 это главная диагональ, при k>0 –
одна из верхних диагоналей, при k<0 – одна из нижних диагоналей.
Пример:
M=
3
8
8
>> triu(M)
ans =
3
0
0
1
3
1
4
2
1
1
3
0
4
2
1
4.2.9. Операции с пустыми матрицами
Возможно создание пустых матриц, например:
>> M=[]
M=
[]
>> whos M
Name
Size
Bytes
Class
M
0x0
0
double array
Grand total is 0 elements using 0 bytes
Из этого примера видно, что пустая матрица создается с помощью пустых
квадратных скобок. С помощью функций zeros, ones, rand или eye также можно
создавать пустые матрицы заданного размера. Например:
>> E=zeros(0,5)
E = Empty matrix: 0-by-5
4.3.2. Тестовые матрицы
Для выполнения ряда вычислений в области линейной алгебры создан ряд специ
альных матриц, именуемых тестовыми матрицами. Такие матрицы создаются
указанными ниже функциями.
[A,B,C,…] = gallery(\'tmfun\',P1,P2,…) – возвращает тестовые матри
цы, определенные строкой tmfun, где tmfun – это имя семейства матриц, вы
С пустыми матрицами можно выполнять некоторые операции, например:
>> sum(M)
ans = 0
208
Операции с векторами и матрицами
Создание и вычисление специальных матриц
209
бранное из списка. P1, P2, … – входные параметры, требуемые для конкретного
семейства матриц. Число используемых параметров P1, P2, … изменяется от мат
рицы к матрице. Функция gallery хранит более 50 различных тестовых матрич
ных функций, полезных для тестирования численных алгоритмов и других целей
(включая матрицы Коши, Чебышева, фон Неймана, Чебышева–Вандермонде,
Вандермонде, Уилкинсона и т.