1、第二講 MATLAB程序設(shè)計1 MATLAB有兩種工作方式：交互式的命令行工作方式M文件的程序工作方式2一. M文件3一. M文件用MATLAB語言編寫的程序，稱為M文件。 M文件有兩類：命令文件和函數(shù)文件。命令文件：沒有輸入?yún)?shù)，也不返回輸出參數(shù)。函數(shù)文件：可以輸入?yún)?shù)，也可返回輸出參數(shù)。4M文件的建立與編輯建立新的M文件： 從MATLAB命令窗口的File菜單中選擇New菜單項，再選擇M-file命令。一. M文件編輯已有的M文件： 從MATLAB命令窗口的Flie菜單中選擇Open M-file命令。5命令文件【例1】 建立一個命令文件將變量a,b的值互換。 一. M文件 e31m文件：

2、 a=1:9; b=11,12,13;14,15,16;17,18,19; c=a;a=b;b=c; a b在MATLAB的命令窗口中輸入e31，將會執(zhí)行該命令文件。7二數(shù)據(jù)的輸入輸出8【例2】 求一元二次方程a2 +bx+c=0的根。二數(shù)據(jù)的輸入輸出 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)將該程序以aa.m文件存盤，然后運行aa.m文件。102 pause函數(shù)：暫停程序的執(zhí)行。調(diào)用格式： pause(延遲秒數(shù))注：如果省略延遲時間，直接使用paus

3、e，則將暫停程序，直到用戶按任一鍵后程序繼續(xù)執(zhí)行。二數(shù)據(jù)的輸入輸出11二數(shù)據(jù)的輸入輸出 3 disp函數(shù)：命令窗口輸出函數(shù)。調(diào)用格式： disp(輸出項)注：輸出項為字符串或矩陣。例如 A=Hello,MATLAB; disp(A)輸出為： Hello,MATLAB12（一）選擇結(jié)構(gòu)選擇結(jié)構(gòu)的語句有if語句和switch語句。1 if語句 格式一： if 條件 語句組 end三程序設(shè)計格式二： if 條件 語句組1 else 語句組2 end14(一)選擇結(jié)構(gòu)格式三： if 條件1 語句組1 elseif 條件2 語句組2 elseif 條件m 語句組m else 語句組m+1 end三程序設(shè)

4、計15【例5】 輸入一個字符，若為大寫字母，則輸出其后繼字符，若為小寫字母，則輸出其前導(dǎo)字符，若為其他字符則原樣輸出。c=input(,s); if c=A & c=a& c=z disp(setstr(abs(c)-1); else disp(c); end三程序設(shè)計172 switch語句switch語句根據(jù)變量或表達(dá)式的取值不同，分別執(zhí)行不同的語句。其格式為： switch 表達(dá)式 case 值1 語句組1 case 值2語句組2 case 值m 語句組m otherwise 語句組m+1 end三程序設(shè)計(一)選擇結(jié)構(gòu)18【例6】 根據(jù)變量 num 的值來決定顯示的內(nèi)容。 num=in

5、put(請輸入一個數(shù)); switch num case -1 disp(I am a teacher.); case 0 disp(I am a student.); case 1 disp(You are a teacher.); otherwise disp(You are a student.); end19（二）循環(huán)結(jié)構(gòu)實現(xiàn)循環(huán)結(jié)構(gòu)的語句：for語句和while語句。1 for語句：格式： for 循環(huán)變量=表達(dá)式1:表達(dá)式2:表達(dá)式3 循環(huán)體語句 end注：其中表達(dá)式1的值為循環(huán)變量的初值，表達(dá)式2的值為步長，表達(dá)式3的值為循環(huán)變量的終值。步長為1時，表達(dá)式2可以省略。三程序設(shè)計2

6、02 while語句格式： while (條件) 循環(huán)體語句 end三程序設(shè)計(二)循環(huán)結(jié)構(gòu)213 循環(huán)的嵌套如果一個循環(huán)結(jié)構(gòu)的循環(huán)體又包括一個循環(huán)結(jié)構(gòu)，就稱為循環(huán)的嵌套，或稱為多重循環(huán)結(jié)構(gòu)。多重循環(huán)的嵌套層數(shù)可以是任意的。可以按照嵌套層數(shù)，分別叫做二重循環(huán)、三重循環(huán)等。處于內(nèi)部的循環(huán)叫作內(nèi)循環(huán)，處于外部的循環(huán)叫作外循環(huán)。三程序設(shè)計(二)循環(huán)結(jié)構(gòu)22四函數(shù)文件24函數(shù)文件是另一種形式的M文件，每一個函數(shù)文件都定義一個函數(shù)。事實上，MATLAB提供的標(biāo)準(zhǔn)函數(shù)大部分都是由函數(shù)文件定義的。四. 函數(shù)文件25【例9】 編寫函數(shù)文件求小于任意自然數(shù)n的Fibonacci數(shù)列各項。 function f=

7、ffib(n) %用于求Fibonacci數(shù)列的函數(shù)文件 %f=ffib(n) %1999年9月30日編 f=1,1; i=1; while f(i)+f(i+1)n f(i+2)=f(i)+f(i+1); i=i+1; end將以上函數(shù)文件以文件名ffib.m存盤，然后在MATLAB命令窗口輸入以下命令，可求小于2 000的Fibonacci數(shù)。 ffib(2000)272函數(shù)調(diào)用 函數(shù)文件編制好后，就可調(diào)用函數(shù)進行計算了。如上面定義ffib函數(shù)后，調(diào)用它求小于2000的Fibonacci數(shù)。函數(shù)調(diào)用的一般格式是： 輸出實參表=函數(shù)名(輸入實參表)四. 函數(shù)文件28【例10】 利用函數(shù)文件，

8、實現(xiàn)直角坐標(biāo)(x,y)與極坐標(biāo)(,)之間的轉(zhuǎn)換。函數(shù)文件tran.m： function gama,theta=tran(x,y) gama=sqrt(x*x+y*y); theta=atan(y/x);調(diào)用tran.m的命令文件main1.m： x=input(Please input x=:); y=input(Please input y=:); gama,theta=tran(x,y); gama theta29【例11】 利用函數(shù)的遞歸調(diào)用，求n！。 function f=factor(n) if n=1 f=1; else f=factor(n-1)*n; end return;

9、%返回 四. 函數(shù)文件在命令文件main2.m中調(diào)用函數(shù)文件factor.m： for i=1:10 fac(i)=factor(i); end fac303 函數(shù)所傳遞參數(shù)的可調(diào)性MATLAB在函數(shù)調(diào)用上有一個與眾不同之處：函數(shù)所傳遞參數(shù)數(shù)目的可調(diào)性。憑借這一點，一個函數(shù)可完成多種功能。 在調(diào)用函數(shù)時，MATLAB用兩個永久變量nargin和nargout分別記錄調(diào)用該函數(shù)時的輸入實參和輸出實參的個數(shù)。只要在函數(shù)文件中包含這兩個變量，就可以準(zhǔn)確地知道該函數(shù)文件被調(diào)用時的輸入輸出參數(shù)個數(shù)，從而決定函數(shù)如何進行處理。四. 函數(shù)文件31【例12】 nargin用法示例函數(shù)文件examp.m： fu

10、nction fout=charray(a,b,c) if nargin=1 fout=a; elseif nargin=2 fout=a+b; elseif nargin=3 fout=(a*b*c)/2; end命令文件mydemo.m： x=1:3;y=1;2;3; examp(x) examp(x,y) examp(x,y,3)執(zhí)行mydemo.m輸出是： ans = 1 2 3 ans = 2 4 6 ans = 2132五全局變量和局部變量33 在MATLAB中，全局變量用命令global定義。函數(shù)文件的內(nèi)部變量是局部的，與其他函數(shù)文件及MATLAB工作空間相互隔離。但是，如果在若

11、干函數(shù)中，都把某一變量定義為全局變量，那么這些函數(shù)將公用這一個變量。全局變量的作用域是整個MATLAB工作空間，即全程有效。所有的函數(shù)都可以對它進行存取和修改。因此，定義全局變量是函數(shù)間傳遞信息的一種手段。四. 函數(shù)文件34四. 函數(shù)文件【例13】 全局變量應(yīng)用示例先建立函數(shù)文件wadd.m，該函數(shù)將輸入的參數(shù)加權(quán)相加。function f=wadd(x,y) %add two variable global ALPHA BETA f=ALPHA*x+BETA*y;在命令窗口中輸入： global ALPHA BETA ALPHA=1; BETA=2; s=wadd(1,2)輸出為： s = 535習(xí) 題1. 輸入20個數(shù)，求其中最大數(shù)和最小數(shù)。要求分別用循環(huán)結(jié)構(gòu)和調(diào)用MATLAB的max函