MATLAB教程2012a-习题解答-第六章-张志涌-北航

PAGEPAGE1M文件和句柄函数习题6及解答请分别写出用for和while循环语句计算的程序。此外，还请写出避免循环的数值、符号计算程序。〖目的〗练习for环和while环的使用方法。避免循环，提高计算速度。计时指令tic,toc的用法。〖解答〗（1）for环tics1=0;fork=0:1e6s1=s1+0.2^(k);end;s1tocs1=1.2500Elapsedtimeis1.453482seconds.（2）while环tics2=1;k=1;whilek<1e6+1 %注意：上限与for环不同 s2=s2+0.2^k; k=k+1;ends2tocs2=1.2500Elapsedtimeis2.716870seconds.（3）数值求和指令tics3=sum(0.2.^(0:1e6))tocs3=1.2500Elapsedtimeis0.626723seconds.（4）符号求和指令ticsymsk;s4=vpa(symsum(0.2^k,0,1e6))tocs4=1.2500000000000000000000000000000Elapsedtimeis4.029501seconds.〖说明〗使用数值求和指令和“数组运算”的计算速度最快。符号计算的精度最高，但速度慢。以上程序运行时间仅供参考。具体时间与所用机器、那程序是否初次运行、在MATLAB指令窗还是在M-book中运行等因素有关。编写一个函数M文件，它的功能：没有输入量时，画出单位圆（见图p6.2-1）；输入量是大于2的自然数N时，绘制正N边形，图名应反映显示多边形的真实边数（见图p6.2-2）；输入量是“非自然数”时，给出“出错提示”。此外，函数M文件应有H1行、帮助说明和程序编写人姓名。〖目的〗较完整形式函数M文件的编写。H1的书写，帮助区内容的组织。指令nargin的用法。〖解答〗（1）函数M文件functionprob_solve602(n)%prob_solve602(n) plotacircleorapolygonwithnedges%prob_solve602 plotacircle%n 应为大于2的自然数%ByZZY,2006-2-15ifnargin==0t=0:pi/100:2*pi;x=exp(i*t);str='Circle';elseif(nargin~=0)&(n<=2)error('输入量应是大于2的自然数')end;ifn-round(n)~=0 %检查非自然数error('输入量应是大于2的自然数')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）各典型运行情况A=SM(3,0,'double') %产生一个3阶“伪特征根”矩阵disp('')pwd %显示当前所在目录，以证实符合题意。which('smoke') %再次检查在当前目录下能否看到smoke.m，起证实作用。disp('')b=C:\MATLAB2006\R2006a\work'smoke'notfound.A=-0.5000+0.8660i1.000000-0.5000-0.8660i1.00001.000001.0000ans=C:\MATLAB2006\R2006a\work'smoke'notfound.（2）在当前目录下，进行直接调用smok