Matlab đại số tuyến tính câu thực hiện trong command window

TRƯỜNG ĐẠI HỌC BÁCH KHOA TP.HCMKHOA KHOA HỌC ỨNG DỤNGMÔN HỌC:ĐẠI SỐ TUYẾN TÍNHBÁO CÁO BÀI TẬP LỚN ĐẠI SỐ TUYẾN TÍNHTÊN ĐỀ TÀI:ĐỀ TÀI 1Giảng viên hướng dẫn : Lê Thị Yến NhiLớp bài tập : L

TRƯỜNG ĐẠI HỌC BÁCH KHOA TP.HCM

KHOA KHOA HỌC ỨNG DỤNG

MÔN HỌC:ĐẠI SỐ TUYẾN TÍNH

BÁO CÁO BÀI TẬP LỚN ĐẠI SỐ TUYẾN TÍNHTÊN ĐỀ TÀI:ĐỀ TÀI 1

Giảng viên hướng dẫn : Lê Thị Yến NhiLớp bài tập : L01

Nhóm :01

Báo cáo tháng 1/2017

Phần tử cơ sở của hàng:là phần tử khác 0 đầu tiên của hàng đó kể từ bên trái sang

Ma trận bậc thang:là ma trận mà phần tử cơ sở của hàng dưới nằm bên phải phần tử cơ sở hàng trên -Hàng toàn số 0 (nếu có) thì nằm dưới.

Hạng của ma trận:

-Mọi ma trận đều có thể đưa về dạng bậc thang bằng các phép biến đổi sơ cấp (theo hàng hoặc cột) -Hạng của ma trận là số hàng khác 0 của ma trận bậc thang

3.Thuật toán

B6 quét ma trận lần lượt từng hàng từ trên xuống nếu gặp phần tử của ma trận khác 0 thì thoát khỏi vòng lặp

B7 nếu i~=k thì ta chuyển chỗ các phần tử của ma trận nằm trên hàng đó với hàng kề với nó

B8:cho i chạy k+1->m

B9:nếu phần tử A(i,h) khác 0 thì ta gán a là thương của phép chia A(i,h)/A(k,h)

B10: cho j chạy từ h->n ,ta lấy lần lượt các giá trị của ma trận trong một hàng trừ đi các phần tử phía trên nó A(i,j)=A(i,j)-a*A(k,j)

B11: cho h=h+1 và lặp lại bước 3B12 hiển thị ma trận bậc thang

B13 đặt i=n+1-s với s chạy từ 1->n để cho i chạy ngược lại so với s

B14: cho j chạy từ 1->m

B15 :quét từ dưới lên từng hàng một từ trái qua phải nếu gặp một phần tử ma trận khác không nào thì hạng là số hàng mà phần tử đó thuộc

B16: hiển thị hạng

function bdsc_mtbt

A=input('Nhap A: ');%nhap giá tri ma tran can tính toán%

[m,n]=size(A);%gán cho m là so hàng n là so cot cua ma tran%

h=1;%gán cho h=1 lát nua se dùng nhý là bien chi cot%

for k=1:m-1;%cho k chay tu 1-> m-1%

for j=h:n;%ung voi moi giá tri k ðang chay thi ta cho j(là bien chi cot) chay tu h=1 ->n

for i=k:m;%ung voi moi giá tri j ta lai cho i(là bien chi hàng)chay tu k=1 ->m

if(A(i,j)~=0)%i và j chay cho ðen khi phan tu ma tran thuoc hàng i cot j khác 0(hieu là:neu quét ma tran tung cot mot tu trên xuong dýoi cho ðen khi gap mot phan tu khác 0 thi thôi)

break;%neu phan tu hàng i cot j khác 0 thi nó se dung và bo lenh 'for'ðang chua nó tuc là lai bat ðau chay tu 'for j=h:n'%

A(i,j)=A(i,j)-a*A(k,j);%sau do phan tu ma tran hang i cot j se nhan duoc gia tri moi bang phep tru tren%

tmp = m; m = n; n = tmp;% doi bien m voi n%

for s=1:n;%g?i s la bien so hang chay tu 1->n%

i=n+1-s; %bien i cung la bien so hang nhung ch?y nguoc lai bien s%

for j=1:m;%j la bien so cot%

if A(i,j)~=0%nghia la lan luot quet ma tran lan luot tu duoi len tren tung hang tu trai sang phai , khi tim duoc phan tu khac 0 thi dung lai%

k=i;%gan cho k tri so cua hang co cac phan tu khac 0%

break;%neu tim duoc mot gia tri cua ma tran A bang cach quet nhu vay thi dung lai

end;

Chương trình hoạt động với mọi ma trận nhập

else Disp(' V1 khong vuong goc V2')