Kn ve Tn İçin İndirgeme Bağıntıları
- Derya PAMUKTULUM
- Kategori: Bilim
- Görüntüleme: 13
Trapez Metodunun Kısa Bir Tarihçesi
Üniversitelerimizin Nümerik Analiz dersinde okutulan Romberg'in yöntemindeki trapez kuralının tarihçesi Newton-Cotes'a kadar gider. Newton-Cotes yöntemine göre f(x)'e n. dereceden polinomal g(x) fonksiyonuyla yapılan yaklaşıma göre şu sonuçlar ortaya çıkar (ki n. dereceden g(x) polinomuyla yaklaşım için f(x)'in üzerinde n+1 tane noktanın alınması gerektiğine dikkat ediniz):
n=1 için Trapez kuralı,
n=2 için Simpson 1/3 kuralı (Kepler kuralı),
n=3 için Simpson 3/8 kuralı,
n=4 için Boole kuralı.
Newton, Roger Cotes hakkında "Eğer o yaşasaydı bazı şeyler hakkında daha fazla bilgi sahibi olabilirdik" demiştir (Bkz. 6.2. Newton Cotes Kuadratür Formülü, S. 74).
İlk kuralda f(x) eğrisinin altında bir trapez (yamuk) çizip, onun alanını bularak f(x)'in altında kalan alana bir yaklaşımda bulunulur. Fakat f(x) üzerindeki noktalar doğru parçalarıyla birleştirildiğinden (ki bu durum g(x)'in lineer bir polinom olmasından kaynaklanır) bu yöntemdeki hata polinomal eğriyle (1<n için g(x) polinomuyla) yapılan yaklaşıma göre daha fazla olmaktadır. Bu nedenle diğer kuralların ortaya çıkması bir zaruret halini almıştır ve 2. kural Johannes Kepler (ki 1615'te şarap fıçılarının hacimlerini ölçmek için kullanmıştı. Bkz. Keplersche Fassregel) ve Thomas Simpson ve 3. kural Thomas Simpson tarafından bulunmuştur.
Bu kurallar hakkında Romberg, gözlemleri bir tablo halinde vermiştir!
Romberg'in Tablosu Hakkında
Romberg, tezinin başında yukarıdaki kurallara ilişkin bir tablo verir. Bu tablodan önce şunları söyler: "En yaygın olarak trapez formülü ve Simpson formülü kullanılır. Gauss, ayrıca yüksek dereceli polinomları Fk noktalarına yerleştirdi ve ilgili katsayılar dizilerini hesapladı. Bunları burada yalnızca 4 ve 8 aralıklarla çoğaltıyoruz (Am häufigsten benutz man die Trapez Formel und die Simpson Formel. Gauss hat auch Polynome höheren Grades durch die Punkte Fk gelegt und die zugehörigen Koeffizientenfolgen berechnet. Wir geben hier nur diejenigen für 4 und 8 Intervalle wieder)":
Romberg'in tablosu, DKNVS 1955. Bu formüllerin Romberg tablosundaki köşegenel elemanlardan geldiklerine dikkat ediniz!
Tablonun altında da şu kritik bilgileri verir: "Gauss formülleri, daha rahatsız edici ağırlıkları nedeniyle nadiren kullanılır. Biraz farklı bir hesaplama yöntemiyle, yukarıda belirtilen çeşitli katsayıları kullanmadan S ve G1'in yanısıra I için daha da iyi bir yaklaşık değer hesaplayabiliriz (Die Gauss’schen Formeln werden, wegen ihrer unbequemeren Gewichte, selten benutzt. Durch eine etwas abgeänderte Berechnungsweise können wir aber ausser S und G1 auch einen noch besseren Näherungswert für I berechnen, ohne die oben genannten verschiedenen Koeffizienten zu benutzen)."
Şimdi Romberg'in tablosunu burada doğrultarak (ki son satırda bir hata olmuş) bu formüllerin, dolayısıyla kuralların ne anlama geldiğini açıklayayım:
1. Bu tablodaki ilk satırdan da görüldüğü üzere trapezin katsayıları "1, 1"dir ve bu formül K0'dır.
2. İkinci satırdaki formül, 3 nokta için Simpson 1/3 kuralıdır ve katsayıları "1, 4, 1"dir (Bkz. (16) no'lu formüle). Bu konuda Romberg, tezinin 2. sayfasının başındaki formüllerden sonra şu müthiş tespiti yapar: "S'nin Simpson S formülüne (Simpson 1/3 kuralındaki formül) göre tam olarak yaklaşıklığı temsil ettiği fark edilmemiş görünür (Man scheint nicht bemerkt zu haben, dass die S genau die Näherung nach der Simpson Formel S darstellen)", Bkz. RİK 4, S. 5. Yani Romberg diyor ki: "Simpson S formülü ya da Simpson 1/3 kuralı, Huygens'in algoritmasının trapez formülüne uygulamasından ibarettir". Bu sonuçla Simpson 1/3 kuralının trapez formülünün Richardson ekstrapolasyonundaki ilk mertebeden elde edilen sonuç olduğunu ve bunun kendisinden önce fark edilmediğini söyler!
3. Üçüncü satırdaki formül 5 nokta için Boole kuralı olarak geçer ve katsayıları "7, 32, 12, 32, 7"dir (Bkz. (18) no'lu formüle). Fakat Romberg, "G1" ile gösterdiği bu formülü "1. Gauss Formülü" olarak anar. Acaba bu formülün Gauss'un olduğuna ilişkin Heidelberg Üniversitesi'nde bir kanıt var mıdır?
4. Romberg, tablosunun son satırında "G11" ile gösterdiği formülün "2. Gauss Formülü" ve katsayılarının "989, 5888, -929, 10496, -4540, 10496, -929, 5888, 989" olduğunu ve 2. sayfadaki formüllerden sonraki 2. paragrafta da şunları söyler: "U4'te F'ye aynı ağırlıkları verdiğimizden, yukarıdaki G11 yaklaşımı Q ve benzeri bir yaklaşım değildir (Die oben angeführte Näherung G11 ist nicht in unseren Näherung Q usw. Enthalten, da wir die F in U4 mit gleichen Gewichten versehen)". Yani Romberg'in, metodundaki yaklaşımlarına dahil etmediği bu formül yerine "217, 1024, 1024, 352, 436, 352, 1024, 1024, 217" katsayıları olan yukarıdaki 4. satırdaki formül gelecektir (ki bunun için her satırdaki formüldeki katsayıların pozitif olduğuna dikkat ediniz). Bununla birlikte, söz konusu bu formül (22) no'lu formülde 9 nokta için 1997'de Ueberhuber tarafından verilmiş olduğu geçer. Yani bunun da Gauss'a ait olduğunun kanıtlanması gerekiyor.
Özetle, bu formüllerden anlaşılacağı üzere f(x)'e diğer yöntemlerdeki gibi polinomal bir eğriyle yaklaşmak yerine uygun bir ekstrapolasyon altında trapez formülüyle yaklaşmak yeterlidir. Bu durum 18 yy.'ın sonlarından 20 yy.'ın ortasına kadar anlaşılamadı. Romberg, 1955'te bunun farkına vardı ve kendi tablosunda örneklerini verdi. 10.02.2020, 01:33'te RİK 3 (ki gerçekte 14.05.2003'te) ve 21.04.2021, 04:45'teki RİK 4'te (ki gerçekte 19.02.2017, 03:09:48'de) ise bu durum tarafımca tamamen aydınlatılmış oldu ve böylece diğer tüm yöntemler tarihin çöp sepetine atılarak orijinale (kaynağa. Bkz. Matrix Reloaded) yani trapez formülüne dönüş sağlanmış oldu!
Şimdi trapez kuralına geçmeden önce Simpson'un kısaca tanıttıktan sonra 576 sayfalık "The Doctrine and Application of Fluxions" adlı 2 ciltlik kitabı hakkındaki gözlemlerimi vermeme izin veriniz.
Thomas Simpson (1710-1761), babası bir dokumacı olduğu için kendisi de başlangıçta dokumacı olarak kendini yetiştirmişken daha sonra kendi kendine Matematik üzerine çalışmıştır. 1750'de 2 cilt "The Doctrine and Application of Fluxions" adlı Hesap (Calculus) kitabının basılmasına rağmen, onun ilgi alanı Olasılık Teorisi olmuştur. Bu kitaplarda matematiğin çeşitli branşlarına ait problemlerin çözümleri vardır. Yukarıda anılan kurallar "SECTION VII of the Use of Fluxions in findings the Areas of Curves (Eğrilerin Alanlarının Bulunmasında Fluksiyonların Kullanımı)" parçasındaki bazı problemlerin çözümünde kullanılmış olduğu teklif edilir. Bu problemlerin çözümlerini ayrıntılı bir şekilde incelemedim ama yüzeysel bir bakışla şunları tespit ettim: Öncelikle bu parçadaki (S. 121-155) problemlerde geçen eğrilerin hepsi özeldir, dolayısıyla ilgili problemlerin çözümlerinden trapez kuralını (ve bu arada orta nokta formülünü) ve adıyla anılan kuralları çıkartabilmek mümkün değildir. Bunları bırakın, günümüzdeki koordinat bilgisini bulabilmek bile mümkün değildir. Bu nedenle yukarıda anılan 2 ve 3. kuralların "Simpson" adıyla anılmasını doğru değil; olsa olsa öncüsü olabilir. Çünkü eğer biz, bu kurallara "Simpson" dersek, o zaman antik dönemdekiler mezarlarından kalkarlar ve her şeyi alıp götürürler, elimizde Matematik adına bir şey kalmaz, ya! Örneğin bu kuralları "Simpson" adıyla anabilmemiz için "1.6.1. Romberg Metodu'nun Deşifrasyonu"ndaki gibi bir ispata ihtiyaç vardır. Biz, işte bu yüzden "Romberg İntegrali Yöntemi" adıyla anıyoruz. Neden? Çünkü kapı gibi bir ispat var ve bu, günümüzdekini karşılar. Bu arada, Romberg, 94 yaşına kadar yaşamasına rağmen, dolayısıyla tezini her ne şekilde anlattıysa Almanlar'ın elinde böyle bir ispat var mıydı? Yani Romberg'in tezindeki metodun günümüzdekini karşıladığını biliyorlar mıydı? Hiç zannetmiyorum. Bence her 2 taraf da açmazdaydı. Ama bu ispatı o amaçla yapmadım; sadece Romberg'in yöntemini günümüze taşımak istedim ve taşıdım. Bu sonuçla Romberg'in adını günümüzde de tescilleyerek 1955'teki patent hakkının geçerli olduğunu ispatlamış (göstermiş) olduk!
İlk İndirgeme Bağıntısının Genelleştirilmesi Hakkında
Bilindiği üzere bileşik trapez yöntemindeki trapez ve orta nokta formülleri Önsöz'deki (3) ve (5)'teki şekliyle 19. yy.'dan günümüze kadar hiç değişmeden geldi. Bu formüller Sheppard tarafından 1900'de verilen (11) ve (23)'teki formüller olup, bunlar arasındaki aritmetik ortalama bağıntıyı da (24)'te vermiştir (Bkz. Survey of Extrapolation Processes In Numerical Analysis). İşte Sheppard'ın trapez ve orta nokta formülleri için verdiği (11) ve (23) no'lu formüller günümüze kadar hiç değişmeden geldi ve kimse de bunları değiştirmeye yeltenmedi. Bu bakımdan ilk indirgeme bağıntısı olan (1) no'lu formül bir devrimdir. Bu formülü daha projeye başlamadan ilk gözlem anımda görmüştüm (bkz. RİK 1, (16) no'lu formül) ve ilk kez 1. Yol'da ve 2.2. Yol'da da Teorem 1.3'teki diğer bağıntıyı görmüştüm (18.11.2016, 17:00). Fakat bunlardan ilk indirgeme bağıntısı Sheppard'ın (11)'deki trapez formülünden görülmesine karşın şimdiye kadar "indirgeme" anlamında hiç gösterilmedi ve kullanılmadı!
Peki bu formülü indirgemeli şekilde gösterince ne kazanmıştım?
Ben sadece projenin başında bir nümerik gözlem olarak not etmiştim bu bağıntıyı. Fakat 25.04.2021, 14:29:46'da (9) no'lu formülü görünce, eğer A-ha'dan Pål Waaktaar-Savoy gibi ifade edersem, her şey ilginç olmaya başladı. Çünkü (9)'daki 2. terim yine K0 idi ve bu, (1)'deki ilk bulduğum bağıntıyla ilişkiliydi ve bana yeni bir kapı açmıştı. Hemen başladım, Pål gibi (9) nasıl çalışıyor, (1)'dekinden farkı ne? gibi sorularının yanıtlarını aramaya. Araştırmamdan çıkan sonuca göre (1) ile (9) formüllerinin birbirinden farklı ama aynı sonuçları ürettiklerini, dolayısıyla (9)'un (1)'dekinden farklı bir formda olduğunu gördüm. Tabii o zaman (9) no'lu formüle "Kn'nin Simetrik Açılımı" diyordum!
Şimdi çok sıkı durun: A-ha'nın söz yazarı ve gitaristi olan Pål'un "Düşünce Treni (Train of Thought)" adlı parçasında Knut Hamsun'u andığını biliyordum ama hayatını ve şarkılarını anlattığı"Tears from a Stone" adlı kitabında anne ve babasının Knut Hamsun'un kitaplarını daima okuduklarını ve Pål'un de onların etkisinde kalarak bu parçayı yazmış olabileceği aklımın ucundan bile geçmiyordu. Anladığım kadarıyla anne ve babası, yaşamlarını Hamsun'un kitaplarıyla doldurmuşlar ve buradan yarınki Babalar Günü nedeniyle Pål'un ve tabii ki herkesinkini şimdiden kutlarım. Ama en çok annesini sevdiğini biliyorum (ki soyadı annesinden gelir).
Bu arada hemen bir karşılaştırma yapayım. A-ha'yı ilk çıktığı günden yani liseden beri bilirim. O zamanlar gençtik ve dinliyorduk, ama şimdi mümkün değil. Çünkü müzik anlayışı, daha doğrusu her alanda anlayışlar değişti. Yine de nostaljik olarak bakarım. Aradan geçen 40 yıla yakın sürede bende şu sonuç çıktı: Tamam, zamanla her şey değişiyor (Bkz. "Bir Entellektüel Olarak Mustafa Kemal Atatürk" makalesindeki "Zaman süratle ilerliyor..." ile başlayan sözüne). Örneğin 80'lerdeki müziğin günümüzde bir anlamı yok ama, Ole Amble'ın 1952'deki ve Romberg'in 1955'teki tezlerinde beni çeken şey ne idi? Tek kelimeyle yanıt vermem gerekirse, bu çalışmaların birer klasik olmasıdır. Yani bu çalışmalarda hala ilk günkü heyecanı hissedebiliyorsunuz, ama diğer şeylerde öyle olmuyor. İşte Hamsun'un neden bu kadar çok sevildiği tam bu noktada ortaya çıkıyor.
Daha sonra Kn için bulduğum formülleri (10)'da toplayıp karşılaştırdım ve bunlardan hareketle (1.14)'teki formüllere göre Tn için (11)'deki formülleri buldum. Ama diğer taraftan (9)'daki yöntemi (2) no'lu orta nokta formülünde kullanınca Tn için (12)'deki formül ortaya çıktı. İşte 5 gündür bu formüllerle oyalandım, çünkü arayışımın bittiğini zannediyordum. Bu formüller Kn için (1)'deki ve Tn için (2)'deki formüllerle aynı sonuçları veriyordu ama formları farklıydı!
Pål'un dediği gibi Kn ve Tn için farklı formlar söz konusu olduğunda, bunların o formlarda yazılmasındaki matematiği tanımam gerekiyordu ve her şey bu noktadan itibaren ilginçleşmeye başladı. Çünkü 30.04.2021, 19:55:47'de Kn yi (1.13)'te 2. kez genelleştirdim ve arkasından da (1.17) geldi. Bundan sonrası çorap söküğü gibi geldi ve Kn ve Tn'yi (1.18)&(1.19)'da tam anlamıyla genelleştirmiş oldum. Bu, Teorem 1.2'dir. Bunun arkasından 3 teorem daha verdim. Tabii ki hepsinin ispatlarını ve sonuçlarını en ince ayrıntısına kadar verdim.
Norveç'in İlk Veri Yöneticisi: Ole Amble
Fakat hiçbirinde Önsöz ve EK 1'deki kadar heyecan duymadım. Çünkü bu teoremler yeni bir şey keşfetmenin ötesine gitmedi. Ama içinde insan hikayesi barındıran çalışmalar ister istemez sizi etkiliyor, kayıtsız kalamıyorsunuz. Norveç'in Bilgi İşlem Sistemi'nin kurucusu ve ilk Veri Yöneticisi olan Ole Amble böyle biridir. O, "Ole Amble Algoritması, 01.02.2019" adlı makalemin 100. sayfasında gördüğünüz gibi babasının Tröndheim'daki ölümü nedeniyle oradaki Norveç Bilim ve Teknoloji Üniversitesi'nde öğretim görevlisi olarak çalışır. Sonra Norveç dijitalleşmeye başlayınca Oslo Üniversitesi'nin Matematik Bölümü'ndeki Veri İşleme Bölümü'ne geçer ve emekli olana kadar da orada çalışır. Hayat hikayesi bayağı acıklıdır. Pål'unkü de öyledir. Ne yazık ki Ole Amble'ın 1952'deki tezi için Byrunjulf Owren ile birlikte araştırma yaptık ama bulamadık. Bu nedenle 16. sayfasında (PDF'de 26) Norveç Kralı V. HARALD ekselanslarına hem taziyemi sundum, hem de himayesindeki DKNVS'deki bu tezin akıbeti için bir çağrıda bulundum. Çünkü Ole Amble sıradan biri değildi. Bilmiyorum, ben buradan Ole Amble'ın tezi için elimden geleni yaptım ve bunları burada makaleler halinde yayımlıyorum. Bakalım, Norveç tarafında ne yapılabilir? Ama Ole Amble'ın tezini bulamasalar bile makalelerimdeki çalışmaların yeterli olduğunu ve hatta ötesine geçtiğimi bildiririm. Yani aslında Norveç Kralı'na taziyemi sunarken asıl amacım Ole Amble'ın tezini istemek değildi, çünkü makalelerimde o tezin ötesine geçtik bile; Türk-Norveç Dostluğu'nda bir adım atarak 2 ülke arasındaki ilişkilerin geliştirilmesinden ibarettir. Bu arada Norveç Kralı'nın çok iyi bir okuyucu olduğunu belirtmem gerekiyor. Buradan Norveç Kralı'nın Babalar Günü'nü kutlar ve bu vesileyle acı kaybından dolayı bir kez daha başsağlığı dilerim!
Şimdi sözü fazla uzatmadan hemen çalışmalara geçelim:
1. Kn ve Tn İçin İndirgeme Bağıntıları, 18.06.2021, 23:11,
1.1. Bileşik trapez yöntemine ait tam boy Geogebra grafiği,
2. Ole Amble Algoritması, 01.02.2019, 20:21,
3. Mathematica'da Lineer E-ATA M Algoritması Üzerindeki Ole Amble Algoritmasının Romberg Örneğine Uygulanmasıyla Elde Edilen Yaklaşıklıklara Sembolik Bir Yaklaşım:
3.1. HTML dosyası,
3.2. Mathematica dosyası.