حلول عددية للمعادلات التفاضلية

[يسار|180px|تصغير|Illustration of numerical integration for the differential equation <math>y'=y, y(0)=1.</math> باللون الأزرق: the [طريقة أويلر, باللون الأخضر: طريقة النقطة المنتصف, باللون الأحمر: الحلحلة الدقيقة, y=e^t. طول الخطوة هو h=1.0.]]

تعتبر المعادلات التفاضلية من الأدوات الرياضية الهامة في فهم العديد من المسائل الفيزيائية والهندسية والاجتماعية وقد امتدت أهميتها مؤخرا إلى حقول العلوم الاقتصادية وظهر ما يسمى بالنمذجة الرياضية.

المعادلة التفاضلية هي معادلة تحتوي على مشتقات وتفاضلات لبعض الدوال الرياضية وتظهر بشكل متغيرات المعادلة.

## تصنيف المعادلات

# معادلات تفاضلية جزئية
# معادلات تفاضلية عادية : هي معادلة تفاضلية تحتوي على متغير مستقل واحد.

## تقسم مسائل المعادلات التفاضلية العادية إلى
# مسائل قيم ابتدائية.
# مسائل قيم حدية.
#
مسائل القيم الابتدائية:

في هذا النوع من المسائل يكون للمعادلة التفاضلية شرط ابتدائي للمتغيرات المستقلة

الشرط الابتدائي يمثل النقطة الابتدائية للدالة التي تمثل حل المعادلة التفاضلية.
,y' = f(x,y),\,\,\,y(a) =\alpha

فيما يلي بعض الطرق العددية لإيجاد الحل العددي للمعادلة التفاضلية.

قبل استعراض الحل التقريبي للطرق العددية لمسالة القيمة الابتدائية نحتاج لبعض التعاريف والنظريات .

**تعريف**

لتكن f(t,y)داله نقول أن f دالة من النوع ليبتشيز في متغير y على المجموعة D ,

R^2\subset D 

إذا وجد L>0

|f(t,y_1)-f(t,y_2)|\le|y_1-y_2|

كلما أخذنا (t,y_1),(t,y_2) في D , يسمى L ثابت ليبتشيز للدالة f

**مثال**

أثبتي أن الدالة |f(t,y)=t|y تحقق شرط ليبتشيز في الفترة

D =\left\{(t,y):{ 1\le t\le 2},{ -3\le y\le 4}\right\}

**الحل :**

لأي زوج من النقاط (t,y_1),(t,y_2) في D

|f(t,y_1)-f(t,y_2)| =|t|y_1|-t|y_2|| =|t|||y_1|-|y_2||\le 2|y_1-y_2|

إذا ƒ تحقق شرط ليبتشيز L =2

## النظرية (1)

نفرض أن D =\left\{(t,y): { a\le t\le b},{ -\infty

بحيث f(t,y) دالة متصلة على D إذا كانت f دالة تحقق شرط ليبتشيز في متغير y على المجموعة D

فإن مسألة القيمة الابتدائية

y'(t) =f(t,y),\,\,\,{ a\le t\le b},\,\,\,{y(a) = \alpha} 

y(t) لها حل وحيد في الفترة {a\le t\le b}

**مثال **

أستخدم نظرية(1) لإثبات أنه يوجد حل وحيد لمسألة القيمة الابتدائية

y' =1+t\sin(ty),\,\,\,{0\le t\le 2},\,\,\,y(0)=0

**الحل :**

t ثابت نطبق نظرية القيمة المتوسطة للدالة

f(t,y) =1+t\sin(ty)

نجد أنه عندما y , و \xi  عدد موجود في (y_1,y_2)

\frac{f(t,y_1)-f(t,y_2)}{(y_1,y_2)}\ ={\frac{ \partial }{ \partial y}\left(f(t, \xi )\right) =t^2 \cos ( \xi t)}

|f(t,y_2)-f(t,y_1)| =|y_2-y_1||t^2\cos ( \xi t)|\le 4|y_2-y_1|

f تحقق شرط ليبتشيز والثابت هو L=4

وبالإضافة إلى ذلك f(t,y)متصلة عند {0\le t\le 2},\,\,\,{ -\infty

من نظرية (1) إنه يوجد حل وحيد لمسألة القيمة الابتدائية

## الطرق العددية

# طريقة اويلر
# طريقة تايلور

y' =f(x,y),\,\,\,{a\le t\le b},\,\,\,y(a) =\alpha

### طريقة تايلور من الرتب العليا

w_i =\alpha 

w_{i+1} =w_i +h T^{n}(t_i,w_i),\,\,\,{i =0,1,2,...N-1}

 h ={\frac {(b-a)}{n}},\,\,\,w_i \cong y(t_i) 

### طريقة اويلر

هي تايلور من الرتبة الأولى

y' =f(x,y),\,\,\,{a\le t\le b},\,\,\,y(a) =\alpha

y_{i+1}=y_i +hf(t_i,y_i)

h = {\frac {(b-a)}{n}},\,\,\,w_i\approx y(t_i)

## خطأ أويلر

### نظرية

إذا كانت f متصلة وتحقق شرط ليبتشيز على

D =\left\{(t,y):{ a\le t \le b},{ -\infty

ويوجد ثابت M بحيث يحقق

{|y''(t)|\le M},\,\,\,{t\in[a,b]} 

عندما تكون y(t) لها حل وحيد لمسألة القيمة الابتدائية

y' =f(x,y),\,\,\,{a\le t\le b},\,\,\,y(a) =\alpha

نفرض أن { w_1,w_2,w_3,...} بحيث أن الخطأ المتولد من طريقة ايلور لأي عدد موجب N

{\forall i =0,1,....,N}

|(y(t_i)-w_i)|\le {\frac {hM}{2L}}[e^{L(t_i-a)} -1]

## ملاحظات
# الخطأ النسبي لطريقة اويلر ناتج عن اختيار أول حدين من متسلسلة تايلور وحذف باقي الحدود.
# كلما زادت الرتبة في طريقة تايلور فإن الدقة تكون أفضل.

## مراجع

* numerical analysis / Richard . Burden /J.Douglas Faires /ninth edition
* مقدمة في التحليل العددي د.مجدي الطويل

تصنيف:حساب التفاضل والتكامل
تصنيف:معادلات تفاضلية عادية
تصنيف:معادلات تفاضلية عددية