使 pde 解决方案适合数据

Question

我已经使用 fipy 解决了 pde（非均匀扩散率的第二菲克定律）。

我现在想通过将 pde 的解拟合到测量值来估计其中一层的扩散系数。

您对如何进行有什么建议吗？

非常感谢！

Answer 1

有趣的问题。正如@wd15 建议的那样，scipy.optimize 有有用的例程。

我整理了一个notebook，它将scipy.optimize.curve_fit 中的示例代码与examples.diffusion.mesh1D 中具有非均匀扩散率的稳态扩散示例相结合。

您的GitHub issue 可能更适合对此问题进行任何扩展讨论。