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

BÀI TẬP LỚN MÔN

CÔNG CỤ MÔ HÌNH HÓA BẰNG SỐ

GIẢNG VIÊN: DƯƠNG THANH QUỲNH

Họ Và Tên : KHÚC VĂN TÔ

Lớp: 63CLC1 MSSV : 206263

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

function Bai11

disp('Bai1.1')

disp('Phan b')

6-2/5+7^2-1

disp('ket thuc bai 1')

Đáp Án

>> Bai11

Bai1.1

Phan b

ans =

53.6000

ket thuc bai 1

II Vector

Bài 2.2: Cho x= [3 1 5 7 9 2 6] dự đoán kết quả các dòng lệnh sau và thử lại bằng matlab

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

x(1:7)

disp('ket thuc bai 2')

ĐÁP ÁN

>> Bai22

Bai2.2

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

x =

3 1 5 7 9 2 6

phan b

ans =

3 1 5 7 9 2 6

ket thuc bai 2

Bài3: Ma Trận

function Bai35

disp( 'Bai 3.5:Cho vecto x=[2 4 1 6],y=[5 9 1 0]' )

disp( 'Cau c: Hay tao ra ma tran C co tinh chat: cot 1 va 3 co gia tri la vecto

x, cot 2 va 4 co gia tri la vecto y' )

x=[2 4 1 6];

y=[5 9 1 0];

M=[x;y;x;y];

disp( 'Ket qua:' )

disp(M.')

end

Đáp án:

>> bai3

Bai 3.5:Cho vecto x=[2 4 1 6],y=[5 9 1 0]

Cau c: Hay tao ra ma tran C co tinh chat: cot 1 va 3 co gia tri la vecto x, cot 2 va 4 co gia tri la vecto y

Ket qua:

2 5 2 5

4 9 4 9

1 1 1 1

6 0 6 0

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 bai55

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

b=X(X<0)

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

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

Bài 6.1: Câu lệnh IF

function bai61

disp( 'bai6.1' )

disp( ' Du doan ket qua va thu lai bang matlab' )

disp( ' tinh m:' )

n=input( ' nhap gia tri cua n: n=' )

if n>1

m=n+1

else

m=n-1

end

disp( ' an tiep phim bat ky de lam phan b' )

pause

n=input( ' nhap gia tri cua n: n=' )

if n>1

m=n+1

else

m=n-1

end

disp( 'an tiep phim bat ky de lam phan c' )

pause

n=input( ' nhap gia tri cua n: n=' )

if n>1

m=n+1

else

m=n-1

end

end

ĐÁP ÁN

>> bai61

bai6.1

Du doan ket qua va thu lai bang matlab tinh m:

nhap gia tri cua n: n=4

n =

4

m =

5

an tiep phim bat ky de lam phan b nhap gia tri cua n: n=

n =

[]

m =

[]

an tiep phim bat ky de lam phan c

nhap gia tri cua n: n=6

n =

6

m =

7

Bài 7.2 : Phần lập trình

Function bai72

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 2 phan tu la C=8.220635e-19

Bai 8 : Lập trình

function bai81

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:

> bai81

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

Bài 9:Bài tập SYMBOLIC MATH

function Bai95

disp( 'Bai 9.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

ĐÁP ÁN

>Bai95

Bai9.5 Hinh cau tam O ban kinh bang 3: