ĐI HC QUC GIA TP.HCM
TRƯNG ĐI HC BCH KHOA
BO CO BI TP LN
MATLAB
MÔN: ĐI S TUYN TNH
GVHD: Phan Th Khnh Vân
Đ ti 3:
Cho mt h vc tơ M. Kim tra M c đc lp tuyn tnh hay không? Nu c hy
trc chun h vc tơ M.
Không đưc dng lnh “qr”.
STT
H V TÊN
MSSV
LP
1
Đon Lê Mnh Cưng
1510353
DC1504
2
Hong Văn Long
1511804
DC1504
3
H Văn Thng
1513121
DC1504
4
Lê Văn Đnh
1510755
DC1504
5
Nguyn Ngô Long
1511814
DC1504
6
Nguyn Thi Quang
1413098
DC1504
7
Nguyn Văn Thnh
1513055
DC1504
8
Nguyn Văn Trng
1513704
DC1504
9
Nguyn Xuân Trc
1513804
DC1504
TP.HCM, Ngy 1, Thng 6, 2016
1 | P a g e
BO CO MATLAB ĐI S
Mc lc
Tiêu đ Trang
I. Li m đu ……………………………………………………………..2
II. Đ ti v phân tch……………………………………………………...3
III. Bo co đ ti …………………………………………………………..4
IV. Nhn xt ca ging viên hưng dn ……………………………..........6
2 | P a g e
BO CO MATLAB ĐI S
I. LI M ĐU…
Ngy nay khoa hc ngy cng pht trin, vi đ pht trin ny vic ng dng
khoa hc v sng ch khoa hc trưng hc l rt thit thc v quan trng. Chnh
v vy, ngay t năm đu cc ging viên trưng ĐH Bch Khoa TP.HCM đ gip cho
cc sinh viên ngnh k thut lm quen vi cc ng dng lp trnh, v d như Chương
trnh Matlab.
MATLAB l mt môi trưng tnh ton s v lp trnh cho php tnh ton s
vi ma trn, v đ th hm s hay biu din thông tin, thc hin thut ton, to cc
giao din ngưi dng v liên kt vi nhiu my tnh vit trên nhiu ngôn ng lp
trnh khc. Vi thư vin Toobox, MATLAB cho php mô phng tnh ton, thc
nghim nhiu mô hnh trong thc t v k thut. Vi hơn 40 năm hnh thnh v pht
trin, ngy nay vi thit k s dng tương đi đơn gin v ph thông, MATLAB l
công c tnh ton hu hiu đ gii quyt cc bi ton k thut.
V vy, đi vi nhng bi ton trong môn Đi s, đc bit l nhng bi ton
Ma trn, h phương trnh, nh x tuyn tnh,… MATLAB c th gip ta gii quyt
nhng bi ton đ mt cch đơn gin v hiu qu, phn mm tr thnh mt tr th
đc lc cho c ging viên v sinh viên trong ging dy v hc tp.
3 | P a g e
BO CO MATLAB ĐI S
II. Đ TI V PHÂN TCH
Đ ti 3:
Cho mt h vc tơ M. Kim tra M c đc lp tuyn tnh hay không? Nu c hy
trc chun h vc tơ M.
Không đưc dng lnh “qr”.
Phân tch:
Input: Cho php nhp vo mt h vc tơ M.
Output: Nu M đc lp tuyn tnh th trc chun h vc tơ M
Cơ s l thuyt :
Kim tra cc vc tơ x1, x2,…, xm c đc lp tuyn tnh hay không?
Đt A = ( x𝑇
1 x𝑇
2 … x𝑇
𝑚 ) v xc đnh r(A)
Nu r(A) = m th x1, x2,…, xm đc lp tuyn tnh.
Hai vc tơ x1 , x2 trc chun nu {<𝑥
1 ,𝑥
2> =0
‖𝑥
1‖= ‖𝑥
2‖=1
Thut ton:
Nhp h vc tơ M t bàn phím.
Kim tra M c đc lp tuyn tính.
Trc chun h vector M (nu M đc lp tuyn tính)
4 | P a g e
BO CO MATLAB ĐI S
III. BO CO Đ TI
Đ ti 3:
Cho mt h vc tơ M. Kim tra M c đc lp tuyn tnh hay không? Nu c hy
trc chun h vc tơ M.
Không đưc dng lnh “qr”.
Đon CODE:
function tracchuan
f=input('nhap ho vecto M: ');
[m,~]=size(f);
if rank(f)~=m
disp('M phu thuoc tuyen tinh')
return
end
fc=zeros(size(f));
fc(1,:)=f(1,:);
if m==1
fc=fc/(sqrt(f*f'));
else
alpha=zeros(m-1,1);
for i=1:m-1
fc(i+1,:)=f(i+1,:);
for j=1:i
alpha(j,1)=(f(i+1,:)*fc(j,:)')/(fc(j,:)*fc(j,:)');
fc(i+1,:)=fc(i+1,:)-alpha(j,1)*fc(j,:);
end
end
for i=1:m
fc(i,:)=fc(i,:)/(sqrt(fc(i,:)*fc(i,:)'));
end
end
disp('M doc lap tuyen tinh,truc chuan M:')
disp(fc)
end