MATLAB教程a第章習(xí)題解答張志涌

MATLAB教程ａ第章習(xí)題解答-張志涌————————————————————————————————作者:————————————————————————————————日期: M文件和句柄函數(shù)習(xí)題6及解答請分別寫出用ｆoｒ和while循環(huán)語句計算的程序。此外，還請寫出避免循環(huán)的數(shù)值、符號計算程序?！冀獯稹?1)fｏr環(huán)tics1＝0；fｏrk=0:1e６ｓ1=ｓ1+0．2＾(k)；end;ｓ１tocs1＝1.2500Elａｐseｄtｉmeｉｓ1.453482ｓeconds．(２)wｈilｅ環(huán)tics２=１;k=1;ｗhilｅk<1e6+1 ??％注意:上限與for環(huán)不同 ｓ2=ｓ2+０．2＾k;?k＝ｋ＋１;eｎｄs２tｏcs2=１．2500Elapsedtiｍeis2.71６８70sｅconｄs.（３）數(shù)值求和指令tics3=suｍ(0.2.＾(０:1e６)）ｔocs３＝１.2500Ｅｌａｐsedtiｍeis0.62６723seconds.(4）符號求和指令ｔicｓymsｋ;s4=vpａ(sｙmsum(0.2^k,０,1e６)）tocs4=1.250０0000000０0０000０00000000０0000Eｌaｐsedｔｉmeiｓ4.029501ｓecｏnds．〖說明〗使用數(shù)值求和指令和“數(shù)組運(yùn)算”的計算速度最快。符號計算的精度最高,但速度慢。以上程序運(yùn)行時間僅供參考。具體時間與所用機(jī)器、那程序是否初次運(yùn)行、在MAＴLAB指令窗還是在M-book中運(yùn)行等因素有關(guān)。編寫一個函數(shù)M文件，它的功能：沒有輸入量時，畫出單位圓（見圖p6.2-1)；輸入量是大于２的自然數(shù)N時,繪制正N邊形,圖名應(yīng)反映顯示多邊形的真實(shí)邊數(shù)(見圖p6.2-２）；輸入量是“非自然數(shù)”時，給出“出錯提示”。此外，函數(shù)M文件應(yīng)有H1行、幫助說明和程序編寫人姓名。〖解答〗（1）函數(shù)M文件funcｔionｐrob_solve6０2(n)%prob_solvｅ６02(n） ｐloｔaｃircleorａpolｙgonwｉthnedges%prob_sｏlｖe60２?ｐｌotａciｒclｅ%ｎ ? 應(yīng)為大于2的自然數(shù)%ByZZY,2006-2-15ifnargｉn＝=0t＝0：pi／100:2＊ｐi;x=ｅxp(i*ｔ);str=＇Cｉrcle';elseif(naｒgｉn～=0）＆(n<=2)erｒor('輸入量應(yīng)是大于２的自然數(shù)')end;ifn-rｏund(ｎ)~=0 ?? ?%檢查非自然數(shù)ｅrror（'輸入量應(yīng)是大于2的自然數(shù)')eｎd；ｔ=(0：n）/n*2*ｐｉ;x=exｐ(i*t);str＝［'Ｐolyｇｏnｗitｈ＇，iｎｔ2ｓtr（n），'ｅdges＇];? %合成字符串enｄplｏt(real(x),imａg(x),'r＇,'ＬineWｉdｔh＇,４）title（str)aｘissqｕareimageoｆfshｇ（2)各典型運(yùn)行情況prob_solｖe602圖p６.2-１prob＿solｖe602（６)圖p6.２-2pｒｏb＿solｖe6０2（２)?？？Erｒorusing==>prob_ｓolve602輸入量應(yīng)是大于2的自然數(shù)proｂ_solvｅ60２(7．3)？？?Eｒｒoruｓiｎg==>prｏb_solve602輸入量應(yīng)是大于2的自然數(shù)用泛函指令fminbnd求在x=０附近的極小值。fｍiｎbnd的第一個輸入量要求使用匿名函數(shù)表達(dá)?！冀獯稹?1）在0附近的較小區(qū)間［-0．５,0.5]內(nèi)搜索y=@(ｘ）-ｅxｐ（-x)．*abs(sin(cos（x)));[ｘ１,y1］=ｆｍinｂnｄ(y，-0.5,0．5)x１＝-0.4999y1＝-1.２６81(２）在包含0的更大區(qū)間內(nèi)搜索 在［-0．5,０．5］區(qū)間搜索到的極小值點(diǎn)在該搜索區(qū)間的邊界上,所以搜索區(qū)間應(yīng)該進(jìn)一步擴(kuò)大。[x2,y２]=fmiｎbnd（ｙ,－１,1)?x2=－0．８63４y2=－1.4348〖說明〗在［-1,1］區(qū)間的搜索結(jié)果是區(qū)間的內(nèi)點(diǎn)。進(jìn)一步擴(kuò)大搜索區(qū)間,只可能找到更遠(yuǎn)離0的極小值點(diǎn),所以沒有必要進(jìn)一步搜索。在matlab的\toｏlbｏx\ｍatlaｂ＼elｍａt(yī)\ｐriｖate文件夾上有一個“煙圈矩陣”發(fā)生函數(shù)smｏkｅ.ｍ。運(yùn)行指令smoke(3,0,＇double'),將生成一個3階偽特征根矩陣如下A=－0.50０0＋０．８660i1.0００0０0-0.5000-０.8660i1.00001.0０００01.0000現(xiàn)在的問題是：在ＭＡTLAB當(dāng)前目錄為＼work情況下,如何利用函數(shù)句柄調(diào)用smｏｋe．m函數(shù),產(chǎn)生3階偽特征根矩陣。請寫出相應(yīng)的程序或操作步驟。ｃleａrb=ｐwd? ? ?? ％獲取當(dāng)前目錄名字符串ｗｈich（'smoｋｅ') ?? %檢查在當(dāng)前目錄下能否看到sｍokｅ.mb_d=b;b_ｄ(ｅnｄ-4:ｅｎｄ)=［］; ? ％在b字符串中去除最后的四個字符,即worｋ。ｓtｒ=［b＿d,'\tooｌbox＼matｌab\elｍat＼privａt(yī)e'];?％構(gòu)造smokｅ所在的目錄字符串ｃd(sｔr)?? ? ％把smｏｋe．ｍ所在目錄設(shè)置為當(dāng)前目錄。SM＝@ｓmoke;? ?? %創(chuàng)建smokｅ．m的函數(shù)句柄cd(b） ?? ? ?%使worｋ恢復(fù)為當(dāng)前目錄diｓp(＇'）?? ? %為顯示清楚起見，產(chǎn)生一個顯示空行。A=SM(３,０，'doｕｂle')???? ％產(chǎn)生一個3階“偽特征根”矩陣ｄｉsp(＇'）pwｄ ??????%顯示當(dāng)前所在目錄，以證實(shí)符合題意。wｈｉcｈ('smoｋe＇) ? ?%再次檢查在當(dāng)前目錄下能否看到ｓｍoke.m,起證實(shí)作用。dｉsp（'')ｂ=C：\MAＴLAB2006\R2０06ａ\wｏrk'ｓｍokｅ'ｎotｆounｄ.Ａ=-0.5０00+０.８6６0i1.000００0-０.５00０－0.8660ｉ１.０0００1.０0000１．０00０anｓ=Ｃ:\M