求助：用相位展开消除复方程分支切割引发的等值线图不连续性

太懂这种卡壳的感觉了！水平方向的分支切割用unwrap一键搞定，结果碰到垂直方向的就直接哑火，确实头疼。不过别慌，这里有几个实用的办法帮你解决：

• 转置矩阵适配unwrap的默认逻辑
  unwrap默认是沿着数组的列（水平方向）处理相位跳变的，那咱们就把垂直方向的问题转成它擅长的水平方向：先把你的相位矩阵转置，用unwrap处理完再转回来就行。代码示例：

  % 假设phase_mat是你的原始相位矩阵，y轴方向存在相位跳变
  phase_unwrapped = unwrap(phase_mat')';
  

• 直接指定unwrap的展开维度
  如果你的MATLAB版本是R2016b及以后，其实unwrap支持指定展开维度！不用转来转去，直接告诉它要沿着垂直方向（行维度）处理：

  % 维度1对应行方向（也就是你的y轴方向），维度2对应列方向
  phase_unwrapped = unwrap(phase_mat, [], 1);
  

  这个方法最直接，省了转置的步骤，效率也更高。

• 手动遍历修正（兼容旧版本）
  要是你的MATLAB版本比较老，不支持指定维度参数，那咱们就手动来处理垂直方向的每一列：遍历每一列的相邻行，检查相位差，超过π就加减2π来修正跳变。代码示例：

  [rows, cols] = size(phase_mat);
  phase_unwrapped = phase_mat;
  % 遍历每一列
  for c = 1:cols
      % 从第二行开始逐行检查
      for r = 2:rows
          phase_diff = phase_unwrapped(r,c) - phase_unwrapped(r-1,c);
          % 相位差超过π，就给后面的行统一调整
          if phase_diff > pi
              phase_unwrapped(r:end,c) = phase_unwrapped(r:end,c) - 2*pi;
          elseif phase_diff < -pi
              phase_unwrapped(r:end,c) = phase_unwrapped(r:end,c) + 2*pi;
          end
      end
  end
  

  这个方法兼容性拉满，就是代码稍微繁琐一点，不过对付垂直跳变完全够用。

处理完之后记得画个等值线图对比一下，看看分支切割是不是已经消失啦。要是你的分支切割不是严格垂直，可能需要再结合方向判断微调，但大部分垂直跳变的情况，上面前两种方法就能轻松搞定。

内容的提问来源于stack exchange，提问作者Tina Joseph