使PDE解决方案适合数据

问题描述

我已经使用fipy解决了pde(扩散率不均匀的第二菲克定律)。

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

您对如何进行有任何建议吗?

非常感谢!

解决方法

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

我将notebook的示例代码与我们的稳态扩散示例(在scipy.optimize.curve_fit中具有非均匀扩散性)组合在一起。

您的examples.diffusion.mesh1D对于此问题的扩展讨论可能会更好。