Gauss-Legendre Kuadratür Hesaplayıcı nedir?
Bu araç, referans aralığı olan [-1, 1] üzerinde n noktalı Gauss-Legendre kuadratür kuralının düğümlerini (apsisleri) ve ağırlıklarını hesaplar. Gauss-Legendre kuadratürü, belirli bir integrali fonksiyon değerlerinin ağırlıklı toplamı olarak yaklaşık hesaplayan sayısal bir integral yöntemidir: f(x) fonksiyonunun -1'den 1'e integrali, yaklaşık olarak i üzerinden \(w_i\) çarpı \(f(x_i)\) toplamına eşittir. Yalnızca n nokta kullanarak \(2n-1\) dereceye kadar olan tüm polinomları tam olarak integre eder; bu da onu yamuk (trapez) ya da Simpson gibi eşit aralıklı kurallara kıyasla çok daha hassas kılar.
$$\int_{-1}^{1} f(x)\,dx \approx \sum_{i=1}^{n} w_i\, f(x_i)$$Nasıl Kullanılır?
Mertebe n değerini (nokta sayısı, 2 ile 100 arası) seçin, isterseniz görüntülenecek anlamlı basamak sayısını da belirleyin. Hesaplayıcı, her satırında bir \(x_i\) düğümü ve buna karşılık gelen \(w_i\) ağırlığı bulunan n satırlık bir tablo döndürür. Düğümler 0 etrafında simetriktir ve tamamı (-1, 1) aralığının kesinlikle içinde yer alır; ağırlıkların tümü pozitiftir ve toplamları, aralığın uzunluğu olan tam olarak 2'ye eşittir. İntegrali herhangi bir [a, b] aralığı üzerinde almak için her düğümü \(t_i = \frac{b-a}{2} x_i + \frac{a+b}{2}\) dönüşümüyle eşleyin ve her ağırlığı \(\frac{b-a}{2}\) ile çarpın.
Formülün Açıklaması
Düğümler, Bonnet özyinelemesiyle oluşturulan \(P_n(x)\) Legendre polinomunun n köküdür: \(P_0=1\), \(P_1=x\) ve \(P_k = \frac{(2k-1)\, x\, P_{k-1} - (k-1)\, P_{k-2}}{k}\). Her ağırlık $$w_i = \frac{2}{\left(1 - x_i^{2}\right)\left[P_n^{\prime}(x_i)\right]^{2}}$$ formülüyle bulunur; burada türev \(P_n^{\prime}(x) = \frac{n\,(x\, P_n(x) - P_{n-1}(x))}{x^2 - 1}\) şeklindedir. Kökler, \(x = \cos\!\left(\frac{\pi (i - 0.25)}{n + 0.5}\right)\) başlangıç değeriyle başlatılan Newton yöntemiyle bulunur ve bu yöntem birkaç yinelemede yakınsar.
Çözümlü Örnek (n = 3)
\(P_3\)'ün kökleri \(x = 0\) ve \(x = \pm\sqrt{3/5} = \pm 0.7745966692\) değerleridir. \(x = 0\) noktasındaki ağırlık \(\frac{8}{9} = 0.888888889\), \(x = \pm 0.7745966692\) noktalarındaki ağırlık ise \(\frac{5}{9} = 0.555555556\)'dır. Ağırlıkların toplamı $$\frac{5}{9} + \frac{8}{9} + \frac{5}{9} = 2$$ olur ve bu 3 noktalı kural, 5. dereceye kadar olan polinomları tam olarak integre eder.
Sıkça Sorulan Sorular
Ağırlıkların toplamı neden 2'dir? Sabit fonksiyon \(f(x) = 1\)'in [-1, 1] üzerindeki integrali 2'dir ve kuadratür kuralının sabitleri tam olarak üretmesi gerektiğinden ağırlıkların toplamı aralığın uzunluğuna eşit olmalıdır.
Değerler ne kadar hassas? Bu, yaklaşık 15 anlamlı basamak veren çift duyarlıklı (double precision) bir hesaplamadır. Daha fazla basamak isteyen görüntüleme seçenekleri, çift duyarlığın temsil edebildiği değere yuvarlanır.
Tam olarak integre edilen maksimum derece nedir? n noktalı bir kural, \(2n-1\) dereceye kadar olan tüm polinomlar için tamdır.