Bài tập lớn công cụ mô hình hóa bằng số (5)

XÂY DỰNG TRƯỜNG ĐẠI HỌC BAN ĐÀO TẠO KỸ SƯ CHẤT LƯỢNG CAO BÀI TẬP MÔN CÔNG CỤ MÔ HINH HOÁ BẰNG SỐ GIẢNG VIÊN: DƯƠNG THANH QUỲNH SINH VIÊN:LẠI MINH ĐỨC LỚP:63CLC1 MSSV:49863 Bài 1.1:Làm quen Matlab 1.1 Tính tay biểu thức sau, thử lại Matlab: a / * b – / + ^ – c 10 / \ – + * d ^ / e ^ ^ f + round(6 / + * 2) / – g + floor(6 / + * 2) / – h + ceil(6 / + * 2) / – k fix(4/9)+fix(3*(5/6)) GIẢI function Bai11 disp('++++++++++++++++++++++++++++++++++++++++++++++++++++++++') disp(' LAM QUEN VOI MATLAP ') disp('+++++++++++++++++++++++++++++++++++++++++++++++++++++++++') disp('Bai1.1') disp('phan a') 2/2*3 disp('phan b') 6-2/5+7^2-1 disp('phan c') 10/2\5-3+2*4 disp('phan d') 3^2/4 disp('phan e') 3^2^2 disp('phan f') 2+round(6/9+3*2)/2-3 disp('phan g') 2+floor(6/9+3*2)/2-3 disp('phan h') 2+ceil(6/9+3*2)/2-3 disp('phan k') fix(4/9)+fix(3*(5/6)) disp(' ket thuc bai ') end ++++++++++++++++++++++++++++++++++++++++++++++++++++++++ LAM QUEN VOI MATLAP +++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Bai1.1 phan a ans = phan b ans = 53.6000 phan c ans = phan d ans = 2.2500 phan e ans = 81 phan f ans = 2.5000 phan g ans = phan h ans = 2.5000 phan k ans = ket thuc bai Bài 2.2: Cho x = [3 6], dự đốn kết dịng lệnh sau thử lại Matlab: a x(3) b x(1:7) c x(1:end) d x(1:end-1) e x(6:-2:1) f x([1 1]) g sum(x) Giải function bai22 disp('=====================================================') disp('| VECTO |') disp('=====================================================') disp('Bai 2.2') disp('cho x=[3 6]') x=[ ] disp('phan a') x(3) disp('phan b') x(1:7) disp('phan c') x(1:end) disp('phan d') x(1:end-1) disp('phan e') x(6:-2:1) disp('phan f') x([1 1]) disp('phan g') sum(x) end >> bai2 ===================================================== | VECTO | ===================================================== Bai 2.2 cho x=[3 6] x= 9 phan a ans = phan b ans = phan c ans = phan d ans = 3 phan e ans = phan f ans = phan g ans = 33 >> Bài 3.2: cho A = [2 ; ; 5], dự đốn kết quả, giải thích; thử lại Matlab: a A’ b A(:,[1 4]) c A([2 3], [3 1]) d reshape(A, 2, 6) e A(:) f flipud(A) g fliplr(A) h [A; A(end,:)] i A(1:3,:) j [A; A(1:2, :)] k sum(A) l sum(A’) m sum(A, 2) n [ [ A ; sum(A) ] [ sum(A,2) ; sum(A(:)) ] ] Giải function Bai32 disp('=====================================') disp('| MA TRAN |') disp('=====================================') disp('Bai3.2') disp('cho ma tran A= [2 ; ; 5]') A= [2 ; ;8 5] disp('phan a') A' disp('phan b') A(:,[1 4]) disp('phan c') A([2 3],[3 1]) disp('phan d') reshape(A,2,6) disp('phan e') A(:) disp('phan f') flipud(A) disp('phan g') fliplr(A) disp('phan h') [A;A(end,:)] disp('phan i') A(1:3,:) disp('phan j') [A;A(1:2,:)] disp('phan k') sum(A) disp('phan l') sum(A') disp('phan m') sum(A,2) disp('phan n') [[A;sum(A)] [sum(A,2);sum(A(:))]] end >> bai3 ===================================== | MA TRAN | ===================================== Bai3.2 cho ma tran A= [2 ; ; 5] A= 8 1 phan a ans = phan b ans = phan c ans = phan d ans = 7 9 phan j ans = phan k ans = 13 16 18 phan l ans = 25 15 16 phan m ans = 25 15 16 phan n ans = 25 15 16 13 16 18 56 >> Bài 4.3: Vẽ đồ thị hàm số f(x) = sin(1/x) với 0.01 < x < 0.1 function Bai43 disp('==========================================') disp('| VE DO THI |') disp('==========================================') disp('Bai43') disp('ve thi ham so f(x)=sin(1/x) voi 0.01