Bu hesaplayıcı ne işe yarar?
Bu araç, n noktalı genelleştirilmiş Gauss-Laguerre kuadratürü için düğümleri (apsisler) \(x_i\) ve ağırlıkları \(w_i\) hesaplar. Saf matematiğe dayanan bir sayısal integrasyon aracıdır ve her yerde aynı şekilde çalışır. Kural, \(x^{\alpha}e^{-x}\) ağırlık fonksiyonunu taşıyan [0, ∞) yarı sonsuz aralığındaki integralleri yaklaşık olarak hesaplar:
$$\int_{0}^{\infty} x^{\alpha} e^{-x} f(x)\, dx \approx \sum_{i=1}^{n} w_i\, f(x_i).$$Düğümler, genelleştirilmiş Laguerre polinomu \(L_{n}^{(\alpha)}(x)\)'in pozitif kökleridir ve kural, f en fazla \(2n-1\) dereceli bir polinom olduğunda tam sonuç verir.
Nasıl kullanılır?
n mertebesini (nokta sayısı, 2 ile 100 arası) seçin, üs parametresi \(\alpha\)'yı girin (−1'den büyük herhangi bir gerçek sayı; klasik Gauss-Laguerre kuralı \(\alpha = 0\) kullanır) ve kaç anlamlı görüntüleme basamağı istediğinizi belirleyin. Sonuç, her düğümü ve ona eşlenen ağırlığı \(x_i\) değerine göre artan sırada listeler; ayrıca yerleşik bir öz-denetim de sunar.
Formül ve yöntem
Her ağırlık şu kapalı biçime uyar: \(w_i = \Gamma(n+\alpha+1)\cdot x_i / (n!\cdot[(n+1)L_{n+1}^{(\alpha)}(x_i)]^{2})\). Dahili olarak buna eşdeğer ve sayısal açıdan kararlı olan Golub-Welsch yöntemini kullanırız: köşegeni \(a_k = 2k+\alpha+1\) ve köşegen dışı elemanları \(b_k = \sqrt{k(k+\alpha)}\) olan simetrik üç köşegenli Jacobi matrisini kurarız. Bu matrisin özdeğerleri düğümleri verir; her ağırlık ise \(\mu_0\cdot(\text{ilk özvektör bileşeni})^2\) değerine eşittir; burada \(\mu_0 = \Gamma(\alpha+1)\) sıfırıncı momenttir. Bu yaklaşım, büyük faktöriyellerden kaynaklanan taşmayı önler.
Çözümlü örnek
\(n = 2\), \(\alpha = 0\) için: \(L_{2}^{(0)}(x) = (x^2-4x+2)/2\) olduğundan kökler \(x = 2 \pm \sqrt{2}\)'dir; bu da \(x_1 = 0{,}5857864\) ve \(x_2 = 3{,}4142136\) verir. Ağırlıklar $$w_1 = \frac{2+\sqrt{2}}{4} = 0{,}8535534 \quad \text{ve} \quad w_2 = \frac{2-\sqrt{2}}{4} = 0{,}1464466.$$ Toplamları \(1 = \Gamma(1)\) olup sonucu doğrular.
Sıkça sorulan sorular
\(\alpha\) ne işe yarar? \(x^{\alpha}\) ağırlığını belirler; \(\alpha = 0\) standart Gauss-Laguerre kuralını verirken \(\alpha > 0\) ağırlığı orijinden uzağa kaydırır. −1'den büyük olmalıdır.
Ne kadar doğru? n noktalı kural, \(2n-1\) dereceye kadar olan polinomları tam olarak integre eder; düzgün (pürüzsüz) fonksiyonlarda yakınsama hızlıdır.
Sonucu nasıl doğrulayabilirim? Tüm ağırlıkların toplamı her zaman \(\Gamma(\alpha+1)\)'e eşittir ve sıfırıncı moment satırında gösterilir.
