83
for i=1:3 %do each color band in turn
band=lower(abc(i));
%get (i)th input and make lower case
if strncmp(band,'bla',3) % black (compare min # of)
vals(i)=0; % chars for unique match)
elseif strncmp(band,'br',2) %brown
vals(i)=1;
elseif strncmp(band,'r',1) %red
vals(i)=2;
elseif strncmp(band,'o',1) %orange
vals(i)=3;
elseif strncmp(band,'y',1) %yellow
vals(i)=4;
elseif strncmp(band,'gre',3) %green
vals(i)=5;
elseif strncmp(band,'blu',3) %blue
vals(i)=6;
elseif strncmp(band,'v',1) %violet
vals(i)=7;
elseif strncmp(band,'gra',3) %gray
vals(i)=8;
elseif strncmp(band,'w',1) %white
vals(i)=9;
else
error(['Unknown Color Band.'])
end
end
if vals(1)==0
error('First Color Band Cannot Be Black.')
end
r=(10*vals(1)+vals(2))*10^vals(3);
Sö dông hµm nµy cho mét vµi vÝ dô:
>> resistor('brown', 'black', 'red')
ans=
1000
------------------------oOo-----------------------
ch−¬ng 13
PH¢N TÝCH D÷ LIÖU
Bëi v× MATLAB lµ mét øng dông h−íng ma trËn nªn nã dÔ dµng thùc hiÖn c¸c ph©n tÝch thèng kª
trªn c¸c tËp d÷ liÖu, trong khi theo mÆc ®Þnh MATLAB coi c¸c tËp d÷ liÖu ®−îc l−u tr÷ trong c¸c
m¶ng cét, viÖc ph©n tÝch d÷ liÖu cã thÓ thùc hiÖn theo bÊt cø chiÒu nµo. §ã lµ trõ khi ®−îc chØ ®Þnh
theo mét c¸ch kh¸c, c¸c cét cña mét m¶ng d÷ liÖu thÓ hiÖn c¸c th«ng sè ®o kh¸c nhau, mçi hµng thÓ
84
hiÖn mét gi¸ trÞ mÉu cña c¸c th«ng sè ®o ®ã. VÝ dô gi¶ sö nhiÖt ®é ban ngµy (tÝnh theo ®é C) cña 3
thµnh phè tÝnh trong mét th¸ng (31 ngµy ®−îc ghi l¹i vµ g¸n cho mét biÕn lµ temps trong mét script
M_file, khi ch¹y M_file th× gi¸ trÞ cña temps ®−îc ®a vµo m«i tr−êng MATLAB, thùc hiÖn c«ng viÖc
nµy, biÕn temps chøa:
>> temps
temps=
12 8 18
15 9 22
12 5 19
14 8 23
12 6 22
11 9 19
15 9 15
8 10 20
19 7 18
12 7 18
14 10 19
11 8 17
9 7 23
8 8 19
15 8 18
8 10 20
10 7 17
12 7 22
9 8 19
12 8 21
12 8 20
10 9 17
13 12 18
9 10 20
10 6 22
14 7 21
12 5 22
13 7 18
15 10 23
13 11 24
12 12 22
Mçi hµng chøa nhiÖt ®é cña mét ngµy nµo ®ã, cßn mçi cét chøa nhiÖt ®é cña mét thµnh phè. §Ó cho
d÷ liÖu trë lªn dÔ dµng h¬n, h·y gâ vµo nh sau:
>> d=1:31; % number the days of the month
>> plot(d,temps)
>> xlabel('Day of month')
>> ylabel('Celsius')
>> title('Daily High Tempratures in three Cities')
85
H×nh 13.1
LÖnh
plot
võa dïng trªn ®©y minh ho¹ thªm mét c¸ch sö dông. BiÕn d lµ mét vector dµi 31,
trong khi biÕn temps lµ mét ma trËn 31x3. Cho tr−íc nh÷ng d÷ liÖu nµy, lÖnh
plot
sÏ trÝc mçi cét
cña biÕn temps cho vµo d.
§Ó minh ho¹ mét vµi kh¶ n¨ng ph©n tÝch d÷ liÖu cña MATLAB, h·y xÐt c¸c lÖnh sau, dùa trªn d÷
liÖu vÒ nhiÖt ®é ®· cho:
>> avg_temp = mean(temps)
avg_temp=
11.9677 8.2258 19.8710
VÝ dô trªn chØ ra r»ng thµnh phè thø 3 lµ cã nhiÖt ®é trung b×nh cao nhÊt, ë ®©y MATLAB ®· tÝnh
nhiÖt ®é trung b×nh cña mçi cét mét c¸ch riªng rÏ. NÕu tÝnh trung b×nh ë c¶ 3 thµnh phè th×:
>> avg_avg = mean(avg_temp)
avg_avg=
13.3548
Khi mµ c¸c gi¸ trÞ ®Çu vµo trong mét hµm ph©n tÝch d÷ liÖu lµ mét vector hµng hay cét th×
MATLAB chØ ®¬n gi¶n lµ tiÕn hµnh c¸c phÐp to¸n trªn vector vµ tr¶ vÒ gi¸ trÞ sè.
B¹n còng cã thÓ dïng m¶ng ®Ó thùc hiÖn c«ng viÖc nµy:
>> avg_temp = mean(temps,1) % Gièng nh− trªn, tÝnh cho c¸c cét
avg_temp =
11.9677 8.2258 19.8710
86
>> avr_tempr = mean(temps,2) % TÝnh cho mçi hµng
avr_tempr =
12.6667
15.3333
12.0000
15.0000
13.3333
13.0000
13.0000
12.6667
14.6667
12.3333
14.3333
12.0000
13.0000
11.6667
13.6667
12.3333
11.3333
13.6667
12.0000
13.6667
13.3333
12.0000
14.3333
13.0000
12.6667
14.0000
13.0000
12.6667
16.0000
16.0000
15.3333
§©y lµ gi¸ trÞ nhiÖt ®é trung b×nh ë c¶ ba thµnh phè trong tõng ngµy.
XÐt bµi to¸n t×m sù chªnh lÖch nhiÖt ®é cña mçi thµnh phè so víi gi¸ trÞ trung b×nh, cã nghÜa lµ
avg_temp(i) ph¶i bÞ trõ ®i bëi cét thø
i
cña biÕn temps. B¹n kh«ng thÓ ra mét c©u lÖnh nh− sau:
>> temps-avg_temp
??? Error using ==> -
Matrix dimensions must agree.
Bëi v× thao t¸c nµy kh«ng ph¶i lµ c¸c thao t¸c ®· ®Þnh nghÜa trªn m¶ng (temps lµ mét m¶ng 31x3,
cßn avg_temp lµ mét m¶ng 1x3). Cã lÏ c¸ch dïng vßng lÆp
for
lµ ®¬n gi¶n nhÊt:
>> for i = 1:3
tdev(:,i) = temps(:,i)- avg_temp(i);
end
>> tdev
tdev =
0.0323 -0.2258 -1.8710
87
3.0323 0.7742 2.1290
0.0323 -3.2258 -0.8710
2.0323 -0.2258 3.1290
0.0323 -2.2258 2.1290
-0.9677 0.7742 -0.8710
3.0323 0.7742 -4.8710
-3.9677 1.7742 0.1290
7.0323 -1.2258 -1.8710
0.0323 -1.2258 -1.8710
2.0323 1.7742 -0.8710
-0.9677 -0.2258 -2.8710
-2.9677 -1.2258 3.1290
-3.9677 -0.2258 -0.8710
3.0323 -0.2258 -1.8710
-3.9677 0.7742 0.1290
-1.9677 -1.2258 -2.8710
0.0323 -1.2258 2.1290
-2.9677 -0.2258 -0.8710
0.0323 -0.2258 1.1290
0.0323 -0.2258 0.1290
-1.9677 0.7742 -2.8710
1.0323 3.7742 -1.8710
-2.9677 1.7742 0.1290
-1.9677 -2.2258 2.1290
2.0323 -1.2258 1.1290
0.0323 -3.2258 2.1290
1.0323 -1.2258 -1.8710
3.0323 1.7742 3.1290
1.0323 2.7742 4.1290
0.0323 3.7742 2.1290
Khi thùc hiÖn ph−¬ng ph¸p nµy ta thÊy nã chËm h¬n so víi c¸c c©u lÖnh ®−îc MATLAB thiÕt kÕ
riªng ®Ó dïng cho m¶ng. Khi ta nh©n b¶n biÕn avg_temp ®Ó kÝch th−íc cña nã b»ng víi kÝch th−íc
cña temps. Sau ®ã thùc hiÖn phÐp trõ th× sÏ nhanh h¬n rÊt nhiÒu:
>> tdev = temps - avg_temp(ones(31,1),:)
tdev =
0.0323 -0.2258 -1.8710
3.0323 0.7742 2.1290
0.0323 -3.2258 -0.8710
2.0323 -0.2258 3.1290
0.0323 -2.2258 2.1290
-0.9677 0.7742 -0.8710
3.0323 0.7742 -4.8710
-3.9677 1.7742 0.1290
7.0323 -1.2258 -1.8710
0.0323 -1.2258 -1.8710
2.0323 1.7742 -0.8710
-0.9677 -0.2258 -2.8710
-2.9677 -1.2258 3.1290