Legendre Polynomials and Their Use for Karhunen-Loève Expansion

Fuente: arXiv
Saved in:
Bibliographic Details
Main Author: Béreš, Michal
Format: Preprint
Published: 2025
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866918090991730688
author Béreš, Michal
author_facet Béreš, Michal
contents This paper makes two main contributions. First, we present a pedagogical review of the derivation of the three-term recurrence relation for Legendre polynomials, without relying on the classical Legendre differential equation, Rodrigues' formula, or generating functions. This exposition is designed to be accessible to undergraduate students. Second, we develop a computational framework for Karhunen-Loève expansions of isotropic Gaussian random fields on hyper-rectangular domains. The framework leverages Legendre polynomials and their associated Gaussian quadrature, and it remains efficient even in higher spatial dimensions. A covariance kernel is first approximated by a non-negative mixture of squared-exponentials, obtained via a Newton-optimized fit with a theoretically informed initialization. The resulting separable kernel enables a Legendre-Galerkin discretization in the form of a Kronecker product over single dimensions, with submatrices that exhibit even/odd parity structure. For assembly, we introduce a Duffy-type transformation followed by quadrature. These structural properties significantly reduce both memory usage and arithmetic cost compared to naive approaches. All algorithms and numerical experiments are provided in an open-source repository that reproduces every figure and table in this work.
format Preprint
id arxiv_https___arxiv_org_abs_2507_09825
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Legendre Polynomials and Their Use for Karhunen-Loève Expansion
Béreš, Michal
Computational Engineering, Finance, and Science
Probability
65C05, 86-08, 82-08, 65C60, 60-08
This paper makes two main contributions. First, we present a pedagogical review of the derivation of the three-term recurrence relation for Legendre polynomials, without relying on the classical Legendre differential equation, Rodrigues' formula, or generating functions. This exposition is designed to be accessible to undergraduate students. Second, we develop a computational framework for Karhunen-Loève expansions of isotropic Gaussian random fields on hyper-rectangular domains. The framework leverages Legendre polynomials and their associated Gaussian quadrature, and it remains efficient even in higher spatial dimensions. A covariance kernel is first approximated by a non-negative mixture of squared-exponentials, obtained via a Newton-optimized fit with a theoretically informed initialization. The resulting separable kernel enables a Legendre-Galerkin discretization in the form of a Kronecker product over single dimensions, with submatrices that exhibit even/odd parity structure. For assembly, we introduce a Duffy-type transformation followed by quadrature. These structural properties significantly reduce both memory usage and arithmetic cost compared to naive approaches. All algorithms and numerical experiments are provided in an open-source repository that reproduces every figure and table in this work.
title Legendre Polynomials and Their Use for Karhunen-Loève Expansion
topic Computational Engineering, Finance, and Science
Probability
65C05, 86-08, 82-08, 65C60, 60-08
url https://arxiv.org/abs/2507.09825