我在对流扩散 PDE 中确实有一个扩散系数,它是不连续的,看起来像(x 轴上的浓度):

出于数字原因,我使用综合形式:

我预先计算了函数,在计算过程中,我根据每个网格点的实际浓度值对扩散系数进行插值。请注意,函数在 x = 0.1 和 x = 0.5 时是不连续的
如果我可以为数据拟合一个分析函数,我可以避免在每个时间步进行插值,并可以加快计算速度。这可能吗?这有意义吗?哪个功能有意义?就像是:
myfun = somefunction*((u>=0.1) && (u<=0.5)
我在对流扩散 PDE 中确实有一个扩散系数,它是不连续的,看起来像(x 轴上的浓度):

出于数字原因,我使用综合形式:

我预先计算了函数,在计算过程中,我根据每个网格点的实际浓度值对扩散系数进行插值。请注意,函数在 x = 0.1 和 x = 0.5 时是不连续的
如果我可以为数据拟合一个分析函数,我可以避免在每个时间步进行插值,并可以加快计算速度。这可能吗?这有意义吗?哪个功能有意义?就像是:
myfun = somefunction*((u>=0.1) && (u<=0.5)
如果您的数据在一个统一的网格上可用,那么您不太可能加快计算速度:您目前所要做的只是一次除法找出你所在的区间,以及数据点之间的线性插值和. 几乎任何不是低阶多项式的函数都将比这些操作花费更多。
但是,如果您坚持要为您的数据拟合一个函数,那么它的来源会很有帮助。曲线看起来像高斯曲线,但它是否真的是当然取决于它的模型。