【正文】 以下代碼，z=peaks。subplot(2,2,1)contour(z)subplot(2,2,2)[c,h]=contour(z,[,,])。clabel(c,h)。subplot(2,2,3)[c,h]=contour(z,6)。clabel(c,h)。subplot(2,2,4)[c,h]=contourf(z,8)。clabel(c,h)。函數(shù)M文件需要接受參數(shù)的輸入和輸出，只處理輸入?yún)?shù)傳遞的數(shù)據(jù)，并把處理結(jié)果作為函數(shù)輸出參數(shù)返回給指定接收變量。matlab提供的很多庫(kù)函數(shù)就是這種文件。一般要求文件名和函數(shù)名一致，才能調(diào)用該文件，否則會(huì)很麻煩。例2：編寫(xiě)函數(shù)實(shí)現(xiàn)求給定矩陣的平均值。%功能：求一個(gè)矩陣所有元素的平均值function average=junzhi(A)s=sum(sum(A))。%求所有元素的和[m,n]=size(A)。%求一共有多少個(gè)元素average=s/m/n。%求平均值二、程序流程控制順序語(yǔ)句見(jiàn)例1條件語(yǔ)句（1）if語(yǔ)句（a）格式1if 邏輯表達(dá)式 執(zhí)行語(yǔ)句end例3：根據(jù)嘴型判斷雞鴨。function ji_ya(s)if s==39。尖嘴39。 disp(39。您描述的是雞雞雞！39。)。endif s==39。扁嘴39。 disp(39。您描述的是鴨鴨鴨!39。)。end 在命令行窗口輸入： s=39。扁嘴39。 ji_ya(s)（b）格式2if 邏輯表達(dá)式 執(zhí)行語(yǔ)句1else 執(zhí)行語(yǔ)句1end例4：根據(jù)嘴型判斷是否是雞。function ji(s)if s==39。尖嘴39。 disp(39。您描述的可能是雞雞雞！39。)。else disp(39。您描述的肯定不是雞!39。)。end在命令行窗口輸入： s=39。扁嘴39。 ji_ya(s)（c）格式3if 邏輯表達(dá)式1 執(zhí)行語(yǔ)句1elseif邏輯表達(dá)式2 執(zhí)行語(yǔ)句2elseif邏輯表達(dá)式3 執(zhí)行語(yǔ)句3…else 執(zhí)行語(yǔ)句end例5：編寫(xiě)M文件求函數(shù)的函數(shù)值function f=hanshu1(x)if x0 f=x^2。elseif 0=xamp。x2 f=23*x。elseif 2=xamp。x4 f=sin(2*x)。else f=exp(2*x)。end（d）if語(yǔ)句嵌套結(jié)構(gòu)請(qǐng)同學(xué)們考慮用嵌套結(jié)構(gòu)編寫(xiě)例5（2）switch語(yǔ)句格式：switch 表達(dá)式case1 表達(dá)式1 語(yǔ)句1case2 表達(dá)式2 語(yǔ)句2…otherwise 語(yǔ)句end例6：function ff(method)switch method case {39。linear39。} disp(39。Method is linear39。) case 39。cubic39。 disp(39。Method is cubic39。) case 39。nearest39。 disp(39。Method is nearest39。) otherwise disp(39。Unknown method.39。)end在matlab命令窗口輸入： method=39。linear39。 ff(method)例7：為了促銷(xiāo)某種商品，銷(xiāo)量在1000以下，按原價(jià)收款，1000至2000按原價(jià)的90%收款，以后每增加1000，少收10%，到5折以后不再增加折扣，編寫(xiě)程序計(jì)算應(yīng)收貨款。function m=money(p,q)%p價(jià)格，q銷(xiāo)售量switch q case q0amp。q1000 money=p*q。 case q=1000amp。q2000 money=p*q*。 case q=2000amp。q3000 money=p*q*。 case q=3000amp。q4000 money=p*q*。 case q=4000amp。q5000 money=p*q*。 otherwise money=p*q*。end循環(huán)語(yǔ)句（1）for循環(huán)格式：for 表達(dá)式 循環(huán)體end嵌套格式：for 表達(dá)式 循環(huán)體 for 表達(dá)式 循環(huán)體 …. end ….end例8：編寫(xiě)程序求小于自然數(shù)n的能被m整除的自然數(shù)。function t=zc(n,m)if mn disp(39。輸入?yún)?shù)m過(guò)大，請(qǐng)重新輸入！39。)else t=[]。 for i=1:n if rem(i,m)==0 t=[t,i]。 end endend（2）while語(yǔ)句格式：while 邏輯表達(dá)式 循環(huán)體end例9：一球從100米高度自由落下，每次落地后又反彈回原高度的一半再次落下，是編寫(xiě)程序?qū)崿F(xiàn)：（1）輸入次數(shù)，求出后一次反彈的高度和經(jīng)過(guò)的總路程；（2）求經(jīng)過(guò)的總路程。function [h1,s1,s2,nn]=zylt(H,n)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%h1最后一次彈起的高度；%s1落下n次后經(jīng)過(guò)的路程%s2經(jīng)過(guò)的總路程%nn總共彈起的次數(shù)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%計(jì)算落下n次后彈起的高度和經(jīng)過(guò)的路程s1=0。h=H。for i=1:n s1=s1+h。 h=h/2。 s1=s1+h。endh1=h。%計(jì)算總路程s2=0。m=0。hh=H。while hheps s2=s2+hh。 hh=hh/2。 s2=s2+hh。 m=m+1。endnn=m。交互語(yǔ)句（1）輸入輸出控制、等待用戶(hù)響應(yīng)例10：輸入線性方程組的系數(shù)矩陣和常數(shù)向量，并求解方程組。：A=input(39。請(qǐng)輸入線性方程組的系數(shù)矩陣：39。)。b=input(39。請(qǐng)輸入常數(shù)向量：39。)。x=A\b。disp(39。按回車(chē)鍵查看結(jié)果！39。)pause%熟悉pause的用法disp(39。您所求解的方程組的解是：39。)x輸入：A=[1 2 2。3 2 1。2 4 3],b=[3。2。6],觀察運(yùn)行結(jié)果。（2）程序中斷、錯(cuò)誤與警告例11：驗(yàn)證輸入的數(shù)是否是素?cái)?shù).function sushu(m)%驗(yàn)證輸入的m是否是素?cái)?shù)if nargin~=1 error(39。輸入?yún)?shù)個(gè)數(shù)應(yīng)為1！39。)endif m=0|mround(m)~=0 warning(39。數(shù)據(jù)輸入錯(cuò)誤，輸入?yún)?shù)應(yīng)是自然數(shù)！39。) returnendk=sqrt(m)。i=2。while i if rem(m,i)==0 disp(39。輸入的數(shù)不是素?cái)?shù)！39。) break else i=i+1。 end if ik disp(39。輸入的是素?cái)?shù)!39。) return endend三、應(yīng)用舉例第六講：Matlab的GUI設(shè)計(jì)