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

BÀI TẬP LỚN MƠN CƠNG CỤ MƠ HÌNH HĨA BẰNG SỐ GIẢNG VIÊN: DƯƠNG THANH QUỲNH Họ Và Tên : KHÚC VĂN TÔ Lớp: 63CLC1 MSSV : 206263 Bài 1.1: Tính tay biểu thức sau thử lại matlab function Bai11 disp('Bai1.1') disp('Phan b') 6-2/5+7^2-1 disp('ket thuc bai 1') Đáp Án >> Bai11 Bai1.1 Phan b ans = 53.6000 ket thuc bai II Vector Bài 2.2: Cho x= [3 6] dự đốn kết dịng lệnh sau thử lại matlab function Bai22 disp('Bai2.2') disp('cho x=[ 6]') x=[ 6] disp('phan b') x(1:7) disp('ket thuc bai 2') ĐÁP ÁN >> Bai22 Bai2.2 cho x=[ 6] x = 9 phan b ans = ket thuc bai Bài3: Ma Trận function Bai35 disp('Bai 3.5:Cho vecto x=[2 6],y=[5 0]') disp('Cau c: Hay tao ma tran C co tinh chat: cot va co gia tri la vecto x, cot va co gia tri la vecto y') x=[2 6]; y=[5 0]; M=[x;y;x;y]; disp('Ket qua:') disp(M.') end Đáp án: >> bai3 Bai 3.5:Cho vecto x=[2 6],y=[5 0] Cau c: Hay tao ma tran C co tinh chat: cot va co gia tri la vecto x, cot va co gia tri la vecto y Ket qua: 5 9 1 1 6 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 bai55 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> bai55 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 Bài 6.1: Câu lệnh IF function bai61 disp('bai6.1') disp(' Du doan ket qua va thu lai bang matlab') disp(' tinh m:') n=input(' nhap gia tri cua n: n=') if n>1 m=n+1 else m=n-1 end disp(' an tiep phim bat ky de lam phan b') pause n=input(' nhap gia tri cua n: n=') if n>1 m=n+1 else m=n-1 end disp('an tiep phim bat ky de lam phan c') pause n=input(' nhap gia tri cua n: n=') if n>1 m=n+1 else m=n-1 end end ĐÁP ÁN >> bai61 bai6.1 Du doan ket qua va thu lai bang matlab tinh m: nhap gia tri cua n: n=4 n = m = an tiep phim bat ky de lam phan b nhap gia tri cua n: n= n = [] m = [] an tiep phim bat ky de lam phan c nhap gia tri cua n: n=6 n = m = Bài 7.2 : Phần lập trình Function bai72 r=input('to hop chap r, nhap r='); n=input('cua n phan tu, nhap n='); C=1; for i=1:n C=C*i; end for i=1:r C=C/i; end for i=1:(n-r) C=C/i; end; fprintf('to hop chap %d cua %d phan tu la C=%d\n',r,n,C) PHẦN ĐÁP ÁN >> bai72 to hop chap r, nhap r=20 cua n phan tu, nhap n=2 to hop chap 20 cua phan tu la C=8.220635e-19 Bai : Lập trình function bai81 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 Bai11 Bai1.1 Phan b ans = 53.6000 ket thuc bai II Vector Bài 2.2: Cho x= [3 6] dự đốn kết dịng lệnh sau thử lại matlab function Bai22 disp('Bai2.2') disp('cho... x(1:7) disp('ket thuc bai 2') ĐÁP ÁN >> Bai22 Bai2.2 cho x=[ 6] x = 9 phan b ans = ket thuc bai Bài3 : Ma Trận function Bai35 disp('Bai 3.5:Cho vecto x=[2 6],y=[5 0]') disp('Cau c: Hay tao ma tran