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

BÀI TẬP MÔN CÔNG CỤ MÔHÌNH HÓA BẰNG SỐ GIẢNG VIÊN:DƯƠNG THANH QUỲNH SINH VIÊN:NGHIÊM XUÂN HIẾU LỚP 63CLC1_71563 1.1:Tính bằng tay các biểu thức sau, rồi thử lại bằng Matlab: c... Những p

BÀI TẬP MÔN CÔNG CỤ MÔ

HÌNH HÓA BẰNG SỐ

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

SINH VIÊN:NGHIÊM XUÂN HIẾU

LỚP 63CLC1_71563

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

c 10 / 2 \ 5 – 3 + 2 * 4

function phan1

disp('Bai1.1')

disp('phan c')

10/2\5-3+2*4

end

ĐÁP ÁN

>> phan1

Bai1.1

phan c

ans =

6

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:

c x(1:end)

function phan1

disp('===================================================') disp('| VECTO |') disp('===================================================') disp('Bai 2.2')

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

x=[ 3 1 5 7 9 2 6 ]

disp('phan c')

x(1:end)

end

ĐÁP ÁN

phan1

=====================================================

| VECTO |

=====================================================

Bai 2.2

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

x =

3 1 5 7 9 2 6

phan c

ans =

3 1 5 7 9 2 6

Bài 3.2: cho A = [2 7 9 7 ; 3 1 5 6 ; 8 1 2 5], dự đoán kết quả, giải thích; rồi thử lại bằng Matlab:

function phan1

disp('=====================================')

disp('| MA TRAN |')

disp('=====================================')

disp('Bai3.2')

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

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

disp('phan a')

A'

End

ĐÁP ÁN

phan1

=====================================

| MA TRAN |

=====================================

Bai3.2

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

A =

2 7 9 7

3 1 5 6

8 1 2 5

phan a

ans =

2 3 8

7 1 1

9 5 2

7 6 5

Bài 4.3: Vẽ đồ thị hàm số f(x) = sin(1/x) với 0.01 < x

< 0.1

function phan1

disp('==========================================')

disp('| VE DO THI |')

disp('==========================================')

disp('Bai43')

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 thhi ham so f(x) la :')

plot(x,sin(1./x))

end

ĐÁP ÁN

Bài 5.5: Cho ma trận X=[1 -2 8 9 ; 5 6 0 25 ; 1 -11 3 7 ; 0 4 -3 2] , tìm:

a Những phần tử dương của ma trận.

b Những phần tử âm của ma trận.

c Tổng các phần tử lẻ của ma trận.

disp('========================================') disp('| BIEU THUC LOGIC |') disp('=========================================') 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('phan a')

disp('nhung phan tu duong cua ma tran la : ')

X(X>0)

disp('phan b')

disp('nhung ohan tu am cua ma tran la : ')

X(X<0)

disp('phan c')

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

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

K=setdiff(D,G);

sum(K)

disp('phan d')

end

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

max(setdiff(D,J))

disp('phan e')

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

X(X==0)=1

ĐÁP ÁN

========================================

| BIEU THUC LOGIC |

=========================================

Bai5.5

cho ma tran X

X =

1 -2 8 9

5 6 0 25

1 -11 3 7

0 4 -3 2

phan a

nhung phan tu duong cua ma tran la :

ans =

1

5

1

6

4

8

3

9

25

7

2

phan b

nhung ohan tu am cua ma tran la :

ans =

-2

-11

-3

phan c

tong cac phan tu le cua ma tran la :

ans =

36

phan d

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

ans =

8

phan e

gan nhung phan tu 0 thanh 1

X =

1 -2 8 9

5 6 1 25

1 -11 3 7

1 4 -3 2

Bài 6.1: Dự đoán kết quả và thử lại bằng Matlab:

if n > 1

m = n+1

else

m = n – 1

end

a n = 7 m =

?

b n = 0 m =

?

c n = -10 m

= ? function phan1

disp('==========================================')

disp('| CAU LENH IF |') disp('==========================================') disp('Bai6.1')

disp(' du doan ket qua va thu lai bang matlab') disp('tinh m:')

disp('phan a')

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

if n>1

m=n+1

else

m=n-1

end

disp('an phim bat ki de lam phan c')

pause

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

if n>1

m=n+1

else

m=n-1

end

ĐÁP ÁN

phan1

==========================================

| CAU LENH IF |

==========================================

Bai6.1

du doan ket qua va thu lai bang matlab

tinh m:

phan a

nhap gia tri cua n: n=3

n =

3

m =

4

an phim bat ki de lam phan c

nhap gia tri cua n: n=5

n =

5

m =

6

7.6 Cho x = [4 1 6] và y = [6 2 7], tính các mảng/vector/ma trận sau

a aij = xiyj

function phan1

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

phan1

Bai7.6

cau a; a(ij)=x(i)y(j)

ket qua

a(ij))

24 8 28

6 2 7

36 12 42

Bài 8.5: Dãy số Fibonacci có mối quan hệ sau: Fn = Fn-1 + Fn-2, với F0 = F1 = 1.

a Tính 10 số Fibonacci đầu tiên.

b Trong 50 số Fibonacci đầu tiên, tính tỷ lệ Fn / Fn-1 N hận xét gì các giá

tỷ lệ này so với “tỷ lệ

vàng” ( (1 + sqrt(5))/2 ).

disp(' VONG LAP ') disp('===================================================')

n=input('nhap do lon:');

s=[]

if n==1|n==2

s=1;

else

s=[1 1]

for i=3:n

s= [s (i-1)+s(i-2)];

end

end

disp(' 10 so Fibonacci dau tien co gia tri lan luot la:') disp('an phim bat ki de lam tiep phan b')

pause

m=input('nhap do lon:');

a=[];

if m==1|m==2

a=1;

else

a=[1 1]

for i=3:m

a= [a (a(i-1)+a(i-2))];

end

end

disp(' 50 so Fibonacci dau tien co gia tri lan luot la :') a

disp('ti so giua Fn/Fn-1 cua day la :')

for i=2:m

A=a(i)/a(i-1)

end

A

disp('nhan xet:ta thay tu so A thu 13 tro di co gia tri bang

ty le vang')

ĐÁP ÁN

phan1

VONG LAP

========================================================

Bai8.5

nhap do lon:5

s =

[]

s =

1 1

10 so Fibonacci dau tien co gia tri lan luot la:

an phim bat ki de lam tiep phan b

nhap do lon:10

a =

1 1

50 so Fibonacci dau tien co gia tri lan luot la :

a =

1 1 2 3 5 8 13 21 34 55

ti so giua Fn/Fn-1 cua day la :

A =

1

A =

2

A =

1.5000

A =

1.6667

A =

1.6000

A =

1.6250

A =

1.6154

A =

1.6190

A =

1.6176

A =

1.6176

nhan xet:ta thay tu so A thu 13 tro di co gia tri bang ty le vang

Bài9.5: Vẽ hình cầu tâm O bán kính 3 bằng phương trình:

x=3*sin(a)*cos(b);

y=3*sin(a)*sin(b);

z=3*cos(a);

disp('====================================================='

)

disp('| Bai Tap symbolic math

|')

disp('====================================================='

)

disp('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;

ĐÁP ÁN