¿Cuál es el significado de los polinomios de Legendre? ¿Cuáles son algunas de sus aplicaciones más comunes?

En física, los polinomios de Legendre se ven a menudo cuando se trabaja con expansiones multipolares para potenciales [matemáticos] 1 / r [/ matemáticos], por ejemplo, Coulombic, gravitacional, campos magnéticos.

[matemáticas] \ frac {1} {| \ mathbf {r} – \ mathbf {r} ‘|} = \ sum_ {l = 0} ^ \ infty \ frac {r _ {} ^ {l + 1}} P_l \ bigl (\ hat {\ mathbf {r}} \ cdot \ hat {\ mathbf {r}} ‘\ bigr) [/ math]

donde [math] r _> = \ max (r, r ‘) [/ math] y [math] r _ <= \ min (r, r') [/ math].

Trabajar con expansiones multipolares de este tipo ofrece ventajas teóricas (al truncar la expansión para obtener rápidamente efectos de primer o segundo orden) y permite el cálculo eficiente de cantidades directas o derivadas (por ejemplo, derivadas).

Otra aplicación en física se deriva del hecho de que los polinomios de Legendre son funciones propias del siguiente operador diferencial hermitiano, que a menudo se encuentra al resolver PDE que involucran al laplaciano 3D:

De la definición de un polinomio Legendre:

[matemáticas] \ frac {d} {dx} \ biggl [(1-x ^ 2) \ frac {d} {dx} P (x) \ biggr] = – \ lambda P (x), [/ matemáticas]

donde los valores propios [math] \ lambda [/ math] corresponden a [math] n (n + 1) [/ math].

Una solución [matemática] X [/ matemática] del laplaciano 3D en coordenadas polares satisface:

[matemáticas] \ frac {d} {d \ mu} \ biggl [(1- \ mu ^ 2) \ frac {d X} {d \ mu} \ biggr] + n (n + 1) X = 0 [/ matemáticas]

donde [matemáticas] \ mu = \ cos (\ theta) [/ matemáticas]. Debería ser bastante obvio que [matemáticas] X [/ matemáticas] es un polinomio de Legendre, es decir

[matemáticas] X = P_n (\ mu) = P_n (\ cos \ theta) [/ matemáticas]

El operador laplaciano [matemático] \ Delta f = \ nabla ^ 2 f [/ matemático] es un componente central de muchos problemas físicos, por ejemplo, la ecuación de onda y la ecuación de Schrodinger.

Por algo completamente diferente …

Al intentar encontrar una aproximación polinómica a una función, o ajustar puntos de datos para minimizar el error cuadrado, es decir, el ajuste de mínimos cuadrados, necesitamos resolver una ecuación lineal que involucre la matriz de Hilbert, cuyas entradas son los productos internos del polinomio aproximado. . Por ejemplo, si buscamos coeficientes [matemática] a_0, a_1, \ puntos, a_N [/ matemática] para [matemática] f (x) = a_0 + a_1x + a_2x ^ 2 + \ puntos + a_Nx ^ N [/ matemática] , entonces la matriz de Hilbert correspondiente es

[matemáticas] A_ {m, n} = \ langle x ^ m, x ^ n \ rangle = \ int_0 ^ 1 x ^ {m + n} dx = \ frac {1} {m + n + 1} [/ math ]

El número de condición en esta matriz explota con el número de puntos de interpolación [matemática] N [/ matemática], es decir, el sistema lineal se vuelve cada vez más mal acondicionado.

Una solución a este problema de acondicionamiento es usar una base ortonormal. Convenientemente, los polinomios de Legendre no solo minimizan la norma [matemática] L ^ 2 [/ matemática], sino que también forman una base ortonormal, es decir, la matriz de Hilbert es diagonal y está mucho mejor condicionada

[matemáticas] A_ {m, n} = \ int_0 ^ 1 P_m (x) P_n (x) dx = \ frac {\ delta_ {mn}} {2m + 1}, [/ matemáticas]

dado que el delta de Kronecker hace que los elementos fuera de la diagonal se vuelvan cero.

Resolver un sistema lineal diagonal, por ejemplo, mediante sustitución posterior es más estable numéricamente que intentar resolver un sistema mal acondicionado.

Bueno, la mejor manera de encontrar una respuesta es averiguar por qué los matemáticos y los físicos incluso se molestaron en nombrarlos. ¿Por qué Legendre necesitaba estos polinomios?

Puede encontrar la respuesta en el artículo de Wikipedia: polinomios de Legendre

Esencialmente, aparecen en los coeficientes de la expansión en serie del potencial newtoniano:

[matemática] \ frac {1} {\ | \ boldsymbol {r_1} – [/ matemática] [matemática] \ boldsymbol {r_2} \ |} = \ sum_ {k = 0} ^ \ infty \ frac {\ | \ boldsymbol {r_2} \ | ^ k} {\ | \ boldsymbol {r_1} \ | ^ {k + 1}} P_k (\ cos (\ alpha)) [/ math]

Donde [math] \ alpha [/ math] es el ángulo entre [math] \ boldsymbol {r_1} [/ math] y [math] boldsymbol {r_2} [/ math] y [math] P_k [/ math] es el kth Leyenda polinomial.

También aparecen en otras series de expansión de soluciones de ecuaciones expresadas en coordenadas polares como la ecuación de Laplace o de Schrödinger.

Estas expansiones de series nos permiten hacer buenas aproximaciones numéricas de las funciones subyacentes.

La respuesta corta es no, no puede aumentar la precisión manipulando la ubicación de los puntos de cuadratura.

La precisión del cálculo de un elemento finito depende de cuán pequeñas son las celdas de malla y del orden de los polinomios interpoladores. Una vez que la malla y el orden polinómico están fijos, la matriz de rigidez se define de manera única. El trabajo de la rutina de ensamblaje es calcularlo, generalmente usando cuadratura. Mientras la cuadratura sea lo suficientemente precisa, obtendrá exactamente la misma matriz de rigidez, sea cual sea la fórmula de cuadratura que utilice.

Sin embargo, se ha trabajado en estrategias de posprocesamiento para aumentar la precisión de, por ejemplo, las evaluaciones de estrés. Es decir, en lugar de simplemente leer los valores de puntos de los valores de tensión de la solución final, puede ser mejor formar varios promedios ponderados de la solución en regiones más grandes. Este procedimiento puede aumentar la precisión para que las evaluaciones de tensión sean del mismo orden de precisión que la tensión. Google Babuska, Miller, “El enfoque de procesamiento posterior en el método de elementos finitos” si tiene curiosidad.

Ayuda con la integración numérica.
La integración numérica es útil para la computadora porque una computadora realmente no puede integrar una función. Solo puede evaluarlo en algunos puntos.
La integración de Gauss consiste básicamente en cambiar la integral de una función en una suma de pocos valores de la función, pero esos pocos valores deben elegirse bien. Resulta que si tomas esos valores como las raíces de los polinomios de Legendre, tienes algo que funciona bastante bien.
No sé cuánto se usa en el software moderno para realizar integrales (¿tal vez alguien en quora podría ayudar con esto?), Pero si necesita una forma simple y rápida de hacer una integración numérica, esta es seguramente una buena manera.

Debido a la naturaleza ortogonal de los polinomios de Legendre, son una base conveniente para ciertos cálculos matemáticos. Por ejemplo, en el método discontinuo de Galerkin (un tipo de método de elementos finitos), a menudo se usan como base, ya que simplifican la teoría y el cálculo de las integrales. Por ejemplo, cuando asigna algunas de las integrales definidas al intervalo [-1,1] y utiliza una cuadratura gaussiana para la integración numérica, parte del cálculo puede reducirse a fórmulas muy simples. El método de elementos finitos se utiliza en muchas aplicaciones de ingeniería de la vida real.


Los polinomios de Legendre aparecen en armónicos esféricos, que son una función propia del operador de Laplace sobre la esfera. Esta descomposición se usa realmente en modelos espectrales meteorológicos, por ejemplo, el Sistema de pronóstico global y el Sistema de pronóstico integrado.

Los polinomios de Legendre son las funciones propias latitudinales naturales en una esfera. Por lo tanto, puede descomponer una función en coordenadas ortogonales tal como lo haría con una serie de Fourier. También puede calcular una serie de problemas relacionados con términos polinomiales utilizando los coeficientes de Clebsch-Gordan.

Mientras lucho por encontrar aplicaciones, que no sean artísticas o estadísticas para matemáticas avanzadas, los trinomios son fáciles. Tienes una receta para un pastel de 6X6 pulgadas, pero tu sartén de 6X6 pulgadas se arruinó. Todo lo que tienes es una sartén de 9X12. ¿Cuánto aumentas tu receta para cubrir la sartén? Sustituye x por 6 en X ^ 2 + 9x + 18.

Además, aunque no tengo solicitudes diarias para matemáticas avanzadas, mi esposo, que es estadístico y modelador matemático, tiene muchas.

En general, si está tratando de descubrir cómo cambiar sus dimensiones afectaría el área de una forma bidimensional que puede encontrar polinomios que contienen cuadrados, si está tratando con volúmenes, puede encontrar polinomios que contienen cubos.

Esencialmente, cada vez que desee calcular de manera realista una función determinada, utiliza un polinomio de Taylor (consulte el teorema de Taylor) o alguna aproximación más avanzada que utiliza polinomios o funciones racionales (proporciones de funciones polinómicas).

Los polinomios son geniales porque son fáciles de calcular y con suficientes potencias y los coeficientes correctos, puede hacer que se aproximen a casi cualquier cosa al nivel deseado de precisión.

Bueno, ciertamente los polinomios tienen una aplicación poderosa en el modelado estadístico no lineal. Cuando uno encuentra que el comportamiento de ciertas variables explicativas y una variable de respuesta se explican por una parábola.

Una aplicación interesante, bastante práctica y divertida es quizás usar un polinomio de grado 2 o 3 para calcular el volumen y la superficie de un tazón. Si puede encontrar una función que se ajuste mejor al contorno de la base, tiene un buen comienzo. Encontrará que esa función será un polinomio.

¡¡Buena pregunta!!