Hinh cau tam O ban kinh bang 3:.
Trang 1TRƯỜNG ĐẠI HỌC XÂY DỰNG
BAN KỸ SƯ CHẤT LƯỢNG CAO
BÀI TẬP MÔN CÔNG CỤ MÔ HÌNH HÓA BẰNG SỐ
GIẢNG VIÊN: DƯƠNG THANH QUỲNH
Họ Tên: Lưu Văn Kiên Lớp: 58CLC1 MSSV:3123558
Bài 1.1: Tính bằng tay các biểu thức sau rồi thử lại bằng matlab
function Bai11
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
Trang 2disp('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')
ĐÁP ÁN
================================================================ LAM QUEN VOI MATLAB
================================================================ Bai1.1
Phan a
ans =
1.3333
Phan b
ans =
53.6000
Phan c
ans =
6
Phan d
ans =
2.2500
Phan e
ans =
81
Phan f
ans =
2.5000
Phan g
Trang 3ans =
2
Phan h
ans =
2.5000
Phan k
ans =
2
ket thuc bai 1
Bài 2.2: Cho x= [3 1 5 7 9 2 6] dự đoán kết quả các dòng lệnh sau và thử lại bằng matlab
function Bai22
disp('===================================================================') disp(' VECTOR ') disp('===================================================================') 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
======================================================= VECTOR
======================================================= Bai2.2
cho x=[ 3 1 5 7 9 2 6]
Trang 4x =
3 1 5 7 9 2 6 phan a
ans =
5
phan b
ans =
3 1 5 7 9 2 6 phan c
ans =
3 1 5 7 9 2 6 phan d
ans =
3 1 5 7 9 2
phan e
ans =
2 7 1
phan f
ans =
3 2 1 3 3
Trang 5phan g
ans =
33
Bài 3.3: Ma Trận
function Bai33
disp('===================================================================') disp(' MA TRAN ') disp('===================================================================') disp('Bai3.3')
disp('Cho ma tr?n A = [2 4 1 ; 6 7 2 ; 3 5 9]')
A = [2 4 1 ; 6 7 2 ; 3 5 9]
disp(' Gán cho vector x là dòng th? nh?t c?a A')
x=A(1,1:3)
disp('Gán cho ma tr?n y là hai dòng còn l?i (cu?i) c?a A')
y=A(2:3,:)
disp('Tính t?ng theo dòng ma tr?n A ')
sum(A)
disp('Tính t?ng theo c?t ma tr?n A ')
sum(A')
disp('gia tri lon nhat')
max(max(A))
disp('gia tri nho nhat')
min(min(A))
disp('Tính t?ng các ph?n t? c?a A ')
sum(sum(A))
ĐÁP ÁN
===================================================================
MA TRAN
===================================================================
Bai3.3
Cho ma tran A = [2 4 1 ; 6 7 2 ; 3 5 9]
A =
2 4 1
6 7 2
3 5 9
Gán cho vector x là dòng thu nhat cua A
x =
2 4 1
Gán cho ma tran y là hai dòng còn lai (cuoi) cua A
y =
Trang 66 7 2
3 5 9
Tính t?ng theo dòng ma tr?n A
ans =
11 16 12
Tính t?ng theo c?t ma tr?n A
ans =
7 15 17
gia tri lon nhat
ans =
9
gia tri nho nhat
ans =
1
Tính t?ng các ph?n t? c?a A
ans =
39
Bài 4.1: Vẽ đồ thị hàm số f(x)=sin(1/x) với 0.01<x<0.1
function Bai41
disp('===================================================================') disp(' VE DO THI ') disp('===================================================================') disp('Bai4.1')
disp(' Ve ?o 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
===================================================================
VE DO THI
===================================================================
Trang 7Ve ?o 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 Bai55
disp('===================================================================') disp(' BIEU THUC LOGIC ') disp('===================================================================') 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))
Trang 8disp('gan nhung phan tu 0 thanh 1:')
X(X==0)=1
=================================================================== BIEU THUC LOGIC
=================================================================== 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 =
Trang 91 -2 8 9
5 6 1 25
1 -11 3 7
1 4 -3 2
Bài 6.1: Dự đoán kết quả rồi thử lại bằng matlab
function Bai61
disp('===================================================================') disp(' CAU LENH IF ') disp('===================================================================') disp('Bai6.1')
disp(' Du doan ket qua va thu lai bang matlab')
disp(' tinh m:')
n=input(' nhap gia tri cua n: n=')
if n>1
m=n+1
else
m=n-1
end
disp(' an tiep phim bat ky de lam phan b')
pause
n=input(' nhap gia tri cua n: n=')
if n>1
m=n+1
else
m=n-1
end
disp('an tiep phim bat ky de lam phan c')
pause
n=input(' nhap gia tri cua n: n=')
if n>1
m=n+1
else
m=n-1
end
end
ĐÁP ÁN
=============================================================== CAU LENH IF
=============================================================== Bai6.1
Du doan ket qua va thu lai bang matlab
tinh m:
nhap gia tri cua n: n=3
n =
3
m =
Trang 104
an tiep de lam phan b
nhap gia tri cua n: n=3
n =
3
m =
4
an tiep de lam phan c
nhap gia tri cua n: n=6
n =
6
m =
7
Bài 7.8
disp('===================================================================') disp(' VONG LAP ') disp('===================================================================') 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
end
if 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')
Trang 11disp(' tri trung binh cua vector ?o ')
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
===================================================================
VONG LAP
===================================================================
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 ?o
x =
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 =
Trang 12so 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.5:
disp('===================================================================') disp(' LAP TRINH ') disp('===================================================================') disp('Bai8.5')
n=input('nhap Do lon:');
s=[];
if n==1|n==2
s=1;
else
s=[1 1]
for i=3:n
s= [s (s(i-1)+s(i-2))];
end
end
disp(' 10 so Fibonacci dau tien co gia tri lan luot la:')
s
disp('an phim bat ky de lam tiep phan b')
pause
m=input('nhap Do lon:');
a=[];
if m==1|m==2
a=1;
else
a=[1 1]
for i=3:m
a= [a (a(i-1)+a(i-2))];
end
end
disp(' 50 so Fibonacci dau tien co gia tri lan luot la:')
a
disp('ti so giua Fn/Fn-1 cua day so tren la:')
for i=2:m
A=a(i)/a(i-1)
end
A
disp('Nhan xet: Ta thay tu so A thu 13 tro di co gia tri bang voi ty le vang') ĐÁP ÁN
Trang 13LAP TRINH
===================================================================
Bai8.5
nhap Do lon:10
s =
1 1
10 so Fibonacci dau tien co gia tri lan luot la:
s =
1 1 2 3 5 8 13 21 34 55
an phim bat ky de lam tiep phan b
nhap Do lon:50
a =
1 1
50 so Fibonacci dau tien co gia tri lan luot la:
a =
1.0e+010 *
Columns 1 through 12
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
Columns 13 through 24
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
Columns 25 through 36
0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0001 0.0002 0.0004 0.0006 0.0009 0.0015
Columns 37 through 48
0.0024 0.0039 0.0063 0.0102 0.0166 0.0268 0.0433 0.0701 0.1135 0.1836 0.2971 0.4808
Columns 49 through 50
0.7779 1.2586
ti so giua Fn/Fn-1 cua day so tren la:
A =
Trang 141
A =
2
A =
1.5000
A =
1.6667
A =
1.6000
A =
1.6250
A =
1.6154
A =
1.6190
A =
1.6176
A =
1.6182
A =
1.6180
A =
1.6181
Trang 15A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
Trang 161.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
Trang 17A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
Trang 18A =
1.6180
A =
1.6180
A =
1.6180
A =
1.6180
Nhan xet: Ta thay tu so A thu 13 tro di co gia tri bang voi ty le vang
>>
Bai 9.5
function Bai95
disp('===================================================================') disp(' BAI TAP SYMBOLIC MATH ') disp('===================================================================') 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
===================================================================
BAI TAP SYMBOLIC MATH
===================================================================
Bai 9.5
Trang 19Hinh cau tam O ban kinh bang 3: