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

PAGEPAGE1M文件和句柄函數(shù)習(xí)題6及解答請分別寫出用for和while循環(huán)語句計算的程序。此外，還請寫出避免循環(huán)的數(shù)值、符號計算程序。〖解答〗（1）for環(huán)tics1=0;fork=0:1e6s1=s1+0.2^(k);end;s1tocs1=1.2500Elapsedtimeis1.453482seconds.（2）while環(huán)tics2=1;k=1;whilek<1e6+1 %注意：上限與for環(huán)不同 s2=s2+0.2^k; k=k+1;ends2tocs2=1.2500Elapsedtimeis2.716870seconds.（3）數(shù)值求和指令tics3=sum(0.2.^(0:1e6))tocs3=1.2500Elapsedtimeis0.626723seconds.（4）符號求和指令ticsymsk;s4=vpa(symsum(0.2^k,0,1e6))tocs4=1.2500000000000000000000000000000Elapsedtimeis4.029501seconds.〖說明〗使用數(shù)值求和指令和“數(shù)組運(yùn)算”的計算速度最快。符號計算的精度最高，但速度慢。以上程序運(yùn)行時間僅供參考。具體時間與所用機(jī)器、那程序是否初次運(yùn)行、在MATLAB指令窗還是在M-book中運(yùn)行等因素有關(guān)。編寫一個函數(shù)M文件，它的功能：沒有輸入量時，畫出單位圓（見圖p6.2-1）；輸入量是大于2的自然數(shù)N時，繪制正N邊形，圖名應(yīng)反映顯示多邊形的真實邊數(shù)（見圖p6.2-2）；輸入量是“非自然數(shù)”時，給出“出錯提示”。此外，函數(shù)M文件應(yīng)有H1行、幫助說明和程序編寫人姓名?！冀獯稹剑?）函數(shù)M文件functionprob_solve602(n)%prob_solve602(n) plotacircleorapolygonwithnedges%prob_solve602 plotacircle%n 應(yīng)為大于2的自然數(shù)%ByZZY,2006-2-15ifnargin==0t=0:pi/100:2*pi;x=exp(i*t);str='Circle';elseif(nargin~=0)&(n<=2)error('輸入量應(yīng)是大于2的自然數(shù)')end;ifn-round(n)~=0 %檢查非自然數(shù)error('輸入量應(yīng)是大于2的自然數(shù)')end;t=(0:n)/n*2*pi;x=exp(i*t);str=['Polygonwith',int2str(n),'edges']; %合成字符串endplot(real(x),imag(x),'r','LineWidth',4)title(str)axissquareimageoffshg（2）各典型運(yùn)行情況prob_solve602ans=C:\MATLAB2006\R2006a\work'smoke'notfound.（2）在當(dāng)前目錄下，進(jìn)行直接調(diào)用smoke.m的嘗試B=smoke(3,0,'dou