MATLAB 循环：for、while

在 MATLAB 中，循环用于重复执行一段代码。最常用的两种循环是 for 和 while。掌握它们的写法和适用场景，能让你高效处理数组、迭代计算或条件驱动的任务。

一、for 循环：已知次数的重复操作

for 循环适用于事先知道要执行多少次的情况，比如遍历一个数组的所有元素。

基本语法

for 变量 = 起始值:步长:结束值
    % 循环体（要重复执行的代码）
end

• 变量：每次循环自动更新。
• 起始值:步长:结束值：定义循环范围。若省略步长，默认为 1。

操作步骤

1. 创建一个新脚本文件（例如命名为 example_for.m）。
2. 输入以下代码：

   for i = 1:5
       disp(i)
   end

3. 运行脚本，控制台将依次显示数字 1 到 5。

实用技巧

• 遍历数组：直接使用数组作为循环范围。

  arr = [10, 20, 30];
  for val = arr
      disp(val)
  end

  输出结果为 10、20、30。

• 反向循环：设置负步长。

  for i = 5:-1:1
      disp(i)
  end

  输出从 5 递减到 1。

注意事项

• 循环变量在循环结束后仍保留最后一次的值。
• 避免在循环体内修改循环变量（如 i = i + 1），这可能导致逻辑混乱。

二、while 循环：条件满足时持续执行

while 循环适用于不知道具体执行次数，但有一个明确的终止条件的情况，比如等待用户输入有效数据或逼近某个数值解。

基本语法

while 条件表达式
    % 循环体
end

• 条件表达式：只要结果为“真”（非零），就继续循环。
• 必须确保循环体内有改变条件的语句，否则会陷入死循环。

操作步骤

1. 创建新脚本（例如 example_while.m）。
2. 输入以下代码：

   x = 1;
   while x < 5
       disp(x);
       x = x + 1;
   end

3. 运行脚本，输出 1 到 4。

实用场景：数值逼近

假设你想找到最小的整数 n，使得 $2^n > 1000$：

n = 0;
while 2^n <= 1000
    n = n + 1;
end
disp(n);

运行后输出 10，因为 $2^{10} = 1024 > 1000$。

安全措施

• 始终初始化循环变量（如上面的 x = 1 或 n = 0）。
• 在循环体内更新影响条件的变量（如 x = x + 1）。
• 若不确定是否会无限循环，可在循环中加入计数器强制退出：

  count = 0;
  while condition && count < 1000
      % 执行操作
      count = count + 1;
  end

三、for 与 while 的选择指南

根据任务特点选择合适的循环结构，可大幅提升代码清晰度和效率。

示例对比：计算阶乘

• 用 for 循环（已知从 1 乘到 n）：

  n = 5;
  result = 1;
  for i = 1:n
      result = result * i;
  end
  disp(result); % 输出 120

• 用 while 循环（通过递减计数器实现）：

  n = 5;
  result = 1;
  temp = n;
  while temp > 0
      result = result * temp;
      temp = temp - 1;
  end
  disp(result); % 同样输出 120

虽然两者都能完成任务，但 for 版本更直观，代码行数更少。

四、常见错误与调试方法

错误1：死循环

现象：程序卡住无响应。
原因：while 循环条件永远为真。
解决：

• 检查循环体内是否有更新条件变量的语句。
• 临时添加显示语句观察变量变化：

  while x < 10
      disp(['当前 x = ', num2str(x)]);
      % 忘记写 x = x + 1; 就会导致死循环
  end

错误2：索引越界

现象：报错 “Index exceeds matrix dimensions”。
原因：for 循环中索引超出数组长度。
解决：

• **使用 length 或 size 动态获取数组长度：

  arr = [1, 2, 3];
  for i = 1:length(arr)
      disp(arr(i));
  end

错误3：浮点数比较陷阱

现象：while 条件看似应终止却未终止。
原因：浮点运算存在精度误差，如 0.1 + 0.2 ~= 0.3。
解决：

• 避免直接用 == 比较浮点数，改用误差容忍：

  tolerance = 1e-6;
  while abs(x - target) > tolerance
      % 迭代更新 x
  end

五、性能优化建议

MATLAB 对向量化操作高度优化，优先考虑不用循环。但在无法向量化时，注意以下几点：

1. 预分配数组：在循环前用 zeros 或 ones 创建结果数组，避免动态扩容。

   % 不推荐：每次循环都扩展数组
   result = [];
   for i = 1:1000
       result = [result, i^2];
   end
   
   % 推荐：预先分配
   result = zeros(1, 1000);
   for i = 1:1000
       result(i) = i^2;
   end

2. 避免在循环内调用耗时函数：如绘图、文件读写等，尽量移到循环外。

3. 用内置函数替代简单循环：例如用 sum(A) 代替手动累加。

点击 MATLAB 主页的“新建脚本”按钮，复制上述任意代码段，运行验证效果。熟悉 for 和 while 的差异后，你就能根据问题本质选择最合适的循环结构。