[論文レビュー] Singular-value decomposition of solution operators to model evolution equations

Alexandru Aleman, Joe Viola|arXiv (Cornell University)|Sep 3, 2014

Stability and Controllability of Differential Equations参考文献 3被引用数 4

ひとこと要約

本稿では、$L^2(\mathbf{R}^n)$ 上の二次型作用素によって生成される発展方程式の解作用素に対する特異値分解を、ユニタリ変換を用いて導出し、解を調和振動子半群の和に還元する。主な結果は、楕円性の仮定なしに、$e^{-tM}$ の特異値を用いて減衰と正則化を正確に特徴づけることであり、長時間行動は $M$ の固有値スペクトルに、短時間行動は $e^{-tM}$ のノルムに支配されることを明らかにする。

ABSTRACT

We consider evolution equations generated by quadratic operators admitting a decomposition in creation-annihilation operators without usual ellipticity-type hypotheses; this class includes hypocoercive model operators. We identify the singular value decomposition of their solution operators with the evolution generated by an operator of harmonic oscillator type, and in doing so derive exact characterizations of return to equilibrium and regularization for any complex time.

研究の動機と目的

• 楕円性や強 coercivity を仮定しない二次型作用素によって生成される発展方程式の解の長時間行動を特徴づけること。
• 非楕円的で、hypocoercive 型の作用素に対しても、平衡への回帰と正則化の理論を拡張すること。
• 複素時間 $t$ に対して有効な、$e^{-tM}$ の特異値分解を用いた解作用素ノルムの正確かつ明確な記述を提供すること。
• 解作用素がユニタリ同値により調和振動子半群の和に変換可能であることを示し、スペクトル解析を可能にする。

提案手法

• Bargmann 変換を用いて、解作用素を調和振動子型作用素と関連付けること。
• 二次型作用素 $P$ を生成・消滅作用素 $A_j^*, A_j$ に分解し、Hermite 関数によるスペクトル解析を可能にする。
• 解作用素 $\exp(-tP)$ が、固有値が $\sum_j (\log \sigma_j) \alpha_j$ である自己共役作用素 $Q$ に対して $\exp(Q)$ とユニタリ同値であることを特定すること。
• $e^{-tM}$ の特異値 $\sigma_j$ を計算し、これが解作用素の減衰率を決定すること。
• ユニタリ変換 $\mathcal{U}_1, \mathcal{U}_2$ を用いて、解作用素をエネルギー部分空間 $E_m$ を保存する1次元調和振動子半群の和に還元すること。
• $\|e^{-tM}\|$ の解析により、非楕円的摂動のもとでも、短時間における非有界性と長時間における指数的減衰を区別すること。

実験結果

リサーチクエスチョン

• RQ1 楕円性や強 coercivity を欠く場合に、二次発展方程式の解作用素をどのように特徴づけられるか？
• RQ2 標準的な散逸性条件が成り立たない状況下で、このような解の長時間減衰および正則化特性は、何によって決定されるか？
• RQ3 解作用素がユニタリ同値により調和振動子半群の和に分解可能であるか？ もしそうならば、パラメータはどのように特定されるか？
• RQ4 複素時間 $t$ に対して、$e^{-tM}$ の特異値が $\exp(-tP)$ の振る舞いをどのように支配するか？ これにより、有界性およびコンパクト性にどのような意味が生じるか？
• RQ5 作用素が正規でも楕円的でもない場合でも、$M$ の固有値が $\exp(-tP)$ の長時間行動をどの程度決定するか？

主な発見

• 解作用素 $\exp(-tP)$ は、$\sigma_j$ を $e^{-tM}$ の特異値とする $Q = \frac{1}{2}\sum_{j=1}^n (\log \sigma_j)(x_j^2 - \partial_{x_j}^2 - 1)$ に対して $\exp(Q)$ とユニタリ同値であり、完全なスペクトル分解が得られる。
• ノルム $\|\exp(-tP)\|$ は特異値 $\sigma_j$ によって正確に決定され、$\|\exp(-tP)\| = \sup_{\alpha \in \mathbb{N}^n} \exp\left(\sum_{j=1}^n (\log \sigma_j) \alpha_j \right)$ と表され、$\sigma_j < 1$ ならば $t \to \infty$ で指数的に減衰する。
• Fokker-Planck モデル $P_a$ に対しては、$\|e^{-tM_a}\| = 1 - \frac{a^2}{12}t^3 + \mathcal{O}(t^4)$ が成り立ち、初期段階の緩やかな減衰を示すが、$t$ が大きいと指数的減衰を示す。
• 非楕円的摂動 $\tilde{P}_{a,b} = P_a - b(x_1^2 - \partial_{x_1}^2)$ ($b>0$) に対しても、$M_{a,b}$ のすべての固有値 $\lambda$ に対して $\operatorname{Re}(\lambda) > 0$ ならば、$\|e^{-tM_{a,b}}\| \to 0$ と指数的に減少し、長時間にわたる平衡状態への回帰が保証される。
• 解作用素 $\exp(-t\tilde{P}_{a,b})$ は、$t>0$ が大きい場合には有界であり、コンパクト性と正則化を示すが、小刻みな $t$ では非有界である。有界性領域は $\log\|e^{-tM}\|$ の等高線図から明確に観察できる。
• $\exp(-tP)$ の長時間行動は $M$ のスペクトルに支配され、短時間行動は $\|e^{-tM}\|$ に支配される。2つの領域は、前者が減衰を決定し、後者が小刻みな $t$ における非有界性を決定するという意味で、独立している。