Đang tải... (xem toàn văn)
Báo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tínhBáo cáo bài tập lớn Đại số tuyến tính
I H C QU C GIA THÀNH PH H CHÍ MINH I H C BÁCH KHOA BÁO CÁO BÀI T P L N IS TUY N TÍNH L P CK16CK04 NHÓM ng d n: Nguy n Xuân M c: 2016-2017 DANH SÁCH THÀNH VIÊN Ph t Ph m H u B o H Tr Phan Nguy n Nh t Huy Nguy n Tr Nguy ng c Ki m 10 Lê Minh Hi u 11 MSSV 1612625 1614091 1610194 1612199 1611317 1610509 1611948 1613846 1611691 1611015 1611675 M CL C Ph n S PH C TRONG MATLAB L nh real, imag 1.1 1.2 Cú pháp 1.3 Ví d L nh abs 2.2 Cú pháp 2.3 Ví d L nh angle 3.2 Cú pháp 3.3 Ví d L nh conj 4.1 Ý 4.2 Cú pháp 4.3 Ví d Ph n MA TR N TRONG MATLAB L nh numel(A) L nh size Ghép ma tr n theo c t L nh inv(A) Tìm Ak Câu l nh [A B] Câu l nh A(:,n)=[ ] Câu l nh A(:,n:end) Câu l nh A(n,:)=[ ] 10.L nh zeros(n) 11.L nh eye(n) 12.L nh ones(n) 13.L nh rank(A) 14.L nh trace(A) 15.L 16.L nh det(A) 17.L nh tril(T) 18.L nh triu(T) 19.L nh reshape(A,m,n) 20.L nh A\b 21.L nh [Q,R]=qr(Y) ho c [L,U]=lu(Y) 22.L nh A[] 23.L nh A(i,j) 24.L nh A(i,:), A(:,j) 25.L nh A(i:k,:), A(:,j:k) 26.L nh rref(A) 27.L nh fliplr 28.L nh flipud 29.L nh magic 30.L nh pascal 31.L nh rand 32.L nh rot(90) 33.L nh isemty 34.L nh diag Ph n M T S L NH L N KHÔNG GIAN VECTOR, KHÔNG GIAN EUCLIDE, TR GIÊNG L L L L L nh dot nh cross nh length nh norm nh qr L nh [P,D]=eig(A) L nh eig(H) L nh max(X), min(X) Ph N XÉT C A GIÁO VIÊN Ph n S PH C TRONG MATLAB ph c >> i2 ans=-1 L nh real, imag 1.1 Real: l y ph n th c c a s ph c Imag: l y ph n o c a s ph c 1.2 Cú pháp phanthuc= real(z) phanao= imag(z) 1.3 Ví d >>z=5+6i >>phanthuc=real(z) phanthuc= >>phanao=imag(z) phanao= L nh abs 2.1 a s ph c 2.2 Cú pháp: y=abs(z) 2.3 Ví d >>z=3+4i z= 3.000 + 4.000i >> Modul=abs(z) Modul= L nh angle 3.1 3.2 Cú pháp: y=angle(z) 3.3 Ví d >> z= 3+4i z= 3.0000 + 4.0000i >> agumen=angle(z) a s ph c v agumen = 0.9273 L nh conj y s ph c liên h p c a s ph c 4.2 Cú pháp: y= conj(z) 4.3 Ví d >> z=3+4iz = 3.0000 + 4.0000i >> conj(z) radian ans = 3.0000 - 4.0000i Ph n MA TR N TRONG MATLAB L nh numel(A): m s ph n t c a a Ví d » A = [01 09 77; 20 04 2001 ] A= 77 20 2001 »u=numel(A) u=6 L nh size: Cho bi t s dòng c t c a m t ma tr n Ví d >> A= [1 3;4 5;2 6] A= fx >> size (A) ans = fx >> size (A,1) ans = fx >> size (A,2) ans = Ghép ma tr n theo c t: a, L nh: C=[A;B] V i: - A,B ma tr c C ma tr n c n tìm b, Ví d : Cho ma tr n A= [1 2], B=[ 7] , Ghép ma tr n A,B theo c t >> C =[1 2;3 4] C= 4 L nh inv(A): Tìm ma tr n ngh Ví d : >> A=[1 2;2 5] A= 2 >> inv(A) ans = -2 -2 o c a ma tr n L nh Ak: V i: - A ma tr s n tính Ví d : Bài 1: Cho ma tr n A=[1 2;2 4] >>A = 2 >> B=A3 B= 25 50 50 100 Câu l nh [A B] :Ghép ma tr n theo hàng -Cú pháp:[A B] -Ví d : >> A=[1 4;5 5;3 1;6 1] A= >> B=[3 5;6 2;3 3;5 9] B= 6 5 >> [A B] ans = 4 5 6 8 Câu l nh A(:,n)=[ ] : Xóa c t th n c a ma tr n A >> B=[3 5;6 2;3 3;5 9] B= 6 5 >> B(:,1)=[] B= Câu l nh A(:,n:end) : Cho phép l y t c t th -Cú pháp:A(:,n:end) n c t cu i c a ma tr n A= >> A(2,3) ans = 24.L nh A(i,:) A(:,j) : tham chi u dòng i tham chi u c t j Ví d : Cho A=[ 5; 6; 3] A= 5 8 >> A(2,:) ans = >> A(:,3) ans = 25 L nh A(i :k, :) A( :,j :k) : Tham chi u t dòng i d n dòng k Tham chi u t c nc tk Ví d A=[ 2; 8; 3] A= 5 >> A(1:2,:) ans = 5 >> A(:,1:2) ans = 5 26 L nh rref(A) : T o ma tr n b c thang t A Ví d A=[ 2; 1; 5] A= 2 >> rref(A) ans = 0 0 27.L nh FLIPLR : Chuy n ph n t c a ma tr n theo th t c - Cú pháp: b = fliplr(a) - Gi i thích: b: tên ma tr c chuy a: tên ma tr n c n chuy Ví d : i i c l i a= 9 b = fliplr(a) 28 L nh FLIPUD : Chuy n ph n t c a ma tr n theo th t - Cú pháp: b = flipud(a) - Gi i thích: b: tên ma tr c chuy a: tên ma tr n c n chuy - Ví d : >>a= [ 4; 5; 6] a= >>b = flipud(a) b= i i c l i 29 L nh MAGIC : T o ma tr n vng có t ng c a ph n t hàng, c t ho ng chéo b ng - Cú pháp: Tên ma tr n = magic(n) - Gi i thích: c ma tr n Giá tr c a m i ph n t ma tr n m t dãy s nguyên liên t c t 2n T ng hàng, c u b ng Ví d : >>tmt = magic(3) tmt = 30 L nh PASCAL :T o ma tr n theo quy lu n tam giác Pascal - Cú pháp:pascal (n) - Gi i thích:n: s hàng (c t) Ví d : pascal(4) ans = 1 1 n 3 10 10 20 31 L nh RAND : T o ma tr n mà k t mà giá tr c a ph n t ng u nhiên - Cú pháp: y = rand(n) y = rand(m,n) Gi i thích: - y: tên ma tr n -n: t o ma tr n có n hàng, n c t -m, n: t o ma tr n có m hàng, n c t - Giá tr c a ph n t n m kho ng [0 1] Ví d : >>y = rand(3) y= 0.9340 0.0920 0.7012 0.8462 0.6539 0.7622 0.5269 0.4160 0.7622 >> y = rand(3,5) y= 0.2625 0.3282 0.9910 0.9826 0.6515 0.0475 0.6326 0.3653 0.7227 0.0727 0.7361 0.7564 0.2470 0.7534 0.6316 32.L nh ROT90 : Xoay ma tr n 900 - Cú pháp: b = rot90(a) - Gi i thích: b: ma tr c xoay 900 a: ma tr n c n xoay Ví d : >>a=[1 3; 6; 9] a= >> b = rot90(a) b= 33.L nh isempty : Ki m tra xem ma tr n có ma tr n r ng không - Cú pháp : isempty(A) - Gi i thích : A ma tr =1 ma tr - c,n i ma tr n r ng, n u ans n r ng Ví d >> B= zeros(4) B= 0 0 0 0 0 0 0 0 >>isempty(B) ans = % B không ma tr n r ng >> A= [] A= [] >>isempty(A) ans = % A ma tr n r ng 34.L nh DIAG : T o ma tr n m i x - Cú pháp: v = diag(x) v = diag(x,k) c > Gi i thích: - x: vector có n ph n t - v: ma tr xn - k: tham s -N c t o t x theo quy t c: s hàng b ng s c t ph n t c a ng chéo c a v nh d ng cho v, s hàng c t c a v = n + abs(k) ng chéo c a v ph n t c a x - N u k > ph n t c a x n ng chéo v - N u k < ph n t c a x n ng chéo v Ví d : >>x = [ 4]; v = diag(x) v= 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 >>v1 = diag(x,2) v1 = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 >>v2 = diag(x,0) v2 = >>v3 = diag(x,-2) v3 = Ph n M T S L NH L N KHÔNG GIAN VECTOR, KHÔNG GIAN EUCLIDE, TR GIÊNG L nh dot -Cú pháp :dot(u,v) - c Ví d : >> u=[1 3] u= >> v=[3 5] v= >> dot(u,v) ans = 26 L nh cross(u,v) : Tích h Ví d : u= [ 3] u= >> v=[ 8] v= >> cross(u,v) ans = ng c a u, v 25 -10 L nh length : Tính chi u dài c Ví d x = [0 9] l = length (x) l = 10 » x = [01 09 77; 20 04 2001 ] x= 77 20 2001 » l=length(x) l= 4.L nh norm -Cú pháp: norm(u) ; u m t kì) t kì Ví d : u=[1 3] u= >> norm(u) ans = 3.7417 5.L nh qr:(tr c chu n hóa h t A) [P,]=qr(A) ; A ma tr n c ct ob ih Ví d c chu n véc >> A=[1 1] A= 1 >> qr(A) ans = 1.7321 -1.3660 -1.3660 6.L nh [P,D]=eig(A) : Chéo hóa Ví d A=[0 -8 6;-1 -8 7;1 -14 11] A= -8 -1 -8 -14 11 >> [P,D]=eig(A) P= -0.5774 -0.3244 0.2673 -0.5774 -0.4867 0.5345 -0.5774 -0.8111 0.8018 D= -2.0000 0 3.0000 0 2.0000 L nh eig(H) : Tr riêng Ví d H=[1 2;3 4] H= >> eig(H) ans = -0.3723 5.3723 L nh max(X), min(X) : Tr v giá tr l n nh t nh nh t vector X Ví d >> X=[ 9; 3; 4] X= 4 >> max(X) ans = >> min(X) ans = Ph n N XÉT C A GIÁO VIÊN ... Ví d : >> A=[1 2;2 5] A= 2 >> inv(A) ans = -2 -2 o c a ma tr n L nh Ak: V i: - A ma tr s n tính Ví d : Bài 1: Cho ma tr n A=[1 2;2 4] >>A = 2 >> B=A3 B= 25 50 50 100 Câu l nh [A B] :Ghép ma tr... >> Q=ones(4) Q= 1 1 1 1 1 1 1 1 13.L nh rank(A) : Tính h ng c a ma tr n: Ví d >> A=[1 1;2 2;3 3] A= 1 2 3 >> rank(A) ans = 14 L nh trace(A) : Tính v t c a ma tr n: Ví d >> A=[1 1;2 2;3 3] A= 1... cross(u,v) : Tích h Ví d : u= [ 3] u= >> v=[ 8] v= >> cross(u,v) ans = ng c a u, v 25 -10 L nh length : Tính chi u dài c Ví d x = [0 9] l = length (x) l = 10 » x = [01 09 77; 20 04 2001 ] x= 77 20 2001