Quadratura de Clenshaw-Curtis: diferència entre les revisions
Contingut suprimit Contingut afegit
m Afegir el nombre d'operacions per calcular els pesos. |
m Avanços en la traducció |
||
Línia 1:
{{traducció|en|Clenshaw-Curtis quadrature|}}
La '''quadratura de Clenshaw–Curtis''' i la '''quadratura de Fejér''' són mètodes d'[[integració numèrica]] o quadratures, basades en l'expansió de l'integrant en termes dels [[polinomis de Txebixev]]. Un resum breu de l'algoritme és el següent: la [[funció]] <math>f(x)</math> que s'ha d'integrar és avaluada als <math>N</math> extrems o arrels dels polinomis de Txebixev i aquests valors es fan servir per construir una aproximació polinòmica de la funció; aquesta és integrada exactament per donar una aproximació de la integral exacta que busquem. El càlcul dels pesos d'integració es pot fer mitjançant una [[Transformada cosinus discreta|DCT]]
<!--▼
==Mètode general==
A simple way of understanding the algorithm is to realize that Clenshaw–Curtis quadrature (proposed by those authors in 1960)<ref name=Clenshaw60>C. W. Clenshaw and A. R. Curtis "[http://www.digizeitschriften.de/resolveppn/GDZPPN001163442 A method for numerical integration on an automatic computer] ''Numerische Mathematik'' '''2''', 197 (1960).</ref> amounts to integrating via a [[change of variables|change of variable]] ''x'' = cos(θ). The algorithm is normally expressed for integration of a function ''f''(''x'') over the interval [-1,1] (any other interval can be obtained by appropriate rescaling). For this integral, we can write:▼
:<math>\int_{-1}^1 f(x)\, dx = \int_0^\pi f(\cos \theta) \sin(\theta)\, d\theta . </math>▼
That is, we have transformed the problem from integrating <math>f(x)</math> to one of integrating <math>f(\cos \theta) \sin \theta</math>. This can be performed if we know the [[cosine series]] for <math>f(\cos \theta)</math>:▼
:<math>f(\cos \theta) = \frac{a_0}{2} + \sum_{k=1}^\infty a_k \cos (k\theta)</math>▼
:<math>\int_0^\pi f(\cos \theta) \sin(\theta)\, d\theta = a_0 + \sum_{k=1}^\infty \frac{2 a_{2k}}{1 - (2k)^2} .</math>▼
▲
▲
▲
de manera que la integral esdevé:
▲
▲<!--
Of course, in order to calculate the cosine series coefficients
:<math>a_k = \frac{2}{\pi} \int_0^\pi f(\cos \theta) \cos(k \theta)\, d\theta</math>
|