Уравнение с вырожденной матрицей

Форум пользователей пакета Maple

Модератор: Admin

KateXX
Сообщения: 71
Зарегистрирован: Сб фев 20, 2010 12:00 pm

Уравнение с вырожденной матрицей

Сообщение KateXX » Чт июн 10, 2010 1:36 pm

Здравствуйте,
Методика вычисления нормального решения СЛАУ с вырожденной матрицей приведена (извиняюсь за качество изображения). Не понимаю сути [A'B], где A' - транспонированная матрица А...
Теоретически по методу должно быть так:
restart:
with(student):
with(linalg):
# Система уравнений с вырожденой матрицей
C := array([[2, -1, sqrt(2)], [3, 2, -3], [3, sqrt(2), -15/7]]):
D := array([5 + 7*sqrt(2), -24, -12 - 3*sqrt(2)]):
D := `<,>`(D):
print(C):
print(D):
# Параметр приближения
t := 1*10^(-6):
alpha := 0.5*sqrt(3*t):
E := array([[1, 0, 0], [0, 1, 0], [0, 0, 1]]):
# Сведение к приближенному решению
A := array(1..3, 1..3):
A := evalm((multiply(transpose(C),C)) + (evalm(alpha*E))):
B := multiply(transpose(C),D):
X := evalf(multiply(inverse(A), B)):
print(X):
Но что-то не то. И не знаю, что именно...
Изображение
Изображение