我正在尝试根据其拉普拉斯算子重建图像，这会产生泊松方程，并且我正在使用 Neumann 边界条件（边界处的导数 = 0）。

我所拥有的是拉普拉斯算子（f，左侧）和图像边界（边界）。

如何在 Matlab 或 Octave 中解决这个问题？我找到了poisolv()andassempde()函数，但我不太明白它们是如何工作的。