Phương pháp tính: chiến lược chọn trụ trong phương pháp khử gauss
Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (24.46 KB, 2 trang )
Chiến lược chọn trụ trong phương pháp khử
Gauss
Nguyễn Quốc Lân - BM Toán Ứng Dụng, DHBK TP. HCM
> with(linalg):
> Digits:=4;
Khai báo ma trận AA: ma trận mở rộng của hệ Ax = b
AA:=matrix(2,3,[0.003,59.14,59.17,5.291,-6.130,46.78]);
Khử cột 1: phần tử trụ khác 0, tính hệ số khử
m21:= AA[2,1]/AA[1,1];
Biến đổi sơ cấp trên hàng
AA:=addrow(AA,1,2,-m21);
Ý nghóa của giải thuật khử Gauss: Phần tử được khử phải triệt
tiêu
AA[2,1]:=0;
> evalm(AA);
Giải lùi
x:=backsub(AA);
> B:=matrix(2,3,[0.003,59.14,59.17,5.291,-6.130,46.78]);
Hoán vò hàng để tìm ra trụ tối đại
BB := swaprow(B,1,2);
> m21:= BB[2,1]/BB[1,1];
> BB:=addrow(BB,1,2,-m21);
> x:=backsub(BB);
Trường hợp thất bại của phương pháp trụ tối đại
C:=matrix(2,3,[30.00,591400.,591700.,5.291,-6.130,46.78]);
> m21:=C[2,1]/C[1,1];
> C:=addrow(C,1,2,-m21);
> C[2,1]:=0;
> evalm(C);
> x:=backsub(C);
>
