matlab課件編程mbook2008a第6章文件和函數句柄

1、第 6M 文件和函數句的解算及仿真就需要借助較復雜的 MATLAB 程序進行。這種 MATLAB 程序由較多的 MATLAB MATLAB 表達式組成，并循著一定的執(zhí)行次序運行。這種程m。這就是腳本M M 文件。第 6M 文件和函數句的解算及仿真就需要借助較復雜的 MATLAB 程序進行。這種 MATLAB 程序由較多的 MATLAB MATLAB 表達式組成，并循著一定的執(zhí)行次序運行。這種程m。這就是腳本M M 文件。作為一種計算機編程語言， MATLAB 提供了多種用于程序流控制的描述關鍵詞 MATLAB給定的描述關鍵詞，對這四種指令進行簡要的說明。if-else-end 條件控 6.1-

2、 x 1x 1xy1 x xfunction% Functioncalculateofexample6.1-for k=1:nif x(k)=1 elseif expr給出“邏輯 expr 1”該使用方法常被switch-casexyxyswitch-case 控制結6.1-2 switch-case6.1-2 for%forswitchA2,k case 100case acase bcase cswitch expr case value_1case value_2 case value_k AA AA for 循環(huán)和while 循6.1-6.1-3】創(chuàng)建Hilbert1i j（1）Hil

3、bert矩陣a(i, j) A=zeros(K,K); for m = 1:K。forn=A(m,n)=1/(m+n-Ashort1form=forn=A1(m,n)=1/(m+n-t1 =for whilefor ix=array whileA2=zeros(K,K); for m = 1:Kforn=A2(m,n)=t2 =A2=zeros(K,K); for m = 1:Kforn=A2(m,n)=t2 =t3 N116.1-4S nNkn1 kkkfunction% % % % CalculatethesumofaspecialseriesS=1+1/(1+2)+1/(1+2+N) S

4、um of a special seriesTheminimumamongallnumberstohave1/sum(1:N)epsilon S =N = syms k n N SINF =2SN = 1 到nsyms k n N SINF =2SN = 1 到n2 n(n2 2% % % n while1n=input(4n if mod(n,4)=0disp(int2str(n),階魔方矩陣的標稱和是 if NS=Ns012階魔方矩陣的標稱和是經驗證，A是魔方矩陣。控制程序流的其它常用指6.14 M 12階魔方矩陣的標稱和是經驗證，A是魔方矩陣?？刂瞥绦蛄鞯钠渌Ｓ弥?.14 M 算就顯得

5、煩瑣。而此時腳本文件最為適宜。“腳本”本身反映這樣一個事實：MATLAB只M 1.8 節(jié)。它只是一串按用戶意圖排列而成的（包括控制流向指令在內的）MATLAB指令集合。clear MATLAB 指令窗不關閉，這些變量將一直保存在基本工作空間中?；究臻g隨 MATLAB 的啟動而產生；只有當關閉 MATLAB時，該基本空間才被刪除。M 與腳本文件不同，函數文件（Function file）猶如一個“黑箱”。從外界只看到：傳給局部變量和全局變global指令，MATLAB也允許幾個不同的函數空間以及基本工作空間共享同一個使 用 說 明MATLAB。message 是提示用的字符串。keyboard

6、 時，將“控制權”交給鍵盤，用戶可以從鍵盤輸入MATLAB 指令。breakwhileforifend, switch - case, try catch 6.1-5。6.1-5。n 秒后，再繼續(xù)執(zhí)行。M 函數文件的一般結6.2-1M M 函數文件的一般結6.2-1M ：（A）根據指定的半徑，畫出藍色圓function% Theareaandperimeterofaregularpolygon正多邊形的面積和周長) The number of sidesTheAlinespecificationtodeterminelinetype/color The area of the regular

7、polygonTheperimeteroftheregular11N% % % % % % RNstrRN邊形給出多邊形面積S % ZhangZhiyong2006-1-switchnargin case 0casecase 2case 3;if nargout=0elseifnargout2axisequalsquare box onS = L = 6.2-1M 觀察，MATLAB M M M MATLABFunction又被細分為：主函數；子函數；嵌套函數；私用函數；匿名函主S = L = 6.2-1M 觀察，MATLAB M M M MATLABFunction又被細分為：主函數；子函數

8、；嵌套函數；私用函數；匿名函主函數（Primary 子函數6.3-1】編寫一個內含子函數的MfunctionHr=exm060301(flag% Demoforhandlesofprimary functionsand%linecircle。 function% %linecirclet產生的橫坐標變量tswitchwd case 10-plot(t,x,b,t,y,r,LineWidth,plot(t,x,b,t,y,r,LineWidth,2) case circleplot(x,y,-g,LineWidth,8), axis square offerror(linecircle) HH

9、=6.3-1綠色圓周t=0:2*pi/5:2*pi;x=cos(t);y=sin(t)6.3-2 匿名函數（Anonymous FH=(arglist)feval(FH,arglist)函數句柄（Function handle）MATLABfeval(FH,arglist)函數句柄（Function handle）MATLAB函數句柄的創(chuàng)建和觀hm =isa(hm,ans = ans =1CC =M1 =594273681M2 =594273681函數句柄的基本用6.3-36.4-1 節(jié)中，已經簡單介紹和演示了函數句柄兩種調用格式的簡單形argout1,argout2,argoutn=FunN

10、ame(argin1,函數句柄的基本用6.3-36.4-1 節(jié)中，已經簡單介紹和演示了函數句柄兩種調用格式的簡單形argout1,argout2,argoutn=FunName(argin1,argin2,argout1, argout2, , argoutn=Hfun(argin1, argin2, , arginn) argout1,argout2,argoutn=feval(Hfun,argin1,argin2,arginn)6.4-26.2-1建立的函數文件exm060201.mS =L6.4-1 直接調用函數exm060201所畫的紅色等邊三角Hexm =exm060201notfo

11、und.?Undefinedcommand/functionS =L = 10-6.4-2 在新視野中里用函數句柄Hexm所畫的紅色等邊三角6.4-36.2-1建立的函數文件exm060201.m6.4-2 在新視野中里用函數句柄Hexm所畫的紅色等邊三角6.4-36.2-1建立的函數文件exm060201.mHexm_2ans =ans =1?Undefinedfor while K 0.2i 10.20.220.21000000s1 = 是大于 2 的自然數N 繪制正 N 應反映顯示多邊形的真實邊（見圖 p6-2 10-p6-p6-function% % p6-p6-function% % % plotacircleorapolygonwithnedges plot a circle%ByZZY,2006-2-if nargin=0 ififn-str=Polygonwith,int2str(n),edges; Polygonwith6axissquareimageoff fminbnd y(x) axissquareimageoff fminbnd