Họ tên : Nguyễn Trung Đức Lớp : 63CLC1 MSSV : 52463 BÀI TẬP LỚN CƠNG CỤ MƠ HÌNH HOÁ BẰNG SỐ I Làm quen Matlab Bài 1.1 function phan1 disp ('PHAN I: LAM QUEN VOI MATLAB') disp ('Bai 1.1') disp ('Cau a') 2/2*3 ĐÁP ÁN >> Bai11 PHAN I: LAM QUEN VOI MATLAB Bai 1.1 Cau a ans = II Vector Bài 2.4 function Phan2 disp ('PHAN 2: VECTO) disp('Bai2.4') x=[3;2;6;8] y=[4;1;3;5] disp('Cau c: Chia cac phan tu cua y tuong ung voi cac phan tu cua x') disp('Ket qua:') disp(y./x) end ĐÁP ÁN >> Phan2 Bai2.4 x= y= Cau c: Chia cac phan tu cua y tuong ung voi cac phan tu cua x Ket qua: 1.3333 0.5000 0.5000 0.6250 III Ma trận Bài 3.3 function PHAN3 disp('Bai3.3') disp('Cho ma tran A = [2 ; ; 9]') A = [2 ; ; 9] disp(' Gan cho vector x la dong thu nhat cua A') x=A(1,1:3) disp('Gan cho ma tran y la hai dong lai (cuoi) cua A') y=A(2:3,:) ĐÁP ÁN >> PHAN3 Bai3.3 Cho ma tran A = [2 ; ; 9] A= Gan cho vector x la dong thu nhat cua A x= Gan cho ma tran y la hai dong lai (cuoi) cua A y= IV Vẽ đồ thị Bài 4.1 function PHAN4 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 >> PHAN4 Bai4.1 Ve thi ham so f(x) = sin(1/x) voi 0.01 < x < 0.1 Do thi ham so f(x) la: V BIỂU THỨC LOGIC Bài 5.5 function Phan5 disp('Bai 5.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 la:') a=X([X]>0) disp('Nhung phan tu am cua ma tran la:') b=X(X> Phan5 Bai 5.5 Cho ma tran X: X= -2 25 -11 -3 Nhung phan tu duong cua ma tran la: a= 25 Nhung phan tu am cua ma tran la: b= -2 -11 -3 Tong cac phan tu le cua ma tran la: S= 37 Gia tri lon nhat cua cac phan tu chan cua ma tran la: Max = Nhan phim bat ki de tiep tuc e= -2 25 -11 -3 Ma tran X sau gan nhung phan tu bang la: e= -2 25 -11 -3 VI CÂU LỆNH IF Bài 6.5 function Phan6 disp('Bai 6.5') T=input('Nhap vao gia tri cua T: '); if (T>0 & 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') end ĐÁP ÁN >> Phan6 Bai 6.5 Nhap vao gia tri cua T: Gia tri cua h(T)la: ans = -2 VII VÒNG LẶP Bài 7.6 function Phan7 disp('Bai 7.6') 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) end ĐÁP ÁN >> Phan7 Bai 7.6 Cau a: a(ij)=x(i)y(j) Ket qua: a(ij)= 24 28 36 12 42 VIII 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 Phan9 Bai 9.5 Hinh cau tam O ban kinh bang 3: ... 86900 So tien ho thu hai phai tra: 4.3450e+05 So tien ho thu ba phai tra: 7.8650e+05 IX Bài tập Symbolic math Bài 9.5 function Phan9 disp('Bai 9.5') disp('Hinh cau tam O ban kinh bang 3:') syms a... vector x la dong thu nhat cua A x= Gan cho ma tran y la hai dong lai (cuoi) cua A y= IV Vẽ đồ thị Bài 4.1 function PHAN4 disp('Bai4.1') disp(' Ve thi ham so f(x) = sin(1/x) voi 0.01 < x < 0.1')... Bai4.1 Ve thi ham so f(x) = sin(1/x) voi 0.01 < x < 0.1 Do thi ham so f(x) la: V BIỂU THỨC LOGIC Bài 5.5 function Phan5 disp('Bai 5.5') disp('Cho ma tran X:') X=[1 -2 9;5 25;1 -11 7;0 -3 2] disp('Nhung