Giáo trình mathlab toàn tập - Chương 13

PHÂN TíCH Dữ LIệUBở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ố...
Nội dung trích xuất từ tài liệu:
Giáo trình mathlab toàn tập - Chương 13 83for 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.]) endendif vals(1)==0 error(First Color Band Cannot Be Black.)endr=(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¸cm¶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Ø ®Þnhtheo 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Ó 84hiÖ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 3thµ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 scriptM_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Öcnµy, biÕn temps chøa:>> tempstemps= 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 22Mç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è. §Ó chod÷ 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étcñ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ÝnhnhiÖ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étavg_temp = 11.9677 8.2258 19.8710 86 % TÝnh cho mçi hµng>> avr_tempr = mean(temps,2)avr_tempr = 12.6667 15.3333 12.0000 15.0000 ...