Giáo trình Matlab V.5.3
Matlab là một phần mềm toán học của hãng Mathworks để tính toán trên các số và có tính trực quan rất cao. Matlab đã qua nhiều phiên bản, giáo trình này giới thiệu phiên bản 7.0 (release 14). Matlab là...Quy trình cài đặt Matlab cũng tương tự như việc cài đặt các chương trình phần mềm khác, chỉ cần theo các hướng dẫn và bổ sung thêm các thông số cho phù hợp.
www.vietsupport.com -visit- www.dientuvietnam.net
[email protected] CM Soft 70 NCT F2 Q10
PhÇn 1
c¬ së Mat lab
Ch−¬ng 1:
Cμi ®Æt matlab
1.1.Cμi ®Æt ch−¬ng tr×nh:
Qui tr×nh cμi ®Æt Matlab còng t−¬ng tù nh− viÖc cμi ®Æt c¸c ch−¬ng tr×nh phÇn
mÒm kh¸c, chØ cÇn theo c¸c h−íng dÉn vμ bæ xung thªm c¸c th«ng sè cho phï hîp.
1.1.1 Khëi ®éng windows.
1.1.2 Do ch−¬ng tr×nh ®−îc cÊu h×nh theo Autorun nªn khi g¾n dÜa CD vμo æ ®Üa th×
ch−¬ng tr×nh tù ho¹t ®éng, cöa sæ ®ãn mõng hiÖn lªn trong gi©y l¸t. KÝch vμo Next
chuyÓn sang cöa sæ cμi ®Æt kÕ tiÕp.
Trang 1
[email protected] CM Soft 70 NCT F2 Q10
1.1.3 Cöa sæ thø hai chØ c¸c th«ng tin vÒ b¶n quyÒn cña ch−¬ng tr×nh. KÝch Yes ®Ó sang
cöa sæ cμi ®Æt kÕ tiÕp.
1.1.4 Trong cöa sæ Costumer Information, nÕu cÇn thiÕt ®¸nh tªn hä, c«ng ty. NhËp m·
kho¸ cña ch−¬ng tr×nh vμo khung Persnal license password. KÝch tiÕp Next.
1.1.5 Trong cöa sæ select Matlab Components, kÝch xo¸ nh÷ng thμnh phÇn kh«ng cÇn
thiÕt sö dông trong ch−¬ng tr×nh ®Ó dμnh thªm dung l−îng trèng. Muèn kiÓm tra dung
l−îng cña ch−¬ng tr×nh , kÝch vμo Disk space ®Ó quan s¸t.
1.1.6 Sau khi kÝch Next, mμn h×nh hiÖn khung th«ng tin Setup vμ ch−¬ng tr×nh b¾t ®Çu cμi
®Æt.
1.1.7 Sau khi hoμn thμnh Setup, mμn h×nh hiÖn khung tho¹i Setup complete. KÝch Finish
®Ó kÕt thóc qui tr×nh cμi ®Æt.
Trang 2
[email protected] CM Soft 70 NCT F2 Q10
1.1.8 Sau khi kÕt thóc cμi ®Æt cöa sæ kÕ tiÕp lμ Internet explorer. KÝch Close ®Ó vÒ cöa sæ
nÒn windows, biÓu t−îng tù ®éng g¸n trong nÒn windows.
1.2. Khëi ®éng vμ tho¸t khái ch−¬ng tr×nh
1.2.1 Khëi ®éng MATLAB: KÝch vμo biÓu
t−îng Matlab trªn mμn h×nh Windows hoÆc
kÝch:
Start/ Programs/ Matlab/ Matlab 5.3, sau
®ã hiÖn cöa sæ Matlab vμ cöa sæ lÖnh (
MATLAB Command Widowd )
Còng nh− c¸c ch−¬ng tr×nh ch¹y trong m«i tr−êng windows kh¸c Matlab còng cã
nh÷ng thμnh phÇn giao diÖn cña ch−¬ng tr×nh.
• Dßng trªn cïng lμ thanh tiªu ®Ò gåm:
+ TËn tr¸i lμ biÓu t−îng ch−¬ng tr×nh. Khi kÝch vμo biÓu
t−îng nμy Matlab hiÖn menu xæ chøa c¸c lÖnh liªn quan tíi
viÖc xö lý khung cöa sæ ch−¬ng tr×nh còng nh− tho¸t ch−¬ng
tr×nh.
+ KÕ biÓu t−îng lμ tªn ch−¬ng tr×nh cïng víi cöa sæ ch−¬ng tr×nh. TËn cïng lμ ba
biÓu t−îng cã chøc n¨ng phong to, thu nhá vμ tho¸t ch−¬ng tr×nh.
• Dßng thø hai lμ thanh tr×nh ®¬n gåm c¸c menu chÝnh ch−a c¸c lÖnh liªn quan ®Õn
viÖc t¹o, xö lý, g¸n thuéc tÝnh vμo ®å thÞ.
• Dßng thø ba lμ thanh c«ng cô chøa c¸c hép c«ng cô chøa c¸c lÖnh nhanh.
• PhÇn chiÕm gÇn chän mμn h×nh lμ n¬i nhËp c¸c lÖnh .
• Cuèi cïng lμ thanh t¸c vô hiÖn th«ng tin vÒ t×nh tr¹ng ®ang xö lý, thùc hiÖn ®èi
t−îng.
1.2.2 Tho¸t khái MATLAB: Tõ cöa sæ Matlab command winwods ®¸nh lÖnh quit hoÆc
kÝch biÓu t−îng (X) n»m ngay gãc ph¶i thanh tiªu ®Ò Matlab.
Trang 3
[email protected] CM Soft 70 NCT F2 Q10
Ch−¬ng 2
C¸c lÖnh menu cña MATLAB
2.1 File
KÝch vμo File hiÖn menu chøa c¸c lÖnh liªn quan ®Õn viÖc t¹o
míi, qu¶n lý, g¸n thuéc tÝnh cho ®èi t−îng, ®å thÞ cò
2.1.1 New: HiÖn menu chøa ba lÖn ®Ó t¹o ®èi t−îng míi
• M-file: HiÖn cöa sæ
Editor/ Debugger. §©y lμ
m«i tr−êng b¹n t¹o vòng
nh− söa ®æi, gì rèi c¸c tËp
tin ch−¬ng tr×nh nh− M-
file hoÆc MEX-file hoÆc
c¸c ®èi t−îng nμo ®ã ®Ó
thùc hiÖn cho mét nhiÖn
Trang 4
[email protected] CM Soft 70 NCT F2 Q10
vô nμo ®ã. §Ó cã thÓ t¹o ®−îc c¸c tËp tin nμy, b¹n ph¶i dïng ng«n ng÷ lËp tr×nh C
hoÆc FORTRAN
• Figure: §©y lμ m«i tr−êng ®å ho¹ ®Ó b¹n tù vÏ c¸c ®èi t−îng hoÆc ®Ó Matlab vÏ
c¸c ®å thÞ theo hμm lÖnh b¹n nhËp tõ khung cöa sæ lÖnh cña MATLAB hoÆc hiÖn
®å thÞ cò víi lÖnh Open tõ lÖnh cña Matlab.
• Model : HiÖn cöa sæ
Simulink n»m bªn ph¶i vμ cöa
sæ th− viÖn n»m bªn tr¸i ®Ó
ban tù chän mμ thiÕt kÕ ®èi
t−îng.
2.1.2 Open : Më tËp tin ®å thÞ hoÆc h×nh ¶nh trong cöa sæ Figure ®Ó sö lý.
• Trong cöa sæ Matlab
Command windows kÝch:
File/ Open
• Tõ khung tho¹i, kÝch chän
tªn tËp tin cïng thÓ lo¹i
tËp tin ( nÕu cÇn thiÕt )
muèn më vμ kÝch OK.
2.1.3 Open Selection:
Cã chøc n¨ng më nh÷ng ®èi t−îng ®· chän riªng vμ l−u thμnh mét tËp tin riªng
Trang 5
[email protected] CM Soft 70 NCT F2 Q10
2.1.4 Run Script : Ch¹y mét ch−¬ng tr×nh ®· l−u thμnh
mét tËp tin.
• Trong cöa sæ Matlab Command windows kÝch:
File/ Run Script. Mμn h×nh xuÊt hiÖn khung tho¹i
Run Script.
• §¸nh tªn tËp tin vμo trong khung nhËp lÖnh Run vμ
kÝch OK. §Ó t×m kiÕm c¸c tËp tin
®· l−u kÝch vμo Browse.
• Tõ khung tho¹i chän nguån chøa
vμ tªn tËp tin cÇn më råi kÝch vμo
Open ®Ó vÒ khung tho¹i Run Script
víi toμn bé ®−êng dÉn cïng tªn tËp
tin muèn më. KÝch OK.
2.1.4 Load Workspace: T¶i tªn tËp tin . mat lªn mμn h×nh.
• Trong cöa sæ Matlab Command windows
kÝch: File/ Load Workspace. Mμn h×nh
xuÊt hiÖn khung tho¹i Load. Mat file.
• Tõ khung tho¹i, kÝch chän tªn tËp tin vμ kÝch
Open.
B¹n cã thÓ dïng lÖnh Load Workspace ®Ó
du nhËp c¸c tËp tin tõ nh÷ng nguån chøa kh¸c nhau
b»ng c¸ch còng trong khung tho¹i Load.mat file, kÝch vμo mòi tªn cña Look in ®Ó t×m
nguån chøa, chän tËp tin muèn du nhËp.
Ghi chó : §Ó cã thÓ du nhËp c¸c tËp tin kh¸c vμo MTALAB, ®iÒu b¾t buéc lμ c¸c tËp tin
®ã ph¶i ®−îc so¹n theo hÖ m· ASCII.
Trang 6
[email protected] CM Soft 70 NCT F2 Q10
2.1.5 Save Workspace as : L−u c¸c néi dung trong
cöa sæ lÖnh cña MATLAB theo workspace.
• NhËp c¸c lÖnh trong cöa sæ lÖnh cña
MATLAB nh− b×nh th−êng.
• Tõ menu File, chän Save Workspace as. Mμn
h×nh hiÖn khung tho¹i Save Workspace as.
• Tõ khung tho¹i, ®¸nh tªn tËp tin muèn l−u (
nÕu cÇn thiÕt chän nguån chøa víi Save in). KÝch Save. Nªn nhí lμ tËp tin nμy
®−îc l−u theo d¹ng *.mat vμ còng lμ tËp tin ch−¬ng tr×nh gèc.
2.1.5 Show Workspace :
HiÖn néi dung cña tËp tin workspace lªn mμn
h×nh b»ng c¸ch tõ menu File, chän Save Workspace.
Trong cöa sæ Save Workspace sÏ hiÖn néi dung cña
tËp tin Workspace ®· ghi.
2.1.6 Show graphics proprty editor :
HiÖn cöa sæ thay ®æi thuéc tÝnh mÆc
®Þnh cña h×nh ¶nh ®Ó ¸p dông cho tÊt c¶ c¸c
h×nh ¶nh sÏ ®−îc thiÕt kÕ sau nμy hoÆc chØ
riªng theo tõng nhiÖm vô.
• Tõ cöa sæ lÖnh MATLAB kÝch :
File/Show graphics proprty editor.
Mμn h×nh xuÊt hiÖn cöa sæ th«ng tin
Show graphics proprty editor.
• Chän thuéc tÝnh cÇn thay ®æi vμ kÝch
Close. ThÝ dô, thay ®æi thuéc tÝnh
Echo tõ mÆc ®Þnh lμ on thμnh off
Trang 7
[email protected] CM Soft 70 NCT F2 Q10
2.1.7 Show gui layout tool : HiÓn thÞ danh môc c¸c biÓu t−îng c«ng cô trong cöa sæ
Guide
Control
Panel.
Tõ cöa
sæ nμy, b¹n
cã thÓ thªm
c¸c biÓu
t−îng c«ng
cô tù t¹o
míi, g¸n biÓu t−îng biÓu t−îng c«ng cô vμo ®èi t−îng ®Ó thùc hiÖn nhanh b»ng c¸ch kÝch
chän biÓu t−îng trong khung biÓu t−îng Guide tool ®Ó hiÖn th«ng tin vμ tõ ®©y g¸n nh·n
còng nh− c¸c thuéc tÝnh kh¸c cho c«ng cô võa chän vμ kÝch Apply.
2.1.8 Set path: G¸n tËp tin vμo th− môc bÊt kú trong MATLAB.
• Tõ menu File
chän Set Path.
Mμn h×nh hiÖn
cöa sæ Path
Browser. Cöa sæ
nμy liÖt kª c¸c
danh s¸ch c¸c
®−êng dÉn th−
môc n»m trong khung cöa bªn tr¸i vμ danh s¸ch tªn tËp tin n»m bªn ph¶i.
• Chän tªn ®−êng dÉn th− môc trong khung cöa sæ Path vμ chän tªn tËp tin muèn
thay ®æi vÞ trÝ chøa trong khung cöa sæ hoÆc ®¸nh trùc tiÕp vμo khung nhËp lÖnh
Current Directory.
• Tõ menu File cña cöa sæ Path Browser, chän Save ®Ó l−u nh÷nh thay ®æi.
2.1.9 Preferenes: HiÖn khung tho¹i víi ba tuú chän ®Ó ng−êi sö dông cã thÓ thay ®æi
tham sè cho phï hîp nhiÖm vô.
Trang 8
[email protected] CM Soft 70 NCT F2 Q10
• General: HiÖn khung tho¹i ngay khi
kÝch chän lÖnh Preferences tõ menu
File cña khung cöa sæ lÖnh MATLAB.
KÝch chän c¸c lo¹i tham sè muèn thay
®æi hoÆc g¸n thªm vμ kÝch OK.
• Command Windows Font: HiÖn khung danh môc font cïng thuéc tÝnh ®Ó ng−êi
sö dông thay ®æi ph«ng mÆc ®Þnh thμnh font quªn thuéc.
• Copy Options : HiÖn khung th«ng tin ®Ó
ng−êi sö dông tuú chän nh÷ng thμnh phÇn ®èi
t−îng muèn sao chÐp.
2.1.9 Print Setup :
Trang 9
[email protected] CM Soft 70 NCT F2 Q10
HiÖn khung tho¹i Print cïng c¸c chøc n¨ng phô liªn quan ®Õn viÖc in.
2.1.10 Exit MATLAB ( Ctrl + Q ) :
Tho¸t khái ch−¬ng tr×nh MATLAB ®Ó vÒ l¹i mμn h×nh nÒn hÖ ®iÒu hμnh windows.
2.2 Edits : HiÖn menu con chøa c¸c lÖnh kiªn quan ®Õn viÖc xö lý
c¸c ®èi t−îng.
2.2.1 Undo : Huû lÖnh ®· thùc hiÖn tr−íc ®ã. Sau khi chän lÖnh
Undo ®Ó huû lÖnh, lÖnh sÏ ®æi thμnh Redo ®Ó ng−êi sö dông kh«i
phôc nh÷ng g× ®· huû tr−¬c ®ã víi lÖnh Undo.
2.2.2 Cut ( Ctrl +X) : C¾t ®èi t−îng trong khung cöa sæ lÖnh hiÖn
hμnhvμ sai ®ã d¸n vμo vÞ trÝ ®· chän. §èi t−îng ®−îc c¾t cã thÓ lμ
mét c«ng thøc, mét chuçi ký tù, hμm lÖnh . .
• Chän ®èi t−îng muèn c¾t trong khung cöa sæ lÖnh MATLAB vμ Ên Ctrl + C hoÆc
chon Copy tõ menu Edit. ThÝ dô c«ng thøc.
• Chän vÞ trÝ bÊt kú ®Ó d¸n t¹mvμ Ên Ctrl + V hoÆc Paste. §èi t−îng sÏ ®−îc d¸n
vμo vÞ trÝ ®· chän . Nªn nhí ®©y chØ lμ phiªn b¶n ®· ®−îc sao chÐp vμ d¸n.
• Chän vÞ trÝ thùc sù muèn d¸n. KÝch lªnh Cut hoÆc Ên Ctrl +X, ®èi t−îng võa d¸n
lóc tr−íc sÏ t¹m biÕn mÊt t¹i vÞ trÝ hiÖn hμnh. Ên Ctrl + V ®Ó d¸n.
2.2.3 Copy ( Ctrl + C) : Sao chÐp ®èi t−îng trong khung cöa sæ lÖnh MATLAB vμ sau ®ã
d¸n vμo vÞ trÝ ®· chän. §èi t−îng ®−îc copy cã thÓ lμ mét c«ng thøc, mét chuçi ký tù,
hμm lÖnh .
• Chän ®èi t−îng cÇn copy t¹i cöa sæ lÖnh MATLAB vμ sau ®ã Ên Ctrl + V hoÆc
Copy. ThÝ dô , chän c«nh thøc.
• Chän vÞ trÝ cÇn copy ®Õn vμ ©ns Ctrl + V.
2.2.4 Paste ( Ctrl + V) : D¸n ®èi t−îng ®· ®−îc Cut hoÆc Copy vμo vÞ trÝ ®· chän. Ngoμi
ra , b¹n cßn cã thÓ dïng lÖnh Paste cña MATLAB ®Ó d¸n c¸c ®èi t−îng kh¸c vμo
MATLAB.
2.2.5 Clear : Xo¸ ®èi t−îng ®· ®−îc chän trong khung cöa sæ MATLAB.
Trang 10
[email protected] CM Soft 70 NCT F2 Q10
2..2.6 Select Att : Chän toand bé néi dung trong khung cöa sæ lÖnh cña MATLAB.
2.2.7 Clear Session : Xo¸ toμn bé néi dung cña cöa sæ lÖnh MATLAB sau khi chän víi
lÖnh Select Att.
2.3 View : §¸nh dÊu chän hoÆc xo¸ ®Ó hiÖn hay giÊu thanh c«ng cô trong khung cöa sæ
lÖnh MATLAB.
2.4 Window : HiÖn th«ng tin vÒ sè tËp tin ®· thùc hiÖn trong khung cöa sæ lÖnh
MATLAB.
Ch−¬ng 3
C¸c kh¸i niÖm c¬ b¶n
MATLAB ( Matrix laboratory) lμ1 phÇn mÒm øng dông ch¹y trong trong m«i
tr−êng Windows. Dïng ®Ó m« pháng c¸c hÖ thèng ®iÒu khiÓn d−íi d¹ng hÖ ph−¬ng tr×nh
tr¹ng th¸i hoÆc s¬ ®å cÊu tróc. Thùc hiÖn c¸c phÐp to¸n vÒ xö lý ma trËn, xö lý tÝn hiÖu
còng nh− xö lý ®å ho¹
3.1 Mét sè phÝm chuyªn dông vμ lÖnh th«ng dông :
- ↑ hoÆc Ctrl + p : Gäi l¹i c¸c lÖnh ®· thùc hiÖn tr−íc ®ã.
- ↓ hoÆc Ctrl +n : Gäi l¹i lÖnh võa thùc hiÖn tr−íc ®ã.
- → hoÆc Ctrl + f : chuyÓn con trá sang bªn ph¶i 1 ký tù.
- ← hoÆc Ctrl + b: chuyÓn con trá sang tr¸i mét ký tù.
- DÊu (;) ®Ó kÕt thóc mét dßng lÖnh
- ↵ nh¶y xuèng dßng d−íi
- Ctrl + A hoÆc Home : chuyÓn con trá vÒ ®Çu dßng.
- Ctrl + E hoÆc End: ChuyÓn con trá ®Õn cuèi dßng.
- BackSpace: Xo¸ ký tù bªn tr¸i con trá.
- Esc: xo¸ dßng lÖnh.
Trang 11
[email protected] CM Soft 70 NCT F2 Q10
- Ctrl + K : Xo¸ tõ vÞ trÝ con trá ®Õn cuèi dßng.
- Ctrl + C : Dõng ch−¬ng tr×nh ®ang thùc hiÖn.
- Clc : lÖnh xo¸ mμn h×nh.
- Clf: LÖnh xo¸ mμn h×nh ®å ho¹.
- Input: lÖnh nhËp d÷ liÖu vμo tõ bμn phÝm.
- Demo: lÖnh cho phÐp xem c¸c ch−¬ng tr×nh mÉu.
- Help: lÖnh cho phÐp xem phÇn trî gióp.
- Ctrl – c: Dõng ch−¬ng tr×nh khi nã bÞ r¬i vμo tr¹ng th¸i lÆp kh«ng
kÕt thóc.
- Dßng lÖnh dμi: NÕu dßng lÖnh dμi qu¸ th× dïng ↵ ®Ó chuyÓn
xuèng dßng d−íi.
3.2 BiÕn trong MATLAB
3.2.1Tªn: Lμ mét d·y ký tù bao gåm c¸c ch÷ c¸i hay c¸c ch÷ sè hoÆc mét sè
ký tù ®Æc biÖt dïng ®Ó chØ tªn cña biÕn hoÆc tªn cña hμm. Chóng ph¶i ®−îc
b¾t ®Çu b»ng ch÷ c¸i sau ®ã cã thÓ lμ c¸c ch÷ sè hoÆc mét vμi ký tù ®Æc biÖt.
ChiÒu dμi tèi ®a cña tªn lμ 31 ký tù.
B×nh th−êng Matlab cã sù ph©n biÖt c¸c biÕn t¹o bëi ch÷ c¸i th−êng vμ
ch÷ c¸i hoa. C¸c lÖnh cña Matlab nãi chung th−êng sö dông ch÷ c¸i th−êng.
ViÖc ph©n biÖt ®ã cã thÓ ®−îc bá qua nÕu chóng ta thùc hiÖn lÖnh : >>
casensen off
3.2.2 Mét sè lÖnh víi biÕn:
- clear: lÖnh xo¸ tÊt c¶ c¸c biÕn ®· ®−îc ®Þnh nghÜa tr−íc trong
ch−¬ng tr×nh .
- clear biÕn1, biÕn 2... : xo¸ c¸c biÕn ®−îc liÖt kª trong c©u lÖnh.
- Who: hiÓn thÞ c¸c biÕn ®· ®−îc ®Þnh nghÜa trong ch−¬ng tr×nh.
Trang 12
[email protected] CM Soft 70 NCT F2 Q10
- Whos: hiÓn thÞ c¸c biÕn ®· ®−îc ®Þnh nghÜa trong ch−¬ng tr×nh
cïng víi c¸c th«ng sè vÒ biÕn.
- Size ( tªn biÕn ®· ®−îc ®Þnh nghÜa ): cho biÐt kÝch cì cña biÕn d−íi
d¹ng ma trËn víi phÇn tö thø nhÊt lμ sè hμng cña ma trËn, phÇn tö
thø 2 lμ sè cét cña ma trËn.
- Save: L−u gi÷ c¸c biÕn vμo mét File cã tªn lμ Matlab. mat.
- Load: T¶i c¸c biÕn ®· ®−îc l−u d÷ trong mét File ®−a vμo vïng lμm
viÖc.
3.2.3 Mét sè biÕn ®· ®−îc ®Þnh nghÜa tr−íc:
- ans: tù ®éng g¸n tªn nμy cho kÕt qu¶ cña mét phÐp tÝnh mμ ta
kh«ng ®Æt tªn.
VD >> [ 1 2] ↵
ans =
2
- pi
- realmax: ®−a ra gi¸ trÞ cña sè lín nhÊt mμ m¸y tÝnh cã thÓ tÝnh
to¸n ®−îc.
- realmin: ®−a ra gi¸ trÞ cña sè nhá nhÊt mμ m¸y tÝnh cã thÓ tÝnh
to¸n ®−îc.
3.3 C¸c phÐp to¸n v« h−íng:
3.3.1 BiÓu thøc sè häc:
Matlab cã thÓ thùc hiÖn tÊt c¶ c¸c biÓu thøc to¸n häc th«ng qua c¸c
to¸n tö sau:
+ phÐp céng.
- PhÐp trõ.
* phÐp nh©n.
Trang 13
[email protected] CM Soft 70 NCT F2 Q10
/ PhÐp chia ph¶i
\ phÐp chia tr¸i
^ phÐp luü thõa.
= ®Ó g¸n mét gi¸ trÞ nμo ®ã cho mét biÕn.
3.3.2 Thø tù −u tiªn trong phÐp to¸n sè häc:
ngoÆc ®¬n.
luü thõa
nh©n, chia.
Céng ,trõ.
3.3.3 C¸c vÝ dô:
• VÝ dô 1: Gi¶i ph−¬ng tr×nh bËc hai ax2 +bx +c = 0
Ta biÕt c¸c nghiÖm cña ph−¬ng tr×nh nμy cã d¹ng:
− b ± b 2 − 4ac
x =
2a
V× Matlab lμ mét ch−¬ng tr×nh tÝnh to¸n sè nªn chóng ta ph¶i x¸c ®Þnh
c¸c gi¸ trÞ a, b, c.
DÊu = ®−îc sö dông ®Ó g¸n gi¸ trÞ cña a, b, c nh− sau ( gâ phÝm Enter ë cuèi
mçi hμng)
>>a = 2
a=
2
>>b = 5;
>>c = -3;
⇒ DÊu (;) ë cuèi dßng th× Matlab sÏ kh«ng hiÓn thÞ l¹i gi¸ trÞ võa nhËp.
>> x1= (-b + sqrt(b^2- 4*a*c))/(2*a)
x1 =
0.5000
Trang 14
[email protected] CM Soft 70 NCT F2 Q10
>> x2= (-b - sqrt(b^2- 4*a*c))/(2*a)
x2 =
-3
• VÝ dô 2: TÝnh gi¸ trÞ cña ®a thø.
>> a = x^3 -2*x^2 - 6;
>>b = x^2 + 5*x -7;
>>x=3;
>> w = a/b
w=
0.1765
3.4 Sö dông lÖnh tõ c¸c file lÖnh:
Nh÷ng lÖnh cña Matlab cã thÓ ®−îc ®−a vμo mét file. Sau ®ã b¹n sÏ
h−íng dÉn Matlab lμm viÖc víi c¸c lÖnh ®ã. B©y giê, víi vÝ dô 1, chóng ta sÏ
®−a toμn bé c¸c lÖnh trªn vμo mét file lÊy tªn lμ vidu.m. Tªn cña file ph¶i
®−îc b¾t ®Çu b»ng mét ký tù vμ phÇn më réng lμ .m. C¸c b−íc nh− sau:
B−íc 1: File / New/ M-file, M«i tr−êng so¹n th¶o Editor / Debugger sÏ xuÊt
hiÖn
B−íc 2: Trªn mμn h×nh so¹n th¶o, ta gâ c¸c lÖnh cña Matlab.
a = 2;
b = 5;
c=-3;
x1= (-b + sqrt(b^2- 4*a*c))/(2*a)
x2= (-b - sqrt(b^2- 4*a*c))/(2*a)
B−íc 3: Ghi l¹i néi dung tËp tin víi tªn vidu.m råi tho¸t khái m«i tr−êng
so¹n th¶o ®Ó trë vÒ cöa sæ lÖnh ( Matlab Command window.
B−íc 4: T¹i cöa sæ lÖnh ta gâ tªn tËp tin.
>>vidu.m ↵
Trang 15
[email protected] CM Soft 70 NCT F2 Q10
• C¸c file *.m cã 2 lo¹i:
- Script file: lμ mét ch−¬ng tr×nh con kh«ng cã gi¸ trÞ tr¶ vÒ.
- Function file: lμ mét ch−¬ng tr×nh con cã gi¸ trÞ tr¶ vÒ.
3.5 Dßng nh¾c g¸n gi¸ trÞ c¸c biÕn:
§Ó thay ®æi c¸c gi¸ trÞ a,b,c ta ph¶i so¹n th¶o l¹i file vidu.m råi ch¹y
l¹i. Ta söa l¹i ch−¬ng tr×nh ®Ó cã dßng nh¾c nhËp a,.b,c víi c¸c lÇn ch¹y
ch−¬ng tr×nh kh¸c nhau.
B−íc 1: File / New/ M-file, M«i tr−êng so¹n th¶o Editor / Debugger sÏ xuÊt
hiÖn
B−íc 2:
a=input('nhap a= ');
b=input('nhap b= ');
c=input('nhap c= ');
x1= (-b + sqrt(b^2- 4*a*c))/(2*a)
x2= (-b - sqrt(b^2- 4*a*c))/(2*a)
B−íc 3: Ghi l¹i néi dung tËp tin víi tªn vidu.m.
Khi ch¹y ch−¬ng tr×nh ta thu ®−îc kÕt qu¶ nh− sau:
B−íc 4: T¹i cöa sæ lÖnh ta gâ tªn tËp tin.
>>vidu ↵
nhap vao a= 1
nhap vao b= -2
nhap vao c= 5
x1 =
1.0000 + 2.0000i
x2 =
1.0000 – 2.0000i
Hai nghiÖm øng víi c¸c gia trÞ a,b,c võa nhËp vμo vμ lμ nghiÖm ¶o.
VD vÒ Script file: Gi¶i bμi tËp m¹ch: cho m¹ch ®iÖn nh− h×nh vÏ
Trang 16
[email protected] CM Soft 70 NCT F2 Q10
R L C
u
H·y tÝnh dßng trong m¹ch vμ ®iÖn ¸p trªn tõng phÇn tö .
H·y vμo cöa sæ so¹n th¶o vμ trong cöa sæ nμy ta viÕt ch−¬ng tr×nh nh− sau:
R=input( 'nhap gia tri cho R = ')
C=input( 'nhap gia tri cho C = ')
L=input( 'nhap gia tri cho L = ')
U=input( 'nhap gia tri cho U = ')
ZL=2*50*pi*L*i
ZC=1/(2*50*pi*C*i)
Z=R+ZL+ZC
i= U/Z
UR=i*R
UL=i*ZL
UC=i*ZC
Sau khi viÕt xong ch−¬ng tr×nh ta kÝch vμo biÓu t−îng save trong cöa sæ so¹n th¶o vμ tªn
lμ vd1.
Muèn ch¹y ta trë l¹i cöa sæ MATLAB command Window vμ tõ dÊu nh¾c lÖnh:
>> vd1
nhap gia tri cho R = 1000
R= 1000
nhap gia tri cho C = 0.1
C= 0.1000
nhap gia tri cho L = 0.2
L = 0.2000
nhap gia tri cho U = 220
U = 220
ZL = 0 +62.8319i
ZC = 0 - 0.0318i
Trang 17
[email protected] CM Soft 70 NCT F2 Q10
Z = 1.0000e+003 +6.2800e+001i
i = 0.2191 - 0.0138i
UR = 2.1914e+002 -1.3762e+001i
UL = 0.8647 +13.7687i
UC = -0.0004 - 0.0070i
3.6 C¸ch t¹o mét hμm:
Mçi mét file hμm cña Matlab ( M - file) ®Òu ®−îc khai b¸o nh− sau:
Function [ Tªn kÕt qu¶] = tªn hμm (danh s¸ch c¸c biÕn).
PhÇn th©n cña ch−¬ng tr×nh trong hμm lμ c¸c lÖnh cña Matlab thùc hiÖn
viÖc tÝnh to¸n gi¸ trÞ cña ®¹i l−îng ®−îc nªu trong phÇn tªn kÕt qu¶ theo c¸c
biÕn ®−îc nªu tronhg phÇn danh s¸ch biÕn. C¸c biÕn chØ cã t¸c dông néi
trong hμm võa ®−îc khai b¸o Tªn cña c¸c biÕn d−îc c¸ch nhau b»nh dÊu
phÈy (,).
VÝ dô ta thμnh lËp hμm ®æi tõ ®é sang radian:
function rad = change(do)
rad = do*pi/180; % doi do sang radian
Trong Matlab c¸c dßng ghi chó sau dÊu % kh«ng cã t¸c dông thùc thi,
chóng ®¬n gi¶n lμ nh÷ng dßng nh¾c ®Ó ng−êi ®äc ch−¬ng tr×nh dÔ hiÓu mμ
th«i.
File.m th−êng lÊy tªn lμ tªn cña hμm, ta ®Æt tªn file hμm võa lËp lμ change.m.
NÕu muèn ®æi 450 sang radian, chØ cÇn gâ:
>>rad = change(45)
rad =
0.7854
VÝ du: t¹o hμm gi¶i ph−¬ng tr×nh bËc hai, tªn tËp tin ®−îc ®Æt lμ bachai.m.
function [x1,x2] = bachai(a,b,c)
delta = b^2-4*a*c;
x1 = (-b + sqrt(delta))/(2*a);
Trang 18
[email protected] CM Soft 70 NCT F2 Q10
x2 = (-b - sqrt(delta))/(2*a);
>>[x1,x2]=bachai1(4,6,-7)
x1 =
0.77707
x2 =
-2.2707
3.7 VÏ c¸c hμm .
§å ThÞ BËc Hai
Khi muèn vÏ mét hμm nμo 200
®ã, ph¶i x¸c ®Þnh hμm ®ã 150
trong mét file.m sau ®ã sö
100
dông lÖnh :
Fplot(‘tªn hμm’,[khoange vÏ]) 50
2
VÝ dô vÏ hμm y = 4x +6x-7
0
trong ®o¹n [-6, 6], ta lËp file
bachai1.m. -50
-6 -4 -2 0 2 4 6
function y = bachai1(x)
a = 4;b = 6; c = -7;
y =a*x^2 + 6*x + c;
>>fplot(‘bachai1’, [-6,6])
3.8 L−u vμ lÊy d÷ liÖu:
Ta cã thÓ cã thÓ t¹o lËp mét file d÷ liÖu, sau ®ã khi cÇn dïng th× lÊy ra. VÝ
dô t¹o lËp mét ma trËn A:
1 1 .1 1 .2
2 2 .1 2 .2
3 3 .2 3 .2
Sau ®ã ta l−u ma trËn vμo mét file cã tªn lμ dulieu1.
>>A = [1 1.1 1.2;2 2.1 2.2;3 3.2 3.2]
A=
Trang 19