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 (246.69 KB, 13 trang )
TRƯỜNG ĐẠI HỌC XÂY DỰNG
Ban Đào Tạo Và Quản Lý 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áo viên hướng dẫn : Dương Thanh Quỳnh
Sinh viên thực hiện : Bùi Thị Thùy Dương
Mã số sinh viên
: 1509063
Lớp quản lý
: 63CLC1
1.Làm quen Matlab
1.1 Tính bằng tay các biểu thức sau, rồi thử lại bằng matlab
function phan1
disp ('PHAN I: LAM QUEN VOI MATLAB')
disp ('Bai 1.1: Tinh roi thu lai voi Matlab')
disp ('Cau a')
2/2*3
disp ('Cau b')
6-2/5+7^2-1
disp ('Cau c')
10/2\5-3+2*4
disp ('Cau d')
3^2/4
disp ('Cau e')
3^2^2
disp ('Cau f')
2+round(6/9+3*2)/2-3
disp ('Cau g')
2+floor(6/9+3*2)/2-3
disp ('Cau h')
2+ceil(6/9+3*2)/2-3
disp ('Cau k')
fix(4/9)+fix(3*(5/6))
ĐÁP ÁN
PHAN I: LAM QUEN VOI MATLAB
Bai 1.1: Tinh roi thu lai voi Matlab
Cau a
ans =
3
Cau b
ans =
53.6000
Cau c
ans =
6
Cau d
ans =
2.2500
Cau e
ans =
81
Cau f
ans =
2.5000
Cau g
ans =
2
Cau h
ans =
2.5000
Cau k
ans =
2
2.
VECTOR
Bài 2.2
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 a')
x(3)
disp('phan b')
x(1:7)
disp('phan c')
x(1:end)
disp('phan d')
x(1:end-1)
disp('phan e')
x(6:-2:1)
disp('phan f')
x([1 6 2 1 1])
disp('phan g')
sum(x)
ĐÁP ÁN
Bai2.2
cho x=[ 3 1 5 7 9 2 6]
x=
3
1
5
7
9
2
6
1
5
7
9
2
6
phan a
ans =
5
phan b
ans =
3
phan c
ans =
3
1
5
7
9
2
1
5
7
9
2
7
1
2
1
3
3
6
phan d
ans =
3
phan e
ans =
2
phan f
ans =
3
phan g
ans =
33
3. Ma trận
Bài 3.3
function Bai33
disp('Bai3.3')
disp('Cho ma tran A = [2 4 1 ; 6 7 2 ; 3 5 9]')
A = [2 4 1 ; 6 7 2 ; 3 5 9]
disp(' Gan cho vector x là dong thu nhat cua A')
x=A(1,1:3)
disp('Gan cho ma tran y là hai dong con lai (cuoi) cua A')
y=A(2:3,:)
disp('Tinh tong theo dong ma tran A. ')
sum(A)
disp('Tinh tong theo cot ma tran A. ')
sum(A')
disp('gia tri lon nhat')
max(max(A))
disp('gia tri nho nhat')
min(min(A))
disp('Tính tong các phan tu cua A ')
sum(sum(A))
ĐÁP ÁN
Bai3.3
Cho ma tran A = [2 4 1 ; 6 7 2 ; 3 5 9]
A=
2
6
3
4
7
5
1
2
9
Gan cho vector x là dong thu nhat cua A
x=
2
4
1
Gan cho ma tran y là hai dong con lai (cuoi) cua A
y=
6
3
7
5
2
9
Tinh tong theo dong ma tran A.
ans =
11
16
12
Tinh tong theo cot ma tran A.
ans =
7
15
17
gia tri lon nhat
ans =
9
gia tri nho nhat
ans =
1
Tính tong các phan tu cua A
ans =
39
4. Vẽ đồ thị
Bài 4.1
function Bai41
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
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:
5. Biểu thức logic
Bài 5.5
function Bai55
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(' nhung phan tu duong cua ma tran:')
X(X>0)
disp(' nhung phan tu am cua ma tran:')
X(X<0)
disp(' tong cac phan tu le cua ma tran:')
D(1:16)=X(1:4,1:4);
G=-100:2:100;
K=setdiff(D,G);
sum(K)
disp('gia tri lon nhat cac phan tu chan cua ma tran:')
J=-101:2:101;
max(setdiff(D,J))
disp('gan nhung phan tu 0 thanh 1:')
X(X==0)=1
ĐÁP ÁN
Bai5.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:
ans =
1
5
1
6
4
8
3
9
25
7
2
nhung phan tu am cua ma tran:
ans =
-2
-11
-3
tong cac phan tu le cua ma tran:
ans =
36
gia tri lon nhat cac phan tu chan cua ma tran:
ans =
8
gan nhung phan tu 0 thanh 1:
X=
1 -2 8 9
5 6 1 25
1 -11 3 7
1 4 -3 2
6.
Câu lệnh if
Bài 6.5
function Bai65
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')
ĐÁP ÁN
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
7.
Vòng lặp
Bài 7.6
function Bai 7.6
disp('Bai 7.6 : Cho x=[4 1 6] va y=[6 2 7], tinh cac
mang/vecto/ma tran sau')
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)
ĐÁP ÁN
Bai 7.6 : Cho x=[4 1 6] va y=[6 2 7], tinh cac mang/vecto/ma tran sau
Cau a: a(ij)=x(i)y(j)
Ket qua:
a(ij)=
24 8 28
6 2 7
36 12 42
8. Lập trình
Bài 8.1
function Phan8
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
elseif (200
elseif x>300
td=(550*100+1200*100+1800*100+(x-300)*2000)*1.1;
end
end
ĐÁP ÁN
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
9.
Bài tập Symbolic math
Bài 9.5
function Phan9
disp('Bai 9.5: 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
Bai 9.5: Hinh cau tam O ban kinh bang 3
