تقريب الخط المماس
الفكرة الهندسية
معادلة تفاضلية عادية 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 (خطوتان):
الخطوة 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 مقابل القيمة الحقيقية 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 }
بالنسبة لـ λ الحقيقية والسالبة (معادلة تفاضلية متحللة مثل dy/dx = -2y): يجب أن تقع hλ في الفترة (-2, 0] على المحور الحقيقي. مع λ = -2 و h = 0.5: hλ = -1. هذا بالضبط على حدود الاستقرار — الطريقة مستقرة بشكل هامشي، مما يفسر الفشل النوعي في المثال السابق.
مع h = 1 و λ = -2: hλ = -2، ما يضعنا خارج منطقة الاستقرار. الحل يتذبذب مع اتساع متزايد.
إيجاد حد الاستقرار
Runge-Kutta 4 (RK4) له منطقة استقرار أكبر من أويلر، وهذا أحد الأسباب لتفضيله لمعظم المشاكل.
بالنسبة لـ λ السالبة الحقيقية، RK4 يسمح بـ hλ حتى حوالي -2.785 على المحور الحقيقي (مقابل حد -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 — توازن المعادلة التفاضلية.
استقرار نقطة ثابتة: إذا كان |g'(y)| < 1، التكرارات القريبة تتقارب نحو y. إذا كان |g'(y*)| > 1، تتباعد.
g'(y) = 1 + h·f'(y). عند نقطة ثابتة y: |1 + h·f'(y)| < 1 للاستقرار.
هذا هو بالضبط شرط استقرار أويلر مع λ = f'(y*) — التخطيط الخطي للمعادلة التفاضلية عند التوازن.
حوض الجذب: مجموعة الشروط الأولية التي تتقارب نحو y* تحت خريطة أويلر. بالنسبة للأنظمة غير الخطية، حد حوض الجذب يحدد أين ستتابع المحاكاة بموثوقية توازن المعادلة التفاضلية مقابل التباعد نحو جاذب آخر.
حلقة المحاكاة هي نظام ديناميكي منفصل. سلوكها النوعي — التقارب، التذبذب، التباعد — يعتمد على حجم الخطوة h بالنسبة لهندسة حقل اتجاه المعادلة التفاضلية.
ربط الهندسة بتصميم المحاكاة
هندسة المحاكاة الرقمية تنقسم إلى ثلاث أسئلة:
1. أين منطقة الاستقرار؟ بالنسبة لأويلر: القرص |1 + hλ| ≤ 1. أكبر بالنسبة لـ RK4. غير محدودة (نصف المستوى الأيسر بالكامل) بالنسبة للطرق الضمنية.
2. أين قيم ذاتية المعادلة التفاضلية؟ القيم الذاتية λ لـ jacobian f في كل نقطة تحدد أي منطقة استقرار يجب أن تحتوي على hλ.
3. أي h يحتفظ بكل hλ داخل المنطقة؟ الحد الأقصى المسموح به h = (نصف قطر منطقة الاستقرار) / max|λ|.
بالنسبة للأنظمة الصلبة: max|λ| ضخم، مما يفرض h صغير جداً بالنسبة للطرق الصريحة. الطرق الضمنية مكلفة لكل خطوة لكنها تسمح بـ h كبير — تقايض تكلفة لكل خطوة مقابل الاستقرار.
رؤية Hamming تُترجم: اختيار الطريقة الرقمية يكود رهان حول هندسة طيف القيم الذاتية للمعادلة التفاضلية. جعل هذا الرهان صريحاً هو أول قرار تصميم في أي محاكاة.