Gauss-Kronrod Sayısal İntegral Hesaplama Aracı nedir?
Bu araç, bir \(f(x)\) fonksiyonunun sonlu bir \([a, b]\) aralığındaki belirli integralini Gauss tipi kuadratür yöntemiyle hesaplar. İntegrandı özenle seçilmiş düğüm noktalarında değerlendirir, bu değerleri önceden hesaplanmış ağırlıklarla çarpar ve sonuçları toplar. Yüksek mertebeli \(K\) tahmini, daha düşük mertebeli Gauss-Legendre tahmini \(G\) ile karşılaştırılarak \(|K - G|\) şeklinde bir hata sınırı elde edilir; böylece sonuca ne kadar güvenebileceğinizi bilirsiniz.
Nasıl kullanılır?
x değişkenine bağlı bir matematik ifadesi girin (örneğin 4/(1+x^2) veya sin(x)), alt sınır a ile üst sınır b değerlerini belirleyin ve kural için nokta sayısı n'yi seçin (tek sayı, 3 ile 99 arası). Düzgün integrandlarda daha fazla nokta daha yüksek doğruluk sağlar. Desteklenen söz dizimi: + - * / ^, parantezler ve sin, cos, tan, asin, acos, atan, exp, ln, log, sqrt, abs, sinh, cosh, tanh fonksiyonları ile pi ve e sabitleri.
Formül
Kuadratür kuralları, \([-1, 1]\) aralığında integrand değerlerinin ağırlıklı toplamı olarak tanımlanır. \([a, b]\) aralığında integral almak için afin bir değişken dönüşümü uygulanır: \(x(t) = \frac{b-a}{2} \cdot t + \frac{a+b}{2}\), burada \(dx = \frac{b-a}{2} \cdot dt\) olur. Dolayısıyla integral, \(\frac{b-a}{2}\) ile \(w_i\) ağırlıklarının eşlenen düğümlerdeki \(f\) değerleriyle çarpımının toplamına eşittir.
$$\int_{a}^{b} f(x)\,dx \;\approx\; \frac{b-a}{2}\sum_{i=1}^{n} w_i\, f\!\left(\frac{b-a}{2}x_i + \frac{a+b}{2}\right)$$Gauss-Legendre düğümleri, Legendre polinomunun kökleridir; burada Newton yöntemiyle bulunur ve ağırlıkları \(w_i = \frac{2}{(1 - t_i^2) \cdot P'_m(t_i)^2}\) ile verilir.
Çözümlü örnek
\(f(x) = \sin(x)\) için \([0, \pi]\) aralığında tam integral $$\left[-\cos(x)\right]_{0}^{\pi} = -\cos(\pi) + \cos(0) = 1 + 1 = 2$$'dir. Hesaplama aracı, çok küçük bir tahmini hatayla yaklaşık \(2\) değerini verir. Benzer şekilde \(f(x) = \frac{4}{1+x^2}\) için \([0, 1]\) aralığında sonuç \(\pi = 3{,}14159265\) olur; çünkü \(4\cdot\arctan(1) = \pi\)'dir.
Sıkça sorulan sorular
n neden tek sayı olmalı? Gömülü Gauss-Kronrod çifti, \(m = \frac{n-1}{2}\) adet Gauss düğümünü yeniden kullanır; bu da toplam düğüm sayısının tek olmasını gerektirir.
Hata tahmini ne anlama gelir? Yüksek mertebeli ve düşük mertebeli tahminler arasındaki mutlak farktır; küçük bir değer yakınsamayı gösterir.
Tekillikler için ne yapılır? İntegrallenebilir uç nokta tekillikleri doğruluğu düşürür. Sonlu olmayan değerlendirmeler atlanır ve \(a = b\) olduğunda sonuç tam olarak \(0\) döner.