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

BÀI TẬP LỚN MƠ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 : TRẦN TIẾN DŨNG Lớp : 63CLC1 MSSV : 31063 LÀM QUEN VỚI MATLAB Bài 1.3 Cho x = 2, y = Dự đoán kết tiếp theo: z=x y = y +z x=y+x–z x+y–z function bai1_3 disp('x=2,y=3') x=2; y=3; z=x; disp(['+) z=x = ',num2str(z)]) y=y+z; disp(['+) y=y+z = ',num2str(y)]) x=y+x-z; disp(['+) x=y+x-z = ',num2str(x)]) g=x+y-z; disp(['+) x+y-z= ',num2str(g)]) >> bai1_3 x=2,y=3 +) z=x = +) y=y+z = +) x=y+x-z = +) x+y-z= Bài 1.5 Cho a =36 b =15 Tính tay biểu thức sau, thử lại Matlab a mod(a,b) = b rem(a,b) = c gcd(a,b) = d lcm(a,b) = 180 function bai1_5 disp('a = 36, b = 15') a=36; b=15; c=mod(a,b); disp(['mod(a,b)= ',num2str(c)]) d=rem(a,b); disp(['rem(a,b)= ',num2str(d)]) e=gcd(a,b); disp(['gcd(a,b)= ',num2str(e)]) f=lcm(a,b); disp(['lcm(a,b)= ',num2str(f)]) >> bai1_5 a = 36, b = 15 mod(a,b)= rem(a,b)= gcd(a,b)= lcm(a,b)= 180 VECTOR 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]) = 33 g sum(x) = 33 function bai2_2 disp('x = [3 6]') x = [3 6]; disp(['a) x(3)= ',num2str(x(3))]) disp(['b) x(1:7)= ',num2str(x(1:7))]) disp(['c) x(1:end)= ',num2str(x(1:end))]) disp(['d) x(1:end-1)= ',num2str(x(1:end-1))]) disp(['e) x(6:-2:1)= ',num2str(x(6:-2:1))]) disp(['f) x([1 1])= ',num2str(x([1 1]))]) disp(['g) sum(x)=',num2str(sum(x))]) >> bai2_2 x = [3 6] 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])= 3 g) sum(x)=33 Bài 2.3 cho x = [ ] a Cộng thêm 16 vào tất phần tử b Cộng thêm vào phần tử vị trí lẻ c Lấy bậc tất phần tử d Bình phương tất phần tử function bai2_3 disp('x = [2 6]') x=[2 6]; disp(['a)x = ',mat2str(x+16)]) y=x; y(1:2:end)=y(1:2:end)+3; disp(['b)x = ',mat2str(y)]) disp(['c)x = ',mat2str(sqrt(x))]) disp(['d)x = ',mat2str(x.*x)]) >> bai2_3 x = [2 6] a)x = [18 21 17 22] b)x = [5 6] c)x = [1.4142 2.2361 2.4495] d)x = [4 25 36] MA TRẬN Bài 3.3 Cho ma trận A = [2 ; ; 9], viết lệnh Matlab để a Gán cho vector x dòng thứ A b Gán cho ma trận y hai dòng lại (cuối) A c Tính tổng theo dịng ma trận A d Tính tổng theo cột ma trận A e Tìm giá trị lớn nhỏ ma trận f Tính tổng phần tử A function bai3_3 A=[2 ; ; 9]; disp(['A = ',mat2str(A)]) x=A(1,:); disp(['a) x = ',mat2str(x)]) y=A(end-1:end,:); disp(['b) y = ',mat2str(y)]) disp(['c) tong theo dong cua ma tran A la: ',mat2str(sum(A,2))]) disp(['d) tong theo cot cua ma tran A la: ',mat2str(sum(A,1))]) disp(['e) gia tri lon nhat cua ma tran A la: ',num2str(max(A(:)))]) disp([' gia tri nho nhat cua ma tran A la: ',num2str(min(A(:)))]) disp(['f) tong cac phan tu cua A la: ',num2str(sum(A(:)))]) >> bai3_3 A = [2 1;6 2;3 9] a) x = [2 1] b) y = [6 2;3 9] c) tong theo dong cua ma tran A la: [7;15;17] d) tong theo cot cua ma tran A la: [11 16 12] e) gia tri lon nhat cua ma tran A la: gia tri nho nhat cua ma tran A la: f) tong cac phan tu cua A la: 39 VẼ ĐỒ THỊ Bài 4.7 Vẽ hai hàm y = x2 y sin3x đồ thị function bai4_7 x = linspace(-2*pi,2*pi); y1 = x.^2; y2 = sin(x); figure plot(x,y1,'rx-',x,y2,'b. ') xlabel('x') ylabel('y') title('do thi ham x^2 va sin(x)') legend ('y1 = x^2','y2 =sinx') axis([-4 -2 15]) BIỂU THỨC LOGIC Bài 5.5 Cho ma trận X = [1 -2 9; 25; -11 7; -3 2] tìm : a Những phần tử dương ma trận b Những phần tử âm ma trận c Tổng phần tử lẻ ma trận d Giá trị lớn phần tử chẵn ma trận e Gán phần tử thành function bai5_5 X=[1 -2 9; 25; -11 7; -3 2]; disp(['X = ',mat2str(X)]) F=X(X>0); disp(['a) nhung phan tu duong cua ma tran X: ',mat2str(F)]) H=X(X> bai5_5 X = [1 -2 9;5 25;1 -11 7;0 -3 2] a) nhung phan tu duong cua ma tran X: [1;5;1;6;4;8;3;9;25;7;2] b) nhung phan tu am cua ma tran X: [-2;-11;-3] c) tong cac phan tu le cua ma tran: 37 d) gia tri lon nhat cac phan tu chan cua ma tran: e) sau gan nhung phan tu 1: X=[1 -2 9;5 25;1 -11 7;1 -3 2] CÂU LỆNH IF Bài 6.5 h(T) = T - 10 = 0.45 T + 900 Thử trường hợp: a T = b T = 110 < T < 100 T > 100 function bai6_5 T=input('nhap gia tri cua T: '); if (T>0)&(T100 h=0.45*T+900; disp(['khi T = ',num2str(T),' thi ham h = ',num2str(h)]) end >> bai6_5 nhap gia tri cua T: T = thi ham h >> bai6_5 nhap gia tri cua T: T = 110 thi ham = -5 110 h = 949.5 VÒNG LẶP Bài 7.3 Tạo ma trận M x N số ngẫu nhiên (dùng lệnh rand) Chuyển giá trị nhỏ 0.2 thành 0, giá trị lớn hay 0.2 thành function bai7_3 M=input('nhap M: '); N=input('nhap N: '); for i=1:M; for j=1:N; X(i,j)=rand; end end disp(['ma tran X ngau nhien ',num2str(M),'x',num2str(N)]) disp(X) X(X=0.2)=1; disp('ma tran X sau chuyen cac gia tri nho hon 0.2 0,') disp('cac gia tri lon hon hay bang 0.2 1:') disp(X) >> bai7_3 nhap M: nhap N: ma tran X ngau nhien 2x5 0.2575 0.8407 0.2543 0.9293 0.3500 0.1966 ma tran X sau chuyen cac 0, cac gia tri lon hon hay bang 1 1 1 0.8143 0.2511 0.2435 0.6160 gia tri nho hon 0.2 0.2 1: 1 LẬP TRÌNH Bài 8.1 Đơn giá tính tiền điện tháng cơng ty điện lực hộ gia đình sau: Nếu lượng điện sử dụng bai8_1 nhap so dien su dung: 120 tien dien phai tra la: 86900 vnd >> bai8_1 nhap so dien su dung: 320 tien dien phai tra la: 434500 vnd >> bai8_1 nhap so dien su dung: 480 tien dien phai tra la: 786500 vnd BÀI TẬP SYMBOLIC MATH Bài 9.5 Vẽ hình cầu tâm O bán kính phương trình: function bai9_5 disp('Hinh cau tam O ban kinh bang 3:') syms a b; x=3*sin(a)*cos(b); y=3*sin(a)*sin(b); z=3*cos(a); ezsurf(x,y,z); axis tight ; colorbar; end >> bai9_5 Hinh cau tam O ban kinh bang 3: ... 434500 vnd >> bai8_1 nhap so dien su dung: 480 tien dien phai tra la: 786500 vnd BÀI TẬP SYMBOLIC MATH Bài 9.5 Vẽ hình cầu tâm O bán kính phương trình: function bai9_5 disp('Hinh cau tam O ban... T: T = 110 thi ham = -5 110 h = 949.5 VÒNG LẶP Bài 7.3 Tạo ma trận M x N số ngẫu nhiên (dùng lệnh rand) Chuyển giá trị nhỏ 0.2 thành 0, giá trị lớn hay 0.2 thành function bai7_3 M=input('nhap... / 1KW Từ 201 đến 300 KW tính 1.800 đ / 1KW Trên 300 KW 2.000 đ / 1KW Tổng số tiền mà hộ gia đình phải trả tổng số tiền sử dụng điện + 10% thuế VAT Lượng điện sử dụng tháng hộ gia đình khu