使用 meshgrid函数绘制网格点坐标的原理与代码实现

使用 `meshgrid` 绘制网格点坐标的原理与代码实现

在 MATLAB 中，meshgrid 是一个常用函数，用于生成二维平面网格点的坐标矩阵。本文将详细介绍如何利用 meshgrid 函数生成的矩阵绘制网格点的坐标，并给出具体的代码实现和原理解析。

实现思路

meshgrid 函数生成的矩阵 wa 和 w 分别代表二维平面网格点的 $x$ 和 $y$ 坐标。我们可以将 wa 和 w 矩阵中的对应元素组合成点的 $(x, y)$ 坐标，并通过 MATLAB 的绘图函数（如 plot）进行可视化。通过这种方式，可以直观地观察网格点在平面上的分布。

代码实现

以下是具体实现的 MATLAB 代码：

% 定义 wa_range 和 w_range
wa_range = 1:5; % 1x5 的行向量
w_range = 10:13; % 1x4 的行向量

% 使用 meshgrid 生成矩阵
[wa, w] = meshgrid(wa_range, w_range);

% 绘制网格点
figure;
plot(wa(:), w(:), 'ro', 'MarkerFaceColor', 'r');
grid on;
xlabel('x 坐标 (对应 wa)');
ylabel('y 坐标 (对应 w)');
title('Meshgrid 生成的网格点坐标');

% 打印生成的矩阵
disp('矩阵 wa 的内容:');
disp(wa);
disp('矩阵 w 的内容:');
disp(w);

运行上述代码后，将会生成一个包含网格点的二维平面图。每个点的横坐标来自矩阵 wa，纵坐标来自矩阵 w。
在这里插入图片描述

结果解析

假设输入为：

wa_range = 1:5; 
w_range = 10:13;

利用 meshgrid 生成的矩阵 wa 和 w 如下：

wa =
     1     2     3     4     5
     1     2     3     4     5
     1     2     3     4     5
     1     2     3     4     5

w =
    10    10    10    10    10
    11    11    11    11    11
    12    12    12    12    12
    13    13    13    13    13

矩阵 wa 的每一行都是 wa_range，表示网格点的 $x$ 坐标。
矩阵 w 的每一列都是 w_range 的转置，表示网格点的 $y$ 坐标。

通过 wa(:) 和 w(:) 可以将上述矩阵展平为列向量：

wa(:)：按列优先排列后的 wa 元素，结果为 [1; 1; 1; 1; 2; 2; ...; 5; 5]。
w(:)：按列优先排列后的 w 元素，结果为 [10; 11; 12; 13; 10; 11; ...; 13]。

plot(wa(:), w(:), ...) 将这两个列向量中对应位置的元素组合为点的 $(x, y)$ 坐标，从而绘制出所有网格点。

`meshgrid` 原理解析

以 w 和 wa 的生成过程为例，假设 wa_range 是 $1\times n$ 的行向量，w_range 是 $1\times k$ 的行向量。

1. 理论推导

生成矩阵 w 时，可以理解为：

将 w_range 转置为 $k\times 1$ 的列向量；
通过水平复制该列向量 $n$ 次，形成 $k\times n$ 的矩阵 w。矩阵 w 的每一列都是 w_range 的转置。

生成矩阵 wa 时，可以理解为：

将 wa_range 保持为 $1\times n$ 的行向量；
通过垂直复制该行向量 $k$ 次，形成 $k\times n$ 的矩阵 wa。矩阵 wa 的每一行都是 wa_range。

2. 手动验证代码

以下代码展示了如何手动生成与 meshgrid 一致的矩阵 w 和 wa：

% 定义 wa_range 和 w_range
wa_range = 1:5;
w_range = 10:13;

% 手动生成矩阵 w 和 wa
w_simulated = repmat(w_range', 1, length(wa_range));
wa_simulated = repmat(wa_range, length(w_range), 1);

% 使用 meshgrid 生成矩阵 w 和 wa
[wa, w] = meshgrid(wa_range, w_range);

% 比较两者是否相同
is_equal_w = isequal(w, w_simulated);
is_equal_wa = isequal(wa, wa_simulated);

disp(['手动生成的矩阵 w 是否与 meshgrid 结果相同: ', num2str(is_equal_w)]);
disp(['手动生成的矩阵 wa 是否与 meshgrid 结果相同: ', num2str(is_equal_wa)]);