[논문 리뷰] Numerical computation of convolutions in free probability theory

Sheehan Olver, Raj Rao Nadakuditi|arXiv (Cornell University)|2012. 03. 08.

Random Matrices and Applications참고 문헌 26인용 수 15

한 줄 요약

이 논문은 자유 확률 이론에서의 자유 복합(덧셈, 곱셈, 압축)을 계산하기 위한 수치적 방법을 제시한다. 이는 적합한 측도의 역카우치 변환을 활용하여 수행되며, 역수 있는 측도(예: 반원형, 가우시안, 제곱근 감쇠형)에 대해 잘 조절된 바르멘드 시스템을 통해 스펙트럴 정확도를 달성한다. 입력 측도에 대해 광범위한 조건 하에서 수렴이 보장된다.

ABSTRACT

We develop a numerical approach for computing the additive, multiplicative and compressive convolution operations from free probability theory. We utilize the regularity properties of free convolution to identify (pairs of) `admissible' measures whose convolution results in a so-called `invertible measure' which is either a smoothly-decaying measure supported on the entire real line (such as the Gaussian) or square-root decaying measure supported on a compact interval (such as the semi-circle). This class of measures is important because these measures along with their Cauchy transforms can be accurately represented via a Fourier or Chebyshev series expansion, respectively. Thus, knowledge of the functional inverse of their Cauchy transform suffices for numerically recovering the invertible measure via a non-standard yet well-behaved Vandermonde system of equations. We describe explicit algorithms for computing the inverse Cauchy transform alluded to and recovering the associated measure with spectral accuracy. Convergence is guaranteed under broad assumptions on the input measures.

연구 동기 및 목표

• 자유 확률 이론에서 자유 복합(덧셈, 곱셈, 압축)을 계산하기 위한 수치적 프레임워크를 개발하기 위해.
• 기호 계산을 초월한 비대수적 측도에 대한 효율적인 수치 도구의 부족을 해결하기 위해.
• 반원형 및 가우시안과 같은 일반적인 분포를 포함한 자유 복합의 고정밀도 계산을 가능하게 하기 위해.
• 역카우치 변환과 측도 복원이 스펙트럴 정확도로 계산될 수 있는 조건을 설정하기 위해.
• 연구자들과 전문가들이 자유 확률 복합을 수치적으로 탐색할 수 있도록 실용적이고 구현 가능한 알고리즘을 제공하기 위해.

제안 방법

• 이 방법은 자유 복합 연산을 입력 측도의 역카우치 변환의 점별 계산으로 환원한다.
• 역카우치 변환을 정확하게 계산할 수 있는 '적합한' 측도(부드럽게 감쇠하는, 자코비형, 반제곱근형, 점 질량 또는 유한 조합)를 식별한다.
• 역수 있는 측도(스바르츠형, 제곱근 감쇠형 또는 반제곱근/스바르츠형)의 경우, 비표준이지만 잘 조절된 바르멘드 시스템을 통해 측도를 복원한다.
• 역카우치 변환은 수치적 적분과 해석적 계속성을 이용하여 계산되며, 카우치 변환의 정규성을 활용한다.
• 스펙트럴 정확도를 달성하기 위해 역카우치 변환을 푸리에 또는 체비셰프 다항식 전개로 표현한다.
• 바르멘드 시스템의 조건수에 대한 경계와 균일한 적분 점 분포, 그리고 기저 함수의 해석성에 기반하여 수렴이 보장된다.

실험 결과

연구 질문

• RQ1 비대수적 확률 측도에 대해 자유 복합 연산을 스펙트럴 정확도로 수치적으로 계산할 수 있는가?
• RQ2 지지가 유한하거나 무한한 측도에 대해 역카우치 변환을 어떻게 정확하게 계산할 수 있는가?
• RQ3 어떤 측도의 클래스에서 역카우치 변환으로부터 출력 측도를 안정적이고 정확하게 복원할 수 있는가?
• RQ4 측도 복원에 사용되는 바르멘드 시스템이 잘 조절된 상태를 유지하는 조건은 무엇인가?
• RQ5 이 방법은 랜덤 매트릭스 이론에서의 실용적 문제, 예를 들어 랜덤 매트릭스의 합 또는 곱의 극한 스펙트럼 분포 계산에 적용될 수 있는가?

주요 결과

• 이 방법은 넓은 범위의 측도에 대해 자유 덧셈, 곱셈, 압축 복합의 스펙트럴 정확도를 달성할 수 있다.
• 반원형 및 가우시안 분포를 포함한 적합한 측도에 대해, 역카우치 변환을 점별로 고정밀도로 계산할 수 있다.
• 바르멘드 시스템을 통한 측도 복원은 미약한 가정 하에 안정적이고 수렴하며, 조건수는 큰 m에 대해 √n + δ로 유계이다.
• 입력 함수가 해석적이며 f(−1) = 0를 만족할 경우, 최소 제곱 근사의 수렴 속도는 초대수적(super-algebraic)이다. 이는 수정된 바르멘드 시스템에서 성립한다.
• 수치 예제를 통해 반원형과 가우시안의 자유 덧셈 복합, 그리고 가우시안의 자유 압축 복합을 고정밀도로 성공적으로 계산하였다.
• 해석 함수의 트렁케이션 오차 감소와 모어-펜로즈 의사역행렬에 대한 경계를 통해 수렴에 대한 이론적 보장을 확립하였다.