Keplers denklemi - Keplers equation - Wikipedia

0 ile 1 arasındaki beş farklı eksantriklik için Kepler'in denklem çözümleri

İçinde yörünge mekaniği, Kepler denklemi bir cismin yörüngesinin çeşitli geometrik özelliklerini bir merkezi kuvvet.

İlk olarak tarafından türetildi Johannes Kepler 1609'da Bölüm 60'ta Astronomia Nova,[1][2] ve onun V. kitabında Kopernik Astronomisinin Özü (1621) Kepler denkleme yinelemeli bir çözüm önerdi.[3][4] Denklem, özellikle klasik olmak üzere hem fizik hem de matematik tarihinde önemli bir rol oynamıştır. gök mekaniği.

Denklem

Kepler denklemi dır-dir

nerede M ... anomali demek, E ... eksantrik anormallik, ve e ... eksantriklik.

'Eksantrik anormallik' E Kepler yörüngesinde hareket eden bir noktanın konumunu hesaplamak için kullanışlıdır. Örneğin, vücut periastronu koordinatlarda geçerse x = a(1 − e), y = 0, zamanda t = t0, sonra herhangi bir zamanda vücudun konumunu bulmak için önce ortalama anormalliği hesaplarsınız M zamandan ve ortalama hareket n formülle M = n(tt0), sonra yukarıdaki Kepler denklemini çözerek Eardından koordinatları şuradan alın:

nerede a ... yarı büyük eksen, b yarı küçük eksen.

Kepler'in denklemi bir aşkın denklem Çünkü sinüs bir aşkın işlev, bunun için çözülemeyeceği anlamına gelir E cebirsel olarak. Sayısal analiz ve dizi genişletmeler genellikle değerlendirmek için gereklidir E.

alternatif formlar

Kepler denkleminin birkaç biçimi vardır. Her form belirli bir yörünge türü ile ilişkilendirilir. Standart Kepler denklemi, eliptik yörüngeler için kullanılır (0 ≤ e <1). Hiperbolik Kepler denklemi, hiperbolik yörüngeler için kullanılır (e > 1). Radyal Kepler denklemi doğrusal (radyal) yörüngeler için kullanılır (e = 1). Barker denklemi parabolik yörüngeler için kullanılır (e = 1).

Ne zaman e = 0, yörünge daireseldir. Artan e dairenin eliptik olmasına neden olur. Ne zaman e = 1, üç olasılık vardır:

  • parabolik bir yörünge,
  • çekim merkezinden çıkan sonsuz bir ışın boyunca giren veya çıkan bir yörünge,
  • veya çekim merkezinden uzaktaki bir noktaya doğru bir çizgi parçası boyunca ileri geri giden bir yörünge.

Hafif bir artış e 1'in yukarısı, 180 derecenin biraz altında dönüş açısına sahip hiperbolik bir yörünge ile sonuçlanır. Daha fazla artış, dönüş açısını azaltır ve e sonsuza gider, yörünge sonsuz uzunlukta düz bir çizgi olur.

Hiperbolik Kepler denklemi

Hiperbolik Kepler denklemi:

nerede H hiperbolik eksantrik anomalidir. Bu denklem, M'nin yeniden tanımlanmasıyla elde edilir. −1'in karekökü eliptik denklemin sağ tarafının çarpımı:

(içinde E artık hayalidir) ve sonra E tarafından iH.

Radyal Kepler denklemi

Radyal Kepler denklemi:

nerede t zamanla orantılıdır ve x ışın boyunca çekim merkezine olan mesafeyle orantılıdır. Bu denklem, Kepler'in denklemini 1/2 ile çarparak ve e 1'e:

ve sonra ikame yapmak

Ters problem

Hesaplanıyor M belirli bir değer için E basittir. Ancak, çözme E ne zaman M verilmesi çok daha zor olabilir. Yok kapalı form çözümü.

Biri yazabilir sonsuz seriler kullanarak Kepler denkleminin çözümü için ifade Lagrange ters çevirme, ancak dizi tüm kombinasyonları için yakınsamıyor e ve M (aşağıya bakınız).

Kepler denkleminin çözülebilirliği konusundaki kafa karışıklığı, dört yüzyıldır literatürde devam ediyor.[5] Kepler'in kendisi genel bir çözüm bulma olasılığından şüphe duyduğunu ifade etti:

Arkın ve sinüsün farklı doğası nedeniyle [Kepler denkleminin] a priori çözülemeyeceğine yeterince tatmin oldum. Ama yanılıyorsam ve biri bana yolu gösterirse, o benim gözümde büyük olacak Apollonius.

— Johannes Kepler[6]

Ters Kepler denklemi

Ters Kepler denklemi, Kepler denkleminin tüm gerçek değerleri için çözümüdür. :

Bu verimi değerlendirmek:

Bu seriler şu şekilde çoğaltılabilir: Mathematica InverseSeries işlemi ile.

InverseSeries[Dizi[M-Günah[M],{M,0,10}]]
InverseSeries[Dizi[M-eGünah[M],{M,0,10}]]

Bu işlevler basittir Maclaurin serisi. Transandantal fonksiyonların bu tür Taylor serisi temsilleri, bu fonksiyonların tanımları olarak kabul edilir. Bu nedenle, bu çözüm ters Kepler denkleminin biçimsel bir tanımıdır. Ancak, E değil tüm işlev nın-nin M belirli bir sıfır olmayan e. Türev

sonsuz karmaşık sayılar kümesinde sıfıra gider e<1. Çözümler var ve bu değerlerde

(ters cosh pozitif olarak kabul edildiğinde) ve dE/dM bu noktalarda sonsuza gider. Bu, Maclaurin serisinin yakınsama yarıçapının ve dizi değerleri için yakınsamayacak M bundan daha büyük. Seri, hiperbolik durum için de kullanılabilir, bu durumda yakınsama yarıçapı şu şekildedir: Ne zaman serisi e = 1 olduğunda yakınsar m <2π.

Bu çözüm, matematiksel anlamda en basit olanı olsa da,[hangi? ], çoğu uygulama için diğer çözümler tercih edilir. Alternatif olarak, Kepler'in denklemi sayısal olarak çözülebilir.

İçin çözüm e ≠ 1, tarafından bulundu Karl Stumpff 1968'de[7] ancak önemi anlaşılmadı.[8][açıklama gerekli ]

Bir Maclaurin serisi de yazılabilir. e. Bu seri ne zaman birleşmez e daha büyük Laplace sınırı (yaklaşık 0,66), değerine bakılmaksızın M (sürece M katları ), ancak hepsi için birleşir M Eğer e Laplace sınırının altında. Birincisi dışındaki serideki katsayılar (basitçe M), bağlıdır M dönem ile periyodik olarak .

Ters radyal Kepler denklemi

Ters radyal Kepler denklemi (e = 1) şu şekilde de yazılabilir:

Bu verimi değerlendirmek:

Bu sonucu kullanarak elde etmek için Mathematica:

InverseSeries[Dizi[ArcSin[Sqrt[t]]-Sqrt[(1-t)t],{t,0,15}]]

Ters problemin sayısal yaklaşımı

Çoğu uygulama için, ters problem sayısal olarak hesaplanabilir. kök fonksiyonun:

Bu, aracılığıyla yinelemeli olarak yapılabilir Newton yöntemi:

Bunu not et E ve M Bu hesaplamada radyan birimleri cinsindendir. Bu yineleme, istenen doğruluk elde edilinceye kadar tekrarlanır (ör. f(E) E0 = M(t) yeterlidir. Olan yörüngeler için e > 0.8, başlangıç ​​değeri E0 = π kullanılmalıdır. Eğer e aynı 1, sonra türevi fNewton yönteminin paydasında bulunan, sıfıra yaklaşarak Newton-Raphson, secant veya regula falsi gibi türev tabanlı yöntemleri sayısal olarak kararsız hale getirebilir. Bu durumda, ikiye bölme yöntemi, özellikle çözüm küçük bir başlangıç ​​aralığında sınırlanabildiğinden, garantili yakınsama sağlayacaktır. Modern bilgisayarlarda, 17 ila 18 yinelemede 4 veya 5 basamaklı doğruluk elde etmek mümkündür.[9] Kepler denkleminin hiperbolik formu için benzer bir yaklaşım kullanılabilir.[10]:66–67 Parabolik bir yörünge olması durumunda, Barker denklemi kullanıldı.

Sabit nokta yineleme

İlgili bir yöntem şunu belirterek başlar: . Sağdaki ifadeyi tekrar tekrar yerine koymak sağdaki basit bir sabit nokta yineleme değerlendirme algoritması . Bu yöntem, Kepler'in 1621 çözümüyle aynıdır.[4]

işleviE(e,M,n)E=Miçink=1-enE=M+e*günahESonrakikdönüşE

Yineleme sayısı, değerine bağlıdır . Hiperbolik form benzer şekilde .

Bu yöntem, Newton yöntemi yukarıdaki çözüm

Küçük miktarlarda ilk sipariş ve ,

.

Ayrıca bakınız

Referanslar

  1. ^ Kepler, Johannes (1609). "LX. Methodus, ex hac Physica, hoc est genuina & verissima hypothesi, extruendi utramque partem æquationis, & distantias genuinas: quorum utrumque simul per vicariam fieri hactenus non potuit. Argumentum falsæ hypotheseos". Astronomia Nova Aitiologētos, Seu Physica Coelestis, tradita commentariis De Motibus Stellæ Martis, Ex observationibus G.V. Tychonis Brahe (Latince). s. 299–300.
  2. ^ Aaboe, Asger (2001). Astronominin Erken Tarihinden Bölümler. Springer. s. 146–147. ISBN  978-0-387-95136-2.
  3. ^ Kepler, Johannes (1621). "Libri V. Pars altera." Epitome astronomiæ Copernicanæ usitatâ formâ Quæstionum & Responsionum conscripta, inq; VII. Libros digesta, quorum tres hi priores sunt de Doctrina Sphæricâ (Latince). s. 695–696.
  4. ^ a b Swerdlow, Noel M. (2000). "Kepler'in Kepler Denklemine Yinelemeli Çözümü". Astronomi Tarihi Dergisi. 31: 339–341. Bibcode:2000JHA .... 31..339S. doi:10.1177/002182860003100404.
  5. ^ Kepler'in denkleminin "analitik olarak çözülemeyeceği" sıklıkla iddia edilir; örneğin bakınız İşte. Bunun doğru olup olmadığı, sonsuz bir serinin (veya her zaman yakınsamayan) analitik bir çözüm olarak kabul edilmesine bağlıdır. Diğer yazarlar bunun çözülemeyeceği saçma iddiasında bulunurlar; örneğin bkz. Madabushi V. K. Chari; Sheppard Joel Salon; Elektromanyetizmada Sayısal Yöntemler, Academic Press, San Diego, CA, ABD, 2000, ISBN  0-12-615760-X, s. 659
  6. ^ "Mihi ſufficit credere, ſolvi a priori non poſſe, propter arcus ve ſinus ετερογενειαν. Erranti mihi, quicumque viam monſtraverit, erit mihi magnus Apollonius'dur." Hall, Asaph (Mayıs 1883). "Kepler'in Sorunu". Matematik Yıllıkları. 10 (3): 65–66. doi:10.2307/2635832.
  7. ^ Stumpff, Karl (1 Haziran 1968). "Lie serisinin gök mekaniğinin problemlerine uygulanması üzerine". NASA Teknik Notu D-4460. Alıntı dergisi gerektirir | günlük = (Yardım)
  8. ^ Colwell, Peter (1993). Üç Yüzyılda Kepler Denklemini Çözme. Willmann-Bell. s. 43. ISBN  0-943396-40-9.
  9. ^ Keister, Adrian. "Dairesel Bir Parçanın Yüksekliğini Bulmanın Sayısal Analizi". Wineman Teknolojisi. Wineman Technology, Inc. Alındı 28 Aralık 2019.
  10. ^ Pfleger, Thomas; Montenbruck, Oliver (1998). Kişisel Bilgisayarda Astronomi (Üçüncü baskı). Berlin, Heidelberg: Springer. ISBN  978-3-662-03349-4.

Dış bağlantılar