حلّ مسألة القيمة الابتدائية y′ + 2y/x = x2 لـ x > 0، مع y(1) = 3. أوجد y(2).
y′ + 2y/x = x2y(1) = 3
y(2) =
3المحاليل
1.الطرائق العددية
الإجابة: T = 0,4746
التكامل العددي: قسّم [a, b] إلى n شريحة عرض كل منها h = (b − a)/n بعقد x₀ … xₙ. قاعدة شبه المنحرف: T = (h/2)[f₀ + 2f₁ + 2f₂ + … + 2fn−1 + fₙ]. قاعدة سيمبسون (n زوجي): S = (h/3)[f₀ + 4f₁ + 2f₂ + 4f₃ + … + 4fn−1 + fₙ]. يتناقص خطأ شبه المنحرف مثل h² (الحد (b − a)h²·M/12 مع M = max|f″|)، وخطأ سيمبسون مثل h⁴، لذا يكون سيمبسون عادةً أدق بكثير للعقد نفسها.
العقد والقيم
h = 1,5
Σ wᵢfᵢ = 0 + 2·0,273278 + 0,0862337 = 0,632789
قاعدة شبه المنحرف
T = (1,5/2)·0,632789 = 0,4746
الإجابة
T = 0,4746
ليس لهذا التكامل دالة أصلية أولية، وهذا هو بالضبط الوقت الذي تلزم فيه قاعدة عددية.
2.الطرائق العددية
الإجابة: y4 = −0,8844
الحل العددي لـ y′ = f(t, y): تسير طريقة أويلر على الميل عند بداية كل خطوة: yn+1 = yₙ + h·f(tₙ, yₙ). وتأخذ طريقة هوين (أويلر المحسّنة) متوسط الميل عند البداية وعند تنبؤ أويلر: k₁ = f(tₙ, yₙ)، k₂ = f(tₙ + h, yₙ + h·k₁)، yn+1 = yₙ + (h/2)(k₁ + k₂). وتجمع طريقة رونج–كوتا الكلاسيكية (RK4) أربعة ميول: yn+1 = yₙ + (h/6)(k₁ + 2k₂ + 2k₃ + k₄). الخطأ الكلي لأويلر يتناسب مع h، ولهوين مع h²، ولـ RK4 مع h⁴.
خطوات أويلر
y₁ = −1 + 0,1·(0) = −1
y₂ = −1 + 0,1·(0,2) = −0,98
y₃ = −0,98 + 0,1·(0,392) = −0,9408
y₄ = −0,9408 + 0,1·(0,56448) = −0,884352
الرسم البياني
y4 = −0,8844
الإجابة
y4 = −0,8844
القيمة الدقيقة هي y(0,4) = −0,852144.
3.المتتاليات والمتسلسلات
الإجابة: A: Σ 1/k³
التقارب: إذا لم تؤول الحدود إلى 0 فالمتسلسلة متباعدة (اختبار الحد النوني). Σ1/kᵖ تتقارب فقط عندما p > 1، و Σrᵏ فقط عندما |r| < 1. قارن متسلسلة بإحداهما عندما تسلك حدودها سلوكهما لـ k كبيرة. اختبار النسبة: L = lim |ak+1/ak|؛ اختبار الجذر: L = lim |ak|1/k. L < 1 تتقارب، L > 1 تتباعد، L = 1 لا يحسم شيئًا. متسلسلة القوى Σck(x − x₀)ᵏ تتقارب لـ |x − x₀| < R مع R = lim |ck/ck+1|؛ والمتسلسلة المتناوبة التي تتناقص حدودها إلى 0 تتقارب.
C: Σ (1 + 1/k)ᵏ → تؤول الحدود إلى e لا إلى 0: تتباعد باختبار الحد النوني.
D: Σ (5/4)ᵏ → هندسية مع r = 5/4 > 1: تنمو الحدود، فتتباعد.
الجواب
A: Σ 1/k³
الإجابة
A: Σ 1/k³
4.المتتاليات والمتسلسلات
الإجابة: |R₅| = 0,0709
كثيرات حدود تايلور: قرب x = a تُقرَّب دالة ملساء بـ Pₙ(x) = f(a) + f′(a)(x − a) + f″(a)(x − a)²/2! + … + f⁽ⁿ⁾(a)(x − a)ⁿ/n!؛ وعندما a = 0 تكون هذه كثيرة حدود ماكلورين. المتسلسلات المعيارية: eˣ = Σxᵏ/k!، sin x = x − x³/3! + x⁵/5! − …، cos x = 1 − x²/2! + x⁴/4! − …، ln(1 + x) = x − x²/2 + x³/3 − …، 1/(1 − x) = Σxᵏ. حد الخطأ للاغرانج: |f(x) − Pₙ(x)| ≤ M·|x − a|ⁿ⁺¹/(n + 1)!، حيث M أكبر قيمة لـ |f⁽ⁿ⁺¹⁾| بين a و x.
حد المشتقة
M = e1,5 = 4,482
حد لاغرانج
|R₅| ≤ 4,482·1,5⁶/6! = 0,0709
الإجابة
|R₅| = 0,0709
الخطأ الحقيقي أصغر: |f(1,5) − P₅(1,5)| = 0,01997.
5.المصفوفات والمحددات
الإجابة: c12 = −2
المصفوفات: محدد مصفوفة 2×2 هو ad − bc. وللمصفوفة 3×3، انشر على صف: كل عنصر مضروب في محدد المصفوفة 2×2 المتبقية بعد حذف صفه وعموده، مع إشارات متناوبة + − +. القيم الذاتية λ تحقق det(A − λI) = 0؛ وللمصفوفة 2×2 يكون ذلك λ² − (الأثر)·λ + det = 0.
صف في عمود
c12 = 2·(−1) + 5·0 = −2
الإجابة
c12 = −2
6.المصفوفات والمحددات
الإجابة: λ = 6
المصفوفات: محدد مصفوفة 2×2 هو ad − bc. وللمصفوفة 3×3، انشر على صف: كل عنصر مضروب في محدد المصفوفة 2×2 المتبقية بعد حذف صفه وعموده، مع إشارات متناوبة + − +. القيم الذاتية λ تحقق det(A − λI) = 0؛ وللمصفوفة 2×2 يكون ذلك λ² − (الأثر)·λ + det = 0.
المعادلة المميزة
λ² − (9)λ + (18) = 0
الحل
λ = (9 ± √(9))/2 → λ = 3, 6
λ = 6
الإجابة
λ = 6
7.المعادلات التفاضلية
الإجابة: y(2) = 0,4035
المعادلات الخطية من الرتبة الثانية: للمعادلة y″ + b·y′ + c·y = 0 جرّب y = ert: فتحصل على المعادلة المميِّزة r² + b·r + c = 0. جذران حقيقيان r₁ ≠ r₂ يعطيان y = A·er₁t + B·er₂t (فوق التخميد)؛ وجذر مضاعف r يعطي y = (A + Bt)·ert (تخميد حرج)؛ وجذران مركّبان −α ± jω يعطيان y = e−αt(A cos ωt + B sin ωt) (تحت التخميد، متذبذب). تحدد القيمتان الابتدائيتان y(0) وy′(0) الثابتين A وB. وعند كتابتها y″ + 2ζωny′ + ωn²y = 0 تكون ωn التردد الطبيعي وζ نسبة التخميد: ζ < 1 يعني تذبذبًا. ومع دخل ثابت y″ + by′ + cy = c·K يكون الحل القيمة المستقرة K زائد الاستجابة الحرة.
الجذور المميِّزة
(r + 1)(r + 3) = 0 → r = −1, −3
الثابتان من القيم الابتدائية
A + B = 2, −1A − 3B = 0
A = 3, B = −1
احسب
y(2) = 3·e−2 + (−1)·e−6 = 0,4035
الإجابة
y(2) = 0,4035
8.المعادلات التفاضلية
الإجابة: y(2) = 2,3
المعادلات التفاضلية من الرتبة الأولى: للمعادلة dy/dt = k·y الحل y = y₀·ekt (نمو عندما k > 0، واضمحلال عندما k < 0)؛ وعمر النصف T يعطي k = ln 2/T. قانون نيوتن للتبريد dT/dt = −k·(T − Tₐ) يعطي T(t) = Tₐ + (T₀ − Tₐ)·e−kt. المعادلة الخطية y′ + p(x)·y = q(x) تُحلّ بعامل التكامل μ = e∫p dx: (μy)′ = μq، لذا μy = ∫μq dx + C.