泊松方程序列

Question

我想解决泊松方程的序列

  

Δu= f（k）

表示不同的k值。

我试过

for k=1:10
f='x+y+k'
u(k)=assempde('problg',p,e,t,c,a,f);
end

它返回错误

  

未定义的函数或变量'k'。

无限感激任何帮助！

Answer 1

您可以尝试使用以下代码替换第二行：

f=['x+y+' num2str(k)];

我无法测试它，因为我没有PDE工具箱。 希望这会有所帮助。

一个。

PS：如果你必须解决很多k值，有一些方法可以在给定间隔内的所有k值中一次性解决这个pde。你会得到三个变量的函数：u（x，y，k）（谁说k是参数而不是变量）。好处是内存成本和计算时间仍然会像2D泊松一样扩展，你永远不会在3D网格上解决/存储事物（k是第三维）。 如果你对那些新兴的方法感兴趣，它叫做：正确的广义分解（PGD），你可能会发现一些关于这些方法的论文。