[论文解读] Implicit higher derivatives, and a formula of Comtet and Fiolet

Tom Wilde|ArXiv.org|May 17, 2008

Mathematics and Applications参考文献 1被引用 5

一句话总结

本文更正了 Comtet 和 Fiolet（1974 年）在计算隐函数 $ y = f(x) $ 满足 $ F(x,y) = 0 $ 时高阶导数 $ d^n y / dx^n $ 的有缺陷的组合公式。通过拉格朗日反演和数学归纳法，作者推导出一个更正后的公式，该公式以满足 $ (0,1) \notin p $ 的二维整数分拆 $ p \vdash (n, |p|-1) $ 为求和指标，包含系数 $ \alpha_p $ 和 $ F $ 的偏导数，分母为 $ F_y^{|p|} $。主要贡献在于给出了隐函数高阶导数的精确、无误表达式，以及对公式中项数的更正计数。

ABSTRACT

Let F(x,y) be a function of two variables, and suppose y = f(x) satisfies F(x,y)=0 in some range. Then dy/dx = -Fx/Fy, where Fx and Fy denote the partial derivatives of F with respect to x and y. It is natural to seek a general expression for the higher derivatives d^ny/dx^n, in terms of partial derivatives of F, and such an expression was given in 1974 by L. Comtet and M. Fiolet. Their formula, however, contains some errors. In this note, we give a corrected expression. We give a derivation using Lagrange inversion and also an elementary proof by induction. We further correct a minor error in Comtet and Fiolet's expression for the number of terms in their formula.

研究动机与目标

• 更正 Comtet 和 Fiolet 于 1974 年发表的关于隐函数 $ y = f(x) $ 满足 $ F(x,y) = 0 $ 时 $ d^n y / dx^n $ 的错误组合公式。
• 通过两种独立方法——拉格朗日反演与数学归纳法——严格推导更正后的公式。
• 更正 Comtet 和 Fiolet 对其公式中项数表达的错误，该错误涉及对分拆 $ p \vdash (n, |p|-1) $ 的计数，且排除 $ (0,1) $。
• 建立更正后项数 $ a(n) $ 的生成函数，即满足条件的此类分拆的数量，并计算其在 $ n = 24 $ 以内的值。

提出的方法

• 通过一种拉格朗日反演的形式推导更正后的公式，将 $ d^n y / dx^n $ 表示为满足 $ (0,1) \notin p $ 的二维整数分拆 $ p \vdash (n, |p|-1) $ 的求和。
• 应用数学归纳法验证更正后公式对 $ d^n y / dx^n $ 的正确性，确保其在各阶导数下的一致性。
• 引入二维整数分拆 $ p \vdash (m,n) $ 的记号，采用字典序排列，并定义多重性计数 $ e_{p,k,l} $。
• 定义系数 $ \alpha_p = \frac{n! m!}{\prod_i p_{i1}! p_{i2}! \prod_{k,l} e_{p,k,l}!} $，其形式类似于 Faà di Bruno 公式，用于加权每个分拆项。
• 使用生成函数 $ \prod_{(i,j) \in E} \frac{1}{1 - t^i u^{i+j-1}} $，其中 $ E = \mathbb{N} \times \mathbb{N} \setminus \{(0,0),(0,1)\} $，将项数 $ a(n) $ 计算为 $ t^n u^{n-1} $ 的系数。
• 通过递推公式 $ p_n(t) = \frac{1}{n} \sum_{s=1}^n s q_s(t) p_{n-s}(t) $ 实现高效计算，其中 $ q_m(t) $ 由对数生成函数导出，从而高效计算 $ a(n) $。

实验结果

研究问题

• RQ1 当 $ y = f(x) $ 由 $ F(x,y) = 0 $ 隐式定义时，$ n $ 阶导数 $ d^n y / dx^n $ 的正确组合公式是什么？
• RQ2 Comtet 和 Fiolet 的原始公式为何存在错误？其推导过程中的错误根源是什么？
• RQ3 公式中不同项数 $ a(n) $ 的正确表达式是什么？如何计算该值？
• RQ4 如何通过拉格朗日反演与数学归纳法推导更正后的公式？二维整数分拆在该推导过程中起到什么作用？

主要发现

• 更正后的 $ d^n y / dx^n $ 公式为 $ \sum_{p: p \vdash (n, |p|-1), (0,1) \notin p} (-1)^{|p|} \alpha_p \frac{F_{x,y,p}}{F_y^{|p|}} $，其中 $ \alpha_p $ 通过二维整数分拆上的多项式系数定义。
• 当 $ n = 2 $ 时，公式给出 $ \frac{d^2 y}{dx^2} = -\frac{F_{xx}}{F_y} + 2\frac{F_x F_{xy}}{F_y^2} - \frac{F_x^2 F_{yy}}{F_y^3} $，与从一阶导数直接微分的结果一致。
• Comtet 和 Fiolet 声称的项数 $ a(n) $ 表达式有误；正确计数应为生成函数 $ \prod_{(i,j) \in E} \frac{1}{1 - t^i u^{i+j-1}} $ 中 $ t^n u^{n-1} $ 的系数，而非 $ \prod_{(i,j) \in E} (1 - t^i u^j)^{-1} $ 中的系数。
• 序列 $ a(n) $ 从 $ n = 1 $ 到 $ 24 $ 的值为 $ 1, 3, 9, 24, 61, 145, 333, 732, 1565, 3247, 6583, 13047, 25379, 48477, 91159, 168883, 308736, 557335, 994638, 1755909, 3068960, 5313318, 9118049, 15516710 $，与 Comtet 和 Fiolet 的值在 $ n = 23 $ 之前一致，但在 $ n = 24 $ 时出现差异。
• 通过递推公式 $ p_n(t) = \frac{1}{n} \sum_{s=1}^n s q_s(t) p_{n-s}(t) $，其中 $ q_m(t) $ 由 $ \log F(u,t) $ 导出，可高效计算 $ a(n) $ 至 $ n = 50 $。