matlab教程2012a习题解答张志涌.docVIP

M文件和句柄函数 习题6及解答 请分别写出用for 和while 循环语句计算的程序。此外，还请写出避免循环的数值、符号计算程序。 〖解答〗 （1）for 环 tic s1=0; for k=0:1e6 s1=s1+0.2^(k); end; s1 toc s1 = 1.2500 Elapsed time is 1.453482 seconds. （2）while环 tic s2=1;k=1; while k1e6+1 %注意：上限与for环不同 s2=s2+0.2^k; k=k+1; end s2 toc s2 = 1.2500 Elapsed time is 2.716870 seconds. （3）数值求和指令 tic s3=sum(0.2.^(0:1e6)) toc s3 = 1.2500 Elapsed time is 0.626723 seconds. （4）符号求和指令 tic syms k; s4=vpa(symsum(0.2^k,0,1e6)) toc s4 = 1.2500000000000000000000000000000 Elapsed time is 4.029501 seconds. 〖说明〗 使用数值求和指令和“数组运算”的计算速度最快。 符号计算的精度最高，但速度慢。 以上程序运行时间仅供参考。具体时间与所用机器、那程序是否初次运行、在MATLAB指令窗还是在M-book中运行等因素有关。 编写一个函数M文件，它的功能：没有输入量时，画出单位圆（见图p6.2-1）；输入量是大于2的自然数N时，绘制正N边形，图名应反映显示多边形的真实边数（见图p6.2-2）；输入量是“非自然数”时，给出“出错提示”。此外，函数M文件应有H1行、帮助说明和程序编写人姓名。 〖解答〗 （1）函数M文件 function prob_solve602(n) % prob_solve602(n) plot a circle or a polygon with n edges % prob_solve602 plot a circle % n 应为大于2的自然数 % By ZZY, 2006-2-15 if nargin==0 t=0:pi/100:2*pi; x=exp(i*t); str=Circle; else if (nargin~=0)(n=2) error(输入量应是大于2的自然数) end; if n-round(n)~=0 %检查非自然数 error(输入量应是大于2的自然数) end; t=(0:n)/n*2*pi; x=exp(i*t); str=[Polygon with , int2str(n), edges]; % 合成字符串 end plot(real(x),imag(x),r,LineWidth,4) title(str) axis square image off shg （2）各典型运行情况 prob_solve602 图 p6.2-1 prob_solve602(6) 图 p6.2-2 prob_solve602(2) ??? Error using == prob_solve602 输入量应是大于2的自然数 prob_solve602(7.3) ??? Error using == prob_solve602 输入量应是大于2的自然数 用泛函指令fminbnd求在x=0附近的极小值。fminbnd的第一个输入量要求使用匿名函数表达。 〖解答〗 （1）在0附近的较小区间 [-0.5, 0.5]内搜索 y=@(x)-exp(-x).*abs(sin(cos(x))); [x1,y1]=fminbnd(y,-0.5,0.5) x1 = -0.4999 y1 = -1.2681 （2）在包含0的更大区间内搜索 在 [-0.5, 0.5] 区间搜索到的极小值点在该搜索区间的边界上，所以搜索区间应该进一步扩大。 [x2,y2]=fminbnd(y,-1,1) x2 = -0.8634 y2 = -1.4348 〖说明〗 在[-1, 1]区间的搜索结果是区间的内点。进一步扩大搜索区间，只可能找到更远离0的极小值点，所以没有必要进一步搜索。 在matlab的 \toolbox\matlab\elmat\private文件夹上有一个“烟圈矩阵”发生函数smoke.m。运行指令smoke