[논문 리뷰] Integral equations, fractional calculus and shift operator

D. Babusci, G. Dattoli외 1명|arXiv (Cornell University)|2010. 07. 29.

Fractional Differential Equations Solutions참고 문헌 4인용 수 4

한 줄 요약

이 논문은 분수 도함수와 일반화된 이동 연산자—특히 $x\partial_x$ 기반 확대 연산자—를 사용하는 새로운 운영 해석학 프레임워크를 개발하여 일반화된 적분 방정식, 즉 일반화된 람브-배트만 및 베크르 유형 방정식을 해결한다. 주요 기여는 $x\partial_x$의 분수 거듭제곱을 통해 표현되는 통합된 방법으로, 특수 함수와 적분 표현을 통해 해석적 해를 도출할 수 있게 한다.

ABSTRACT

We present an extension of a previously developed method employing the formalism of the fractional derivatives to solve new classes of integral equations. This method uses different forms of integral operators that generalizes the exponential shift operator.

연구 동기 및 목표

• 기존의 고전적 람브-배트만 방정식을 초월하여 분수 미적분학과 이동 연산자 체계를 이용해 새로운 종류의 적분 방정식을 해결하는 데 그 목적이 있다.
• 비정수 거듭제곱의 연산자 $x\partial_x$ 를 포함하는 적분 방정식을 체계적으로 해결하는 방법을 개발하는 데 그 목적이 있다.
• 비선형 스케일링 커널을 가진 적분 방정식을 해결하기 위해 지수 이동 연산자를 확대 및 일반화된 기능적 연산자로 일반화하는 데 그 목적이 있다.
• $\Gamma(x\partial_x + 1)$ 과 그 역의 운영 항등식을 수립하여 지수 및 거듭제곱 법칙 커널을 가진 방정식의 해를 구하는 데 그 목적이 있다.
• 비정수 거듭제곱의 $x\partial_x$ 를 위한 스틸링 수 항등식의 분수 확장 형태를 추측하고 수치적으로 검증하는 데 그 목적이 있다.

제안 방법

• 방법론은 $a^{-\nu} = \frac{1}{\Gamma(\nu)}\int_0^\infty e^{-sa} s^{\nu-1} ds$ 를 이용하여 $ (x\partial_x)^{-\nu} $ 를 지수 확대 연산자 $ \hat{E}(-s) = \exp(-s x\partial_x) $ 의 적분으로 표현한다.
• 확대 연산자 $ \hat{E}(\lambda) = \exp(\lambda x\partial_x) $ 는 $ \hat{E}(\lambda)g(x) = g(e^\lambda x) $ 로 작용하여 적분 방정식을 운영 형태로 변환한다.
• 가우시안 또는 지수 커널을 가진 경우, $ \int_0^\infty e^{-y^2} \hat{E}(-y^2) dy $ 와 같은 적분을 평가하여 $ \sqrt{x\partial_x} $ 와 같은 분수 연산자를 도출한다.
• 방정식 $ \int_0^\infty dy\, u(e^{-y^2}x) = f(x) $ 의 해는 $ u(x) = \frac{2}{\sqrt{\pi}} \sqrt{x\partial_x} f(x) $ 로 유도되며, 이는 $ (x\partial_x)^{-1/2} $ 를 통해 적분 형태로 표현된다.
• 비선형 스케일링을 위한 일반화된 이동 연산자 항등식 $ \exp(\lambda q(x)\partial_x)g(x) = g\left[F^{-1}(\lambda + F(x))\right] $ 는 $ F(x) = \int^x d\xi / q(\xi) $ 를 사용하여 유도되며, 비선형 스케일링을 가진 방정식을 해결하는 데 사용된다.
• 이 방법은 $ \int_0^\infty dt\, f(x g(t)) $ 와 같은 형태의 적분으로 확장되며, $ \hat{O}(x\partial_x)f(x) $ 로 재구성되며, 여기서 $ \hat{O}(\mu) = \int_0^\infty [g(t)]^\mu dt $ 이다.

실험 결과

연구 질문

• RQ1 분수 도함수와 일반화된 이동 연산자를 사용하여 비정수 거듭제곱의 $ x\partial_x $ 를 포함하는 적분 방정식을 어떻게 해결할 수 있는가?
• RQ2 커널이 $ e^{-y^2} $, $ e^{-y} $, 또는 $ (1+y^2)^{-\nu} $ 와 같은 경우, 해의 운영 형태는 어떻게 되는가?
• RQ3 비선형 스케일링 변환, 예를 들어 $ F^{-1}(-y^2 + F(x)) $ 를 포함하는 방정식에 대해 일반화된 이동 연산자 체계를 적용할 수 있는가?
• RQ4 지수 커널을 가진 방정식을 해결하기 위해 $ \Gamma(x\partial_x + 1) $ 의 역을 어떻게 운영적으로 표현할 수 있는가?
• RQ5 $ (x\partial_x)^\nu $ 의 구조는 비정수 $ \nu \in (0,1) $ 에 대해 어떻게 되며, 일반화된 스틸링 수를 사용해 $ x^k \partial_x^k $ 로 전개할 수 있는가?

주요 결과

• 수정된 람브-배트만 방정식 $ \int_0^\infty dy\, u(e^{-y^2}x) = f(x) $ 의 해는 $ u(x) = \frac{2}{\sqrt{\pi}} \sqrt{x\partial_x} f(x) $ 로 주어지며, 이는 $ u(x) = \frac{2}{\pi} \int_0^x \frac{\partial_\xi f(\xi)}{\sqrt{\ln(x/\xi)}} d\xi $ 로 단순화된다.
• 방정식 $ \int_0^\infty dy\, e^{-y} u(y x) = f(x) $ 의 해는 $ u(x) = \left\{ \Gamma(x\partial_x + 1) \right\}^{-1} f(x) $ 이며, $ f(x) = e^{-x} $ 일 경우 $ u(x) = J_0(2\sqrt{x}) $, 즉 0차 베셀 함수가 된다.
• 만약 $ f(x) = e^{-x} $ 이고 커널이 $ u(y^m x) $ 라면, 해는 $ u(x) = W_0(x|m) = \sum_{k=0}^\infty \frac{(-1)^k}{k! (m k)!} x^k $ 로 주어지며, 이는 0차 베셀-라이트 함수이다.
• 방정식 $ \int_0^\infty dy\, u(\sqrt{x^2 + 2y^2}) = f(x) $ 의 해는 $ u(x) = -\frac{2}{\sqrt{\pi}} \sqrt{ -\frac{1}{x}\partial_x } f(x) $ 로 유도되며, 항등식 $ \exp(\lambda \frac{1}{x}\partial_x) g(x) = g(\sqrt{x^2 + 2\lambda}) $ 를 사용한다.
• 방정식 $ \int_0^\infty dt\, \frac{1}{(1+t^2)^{\nu x\partial_x}} \sin x $ 의 해는 $ F(x;\nu) = \frac{\sqrt{\pi}}{2} \sum_{n=0}^\infty \frac{(-1)^n}{(2n+1)!} \frac{\Gamma((2n+1)\nu - 1/2)}{\Gamma((2n+1)\nu)} x^{2n+1} $ 로 주어지며, 수치적으로 확인되었다.
• 논문은 $ 0 < \nu < 1 $ 에 대해 $ (x\partial_x)^\nu = \sum_{k=0}^\infty S(k,\nu) x^k \partial_x^k $ 라는 추측을 제시하며, 수치적 증거는 이를 지지하지만 엄밀한 증명은 아직 남아 있다.