最近要用matlab画二维网格，网上查阅和实践后，做一下总结。
一、利用meshgrid绘制二维网格的正确做法
xtick = linspace(10,20,51);
ytick = linspace(0,2*pi,21);
[X,Y] = meshgrid(xtick,ytick);
plot(X, Y, 'k', X', Y', 'k');
axis([10 20 0 2*pi]);
绘图结果如下：
到这儿就结束了。接下来的内容是我走过的坑和原因分析，不感兴趣的可以直接跳过。
二、错误示例
一开始经过网上搜索，发现有二维网格绘制的相关例子，一个简单例子如下：
xtick = linspace(0,10,21);
ytick = linspace(0,10,21);
[X,Y] = meshgrid(xtick,ytick);
plot(X,Y,'k',Y,X,'k');
绘图如下：
一看就是我想要的图，直接根据这个例子，开始动手：
xtick = linspace(10,20,51);
ytick = linspace(0,2*pi,21);
[X,Y] = meshgrid(xtick,ytick);
plot(X, Y, 'k', Y, X, 'k');
画出来是这样的，傻眼了。。。网格没有画到一起，和我想的不一样。。。
经过仔细研究meshgrid这个函数的输出，以及画出来的图，又重新修改了程序，如第一部分的正确示例所示（只有plot命令做了修改）。
三、原因
首先研究一下Meshgrid的输入输出
直接截取官方说明：
[X,Y] = meshgrid(x,y) 基于向量 x 和 y 中包含的坐标返回二维网格坐标。X 是一个矩阵，每一行是 x 的一个副本；Y 也是一个矩阵，每一列是 y 的一个副本。坐标 X 和 Y 表示的网格有 length(y) 个行和 length(x) 个列。
举个官方的例子：
x=1:3;
y=1:5;
[X,Y]=meshgrid(x,y);
输出X,Y看看长什么样
先plot(X,Y)，查看是如何画图的，如下图：有三条竖直线，分别是x=1，x=2和x=3。其画图次序为：X第1列和Y的第1列（红线），X第2列和Y第2列（绿线），以及X第3列和Y的第3列（蓝线）。
接着plot(Y,X), 继续查看绘图情况，如下图：有三条水平线，分别是y=1，y=2和y=3。其画图次序为：Y第1列和X的第1列（红线），Y第2列和X第2列（绿线），以及Y第3列和X的第3列（蓝线）。
可以发现，利用plot命令画矩阵时，它会按照对应列进行绘图。且，第一个矩阵会作为横坐标值、第二个矩阵作为纵坐标值进行绘图。
接着我们再研究一下，如何在plot(X,Y)的基础上继续把我们想要的水平线画上，即y=1，y=2,…,y=5。
先看一下下图中的y=1的水平线，其坐标点分别为(1,1),(2,1)和(3,1)，这三个坐标值分别对应X矩阵的第1行和Y矩阵的第1行。这样我们就可以把X和Y矩阵转置再plot(X',Y')就可以得到5条水平线。
最后，plot(X,Y,'k',X',Y','k');可以得到下图。
总结：出错的原因在于没有真正理解用plot命令画两个矩阵时程序是如何执行的。
注：如果横坐标和纵坐标的坐标范围和格点步长是一样的，那么X’和Y是相同的，Y’和X是相同的，所以用plot(X,Y,'k',Y,X,'k')不会出现问题。
看来还是不能偷懒-_-!
(原创-转载请注明出处)
参考：
https://ww2.mathworks.cn/help/matlab/ref/meshgrid.html
https://blog.csdn.net/zhanghanqmx/article/details/121891071