Линейная алгебра. Теоремы и алгоритмы. Яцкин Н.И. - 541 стр.

UptoLike

ВУЗ: 

ИвГУ | Иваново

Составители: 

Рубрика: 

Прил. 1 Коды Maple-процедур 541
if ms<n then
exist_jbas:=false;
print(`Не существует полного жорданова базиса,
только - частичный.`);
else
print(`Полный жорданов базис существует.`);
fi;
print(str);
# Заготовки для рабочих массивов B,BGJ,d,Nbas,F,l.
# Элементы всех этих массивов (кроме последнего)
# сами являются индексированными переменными.
# Массив степеней матриц B[i]=A-lambda[i]*E:
# B[i][k]=(A-lambda[i]*E)^k,
# где i - номер собственного значения (i=1,...,s),
# k - номер итерации (0<=k<=m[i]).
B:=array(1..s);
# Те же матрицы, приведенные к виду Жордана-Гаусса.
BGJ:=array(1..s);
# Массив итерированных дефектов: d[i][k]=dfc(B[i][k]);
d:=array(1..s);
# Массив базисов в итерированных ядрах:
# множество Nbas[i][k] должно содержать
# необработанный базис
# в итерированном ядре N[i][k]=Ker(B[i][k]).
Nbas:=array(1..s);
# Массив фундаментальных матриц:
# матрица F[i][k] формируется
# из векторов множества Nbas[i][k].
F:=array(1..s);
# Массив показателей стабилизации:
# l[i] - показатель стабилизации для матрицы B[i].
l:=array(1..s);
# Заполнение массивов B,BGJ,d,Nbas,F,l.
for i from 1 to s do
# Начальные значения для
# индексированных переменных - элементов массивов B,d.
B[i][0]:=E:d[i][0]:=0:
# Начало вычисления показателей стабилизации.
l[i]:=1: 
Прил. 1                 Коды Maple-процедур                541

if ms

Страницы