تقريب خط الطلاء
الفكرة الجبرية
معادلة التفاضل العددي dy/dx = f(x,y) تمنح كتلة من الميل لكل نقطة في المستوى (x,y) — حقل اتجاه. الحل الصحيح y(x) هو مجرة تتبع دائمًا تلك الميلات المخصصة.
طريقة أيلر تحول الحقل المتناوب ongoing إلى مشي مستمر:
>(xₙ, yₙ) → (xₙ + h, yₙ + h·f(xₙ, yₙ))
من نقطة (xₙ, yₙ)، اذهب مسافة h بتangent line. وصل إلى نقطة تقريبية التالية. كرر.
الخطأ الجبري: خط الطلاء في (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 مقابل القيمة الحقيقية 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. هذا هوexact على الحدود الاستقرار - التفسير النوعي الفاشل في المثال السابق.
مع h = 1 و λ = -2: hλ = -2، مما يجعلنا خارج منطقة الاستقرار. الحل يتردد مع تăng في الأмпليود.
إيجاد الحدود الاستقرار
طريقة رونغ-كوتا 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 relative إلى هندسة مجال الاتجاه للمعادلة التفاضلية.
ربط الهندسة مع تصميم التمثيل المخطط
جبر النمذجة الرقمية تنتهي إلى ثلاث أسئلة:
1. أين منطقة الاستقرار؟ بالنسبة لإولر: الدائرة |1 + hλ| ≤ 1. أكبر بالنسبة لـ RK4. غير محدودة (نصف الكرة الأيسر في المستوى المثلثي) للأساليب الضمنية.
2. أين قيم الأبعاد للمسلمات؟ قيم الأبعاد λ من جاكوبيان f في كل نقطة تحدد منطقة الاستقرار التي يجب أن تحتوي hλ.
3. ما الذي يحافظ على جميع hλ داخل المنطقة؟ الحد الأقصى للسماح ب h = (راديوس منطقة الاستقرار) / max|λ|.
في الأنظمة الصلبة: max|λ| ضخم، مما يضطر إلى استخدام h صغير للأساليب الواضحة. الأساليب الضمنية مكلفة لكل خطوة ولكنها تسمح بالحصول على h كبير - وهي تتبادل تكلفة الخطوة لكل خطوة من أجل الاستقرار.
فكرة هامنج تتtranslate: اختيار طريقة النمذجة يحتوي على تحديق حول الهندسة للمسلمات الموجودة في طيف الأبعاد. جعل هذا التحديق واضحًا هو القرار التصميم الأول في أي نمذجة.