TRƯỜNG ĐẠI HỌC XÂY DỰNG BAN 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ẢNG VIÊN : DƯƠNG THANH QUỲNH
Họ Tên : Nguyễn Tiến Long
Lớp : 63CLC1
MSSV : 130063
Bài 1.1
function Bai11
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')
Trang 2Đáp Án
>> Bai11
===================================================================
| 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
Trang 3ans =
81
Phan f
ans =
2.5000 Phan g ans =
2
Phan h ans =
9.5000 Phan k ans =
2
ket thuc bai 1
Bài 2.2
Trang 4function Bai22
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
>> Bai22
===================================================================
| Bai 2: VECTOR |
===================================================================
cho x=[ 3 1 5 7 9 2 6]
x =
3 1 5 7 9 2 6
Phan a
ans =
5
Phan b
Trang 5ans =
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
Phan g
ans =
Trang 633
Bài 3.4
function Bai34
disp('Tao ra ma tran 4x4 co gia tri nguyen nam trong
khoang[-10,10]')
A=round(rand(4,4)*(10 -(-10))+ (-10))
disp('Cong moi phan tu cua ma tran cho 15')
A+15
disp('Binh phuong moi phan tu cua ma tran')
A.^2
disp('Cong them 10 vao cac phan tu o dong 1 va dong 2') A([1 2],:)=A([1 2],:)+10
disp('Cong them 10 vao cac phan tu o cot 1 va cot 4')
A(:,[1 4])=A(:,[1 4])+10
Đáp Án
>> Bai34
===================================================================
| Bai 3: Ma Tran |
===================================================================
Tao ra ma tran 4x4 co gia tri nguyen nam trong khoang[-10,10]
A =
6 3 9 9
8 -8 9 0
-7 -4 -7 6
8 1 9 -7
Cong moi phan tu cua ma tran cho 15
ans =
Trang 721 18 24 24
23 7 24 15
8 11 8 21
23 16 24 8
Binh phuong moi phan tu cua ma tran
ans =
36 9 81 81
64 64 81 0
49 16 49 36
64 1 81 49
Cong them 10 vao cac phan tu o dong 1 va dong 2
A =
16 13 19 19
18 2 19 10
-7 -4 -7 6
8 1 9 -7
Cong them 10 vao cac phan tu o cot 1 va cot 4
A =
26 13 19 29
28 2 19 20
3 -4 -7 16
18 1 9 3
Trang 8Bài 4.7
function Bai47
disp(' Ve hai ham y=x^2 va y =sin3x tren cung mot do thi')
x = 0:pi/100:2*pi
plot(x,x.^2)
hold on
plot(x,sin(3.*x))
Bài 5.7
function Bai57
disp('Tao vector randperm(35)')
x= randperm(35)
(x<6)=2
y((x>=6)&(x<20))=x((x>=6)&(x<20))-4
y((x>=20)&(x<=35))=36-x((x>=20)&(x<=35))
Đáp Án
>> Bai57
===================================================================
| Bai 5: Bieu Thuc Logic |
===================================================================
Trang 9Tao vector randperm(35)
x =
Columns 1 through 17
2 31 25 6 27 15 19 14 12 28 26 29 8 13 22 7 24 Columns 18 through 34
30 10 35 34 23 18 17 21 3 9 4 20 5 16 32 11 33 Column 35
1
ans =
1×35 logical array
Columns 1 through 26
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Columns 27 through 35
0 0 0 0 0 0 0 0 0
y =
Trang 10Columns 1 through 17
0 0 0 2 0 11 15 10 8 0 0 0 4 9 0 3 0
Columns 18 through 33
0 6 0 0 0 14 13 0 0 5 0 0 0 12 0 7
y =
Columns 1 through 17
0 5 11 2 9 11 15 10 8 8 10 7 4 9 14 3 12
Columns 18 through 34
6 6 1 2 13 14 13 15 0 5 0 16 0 12 4 7 3
>>
Bài 6.5
function bai65
T= input('T=')
if (T>0)&(T<100)
disp('Gia tri h(T):')
h=T-10
elseif T>100
disp('Gia tri h(T):')
h=0.45*T+900
end
Đáp Án
Trang 11>> bai65
===================================================================
| Bai 6: Cau Lenh If |
=================================================================== T=80
T =
80
Gia tri h(T):
h =
70
T=200
T =
200
Gia tri h(T):
h =
990
Trang 12Bài 7.8
disp('Co bao nhieu phan tu dung truoc mot phan tu co gia tri nam trong khoang 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('Co bao nhieu phan tu trong khoang[0.01,0.5] ') disp('gia tri trung binh cua vertor ')
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.01 0.5] truoc gia tri dau tien lon hon trung binh')
else
y= x((x(1:n-1)>=0.01)&(x(1:n-1)<=0.5));
if isempty(y)==1
disp('Ko co so nao thuoc khoang [0.01 0.5] truoc gia tri dau tien lon hon trung binh')
else
strcat('co :',num2str(length(y)),' so thuoc khoang [0.01 0.5] truoc gia tri dau tien lon hon trung binh') end
end
Đáp Án
>> bai78
Co bao nhieu phan tu dung truoc mot phan tu co gia tri nam trong khoang 0.8 den 0.85
Trang 13Vo nghiem
Co bao nhieu phan tu trong khoang[0.01,0.5]
gia tri trung binh cua vertor
x =
0.9421 0.9561 0.5752 0.0598 0.2348 0.3532 0.8212 0.0154 0.0430 0.1690
x =
0.9421 0.9561 0.5752 0.0598 0.2348 0.3532 0.8212 0.0154 0.0430 0.1690
TrungBinh =
0.4170
ans =
'so thu :1 la so dau tien co gia tri lon hon trung binh'
Ko co so nao thuoc khoang [0.01 0.5] truoc gia tri dau tien lon hon trung binh
>>
Bài 8.1
function bai8_1
Trang 14d=input('nhap so dien su dung: ');
T=0;
G=0;
if d<=100
T=d*550;
G=T+T*0.1
disp(['tien dien phai tra la: ',num2str(G),' vnd'])
elseif (d>100)&(d<=200)
T=100*550+(d-100)*1200;
G=T+T*0.1;
disp(['tien dien phai tra la: ',num2str(G),' vnd'])
elseif (d>200)&(d<=300)
T=100*550+100*1200+(d-200)*1800;
G=T+T*0.1;
disp(['tien dien phai tra la: ',num2str(G),' vnd'])
else %>300
T=100*550+100*1200+100*1800+(d-300)*2000;
G=T+T*0.1;
disp(['tien dien phai tra la: ',num2str(G),' vnd'])
end
Đáp Án
>> bai8_1
nhap so dien su dung: 60
G = 36300
tien dien phai tra la: 36300 vnd
- nhap so dien su dung: 150
tien dien phai tra la: 126500 vnd
- nhap so dien su dung: 270
tien dien phai tra la: 331100 vnd
- nhap so dien su dung: 360
tien dien phai tra la: 522500 vnd
>>
Bài 9.5
function Bai95
disp('Bai 9.5')
disp('Hinh cau tam O ban kinh bang 3:')
syms a b;
Trang 15x=3*sin(a)*cos(b); y=3*sin(a)*sin(b); z=3*cos(a);
ezsurf(x,y,z);
axis tight ;
colorbar;
end