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

TRƯỜNG ĐẠI HỌC XÂY DỰNG Ban Đào Tạo Và Quản Lý Kỹ Sư Chất Lượng Cao BÀI TẬP LỚN CÔNG CỤ MƠ HÌNH HĨA BẰNG SỐ Giáo viên hướng dẫn : Dương Thanh Quỳnh Sinh viên thực : Bùi Thị Thùy Dương Mã số sinh viên : 1509063 Lớp quản lý : 63CLC1 1.Làm quen Matlab 1.1 Tính tay biểu thức sau, thử lại matlab function phan1 disp ('PHAN I: LAM QUEN VOI MATLAB') disp ('Bai 1.1: Tinh roi thu lai voi Matlab') disp ('Cau a') 2/2*3 disp ('Cau b') 6-2/5+7^2-1 disp ('Cau c') 10/2\5-3+2*4 disp ('Cau d') 3^2/4 disp ('Cau e') 3^2^2 disp ('Cau f') 2+round(6/9+3*2)/2-3 disp ('Cau g') 2+floor(6/9+3*2)/2-3 disp ('Cau h') 2+ceil(6/9+3*2)/2-3 disp ('Cau k') fix(4/9)+fix(3*(5/6)) ĐÁP ÁN PHAN I: LAM QUEN VOI MATLAB Bai 1.1: Tinh roi thu lai voi Matlab Cau a ans = Cau b ans = 53.6000 Cau c ans = Cau d ans = 2.2500 Cau e ans = 81 Cau f ans = 2.5000 Cau g ans = Cau h ans = 2.5000 Cau k ans = 2 VECTOR Bài 2.2 function Bai22 disp('Bai2.2') disp('cho x=[ 6]') x=[ 6] 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) ĐÁP ÁN Bai2.2 cho x=[ 6] x= 9 phan a ans = phan b ans = phan c ans = 9 3 phan d ans = phan e ans = phan f ans = phan g ans = 33 Ma trận Bài 3.3 function Bai33 disp('Bai3.3') disp('Cho ma tran A = [2 ; ; 9]') A = [2 ; ; 9] disp(' Gan cho vector x dong thu nhat cua A') x=A(1,1:3) disp('Gan cho ma tran y hai dong lai (cuoi) cua A') y=A(2:3,:) disp('Tinh tong theo dong ma tran A ') sum(A) disp('Tinh tong theo cot ma tran A ') sum(A') disp('gia tri lon nhat') max(max(A)) disp('gia tri nho nhat') min(min(A)) disp('Tính tong phan tu cua A ') sum(sum(A)) ĐÁP ÁN Bai3.3 Cho ma tran A = [2 ; ; 9] A= Gan cho vector x dong thu nhat cua A x= Gan cho ma tran y hai dong lai (cuoi) cua A y= Tinh tong theo dong ma tran A ans = 11 16 12 Tinh tong theo cot ma tran A ans = 15 17 gia tri lon nhat ans = gia tri nho nhat ans = Tính tong phan tu cua A ans = 39 Vẽ đồ thị Bài 4.1 function Bai41 disp('Bai4.1') disp(' Ve thi ham so f(x) = sin(1/x) voi 0.01 < x < 0.1') x=0.01:0.0001:0.1; disp(' Do thi ham so f(x) la:') plot(x,sin(1./x)) ĐÁP ÁN Bai4.1 Ve thi ham so f(x) = sin(1/x) voi 0.01 < x < 0.1 Do thi ham so f(x) la: Biểu thức logic Bài 5.5 function Bai55 disp('Bai5.5') disp(' Cho ma tran X') X=[1 -2 9;5 25;1 -11 7;0 -3 2] disp(' nhung phan tu duong cua ma tran:') X(X>0) disp(' nhung phan tu am cua ma tran:') X(X0 & T100) disp('Gia tri cua h(T) la:') 0.45*T+900 else disp('Khong xac dinh duoc gia tri cua ham h(T)!Vui long nhap gia tri khac') ĐÁP ÁN Nhap vao gia tri cua T: Gia tri cua h(T)la: ans = -5 Nhap vao gia tri cua T: 110 Gia tri cua h(T) la: ans = 949.5000 Vòng lặp Bài 7.6 function Bai 7.6 disp('Bai 7.6 : Cho x=[4 6] va y=[6 7], tinh cac mang/vecto/ma tran sau') disp('Cau a: a(ij)=x(i)y(j)') x=[4 6]; y=[6 7]; for(i=1:3) for(j=1:3) a(i,j)=x(i)*y(j); end end disp('Ket qua:') disp('a(ij)=') disp(a) ĐÁP ÁN Bai 7.6 : Cho x=[4 6] va y=[6 7], tinh cac mang/vecto/ma tran sau Cau a: a(ij)=x(i)y(j) Ket qua: a(ij)= 24 28 36 12 42 Lập trình Bài 8.1 function Phan8 disp('Bai 8.1:Tien dien') disp('So tien ho thu nhat phai tra:') disp(tiendien(120)) disp('So tien ho thu hai phai tra:') disp(tiendien(320)) disp('So tien ho thu ba phai tra:') disp(tiendien(480)) end function td=tiendien(x) if x