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 : DƯƠNG THANH QUỲNH SINH VIÊN THỰC HIỆN : PHÙNG NGỌC HÀ MY
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Ụ HÓA BẰNG MÔ HÌNH
SỐ
GIÁO VIÊN HƯỚNG DẪN : DƯƠNG THANH
QUỲNH
SINH VIÊN THỰC HIỆN : PHÙNG NGỌC HÀ MY
MÃ SỐ SINH VIÊN : 1532463
Trang 2I LÀM QUEN VỚI MATLAB
Bài 1.1: Tính bằng tay rồi thử lại bằng Matlab
PHẦN LẬP TRÌNH
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))
PHẦN ĐÁP ÁN
>> Bai11
PHAN I: LAM QUEN VOI MATLAB
Bai 1.1: Tinh roi thu lai voi Matlab
Cau a
Trang 3ans = 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
Trang 4II VECTOR
Bài 2.2 : Cho ma trận
PHẦN LẬP TRÌNH
function Bai23
disp ('PHAN II: VECTOR')
disp ('Bai 2.3')
disp ('Nhap ma tran x')
x=[2 5 1 6]
disp ('Cong tat ca cac phan tu trong ma tran x voi 16')
x+16
disp ('Cong cac phan tu le cua ma tran voi 3')
y=[3 0 3 0]
x+y
disp ('Lay can binh phuong tat ca cac phan tu cua ma tran x') sqrt(x)
disp ('Binh phuong tat ca cac phan tu')
x.^2
PHẦN ĐÁP ÁN
>> Bai23
PHAN II: VECTOR
Bai 2.3
Nhap ma tran x
x =
2 5 1 6
Cong tat ca cac phan tu trong ma tran x voi 16
ans =
18 21 17 22
Cong cac phan tu le cua ma tran voi 3
y =
3 0 3 0
ans =
Trang 55 5 4 6
Lay can binh phuong tat ca cac phan tu cua ma tran x
ans =
1.4142 2.2361 1.0000 2.4495
Binh phuong tat ca cac phan tu
ans =
4 25 1 36
III MA TRẬN
Bài 3.1 :
a x + y là dòng lệnh hợp lệ, vì x và y đều là các ma trận 3x1
b x + A là dòng lệnh hợp lệ vì x và A đều là các ma trận có 3 cột
c x’ + y là dòng lệnh hợp lệ vì x’ là ma trận 1x3
d A – [x’ y’] là dòng lệnh không hợp lệ vì A là ma trận 3x2 mà [x’ y’]
là ma trận 2x3
e [x ; y’] là dòng lệnh không hợp lệ vì x là ma trận 3x1 mà y’ là ma trận 1x3 nên không thể ghép dòng
f [x ; y] là dòng lệnh hợp lệ vì đây là lệnh ghép 2 ma trận 3x1
g A – 3 là dòng lệnh hợp lệ vì đây là phép trừ tất cả các phần tử của
ma trận A cho 3
PHẦN LẬP TRÌNH
function Bai31
disp('Nhap ma tran x, y, A')
x=[1 4 8]
y=[2 1 5]
A=[3 1 6 ; 5 2 7]
disp ('Thu xem dong lenh nao la hop le')
x+y
x+A
x'+y
A-[x' y']
[x,y']
[x,y]
A-3
PHẦN ĐÁP ÁN
Trang 6>> Bai31
Nhap ma tran x, y, A
x =
1 4 8
y =
2 1 5
A =
3 1 6
5 2 7
Thu xem dong lenh nao la hop le ans =
3 5 13
ans =
4 5 14
6 6 15
ans =
3 2 6
6 5 9
10 9 13
Matrix dimensions must agree Error in Bai31 (line 10)
A-[x' y']
Error in Bai31 (line 10)
[x,y']
ans =
1 4 8 2 1 5
Trang 7ans =
0 -2 3
2 -1 4
IV VẼ ĐỒ THỊ
Bài 4.4: Vẽ lại đồ thị hàm số ex (0 < x < 4) nhưng với đường đứt
nét, màu đỏ và dấu (marker) là dấu +
PHẦN LẬP TRÌNH
function Bai44
disp ('Ve do thi ham so')
disp ('Khoang gia tri cua x la')
x=0:4
disp ('Do thi ham so can ve la')
plot(x,exp(x),' r+')
PHẦN ĐÁP ÁN
>> Bai44
Ve do thi ham so
khoang gia tri cua x la
x =
0 1 2 3 4
Do thi ham so can ve la
>> Bai44
Ve do thi ham so
khoang gia tri cua x la
x =
0 1 2 3 4
Do thi ham so can ve la
Trang 8V BIỂU THỨC LOGIC Bài 5.5:
PHẦN LẬP TRÌNH
function Bai55
disp( 'Bai 5.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 la:' )
a=X([X]>0)
disp( 'Nhung phan tu am cua ma tran la:' )
b=X(X<0)
disp( 'Tong cac phan tu le cua ma tran la:' )
S=sum(X(mod([X],2)==1))
disp( 'Gia tri lon nhat cua cac phan tu chan cua ma tran la:' ) Max=max(X(mod(X,2)==0))
disp( 'Nhan phim bat ki de tiep tuc' )
pause
e=X
disp( 'Ma tran X sau khi gan nhung phan tu bang 0 thanh 1 la:' ) e(X==0)=1
PHẦN ĐÁP ÁN
>> Bai55
Bai 5.5
Cho ma tran X:
X =
1 -2 8 9
5 6 0 25
Trang 91 -11 3 7
0 4 -3 2
Nhung phan tu duong cua ma tran la:
a =
1
5
1
6
4
8
3
9
25
7
2
Nhung phan tu am cua ma tran la:
b =
-2
-11
-3
Tong cac phan tu le cua ma tran la:
S =
37
Gia tri lon nhat cua cac phan tu chan cua ma tran la: Max =
8
Nhan phim bat ki de tiep tuc
e =
1 -2 8 9
5 6 0 25
1 -11 3 7
0 4 -3 2
Ma tran X sau khi gan nhung phan tu bang 0 thanh 1 la:
e =
1 -2 8 9
5 6 1 25
1 -11 3 7
1 4 -3 2
Trang 10VI CÂU LỆNH IF
Bài 6.5:
PHẦN LẬP TRÌNH
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
PHẦN ĐÁP ÁN
>> Bai65
Nhap vao gia tri cua T: 8
Gia tri cua h(T)la:
ans =
-2
VII VÒNG LẶP
Bài 7.2:
PHẦN LẬP TRÌNH r=input('to hop chap r, nhap r=');
n=input('cua n phan tu, nhap n=');
C=1;
for i=1:n
C=C*i;
end
for i=1:r
C=C/i;
end
for i=1:(nr)
C=C/i;
end;
fprintf('to hop chap %d cua %d phan tu la C=%d\n',r,n,C)
PHẦN ĐÁP ÁN
>> Bai72
to hop chap r, nhap r=20
cua n phan tu, nhap n=2
to hop chap 20 cua 2 phan tu la C=8.220635e-19
Trang 11VIII LẬP TRÌNH
Bài 8.1:
PHẦN LẬP TRÌNH
A=input('nhap vao so dien tieu thu cua cac ho gia dinh:'); T=zeros(1,length(A));
for i=1:length(A)
if A(i)<=100
T(i)=A(i)*550;
elseif A(i)<=200
T(i)=55000+(A(i)-100)*1200;
elseif A(i)<=300
T(i)=55000+ 120000+(A(i)-200)*1800;
else T(i)=55000+ 120000+180000+(A(i)-300)*2000;
end
T(i)=T(i)*1.1;
end
for i=1:length(A)
fprintf('tien dien cua ho gia dinh thu %d la:
%d\n',i,T(i))
end
PHẦN ĐÁP ÁN
>> Bai81
nhap vao so dien tieu thu cua cac ho gia dinh:273
tien dien cua ho gia dinh thu 1 la: 337040
IX BÀI TẬP SYMBOLIC MATH
Bài 9.5
PHẦN LẬP TRÌNH
syms t p
x=3*sin(t)*cos(p);
y=3*sin(t)*sin(p) ;
z=3*cos(t) ;
ezsurf(x,y,z)
PHẦN ĐÁP ÁN