Inverse problems involve making inference about unknown parameters of a physical process using observational data. This paper investigates an important class of inverse problems -- the estimation of the initial condition of a spatio-temporal advection-diffusion process using spatially sparse data streams. Three spatial sampling schemes are considered, including irregular, non-uniform and shifted uniform sampling. The irregular sampling scheme is the general scenario, while computationally efficient solutions are available in the spectral domain for non-uniform and shifted uniform sampling. For each sampling scheme, the inverse problem is formulated as a regularized convex optimization problem that minimizes the distance between forward model outputs and observations. The optimization problem is solved by the Alternating Direction Method of Multipliers algorithm, which also handles the situation when a linear inequality constraint (e.g., non-negativity) is imposed on the model output. Numerical examples are presented, code is made available on GitHub, and discussions are provided to generate some useful insights of the proposed inverse modeling approaches.
翻译:反面问题涉及利用观测数据推断物理过程的未知参数。本文调查了一个重要的反向问题类别 -- -- 利用空间稀少的数据流估计时空对流扩散过程的初始条件。考虑的是三个空间抽样办法,包括非常规、非统一和移动统一抽样办法。不规则抽样办法是一般情况,而光谱域则为非统一和移动统一抽样提供计算高效的解决方案。对于每个抽样办法而言,反面问题是一个正规化的平面优化问题,最大限度地减少前方模型输出和观测之间的距离。优化问题由多动算法的调整方向方法解决,该方法也处理模型输出中强加线性不平等限制(如非集中)时的情况。提出了数字实例,在吉特胡卜提供了代码,并提供了讨论,以对拟议的反型方法产生一些有用的见解。