Đang tải... (xem toàn văn)
Lệnh maxX, minXPhần 4 ĐÁNH GIÁ NHẬN XÉT CỦA GIÁO VIÊNPhần 1 SỐ PHỨC TRONG MATLABTa có i là đơn vị phức >> i2 ans=-11.. Lệnh real, imag1.1 Ý nghĩaReal: lấy phần thực của số phứcImag: lấy
ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA KHOA CƠ KHÍ BÁO CÁO BÀI TẬP LỚN ểu Ti ĐẠI SỐ TUYẾN TÍNH ận lu ôn m LỚP CK16CK04 NHÓM c họ Giáo viên hướng dẫn: Nguyễn Xuân Mỹ Năm học: 2016-2017 DANH SÁCH THÀNH VIÊN ểu Ti Phạm Văn Phú Đinh Hoàng Việt Phạm Hữu Bảo Hồ Trần Trung Nghĩa Phan Nguyễn Nhật Huy Trương Công Duy Nguyễn Trọng Lương Nguyễn Xuân Trường Vy Đức Kiệm 10 Lê Minh Hiếu 11 Phan Duy Khương MSSV ận lu ôn m 1612625 1614091 1610194 1612199 1611317 1610509 1611948 1613846 1611691 1611015 1611675 c họ MỤC LỤC Phần SỐ PHỨC TRONG MATLAB ểu Ti Lệnh real, imag 1.1 Ý nghĩa 1.2 Cú pháp 1.3 Ví dụ Lệnh abs 2.1 Ý nghĩa 2.2 Cú pháp 2.3 Ví dụ Lệnh angle 3.1 Ý nghĩa 3.2 Cú pháp 3.3 Ví dụ Lệnh conj 4.1 Ý nghĩa 4.2 Cú pháp 4.3 Ví dụ ận lu ơn m 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 họ Phần MA TRẬN TRONG MATLAB ểu Ti 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ệnh A’ 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) [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 ận lu ôn m c họ Phần MỘT SỐ LỆNH LẦN KHÔNG GIAN VECTOR, KHÔNG GIAN EUCLIDE, TRỊ GIÊNG Lệnh dot Lệnh cross Lệnh length Lệnh norm Lệnh qr Lệnh [P,D]=eig(A) Lệnh eig(H) Lệnh max(X), min(X) Phần ĐÁNH GIÁ NHẬN XÉT CỦA GIÁO VIÊN ểu Ti PHỨC TRONG MATLAB lu Phần SỐ ận Ta có i đơn vị phức >> i2 ans=-1 ôn m Lệnh real, imag 1.1 Ý nghĩa Real: lấy phần thực số phức Imag: lấy phần ảo 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 Ý nghĩa: tìm modul số phức c họ 2.2 Cú pháp: y=abs(z) 2.3 Ví dụ >>z=3+4i z= ểu Ti 3.000 + 4.000i >> Modul=abs(z) Modul= Lệnh angle 3.1 Ý nghĩa: Tìm agument số phức với đơn vị radian 3.2 Cú pháp: y=angle(z) 3.3 Ví dụ >> z= 3+4i z= 3.0000 + 4.0000i >> agumen=angle(z) ôn họ Lệnh conj m 0.9273 ận lu agumen = c 4.1 Ý nghĩa: Lấy số phức liên hợp số phức 4.2 Cú pháp: y= conj(z) 4.3 Ví dụ >> z=3+4iz = 3.0000 + 4.0000i >> conj(z) ans = 3.0000 - 4.0000i ểu Ti lu TRẬN TRONG MATLAB ận Phần MA ôn m Lệnh numel(A): Đếm số phần tử a c họ 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 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: - C ma trận cần tìm ểu Ti - A,B ma trận cho trước b, Ví dụ: Cho ma trận A= [1 2], B=[ 7] , Ghép ma trận A,B theo cột họ ôn m ận C= lu >> C =[1 2;3 4] c Lệnh inv(A): Tìm ma trận nghịch đảo ma trận Ví dụ : >> A=[1 2;2 5] A= 2 >> inv(A) ans = -2 -2 Lệnh Ak: Với: - A ma trận cho k hệ số mũ cần tính Ví dụ: Bài 1: Cho ma trận A=[1 2;2 4] >>A = 2 >> B=A3 B= 25 50 ểu Ti 50 100 Câu lệnh [A B] :Ghép ma trận theo hàng ận lu -Cú pháp:[A B] -Ví dụ: >> B=[3 5;6 2;3 3;5 9] B= 6 c họ ôn A= m >> A=[1 4;5 5;3 1;6 1] 5 >> [A B] ans = 4 5 6 8 ểu Ti Câu lệnh A(:,n)=[ ] : Xóa cột thứ n ma trận A >> B=[3 5;6 2;3 3;5 9] 6 5 c họ ôn m ận lu B= >> B(:,1)=[] B= Câu lệnh A(:,n:end) : Cho phép lấy từ cột thứ n đến cột cuối ma trận -Cú pháp:A(:,n:end) A= >> A(2,3) ểu Ti ans = lu ận 24.Lệnh A(i,:) A(:,j) : tham chiếu dòng i tham chiếu cột j 5 8 >> A(2,:) ans = c họ ôn A= m Ví dụ: Cho A=[ 5; 6; 3] >> A(:,3) ans = ểu Ti 25 Lệnh A(i :k, :) và A( :,j :k) : Tham chiếu từ dòng i dến dòng k Tham chiếu từ cột j đến cột k ận lu Ví dụ A=[ 2; 8; 3] ôn m A= >> A(1:2,:) ans = 5 >> A(:,1:2) ans = c họ 1 5 26 Lệnh rref(A) : Tạo ma trận bậc thang từ A Ví dụ A=[ 2; 1; 5] A= ận lu ểu Ti 1 0 c họ ôn ans = m >> rref(A) 27.Lệnh FLIPLR : Chuyển phần tử ma trận theo thứ tự cột ngược lại - Cú pháp: b = fliplr(a) - Giải thích: b: tên ma trận chuyển đổi a: tên ma trận cần chuyển đổi Ví dụ: a= 9 b = fliplr(a) 28 Lệnh FLIPUD : Chuyển phần tử ma trận theo thứ tự hàng ngược lại ểu Ti - Cú pháp: b = flipud(a) - Giải thích: lu b: tên ma trận chuyển đổi ận a: tên ma trận cần chuyển đổi Ví dụ: c a= họ >>a= [ 4; 5; 6] ôn m - >>b = flipud(a) b= 29 Lệnh MAGIC : Tạo ma trận vng có tổng phần tử hàng, cột đường chéo - Cú pháp: Tên ma trận = magic(n) - Giải thích: n: kích thước ma trận Giá trị phần tử ma trận dãy số nguyên liên tục từ đến 2n Tổng hàng, cột đường chéo >>tmt = magic(3) ểu Ti Ví dụ: c họ ôn m ận lu 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 3 4 10 10 20 31 Lệnh RAND : Tạo ma trận mà kết mà giá trị phần tử ngẫu nhiên - Cú pháp: y = rand(n) y = rand(m,n) - y: tên ma trận ểu Ti Giải thích: lu -n: tạo ma trận có n hàng, n cột ận -m, n: tạo ma trận có m hàng, n cột m ơn - Giá trị phần tử nằm khoảng [0 1] c >>y = rand(3) họ Ví dụ: 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ận xoay 900 a: ma trận cần xoay >>a=[1 3; 6; 9] ểu Ti Ví dụ: lu a= họ ôn m ận c >> 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ận cho trước , ans =0 nghĩa A ma trận rỗng, ans =1 ma trận cho ma trận rộng - Ví dụ >> B= zeros(4) B= 0 0 0 0 0 0 0 0 % B không ma trận rỗng ôn m >> A= [] ận lu ểu ans = Ti >>isempty(B) A= họ [] c >>isempty(A) ans = % A ma trận rỗng 34.Lệnh DIAG : Tạo ma trận xử lý đường chéo theo quy ước - Cú pháp: v = diag(x) v = diag(x,k) > Giải thích: - x: vector có n phần tử - v: ma trận tạo từ x theo quy tắc: số hàng số cột phần tử x nằm đường chéo v - k: tham số định dạng cho v, số hàng cột v = n + abs(k) - Nếu k = đường chéo v phần tử x - Nếu k > phần tử x nằm phía đường chéo v - Nếu k < phần tử x nằm phía đường chéo v ểu Ti >>x = [ ận lu Ví dụ: 4]; ôn m v = diag(x) v= họ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 c >>v1 = diag(x,2) v1 = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 >>v2 = diag(x,0) v2 = ểu Ti 0 0 lu 0 ận >>v3 = diag(x,-2) ôn m v3 = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 c họ Phần MỘT SỐ LỆNH LẦN KHÔNG GIAN VECTOR, KHÔNG GIAN EUCLIDE, TRỊ GIÊNG Lệnh dot:( tính tích vơ hướng vectơ) -Cú pháp :dot(u,v) - u,v:hai vectơ cho trước Ví dụ: >> u=[1 3] u= ểu v= Ti >> v=[3 5] ận >> dot(u,v) lu m ôn ans = 26 c họ Lệnh cross(u,v) : Tích hữu hướng u, v Ví dụ : u= [ 3] u= >> v=[ 8] v= >> cross(u,v) ans = 25 -10 Lệnh length : Tính chiều dài vectơ Ví dụ x = [0 9] l = length (x) l = 10 » x = [01 09 77; 20 04 2001 ] x= 77 20 2001 lu » l=length(x) ận l= ểu Ti c Ví dụ: u=[1 3] họ -Cú pháp: norm(u) ; u vec tơ ơn m 4.Lệnh norm:(tính độ dài vectơ bất kì) u= >> norm(u) ans = 3.7417 5.Lệnh qr:(trực chuẩn hóa họ vectơ cột A) [P,]=qr(A) ; A ma trận cột tạo họ véc tơ Ví dụ: Trong R^3 cho véc tơ u=(1,1,1).Hãy trực chuẩn véc tơ u >> A=[1 1] A= 1 >> qr(A) ans = 1.7321 ểu Ti -1.3660 -1.3660 lu 6.Lệnh [P,D]=eig(A) : Chéo hóa ận Ví dụ -1 -8 -14 c -8 họ 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= ôn A= m A=[0 -8 6;-1 -8 7;1 -14 11] -2.0000 0 3.0000 0 2.0000 Lệnh eig(H) : Trị riêng Ví dụ H=[1 2;3 4] H= ểu Ti >> eig(H) lu ans = ận -0.3723 m ôn 5.3723 X= 4 >> max(X) ans = >> min(X) c Ví dụ >> X=[ 9; 3; 4] họ Lệnh max(X), min(X) : Trả giá trị lớn nhỏ vector X ans = 2 Phần ĐÁNH GIÁ NHẬN XÉT CỦA GIÁO VIÊN ểu Ti ận lu ôn m c họ