问题描述
我已经使用fipy解决了pde(扩散率不均匀的第二菲克定律)。
我现在想通过将pde的解拟合到测量值来估计其中一层的扩散系数。
您对如何进行有任何建议吗?
非常感谢!
解决方法
有趣的问题。正如@ wd15所建议的那样,scipy.optimize
具有有用的例程。
我将notebook的示例代码与我们的稳态扩散示例(在scipy.optimize.curve_fit
中具有非均匀扩散性)组合在一起。
您的examples.diffusion.mesh1D
对于此问题的扩展讨论可能会更好。