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

8.Lập trìnhfunction Bai8 disp'Bai 8.1:Tien dien' disp'So tien ho thu nhat phai tra:' disptiendien120 disp'So tien ho thu hai phai tra:' disptiendien320 disp'So tien ho thu ba phai tra:'

Bài tập Công cụ mô hình hóa bằng số

Giảng viên : DƯƠNG THANH QUỲNH

Họ và tên: Vũ Đức Khang

Lớp:63CLC1 MSSV:107063

1.Làm quen Matlab

function Bai1

disp('Bai1.1.a')

disp('2/2*3')

disp('Ket qua')

disp(2/2*3)

Đáp án

>> Bai1

Bai1.1.a

2/2*3

Ket qua

3

2.Vec to

function Bai2

disp('Bai 2.2.a')

disp('Cho x=[3 1 5 7 9 2 6]')

disp('Du doan dong lenh sau va thu lai bang Matlab') disp('x(3)')

x=[3 1 5 7 9 2 6];

disp('Dap an')

x(3)

Đáp án

>> Bai2

Bai 2.2.a

Cho x=[3 1 5 7 9 2 6]

Du doan dong lenh sau va thu lai bang Matlab

x(3)

Dap an

ans =

5

3.Ma trận

function Bai3

disp('Bai 3.2.d Du doan ket qua giai thich va thu lai bang Matlab ')

disp('Cho A=[2 7 9 7;3 1 5 6;8 1 2 5]')

disp('reshape(A,2,6)')

disp('Ket qua')

A=[2 7 9 7;3 1 5 6;8 1 2 5]

reshape(A,2,6)

Đáp án

>> Bai3

Bai 3.2.d Du doan ket qua giai thich va thu lai bang Matlab

Cho A=[2 7 9 7;3 1 5 6;8 1 2 5]

reshape(A,2,6)

Ket qua

A =

2 7 9 7

3 1 5 6

8 1 2 5

ans =

2 8 1 9 2 6

3 7 1 5 7 5

4 Vẽ đồ thị

function Bai4

disp('Bai4.1')

disp(' Ve do 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))

disp('ket thuc bai 4')

Đáp án

>> Bai4

Bai4.1

Ve do thi ham so f(x) = sin(1/x) voi 0.01 < x < 0.1

Do thi ham so f(x) la:

Bài 5: Biểu thức logic

function Bai5

disp('Bai5.5')

disp(' Cho ma tran X')

X=[1 -2 8 9;5 6 0 25;1 -11 3 7;0 4 -3 2] disp(' nhung phan tu duong cua ma tran:') X(X>0)

Đáp án

>> Bai5

Bai5.5

Cho ma tran X

X =

1 -2 8 9

5 6 0 25

1 -11 3 7

0 4 -3 2

nhung phan tu duong cua ma tran:

ans =

1

5

1

6

4

8

3

9

25

7

2

6.Câu lệnh If

function Bai6

disp('Bai 6.5')

T=input('Nhap vao gia tri cua T: ');

if (T>0 & T<100)

disp('Gia tri cua h(T)la:')

T-10

elseif (T>100)

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

>> Bai6

Bai 6.5

Nhap vao gia tri cua T: 5

Gia tri cua h(T)la:

ans =

-5

7.Vòng lặp

function Bai7

disp('Bai 7.6 : Cho x=[4 1 6] va y=[6 2 7], tinh cac mang/vecto/ma tran sau')

disp('Cau a: a(ij)=x(i)y(j)')

x=[4 1 6];

y=[6 2 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

>> Bai7

Bai 7.6 : Cho x=[4 1 6] va y=[6 2 7], tinh cac mang/vecto/ma tran sau

Cau a: a(ij)=x(i)y(j)

Ket qua:

a(ij)=

24 8 28

6 2 7

36 12 42

8.Lập trình

function Bai8

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<=100

td=550*x*1.1;

elseif (100<x)&&(x<=200)

td=(550*100+(x-100)*1200)*1.1;

elseif (200<x)&&(x<=300)

td=(550*100+1200*100+(x-200)*1800)*1.1;

elseif x>300

td=(550*100+1200*100+1800*100+(x-300)*2000)*1.1; end

end

Đáp án

>> Bai8

Bai 8.1:Tien dien

So tien ho thu nhat phai tra:

86900

So tien ho thu hai phai tra:

4.3450e+05

So tien ho thu ba phai tra:

7.8650e+05

Đáp án

9.Bài tập Symbolic Math

function Bai9

disp('Bai 9.5: 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

Đáp án

>> Bai9

Bai 9.5: Hinh cau tam O ban kinh bang 3