2쪽
Updatesofts.com Ebook Team 1 Matlab toμn tËpMatlab toμn tËpMatlab toμn tËpMatlab toμn tËp cμi ®Ætcμi ®Ætcμi ®Ætcμi ®Æt 1. MATLAB for WIN Yªu cÇu hÖ thèngYªu cÇu hÖ thèngYªu cÇu hÖ thèngYªu cÇu hÖ thèng • HÖ thèng IBM hoÆc t−¬ng thÝch 100% víi bé vi xö lÝ 486 Intel céng víi bé ®ång xö lÝ to¸n häc 487 ( ngo¹i trõ 486 DX cã bé xö lÝ bªn trong ), Pentium hoÆc Pentium Pro Processor. • Microsoft Window 95 hoÆc Window NT. a) æ CD ROMa) æ CD ROMa) æ CD ROMa) æ CD ROM - Bé ®iÒu phèi ®å ho¹ 8 bit vμ card mμn h×nh ( 256 mμu ®ång thêi ) - Kho¶ng trèng ®Üa ®ñ ®Ó cμi ®Æt vμ ch¹y c¸c tuú chän. Sù yªu cÇu ®Üa cøng thay ®æi tuú theo kÝch cì c¸c partition vμ c¸c tÖp trî gióp help ®−îc cμi ®Æt trùc tiÕp theo tuú chän. Qu¸ tr×nh cμi ®Æt sÏ th«ng b¸o cho b¹n biÕt tØ mØ vÒ dung l−îng ®Üa yªu cÇu. VÝ dô: Partition víi mét liªn cung mÆt 0 cÇn 25 MB cho riªng MATLAB vμ 50 MB cho c¶ MATLAB vμ HELP. Partition víi liªn cung 64 KB cÇn 115 MB cho riªng MATLAB vμ 250 MB cho c¶ MATLAB vμ HELP. b ) Bé nhí.b ) Bé nhí.b ) Bé nhí.b ) Bé nhí. Microsoft WIndow 95: 8 MB tèi thiÓu vμ 16 MB khuyÕn nghÞ. Microsoft WIN NT 3.51 hoÆc 4.0: 12 MB tèi thiÓu vμ 16 MB khuyÕn nghÞ. C¸c khuyÕn nghÞ • Bé nhí phô vμo (Bé nhí bæ sung: additional Memory). • VØ m¹ch t¨ng tèc ®å ho¹ bæ trî cho Microsoft Window. • M¸y in trî gióp cho Microsoft Window. • VØ m¹ch ©m thanh trî gióp cho Microsoft Window. • Microsoft Word 7.0 hoÆc h¬n ( nÕu b¹n cã ý ®Þnh sö dông MATLAB NoteBook ). • Tr×nh biªn dÞch Watcom C, Borland, Microsoft (x©y dùng file MEX). • Netscape Navigator 2.0 hoÆc version cao h¬n hoÆc Microsoft Internet Explorer 3.0 ®Ó ch¹y MATLAB Help Desk. Qu¸ tr×nh cμi ®ÆtQu¸ tr×nh cμi ®ÆtQu¸ tr×nh cμi ®ÆtQu¸ tr×nh cμi ®Æt 1. §Æt ®Üa vμo æ CD. Trªn WIN 95 ch−¬ng tr×nh SETUP b¾t ®Çu ch¹y tù ®éng nÕu nh− MATLAB cha ®−îc cμi tõ tr−íc. Cßn kh«ng, nhÊn ®óp vμo biÓu t−îng setup.exesetup.exe ®Ó b¾t ®Çusetup.exesetup.exe qu¸ tr×nh cμi ®Æt. 2. ChÊp nhËn hay bá ®i nh÷ng khuyÕn c¸o vÒ cÊp ®¨ng kÝ phÇn mÒm trªn mμn h×nh. NÕu chÊp nhËn b¹n míi cã thÓ b¾t ®Çu qu¸ tr×nh cμi ®Æt. 3. Trªn Custumer Information, nhËp vμo tªn b¹n, ®Þa chØ cña b¹n. Tªn kh«ng ®−îc qu¸ 30 kÝ tù. NhÊn nót NEXT. 4. NhÊn vμo c¸c hép trèng thμnh phÇn dÊu ‘v‘ nÕu nh− b¹n muèn tuú chän ®ã vμ nhÊn tiÕp nÕu b¹n cã ý ®Þnh kh«ng muèn tuú chän ®ã ( cã thÓ thªm vμo sau nμy nÕu muèn ThiNganHang.com
3쪽
Updatesofts.com Ebook Team 2 ). Trªn mμn h×nh hiÓn thÞ C:\MATLAB lμ th− môc ®Ých mÆc ®Þnh cña qu¸ tr×nh cμi ®Æt. NÕu b¹n muèn cμi ®Æt vμo th− môc kh¸c hoÆc ®æi tªn th− môc th× b¹n lùa chän Browse. MATLAB cho Macintosh. MATLAB cho m¸y Macintosh ch¹y ®−îc trªn: • Mäi m¸y Macintosh cã cÊu h×nh ®ñ m¹nh ( power Macintosh ). • Mäi Macintosh ®−îc trang bÞ bé vi xö lÝ 68040 ( bé ®ång xö lÝ to¸n häc bªn trong ). • Mäi m¸y Macintosh ®−îc trang bÞ bé vi xö lÝ 68020 hoÆc 68030 vμ bé ®ång xö lÝ to¸n häc 68881 hoÆc 68882. Yªu cÇu tèi thiÓu ®Ó ch¹y MATLAB. • §Üa cøng trèng tèi thiÓu 26 MB, cÇn thªm 60 MB cho hÖ thèng tuú chon HELP trùc tuyÕn. • 16 MB cho ph©n vïng bé nhí. • æ CD ROM. • Color Quick Draw. ---------------oOo---------------- Ch−¬ng1 GIíI THIÖU chung B©ygiê b¹n ®· cμi ®Æt xong, chóng ta h·y xem MATLAB cã thÓ lμm ®−îc nh÷ng g×. Trong phÇn nμy chóng ta sÏ tr×nh bμy mét sè nh÷ng øng dông cña nã; v× ®Ó tr×nh bμy tÊt c¶ nh÷ng øng dông cña MATLAB sÏ rÊt dμi vμ tèn thêi gian. NÕu b¹n ®äc quyÓn h−íng dÉn nμy, b¹n sÏ thÊy MATLAB lμ ng«n ng÷ rÊt m¹nh ®Ó gi¶i quyÕt nh÷ng vÊn ®Ò quan träng vμ khã kh¨n cña b¹n. Nã sÏ rÊt h÷u Ých khi b¹n ®äc phÇn h−íng dÉn c¬ b¶n v× nã sÏ cung cÊp cho b¹n nh÷ng kiÕn thøc c¬ b¶n ®Ó b¹n hiÓu râ MATLAB vμ ph¸t triÓn ®−îc nh÷ng kh¶ n¨ng cña m×nh sau nμy. Cã lÏ c¸ch dÔ nhÊt ®Ó h×ng dung vÒ MATLAB lμ nã cã ®Çy ®ñ c¸c ®Æc ®iÓm cña m¸y tÝnh c¸ nh©n: gièng nh− c¸c m¸y tÝnh c¬ b¶n, nã lμm tÊt c¶ c¸c phÐp tÝnh to¸n häc c¬ b¶n nh− céng, trõ, nh©n, chia; gièng nh− m¸y tÝnh kü thuËt, nã bao gåm: sè phøc, c¨n thøc, sè mò, logarithm, c¸c phÐp to¸n l−îng gi¸c nh− sine, cosine, tang; nã còng gièng nh− m¸y tÝnh cã kh¶ n¨ng lËp tr×nh, cã thÓ l−u tr÷, t×m kiÕm l¹i d÷ liÖu, còng cã thÓ t¹o, b¶o vÖ vμ ghi tr×nh tù c¸c lÖnh ®Ó tù ®éng phÐp to¸n khi gi¶i quyÕt c¸c vÊn ®Ò, b¹n cã thÓ so s¸nh logic, ®iÒu khiÓn thùc hiªn lÖnh ®Ó ®¶m b¶o tÝnh ®óng ®¾n cña phÐp to¸n. Gièng nh− c¸c m¸y tÝnh hiÖn ®¹i nhÊt, nã cho phÐp b¹n biÓu diÔn d÷ liÖu díi nhiÒu d¹ng nh−: biÓu diÔn th«ng th−êng, ma tr©n ®¹i sè, c¸c hμm tæ hîp vμ cã thÓ thao t¸c víi d÷ liÖu th−êng còng nh− ®èi víi ma trËn. Trong thùc tÕ MATLAB cßn øng dông rÊt réng r·i trong nhiÒu lÜnh vùc vμ nã còng sö dông rÊt nhiÒu c¸c phÐp tÝnh to¸n häc. Víi nh÷ng ®Æc ®iÓm ®ã vμ kh¶ n¨ng th©n thiÖn víi ng−êi sö dông nªn nã dÔ dμng sö dông h¬n c¸c ng«n ng÷ kh¸c nh− Basic, Pascal, C. Nã cung cÊp mét m«i tr−êng phong phó cho biÓu diÔn d÷ liÖu, vμ cã kh¶ n¨ng m¹nh mÏ vÒ ®å ho¹, b¹n cã thÓ t¹o c¸c giao diÖn riªng cho ng−êi sö dông(GUIs) ®Ó g¶i quyÕt nh÷ng vÊn ®Ò riªng cho m×nh. Thªm vμo ®ã MATLAB ®−a ra nh÷ng c«ng cô ®Ó gi¶i quyÕt nh÷ng vÊn ®Ò ®Æc biÖt, gäi lμ Toolbox (hép c«ng cô). VÝ dô Student Edition cña MATLAB bao gåm c¶ Toolbox ®iÒu khiÓn hÖ thèng, Toolbox xö lÝ tÝn hiÖu, Toolbox biÓu t−îng to¸n häc. Ngoμi ra b¹n cã thÓ t¹o Toolbox cho riªng m×nh. Víi nh÷ng kh¶ n¨ng m¹nh mÏ, réng lín cña MATLAB nªn nã rÊt cÇn thiÕt cho b¹n b¾t ®Çu tõ phÇn c¬ b¶n. Sau ®©y chóng ta sÏ nghiªn cøu tõng phÇn, vμ cuèn s¸ch nμy sÏ gióp ThiNganHang.com
4쪽
Updatesofts.com Ebook Team 3 b¹n hiÓu ®−îc chóng. Tr−íc tiªn, mét c¸ch ®¬n gi¶n nhÊt lμ chóng ta quan niÖm nh− lμ mét m¸y tÝnh c¬ b¶n, tiÕp theo lμ nh− m¸y tÝnh kü thuËt vμ nh− m¸y tÝnh cã thÓ lËp tr×nh ®−îc, cuèi cïng lμ nh− m¸y tÝnh hiÖn ®¹i nhÊt. B»ng c¸ch quan niÖm nμy b¹n sÏ dÔ dμng hiÓu ®ù¬c nh÷ng c¸ch mμ MATLAB gi¶i quyÕt nh÷ng vÊn ®Ò th«ng th−êng vμ xem MATLAB gi¶i quyÕt nh÷ng vÊn ®Ò vÒ sè phøc mÒm dÎo nh− thÕ nμo. Tuú thuéc vμo kiÕn thøc cña b¹n, b¹n cã thÓ t×m thÊy nh÷ng phÇn trong cuèn s¸ch h−- íng dÉn nμy høng thó hay buån tÎ... Khi b¹n ch¹y ch−¬ng tr×nh MATLAB, nã sÏ t¹o mét hoÆc nhiÒu cöa sæ trªn mμn h×nh cña b¹n, vμ cöa sæ lÖnh (command) lμ cöa sæ chÝnh ®Ó b¹n giao tiÕp víi MATLAB, cöa sæ nμy xuÊt hiÖn nh h×nh díi ®©y. C¸c kÝ tù ‘EDU>>‘ lμ dÊu nh¾c cña MATLAB trong student MATLAB. Trong c¸c version kh¸c cña MATLAB, dÊu nh¾c ®¬n gi¶n chØ lμ ‘>>‘. Khi cöa sæ lÖnh xuÊt hiÖn, lμ cöa sæ ho¹t ®éng, con trá xuÊt hiÖn bªn ph¶i dÊu nh¾c nh− ë h×nh d−íi. Con trá vμ dÊu nh¾c nμy cña MATLAB b¸o r»ng MATLAB ®ang ®îi ®Ó thùc hiÖn lÖnh. H×nh 1.1H×nh 1.1H×nh 1.1H×nh 1.1 Cöa sæ lÖnh cña Student MATLABCöa sæ lÖnh cña Student MATLABCöa sæ lÖnh cña Student MATLABCöa sæ lÖnh cña Student MATLAB 1.1 C¸c phÐp to¸n ®¬n gi¶n1.1 C¸c phÐp to¸n ®¬n gi¶n1.1 C¸c phÐp to¸n ®¬n gi¶n1.1 C¸c phÐp to¸n ®¬n gi¶n Gièng nh− m¸y tÝnh ®¬n gi¶n th«ng th−êng, MATLAB cã thÓ thùc hiÖn c¸c phÐp to¸n ®¬n gi¶n, nh− vÝ dô d−íi ®©y: Mary ®Õn mét cöa hμng v¨n phßng phÈm vμ mua 4 côc tÈy, 25 xu mét côc, 6 tËp vë, 52 xu mét tËp, hai cuén b¨ng ®μi, 99 xu mét cuén. H·y tÝnh xem Mary mua bao nhiªu vËt, vμ tæng sè tiÒn lμ bao nhiªu? NÕu dïng m¸y tÝnh th«ng th−êng, ta vμo c¸c sè: 4 + 6 + 2 = 12 ( vËt) 4x25 + 6x52 + 2x99 = 610 (xu) H×nh 1.2H×nh 1.2H×nh 1.2H×nh 1.2 Cöa sæ lÖnh cña MATLAB version 5.2Cöa sæ lÖnh cña MATLAB version 5.2Cöa sæ lÖnh cña MATLAB version 5.2Cöa sæ lÖnh cña MATLAB version 5.2 Trong MATLAB chóng ta cã thÓ gi¶i quyÕt vÊn ®Ò nμy theo nhiÒu c¸ch. Tr−íc tiªn gièng nh− m¸y tÝnh ë trªn, chóng ta cã thÓ tÝnh: >> 4 + 6 + 2 ans= 12 >> 4*25 + 6*52 + 2*99 ans= 610 Chó ý r»ng MATLAB kh«ng chó ý ®Õn nh÷ng kho¶ng trèng, cho tÊt c¶ c¸c phÇn, vμ phÐp nh©n cã møc ®é −u tiªn cao h¬n phÐp céng. Vμ mét chó ý kh¸c lμ MATLAB gäi kÕt qu¶ ans (viÕt t¾t cña answer) cho c¶ hai phÐp tÝnh. Nh− ®· nãi ë trªn, vÊn ®Ò trªn cã thÓ gi¶i quyÕt b»ng c¸ch chøa c¸c th«ng tin vμo biÕn cña MATLAB: >> erasers = 4 erasers= 4 ThiNganHang.com
5쪽
Updatesofts.com Ebook Team 4 >> pads = 6 pads= 6 >> tape = 2; >> iterms = erases + pads + tape iterms= 12 >> cost = erases*25 + pads*52 + tape*99 cost= 610 ë ®©y chóng ta t¹o 3 biÕn MATLAB: erases, pads, tape ®Ó chøa sè l−îng mçi lo¹i vËt. Sau khi vμo c¸c gi¸ trÞ cho c¸c biÕn nμy, MATLAB hiÓn thÞ kÕt qu¶ ra mμn h×nh, trõ tr−êng hîp biÕn tape. DÊu hai chÊm ®»ng sau c©u lÖnh ‘>> tape = 2;‘ th«ng b¸o cho MATLAB nhËn gi¸ trÞ g¸n nh−ng kh«ng hiÓn thÞ ra mμn h×nh. Cuèi cïng kh¸c víi gäi kÕt qu¶ ans, chóng ta yªu cÇu MATLAB gäi kÕt qu¶ tæng sè c¸c vËt lμ iterms, vμ tæng sè tiÒn lμ cost. T¹i mçi b−íc MATLAB ®Òu ®a ra c¸c th«ng tin. V× cã l−u gi÷ c¸c biÕn nªn chóng ta cã thÓ yªu cÇu MATLAB tÝnh gi¸ trÞ trung b×nh cho mçi vËt: >> everage_cost = cost/iterms everage_cost= 50.8333 Bëi v× everage cost cã hai tõ, mμ MATLAB yªu cÇu biÕn chØ cã mét tõ, nªn chóng ta dïng dÊu g¹ch d−íi ®Ó nèi hai tõ nμy thμnh mét tõ. Ngoμi c¸c phÐp tÝnh trªn, MATLAB cßn cã mét sè phÐp tÝnh c¬ b¶n kh¸c nh− b¶ng d−íi ®©y: PhÐp tÝnh BiÓu t−îng VÝ dô PhÐp céng, a + b + 5 + 3 PhÐp trõ, a - b - 7 - 4 PhÐp nh©n, a.b * 18*24 PhÐp chia, ab / hoÆc \ 56/8 = 8\ 56 PhÐp luü thõa, ab ^ 5^2 Trong c¸c phÐp to¸n trªn cã møc ®é −u tiªn kh¸c nhau, khi tÝnh tõ tr¸i sang ph¶i cña mét dßng gåm nhiÒu lÖnh th× phÐp to¸n luü thõa cã møc ®é −u tiªn cao nhÊt, tiÕp theo lμ phÐp nh©n vμ phÐp chia cã møc ®é −u tiªn b»ng nhau cuèi cïng lμ phÐp céng vμ phÐp trõ còng cã møc ®é −u tiªn b»ng nhau. 1.2 Kh«ng gian lμm viÖc cña MATLAB1.2 Kh«ng gian lμm viÖc cña MATLAB1.2 Kh«ng gian lμm viÖc cña MATLAB1.2 Kh«ng gian lμm viÖc cña MATLAB Còng nh− b¹n lμm viÖc víi cöa sæ LÖnh, MATLAB nhí c¸c lÖnh b¹n gâ vμo còng nh− c¸c gi¸ trÞ b¹n g¸n cho nã hoÆc nã ®−îc t¹o lªn. Nh÷ng lÖnh vμ biÕn nμy ®−îc gäi lμ l−u gi÷ trong kh«ng gian lμm viÖc cña MATLAB, vμ cã thÓ ®−îc gäi l¹i khi b¹n muèn. VÝ dô, ®Ó kiÓm tra gi¸ trÞ cña biÕn tape, tÊt c¶ nh÷ng g× b¹n ph¶i lμm lμ yªu cÇu MATLAB cho biÕt b»ng c¸ch ®¸nh vμo tªn biÕn t¹i dÊu nh¾c: >> tape ThiNganHang.com
6쪽
Updatesofts.com Ebook Team 5 tape= 2 NÕu b¹n kh«ng nhí tªn biÕn, b¹n cã thÓ yªu cÇu MATLAB cho danh s¸ch c¸c biÕn b»ng c¸ch d¸nh lÖnh whowhowhowho tõ dÊu nh¾c lÖnh: >> who Your variables are: ans cost iterms tape average_cost erasers pads Chó ý r»ng MATLAB kh«ng ®−a ra gi¸ trÞ cña tÊt c¶ c¸c biÕn, nÕu b¹n muèn biÕt gi¸ trÞ, b¹n ®¸nh vμo tªn biÕn t¹i dÊu nh¾c lÖnh cña MATLAB. §Ó gäi l¹i c¸c lÖnh b¹n ®· dïng, MATLAB dïng c¸c phÝm mòi tªn (↑ ↓ ) trªn bμn phÝm cña b¹n. VÝ dô ®Ó gäi l¹i lÖnh b¹n gâ vμo lóc gÇn hiÖn t¹i nhÊt, b¹n nhÊn phÝm mòi tªn ↓, tiÕp tôc nhÊn phÝm nμy, nã sÏ l¹i gäi tiÕp lÖnh tr−íc ®ã, NÕu b¹n dïng phÝm mòi tªn ↑ nã sÏ gäi l¹i lÖnh tõ lÖnh ®Çu tiªn cho ®Õn lÖnh gÇn hiÖn t¹i nhÊt. C¸c phÝm mòi tªn ← vμ → cã thÓ dïng ®Ó thay ®æi vÞ trÝ con trá trong dßng lÖnh t¹i dÊu nh¾c cña MATLAB, nh− vËy chóng ta cã thÓ söa dßng lÖnh, thªm n÷a, chóng ta cã thÓ dïng chuét cïng víi bé nhí ®Öm ®Ó c¾t, copy, d¸n, vμ söa v¨n b¶n t¹i dÊu nh¾c cña dßng lÖnh. 1.3 BiÕn1.3 BiÕn1.3 BiÕn1.3 BiÕn Gièng nh− nh÷ng ng«n ng÷ lËp tr×nh kh¸c, MATLAB cã nh÷ng quy ®Þnh riªng vÒ tªn biÕn. Tr−íc tiªn tªn biÕn ph¶i lμ mét tõ, kh«ng chøa dÊu c¸ch, vμ tªn biÕn ph¶i cã nh÷ng quy tu©n thñ nh÷ng quy t¾c sau: Quy ®Þnh vÒ tªn biÕn Tªn biÕn cã ph©n biÖt ch÷ hoa ch÷ th−êng. VÝ dô: Iterms, iterms, itErms, vμ ITERMS lμ c¸c biÕn kh¸c nhau Tªn biÕn cã thÓ chøa nhiÒu nhÊt 31 kÝ tù, cßn c¸c kÝ tù sau kÝ tù thø 31 bÞ lê ®i. VÝ dô: howaboutthisveriablename Tªn biÕn b¾t ®Çu ph¶i lμ ch÷ c¸i, tiÕp theo cã thÓ lμ ch÷ sè, sè g¹ch d−íi VÝ dô: how_about_this_veriable_name, X51483. a_b_c_d_e KÝ tù chÊm c©u kh«ng ®−îc phÐp dïng v× nã cã nh÷ng ý nghÜa ®Æc biÖt Cïng víi nh÷ng quy ®Þnh trªn, MATLAB cã nh÷ng biÕn ®Æc biÖt trong b¶ng sau: C¸c biÕn ®Æc biÖt Gi¸ trÞ ans Tªn biÕn mÆc ®Þnh dïng ®Ó tr¶ vÒ kÕt qu¶ pi = 3.1415.. Eps Sè nhá nhÊt, nh− vËy dïng céng víi 1 ®Ó ®−îcsè nhá nhÊt lín h¬n 1 flops Sè cña phÐp to¸n sè thùc inf §Ó chØ sè v« cïng nh kÕt qu¶ cña 1/0 NaN hoÆc nan Dïng ®Ó chØ sè kh«ng x¸c ®Þnh nh− kÕt qu¶ cña 0/0 i (vμ) j i = j = nargin Sè c¸c ®èi sè ®a vμo hμm ®−îc sö dông narout Sè c¸c ®èi sè hμm ®a ra realmin Sè nhá nhÊt cã thÓ ®−îc cña sè thùc realmax Sè lín nhÊt cã thÓ ®−îc cña sè thùc ThiNganHang.com
7쪽
Updatesofts.com Ebook Team 6 Nh− b¹n cã thÓ t¹o mét biÕn cña MATLAB, vμ b¹n còng cã thÓ g¸n l¹i gi¸ trÞ cho mét hoÆc nhiÒu biÕn. VÝ dô: >> erases = 4; >> pads = 6; >> tape = 2; >> iterms = eases + pads + tape iterms= 12 >> erases = 6 erases= 6 >> iterms iterms= 12 ë ®©y chóng ta sö dông l¹i vÝ dô trªn, chóng ta t×m ®−îc sè vËt mμ Mary ®· mua sau ®ã chóng ta thay ®æi sè côc tÈy lªn 6, gi¸ trÞ nμy sÏ ®Ì lªn gi¸ trÞ tr−íc cña nã lμ 4. Khi b¹n lμm nh− vËy, gi¸ trÞ cña iterms vÉn kh«ng thay ®æi, v× MATLAB kh«ng tÝnh l¹i iterms víi gi¸ trÞ míi cña erases. Khi MATLAB thùc hiÖn mét phÐp tÝnh, nã lÊy gi¸ trÞ cña c¸c biÕn hiÖn thêi, nªn nÕu b¹n muèn tÝnh gi¸ trÞ míi cña iterms, cost, average_cost, b¹n gäi l¹i c¸c lÖnh tÝnh c¸c gi¸ trÞ ®ã. §èi víi c¸c biÕn ®Æc biÖt ë trªn, nã cã s½n gi¸ trÞ, nh− vËy khi b¹n khëi ®éng MATLAB; nÕu b¹n thay ®æi gi¸ trÞ cña nã th× nh÷ng gi¸ trÞ ®Æc biÖt ban ®Çu sÏ bÞ mÊt cho ®Õn khi b¹n xo¸ biÕn ®ã ®i hoÆc khëi ®éng l¹i MATLAB. Do ®ã b¹n kh«ng nªn thay ®æi gi¸ trÞ cña biÕn ®Æc biÖt, trõ khi nã thùc sù cÇn thiÕt. C¸c biÕn trong kh«ng gian lμm viÖc cña MATLAB cã thÓ bÞ xo¸ kh«ng ®iÒu kiÖn b»ng c¸ch dïng lÖnh clearclearclearclear. VÝ dô: >> clear erases chØ xo¸ mét biÕn erases >> clear cost iterms xo¸ c¶ hai biÕn cost vμ iterms >> clear cl* dÊu * ®Ó chØ r»ng xo¸ tÊt c¶ c¸c biÕn b¾t ®Çu b»ng hai kÝ tù cl. >> clear xo¸ tÊt c¶ c¸c biÕn trong kh«ng gian lμm viÖc!. B¹n sÏ kh«ng ®−îc hái ®Ó x¸c nhËn c©u lÖnh nμy vμ tÊt c¶ c¸c biÕn ®· bÞ xo¸ kh«ng thÓ kh«i phôc l¹i. Cã thÓ nãi r»ng dïng lÖnh clearclearclearclear rÊt nguy hiÓm, v× vËy khi dïng lÖnh nμy b¹n nªn dïng ®óng vÞ trÝ. 1.4 C©u gi¶i thÝch (comment) vμ sù chÊm c©u1.4 C©u gi¶i thÝch (comment) vμ sù chÊm c©u1.4 C©u gi¶i thÝch (comment) vμ sù chÊm c©u1.4 C©u gi¶i thÝch (comment) vμ sù chÊm c©u TÊt c¶ c¸c v¨n b¶n ®»ng sau kÝ hiÖu phÇn tr¨m (%) ®Òu lμ c©u gi¶i thÝch. VÝ dô: >> erases = 4 % Sè côc tÈy. erases= ThiNganHang.com
8쪽
Updatesofts.com Ebook Team 7 4 BiÕn erases ®−îc g¸n gi¸ trÞ lμ 4, cßn tÊt c¶ kÝ hiÖu phÇn tr¨m vμ v¨n b¶n ®»ng sau nã ®Òu bÞ lê ®i. §Æc ®iÓm nμy gióp cho chóng ta dÔ theo dâi c«ng viÖc chóng ta ®ang lμm. NhiÒu lÖnh cã thÓ ®Æt trªn cïng mét hμng, chóng c¸ch nhau bëi dÊu phÈy hoÆc dÊu chÊm phÈy, nh−: >> erases = 4, pads = 6; tape = 2 erases= 4 tape= 2 dÊu phÈy ®Ó yªu cÇu MATLAB hiÓn thÞ kÕt qu¶ trªn mμn h×nh; cßn dÊu chÊm phÈy lμ kh«ng hiÓn thÞ kÕt qu¶ trªn mμn h×nh. >> average_cost = cost/ ... iterms average_cost= 50.83333 Nh− vÝ dô trªn, ta cã thÓ dïng dÊu ba chÊm (...) ®Ó chØ c©u lÖnh ®−îc tiÕp tôc ë hμng d−íi, phÐp tÝnh thùc hiÖn ®−îc khi dÊu ba chÊm ng¨n c¸ch gi÷a to¸n tö vμ biÕn, nghÜa lμ tªn biÕn kh«ng bÞ ng¨n c¸ch gi÷a hai hμng: >> average_cost = cost/ it... erms ??? age_cost = cost/iterms Missing operator, coma, or semicolon. gièng nh− vËy, tr¹ng th¸i cña lêi gi¶i thÝch kh«ng thÓ tiÕp tôc: >> % Comments cannot be continued ... >> either ??? Undefined function or variable either. B¹n cã thÓ dõng ch−¬ng tr×nh b»ng c¸ch nhÊn ®ång thêi Ctrl vμ C. 1.5 Sè phøc1.5 Sè phøc1.5 Sè phøc1.5 Sè phøc Mét trong nh÷ng ®Æc ®iÓm m¹nh mÏ nhÊt cña MATLAB lμ lμm viÖc víi sè phøc. Sè phøc trong MATLAB ®−îc ®Þnh nghÜa theo nhiÒu c¸ch, vÝ dô nh− sau: >> c1 = 1 - 2i % ChÌn thªm kÝ tù i vμo phÇn ¶o. c1= 1.0000 - 2.0000i >> c1 = 1 - 2j % j ë ®©y t−¬ng tù nh− i ë trªn. c1= 1.0000 - 2.0000i >> c2 = 3*(2-sqrt(-1)*3) c2= ThiNganHang.com
9쪽
Updatesofts.com Ebook Team 8 6.0000 - 9.0000i >> c3 = sqrt(-2) c3= 0 + 1.4142i >> c4 = 6 + sin(.5)*i c4= 6.0000 + 0.4794i >> c5 = 6 + sin(.5)*j c5= 6.0000 + 0.4794i Trong hai vÝ dô cuèi, MATLAB mÆc ®Þnh gi¸ trÞ cña i = j = dïng cho phÇn ¶o. Nh©n víi i hoÆc j ®−îc yªu cÇu trong tr−êng hîp nμy, sin(.5)i vμ sin(.5)j kh«ng cã ý nghÜa ®èi víi MATLAB. Cuèi cïng víi c¸c kÝ tù i vμ j, nh− ë trong hai vÝ dô ®Çu ë trªn chØ lμm viÖc víi sè cè ®Þnh, kh«ng lμm viÖc ®−îc víi biÓu thøc. Mét sè ng«n ng÷ yªu cÇu sù ®iÒu khiÓn ®Æc biÖt cho sè phøc khi nã xuÊt hiÖn, trong MATLAB th× kh«ng cÇu nh− vËy. TÊt c¶ c¸c phÐp tÝnh to¸n häc ®Òu thao t¸c ®−îc nh− ®èi víi sè thùc th«ng th−êng: >> c6 = (c1 + c2)/c3 % Tõ c¸c d÷ liÖu ë trªn c6= -7.7782 - 4.9497i >> check_it_out = i^2 % B×nh ph−¬ng cña i ph¶i lμ -1 check_it_out= -1.0000 + 0.0000i trong vÝ dô nμy chØ cßn l¹i phÇn thùc, phÇn ¶o b»ng kh«ng. Chóng ta cã thÓ dïng hμm realrealrealreal vμ imagimagimagimag ®Ó kiÓm tra tõng phÇn thùc vμ ¶o. Chóng ta cã thÓ biÓu diÔn sè phøc d¹ng ®é lín vμ gãc (d¹ng cùc): M M.ej = a+bi ë trªn sè phøc ®−îc biÓu diÔn b»ng ®é lín M vμ gãc , quan hÖ gi÷a c¸c ®¹i l−îng nμy vμ phÇn thùc, phÇn ¶o cña sè phøc biÓu diÔn d−íi d¹ng ®¹i sè lμ: M = = tan-1(b/ a) a = Mcos b = Msin Trong MATLAB, ®Ó chuyÓn tõ d¹ng cùc sang d¹ng ®¹i sè, dïng c¸c hμm realrealrealreal, imagimagimagimag, vμ angleangleangleangle: >> c1 % Gäi l¹i c1 c1= 1.0000 - 2.0000i >> M_c1 = abs(c1) % TÝnh argument cña sè phøc M_c1= 2.2361 >> angle_c1 = angle(c1) % TÝnh gãc cña sè phøc theo radian angle_c1= ThiNganHang.com
10쪽
Updatesofts.com Ebook Team 9 -1.1071 >> deg_c1 = angle_c1*180/ pi % ChuyÓn tõ radian sang ®é -63.4349 >> real_c1 = real(c1) % TÝnh phÇn thùc real_c1= 1 >> imag_c1 = imag(c1) % TÝnh phÇn ¶o imag_c1= -2 -------------------oOo------------------ Ch−¬ng2 C¸C §ÆC tÝnh KÜ THUËT Gièng nh− hÇu hÕt c¸c m¸y tÝnh kü thuËt, MATLAB ®a ra rÊt nhiÒu c¸c hμm to¸n häc, kÜ thuËt th«ng dông, ngoμi ra MATLAB cßn cung cÊp hμng tr¨m c¸c hμm ®Æc biÖt vμ thuËt to¸n, nã rÊt h÷u Ých ®Ó gi¶i quyÕt c¸c vÊn ®Ò khoa häc. TÊt c¶ c¸c hμm nμy ®−îc liÖt kª trong online help, cßn ë ®©y chØ ®Ò cËp ®Õn nh÷ng hμm th«ng dông nhÊt. 2.1 C¸c hμm to¸n häc th«ng th−êng2.1 C¸c hμm to¸n häc th«ng th2.1 C¸c hμm to¸n häc th«ng th2.1 C¸c hμm to¸n häc th«ng th−êng−êng−êng C¸c hμm to¸n häc cña MATLAB ®−îc liÖt kª trong b¶ng d−íi ®©y, chóng ®Òu cã chung mét c¸ch gäi hμm nh− vÝ dô d−íi ®©y: >> x = sqrt(2)/2 x= 0.7071 >> y = sin(x) y= 0.7854 >> y_deg = y*180/pi y_deg= 45.0000 Nh÷ng lÖnh nμy ®Ó t×m mét gãc (tÝnh b»ng ®é) khi biÕt gi¸ trÞ hμm sin cña nã lμ / 2. TÊt c¶ c¸c hμm liªn quan ®Õn gãc cña MATLAB ®Òu lμm viÖc víi radian. B¶ng c¸c hμm: C¸c hμm th«ng th−êng ThiNganHang.com
11쪽
Updatesofts.com Ebook Team 10 abs(x) TÝnh argument cña sè phøc x acos(x) Hμm ng−îc cña cosine acosh(x) Hμm ng−îc cña hyperbolic cosine angle(x) TÝnh gãc cña sè phøc x asin(x) Hμm ng−îc cña sine asinh(x) Hμm ng−îc cña hyperbolic sine atan(x) Hμm ng−îc cña tangent atan2(x, y) Lμ hμm arctangent cña phÇn thùc cña x vμ y atanh(x) Hμm ng−îc cña hyperbolic tangent ceil(x) XÊp xØ d−¬ng v« cïng conj(x) Sè phøc liªn hîp cos(x) Hμm cosine cña x cosh(x) Hμm hyperbolic cosine cña x exp(x) Hμm ex fix(x) XÊp xØ kh«ng floor(x) XÊp xØ ©m v« cïng gdc(x, y) ¦íc sè chung lín nhÊt cña hai sè nguyªn xvμ y imag(x) Hμm tr¶ vÒ phÇn ¶o cña sè phøc lcm(x, y) Béi sè chung nhá nhÊt cña hai sè nguyªn x vμ y log(x) Logarithm tù nhiªn log10(x) Logarithm c¬ sè 10 real(x) Hμm tr¶ vÒ phÇn thùc cña x rem(x, y) PhÇn d− cña phÐp chia x/ y round(x) Hμm lμm trßn vÒ sè nguyªn tè sign(x) Hμm dÊu: tr¶ vÒ dÊu cña argument nh−: sign(1.2)=1; sign(-23.4)=-1; sign(0)=0 sin(x) Hμm tÝnh sine cña x sinh(x) Hμm tÝnh hyperbolic sine cña x sqrt(x) Hμm khai c¨n bËc hai tan(x) Tangent tanh(x) Hyperbolic tangent >> 4*atan(1) % Mét c¸ch tÝnh xÊp xØ gi¸ trÞ cña pi ans= 3.1416 >> help atant2 % Yªu cÇu gióp ®ì ®èi víi hμm atan2 ATAN2 four quadrant inverse tangent ATAN2(Y, X) is the four quadrant arctangent of the real parts of the elements of X and Y. -pi <= ATAN2(Y, X) <= pi see also ATAN. >> 180/pi*atan(-2/ 3) ans= -33.69 >> 180/pi*atan2(2, -3) ans= 146.31 ThiNganHang.com
12쪽
Updatesofts.com Ebook Team 11 >> 180/pi*atan2(-2, 3) ans= -33.69 >> 180/pi*atan2(2, 3) ans= 33.69 >> 180/pi*atan2(-2, -3) ans= -146.31 Mét sè vÝ dô kh¸c: >> y = sqrt(3^2 + 4^2) % TÝnh c¹nh huyÒn cña tam gi¸c pitago 3-4-5 y= 5 >> y = rem(23,4) % 23/4 cã phÇn d− lμ 3 y= 3 >> x = 2.6,y1 = fix(x),y2 = floor(x),y3 = ceil(x),y4 = round(x) x= 2.6000 y1= 2 y2= 2 y3= 3 y4= 3 >> gcd(18,81) % 9 lμ −íc sè chung lín nhÊt cña 18 vμ 81 ans= 9 >> lcm(18,81) % 162 lμ béi sè chung lín nhÊt cña 18 vμ 81 ans= 162 VÝ dô: ¦íc l−îng chiÒu cao cña ng«i nhμVÝ dô: ¦íc lVÝ dô: ¦íc lVÝ dô: ¦íc l−îng chiÒu cao cña ng«i nhμ−îng chiÒu cao cña ng«i nhμ−îng chiÒu cao cña ng«i nhμ VÊn ®Ò:VÊn ®Ò:VÊn ®Ò:VÊn ®Ò: Gi¶ thiÕt biÕt kho¶ng c¸ch tõ ngêi quan s¸t ®Õn ng«i nhμ lμ D, gãc tõ ngêi quan s¸t ®Õn ng«i nhμ lμ ; chiÒu cao cña ngêi quan s¸t lμ h. Hái ng«i nhμ cao bao nhiªu? Gi¶i ph¸p:Gi¶i ph¸p:Gi¶i ph¸p:Gi¶i ph¸p: Ta biÓu diÔn kÝch thøc nh− h×nh 2.1: (kh«ng thÊy) H×nh 2.1H×nh 2.1H×nh 2.1H×nh 2.1 Ng«i nhμ cã chiÒu cao lμ H + h, H lμ chiÒu dμi cña mét c¹nh cña tam gi¸c, chiÒu dμi nμy cã thÓ tÝnh ®−îc b»ng c«ng thøc quan hÖ gi÷a gãc vμ c¹nh cña tam gi¸c: tan() = ThiNganHang.com
13쪽
Updatesofts.com Ebook Team 12 Tõ ®ã ta cã chiÒu cao cña ng«i nhμ lμ h + H = h + D.tan() NÕu h =2meters, D =50meters, vμ lμ 60o, MATLAB sÏ ®a ra kÕt qu¶ lμ: >> h = 2 h = 2 >> theta = 60 theta = 60 >> D = 50 D = 50 >> buiding_height = h+D*atan(theta*pi/180) buiding_height = 54.3599 VÝ dô sù suy gi¶m do ph©n r· VÊn ®ÒVÊn ®ÒVÊn ®ÒVÊn ®Ò : Sù ph©n r· ph©n tö polonium cã chu kú ph©n r· lμ 140 ngμy, tøc lμ sau 140 ngμy th× l−îng poloniun cßn l¹i lμ 1/2 l−îng ban ®Çu. Hái nÕu ban ®Çu cã 10 grams polonium, nã sÏ cßn l¹i bao nhiªu sau 250 ngμy? Gi¶i quyÕtGi¶i quyÕtGi¶i quyÕtGi¶i quyÕt: Sau 1 chu kú ph©n r· hoÆc 140 ngμy, cßn l¹i 10x0.5 = 5 grams; sau 2 chu kú ph©n r· hoÆc 280 ngμy, cßn l¹i 5x0.5 = 10x(0.5)2 = 2.5grams, tõ ®ã ta cã kÕt qu¶ n»m trong kho¶ng 5 vμ 2.5 grams, vμ ta cã c«ng thøc tÝnh phÇn cßn l¹i sau kho¶ng thêi gian bÊt kú: khèi l−îng cßn l¹i = khèi l−îng ban ®Çu x(0.5)thêi gian/ chu kú vÝ dô thêi gian lμ 250 ngμy, vμ kÕt qu¶ MATLAB ®a ra lμ: >> initial_amount = 10; % Khèi l−îng ban ®Çu >> half_life = 140; % Chu kú ph©n r· >> time = 250; % Thêi gian tÝnh khèi l−îng >> amount_left = initial_*0.5^(time/half_life) amount_left= 2.9003 VÝ dô tÝnh to¸n vÒ l·i xuÊtVÝ dô tÝnh to¸n vÒ l·i xuÊtVÝ dô tÝnh to¸n vÒ l·i xuÊtVÝ dô tÝnh to¸n vÒ l·i xuÊt VÊn ®ÒVÊn ®ÒVÊn ®ÒVÊn ®Ò: B¹n ®ång ý mua «t« míi víi gi¸ 18,500 dollars. Ng−êi b¸n «t« ®a ra hai gi¶i ph¸p vÒ tμi chÝnh lμ: thø nhÊt, tr¶ 2.9% l·i xuÊt cña sè tiÒn trªn trong vßng 4 n¨m. Thø hai lμ tr¶ 8.9% l·i xuÊt cña sè tiÒn trªn trong vßng 4 n¨m vμ gi¸ b¸n ®−îc gi¶m ®i mét kho¶n lμ 1500 dollars. Hái víi gi¶i ph¸p nμo th× b¹n mua ®−îc «t« víi gi¸ rÎ h¬n? Gi¶i ph¸p:Gi¶i ph¸p:Gi¶i ph¸p:Gi¶i ph¸p: Sè tiÒn tr¶ hμng th¸ng lμ P, trªn tæng sè tiÒn lμ A dollars, tØ sè l·i xuÊt hμng th¸ng lμ R, tr¶ trong M th¸ng: P = A ThiNganHang.com
14쪽
Updatesofts.com Ebook Team 13 Tæng sè tiÒn ph¶i tr¶ sÏ lμ: T = PxM Gi¶i ph¸p MATLAB ®a ra lμ: >> format bank % Dïng d¹ng hiÓn thÞ ng©n hμng >> A = 18500; % Tæng sè tiÒn >> M = 12*4; % Sè th¸ng ph¶i tr¶ l·i >> FR = 1500; % TiÒn gi¶m gi¸ cña nhμ m¸y >> % Gi¶i ph¸p thø nhÊt >> R = (2.9/100)/12; % TØ lÖ l·i xuÊt hμng th¸ng >> P = A*(R*(1+R)^M/((1+R)^M - 1)) % Kho¶n tiÒn ph¶i tr¶ hμng th¸ng P= 408.67 >> T1 = P*M % Tæng gi¸ trÞ cña «t« T1= 19616.06 >> % Gi¶i ph¸p thø hai >> R = (8.9/100)/12; % TØ lÖ l·i xuÊt hμng th¸ng >> P = (A-FR)*(R*(1 + R)^M/((1+R)^M - 1)) % TiÒn ph¶i tr¶ hμng th¸ng P= 422.24 >> T2 = P*M % Tæng gi¸ trÞ cña «t« T2= 20267.47 >> Diff = T2 - T1 Diff= 651.41 Nh− vËy ta cã gi¶i ph¸p thø nhÊt gi¸ rÎ h¬n gi¶i ph¸p thø hai. VÝ dô: VÊn ®Ò nång ®é acidVÝ dô: VÊn ®Ò nång ®é acidVÝ dô: VÊn ®Ò nång ®é acidVÝ dô: VÊn ®Ò nång ®é acid VÊn ®Ò:VÊn ®Ò:VÊn ®Ò:VÊn ®Ò: Nh− mét phÇn cña qu¸ tr×nh s¶n xuÊt bé phËn cña vËt ®óc t¹i mét nhμ m¸y tù ®éng, bé phËn ®ã ®−îc nhóng trong n−íc ®Ó lμm nguéi, sau ®ã nhóng trong bån ®ùng dung dÞch acid ®Ó lμm s¹ch. Trong toμn bé cña qu¸ tr×nh nång ®é acid gi¶m ®i khi c¸c bé phËn ®- −îc lÊy ra khái bån acid v× khi nhóng bé phËn cña vËt ®óc vμo bån th× mét l−îng n−íc cßn b¸m trªn vËt ®óc khi nhóng ë bÓ tr−íc còng vμo theo vμ khi nhÊc ra khái bån mét l−îng acid b¸m theo vËt. §Ó ®¶m b¶o chÊt l−îng th× nång ®é acid ph¶i kh«ng ®−îc nhá h¬n mét l−îng tèi thiÓu. B¹n h·y b¾t ®Çu víi nång ®é dung dÞch lμ 90% th× nång ®é tèi thiªu ph¶i lμ 50%. L−îng chÊt láng thªm vμo vμ lÊy ®i sau mçi lÇn nhóng dao ®éng trong kho¶ng tõ 1% ®Õn 10%. Hái bao nhiªu bé phËn cã thÓ nhóng vμo bÓ dung dÞch acid tr−íc khi nång ®é cña nã gi¶m xuèng d−íi møc cho phÐp? Gi¶i ph¸p:Gi¶i ph¸p:Gi¶i ph¸p:Gi¶i ph¸p: Ban ®Çu nång ®é acid lμ initial_con = 90% = acid/ (acid + water) sau lÇn nhóng thø nhÊt nång ®é acid cßn: con = = = ThiNganHang.com
15쪽
Updatesofts.com Ebook Team 14 = ‘acid‘ lμ l−îng acid ban ®Çu trong dung dÞch, ‘water‘ lμ l−îng n−íc ban ®Çu trong dung dÞch, ‘lost‘ lμ l−îng phÇn tr¨m n−íc thªm vμo. Sè acid cßn l¹i trong dung dÞch sau lÇn nhóng thø nhÊt lμ: acid_left = NghÜa lμ, khi nhóng lÇn thø hai nång ®é dung dÞch sÏ lμ: con = = = TiÕp tôc qu¸ tr×nh nμy, sau n lÇn nhóng, nång ®é acid lμ: con = NÕu nång ®é acid cßn l¹i lμ møc tèi thiÓu chÊp nhËn ®−îc, sè lÇn nhóng cùc ®¹i sÏ lμ mét sè nguyªn b»ng hoÆc nhá h¬n n: n = Trong MATLAB gi¶i ph¸p sÏ lμ: >> initial_con = 90 initial_con= 90 >> min_con = 50 min_con= 50 >> lost = 0.01; >> n = floor(log( initial_con/min_con)/log(1+lost)) n= 59 Nh− vËy cã thÓ nhóng 59 lÇn tr−íc khi nång ®é acid gi¶m xuèng d−íi 50%. Chó ý hμm floorfloorfloorfloor dïng ®Ó lμm trßn sè n xuèng sè nguyªn gÇn nhÊt, vμ ë d©y ta còng cã thÓ dïng hμm logarithmlogarithmlogarithmlogarithm c¬ sè 10 vμ logarithmlogarithmlogarithmlogarithm c¬ sè 2 thay cho hμm logarithmlogarithmlogarithmlogarithm tù nhiªn ë trªn. -------------------oOo---------------------- ThiNganHang.com
16쪽
Updatesofts.com Ebook Team 15 ch−¬ng 3 NH÷NG §ÆC §IÓM CñA CöA Sæ LÖNH Cöa sæ lÖnh (comand) cña MATLAB cã rÊt nhiÒu nh÷ng ®Æc ®iÓm cÇn chó ý, mét sè chóng ®· ®−îc giíi thiÖu ë ch−¬ng tr−íc, vμ sau ®©y chóng ta t×m hiÓu râ h¬n vÒ chóng. 3.1 Qu¶n lÝ kh«ng gian lμm viÖc cña MATLAB3.1 Qu¶n lÝ kh«ng gian lμm viÖc cña MATLAB3.1 Qu¶n lÝ kh«ng gian lμm viÖc cña MATLAB3.1 Qu¶n lÝ kh«ng gian lμm viÖc cña MATLAB C¸c d÷ liÖu vμ biÕn ®−îc t¹o lªn trong cöa sæ lÖnh, ®−îc l−u trong mét phÇn gäi lμ kh«ng gian lμm viÖc cña MATLAB. Muèn xem tªn biÕn trong kh«ng gian lμm viÖc cña MATLAB ta dïng lÖnh who: >> who Your variables are: D h buiding_height theta C¸c biÕn nμy ®−îc dïng trong vÝ dô −íc l−îng chiÒu cao ng«i nhμ. §Ó xem chi tiÕt h¬n vÒ c¸c biÕn ta dïng lÖnh whos: >> whos Name Size Bytes Class D 1x1 8 double array buiding_height 1x1 8 double array h 1x1 8 double array theta 1x1 8 double array Grand total is 4 elements using 32 bytes Mçi biÕn ®−îc liÖt kª víi kÝch cì cña nã, sè bytes sö dông, vμ c¸c líp cña chóng (class), trong vÝ dô ®Æc biÖt nμy, c¸c biÕn ®Òu lμ sè ®¬n, cã ®é chÝnh x¸c hai sè sau dÊu phÈy. LÖnh whoswhoswhoswhos ®Æc biÖt cã Ých khi nghiªn cøu ®Õn phÇn m¶ng vμ c¸c kiÓu d÷ liÖu kh¸c. Ngoμi c¸c hμm nμy, trong môc Show WorkspaceShow Workspace trong b¶ng chän fileShow WorkspaceShow Workspace file t¹o ra cöa sæfilefile GUI gäi lμ Workspace BrowserWorkspace Browser, nã chøa c¸c th«ng tin t−¬ng tù nh− lÖnhWorkspace BrowserWorkspace Browser whoswhoswhoswhos. Thªm n÷a nã t¹o cho b¹n kh¶ n¨ng xo¸, lμm s¹ch c¸c biÕn mμ b¹n chän. Cöa sæ nμy còng cã thÓ t¹o b»ng c¸ch nhÊn nót Workspace BrowserWorkspace Browser, trªn thanh c«ng cô cña cöa sæ lÖnh.Workspace BrowserWorkspace Browser Nh− ®· tr×nh bμy ë trªn, lÖnh clearclearclearclear cã thÓ xo¸ biÕn tõ kh«ng gian lμm viÖc cña MATLAB. VÝ dô: >> clear h D % Xo¸ c¸c biÕn h vμ D >> who Your variables are: buiding_height theta ThiNganHang.com
17쪽
Updatesofts.com Ebook Team 16 C¸c tuú chän kh¸c cña hμm clear chóng ta cã thÓ t×m hiÓu thªm b»ng lÖnh help: >> help clear CLEAR Clear variables and functions from memory. CLEAR removes all variables from the workspace. CLEAR VARIABLES does the same thing. CLEAR GLOBAL removes all global variables. CLEAR FUNCTIONS removes all compiled M-functions. CLEAR MEX removes all links to MEX-files. CLEAR ALL removes all variables, globals, functions and MEX links. CLEAR VAR1 VAR2 ... clears the variables specified. The wildcard character '*' can be used to clear variables that match a pattern. For instance, CLEAR X* clears all the variables in the current workspace that start with X. If X is global, CLEAR X removes X from the current workspace, but leaves it accessible to any functions declaring it global. CLEAR GLOBAL X completely removes the global variable X. CLEAR FUN clears the function specified. If FUN has been locked by MLOCK it will remain in memory. CLEAR ALL also has the side effect of removing all debugging breakpoints since the breakpoints for a file are cleared whenever the m-file changes or is cleared. Use the functional form of CLEAR, such as CLEAR('name'), when the variable name or function name is stored in a x©u. See also WHO, WHOS, MLOCK, MUNLOCK. Cuèi cïng, khi lμm viÖc trong kh«ng gian lμm viÖc cña MATLAB, nã th−êng thuËn tiÖn ®Ó ghi hoÆc in mét b¶n sao c«ng viÖc cña b¹n, lÖnh diarydiarydiarydiary ghi d÷ liÖu ng−êi dïng ®−a vμo vμ cöa sæ lÖnh vμ ®−a ra file v¨n b¶n d¹ng m· ASCII cã tªn lμ diary trong th− môc hiÖn t¹i. >> diary frame % ghi d÷ liÖu vao file frame >> diary off % kÕt thóc lÖnh diary vμ ®ãng file Khi cöa sæ lÖnh ®−îc chän, chän printprint...printprint tõ b¶ng chän filefile ®Ó in mét b¶n cña cöa sæfilefile lÖnh, b¹n cã thÓ dïng chuét ®Ó lùa chän phÇn m×nh muèn ghi, chän Pint SelectionPint Selection... tõ b¶ngPint SelectionPint Selection chän filefile, ®Ó in mét phÇn v¨n b¶n ®· lùa chän.filefile 3.2 Ghi vμ phôc håi d÷ liÖu3.2 Ghi vμ phôc3.2 Ghi vμ phôc3.2 Ghi vμ phôc håi d÷ liÖuhåi d÷ liÖuhåi d÷ liÖu §Ó nhí c¸c biÕn MATLAB cã thÓ ghi vμ gäi l¹i d÷ liÖu tõ file trong m¸y tÝnh cña b¹n. Môc Workspace asWorkspace as... trong b¶ng chän fileWorkspace asWorkspace as file më hép chuÈn héi tho¹i ®Ó ghi tÊt c¶ c¸c biÕnfilefile ThiNganHang.com
18쪽
Updatesofts.com Ebook Team 17 hiÖn t¹i. Gièng nh− vËy, trong môc Load WorkspaceLoad Workspace trong b¶ng chän fileLoad WorkspaceLoad Workspace file më hép héi tho¹ifilefile ®Ó gäi l¹i tÊt c¶ c¸c biÕn mμ ta ®· ghi l¹i tõ kh«ng gian lμm viÖc tr−íc, nã kh«ng lμm mÊt c¸c biÕn nμy trong kh«ng gian lμm viÖc hiÖn t¹i. Khi ta gäi l¹i c¸c biÕn, mμ c¸c biÕn nμy trïng tªn víi c¸c biÕn trong kh«ng gian lμm viÖc cña MATLAB, nã sÏ thay ®æi gi¸ trÞ cña c¸c biÕn theo gi¸ trÞ cña c¸c biÕn gäi ra tõ file. NÕu b¶ng chän file kh«ng thuËn tiÖn hoÆc kh«ng ®¸p øng ®−îc nh÷ng yªu cÇu cña b¹n, MATLAB cung cÊp hai lÖnh savesavesavesave vμ loadloadloadload, nã thùc hiÖn mét c¸ch mÒm dÎo h¬n, trong tr−êng hîp ®Æc biÖt, lÖnh savesavesavesave cho phÐp b¹n ghi mét hoÆc nhiÒu h¬n mét biÕn tuy theo sù lùa chon cña b¹n. VÝ dô: >> save Chøa tÊt c¶ c¸c biÕn trong MATLAB theo kiÓu nhÞ ph©n trong file MATLAB.mat >> save data chøa tÊt c¶ c¸c biÕn trong MATLAB theo kiÓu nhÞ ph©n trong fle data.mat. >> save data erasers pads tape -ascii Ghi c¸c biÕn erasers, pads, tape trong d¹ng m· ASCII 8 sè trong file data. File d¹ng m· ASCII cã thÓ söa ®æi b»ng bÊt cø ch−¬ng tr×nh so¹n th¶o v¨n b¶n nμo, chó ý r»ng file ASCII kh«ng cã phÇn më réng .mat. >> save data erasers pads tape -ascii -double Ghi c¸c biÕn erasers, pads, tape d¹ng ASCII 16 sè trong file data. LÖnh loadloadloadload còng dïng víi có ph¸p t−îng tù. 3.3 Khu«n d¹ng hiÓn thÞ sè3.3 Khu«n d¹ng hiÓn thÞ sè3.3 Khu«n d¹ng hiÓn thÞ sè3.3 Khu«n d¹ng hiÓn thÞ sè Khi MATLAB hiÓn thÞ kÕt qu¶ d¹ng sè, nã tu©n theo mét sè quy ®Þnh sau: MÆc ®Þnh, nÕu kÕt qu¶ lμ sè nguyªn th× MATLAB hiÓn thÞ nã lμ mét sè nguyªn, khi kÕt qu¶ lμ mét sè thùc th× MATLAB hiÓn thÞ sè xÊp xØ víi bèn ch÷ sè sau dÊu phÈy, cßn c¸c sè d¹ng khoa häc th× MATLAB hiÓn thÞ còng gièng nh− trong c¸c m¸y tÝnh khoa häc. B¹n cã thÓ kh«ng dïng d¹ng mÆc ®Þnh, mμ t¹o mét khu«n d¹ng riªng tõ môc Preferences, trong b¶ng chän filePreferencesPreferencesPreferences file, cã thÓ mÆc ®Þnh hoÆc ®¸nh d¹ng xÊp xØ t¹i dÊu nh¾c.filefile Chóng ta dïng biÕn average_cost ( trong vÝ dô tr−íc) lμm vÝ dô, d¹ng sè nμy lμ: LÖnh cña MATLAB Average_cost Chó thÝch format short 50.833 5 sè format long 50.83333333333334 16 sè format short e 5.0833e+01 5 sè víi sè mò format long e 5.083333333333334e+01 16 sè víi sè mò format short g 50.833 chÝnh x¸c h¬n format short hoÆc format short e format long g 50.83333333333333 chÝnh x¸c h¬n format long hoÆc format long e format hex 40496aaaaaaaaaab hÖ c¬ sè 16 ThiNganHang.com
19쪽
Updatesofts.com Ebook Team 18 format bank 50.83 hai sè hÖ 10 format + + d−¬ng, ©m hoÆc b»ng kh«ng format rat 305/ 6 d¹ng ph©n sè Mét chó ý quan träng lμ MATLAB kh«ng thay ®æi sè khi ®Þnh l¹i khu«n d¹ng hiÓn thÞ ®−îc chän, mμ chØ thay ®æi mμn h×nh thay ®æi. ---------------oOo----------------- Ch−¬ng 4 Script M_files Mét vÊn ®Ò ®¬n gi¶n lμ, yªu cÇu cña b¹n t¹i dÊu nh¾c cña MATLAB trong cöa sæ lÖnh lμ nhanh vμ hiÖu qu¶. Tuy nhiªn v× sè lÖnh t¨ng lªn, hoÆc khi b¹n muèn thay ®æi gi¸ trÞ cña mét hoÆc nhiÒu biÕn vμ thùc hiÖn l¹i mét sè lÖnh víi gi¸ trÞ míi, nÕu cø ®¸nh lÆp l¹i t¹i dÊu nh¾c cña MATLAB th× sÏ trë lªn buån tÎ, do vËy MATLAB cung cÊp mét gi¶i ph¸p cho vÊn ®Ò nμy lμ: nã cho phÐp b¹n thay thÕ c¸c lÖnh cña MATLAB b»ng mét file v¨n b¶n ®¬n gi¶n, vμ yªu cÇu MATLAB më file vμ thùc hiÖn lÖnh chÝnh x¸c nh− lμ ®¸nh t¹i dÊu nh¾c cña MATLAB t¹i cöa sæ lÖnh, nh÷ng file nμy gäi lμ script filescript file, hoÆc ®¬n gi¶n lμ M_filescript filescript file M_file. Danh tõM_fileM_file "script" ®Ó chØ r»ng thùc tÕ MATLAB ®äc tõ file kÞch b¶n t×m thÊy trong file. Danh tõ "M_file" ®Ó chØ r»ng tªn script file ®ã ph¶i kÕt thóc b»ng phÇn më réng lμ '.m' nh vÝ dô example1.m. §Ó t¹o mét script M_file, chän NewNew trong b¶ng chän fileNewNew file vμ chän M_filefilefile M_file. Thñ tôc nμyM_fileM_file sÏ t¹o ra mμn h×nh so¹n th¶o, vμ b¹n cã thÓ ®¸nh ®−îc c¸c lÖnh cña MATLAB trong ®ã. VÝ dô d−íi ®©y lμ c¸ch lÖnh trong vÝ dô −íc l−îng chiÒu cao ng«i nhμ ë tr−íc: function example1 % example1.m VÝ dô −íc l−îng chiÒu cao ng«i nhμ h = 2 theta = 60 D = 50; building_height = h + D*tan(theta*pi/180) B¹n cã thÓ ghi vμ l−u gi÷ file nμyb»ng c¸ch chän SaveSave tõ b¶ng chän fileSaveSave file. Khi b¹n ghi lªnfilefile file chó ý ph¶i ®¸nh tªn file trïng víi tªn hμm (example) kh«ng cÇn ®¸nh vμo phÇn më réng, MATLAB tù g¸n vμo cho nã. Khi ®ã tõ dÊu nh¾c ta cã thÓ ®¸nh: >> example1 h= 2 theta= 60 building_height= 54.3599 Khi MATLAB diÔn gi¶i c¸c tr¹ng th¸i cña example1 ë trªn, nã sÏ ®−îc nãi kü h¬n ë ch−¬ng sau, nh−ng mét c¸ch ng¾n gän, MATLAB dïng c¸c tr¹ng th¸i cña biÕn MATLAB hiÖn t¹i vμ t¹o lªn c¸c lÖnh cña nã, b¾t ®Çu b»ng tªn M_file. NghÜa lμ, nÕu example1 kh«ng ph¶i lμ biÕn hiÖn t¹i, hoÆc mét lÖnh MATLAB x©y dùng lªn, MATLAB më file example1.m ThiNganHang.com
20쪽
Updatesofts.com Ebook Team 19 (nÕu nã t×m thÊy) vμ tÝnh gi¸ trÞ c¸c lÖnh t×m thÊy chØ khi chóng ta vμo c¸c th«ng sè chÝnh x¸c t¹i dÊu nh¾c cña cöa sæ lÖnh. Nh− ®· thÊy lÖnh trong M_file truy cËp ®Õn tÊt c¶ c¸c biÕn trong kh«ng gian lμm viÖc cña MATLAB, vμ tÊt c¶ c¸c biÕn trong M_file trë thμnh mét phÇn cña kh«ng gian lμm viÖc. B×nh th−êng c¸c lÖnh ®äc trong M_file kh«ng ®−îc hiÓn thÞ nh− lμ nã ®−îc tÝnh trong cöa sæ lÖnh, nh−ng lÖnh echo onecho onecho onecho on yªu cÇu MATLAB hiÓn thÞ hoÆc lÆp l¹i lÖnh ®èi víi cöa sæ lÖnh nh− chóng ta ®· ®äc vμ tÝnh. TiÕp theo b¹n cã thÓ ®o¸n ®−îc lÖnh echo offecho offecho offecho off lμm g×. Gièng nh− vËy, lÖnh echoechoechoecho lÆp l¹i bëi chÝnh nã lμm thay ®æi chÝnh tr¹ng th¸i cña nã. Víi ®Æc ®iÓm nμy cña M_file b¹n cã thÓ thay ®æi l¹i néi dung cña file, vÝ dô b¹n cã thÓ më M_file example1.m thay ®æi l¹i c¸c gi¸ trÞ cña h, D, hoÆc theta, ghi l¹i file ®ã vμ yªu cÇu MATLAB tÝnh l¹i lÖnh trong file. Thªm n÷a, b»ng c¸ch t¹o M_file, c¸c lÖnh cña b¹n ®−îc l−u trªn ®Üa vμ cã thÓ øng dông vÒ sau khi b¹n cÇn. Nh÷ng øng dông cña chØ dÉn cña MATLAB gióp chóng ta hiÓu ®−îc khi dïng script file nh− trong example1.m, chØ dÉn cho phÐp b¹n l−u gi÷ cïng c¸c lÖnh trong script file, v× vËy b¹n nhí ®−îc nh÷ng lÖnh ®ã lμm g× khi b¹n nh×n l¹i file sau ®Êy. Thªm n÷a, dÊu chÊm phÈy ®»ng sau c©u lÖnh kh«ng cho hiÓn thÞ kÕt qu¶, tõ ®ã b¹n cã thÓ ®iÒu chØnh script file ®a ra nh÷ng kÕt qu¶ cÇn thiÕt. V× nh÷ng øng dông cña script file, MATLAB cung cÊp mét sè hμm ®Æc biÖt cã Ých khi b¹n sö dông trong M_file: C¸c hμm M_file disp(ans) HiÓn thÞ c¸c kÕt qu¶ mμ kh«ng hiÖn tªn biÕn echo §iÒu khiÓn cöa sæ lÖnh lÆp l¹i c¸c lÖnh cña script file input Sö dông dÊu nh¾c ®Ó ®a d÷ liÖu vμo keyboard Trao ®iÒu khiÓn t¹m thêi cho bμn phÝm pause Dõng l¹i cho ®Õn khi ng−êi dïng nhÊn mét phÝm bÊt kú pause(n) Dõng l¹i n gi©y waitforbuttonpress Dõng l¹i cho ®Õn khi ng−êi dïng nhÊn chuét hoÆc phÝm. Khi lÖnh cña MATLAB kh«ng kÕt thóc b»ng dÊu chÊm phÈy, kÕt qu¶ cña lÖnh ®−îc hiÓn thÞ trªn cöa sæ lÖnh cïng víi tªn biÕn. §«i lóc nã thuËn tiÖn khi kh«ng cho hiÖn tªn biÕn, trong MATLAB ta dïng lÖnh disp ®Ó thùc hiÖn viÖc nμy: >> h % C¸ch truyÒn thèng ®Ó hiÖn kÕt qu¶ h= 2 >> disp(h) % HiÖn kÕt qu¶ kh«ng cã tªn biÕn 2 §Ó gióp b¹n so¹n th¶o script file khi tÝnh to¸n cho nhiÒu tr−êng hîp, lÖnh inputinputinputinput cho phÐp b¹n t¹o c©u nh¾c ®Ó vμo d÷ liÖu ®−îc an toμn. VÝ dô example1.m víi nh÷ng phÇn ®−- îc söa: function example1 % example1.m VÝ dô −íc l−îng chiÒu cao ng«i nhμ h = 2 theta = 60 ThiNganHang.com
앞 20페이지만 표시됩니다. 전체를 보려면 자료를 다운로드하세요.
댓글 (0)
댓글이 없습니다. 첫 댓글을 남겨보세요!