ტანგენტის ხაზის მიახლოება
გეომეტრიული იდეა
ჩვეულებრივი დიფერენციალური განტოლება dy/dx = f(x,y) ენიჭებს დახრილობას (x,y) სიბრტყის ყველა წერტილში — მიმართულების ველი. ჭეშმარიტი ამოხსნა y(x) არის მრუდი, რომელიც ყველგან მოცემულ დახრილობებს მიჰყვება.
ოილერის მეთოდი უწყვეტ მიმართულების ველს გარდაქმნის დისკრეტულ სამოსახლოში:
> (xₙ, yₙ) → (xₙ + h, yₙ + h·f(xₙ, yₙ))
წერტილი (xₙ, yₙ) დან მოძრაობენ მანძილი h ტანგენტის ხაზის გასწვრივ. ჩვენ მივაღწიეთ მიახლოებით შემდეგ წერტილს. გაიმეორეთ.
გეომეტრიული შეცდომა: (xₙ, yₙ) წერტილში ტანგენტის ხაზის დახრილი არის f(xₙ, yₙ), მაგრამ ჭეშმარიტი მრუდი აქვს განსხვავებული დახრილი [xₙ, xₙ + h] ინტერვალის ყველა წერტილში. ოილერის ნაბიჯი იყენებს ის დახრილს მარცხენა დასკვნის წერტილიდან მთელი გამგეობის მანძილზე — 'ის დახრილი რომელი იყო.' შეცდომა ნაბიჯზე იზრდება h² წესით.
დაგროვილი შეცდომა
N ნაბიჯი ფიქსირებული დასკვნის x = a მიმართულებით, h = a/N ხემპის სიტუაციაში:
- ლოკალური შეკვეცის შეცდომა ნაბიჯზე: O(h²)
- ნაბიჯების რაოდენობა: N = a/h
- გლობალური შეცდომა: O(h²) × (a/h) = O(h) — პირველი რიგის სიზუსტე
ოილერის მეთოდი პირველი რიგის: h-ის განახევება განახევებს გლობალურ შეცდომას.
ოილერის მეთოდის გაშვება
განიხილეთ dy/dx = y, საწყის პირობით y(0) = 1. ჭეშმარიტი ამოხსნა: y(x) = eˣ, ამიტომ y(1) = e ≈ 2.71828.
გამოიყენეთ ოილერის მეთოდი h = 0.5 ხემპით, x = 0 დან x = 1 მდე (2 ნაბიჯი):
ნაბიჯი 1: y₁ = y₀ + h·f(x₀, y₀) = 1 + 0.5·(1) = 1.5. ახალი წერტილი: (0.5, 1.5).
ნაბიჯი 2: y₂ = y₁ + h·f(x₁, y₁) = 1.5 + 0.5·(1.5) = 2.25. ახალი წერტილი: (1.0, 2.25).
ოილერი იძლევა 2.25 vs ჭეშმარიტი მნიშვნელობა 2.71828. შეცდომა: 0.468. ფარდობითი შეცდომა: ~17%.
ოილერის სტაბილურობის რეგიონის წარმოშობა
ტესტის განტოლებაზე dy/dx = λy (სადაც λ არის რთული რიცხვი), ოილერის მეთოდი იძლევა:
> yₙ₊₁ = yₙ + h·λ·yₙ = yₙ·(1 + hλ)
გამტენი ფაქტორი ნაბიჯზე: z = 1 + hλ.
სტაბილურობის პირობა: გამოთვლილი ამოხსნა რჩება შებოჭილი თუ და მხოლოდ თუ |z| ≤ 1, ე.ი., |1 + hλ| ≤ 1.
ეს არის გეომეტრიული პირობა რთული hλ-სიბრტყეში: წერტილი hλ უნდა იკრიბებოდეს დისკის შიგნით სახლის რადიუსი 1 ცენტრით (-1, 0) ზე.
ოილერის სტაბილურობის რეგიონი: { hλ ∈ ℂ : |1 + hλ| ≤ 1 }
რეალური, უარყოფითი λ (ქვემოთ მოშვებული ODE როგორც dy/dx = -2y): hλ უნდა იკრიბებოდეს ინტერვალში (-2, 0] რეალურ ღერძზე. λ = -2 & h = 0.5 ხემპით: hλ = -1. ეს არის ზუსტად სტაბილურობის საზღვარი — მეთოდი მარტივი სტაბილური, რომელი აუხსნის ხარისხობრივი მარცხი წინა მაგალითში.
h = 1 & λ = -2 ხემპით: hλ = -2, დაგვიშვებ სტაბილურობის რეგიონის გარეთ. ამოხსნა რხევა გაზდილი ამპლიტუდით.
სტაბილურობის საზღვრის ძებნა
რუნგე-კუტა 4 (RK4) აქვს უფრო დიდი სტაბილურობის რეგიონი ოილერზე, რომელი არის ერთი მიზეზი ის აირჩიეს უმეტესი პრობლემებით.
რეალური უარყოფითი λ, RK4 იძლევა hλ დაქვემდებარებული თითქმის -2.785 რეალურ ღერძზე (vs ოილერი -2 ლიმიტი).
მხიფე განტოლებებისთვის ხელოვნებრივი λ ძალიან განსხვავებული დიდე — ვთქვი λ₁ = -1 & λ₂ = -1000 — სტაბილურობა მოითხოვს hλ₂ დაკვზე სტაბილურობის რეგიონში. RK4 რეალურ ღერძზე: h·(-1000) ≥ -2.785, ამიტომ h ≤ 0.002785.
ეს ნაკადი ნაბიჯი ზომა, დიქტირებული მხიფე ხელოვნებრივი λ₂, გრაფიკი სიმულაცია ძვირი თუნდაც ნელი კომპონენტი λ₁ შეიძლება გამოიყენოს h = 2.
ფიქსირებული წერტილები & მოზიდვის აუზი
ოილერის მეთოდი გამოყენებული dy/dx = f(y) განსაზღვრებს დისკრეტულ რუკას: yₙ₊₁ = g(yₙ) = yₙ + h·f(yₙ).
ფიქსირებული წერტილი ამ რუკის: y ისეთი რომ g(y) = y. ოილერი dy/dx = f(y), ფიქსირებული წერტილები აკმაყოფილებს f(y) = 0 — ODE-ის ბალანსი.
სტაბილურობა ფიქსირებული წერტილი: თუ |g'(y)| < 1, მახლობელი ანესტეზია კრებულია y მდე. თუ |g'(y*)| > 1, ის განსხვავდება.
g'(y) = 1 + h·f'(y). ფიქსირებული წერტილ y: |1 + h·f'(y)| < 1 სტაბილურობა.
ეს არის ზუსტად ოილერი სტაბილურობის პირობა λ = f'(y*) ხემპით — ხელოვნებრივი ODE-ის ბალანსი.
მოზიდვის აუზი: დაყენებული საწყისი პირობები რომელი კრებულია y* უკან ოილერი რუკის. არაწრფივი სისტემებისთვის, აუზის საზღვარი განსაზღვრებს სად სიმულაცია გაიზრდება საიმედოდ რჩება ODE-ის ბალანსი vs განსხვავდება სხვა მოზიდვის ადგილი.
სიმულაციის ციკლი არის დისკრეტული დინამიური სისტემა. მისი ხარისხობრივი ქცევა — კრებულება, რხევა, განსხვავება — დამოკიდებული ნაბიჯის ზომა h მოცემული ნიმუში ODE-ის მიმართულების ველი.
გეომეტრიის დაკავშირება სიმულაციის დიზაინი
გეომეტრია მათემატიკური სიმულაციის ქვემოთ გამოდის სამი კითხვა:
1. სად არის სტაბილურობის რეგიონი? ოილერი: დისკი |1 + hλ| ≤ 1. უფრო დიდი RK4. უბოჭილი (მთელი მარცხენა ნახევარი-სიბრტყე) იმპილიციტი მეთოდებისთვის.
2. სად არის ODE-ის ხელოვნებრივი? ხელოვნებრივი ყოკოლიანი-ის ბალანსი განსაზღვრებს რომელი სტაბილურობის რეგიონი უნდა შეიცავდეს hλ.
3. რა h ინახებს ყველა hλ რეგიონში? მაქსიმალური დაშვებული h = (სტაბილურობის რეგიონი სახელი) / max|λ|.
მხიფე სისტემებისთვის: max|λ| არის ენორმული, ზეწოვა ნაკადი h გამოკრებული მეთოდებისთვის. იმპილიციტი მეთოდი ძვირი ნაბიჯზე მაგრამ დაშვება დიდი h — ის სავაჭრო ნაბიჯ-ღირებულება სტაბილურობა.
ჰამინგის შემოთავაზება გარდაქმნა: სიმუჩიული მეთოდი არჩეული კოდებს რესურსი მოითხოვს ODE-ის ხელოვნებრივი სპექტრი გეომეტრია. რეკა რესურსი ნათელი არის პირველი დიზაინი გადაწყვეტილება რომელი სიმულაცია.