Автор: zhenya_nrok
Дата сообщения: 06.11.2007 07:12
Добрый день!
Имеется затруднение с разложением матрицы путем LU-факторизации.
Нужно написать программу и проверить путем встроенных функций.
Попытался, не совсем удачно:
A=[13.754 1.71 5.42; -6.39 14.25 1.954; 4.21 7.92 13.41]
b=[8.21; 5.42; 12.29]
%Проверка условий для использования метода LU-факторизации
opr=A(1,1)*A(2,2)*A(3,3)+A(1,2)*A(2,3)*A(3,1)+A(2,1)*A(3,2)*A(1,3)-(A(1,3)*A(2,2)*A(3,1)+A(1,2)*A(2,1)*A(3,3)+A(3,2)*A(2,3)*A(1,1))
d=det(A)
if d==0 'матрица вырожденная'
elseif isequal(d,opr) 'определитель определен правильно'
else 'определитель определен не правильно'
end
for j=1:3
U(1,j)=A(1,j);
end
U
for j=1:3
L(j,1)=A(j,1)/U(1,1);
end
L
for i=2
for j=2:3
for k=1
U(i,j)=A(i,j)-L(i,k)*U(k,j);
end
end
end
U
for k=1
for i=3
for j=2
L(i,j)=1/U(j,j)*(A(i,j)-L(i,k)*U(k,j));
end
end
end
L
for i=3
for j=3
U(i,j)=A(3,3)-(L(3,1)*U(1,3)+L(3,2)*U(2,3));
end
end
U
Подскажите как сделать? ошибки? может у кого есть готовая?
Заранее благодарен.
Прошу извенение за топ не совсем соответствующий тематике.