Python Matplotlib cmap未按预期工作

我为Python3编写了以下代码，它随时间生成一个波函数，并在3D中绘制结果。请注意，

schroedinger1D（…）

函数返回两个numpy数组，每个数组的形状为（361000）

---

输出的X、T和W_real.T（1000,36）的形状正确，但据我所知，曲面的颜色分配错误。我原以为Z轴上的颜色会变化，但在这里我无法判断测量的是什么：

减少网格点数量时，可能更容易理解发生了什么

xl, xr, yb, yt, M, N, xc, wl, s = (-9, 5, 0, 4, 10, 10, -5, 4.0, 3.0)

并且还使用rstride和cstride为1

在这种情况下，情节可能会像这样

现在很容易发现问题：波函数中的振荡频率大于网格的分辨率。这意味着曲面图的一个面片可能从一个非常低的值开始，然后上升到一个非常高的值。在这种情况下，其颜色可以是任何颜色，因为它由面片的一条边的值决定。（如果面片从高值开始并下降到低值，则比从低值开始并上升到高值时更红。）

---

唯一可能的解决方案是使网格比振荡频率更密集。即，尝试可视化变化更平稳的波函数，或仅可视化部分波函数，但网格更密集

你能展示剩下的代码吗？这可能是一个真正的问题，一个bug，或者输入数据的一些问题；没有人知道这个问题。@ImportanceOfBeingErnest添加：）@user2027202827添加。请使用Python3执行。明白了，我将wavefunction更改为一个更平滑的函数，并在一个小的域上执行（只是为了感受更改），然后我开始使用网格，现在我可以正确地看到颜色变化。谢谢

xl, xr, yb, yt, M, N, xc, wl, s = (-9, 5, 0, 4, 10, 10, -5, 4.0, 3.0)