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

TRƯỜNG ĐẠI HỌC XÂY DỰNG BAN QUẢN LÝ VÀ ĐÀO TẠO KỸ SƯ CHẤT LƯỢNG CAO BÀI TẬP LỚN CÔNG CỤ HĨA BẰNG MƠ HÌNH SỐ GIÁO VIÊN HƯỚNG DẪN : QUỲNH SINH VIÊN THỰC HIỆN DƯƠNG THANH : PHÙNG NGỌC HÀ MY MÃ SỐ SINH VIÊN LỚP QUẢN LÝ : 1532463 : 63CLC1 I LÀM QUEN VỚI MATLAB Bài 1.1: Tính tay thử lại Matlab PHẦN LẬP TRÌNH 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)) PHẦN ĐÁP ÁN >> Bai11 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 = II VECTOR Bài 2.2 : Cho ma trận PHẦN LẬP TRÌNH function Bai23 disp ('PHAN II: VECTOR') disp ('Bai 2.3') disp ('Nhap ma tran x') x=[2 6] disp ('Cong tat ca cac phan tu ma tran x voi 16') x+16 disp ('Cong cac phan tu le cua ma tran voi 3') y=[3 0] x+y disp ('Lay can binh phuong tat ca cac phan tu cua ma tran x') sqrt(x) disp ('Binh phuong tat ca cac phan tu') x.^2 PHẦN ĐÁP ÁN >> Bai23 PHAN II: VECTOR Bai 2.3 Nhap ma tran x x= Cong tat ca cac phan tu ma tran x voi 16 ans = 18 21 17 22 Cong cac phan tu le cua ma tran voi y= ans = 5 Lay can binh phuong tat ca cac phan tu cua ma tran x ans = 1.4142 2.2361 1.0000 2.4495 Binh phuong tat ca cac phan tu ans = 25 36 III MA TRẬN Bài 3.1 : a b c d x + y dòng lệnh hợp lệ, x y ma trận 3x1 x + A dịng lệnh hợp lệ x A ma trận có cột x’ + y dịng lệnh hợp lệ x’ ma trận 1x3 A – [x’ y’] dịng lệnh khơng hợp lệ A ma trận 3x2 mà [x’ y’] ma trận 2x3 e [x ; y’] dịng lệnh khơng hợp lệ x ma trận 3x1 mà y’ ma trận 1x3 nên khơng thể ghép dịng f [x ; y] dịng lệnh hợp lệ lệnh ghép ma trận 3x1 g A – dòng lệnh hợp lệ phép trừ tất phần tử ma trận A cho PHẦN LẬP TRÌNH function Bai31 disp('Nhap ma tran x, y, A') x=[1 8] y=[2 5] A=[3 ; 7] disp ('Thu xem dong lenh nao la hop le') x+y x+A x'+y A-[x' y'] [x,y'] [x,y] A-3 PHẦN ĐÁP ÁN >> Bai31 Nhap ma tran x, y, A x= y= A= Thu xem dong lenh nao la hop le ans = 13 14 6 15 6 ans = ans = 10 13 Matrix dimensions must agree Error in Bai31 (line 10) A-[x' y'] Error in Bai31 (line 10) [x,y'] ans = ans = -2 -1 IV VẼ ĐỒ THỊ Bài 4.4: Vẽ lại đồ thị hàm số ex (0 < x < 4) với đường đứt nét, màu đỏ dấu (marker) dấu + PHẦN LẬP TRÌNH function Bai44 disp ('Ve thi ham so') disp ('Khoang gia tri cua x la') x=0:4 disp ('Do thi ham so can ve la') plot(x,exp(x),' r+') PHẦN ĐÁP ÁN >> Bai44 Ve thi ham so khoang gia tri cua x la x= Do thi ham so can ve la >> Bai44 Ve thi ham so khoang gia tri cua x la x= Do thi ham so can ve la V BIỂU THỨC LOGIC Bài 5.5: PHẦN LẬP TRÌNH 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 PHẦN ĐÁP ÁN -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: PHẦN LẬP TRÌNH 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 PHẦN ĐÁP ÁN >> Bai65 Nhap vao gia tri cua T: Gia tri cua h(T)la: ans = -2 VII VỊNG LẶP Bài 7.2: PHẦN LẬP TRÌNH 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 VIII LẬP TRÌNH Bài 8.1: PHẦN LẬP TRÌNH A=input('nhap vao so dien tieu thu cua cac ho gia dinh:'); T=zeros(1,length(A)); for i=1:length(A) if A(i)