[论文解读] D-Modules and Holonomic Functions

Anna-Laura Sattelberger, Bernd Sturmfels|arXiv (Cornell University)|Oct 3, 2019

Polynomial and algebraic computation参考文献 21被引用 9

一句话总结

本文為應用代數幾何中的 D-模與正則函數提供了導論，專注於使用 Weyl 代數中的 Gröbner 基來以消去 D-理想表示正則函數。透過利用正則秩、D-模理論以及有理函數與參數化積分的演算法，該方法實現了統計學中精確的體積計算與概似比推斷。

ABSTRACT

In algebraic geometry, one studies the solutions to polynomial equations, or, equivalently, to linear partial differential equations with constant coefficients. These lecture notes address the more general case when the coefficients are polynomials. The letter D stands for the Weyl algebra, and a D-module is a left module over D. We focus on left ideals, or D-ideals. We represent holonomic functions in several variables by the linear differential equations they satisfy. This encoding by a D-ideal is useful for many problems, e.g., in geometry, physics and statistics. We explain how to work with holonomic functions. Applications include volume computations and likelihood inference.

研究动机与目标

• 為應用代數幾何領域的研究者提供多變數 D-模與正則函數的基礎介紹。
• 展示如何將正則函數以 Weyl 代數中的 D-理想編碼，以利幾何、物理與統計學中的計算應用。
• 透過參數化正則函數，實現對半代數集上積分的精確數值評估。
• 透過正則梯度法與概似函數，建立 D-模理論與統計推斷之間的聯繫。
• 提供一個使用電腦代數系統（如 Singular 和 Mathematica）計算特殊函數與概似函數的消去理想的框架。

提出的方法

• 在 Weyl 代數 $D = \mathbb{C}[x]\langle\partial\rangle$ 中使用 Gröbner 基來計算與操作消去正則函數的 D-理想。
• 應用代數分析基本定理，透過有限維解空間來特徵化正則 D-理想。
• 將積分參數化為新變數的正則函數，隨後推導出消去體積函數的 D-理想。
• 應用正則梯度法與正則梯度下降法，以優化與評估指數族中的概似函數。
• 利用專用電腦代數工具，如 Mathematica 中的 HolonomicFunctions 套件與 Singular 中的 dmod.lib 庫進行理想計算。
• 使用有理 Weyl 代數 $R = \mathbb{C}(x)\langle\partial\rangle$ 來處理有理函數，並計算正則函數之和與積的消去理想。

实验结果

研究问题

• RQ1 如何在 Weyl 代數中有效使用 D-理想來編碼與操作多變數正則函數？
• RQ2 Gröbner 基在 Weyl 代數中於計算特殊函數與積分的消去理想時扮演何種角色？
• RQ3 如何利用正則 D-理想精確評估半代數集上的參數化積分？
• RQ4 D-模理論在離散指數模型中的概似推斷方面，能以何種方式提供增益？
• RQ5 何種結構性質（如正則秩或初始理想的素性）有助於判斷 D-理想是否為最大或完整？

主要发现

• 當參數化後，半代數集的體積函數為正則函數，且其存在可透過 D-模技術實現高度精確的數值評估之 D-理想。
• 正則梯度法與正則梯度下降法提供了有效演算法，可用 D-理想表示來優化與評估統計學中的概似函數。
• 在 Mathematica 的 HolonomicFunctions 套件中，可使用 DFinitePlus 與 DFiniteTimes 指令計算正則函數之和與積的消去理想。
• 對於 $n=3$、$k=4$ 的隨機右刪失模型，概似函數 $f^s$ 存在可透過 Singular 使用 dmod.lib 與 dmodapp.lib 庫計算的 D-理想。
• 概似函數的理想 $\text{Ann}_{D[s]}(f^s)$ 透過適應 [SST00, Algorithm 5.3.15] 至有理函數而計算得出，確認了該模型的正則性。
• 若 $\text{in}_{(0,e)}(I)$ 為素理想，則正則 D-理想 $I$ 為最大理想，此提供了一項驗證消去理想完整性之準則。