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

Vector Bài 2.4 function Phan2... LẬP TRÌNH Bài 8.1 function Phan8 disp'Bai 8.1:Tien dien' disp'So tien ho thu nhat phai tra:' disptiendien120 disp'So tien ho thu hai phai tra:' disptiend

Họ và tên : Nguyễn Trung Đức

Lớp : 63CLC1

MSSV : 52463

BÀI TẬP LỚN CÔNG CỤ MÔ HÌNH HOÁ BẰNG SỐ

I Làm quen Matlab

Bài 1.1

function phan1

disp ('PHAN I: LAM QUEN VOI MATLAB')

disp ('Bai 1.1')

disp ('Cau a')

2/2*3

ĐÁP ÁN

>> Bai11

PHAN I: LAM QUEN VOI MATLAB

Bai 1.1

Cau a

ans =

3

II Vector

Bài 2.4

function Phan2

disp ('PHAN 2: VECTO)

disp('Bai2.4')

x=[3;2;6;8]

y=[4;1;3;5]

disp('Cau c: Chia cac phan tu cua y tuong ung voi cac phan tu cua x') disp('Ket qua:')

disp(y./x)

end

ĐÁP ÁN

>> Phan2

Bai2.4

x =

3

2

6

8

y =

4

1

3

5

Cau c: Chia cac phan tu cua y tuong ung voi cac phan tu cua x

Ket qua:

1.3333

0.5000

0.5000

0.6250

III Ma trận

Bài 3.3

function PHAN3

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 la dong thu nhat cua A')

x=A(1,1:3)

disp('Gan cho ma tran y la hai dong con lai (cuoi) cua A')

y=A(2:3,:)

ĐÁP ÁN

>> PHAN3

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 la dong thu nhat cua A

x =

2 4 1

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

y =

6 7 2

3 5 9

IV Vẽ đồ thị

Bài 4.1

function PHAN4

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

>> PHAN4

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:

V BIỂU THỨC LOGIC

Bài 5.5

function Phan5

disp('Bai 5.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 la:') a=X([X]>0)

disp('Nhung phan tu am cua ma tran la:')

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

S=sum(X(mod([X],2)==1))

disp('Gia tri lon nhat cua cac phan tu chan cua ma tran la:') Max=max(X(mod(X,2)==0))

disp('Nhan phim bat ki de tiep tuc')

pause

e=X

disp('Ma tran X sau khi gan nhung phan tu bang 0 thanh 1 la:') e(X==0)=1

ĐÁP ÁN

>> Phan5

Bai 5.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 la:

a =

1

5

1

6

4

8

3

9

25

7

2

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 =

8

Nhan phim bat ki de tiep tuc

e =

1 -2 8 9

5 6 0 25

1 -11 3 7

0 4 -3 2

Ma tran X sau khi gan nhung phan tu bang 0 thanh 1 la:

e =

1 -2 8 9

5 6 1 25

1 -11 3 7

1 4 -3 2

VI CÂU LỆNH IF

Bài 6.5

function Phan6

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

>> Phan6

Bai 6.5

Nhap vao gia tri cua T: 8

Gia tri cua h(T)la:

ans =

-2

VII VÒNG LẶP

Bài 7.6

function Phan7

disp('Bai 7.6')

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

>> Phan7

Bai 7.6

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

Ket qua:

a(ij)=

24 8 28

6 2 7

36 12 42

VIII 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

>> Phan8

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

IX Bài tập Symbolic math

Bài 9.5

function Phan9

disp('Bai 9.5')

disp('Hinh cau tam O ban kinh bang 3:')

syms ab;

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

>> Phan9

Bai 9.5

Hinh cau tam O ban kinh bang 3: