Matlab đại số tuyến tính câu thực hiện trong command window
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 (133.11 KB, 32 trang )
<span class="text_page_counter">Trang 1</span><div class="page_container" data-page="1">
<b>TRƯỜNG ĐẠI HỌC BÁCH KHOA TP.HCM</b>
<b>KHOA KHOA HỌC ỨNG DỤNG</b>
<b>MƠN HỌC:ĐẠI SỐ TUYẾN TÍNH</b>
<b>BÁO CÁO BÀI TẬP LỚN ĐẠI SỐ TUYẾN TÍNHTÊN ĐỀ TÀI:ĐỀ TÀI 1</b>
<i><b>Giảng viên hướng dẫn : Lê Thị Yến NhiLớp bài tập : L01</b></i>
<i><b>Nhóm :01</b></i>
<b> </b>
<b> </b>
<small>Báo cáo tháng 1/2017</small></div><span class="text_page_counter">Trang 5</span><div class="page_container" data-page="5">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
<b>3.Thuật toán</b>
</div><span class="text_page_counter">Trang 6</span><div class="page_container" data-page="6">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ì thố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
</div><span class="text_page_counter">Trang 7</span><div class="page_container" data-page="7">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
</div><span class="text_page_counter">Trang 8</span><div class="page_container" data-page="8"><small>function bdsc_mtbt</small>
<small> A=input('Nhap A: ');%nhap giá tri ma tran can tính tốn%</small>
<small> [m,n]=size(A);%gán cho m là so hàng n là so cot cua ma tran%</small>
<small> h=1;%gán cho h=1 lát nua se dùng nhý là bien chi cot%</small>
<small> for k=1:m-1;%cho k chay tu 1-> m-1%</small>
<small> 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</small>
<small> for i=k:m;%ung voi moi giá tri j ta lai cho i(là bien chi hàng)chay tu k=1 ->m</small>
<small> 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)</small>
<small> 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'%</small>
<small> 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%</small>
<small> tmp = m; m = n; n = tmp;% doi bien m voi n%</small>
<small> for s=1:n;%g?i s la bien so hang chay tu 1->n%</small>
<small> i=n+1-s; %bien i cung la bien so hang nhung ch?y nguoc lai bien s%</small>
<small> for j=1:m;%j la bien so cot%</small>
<small> 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%</small>
<small> k=i;%gan cho k tri so cua hang co cac phan tu khac 0%</small>
<small> break;%neu tim duoc mot gia tri cua ma tran A bang cach quet nhu vay thi dung lai</small>
<small> end;</small>
</div><span class="text_page_counter">Trang 12</span><div class="page_container" data-page="12">Chương trình hoạt động với mọi ma trận nhập
</div><span class="text_page_counter">Trang 20</span><div class="page_container" data-page="20">else Disp(' V1 khong vuong goc V2')
</div>
