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 (158.04 KB, 11 trang )
Bài tập
Cơng cụ mơ hình hóa bằng số
Giảng viên : DƯƠNG THANH QUỲNH
Họ và tên: Vũ Đức Khang
Lớp:63CLC1
MSSV:107063
1.Làm quen Matlab
function Bai1
disp('Bai1.1.a')
disp('2/2*3')
disp('Ket qua')
disp(2/2*3)
Đáp án
>> Bai1
Bai1.1.a
2/2*3
Ket qua
3
2.Vec to
function Bai2
disp('Bai 2.2.a')
disp('Cho x=[3 1 5 7 9 2 6]')
disp('Du doan dong lenh sau va thu lai bang Matlab')
disp('x(3)')
x=[3 1 5 7 9 2 6];
disp('Dap an')
x(3)
Đáp án
>> Bai2
Bai 2.2.a
Cho x=[3 1 5 7 9 2 6]
Du doan dong lenh sau va thu lai bang Matlab
x(3)
Dap an
ans =
5
3.Ma trận
function Bai3
disp('Bai 3.2.d Du doan ket qua giai thich va thu lai bang
Matlab ')
disp('Cho A=[2 7 9 7;3 1 5 6;8 1 2 5]')
disp('reshape(A,2,6)')
disp('Ket qua')
A=[2 7 9 7;3 1 5 6;8 1 2 5]
reshape(A,2,6)
Đáp án
>> Bai3
Bai 3.2.d Du doan ket qua giai thich va thu lai bang
Matlab
Cho A=[2 7 9 7;3 1 5 6;8 1 2 5]
reshape(A,2,6)
Ket qua
A =
2
7
9
7
3
1
5
6
8
1
2
5
ans =
2
8
1
9
2
6
3
7
1
5
7
5
4. Vẽ đồ thị
function Bai4
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))
disp('ket thuc bai 4')
Đáp án
>> Bai4
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:
Bài 5: Biểu thức logic
function Bai5
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)
Đáp án
>> Bai5
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
6.Câu lệnh If
function Bai6
disp('Bai 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
Đáp án
>> Bai6
Bai 6.5
Nhap vao gia tri cua T: 5
Gia tri cua h(T)la:
ans =
-5
7.Vòng lặp
function Bai7
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)
end
Đáp án
>> Bai7
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
function Bai8
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
>> Bai8
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
Đáp án
9.Bài tập Symbolic Math
function Bai9
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;
end
Đáp án
>> Bai9
Bai 9.5: Hinh cau tam O ban kinh bang 3
