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 QUẢN LÝ VÀ ĐÀO TẠO KỸ SƯ CHẤT LƯỢNG CAO BÀI TẬP LỚN CÔNG CỤ HÓA BẰNG MÔ HÌNH SỐ GIÁO VIÊN HƯỚNG DẪN: SINH VIÊN THỰC HIỆN: MÃ SỐ SINH VIÊN : LỚP QUẢN LÝ :

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

BAN QUẢN LÝ VÀ ĐÀO TẠO KỸ SƯ CHẤT LƯỢNG CAO

BÀI TẬP LỚN

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

GIÁO VIÊN HƯỚNG DẪN:

SINH VIÊN THỰC HIỆN:

MÃ SỐ SINH VIÊN :

LỚP QUẢN LÝ :

DƯƠNG THANH QUỲNH PHẠM MẠNH HÙNG

93363 63CLC1

function bai1_3

disp( '||||||||||||||\\\\\\\\\\\\\\\////////////////|||||||||||||||||||||' )

disp( '| |' )

disp( '| TRUONG DAI HOC XAY DUNG |' )

disp( '| |' )

disp( '| Ho va ten: PHẠM MẠNH HÙNG |' )

disp( '| Lop : 63CLC1 |' )

disp( '| MSSV : 93363 |' )

disp( '| GVHD : DUONG THANH QUYNH |' )

disp( '| Mon thi : Cong cu mo hinh hoa bang so |' )

disp( '| |' )

disp( '||||||||||||||////////////////\\\\\\\\\\\\\\\||||||||||| |||||||||' ) disp( ' ' )

disp( 'Bai 1.3 ' )

disp( 'Cho x=2, y=3,z=x' ) x=2; y=3; z=x disp( 'Ta co y moi: y=y+z' ) y=y+z disp( 'Ta co x moi: x=y+x-z' ) x=y+x-z disp( 'Ket qua phep tinh x+y-z luc nay bang' ) x+y-z end

Bai 1.3

Cho x=2, y=3,z=x z = 2

Ta co y moi: y=y+z y = 5

Ta co x moi: x=y+x-z

x =

5

Ket qua phep tinh x+y-z luc nay bang

ans =

8

function bai2_3

disp( 'Bai 2.3 ' ) disp( 'Cho x=[2 5 1 6]' )

x=[2 5 1 6]

disp( 'a.Cong them 16 vao tat ca cac phan tu' )

x+16

disp( 'b.cong them 3 vao cac phan tu vi tri le' )

x(1:2:end)=x(1:2:end)+3

disp( 'c.lay can bac 2 tung phan tu' )

sqrt(x)

disp( 'd.binh phuong tat ca cac phan tu' )

x.^2

end

>> bai2_3

Bai 2.3

Cho x=[2 5 1 6]

x =

2 5 1 6

a.Cong them 16 vao tat ca cac phan tu

ans =

18 21 17 22

b.cong them 3 vao cac phan tu vi tri le

x =

5 5 4 6

c.lay can bac 2 tung phan tu

ans =

2.2361 2.2361 2.0000 2.4495

d.binh phuong tat ca cac phan tu

ans =

25 25 16 36

function bai3_7

disp( 'Bai 3.7 ' ) disp( 'a.Tu he phuong trinh, ta co ma tran A va B:' )

A=[2 1 5 1; 1 1 -3 -4; 3 6 -2 1; 2 2 2 -3]

B=[5 -1 8 2]'

disp( 'ta tim duoc X1=A^(-1)*B' )

X1=inv(A)*B

disp( 'b.Tu he phuong trinh, ta co ma tran C va D:' )

C=[1 1 1 1;1 2 3 4; 2 3 5 7; 1 1 2 7]

D=[2 2 2 2]'

disp( 'ta tim duoc X2=C^(-1)*D' )

X2=inv(C)*D

end

>> bai3_7

Bai 3.7

a.Tu he phuong trinh, ta co ma tran A va B:

A =

2 1 5 1

1 1 -3 -4

3 6 -2 1

2 2 2 -3

B =

5

-1

8

2

ta tim duoc X1=A^(-1)*B

X1 =

2.0000

0.2000

0

0.8000

b.Tu he phuong trinh, ta co ma tran C va D:

C =

1 1 1 1

1 2 3 4

2 3 5 7

1 1 2 7

D =

2

2

2

2

ta tim duoc X2=C^(-1)*D

X2 =

0

4.5000

-3.0000

0.5000

function bai4_8

x1=linspace(0,20,100); x2=linspace(20,25,10) y1=x1.*sin(x1);

y2=50*x2;

figure

plotyy(x1,y1,x2,y2)

end

function bai5_6

x=[3 15 9 12 -1 0 -12 9 6 1]

disp( 'a.Chuyen cac gia tri duong thanh gia tri am' )

x(x>0)=0

disp( 'b.Chuyen cac boi so cua 3 thanh so 3' )

x(mod(x,3)==0)=3

clear

disp( 'c.Nhan cac gia tri chan cho 5' )

x=[3 15 9 12 -1 0 -12 9 6 1];

x(mod(x,2)==0)=x(mod(x,2)==0)*5

disp( 'd.Gan cho vector y gia tri lon hon 10 cua x' )

clear

x=[3 15 9 12 -1 0 -12 9 6 1];

y=x(x>10)

disp( 'e.Chuyen cac gia tri nho hon trung binh cong thanh gia tri 0' )

clear

x=[3 15 9 12 -1 0 -12 9 6 1]; x(x<mean(x))=0

>> bai5_6

x =

3 15 9 12 -1 0 -12 9 6 1

a.Chuyen cac gia tri duong thanh gia tri am

x =

0 0 0 0 -1 0 -12 0 0 0

b.Chuyen cac boi so cua 3 thanh so 3

x =

3 3 3 3 -1 3 3 3 3 3

c.Nhan cac gia tri chan cho 5

x =

3 15 9 60 -1 0 -60 9 30 1

d.Gan cho vector y gia tri lon hon 10 cua x

y =

15 12

e.Chuyen cac gia tri nho hon trung binh cong thanh gia tri 0

x =

0 15 9 12 0 0 0 9 6 0

function bai6_7

disp( 'ta nhap mot so y' )

y=input( 'y=' )

if y<=10000

disp( 'Khi y<=10000, ta co ket qua t(y) sau' )

t=200

elseif (y>10000)&(y<=20000)

disp( 'Khi 10000<y<=20000, ta co ket qua t(y) sau' ) t=200 + 0.1*(y - 10000)

elseif (y>20000)&(y<=50000)

disp( 'Khi 20000<y<=50000, ta co ket qua t(y) sau' ) t=1200 + 0.15*(y - 20000)

elseif y>50000

disp( 'Khi 50000<y, ta co ket qua t(y) sau' )

t=5700 + 0.25*(y - 50000)

end

>> bai6_7

ta nhap mot so y

y=5000

y =

5000

Khi y<=10000, ta co ket qua t(y) sau

t =

200

>> bai6_7

ta nhap mot so y

y=17000

y =

17000

Khi 10000<y<=20000, ta co ket qua t(y) sau

t =

900

>> bai6_7

ta nhap mot so y

y=25000

y =

25000

Khi 20000<y<=50000, ta co ket qua t(y) sau

t =

1950

>> bai6_7

ta nhap mot so y

y=75000

y =

75000

Khi 50000<y, ta co ket qua t(y) sau

t =

11950

function bai7_6

disp( 'phan a aij=xi*yj' )

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

a

disp( 'phan b bij=xi/yj' )

clear

x = [4 1 6] , y = [6 2 7]

for i=1:3

for j=1:3

b(i,j)=x(i)/y(j);

end

end

b

disp( 'phan c ci=xi*yj' )

clear

x = [4 1 6] , y = [6 2 7]

for i=1:3

c(i)=x(i)*y(i)

disp( 'Tinh tong cac phan tu cua c' ) sum(c)

disp( 'phan d dij=xi*(2+xi+yj)' ) clear

x = [4 1 6] , y = [6 2 7]

for i=1:3

for j=1:3

d(i,j)=x(i)/(2+x(i)+y(j)); end

end

d

>> bai7_6

phan a aij=xi*yj

x =

4 1 6

y =

6 2 7

a =

24 8 28

6 2 7

36 12 42

phan b bij=xi/yj

x =

4 1 6

y =

6 2 7

b =

0.6667 2.0000 0.5714 0.1667 0.5000 0.1429 1.0000 3.0000 0.8571

phan c ci=xi*yj

x =

4 1 6

y =

6 2 7

c =

24

c =

24 2

c =

24 2 42

Tinh tong cac phan tu cua c

ans =

68

phan d dij=xi*(2+xi+yj)

x =

4 1 6

y =

6 2 7

d =

0.3333 0.5000 0.3077

0.1111 0.2000 0.1000

0.4286 0.6000 0.4000

function bai8_1

disp( 'ta nhap so dien y (KW)' )

y=input( 'y=' );

if y<=100

disp( 'Khi y<=10000, ta co ket qua tien dien t(y) sau' ) t=550*y;

t=t+t*(10/100)

elseif (y>100)&(y<=200)

disp( 'Khi 100<y<=200, ta co ket qua tien dien t(y) sau' ) t=550*100 + 1200*(y - 100);

t=t+t*(10/100)

elseif (y>200)&(y<=300)

disp( 'Khi 200<y<=300, ta co ket qua tien dien t(y) sau' ) t= 550*100 + 1200*100 + 1800*(y-200);

t=t+t*(10/100)

elseif y>300

disp( 'Khi 300<y, ta co ket qua tien dien t(y) sau' )

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

t=t+t*(10/100)

end

>> bai8_1

ta nhap so dien y (KW)

y=120

Khi 100<y<=200, ta co ket qua tien dien t(y) sau

t =

86900

>> bai8_1

ta nhap so dien y (KW)

y=320

Khi 300<y, ta co ket qua tien dien t(y) sau

t =

434500

>> bai8_1

ta nhap so dien y (KW)

y=480

Khi 300<y, ta co ket qua tien dien t(y) sau

t =

786500

function bai9_2

syms a b u v A B

A=[sin(a+b) sin(a-b); cos(2*a) sin(2*b)]

B=[u v; v-v u+v]

disp( 'a.Tinh ma tran C=A*B, H=A+B, E=A-B' )

C=A*B

H=A+B

E=A-B

clear A B

A=[sin(a+b) sin(a-b); cos(2*a) sin(2*b)];

B=[u v; v-v u+v];

disp( 'b Tim ma tran nghich dao, chuyen vi cua A va B' ) chuyenviA=A'

chuyenviB=B'

nghichdaoA=inv(A)

nghichdaoB=inv(B)

disp( 'c Khai trien ma tran C' )

C=A*B;

D=det(C)

disp( 'Tính dao ham cua D theo a' )

diff(D,a)

disp( 'Tinh dao ham cua D theo v' )

diff(D,v)

disp( 'Tinh nguyen ham cua D theo a' )

int(D,a)

disp( 'Tinh nguyen ham cua D theo v' )

int(D,v)

disp( 'e Thay so' )

disp( ' Tính A voi a=0, b=pi/2' )

A=subs(A,{a,b},{0,pi/2})

disp( 'Tinh B voi u=3, v=6' )

B=subs(B,{u,v},{3,6})

disp( '||||||||||||\\\\\\\\\\\\\\\////////////////||||||||||||||||||||' ) disp( '| |' ) disp( '| SEE YOU AGAIN |' ) disp( '| |' ) disp( '||||||||||||////////////////\\\\\\\\\\\\\\\||||||||||||||||||||' )

>> bai9_2

A =

[ sin(a + b), sin(a - b)]

[ cos(2*a), sin(2*b)]

B =

[ u, v]

[ 0, u + v]

a.Tinh ma tran C=A*B, H=A+B, E=A-B

C =

[ u*sin(a + b), v*sin(a + b) + sin(a - b)*(u + v)]

[ u*cos(2*a), sin(2*b)*(u + v) + v*cos(2*a)]

H =

[ u + sin(a + b), v + sin(a - b)]

[ cos(2*a), u + v + sin(2*b)]

E =

[ sin(a + b) - u, sin(a - b) - v]

[ cos(2*a), sin(2*b) - v - u]

b Tim ma tran nghich dao, chuyen vi cua A va B

chuyenviA =

[ sin(conj(a) + conj(b)), cos(2*conj(a))]

[ sin(conj(a) - conj(b)), sin(2*conj(b))]

chuyenviB =

[ conj(u), 0]

[ conj(v), conj(u) + conj(v)]

nghichdaoA =

[ sin(2*b)/(sin(2*b)*sin(a + b) - cos(2*a)*sin(a - b)), -sin(a

- b)/(sin(2*b)*sin(a + b) - cos(2*a)*sin(a - b))]

[ -cos(2*a)/(sin(2*b)*sin(a + b) - cos(2*a)*sin(a - b)), sin(a + b)/(sin(2*b)*sin(a + b) - cos(2*a)*sin(a - b))]

nghichdaoB =

[ 1/u, -v/(u*(u + v))]

[ 0, 1/(u + v)]

c Khai trien ma tran C

D =

u^2*sin(2*b)*sin(a + b) - u^2*cos(2*a)*sin(a - b) +

u*v*sin(2*b)*sin(a + b) - u*v*cos(2*a)*sin(a - b)

Tính dao ham cua D theo a

ans =

2*u^2*sin(2*a)*sin(a - b) - u^2*cos(2*a)*cos(a - b) +

u^2*sin(2*b)*cos(a + b) + u*v*sin(2*b)*cos(a + b) -

u*v*cos(2*a)*cos(a - b) + 2*u*v*sin(2*a)*sin(a - b)

Tinh dao ham cua D theo v

ans =

u*sin(2*b)*sin(a + b) - u*cos(2*a)*sin(a - b)

Tinh nguyen ham cua D theo a

ans =

(u*(u + v)*(cos(3*a - b)/2 - (3*cos(a + b))/2 + (3*sin(a - b))/2

- (3*sin(a + 3*b))/2))/3

Tinh nguyen ham cua D theo v

ans =

((u*sin(2*b)*sin(a + b))/2 - (u*cos(2*a)*sin(a - b))/2)*v^2 + (u^2*sin(2*b)*sin(a + b) - u^2*cos(2*a)*sin(a - b))*v

e Thay so

Tính A voi a=0, b=pi/2

A =

[ 1, -1]

[ 1, 0]

Tinh B voi u=3, v=6

B =

[ 3, 6]

[ 0, 9]

||||||||||||\\\\\\\\\\\\\\\////////////////||||||||||||||||||||

| |

| SEE YOU AGAIN |

| |

||||||||||||////////////////\\\\\\\\\\\\\\\||||||||||||||||||||