06-程序流程控制教程.pptVIP

MATLAB 程序设计 MATLAB程序设计入门篇程序流程控制 6-1 循环指令 MATLAB 提供两种循环指令 for 循环（For Loop） while 循环（While Loop） Formats of For Loops Format 1: for 变数 = 向量 表达式 end 在上述语法中，变量的值会被依次设定为向量的每一个元素值，来执行介于 for和 end 之间的表达式。 Format 2: for 变数 = 矩阵 表达式 end 在上述语法中，变量的值会被依次设定为矩阵的每一个列，来执行介于 for 和 end 之间的表达式。 下列 for 循环会产生一个长度为 6 的调和数列（Harmonic Sequence）： x = zeros(1,6); % 变量 x 是一个 1×6 大小的零矩阵 for i = 1:6 x(i) = 1/i; end x % 显示 x x = 1.0000 0.5000 0.3333 0.2500 0.2000 0.1667 在上例中，矩阵 x 最初是一个 1×6 大小的零矩阵，在 for 循环中，变量 i 的值依次是 1 到 6，因此矩阵 x 的第 i 个元素的值依次被设为 1/i。 向量化编程一句可实现 x=1./[1:6] 我们接着可用分数形式来显示此数列： format rat % 使用分数形式来显式数值 disp(x) 1 1/2 1/3 1/4 1/5 1/6 程序流程控制之范例一 for 循环可以是多层或巢状式（Nested）的，在下例中即产生一个 6×6 的Hilbert 矩阵 h，其中为于第 i 列、第 j 行的元素为 ： h = zeros(6); % 变数 x 是一个 6×6 大小的零矩阵 for i = 1:6 for j = 1:6 h(i,j) = 1/(i+j-1); end end format rat % 使用分数形式来显式所有数值 h % 显示 h 1 1/2 1/3 1/4 1/5 1/6 1/2 1/3 1/4 1/5 1/6 1/7 1/3 1/4 1/5 1/6 1/7 1/8 1/4 1/5 1/6 1/7 1/8 1/9 1/5 1/6 1/7 1/8 1/9 1/10 1/6 1/7 1/8 1/9 1/10 1/11 向量化编程： [I,J]=meshgrid(1:6,1:6);h=1./[I+J-1] 程序流程控制之范例二 在下例中，for 循环列出先前产生的 Hilbert 矩阵的每一直行的平方和： format short % 回到预设形式来显式所有数值 for i = h disp(sum(i.^2)); % 返回h每一列的平方和 end 1.4914 0.5118 0.2774 0.1787 0.1262 0.0944 在上例中，由于 h 是一个矩阵，因此每一次 i 的值就是矩阵 h 的一列的内容。 向量化编程：sum(h.^2) 程序流程控制之范例三 若要跳出 for 循环，可用 break 指令。例如，若要找出最小的 n 值，满足 n! 1e100，可输入如下： for i = 1:1000 if prod(1:i) 1e100 fprintf(%g! = %e 1e100\n, i, prod(1:i)); break; % 跳出 for 循环 end end 70! = 1.197857e+100 1e100 向量化编程： find(arrayfun(@(x) prod(1:x),1:1000)1e100,1) 程序流程控制之范例