Bài tập lớn công cụ mô hình hóa bằng số (1)
Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (276.91 KB, 18 trang )
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Ý
:
function bai1_3
DƯƠNG THANH QUỲNH
PHẠM MẠNH HÙNG
93363
63CLC1
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
2
2
1
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
>> 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
d.Gan cho vector y gia tri lon hon 10 cua x
y=
1
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
elseif (y>20000)&(y<=50000)
disp('Khi 20000
elseif y>50000
disp('Khi 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
t =
900
>> bai6_7
ta nhap mot so y
y=25000
y =
25000
Khi 20000
1950
>> bai6_7
ta nhap mot so y
y=75000
y =
75000
Khi 50000
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)
end
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
2
7
y=
6
a=
24
8 28
6
2
7
36 12 42
phan b. bij=xi/yj
x=
4
1
6
2
7
y=
6
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
2
7
y=
6
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
2
7
y=
6
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
t=t+t*(10/100)
elseif (y>200)&(y<=300)
disp('Khi 200
t=t+t*(10/100)
elseif y>300
disp('Khi 300
t=t+t*(10/100)
end
>> bai8_1
ta nhap so dien y (KW)
y=120
Khi 100
t=
86900
>> bai8_1
ta nhap so dien y (KW)
y=320
Khi 300
t=
434500
>> bai8_1
ta nhap so dien y (KW)
y=480
Khi 300
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)/(sin(2*b)*sin(a + b) [ -cos(2*a)/(sin(2*b)*sin(a
+ b)/(sin(2*b)*sin(a + b) -
+ b) - cos(2*a)*sin(a - b)), -sin(a
cos(2*a)*sin(a - b))]
+ b) - cos(2*a)*sin(a - b)), sin(a
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
|
|
||||||||||||////////////////\\\\\\\\\\\\\\\||||||||||||||||||||
