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

TRƯỜNG ĐẠI HỌC XÂY DỰNGBan Đào Tạo Và Quản Lý Kỹ Sư Chất Lượng Cao BÀI TẬP LỚ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 thực hiện : Bùi Thị Thùy Dươ

TRƯỜNG ĐẠI HỌC XÂY DỰNG

Ban Đào Tạo Và Quản Lý Kỹ Sư Chất Lượng Cao

BÀI TẬP LỚ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 thực hiện : Bùi Thị Thùy Dương

Mã số sinh viên : 1509063

Lớp quản lý : 63CLC1

1.Làm quen Matlab

1.1 Tính bằng tay các biểu thức sau, rồi thử lại bằng matlab

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))

ĐÁP ÁN

PHAN I: LAM QUEN VOI MATLAB

Bai 1.1: Tinh roi thu lai voi Matlab

Cau a

ans =

3

Cau b

ans =

53.6000

Cau c

ans =

6 Cau d ans = 2.2500 Cau e ans = 81 Cau f ans = 2.5000 Cau g ans = 2 Cau h ans = 2.5000 Cau k ans = 2

2 VECTOR

Bài 2.2

function Bai22

disp( 'Bai2.2' )

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

x=[ 3 1 5 7 9 2 6]

disp( 'phan a' )

x(3)

disp( 'phan b' )

x(1:7)

disp( 'phan c' )

x(1:end)

disp( 'phan d' )

x(1:end-1)

disp( 'phan e' )

x(6:-2:1)

disp( 'phan f' )

x([1 6 2 1 1])

disp( 'phan g' )

sum(x)

ĐÁP ÁN

Bai2.2

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

x =

3 1 5 7 9 2 6

phan a

ans =

5

phan b

ans =

3 1 5 7 9 2 6

phan c

ans =

3 1 5 7 9 2 6

phan d

ans =

3 1 5 7 9 2

phan e

ans =

2 7 1

phan f

ans =

3 2 1 3 3

phan g

ans =

33

3 Ma trận

Bài 3.3

function Bai33

disp( 'Bai3.3' )

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

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

disp( ' Gan cho vector x là dong thu nhat cua A' )

x=A(1,1:3)

disp( 'Gan cho ma tran y là hai dong con lai (cuoi) cua A' ) y=A(2:3,:)

disp( 'Tinh tong theo dong ma tran A ' )

sum(A)

disp( 'Tinh tong theo cot ma tran A ' )

sum(A')

disp( 'gia tri lon nhat' )

max(max(A))

disp( 'gia tri nho nhat' )

min(min(A))

disp( 'Tính tong các phan tu cua A ' )

sum(sum(A))

ĐÁP ÁN

Bai3.3

Cho ma tran A = [2 4 1 ; 6 7 2 ; 3 5 9]

A =

2 4 1

6 7 2

3 5 9

Gan cho vector x là dong thu nhat cua A

x =

2 4 1

Gan cho ma tran y là hai dong con lai (cuoi) cua A

y =

6 7 2

3 5 9

Tinh tong theo dong ma tran A

ans =

11 16 12

Tinh tong theo cot ma tran A

ans =

7 15 17

gia tri lon nhat

ans =

9

gia tri nho nhat

ans =

1

Tính tong các phan tu cua A

ans =

39

4 Vẽ đồ thị

Bài 4.1

function Bai41

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))

ĐÁP ÁN

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:

5 Biểu thức logic

Bài 5.5

function Bai55

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)

disp( ' nhung phan tu am cua ma tran:' )

X(X<0)

disp( ' tong cac phan tu le cua ma tran:' )

D(1:16)=X(1:4,1:4);

G=-100:2:100;

K=setdiff(D,G);

sum(K)

disp( 'gia tri lon nhat cac phan tu chan cua ma tran:' ) J=-101:2:101;

max(setdiff(D,J))

disp( 'gan nhung phan tu 0 thanh 1:' )

X(X==0)=1

ĐÁP ÁN

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

nhung phan tu am cua ma tran:

ans =

-2

-11

-3

tong cac phan tu le cua ma tran:

ans =

36

gia tri lon nhat cac phan tu chan cua ma tran:

ans =

8

gan nhung phan tu 0 thanh 1:

X =

1 -2 8 9

5 6 1 25

1 -11 3 7

1 4 -3 2

6 Câu lệnh if

Bài 6.5

function Bai65

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' )

ĐÁP ÁN

Nhap vao gia tri cua T: 5

Gia tri cua h(T)la:

ans =

-5

Nhap vao gia tri cua T: 110

Gia tri cua h(T) la:

ans =

949.5000

7 Vòng lặp

Bài 7.6

function Bai 7.6

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)

ĐÁP ÁN

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

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

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

9 Bài tập Symbolic math

Bài 9.5

function Phan9

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;

ĐÁP ÁN

Bai 9.5: Hinh cau tam O ban kinh bang 3