Basin tarik untuk metode Newton

9

Metode Newton untuk menyelesaikan persamaan nonlinier diketahui konvergen secara kuadrat ketika tebakan awal "cukup dekat" dengan solusinya.

Apa itu "cukup dekat"?

Apakah ada literatur tentang struktur cekungan tarik ini?

David Ketcheson
sumber
Root harus diisolasi (bukan multipel). Jika Hessian jelas seragam (cekung ke atas atau ke bawah) di wilayah tersebut, Anda harus bersikap baik. Tentu saja menjamin atau menguji kondisi ini secara empiris biasanya tidak praktis.
hardmath
Saya melihat pertanyaan di NA-Digest tempo hari dan berpikir itu menarik. Tampaknya saya bukan satu-satunya :-)
Wolfgang Bangerth

Jawaban:

8

Untuk persamaan rasional tunggal dalam domain kompleks, cekungan tarik adalah fraktal, kompelemen dari apa yang disebut set Julia. http://en.wikipedia.org/wiki/Julia_set . Untuk teori dengan beberapa angka online yang bagus, lihat, misalnya,
http://mathlab.mathlab.sunysb.edu/~scott/Papers/Newton/Published.pdf
http://hera.ugr.es/doi/15019160.pdf

x31=0

Dengan demikian ada sedikit gunanya dalam menentukan secara rinci apa yang "cukup dekat" dengan solusi. Jika seseorang tahu batas pada turunan kedua, ada teorema Newton - Kantorovich yang memberikan batas lebih rendah pada jari-jari bola di mana metode Newton bertemu, tetapi kecuali dalam 1D, ini cenderung sangat pesimis.

Batas yang bermanfaat secara komputasi dapat diperoleh dengan menggunakan aritmatika interval; lihat, misalnya, makalah saya
Shen Zuhe dan A. Neumaier, operator Krawczyk dan teorema Kantorovich, J. Math. Anal Appl. 149 (1990), 437-443.
http://www.mat.univie.ac.at/~neum/scan/61.pdf

Arnold Neumaier
sumber
x31=0x>0x>1
1
@hardmath: ya, tetapi persamaan kompleks menjadi dua persamaan nyata dalam 2 variabel, yang mana berlaku sama.
Arnold Neumaier
4

"Cukup dekat" sulit untuk dikarakterisasi mengingat hal itu menimbulkan kelas fraktal . Metode Newton dengan strategi globalisasi seperti pencarian garis dan wilayah kepercayaan memperluas daya tarik. Jika struktur masalah tambahan tersedia, seperti dalam optimasi, asumsi yang diperlukan untuk konvergensi dapat semakin melemah.

Jed Brown
sumber
Hanya untuk rasa ingin tahu, apakah Anda memiliki contoh untuk "Jika struktur masalah tambahan tersedia, seperti dalam pengoptimalan, asumsi yang diperlukan untuk konvergensi dapat semakin melemah."?
vanCompute
@vanCompute Lihat contoh ini untuk contoh terkait dengan optimasi, di mana objek fungsional memberikan informasi yang hilang dalam kondisi optimal urutan pertama. Bentuk lain adalah pengetahuan bahwa kelanjutan tertentu (pseudotransient, parameter, kisi, dll) selalu menyatu, tetapi residu mungkin harus meningkat sebelum mencapai solusi jika seseorang mencoba menyelesaikan masalah secara langsung.
Jed Brown
3

Ada beberapa hasil yang bermanfaat untuk metode Newton yang diterapkan pada polinomial kompleks.

f

r=η2d
ηfdf

Batas eksplisit lainnya diberikan oleh Anthony Manning di Bagaimana memastikan menemukan akar polinomial kompleks menggunakan metode Newton (Teorema 1.2).

Lihat juga Cara menemukan semua akar polinomial kompleks dengan metode Newton oleh Hubbard et al.
Menciptakan. Matematika 146 (2001), no. 1, 1–33. pdf

lhf
sumber
Diadaptasi dari math.stackexchange.com/a/1038487/589 .
lhf