




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)
文檔簡介
1、MATLAB編程介紹,MATLAB中各種命令可以完成許多單一的任務(wù),對于某些較為復(fù)雜的問題,僅靠現(xiàn)有的命令或函數(shù)來解決,往往是難以直接達到目的。為此,要運用MATLAB編程語言編制程序,形成M-文件。MATLAB中每一個命令都是一個M-文件。 MATLAB程序(M-文件)有兩種形式,一種是可直接運行的命令文件,另一種是可供調(diào)用的函數(shù)文件,這兩種文件的擴展名相同,均為“.m”,故稱為M-文件。,MATLAB編程介紹,MATLAB程序編輯是在編輯窗口(MATLAB Editor/Debugger)中進行,程序運行結(jié)果或錯誤 信息顯示在命令窗口(Command Windows)中。,1 、編寫M-文
2、件的步驟,(1)打開程序編輯窗口: 點擊MATLAB命令窗口工具條 上的New File圖標,或點擊命令窗口中的File菜單,點 擊選項中的New,選擇M-file, 就可打開MATLAB程序編 輯窗口(MATLAB Editor/Debugger),用戶可在窗口中 編寫、修改程序。,(2)保存編寫的程序:在程序編寫過程中,或程序編寫完成后,點擊工具條上的Save圖標,或按F2鍵,再選擇當前工作的路徑,輸入自己定義的文件名,點擊保存保存文件,保存后文件名后自動加上“.m”后綴,即為M-文件。, 文件名要以字母開頭,中間不含空格和標點符號,且不超過8個 字符;不能用數(shù)字和專用變量名作為文件名,如
3、1,pi,ans,eps等。, 程序通常保存在文件夾:MatlabWork中,此路徑為MATLAB 軟件默認路徑。, 程序若要存放到其它文件夾中,則要添加該路徑到MATLAB的路徑設(shè)置中,否則程序不能連接執(zhí)行。,點擊編輯窗工具條上File菜單,選擇Set Path.,此時打開路徑設(shè)置窗口,再點擊Add Foider.,選擇你所使用的文件夾(若沒有則要先新建),選定后該路徑自動顯示在列表的第一行或最后一行,再點擊 Save 保存即可。,具體操作:,2、 運行M-文件,運行一個已經(jīng)編寫好且存盤的M-文件有兩種方法。,在MATLAB編輯窗口中對正在編輯的M-文件,保存后可直接點擊工具條上的Debug
4、中的Run命令運行或按F5鍵,即可運行當前程序,其運行結(jié)果或信息顯示在命令窗口,此時,用戶可以切換到命令窗口進行查看,之后可回到編輯窗對程序進行修改、重新保存再次運行。,方法一:,在MATLAB命令窗口中直接輸入文件名回車,即可運行相應(yīng)的M-文件。其運行結(jié)果或程序運行中的錯誤信息等顯示在命令窗口,此時,用戶若要修改程序則可以切換到編輯窗口進行修改,修改后重新保存。,方法二:,當對程序進行修改保存時,當前程序?qū)⒏采w原來程序,此時若要想保存為另外的文件,則要點擊工具條上File菜單,選擇Save as.,另取文件名保存即可。,注意:,例. 分別建立命令文件和函數(shù)文件,將華氏溫度f轉(zhuǎn)換 為攝氏溫度c
5、。,首先建立命令文件并以文件名function1.m存盤。 clear; %清除工作空間中的變量 f=input(Input Fahrenheit temperature:); c=5*(f-32)/9 然后在MATLAB的命令窗口中輸入function1,將會執(zhí)行該命令文件,執(zhí)行情況為: Input Fahrenheit temperature:73 c = 22.7778,程序1:,首先建立函數(shù)文件function2.m: function c=function2(f) c=5*(f-32)/9 然后在MATLAB的命令窗口調(diào)用該函數(shù)文件: clear; y=input(Input Fah
6、renheit temperature:); x=function2(y) 輸出情況為: Input Fahrenheit temperature:70 c = 21.1111 x = 21.1111,程序2:,例: 建立同時計算,,,的函數(shù)。, function y1 , y2=fun1(a , b , n), y1=(a+b).n ;, y2=(a-b).n;,編寫完成后用函數(shù)名(fun1)作為文件名存盤, 形成 一個M-文件fun1.m可以進行調(diào)用。,MATLAB程序的基本結(jié)構(gòu),MATLAB提供了三種基本的程序結(jié)構(gòu):順序結(jié)構(gòu)、 循環(huán)結(jié)構(gòu)和選擇結(jié)構(gòu)。一個MATLAB程序就是使計算機 實現(xiàn)某
7、一類功能任務(wù)的命令集。用戶編制程序時,需要 按照解決問題的邏輯步驟和求解算法,用MATLAB的語 句和函數(shù),按照三種基本結(jié)構(gòu)逐條編寫,一般一行一條 語句,便于調(diào)試和查找錯誤。,其中功能處理部分是核心,主要是根據(jù)處理問題的 邏輯步驟,采用循環(huán)迭加、判斷推理、分步計算等 來解決問題。,程序一般包含三個部分:,數(shù)據(jù)輸入、功能處理、結(jié)果輸出,程序編寫完成后,需要存盤,形成永久性文件, 之后可以隨時對它進行調(diào)用或修改。,M-文件中的輸入、輸出方式,1、 輸入方式,s=Any string!; %輸入字符串,(1)程序中直接賦值,x=2; %輸入單個值,a=1,2,3;4,5,6;7,8,0; %輸入矩陣
8、(a為3階方陣),(2)程序運行時賦值,變量 = input( 提示符號串 );,格式:,功能:,從鍵盤接收數(shù)據(jù)并賦予等號左邊變量。其中 提示符號串 是在命令窗口中的提示。,y=input(Please input y=); %提示語句輸入,如:,x=input( x=);,當執(zhí)行到此語句時,暫停程序運行,在命令窗口中光標將停留在提示符號串后閃爍,等待從鍵盤接收數(shù)據(jù)。用戶從鍵盤輸入數(shù)據(jù)并按回車鍵確認,繼續(xù)程序運行。其優(yōu)點是每次運行可以輸入不同數(shù)據(jù),實現(xiàn)人機對話。,如果在input函數(shù)調(diào)用時采用s選項,則允許用戶輸入一個 字符串。例如,想輸入一個人的姓名,可采用命令:,xm=input(What
9、s your name?,s);,2、 輸出方式,(1)直接輸出,x %變量后無分號“;”,變量直接輸出,(2)格式控制輸出(fprintf命令),fprintf(x=%.3f y=%.0fn , a,b); %按指定格式輸出,格式:,其中,單引號內(nèi)為格式說明符:%前為提示符,%.3f為數(shù)據(jù)格式符, 通常與單引號后面的數(shù)據(jù)相對應(yīng), f指顯示一實數(shù),.3表示顯示的數(shù)小數(shù) 點后3位,小數(shù)點前整數(shù)部分不限。 %.0f指顯示一實數(shù)其小數(shù)部分顯示 0位,結(jié)果只顯示整數(shù)部分。n為換行命令符,即一行顯示結(jié)束后,光標 自動換到下一行。a,b為要顯示的變量,它們與前面的格式符相對應(yīng)。,disp(a); % a可
10、以是矩陣、數(shù)組和單個值,例:,fprintf(x=%.5f y=%.0fn, pi, sqrt(2);,結(jié)果為: x=3.14159 y=1,x=input(Input x please.); y=input(Input y please.); z=x; x=y; y=z; disp(x); disp(y);,例:輸入x,y的值,并將它們的值互換后輸出。,a=input(a=); b=input(b=); c=input(c=); d=b*b-4*a*c; x=(-b+sqrt(d)/(2*a),(-b-sqrt(d)/(2*a); disp(x1=,num2str(x(1),x2=,num2
11、str(x(2);,例:求一元二次方程ax2 +bx+c=0的根。,a=1 b=2 c=3 x1=-1+1.4142i,x2=-1-1.4142i,結(jié)果為:,num2str:實數(shù)轉(zhuǎn)換為字符串型,MATLAB循環(huán)結(jié)構(gòu)與應(yīng)用,循環(huán)結(jié)構(gòu)是MATLAB程序中最常用的結(jié)構(gòu)之一,它 是讓計算機按照一定的條件多次重復(fù)執(zhí)行某一命令集, 從而實現(xiàn)諸如累加、迭代、分層計算等功能。MATLAB 軟件提供了兩種循環(huán)結(jié)構(gòu):,有限次循環(huán)for-end結(jié)構(gòu)和條件循環(huán)while-end結(jié)構(gòu)。,1、有限次循環(huán)(for-end)結(jié)構(gòu),for n=n1 : step : n2 commands-1 end commands-2,格
12、式:,其中n為循環(huán)變量(自己定義),n1為起始值,step為步長,n2為結(jié)束值;commands-1為循環(huán)體(循環(huán)模塊),commands-2為后序命令 。,循環(huán)變量n從n1開始,執(zhí)行commands-1,遇到end時n自動增加 step步長,同時與n2比較,當n不超過n2時重復(fù)執(zhí)行commands-1;當 n超過n2時轉(zhuǎn)向執(zhí)行commands-2。當步長為1時,格式中step可以省略, 即為for n=n1 : n2格式。在設(shè)計for-end循環(huán)結(jié)構(gòu)時,應(yīng)在循環(huán)模塊內(nèi)充 分利用循環(huán)變量的變化規(guī)律,開展相關(guān)的運算。另外,結(jié)構(gòu)中for與end 要成對出現(xiàn)。,作用:,例:求正整數(shù)n的階乘:p=12
13、 3 n = n!,,并求出n=20時的結(jié)果。,分析:,這是一個逐次乘積的過程,從1開始存放在積p中, 再乘以下一個正整數(shù)之后再取代積p,依次類推,直到 乘至第n個數(shù)。因此,使用循環(huán)結(jié)構(gòu)來實現(xiàn)。,clear;clc; n=20; %賦值給定正整數(shù) p=1; %設(shè)定存放階乘的變量p并賦初值1 for i=1:n %定義循環(huán)變量i從1到n,以1為步長,即連續(xù)正整數(shù) p=p*i; %先計算右端乘積后再賦給左端的變量p fprintf(i=%.0f, p=%.0fn,i,p) %逐行顯示 end %循環(huán)結(jié)構(gòu)結(jié)束,問題:,是否可以考慮利用input命令對n進行賦值,隨時改變其大小。 如果可以,請修改上述
14、程序并運行。,2、 條件循環(huán)(while-end)結(jié)構(gòu),格式:,while (conditions) commands-1; end commands-2;,當條件 conditions為真時,執(zhí)行commands-1,當遇 到end時,自動檢測條件;當條件 conditions不滿足時,轉(zhuǎn) 向執(zhí)行commands-2。在設(shè)計while-end循環(huán)結(jié)構(gòu)時,應(yīng)在 確保在執(zhí)行了一定次數(shù)之后可以結(jié)束循環(huán);否則,就成了 “死循環(huán)”,即無限次重復(fù)執(zhí)行循環(huán)。另外,結(jié)構(gòu)中while與 end要成對出現(xiàn)。,作用:,關(guān)系運算,MATLAB提供了6種關(guān)系運算符:,(小于),=(小于或等于),(大于),=(大于或等
15、于),=(等于),=(不等于),關(guān)系、邏輯運算,關(guān)系運算符的運算法則為:, 當兩個比較量是標量時,直接比較兩數(shù)的大小。若關(guān)系成立,關(guān)系表達式結(jié)果為1,否則為0。, 當參與比較的量是兩個維數(shù)相同的矩陣時,比較是對兩矩陣相同位置的元素按標量關(guān)系運算規(guī)則逐個進行,并給出元素比較結(jié)果。最終的關(guān)系運算的結(jié)果是一個維數(shù)與原矩陣相同的矩陣,它的元素由0或1組成。, 當參與比較的一個是標量,而另一個是矩陣時,則把標量 與矩陣的每一個元素按標量關(guān)系運算規(guī)則逐個比較,并給出元 素比較結(jié)果。最終的關(guān)系運算的結(jié)果是一個維數(shù)與原矩陣相同 的矩陣,它的元素由0或1組成。, 設(shè)參與邏輯運算的是兩個標量a和b,那么,邏輯運算
16、,MATLAB提供了3種邏輯運算符:當a非零時,運算結(jié)果為0。, 若參與邏輯運算的是兩個同維矩陣,那么運算將對矩陣相同位置上的元素按標量規(guī)則逐個進行。最終運算結(jié)果是一個與原矩陣同維的矩陣,其元素由1或0組成。, 若參與邏輯運算的一個是標量,一個是矩陣,那么運算將在標量與矩陣中的每個元素之間按標量規(guī)則逐個進行。最終運算結(jié)果是一個與矩陣同維的矩陣,其元素由1或0組成。, 邏輯非是單目運算符,也服從矩陣運算規(guī)則。, 在算術(shù)、關(guān)系、邏輯運算中,算術(shù)運算優(yōu)先級最高,邏輯 運算優(yōu)先級最低。,例:根據(jù)e1+1+1/2!+1/3!+1/n! 求e的近似值,要求精確到小數(shù)點后8位。,分析:,這里不能給定n的大小
17、,但注意到第n次近似值和第n-1次近似值的差:,?,F(xiàn)要精確到小數(shù)點后8位,則只需要,因此,可以以此為條件運用條件循環(huán)來實現(xiàn)。,,,注意:在MATLAB中,用科學(xué)記數(shù)法1.0e-8來表示。,clear;clc; p=1; %設(shè)定存放階乘的變量p并賦初值1 s=1; %設(shè)定存放累加和的變量s并賦初值1 r=1; %設(shè)定前后兩次近似值的誤差r并賦初值1 k=0; %設(shè)定構(gòu)造連續(xù)正整數(shù)的變量k賦初值0又為循環(huán)次數(shù) while r=1.0e-8 %根據(jù)近似值的精度r大小來定義條件循環(huán) k=k+1; %轉(zhuǎn)換下一個正整數(shù),或開始統(tǒng)計循環(huán)次數(shù) p=p*k; %計算k的階乘p r=1/p; %計算前后兩次近似值
18、的誤差r s=s+r; %計算e的近似值s fprintf(k=%.0f, s=%.8fn,k,s) %逐行顯示出e的近似值 end %循環(huán)結(jié)構(gòu)結(jié)束,MATLAB選擇結(jié)構(gòu),1、單項選擇判斷(if-end)結(jié)構(gòu),格式:,if (condition) commands-1; end commands-2;,作用:,若條件condition成立,則執(zhí)行commands-1, 再順序執(zhí)行commands-2;否則,跳過commands-1, 直接執(zhí)行commands-2.,2 、多項選擇判斷(if-else-end)結(jié)構(gòu),格式:,if (condition1) commands-1; elseif (
19、condition2) commands-2; else commands-3; end commands-4;,作用 :,若條件condition1成立,則執(zhí)行commands-1,再轉(zhuǎn)向end,順序執(zhí)行 后續(xù)的命令commands-4;否則判斷條件condition2是否成立,若成立,則執(zhí)行commands-2,再轉(zhuǎn)向end執(zhí)行后續(xù)的命令;若條件ondition1和條件condition2均不成立,則執(zhí)行命令集commands-3,再順序向下執(zhí)行。,x=input(請輸入x的值:); if x=0 y= (x+sqrt(pi)/exp(2); else y=log(x+sqrt(1+x*x
20、)/2; end y,例: 計算分段函數(shù)的值。,例:編寫一個函數(shù)將百分制成績轉(zhuǎn)換為優(yōu)(A),良(B),中(C),差(D)四等級.,分析:,按照通常的等級劃分,一般90-100分為優(yōu),78-89分為良,60-77分為中,60分以下的為差,因此,可以用多項選擇判斷結(jié)構(gòu)來實現(xiàn)。,function jb=grade(fs) %定義函數(shù)grade,輸入變量為分數(shù)fs,返回值為等級A,B,C,D if fs=90 %判斷分數(shù)fs是否處在優(yōu)秀級別上 jb= A ; %定義為A級 elseif fs=78 %判斷分數(shù)fs是否處在良好級別上 jb= B ; %定義為B級 elseif fs=60 %判斷分數(shù)fs
21、是否處在合格級別上 jb= C ; %定義為C級 else %分數(shù)fs不處于以上任何級別上 jb= D ; %定義為D級 end %選擇結(jié)構(gòu)結(jié)束,3、switch語句,switch語句根據(jù)表達式的取值不同,分別執(zhí)行不同的 語句,其語句格式為:,switch 表達式 case 表達式1 語句組1 case 表達式2 語句組2 case 表達式m 語句組m otherwise 語句組n end,當表達式的值等于表達式1的值時,執(zhí)行語句組1,當表達式的值等于表達式2的值時,執(zhí)行語句組2,當表達式的值等于表達式m的值時,執(zhí)行語句組m,當表達式的值不等于case所列的表達式的值時,執(zhí)行語句組n。當任意一
22、個分支的語句執(zhí)行完后,直接執(zhí)行switch語句的下一句.,clear n=input(input n=); switch mod(n,5) case 0 fprintf(%d 是5的倍數(shù),n) otherwise fprintf(%d 不是5的倍數(shù),n) end,input n=46 46 不是5的倍數(shù),例: 判斷任意給定的整數(shù)n是否為5的倍數(shù),并輸出判斷結(jié)果。,MATLAB程序流程控制,在編寫程序過程中,根據(jù)問題求解的要求,往往需要在一定條件下跳出當前循環(huán),或終止程序運行,或暫停程序運行等,這也就是在編制程序過程中對程序流程進行控制。MATLAB軟件提供了相應(yīng)流程控制語句:break、return和pause等。,1、 break語句,break 語句導(dǎo)致包含break 指令的最內(nèi)層while或for循環(huán)的終止。通常是根據(jù)循環(huán)內(nèi)部另設(shè)的某種條件是否滿足來決定是否跳出循環(huán),因此,常和if判斷一起使用。
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 二零二五年度水上樂園游泳館場地租賃與水上樂園配套設(shè)施租賃協(xié)議
- 2025年度老舊小區(qū)外墻改造工程安全責(zé)任合同
- 二零二五年度國際貿(mào)易信用證業(yè)務(wù)代理及風(fēng)險管理協(xié)議
- 海洋漁業(yè)資源保護與海產(chǎn)品銷售一體化合同
- 二零二五年度企業(yè)用工協(xié)議與勞動權(quán)益保障與員工激勵機制合同
- 二零二五年度廠房裝修施工安全責(zé)任與綠色施工標準協(xié)議書
- 2025年度酒店與旅游紀念品店合作經(jīng)營合同
- 二零二五年度籃球活動參與者免責(zé)責(zé)任協(xié)議
- 二零二五年度汽車美容店員工勞動爭議解決合同模板
- 二零二五年度農(nóng)村房屋贈與合同附農(nóng)業(yè)保險合作協(xié)議
- 經(jīng)濟學(xué)彈性理論課件
- 公司員工獎懲制度流程
- 星巴克案例分析-星巴克成功之道
- 把未來點亮歌詞打印版
- 危險化學(xué)品建設(shè)項目竣工驗收報告
- 國家中醫(yī)藥管理局第3批24個專業(yè)104個病種中醫(yī)診療方案
- 婦產(chǎn)科學(xué)(第9版)第三章 女性生殖系統(tǒng)生理
- LY/T 2241-2014森林生態(tài)系統(tǒng)生物多樣性監(jiān)測與評估規(guī)范
- GB/T 9086-2007用于色度和光度測量的標準白板
- 2023年山東力明科技職業(yè)學(xué)院高職單招(數(shù)學(xué))試題庫含答案解析
- GB/T 24338.4-2018軌道交通電磁兼容第3-2部分:機車車輛設(shè)備
評論
0/150
提交評論