矢量场：quiver函数如何绘制流体力学中的速度矢量场分布

在流体力学仿真或实验数据处理中，可视化速度分布是个关键环节。但你是否遇到过这样的困扰：用MATLAB的quiver函数画出来的矢量图，箭头方向看着别扭、密度失真，或者箭头的长度根本反映不了真实的速度大小？

别急，这通常不是什么高深的理论问题，根源往往在于两个基础但至关重要的技术细节：坐标网格与速度分量矩阵的维度是否严格匹配，以及缩放参数（scale）设置是否得当。只要把这两点理顺，一张清晰、准确、专业的矢量场图就能轻松呈现。

一、确保坐标网格与速度分量维度严格一致

这是最常踩的“坑”。quiver函数要求起点坐标（X, Y）与速度分量（U, V）必须是尺寸完全相同的矩阵。如果维度对不上，MATLAB要么报错，要么进行隐式扩展，结果就是矢量位置全乱套，图形变得毫无意义。

最稳妥的办法，就是使用meshgrid函数来生成规则网格：

1. 先定义x和y方向的一维坐标向量。步长根据你的数据空间分辨率来定，比如：x = -2:0.4:2; y = -1.5:0.3:1.5;

2. 调用[X, Y] = meshgrid(x, y)生成二维网格矩阵。这里有个细节要注意：在MATLAB中，meshgrid输出的X矩阵，其行数等于y向量的长度，列数等于x向量的长度，别把行列关系搞反了。

3. 计算网格点上对应的速度分量U和V。确保你写的公式最终生成的U矩阵与X同尺寸，V与Y同尺寸。例如：U = cos(Y).*exp(-X.^2); V = sin(X).*exp(-Y.^2);

4. 执行quiver(X, Y, U, V)。如果这时MATLAB弹出“Matrix dimensions must agree”的警告，别犹豫，立刻用size(X)和size(U)检查一下，它们的维度必须百分百相等。

二、使用scale参数校准箭头长度以真实反映速度模值

默认情况下，quiver为了图形美观，会对所有箭头进行统一的自动缩放。这虽然让图看起来清爽，但却牺牲了速度绝对大小的信息，你无法从箭头长度直接判断哪里快、哪里慢。

想让箭头长度真实反映速度大小，关键在于scale参数：

1. 先计算速度模值矩阵：mag = sqrt(U.^2 + V.^2)。这个矩阵后面用来上色或者筛选都很有用。

2. 调用quiver(X, Y, U, V, 1)。把scale参数设为1，意味着每个箭头的显示长度将严格等于其速度模值sqrt(u²+v²)，此时图形具有几何真实性。

3. 如果设成1后，箭头太长导致彼此重叠、难以辨认，可以尝试将scale设为一个小于1的正数（比如0.8或0.5）。这相当于对所有的箭头长度进行线性压缩，公式是：显示长度 = scale × sqrt(u²+v²)。

4. 对于速度特别低的区域（比如mag < 0.1*max(mag(:))），可以考虑不画箭头，或者用特殊的标记（比如点）来表示，避免图上出现一堆几乎看不见的短箭头，影响整体视觉效果。

三、叠加散点图实现速度大小彩色编码

单靠箭头长度，人眼很难分辨细微的速度差异。这时候，给箭头“上色”就成了增强定量辨识能力的利器。结合scatter函数，可以基于速度模值为每个网格点填充颜色。

1. 先用quiver(X, Y, U, V, 1, 'k', 'LineWidth', 0.8)画出黑色的箭头骨架。

2. 紧接着调用hold on命令，保持当前坐标轴，以便叠加图形。

3. 执行scatter(X(:), Y(:), 30, mag(:), 'filled')。这里，30是散点的大小，mag(:)是将速度模值矩阵展平成一列，作为颜色映射的依据。'filled'让散点是实心的。

4. 最后，添加colorbar，并用colormap设置一个好看的色系，比如jet或turbo。**关键一步：用caxis([min(mag(:)) max(mag(:))])固定色标的范围**。这样，当你需要对比多张图时，颜色代表的数值意义是一致的，不会因单张图的数据范围不同而产生误导。

四、控制矢量密度防止视觉过载

处理高分辨率PIV（粒子图像测速）数据时，动辄几千上万个网格点。如果全部画成箭头，结果就是黑压压一片，什么也看不清。对数据进行适当的下采样，保留有代表性的矢量，是绘图前的必要预处理。

1. 根据你的数据量和图形大小，确定x和y方向的下采样步长，比如step_x = 3; step_y = 2;。

2. 对坐标网格和速度分量矩阵进行切片抽取： X_sub = X(1:step_y:end, 1:step_x:end); 同理，提取出对应的Y_sub, U_sub, V_sub。

3. 务必验证一下：size(X_sub)必须等于size(U_sub)，并且它们的行列数应该是原始矩阵按步长下采样后的结果。

4. 用下采样后的数据绘图：quiver(X_sub, Y_sub, U_sub, V_sub, 0.7)。**对于下采样后的图，建议将scale初始值设为0.5到0.8之间**，这样能在箭头长度辨识度和图形密度之间取得一个较好的平衡。

五、在三维流场中使用quiver3绘制截面速度矢量

当你的流场是三维的，包含了z方向的速度分量时，就需要请出quiver3函数了。它的核心要求同样严格：x, y, z, u, v, w这六组数据，必须全部是同型的三维矩阵。通常，我们借助meshgrid来生成基础坐标。

1. 定义三维离散域：x = linspace(-1,1,11); y = x; z = linspace(-0.5,0.5,7);

2. 执行[X,Y,Z] = meshgrid(x,y,z);，这会生成三个尺寸为11×11×7的三维矩阵。

3. 计算三维速度分量U, V, W。确保每个公式算出来的结果都是11×11×7的矩阵。例如：U = Z.*cos(X); V = Z.*sin(Y); W = X.*Y;

4. 调用quiver3(X,Y,Z,U,V,W,'AutoScale','on','MaxHeadSize',0.02)进行绘制。**这里将'AutoScale'设置为'on'非常有用，它能自动适配不同截面或不同区域的速度尺度，避免因局部速度过大或过小导致图形比例失调。**'MaxHeadSize'参数则可以控制箭头头部的大小，让图形更美观。