[논문 리뷰] Implicit higher derivatives, and a formula of Comtet and Fiolet

Tom Wilde|ArXiv.org|2008. 05. 17.

Mathematics and Applications참고 문헌 1인용 수 5

한 줄 요약

이 논문은 Comtet와 Fiolet(1974)이 제안한 고차 도함수 $ d^n y / dx^n $를 계산하기 위한 결함이 있는 조합적 공식을 수정한다. $ y = f(x) $가 $ F(x,y) = 0 $에 의해 암시적으로 정의될 때, 저자들은 라그랑주 역전개와 수학적 귀납법을 사용하여, $ (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.

연구 동기 및 목표

• 암시적으로 정의된 함수 $ y = f(x) $에 대해 $ F(x,y) = 0 $을 만족할 때, Comtet와 Fiolet이 1974년에 발표한 $ d^n y / dx^n $에 대한 잘못된 조합적 공식을 수정하는 것.
• 라그랑주 역전개와 수학적 귀납법이라는 두 가지 다른 방법을 사용하여 수정된 공식을 엄밀하게 유도하는 것.
• Comtet와 Fiolet이 잘못된 표현을 사용한 항의 수를 수정하는 것 — 이는 $ (0,1) $을 제외한 분할 $ p \vdash (n, |p|-1) $을 세는 방식이다.
• 수정된 항의 수 $ a(n) $에 대한 생성함수를 수립하고, $ n = 24 $까지의 값을 계산하는 것.

제안 방법

• 라그랑주 역전개의 한 형태를 사용하여 수정된 공식을 유도하며, $ d^n y / dx^n $를 $ (0,1) \notin p $인 두 차원 분할 $ p \vdash (n, |p|-1) $에 대한 합으로 표현한다.
• 수학적 귀납법을 적용하여 수정된 공식이 차수에 관계없이 일관성을 유지함을 확인한다.
• 두 차원 분할 $ p \vdash (m,n) $에 대한 표기법을 도입하며, 사전순서와 다중도 수 $ e_{p,k,l} $를 포함한다.
• Faà di Bruno의 공식과 유사하게, 각 분할 항의 가중치를 주기 위해 $ \alpha_p = \frac{n! m!}{\prod_i p_{i1}! p_{i2}! \prod_{k,l} e_{p,k,l}!} $를 정의한다.
• 항의 수 $ a(n) $를 계산하기 위해 생성함수 $ \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)\} $이다.
• 로그 생성함수에서 유도된 $ q_m(t) $를 사용하여 $ p_n(t) = \frac{1}{n} \sum_{s=1}^n s q_s(t) p_{n-s}(t) $의 재귀적 계산을 통해 $ a(n) $을 효율적으로 계산한다.

실험 결과

연구 질문

• RQ1 암시적으로 $ F(x,y) = 0 $로 정의된 $ y = f(x) $에 대해 $ n $-번째 도함수 $ d^n y / dx^n $에 대한 올바른 조합적 공식은 무엇인가요?
• RQ2 Comtet와 Fiolet의 원래 공식이 $ d^n y / dx^n $에 대해 오류를 포함하는 이유는 무엇이며, 그 유도 과정에서 오류의 근본 원인은 무엇인가요?
• RQ3 공식의 서로 다른 항의 수 $ a(n) $에 대한 정확한 표현은 무엇이며, 이를 어떻게 계산할 수 있나요?
• RQ4 수정된 공식은 어떻게 라그랑주 역전개와 귀납법을 사용하여 유도될 수 있으며, 두 차원 분할은 이 유도 과정에서 어떤 역할을 하나요?

주요 결과

• 수정된 공식은 $ \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} $의 계수가 아니라.
• $ n = 1 $부터 $ 24 $까지의 수열 $ a(n) $은 $ 1, 3, 9, 24, 61, 145, 333, 732, 1565, 3247, 6583, 13047, 25379, 48477, 91159, 168883, 308736, 557335, 994638, 1755909, 3068960, 5313318, 9118049, 15516710 $으로 계산되며, $ n = 23 $까지는 Comtet와 Fiolet의 값과 일치하지만 $ 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) $에서 유도되며, 이는 $ n = 50 $까지 $ a(n) $를 효율적으로 계산할 수 있게 한다.