1

我u从 Gauss-Siedel（求解 2D 泊松方程）获得了一个解，我想用trisurf. 但是，我似乎无法理解它是如何trisurf工作的（我尝试过使用它，但没有运气）。我还查看了 Matlab 文档，但我不明白峰值是什么。目前，我有这样的东西（但它不起作用）：

ss = a:h:b; %subspace
[X,Y] = meshgrid(ss,ss);
tri = delaunay(X,Y);
Z = peaks(ss);
trisurf(tri,X,Y,Z);

我在哪里插入解决方案值（向量u）？peaks在这种情况下必须采取哪些类型的价值观？

细节：

  • a= 下限
  • b= 上限
  • h=步长
  • u= Gauss-Siedel 的解。它的大小取决于我在网格上评估的点数。解决方案应该收敛，因为它是通过运行i多次迭代获得的。
4

1 回答 1

1

首先，让我给你我的输出whos

  triangle_obj       1x1               540  TriScatteredInterp              
  u                 15x1               120  double                          
  X                 15x1               120  double                          
  Y                 15x1               120  double

试试下面的方法。我认为您可能会发现这对您有用（假设您的所有点都不是共线的）。您可能会因为 TriScatteredInterp 中的 [X, Y] 参数混淆而感到困惑。

triangle_obj = TriScatteredInterp([X, Y], u);
figure
trisurf(triangle_obj, X, Y, u);
于 2013-02-27T23:24:39.300 回答