Code Matlab phân tích A = LU

Dưới đây là Code Matlab phân tích A = LU. Mời bạn đọc đón xem!

74 37 lượt tải Tải xuống
Chia sẻ Báo cáo tài liệu

Môn: Kỹ thuật lập trình

Trường: Công nghệ thông tin

Thông tin:

2 trang 6 tháng trước

Tác giả:

Profile Mai Anh Đào
%Tạo A ở workspace để chương trình nhận dạng
function f(A)
%Input: ma tran vuong A
%Output: ma tran tam giac tren U va ma tran tam giac
duoi don vi L
A = [1 2 3;
3 4 5;
9 21 1];%Nhap ma tran A
[m,n]=size(A);
U=A;
L=zeros(n);
if U(1,1)==0 error('Nhap lai ma tran khac');
end;
%Tim U,L
for i=1:n
L(i,i)=1;
L(i,1)=A(i,1)/A(1,1);
end;
for k=2:n
for i=k:n
for j=k-1:n
U(i,j)=U(i,j)-U(k-1,j)*L(i,k-1);
end;
L(i,k)=U(i,k)/U(k,k);
end;
if U(k,k)==0
disp('U=');disp(U);
disp('i=');disp(k);
error('hoan doi hang i voi hang khac');
end;
end;
format rat;%Hien ket qua duoi dang so thap phan
disp('A=');disp(A);
disp('U=');disp(U);
disp('L=');disp(L);
%Kiem tra lai ket qua
disp('det(A)=');disp(det(A));
disp('det(U)*det(L)=');disp(det(U)*det(L));
| 1/2
| | Tải xuống

Có thể bạn quan tâm:

Preview text:

%Tạo A ở workspace để chương trình nhận dạng function f(A) %Input: ma tran vuong A
%Output: ma tran tam giac tren U va ma tran tam giac duoi don vi L A = [1 2 3; 3 4 5; 9 21 1];%Nhap ma tran A [m,n]=size(A); U=A; L=zeros(n);
if U(1,1)==0 error('Nhap lai ma tran khac'); end; %Tim U,L for i=1:n L(i,i)=1; L(i,1)=A(i,1)/A(1,1); end; for k=2:n for i=k:n for j=k-1:n
U(i,j)=U(i,j)-U(k-1,j)*L(i,k-1); end; L(i,k)=U(i,k)/U(k,k); end; if U(k,k)==0 disp('U=');disp(U); disp('i=');disp(k);
error('hoan doi hang i voi hang khac'); end; end;
format rat;%Hien ket qua duoi dang so thap phan disp('A=');disp(A); disp('U=');disp(U); disp('L=');disp(L); %Kiem tra lai ket qua disp('det(A)=');disp(det(A));
disp('det(U)*det(L)=');disp(det(U)*det(L));