bài tập lớn matlab (3)

Trường Đại Học Xây Dựng

Ban 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ảng viên: DƯƠNG THANH QUỲNH Sinh Viên : NGUYỄN THỊ HỒNG Lớp : 58CLC1

MSSV : 345458

BÀI 1: LÀM QUEN VỚI MATLAB

function bai11

disp('bai tap 1.1');

2/2*3

6-2/5+7^2-1

10/2\5-3+2*4

3^2/4

3^2^2

2 + round(6/9+3*2)/2-3

2 + floor(6/9+3*2)/2-3

2 + ceil(6/9+3*2)/2-3

fix(4/9)+(fix(3*(5/6)))

end

bai tap 1.1

ans =

3

ans =

53.6000

ans =

6

ans =

2.2500

ans =

81

ans =

2.5000

ans =

2

ans =

2.5000

ans =

2

BÀI 2: VECTO

function BAI21

disp('Bai 2.1')

disp('Vector chua cac so nguyen tu 31 den 75 la :')

x1=[31:75]

disp('Vector co 10 phan tu co gia tri nguyen duoc lay ngau nhien trong khoang [0,100] la:')

x2=randi([0,100],[1,10])

disp('Vector co 10 phan tu co gia tri nguyen duoc lay ngau nhien trong khoang [-20,10] la:')

x3=randi([-20,10],[1,10])

end

Bai 2.1

Vector chua cac so nguyen tu 31 den 75 la :

x1 =

Columns 1 through 16

31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46

Columns 17 through 32

47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62

Columns 33 through 45

63 64 65 66 67 68 69 70 71 72 73 74 75

Vector co 10 phan tu co gia tri nguyen duoc lay ngau nhien trong khoang [0,100] la: x2 =

82 91 12 92 63 9 28 55 96 97

Vector co 10 phan tu co gia tri nguyen duoc lay ngau nhien trong khoang [-20,10] la: x3=

-16 10 9 -5 4 -16 -7 8 4 9

BÀI 3: MA TRẬN

function BAI33

disp('Bai 3.3 ');

disp('Cho ma tran A ');

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

disp(' CAU a ');

disp('Gan cho vecto x la dong thu nhat cua a');

disp('Ket qua cau a la:');

x=A(1,:)

disp('Nhan enter de tiep tuc');

disp(' ');

pause;

disp(' Cau b ');

disp('Gan cho ma tran y la 2 dong cuoi cua A');

disp('Ket qua cau b la:');

y=A(2:3,:)

disp('Nhan enter de tiep tuc');

disp(' ');

pause;

disp(' Cau c ');

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

disp('Tinh chuyen vi cua ma tran A');

B=A'

disp('Tinh tong theo cot cua ma tran chuyen vi');

sum(B)

disp('Tong theo dong ma tran A la tong theo cot cua B'); disp('Nhan enter de tiep tuc');

disp(' ');

pause;

disp(' Cau d ');

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

disp('Tong theo cac cot lan luot la:');

sum(A)

disp('Nhan enter de tiep tuc');

disp(' ');

pause;

disp(' Cau e ');

disp('Tim gia tri max va min cua ma tran A');

disp('Tim gia tri max cua A');

disp('Tim gia tri max cua moi cot la:');

B=max(A)

disp('Gia tri lon nhat cua matran B chinh la gia tri max cua A'); max(B)

disp('Tim gia tri min cua A');

disp('Tim gia tri min cua moi cot la:');

B=min(A)

disp('Gia tri nho nhat cua matran B chinh la gia tri min cua A'); min(B)

pause;

disp(' Cau f ');

disp('Tinh tong cac phan tu cua A');

disp('Tinh tong cac cot cua A la B=');

B=sum(A)

disp('Tinh tong cac phan tu cua A');

disp('Tong A la:');

sum(B)

disp('BAI 3.3 KET THUC');

pause;

end

Bai 3.3

Cho ma tran A

A =

2 4 1

6 7 2

3 5 9

CAU a

Gan cho vecto x la dong thu nhat cua a

Ket qua cau a la:

x =

2 4 1

Nhan enter de tiep tuc

Cau b

Gan cho ma tran y la 2 dong cuoi cua A

Ket qua cau b la:

y =

6 7 2

3 5 9

Nhan enter de tiep tuc

Cau c

Tinh tong theo dong ma tran A

Tinh chuyen vi cua ma tran A

B =

2 6 3

4 7 5

1 2 9

Tinh tong theo cot cua ma tran chuyen vi

ans =

7 15 17

Tong theo dong ma tran A la tong theo cot cua B Nhan enter de tiep tuc

Cau d

Tinh tong theo cot ma tran A

Tong theo cac cot lan luot la:

ans =

11 16 12

Nhan enter de tiep tuc

Cau e

Tim gia tri max va min cua ma tran A

Tim gia tri max cua A

Tim gia tri max cua moi cot la:

B =

6 7 9

Gia tri lon nhat cua matran B chinh la gia tri max cua A ans =

9

Tim gia tri min cua A

Tim gia tri min cua moi cot la:

B =

2 4 1

Gia tri nho nhat cua matran B chinh la gia tri min cua A

ans =

1

Cau f

Tinh tong cac phan tu cua A

Tinh tong cac cot cua A la B=

B =

11 16 12

Tinh tong cac phan tu cua A

Tong A la:

ans =

39

BAI 3.3 KET THUC

BÀI 4: VẼ ĐỒ THỊ

Bài 4.1

function BAI41

disp(' Do thi cua cac ham so x^3, e^x va e^(x^2) la ') x=0:0.1:4;

a=x.^3;

b=exp(x);

c=exp(x.^2);

subplot(2,2,1)

plot(x,a)

xlabel('0<=x<=4')

ylabel('a(x)')

subplot(2,2,2)

plot(x,b)

xlabel('0<=x<=4')

ylabel('b(x)')

subplot(2,2,3)

plot(x,c)

xlabel('0<=x<=4')

ylabel('c(x)')

end

Do thi cua cac ham so x^3, e^x va e^(x^2) la

Bài 4.18

function BAI418

disp('Ve bieu do quat cua mot bao cao thi truong'); disp('So nguoi uong coca la 850 nguoi');

disp('So nguoi uong pepsi la 720 nguoi');

disp('So nguoi uong nuoc tinh khiet la 600 nguoi'); disp('So nguoi uong do uong khac la 320 nguoi');

disp('Ta co do thi:');

Y=[850 720 600 320];

pie(Y,{' Coca ',' Pepsi ',' Nuoc tinh khiet ',' Khac '}) end

Ve bieu do quat cua mot bao cao thi truong

So nguoi uong coca la 850 nguoi

So nguoi uong pepsi la 720 nguoi

So nguoi uong nuoc tinh khiet la 600 nguoi

So nguoi uong do uong khac la 320 nguoi

Ta co do thi:

BÀI 5: BIỂU THỨC LOGIC

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

end

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: CÂU LỆNH If

Bài 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

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

BÀI 7: VÒNG LẶP

Bài 7.2

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)

end

To hop chap r, Nhap r=4

Cua n phan tu, Nhap n=8

To hop chap 4 cua 8 phan tu la C=70

BÀI 8: LẬP TRÌNH

Bài 8.1

function BAI81

disp('Tinh tien dien cua cac ho gia dinh')

disp(' BAI LAM ')

i = 1;

while i == 1

k = input('Nhap so KW dien tieu thu trong thang: ');

if k > 0

if k <= 100

t = k*550;

elseif k <= 200

t = 550*100 + 1200*(k - 100);

elseif k < 300

t = 550*100 + 1200*100 + 1800*(k - 200);

else

t = 550*100 + 1200*100 + 1800*100 + 2000*(k - 300);

end

fprintf('So tien dien phai tra trong thang la: %d dong \n' ,t)

fprintf('So tien thue VAT phai tra trong thang la: %d dong \n' ,t*10/100) fprintf('Tong so tien dien phai tra trong thang (gom VAT) la: %d dong \n'

,t*110/100)

else disp('So KW phai la so duong')

i = input('Nhan enter de nhap tiep, nhan phim so bat ky roi an enter de thoat \n');

if i == 1

i = 0;

end

if isempty(i) == 1

i = 1;

end

end

end

Tinh tien dien cua cac ho gia dinh

BAI LAM

Nhap so KW dien tieu thu trong thang: 356

So tien dien phai tra trong thang la: 467000 dong

So tien thue VAT phai tra trong thang la: 46700 dong

Tong so tien dien phai tra trong thang (gom VAT) la: 513700 dong

Nhan enter de nhap tiep, nhan phim so bat ky roi an enter de thoat

BÀI 9: BÀI TẬP SYMBOLC MATH

Bài 9.5

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

Bai 9.5

Hinh cau tam O ban kinh bang 3: