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Ụ MÔ HÌNH HOÁ BẰNG SỐ GIÁO VIÊN HƯỚNG DẪN : DƯƠNG THANH QUỲNH SINH VIÊN THỰC HIỆN : NGUYỄN ĐỨC LONG
Trang 1TRƯỜ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Ụ MÔ HÌNH
HOÁ BẰNG SỐ
GIÁO VIÊN HƯỚNG DẪN : DƯƠNG THANH
QUỲNH
SINH VIÊN THỰC HIỆN : NGUYỄN ĐỨC LONG
MÃ SỐ SINH VIÊN : 1528563
Trang 2Bài 1.1:Làm quen Matlab
GI IẢ
disp( ' _' ) disp( ' LAM QUEN VOI MATLAB ' ) disp( ' -' ) disp( 'Bai1.1' )
disp( 'phan a' )
2/2*3
disp( 'phan b' )
6-2/5+7^2-1
disp( 'phan c' )
10/2\5-3+2*4
disp( 'phan d' )
3^2/4
disp( 'phan e' )
3^2^2
disp( 'phan f' )
2+round(6/9+3*2)/2-3
disp( 'phan g' )
2+floor(6/9+3*2)/2-3
disp( 'phan h' )
2+ceil(6/9+3*2)/2-3
disp( 'phan k' )
fix(4/9)+fix(3*(5/6))
disp( ' ket thuc bai 1 ' )
end
LAM QUEN VOI MATLAB
-Bai1.1
Trang 3phan a ans = 3 phan b ans = 53.6000 phan c ans = 6 phan d ans = 2.2500 phan e ans = 81 phan f ans =
Trang 42.5000
phan g
ans =
2
phan h
ans =
2.5000
phan k
ans =
2
ket thuc bai 1
Bài 2: Vector
function Phan2
disp( '====================================' )
disp( ' 2.Vector ' )
disp( '====================================' )
disp( 'Bai2.4' )
x=[3;2;6;8]
y=[4;1;3;5]
disp( 'Cau c: Chia cac phan tu cua y tuong ung voi cac phan tu cua x' ) disp( 'Ket qua:' )
disp(y./x)
end
Đáp án:
>> Phan2
====================================
2.Vector
Trang 5Bai2.4
x =
3
2
6
8
y =
4
1
3
5
Cau c: Chia cac phan tu cua y tuong ung voi cac phan tu cua x
Ket qua:
1.3333
0.5000
0.5000
0.6250
Bài 3.3: Ma trận
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 la dong thu nhat cua A')
x=A(1,1:3)
disp('Gan cho ma tran y la hai dong con lai (cuoi) cua A') y=A(2:3,:)
disp('ket thuc bai 3')
ĐÁP ÁN
>> Bai33
Bai3.3
Cho ma tran A = [2 4 1 ; 6 7 2 ; 3 5 9]
A =
2 4 1
Trang 66 7 2
3 5 9
Gan cho vector x la dong thu nhat cua A
x =
2 4 1
Gan cho ma tran y la hai dong con lai (cuoi) cua A
y =
6 7 2
3 5 9
ket thuc bai 3
Bài 4.1 Vẽ đồ thị hàm số f(x)=sin(1/x) với 0.01<x<0.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))
disp('ket thuc bai 4')
ĐÁP ÁN
>> Bai41
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:
Trang 7ket thuc bai 4
Bài 5.6 Biểu thức lôgic
function Phan5
disp( '====================================' )
disp( ' 5.Bieu thuc logic ' )
disp( '====================================' )
disp( 'Bai 5.6: Cho x = [3 15 9 12 -1 0 -12 9 6 1]' )
disp( 'Cau b: Viet lenh thuc hien chuyen cac boi so cua so 3 thanh so 3' ) x=[3 15 9 12 -1 0 -12 9 6 1];
for (i=1:length(x))
if x(i)~=0
if mod(x(i),3)==0
x(i)=3;
end
end
end
Trang 8disp( 'x sau khi chuyen la:' )
disp(x)
end
Đáp án:
>> Phan5
====================================
5.Bieu thuc logic
====================================
Bai 5.6: Cho x = [3 15 9 12 -1 0 -12 9 6 1]
Cau b: Viet lenh thuc hien chuyen cac boi so cua so 3 thanh so 3
x sau khi chuyen la:
3 3 3 3 -1 0 3 3 3 1
Bài 6.5 Câu lệnh If
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
>> Bai65
Nhap vao gia tri cua T: 8
Gia tri cua h(T)la:
ans =
-2
Bài 7.8 Vòng lặp
function Bai78
disp('Bai7.8')
disp(' Co bao nhieu phan tu dung truoc 1 phan tu co gia tri 0.8 den 0.85')
x=rand(1,10);
n=1;
while (x(n)<0.8)|(x(n)>0.85)
n=n+1;
if n>10
break
end
Trang 9if n>10
disp('Vo nghiem')
else
ketqua=n-1
end
disp(' an tiep de lam phan b')
pause
disp('co bao nhieu phan tu nam trong khoang [0.01, 0.5] dung truoc gia')
disp(' tri trung binh cua vector to ')
x=rand(1,10)
n=1;
x
TrungBinh=mean(x)
while x(n)<=mean(x)
n=n+1;
end
strcat('so thu : ',num2str(n),' la so dau tien co gia tri lon hon trung binh')
if n==1
disp('Ko co so nao thuoc khoang [0.1 0.4] truoc gia tri dau tien lon hon trung binh')
else
y= x((x(1:n-1)>=0.1)&(x(1:n-1)<=0.4));
if isempty(y)==1
disp('Ko co so nao thuoc khoang [0.1 0.4] truoc gia tri dau tien lon hon trung binh')
else
strcat('co :',num2str(length(y)),' so thuoc khoang [0.1 0.4] truoc gia tri dau tien lon hon trung binh')
end
end
ĐÁP ÁN
>> Bai78
Bai7.8
Co bao nhieu phan tu dung truoc 1 phan tu co gia tri 0.8 den 0.85
Vo nghiem
an tiep de lam phan b
co bao nhieu phan tu nam trong khoang [0.01, 0.5] dung truoc gia tri trung binh cua vector to
Trang 10x =
0.2760 0.6797 0.6551 0.1626 0.1190 0.4984 0.9597 0.3404 0.5853 0.2238
x =
0.2760 0.6797 0.6551 0.1626 0.1190 0.4984 0.9597 0.3404 0.5853 0.2238
TrungBinh =
0.4500
ans =
'so thu :2 la so dau tien co gia tri lon hon trung binh' ans =
'co :1 so thuoc khoang [0.1 0.4] truoc gia tri dau tien lon hon trung binh'
Bài 8.1 Lập trình
function Phan8
disp( '===================================' )
disp( ' 8.Lap trinh ' )
disp( '===================================' )
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<x)&&(x<=200)
td=(550*100+(x-100)*1200)*1.1;
elseif (200<x)&&(x<=300)
td=(550*100+1200*100+(x-200)*1800)*1.1;
elseif x>300
td=(550*100+1200*100+1800*100+(x-300)*2000)*1.1;
Trang 11end
Đáp án:
> Phan8
===================================
8.Lap trinh
===================================
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
Bài 9.5 Bài tập Symbolic Math
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
ĐÁP ÁN
>> Bai95
Bai 9.5
Hinh cau tam O ban kinh bang 3:
Trang 12
-Bài tập lớn của em đến đây là hết rồi ạ Em chúc thầy luôn mạnh
khỏe, có nhiều niềm vui trong công việc và cuộc sống ạ