[论文解读] Properties of the map associated with recovering of the Sturm-Liouville operator by its spectral function. Uniform stability in the scale of Sobolev spaces

А. М. Савчук, Андрей Андреевич Шкаликов|arXiv (Cornell University)|Oct 26, 2010

Spectral Theory in Mathematical Physics参考文献 3被引用 5

一句话总结

本论文为在Sobolev空间中从谱数据恢复Sturm-Liouville势的反问题建立了统一的稳定性估计。它引入了一个非线性映射 $ F: W_2^\theta \to l_D^\theta $，将势的原函数 $ \sigma = \int q $ 映射到正则化谱数据，证明该映射在局部范围内可逆，并给出了显式的Fréchet导数；同时在 $ l_D^\theta $ 的凸子集上建立了统一的双边估计 $ \|\sigma - \sigma_1\|_\theta \asymp \|\mathbf{s} - \mathbf{s}_1\|_\theta $，确保了在全部Sobolev正则度 $ \theta \geq 0 $ 下的鲁棒重构。

ABSTRACT

Denote by $L_D$ the Sturm-Liouville operator $Ly=-y" +q(x)y$ on the finite interval $[0,π]$ with Dirichlet boundary conditions $y(0)=y(π)=0$. Let $\{λ_k\}_1^\infty$ and $\{α_k\}_1^\infty$ be the sequences of the eigenvalues and norming constants of this operator. For all $θ\geqslant 0$ we study the map $F: W_2^θ o l_D^θ$ defined by $F(σ) =\{s_k\}_1^\infty$. Here $σ= \int q $ is the primitive of $q$, $\bold s = \{s_k\}_1^\infty$ be regularized spectral data defined by $s_{2k} =\sqrt{λ_k}-k,\ s_{2k-1}=α_k-π/2$ and $l_D^θ$ are special Hilbert spaces which are constructed in the paper as finite dimensional extensions of the usual weighted $l_2$ spaces. We give a complete characterization of the image of this nonlinear operator, show that it is locally invertible analytic map, find explicit form of its Frechet derivative. The main result of the paper are the uniform estimates of the form $\|σ-σ_1\|_θ\asymp \|\bold s -\bold s_1\|_θ$, provided that the spectral data $\bold s$ and $\bold s_1$ run through special convex sets in the spaces $l_D^θ$.

研究动机与目标

• 解决从谱函数重构Sturm-Liouville势 $ q \in W_2^\alpha[0,\pi] $ 的反问题，扩展至 $ \alpha \geq -1 $。
• 刻画非线性映射 $ F: W_2^\theta \to l_D^\theta $ 的像，该映射将势的原函数 $ \sigma = \int q $ 映射到正则化谱数据 $ \mathbf{s} = \{s_k\} $，其中 $ s_{2k} = \sqrt{\lambda_k} - k $，$ s_{2k-1} = \alpha_k - \pi/2 $。
• 在 $ l_D^\theta $ 的凸子集上建立统一的双边稳定性估计 $ \|\sigma - \sigma_1\|_\theta \asymp \|\mathbf{s} - \mathbf{s}_1\|_\theta $，适用于 $ \theta > 0 $，确保重构的鲁棒性。
• 提供 $ F $ 及其逆的显式公式，实现局部解析反演。
• 通过非线性插值理论将框架扩展至完整的Sobolev空间尺度 $ W_2^\theta $，$ \theta \geq 0 $，并构造Hilbert空间 $ l_D^\theta $ 作为加权 $ \ell_2 $ 空间的有限维扩展。

提出的方法

• 通过将势 $ q \in W_2^\theta $ 的原函数 $ \sigma = \int q $ 关联到正则化谱数据 $ \mathbf{s} = \{s_k\} $ 来定义映射 $ F: W_2^\theta \to l_D^\theta $，其中 $ s_{2k} = \sqrt{\lambda_k} - k $，$ s_{2k-1} = \alpha_k - \pi/2 $，$ \lambda_k $ 和 $ \alpha_k $ 分别为Sturm-Liouville算子 $ L_D = -y'' + q(x)y $ 在区间 $[0,\pi] $ 上带Dirichlet边界条件的特征值与归一化常数。
• 将Hilbert空间 $ l_D^\theta $ 构造为加权 $ \ell_2 $ 空间的有限维扩展，以容纳正则化谱数据，并确保映射 $ F $ 的完备性与有界性。
• 证明 $ F $ 是从Sobolev空间 $ W_2^\theta $ 到空间 $ l_D^\theta $ 的局部可逆解析映射，并显式计算其Fréchet导数。
• 通过使用凸集 $ \Omega^\theta(r,h) \subset l_D^\theta $ 和球 $ \mathcal{B}_\Gamma^\theta(R) \subset W_2^\theta $ 建立统一的稳定性估计，表明势的差的 $ \theta $-范数与谱数据差的 $ \theta $-范数等价。
• 应用非线性插值理论，统一分析整个Sobolev空间尺度 $ W_2^\theta $，$ \theta \geq 0 $ 上的稳定性，确保估计具有鲁棒性且尺度不变。
• 通过仿射平移（添加 $ cx $）实现 $ \mathcal{B}_\Gamma^\theta(R) $ 与 $ \mathcal{B}_\mathbb{R}^\theta(R) $ 之间的等价性，简化分析并将其结果扩展至标准Sobolev球。

实验结果

研究问题

• RQ1 映射 $ F: W_2^\theta \to l_D^\theta $ 将势的原函数 $ \sigma = \int q $ 映射到正则化谱数据的像的精确刻画是什么？
• RQ2 映射 $ F $ 的Fréchet导数如何显式计算？其结果揭示了反问题的局部可逆性有何特征？
• RQ3 在 $ W_2^\theta $-范数下势的差与 $ l_D^\theta $-范数下其正则化谱数据差之间可建立何种统一的双边估计？
• RQ4 估计中的稳定性常数如何依赖于数据集的大小？当半径 $ R \to \infty $ 或 $ h \to 0 $ 时，情况如何？
• RQ5 能否将结果从特殊集合 $ \mathcal{B}_\Gamma^\theta(R) $ 扩展至标准Sobolev球 $ \mathcal{B}_\mathbb{R}^\theta(R) $？谱数据定义需要进行何种修改？

主要发现

• 映射 $ F: W_2^\theta \to l_D^\theta $，定义为 $ F(\sigma) = \{s_k\} $，其中 $ s_{2k} = \sqrt{\lambda_k} - k $，$ s_{2k-1} = \alpha_k - \pi/2 $，是Sobolev空间 $ W_2^\theta $ 与所构造的Hilbert空间 $ l_D^\theta $ 之间的局部可逆解析映射，确保了反问题在局部范围内的适定性。
• 映射 $ F $ 的Fréchet导数被显式计算，并证明为有界线性算子，从而可在Sobolev空间尺度上应用隐函数定理技术。
• 在 $ \theta > 0 $ 条件下，若谱数据 $ \mathbf{s}, \mathbf{s}_1 $ 属于凸子集 $ \Omega^\theta(r,h) \subset l_D^\theta $，则建立了统一的双边估计 $ \|\sigma - \sigma_1\|_\theta \asymp \|\mathbf{s} - \mathbf{s}_1\|_\theta $，其中常数 $ C_1, C_2 $ 仅依赖于 $ r $ 和 $ h $。
• 逆映射 $ F^{-1} $ 在 $ \Omega^\theta(r,h) $ 上一致有界，满足 $ \|(F^{-1})'\| \leq C $，其中 $ C $ 仅依赖于 $ r $ 和 $ h $，确保了重构过程的稳定性。
• 稳定性估计具有鲁棒性：当 $ R \to \infty $ 时，常数 $ r $、$ h^{-1} $、$ C_2 $ 和 $ C_1^{-1} $ 虽然增大，但估计在指定的凸集上仍保持统一。
• 通过使用仿射平移 $ \sigma \mapsto \sigma + cx $，将结果从特殊集合 $ \mathcal{B}_\Gamma^\theta(R) $ 扩展至标准Sobolev球 $ \mathcal{B}_\mathbb{R}^\theta(R) $，并重新定义谱数据以包含常数 $ c $，使得 $ c \geq -\lambda_1(\sigma) - 1 $，从而保证了估计中范数的等价性。