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

حركة المقذوفات المكافئة موضحة متجه السرعة
يتبع مسار المقذوف الذي يتم إطلاقه من المدفع منحنى يتم تحديده بواسطة معادلة تفاضلية عادية مشتقة من قانون نيوتن الثاني.

في الرياضيات، المعادلة التفاضلية العادية ( ODE ) هي معادلة تفاضلية (DE) تعتمد على متغير مستقل واحد فقط . وكما هو الحال مع أي معادلة تفاضلية أخرى، فإن مجهولاتها تتكون من دالة واحدة (أو أكثر) وتتضمن مشتقات تلك الدوال. [ 1 ] يُستخدم مصطلح "عادية" تمييزًا لها عن المعادلات التفاضلية الجزئية (PDEs) التي قد تعتمد على أكثر من متغير مستقل، [ 2 ] وبشكل أقل شيوعًا، تمييزًا لها عن المعادلات التفاضلية العشوائية (SDEs) حيث تكون العملية المُمَثَّلة عشوائية. [ 3 ]

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

المعادلة التفاضلية الخطية هي معادلة تفاضلية تُعرَّف بواسطة متعددة حدود خطية في الدالة المجهولة ومشتقاتها، أي معادلة على الصورة التالية:

أ0(x)y+أ1(x)y+أ2(x)y"++أن(x)y(ن)+ب(x)=0،{\displaystyle a_{0}(x)y+a_{1}(x)y'+a_{2}(x)y''+\cdots +a_{n}(x)y^{(n)}+b(x)=0,}

أينأ0(x)،...،أن(x){\displaystyle a_{0}(x),\ldots ,a_{n}(x)}وب(x){\displaystyle b(x)}هي دوال قابلة للتفاضل بشكل عشوائي ، ولا يشترط أن تكون خطية، و y،...،y(ن){\displaystyle y',\ldots ,y^{(n)}}هي المشتقات المتتالية للدالة المجهولةy{\displaystyle y}من المتغيرx{\displaystyle x}[ 4 ]

من بين المعادلات التفاضلية العادية، تحتل المعادلات التفاضلية الخطية مكانة بارزة لعدة أسباب. فمعظم الدوال الأساسية والخاصة التي تُصادف في الفيزياء والرياضيات التطبيقية هي حلول لمعادلات تفاضلية خطية (انظر: الدالة الهولونومية ). وعندما تُنمذج الظواهر الفيزيائية بمعادلات غير خطية، تُقارب عادةً بمعادلات تفاضلية خطية لتسهيل الحل. أما المعادلات التفاضلية العادية غير الخطية القليلة التي يمكن حلها صراحةً، فتُحل عادةً بتحويل المعادلة إلى معادلة تفاضلية عادية خطية مكافئة (انظر، على سبيل المثال، معادلة ريكاتي ). [ 5 ]

يمكن حل بعض المعادلات التفاضلية العادية بشكل صريح بدلالة الدوال والتكاملات المعروفة . وعندما يتعذر ذلك، قد يكون من المفيد استخدام معادلة حساب متسلسلة تايلور للحلول. أما بالنسبة للمسائل التطبيقية، فيمكن للطرق العددية للمعادلات التفاضلية العادية أن توفر تقريبًا للحل. [ 6 ]

خلفية

تظهر المعادلات التفاضلية العادية في العديد من سياقات الرياضيات والعلوم الاجتماعية والطبيعية . وتستخدم الأوصاف الرياضية للتغير التفاضلات والمشتقات. وترتبط مختلف التفاضلات والمشتقات والدوال عبر المعادلات، بحيث تكون المعادلة التفاضلية نتيجة تصف الظواهر المتغيرة ديناميكيًا، والتطور، والتغير. غالبًا ما تُعرَّف الكميات على أنها معدل تغير كميات أخرى (على سبيل المثال، مشتقات الإزاحة بالنسبة للزمن)، أو تدرجات الكميات، وهي الطريقة التي تدخل بها في المعادلات التفاضلية. [ 7 ]

تشمل المجالات الرياضية المحددة الهندسة والميكانيكا التحليلية . وتشمل المجالات العلمية جزءًا كبيرًا من الفيزياء وعلم الفلك (الميكانيكا السماوية)، والأرصاد الجوية (نمذجة الطقس)، والكيمياء (معدلات التفاعل)، [ 8 ] وعلم الأحياء (الأمراض المعدية، والتنوع الجيني)، وعلم البيئة ونمذجة السكان (التنافس السكاني)، والاقتصاد (اتجاهات الأسهم، وأسعار الفائدة، وتغيرات أسعار توازن السوق).

لقد درس العديد من علماء الرياضيات المعادلات التفاضلية وساهموا في هذا المجال، بما في ذلك نيوتن ، وليبنيز ، وعائلة برنولي ، وريكاتي ، وكليرو ، ودالمبير ، وأويلر .

ومن الأمثلة البسيطة على ذلك قانون نيوتن الثاني للحركة العلاقة بين الإزاحةx{\displaystyle x}والوقتت{\displaystyle t}جسم تحت تأثير القوةF{\displaystyle F}، ويتم إعطاؤها بواسطة المعادلة التفاضلية

مد2x(ت)دت2=F(x(ت)){\displaystyle m{\frac {\mathrm {d} ^{2}x(t)}{\mathrm {d} t^{2}}}=F(x(t))\,}

مما يقيد حركة جسيم ذي كتلة ثابتةم{\displaystyle m}. على العموم،F{\displaystyle F}هي دالة للموقعx(ت){\displaystyle x(t)}من الجسيم في وقتت{\displaystyle t}الدالة المجهولةx(ت){\displaystyle x(t)}يظهر على جانبي المعادلة التفاضلية، ويُشار إليه في الترميزF(x(ت)){\displaystyle F(x(t))}[ 9 ] [ 10 ] [ 11 ] [ 12 ]

التعريفات

فيما يلي،y{\displaystyle y}هو متغير تابع يمثل دالة غير معروفةy=و(x){\displaystyle y=f(x)}المتغير المستقلx{\displaystyle x}تختلف رموز التفاضل باختلاف المؤلف، وباختلاف الرمز الأكثر فائدة للمهمة المطلوبة. في هذا السياق، يُعدّ رمز لايبنتز...دyدx،د2yدx2،...،دنyدxن{\displaystyle {\frac {dy}{dx}},{\frac {d^{2}y}{dx^{2}}},\ldots ,{\frac {d^{n}y}{dx^{n}}}}يُعدّ أكثر فائدةً للتفاضل والتكامل ، بينما تُعدّ رموز لاغرانجy،y"،...،y(ن){\displaystyle y',y'',\ldots ,y^{(n)}}يُعدّ أكثر فائدة لتمثيل المشتقات من الرتب العليا بشكل مختصر، ورمز نيوتن(y˙،y¨،y...){\displaystyle ({\dot {y}},{\ddot {y}},{\overset {...}{y}})}يُستخدم غالبًا في الفيزياء لتمثيل المشتقات ذات الرتبة المنخفضة بالنسبة للزمن.

التعريف العام

منحF{\displaystyle F}، وهي وظيفة منx{\displaystyle x}،y{\displaystyle y}ومشتقاتy{\displaystyle y}ثم معادلة على الصورة

F(x،y،y،...،y(ن-1))=y(ن){\displaystyle F\left(x,y,y',\ldots ,y^{(n-1)}\right)=y^{(n)}}

تُسمى معادلة تفاضلية عادية صريحة من الرتبةن{\displaystyle n}[ 13 ] [ 14 ]

وبشكل أعم، معادلة تفاضلية عادية ضمنية من الرتبةن{\displaystyle n}يأخذ الشكل التالي: [ 15 ]

F(x،y،y،y"، ...، y(ن))=0{\displaystyle F\left(x,y,y',y'',\ \ldots ,\ y^{(n)}\right)=0}

وهناك تصنيفات أخرى:

مستقل
تكون المعادلة التفاضلية مستقلة إذا لم تعتمد على المتغير x .
خطي
تكون المعادلة التفاضلية خطية إذاF{\displaystyle F}يمكن كتابتها كمزيج خطي من مشتقاتy{\displaystyle y}أي أنه يمكن إعادة كتابتها على النحو التالي:
y(ن)=أنا=0ن-1أأنا(x)y(أنا)+ر(x){\displaystyle y^{(n)}=\sum _{i=0}^{n-1}a_{i}(x)y^{(i)}+r(x)}

أينأأنا(x){\displaystyle a_{i}(x)}ور(x){\displaystyle r(x)}هي دوال متصلة لـx{\displaystyle x}[ 13 ] [ 16 ] [ 17 ]

الوظيفةر(x){\displaystyle r(x)}يُطلق عليه اسم المصطلح المصدر ، مما يؤدي إلى مزيد من التصنيف. [ 16 ] [ 18 ]
متجانس
تكون المعادلة التفاضلية الخطية متجانسة إذار(x)=0{\displaystyle r(x)=0}في هذه الحالة، يوجد دائمًا " الحل البسيط " .y=0{\displaystyle y=0}.
غير متجانس (أو غير متجانس)
تكون المعادلة التفاضلية الخطية غير متجانسة إذار(x)0{\displaystyle r(x)\neq 0}.
غير خطي
معادلة تفاضلية غير خطية.

نظام المعادلات التفاضلية العادية

يشكل عدد من المعادلات التفاضلية المترابطة نظامًا من المعادلات. إذاy{\displaystyle \mathbf {y} }هو متجه عناصره عبارة عن دوال؛y(x)=[y1(x)،y2(x)،...،yم(x)]{\displaystyle \mathbf {y} (x)=[y_{1}(x),y_{2}(x),\ldots ,y_{m}(x)]}، وF{\displaystyle \mathbf {F} }هي دالة ذات قيم متجهة لـy{\displaystyle \mathbf {y} }ومشتقاتها، إذن

y(ن)=F(x،y،y،y"،...،y(ن-1)){\displaystyle \mathbf {y} ^{(n)}=\mathbf {F} \left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^{(n-1)}\right)}

هو نظام صريح من المعادلات التفاضلية العادية من الرتبةن{\displaystyle n}والأبعادم{\displaystyle m}. في شكل متجه عمودي :

(y1(ن)y2(ن)yم(ن))=(و1(x،y،y،y"،...،y(ن-1))و2(x،y،y،y"،...،y(ن-1))وم(x،y،y،y"،...،y(ن-1)))// {y} '',\ldots ,\mathbf {y} ^{(n-1)}\right)\\f_{2}\left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^ {(n-1)}\right)\\\vdots \\f_{m}\left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^{(n-1)}\right)\end{pmatrix}}}

هذه ليست بالضرورة خطية. والنظير الضمني لها هو:

F(x،y،y،y"،...،y(ن))=0{\displaystyle \mathbf {F} \left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^{(n)}\right)={\boldsymbol {0}}}

أين0=(0،0،...،0){\displaystyle {\boldsymbol {0}}=(0,0,\ldots ,0)}هو المتجه الصفري . في شكل مصفوفة

(و1(x،y،y،y"،...،y(ن))و2(x،y،y،y"،...،y(ن))وم(x،y،y،y"،...،y(ن)))=(000){\displaystyle {\begin{pmatrix}f_{1}(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^{(n)})\\f_{2}(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^{(n)})\\\vdots \\f_{m}(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^ {(ن)})\end{pmatrix}}={\begin{pmatrix}0\\0\\\vdots \\0\end{pmatrix}}}

بالنسبة لنظام من الشكلF(x،y،y)=0{\displaystyle \mathbf {F} \left(x,\mathbf {y} ,\mathbf {y} '\right)={\boldsymbol {0}}}تتطلب بعض المصادر أيضًا أن تكون مصفوفة جاكوبيF(x،u،v)v{\displaystyle {\frac {\partial \mathbf {F} (x,\mathbf {u} ,\mathbf {v} )}{\partial \mathbf {v} }}}يُشترط أن تكون المصفوفة غير منفردة لكي يُطلق على هذا النظام اسم نظام معادلات تفاضلية عادية ضمني؛ ويمكن تحويل نظام المعادلات التفاضلية العادية الضمني الذي يحقق شرط عدم التفرد في مصفوفة جاكوبي إلى نظام معادلات تفاضلية عادية صريح. في المصادر نفسها، تُسمى أنظمة المعادلات التفاضلية العادية الضمنية ذات مصفوفة جاكوبي المنفردة بالمعادلات التفاضلية الجبرية . هذا التمييز ليس مجرد تمييز مصطلحي؛ فالمعادلات التفاضلية الجبرية لها خصائص مختلفة جوهريًا، وعادةً ما يكون حلها أكثر تعقيدًا من حل أنظمة المعادلات التفاضلية العادية (غير المنفردة). [ 19 ] [ 20 ] [ 21 ] يُفترض، بالنسبة للمشتقات الإضافية، أن مصفوفة هيسيان وما إلى ذلك تُفترض أيضًا غير منفردة وفقًا لهذا المخطط، مع ملاحظة أنه يمكن (ويُعاد عادةً) إعادة كتابة أي معادلة تفاضلية عادية من رتبة أكبر من واحد كنظام معادلات تفاضلية عادية من الرتبة الأولى ، [ 22 ] مما يجعل معيار تفرد مصفوفة جاكوبي كافيًا لشمولية هذا التصنيف لجميع الرتب.

يمكن تصور سلوك نظام المعادلات التفاضلية العادية من خلال استخدام مخطط الطور .

الحلول

بفرض معادلة تفاضلية

F(x،y،y،...،y(ن))=0{\displaystyle F\left(x,y,y',\ldots ,y^{(n)}\right)=0}

دالةu:أناRR{\displaystyle u:I\subset \mathbb {R} \to \mathbb {R} }، أينأنا{\displaystyle I}إذا كانت فترة، فإنها تسمى حلاً أو منحنى تكاملي لـF{\displaystyle F}، لوu{\displaystyle u}يكونن{\displaystyle n}قابلة للتفاضل مرات علىأنا{\displaystyle I}، و

F(x،u،u، ...، u(ن))=0xأنا.{\displaystyle F(x,u,u',\ \ldots ,\ u^{(n)})=0\quad x\in I.}

بافتراض وجود حلينu:جRR{\displaystyle u:J\subset \mathbb {R} \to \mathbb {R} }وv:أناRR{\displaystyle v:I\subset \mathbb {R} \to \mathbb {R} }،u{\displaystyle u}يُطلق عليه اسم امتداد لـv{\displaystyle v}لوأناج{\displaystyle I\subset J}و

u(x)=v(x)xأنا.{\displaystyle u(x)=v(x)\quad x\in I.\,}

يُطلق على الحل الذي ليس له امتداد اسم الحل الأقصى . الحل المعرف على جميعR{\displaystyle \mathbb {R} }يُطلق عليه اسم الحل العالمي .

حل عام لـن{\displaystyle n}المعادلة من الرتبة n هي حل يحتوي علىن{\displaystyle n}ثوابت تكامل مستقلة اختيارية . يُشتق حل خاص من الحل العام بتحديد قيم معينة للثوابت، تُختار غالبًا لتحقيق مجموعة من " الشروط الابتدائية أو الشروط الحدية ". [ 23 ] الحل الشاذ هو حل لا يمكن الحصول عليه بتعيين قيم محددة للثوابت الاختيارية في الحل العام. [ 24 ]

في سياق المعادلات التفاضلية الخطية، يُمكن أن يُشير مصطلح " الحل الخاص" إلى أي حل للمعادلة التفاضلية (حتى لو لم يُحقق الشروط الابتدائية)، والذي يُضاف إلى الحل المتجانس (الحل العام للمعادلة التفاضلية المتجانسة)، ليُشكّل في النهاية الحل العام للمعادلة التفاضلية الأصلية. هذا هو المصطلح المُستخدم في قسم طريقة التخمين في هذه المقالة، ويُستخدم بكثرة عند مناقشة طريقة المعاملات غير المُحددة وتغيير المعاملات .

حلول ذات مدة زمنية محدودة

بالنسبة للمعادلات التفاضلية العادية غير الخطية المستقلة، من الممكن، في ظل شروط معينة، إيجاد حلول ذات مدة زمنية محدودة [ 25 أي أن النظام، بحكم ديناميكيته، سيصل إلى القيمة صفر عند لحظة زمنية محددة، ويبقى عند الصفر إلى الأبد. لا يمكن أن تكون هذه الحلول ذات المدة الزمنية المحدودة دوالًا تحليلية على كامل خط الأعداد الحقيقية، ولأنها ستكون دوالًا غير ليبشيتزية عند لحظة نهايتها، فإنها لا تدخل ضمن نظرية تفرد حلول المعادلات التفاضلية الليبشيتزية.

على سبيل المثال، المعادلة

y=-علامة(y)|y|،y(0)=1{\displaystyle y'=-{\text{sgn}}(y){\sqrt {|y|}},\,\,y(0)=1}

يقبل الحل ذو المدة المحدودة:

y(x)=14(1-x2+|1-x2|)2{\displaystyle y(x)={\frac {1}{4}}\left(1-{\frac {x}{2}}+\left|1-{\frac {x}{2}}\right|\right)^{2}}

النظريات

حلول فريدة

كانت نظرية الحلول الشاذة للمعادلات التفاضلية العادية والجزئية موضوعًا للبحث منذ عهد لايبنتز، لكنها لم تحظَ باهتمام خاص إلا منذ منتصف القرن التاسع عشر. ومن الأعمال القيّمة، وإن كانت غير معروفة على نطاق واسع، في هذا المجال عمل هوتان (1854). وكان داربو (منذ عام 1873) رائدًا في هذه النظرية، وفي التفسير الهندسي لهذه الحلول، فتح آفاقًا جديدةً أمام العديد من الباحثين، ولا سيما كاسوراتي وكايلي . ويعود الفضل إلى كايلي (1872) في وضع نظرية الحلول الشاذة للمعادلات التفاضلية من الرتبة الأولى، كما تم قبولها في حوالي عام 1900 .

الاختزال إلى التربيعات

كانت المحاولة البدائية في معالجة المعادلات التفاضلية تهدف إلى اختزالها إلى تكاملات تربيعية ، أي التعبير عن الحلول بدلالة دوال معلومة وتكاملاتها. كان هذا ممكنًا للمعادلات الخطية ذات المعاملات الثابتة، ولكن تبين في القرن التاسع عشر أنه مستحيل عمومًا في الحالات الأخرى. ومن ثم، بدأ المحللون بدراسة الدوال التي تمثل حلولًا للمعادلات التفاضلية (لأغراضهم الخاصة)، مما فتح مجالًا جديدًا وخصبًا. وكان كوشي أول من أدرك أهمية هذا المنظور. بعد ذلك، لم يعد السؤال الحقيقي هو ما إذا كان الحل ممكنًا بالتكاملات التربيعية، بل ما إذا كانت معادلة تفاضلية معينة كافية لتعريف دالة، وإذا كان الأمر كذلك، فما هي الخصائص المميزة لهذه الدوال.

نظرية فوكس

ألهمت مذكرتان لفوش [ 26 ] منهجًا جديدًا، طوّره لاحقًا كلٌ من تومي وفروبينيوس . (من؟) كان كوليه مساهمًا بارزًا بدءًا من عام 1869. وقد أبلغ برتراند بطريقته لتكامل نظام غير خطي في عام 1868. هاجم كليبش (1873) النظرية بخطوط موازية لتلك الموجودة في نظريته عن التكاملات الأبيلية . ولأن الأخيرة يمكن تصنيفها وفقًا لخصائص المنحنى الأساسي الذي يظل دون تغيير تحت تحويل كسري، فقد اقترح كليبش تصنيف الدوال المتسامية المعرفة بواسطة المعادلات التفاضلية وفقًا للخصائص الثابتة للأسطح المقابلة.و=0{\displaystyle f=0}في ظل التحويلات النسبية أحادية التقابل.

نظرية لي

منذ عام 1870، أرست أعمال سوفوس لي أسسًا أكثر متانة لنظرية المعادلات التفاضلية. فقد بيّن أن نظريات التكامل التي وضعها علماء الرياضيات الأقدم يمكن ربطها، باستخدام زمر لي ، بمصدر مشترك، وأن المعادلات التفاضلية العادية التي تقبل نفس التحويلات المتناهية الصغر تُظهر صعوبات تكاملية مماثلة. كما ركّز على موضوع تحويلات التماس .

لقد تم إثبات صحة نظرية لي للمجموعات في المعادلات التفاضلية، وتحديداً: (1) أنها توحد العديد من الطرق المخصصة المعروفة لحل المعادلات التفاضلية، و(2) أنها توفر طرقاً جديدة فعالة لإيجاد الحلول. وللنظرية تطبيقات على كل من المعادلات التفاضلية العادية والجزئية. [ 27 ]

يعتمد منهج الحل العام على خاصية التناظر للمعادلات التفاضلية، والتحويلات المتناهية الصغر المستمرة للحلول إلى حلول أخرى ( نظرية لي ). وتُستخدم نظرية الزمر المستمرة ، وجبر لي ، والهندسة التفاضلية لفهم بنية المعادلات التفاضلية الخطية وغير الخطية (الجزئية) لتوليد معادلات قابلة للتكامل، ولإيجاد أزواج لاكس ، ومؤثرات التكرار، وتحويل باكلوند ، وأخيرًا إيجاد حلول تحليلية دقيقة للمعادلات التفاضلية.

تم تطبيق أساليب التناظر على المعادلات التفاضلية التي تنشأ في الرياضيات والفيزياء والهندسة وغيرها من التخصصات.

نظرية ستورم-ليوفيل

نظرية ستورم-ليوفيل هي نظرية لنوع خاص من المعادلات التفاضلية الخطية العادية من الرتبة الثانية. تعتمد حلولها على القيم الذاتية والدوال الذاتية المناظرة لها لمؤثرات خطية مُعرَّفة عبر معادلات خطية متجانسة من الرتبة الثانية . تُعرف هذه المسائل بمسائل ستورم-ليوفيل (SLP)، وهي نسبةً إلى يوهان  سي.  إف.  ستورم وجون ليوفيل اللذين درساها في منتصف القرن التاسع عشر. تتميز مسائل ستورم- ليوفيل  بعدد لا نهائي من القيم الذاتية، وتشكل الدوال الذاتية المناظرة لها مجموعة كاملة ومتعامدة، مما يُتيح إجراء عمليات توسيع متعامدة. تُعد هذه فكرة أساسية في الرياضيات التطبيقية والفيزياء والهندسة. [ 28 ] كما تُستخدم مسائل ستورم-ليوفيل في تحليل بعض المعادلات التفاضلية الجزئية.

وجود الحلول وتفردها

توجد عدة نظريات تثبت وجود حلول فريدة لمسائل القيم الابتدائية التي تتضمن معادلات تفاضلية عادية، سواء على المستوى المحلي أو العالمي. وأهم نظريتين هما:

نظريةافتراضخاتمة
نظرية وجود بيانوF{\displaystyle F}مستمرالوجود المحلي فقط
نظرية بيكارد-ليندلوفF{\displaystyle F}ليبشيتز مستمرالوجود المحلي والتفرد

في شكلها الأساسي، تضمن كلتا النظريتين نتائج محلية فقط، على الرغم من أنه يمكن توسيع الأخيرة لإعطاء نتيجة عالمية، على سبيل المثال، إذا تم استيفاء شروط متباينة غرونوال .

كذلك، فإن نظريات التفرد مثل نظرية ليبشيتز المذكورة أعلاه لا تنطبق على أنظمة المعادلات التفاضلية الجبرية ، والتي قد يكون لها حلول متعددة ناتجة عن الجزء الجبري (غير الخطي) وحده. [ 29 ]

تبسيط نظرية الوجود المحلي والوحدانية

يمكن صياغة النظرية ببساطة على النحو التالي. [ 30 ] بالنسبة للمعادلة ومسألة القيمة الابتدائية: y=F(x،y)،y0=y(x0){\displaystyle y'=F(x,y)\,,\quad y_{0}=y(x_{0})} لوF{\displaystyle F}وF/y{\displaystyle \partial F/\partial y}متصلة في مستطيل مغلق R=[x0-أ،x0+أ]×[y0-ب،y0+ب]{\displaystyle R=[x_{0}-a,x_{0}+a]\times [y_{0}-b,y_{0}+b]} فيx-y{\displaystyle x-y}طائرة، حيثأ{\displaystyle a}وب{\displaystyle b}حقيقية (رمزياً :أ،بR{\displaystyle a,b\in \mathbb {R} }) و ×{\displaystyle \times }يرمز إلى الضرب الديكارتي ، وترمز الأقواس المربعة إلى فترات مغلقة ، إذن توجد فترة أنا=[x0-ح،x0+ح][x0-أ،x0+أ]{\displaystyle I=[x_{0}-h,x_{0}+h]\subset [x_{0}-a,x_{0}+a]} بالنسبة للبعضحR{\displaystyle h\in \mathbb {R} }حيث يمكن إيجاد حل المعادلة المذكورة أعلاه ومسألة القيمة الابتدائية. أي أن هناك حلاً وهو حل وحيد. بما أنه لا يوجد قيد علىF{\displaystyle F}لكي تكون خطية، ينطبق هذا على المعادلات غير الخطية التي تأخذ الشكل التالي:F(x،y){\displaystyle F(x,y)}ويمكن تطبيقه أيضاً على أنظمة المعادلات.

التفرد العالمي وأقصى نطاق للحل

عندما تتحقق فرضيات نظرية بيكارد-ليندلوف، يمكن تعميم الوجود المحلي والوحدانية ليشمل نتيجة شاملة. بتعبير أدق: [ 31 ]

لكل شرط ابتدائي(x0،y0){\displaystyle (x_{0},y_{0})}توجد فترة مفتوحة قصوى فريدة (ربما لا نهائية)

أناالأعلى=(x-،x+)،x±R{±}،x0أناالأعلى{\displaystyle I_{\max }=(x_{-},x_{+}),x_{\pm }\in \mathbb {R} \cup \{\pm \infty \},x_{0}\in I_{\max }}

بحيث يكون أي حل يحقق هذا الشرط الابتدائي بمثابة تقييد للحل الذي يحقق هذا الشرط الابتدائي ذي المجالأناالأعلى{\displaystyle I_{\max }}.

في حالة أنx±±{\displaystyle x_{\pm }\neq \pm \infty }هناك احتمالان فقط

  • انفجار في وقت محدود:ليم سوبxx±y(x){\displaystyle \limsup _{x\to x_{\pm }}\|y(x)\|\to \infty }
  • نطاق التعريف:ليمxx±y(x) Ω¯{\displaystyle \lim _{x\to x_{\pm }}y(x)\ \in \partial {\bar {\Omega }}}

أينΩ{\displaystyle \Omega }هي المجموعة المفتوحة التيF{\displaystyle F}يتم تعريفها، وΩ¯{\displaystyle \partial {\bar {\Omega }}}هي حدودها.

لاحظ أن المجال الأقصى للحل

  • هي دائماً فترة (لتحقيق التفرد)
  • قد يكون أصغر منR{\displaystyle \mathbb {R} }
  • قد يعتمد ذلك على الاختيار المحدد لـ(x0،y0){\displaystyle (x_{0},y_{0})}.
مثال.
y=y2{\displaystyle y'=y^{2}}

هذا يعني أنF(x،y)=y2{\displaystyle F(x,y)=y^{2}}، وهوج1{\displaystyle C^{1}}وبالتالي فهي متصلة محليًا وفقًا لشرط ليبشيتز، مما يحقق نظرية بيكارد-ليندلوف.

حتى في مثل هذه الحالة البسيطة، لا يمكن أن يكون نطاق الحل الأقصى هو كل شيءR{\displaystyle \mathbb {R} }لأن الحل هو

y(x)=y0(x0-x)y0+1{\displaystyle y(x)={\frac {y_{0}}{(x_{0}-x)y_{0}+1}}}

الذي يمتلك أقصى مدى:

{Ry0=0(-،x0+1y0)y0>0(x0+1y0،+)y0<0{\displaystyle {\begin{cases}\mathbb {R} &y_{0}=0\\[4pt]\left(-\infty ,x_{0}+{\frac {1}{y_{0}}}\right)&y_{0}>0\\[4pt]\left(x_{0}+{\frac {1}{y_{0}}},+\infty \right)&y_{0}<0\end{cases}}}

يُظهر هذا بوضوح أن الحد الأقصى للفترة الزمنية قد يعتمد على الشروط الابتدائية. مجالy{\displaystyle y}يمكن اعتبار ذلك بمثابةR(x0+1/y0)،{\displaystyle \mathbb {R} \setminus (x_{0}+1/y_{0}),}لكن هذا سيؤدي إلى مجال ليس فترة، بحيث يكون الجانب المقابل للشرط الأولي منفصلاً عن الشرط الأولي، وبالتالي لا يتم تحديده بشكل فريد من خلاله.

المجال الأقصى ليسR{\displaystyle \mathbb {R} }لأن

ليمxx±y(x)،{\displaystyle \lim _{x\to x_{\pm }}\|y(x)\|\to \infty ,}

وهو أحد الحالتين المحتملتين وفقًا للنظرية المذكورة أعلاه.

تخفيض الطلب

عادة ما يكون حل المعادلات التفاضلية أسهل إذا أمكن تقليل رتبة المعادلة.

الاختزال إلى نظام من الدرجة الأولى

أي معادلة تفاضلية صريحة من الرتبةن{\displaystyle n}،

F(x،y،y،y"، ...، y(ن-1))=y(ن){\displaystyle F\left(x,y,y',y'',\ \ldots ,\ y^{(n-1)}\right)=y^{(n)}}

يمكن كتابتها كنظام منن{\displaystyle n}المعادلات التفاضلية من الدرجة الأولى عن طريق تعريف عائلة جديدة من الدوال المجهولة

yأنا=y(أنا-1).{\displaystyle y_{i}=y^{(i-1)}.\!}

لأنا=1،2،...،ن{\displaystyle i=1,2,\ldots ,n}. الن{\displaystyle n}ثم يكون النظام ذو الأبعاد n من المعادلات التفاضلية المقترنة من الدرجة الأولى هو

y1=y2y2=y3yن-1=yنyن=F(x،y1،...،yن).{\displaystyle {\begin{array}{rcl}y_{1}'&=&y_{2}\\y_{2}'&=&y_{3}\\&\vdots &\\y_{n-1}'&=&y_{n}\\y_{n}'&=&F(x,y_{1},\ldots ,y_{n}).\end{array}}}

بشكل أكثر إيجازًا في تدوين المتجهات:

y=F(x،y){\displaystyle \mathbf {y} '=\mathbf {F} (x,\mathbf {y} )}

أين

y=(y1،...،yن)،F(x،y1،...،yن)=(y2،...،yن،F(x،y1،...،yن)).{\displaystyle \mathbf {y} =(y_{1},\ldots ,y_{n}),\quad \mathbf {F} (x,y_{1},\ldots ,y_{n})=(y_{2},\ldots ,y_{n},F(x,y_{1},\ldots ,y_{n})).}

ملخص الحلول الدقيقة

بعض المعادلات التفاضلية لها حلول يمكن كتابتها بصيغة دقيقة ومغلقة. وقد تم تقديم عدة فئات مهمة منها هنا.

في الجدول أدناه،P(x){\displaystyle P(x)}،سؤال(x){\displaystyle Q(x)}،P(y){\displaystyle P(y)}،سؤال(y){\displaystyle Q(y)}، وم(x،y){\displaystyle M(x,y)}،شمال(x،y){\displaystyle N(x,y)}هل توجد أي وظائف قابلة للتكامل منx{\displaystyle x}،y{\displaystyle y}؛ب{\displaystyle b}وج{\displaystyle c}هي ثوابت حقيقية معطاة؛ج1،ج2،...{\displaystyle C_{1},C_{2},\ldots }هي ثوابت اختيارية ( مركبة بشكل عام). المعادلات التفاضلية في صيغها المكافئة والبديلة التي تؤدي إلى الحل من خلال التكامل.

في الحلول التكاملية،λ{\displaystyle \lambda }وε{\displaystyle \varepsilon }هي متغيرات وهمية للتكامل (النظائر المتصلة للمؤشرات في الجمع )، والترميزxF(λ)دλ{\displaystyle \int ^{x}F(\lambda )\,d\lambda }يعني ببساطة الاندماجF(λ){\displaystyle F(\lambda )}بالنسبة إلىλ{\displaystyle \lambda }ثم بعد التكامل، يتم استبدالλ=x{\displaystyle \lambda =x}، دون إضافة ثوابت (مذكورة صراحةً).

المعادلات القابلة للفصل

معادلة تفاضليةطريقة الحلحل عام
من الدرجة الأولى، قابلة للفصل فيx{\displaystyle x}وy{\displaystyle y}(الحالة العامة، انظر أدناه للحالات الخاصة) [ 32 ]

P1(x)سؤال1(y)+P2(x)سؤال2(y)دyدx=0P1(x)سؤال1(y)دx+P2(x)سؤال2(y)دy=0{\displaystyle {\begin{aligned}P_{1}(x)Q_{1}(y)+P_{2}(x)Q_{2}(y)\,{\frac {dy}{dx}}&=0\\P_{1}(x)Q_{1}(y)\,dx+P_{2}(x)Q_{2}(y)\,dy&=0\end{aligned}}}

فصل المتغيرات (القسمة علىP2سؤال1{\displaystyle P_{2}Q_{1}}).xP1(λ)P2(λ)دλ+yسؤال2(λ)سؤال1(λ)دλ=ج{\displaystyle \int ^{x}{\frac {P_{1}(\lambda )}{P_{2}(\lambda )}}\,d\lambda +\int ^{y}{\frac {Q_{2}(\lambda )}{Q_{1}(\lambda )}}\,d\lambda =C}
من الدرجة الأولى، قابلة للفصل فيx{\displaystyle x}[ 30 ]

دyدx=F(x)دy=F(x)دx{\displaystyle {\begin{aligned}{\frac {dy}{dx}}&=F(x)\\dy&=F(x)\,dx\end{aligned}}}

التكامل المباشر.y=xF(λ)دλ+ج{\displaystyle y=\int ^{x}F(\lambda )\,d\lambda +C}
من الدرجة الأولى، مستقلة، قابلة للفصل فيy{\displaystyle y}[ 30 ]

دyدx=F(y)دy=F(y)دx{\displaystyle {\begin{aligned}{\frac {dy}{dx}}&=F(y)\\dy&=F(y)\,dx\end{aligned}}}

فصل المتغيرات (القسمة علىF{\displaystyle F}).x=yدλF(λ)+ج{\displaystyle x=\int ^{y}{\frac {d\lambda }{F(\lambda )}}+C}
من الدرجة الأولى، قابلة للفصل فيx{\displaystyle x}وy{\displaystyle y}[ 30 ]

P(y)دyدx+سؤال(x)=0P(y)دy+سؤال(x)دx=0{\displaystyle {\begin{aligned}P(y){\frac {dy}{dx}}+Q(x)&=0\\P(y)\,dy+Q(x)\,dx&=0\end{aligned}}}

دمجها في جميع أنحاء المكان.yP(λ)دλ+xسؤال(λ)دλ=ج{\displaystyle \int ^{y}P(\lambda )\,d\lambda +\int ^{x}Q(\lambda )\,d\lambda =C}

المعادلات العامة من الدرجة الأولى

معادلة تفاضليةطريقة الحلحل عام
من الدرجة الأولى، متجانس [ 30 ]

دyدx=F(yx){\displaystyle {\frac {dy}{dx}}=F\left({\frac {y}{x}}\right)}

ضع y = ux ، ثم قم بالحل عن طريق فصل المتغيرات في u و x .ln(جx)=y/xدλF(λ)-λ{\displaystyle \ln(Cx)=\int ^{y/x}{\frac {d\lambda }{F(\lambda )-\lambda }}}
من الدرجة الأولى، قابلة للفصل [ 32 ]

yم(xy)+xشمال(xy)دyدx=0yم(xy)دx+xشمال(xy)دy=0{\displaystyle {\begin{aligned}yM(xy)+xN(xy)\,{\frac {dy}{dx}}&=0\\yM(xy)\,dx+xN(xy)\,dy&=0\end{aligned}}}

فصل المتغيرات (القسمة علىxy{\displaystyle xy}).

ln(جx)=xyشمال(λ)دλλ[شمال(λ)-م(λ)]{\displaystyle \ln(Cx)=\int ^{xy}{\frac {N(\lambda )\,d\lambda }{\lambda [N(\lambda )-M(\lambda )]}}}

لوشمال=م{\displaystyle N=M}الحل هوxy=ج{\displaystyle xy=C}.

التفاضل التام ، من الدرجة الأولى [ 30 ]

م(x،y)دyدx+شمال(x،y)=0م(x،y)دy+شمال(x،y)دx=0{\displaystyle {\begin{aligned}M(x,y){\frac {dy}{dx}}+N(x,y)&=0\\M(x,y)\,dy+N(x,y)\,dx&=0\end{aligned}}}

أينمy=شمالx{\displaystyle {\frac {\partial M}{\partial y}}={\frac {\partial N}{\partial x}}}

دمجها في جميع أنحاء المكان.F(x،y)=xم(λ،y)دλ+yY(λ)دλ=yشمال(x،λ)دλ+xX(λ)دλ=ج{\displaystyle {\begin{aligned}F(x,y)&=\int ^{x}M(\lambda ,y)\,d\lambda +\int ^{y}Y(\lambda )\,d\lambda \\&=\int ^{y}N(x,\lambda )\,d\lambda +\int ^{x}X(\lambda )\,d\lambda =C\end{aligned}}}

أين Y(y)=شمال(x،y)-yxم(λ،y)دλ{\displaystyle Y(y)=N(x,y)-{\frac {\partial }{\partial y}}\int ^{x}M(\lambda ,y)\,d\lambda }و X(x)=م(x،y)-xyشمال(x،λ)دλ{\displaystyle X(x)=M(x,y)-{\frac {\partial }{\partial x}}\int ^{y}N(x,\lambda )\,d\lambda }

التفاضل غير الدقيق ، من الدرجة الأولى [ 30 ]

م(x،y)دyدx+شمال(x،y)=0م(x،y)دy+شمال(x،y)دx=0{\displaystyle {\begin{aligned}M(x,y){\frac {dy}{dx}}+N(x,y)&=0\\M(x,y)\,dy+N(x,y)\,dx&=0\end{aligned}}}

أينمyشمالx{\displaystyle {\frac {\partial M}{\partial y}}\neq {\frac {\partial N}{\partial x}}}

عامل التكاملμ(x،y){\displaystyle \mu (x,y)}مُرضٍ

(μم)y=(μشمال)x{\displaystyle {\frac {\partial (\mu M)}{\partial y}}={\frac {\partial (\mu N)}{\partial x}}}

لوμ(x،y){\displaystyle \mu (x,y)}يمكن إيجادها بطريقة مناسبة، ثم

F(x،y)=xμ(λ،y)م(λ،y)دλ+yY(λ)دλ=yμ(x،λ)شمال(x،λ)دλ+xX(λ)دλ=ج{\displaystyle {\begin{aligned}F(x,y)=&\int ^{x}\mu (\lambda ,y)M(\lambda ,y)\,d\lambda +\int ^{y}Y(\lambda )\,d\lambda \\=&\int ^{y}\mu (x,\lambda )N(x,\lambda )\,d\lambda +\int ^{x}X(\lambda )\,d\lambda =C\end{aligned}}}

أينY(y)=شمال(x،y)-yxμ(λ،y)م(λ،y)دλ{\displaystyle Y(y)=N(x,y)-{\frac {\partial }{\partial y}}\int ^{x}\mu (\lambda ,y)M(\lambda ,y)\,d\lambda } وX(x)=م(x،y)-xyμ(x،λ)شمال(x،λ)دλ{\displaystyle X(x)=M(x,y)-{\frac {\partial }{\partial x}}\int ^{y}\mu (x,\lambda )N(x,\lambda )\,d\lambda }

المعادلات العامة من الدرجة الثانية

معادلة تفاضليةطريقة الحلحل عام
من الدرجة الثانية، مستقلة [ 33 ]

د2yدx2=F(y){\displaystyle {\frac {d^{2}y}{dx^{2}}}=F(y)}

اضرب طرفي المعادلة في2دyدx{\displaystyle 2{\frac {dy}{dx}}}، بديل2دyدxد2yدx2=ددx(دyدx)2=2دyدxF(y){\displaystyle 2{\frac {dy}{dx}}{\frac {d^{2}y}{dx^{2}}}={\frac {d}{dx}}\left({\frac {dy}{dx}}\right)^{2}=2{\frac {dy}{dx}}F(y)}ثم قم بإجراء التكامل مرتين.x=±yدλ2λF(ε)دε+ج1+ج2{\displaystyle x=\pm \int ^{y}{\frac {d\lambda }{\sqrt {2\int ^{\lambda }F(\varepsilon )\,d\varepsilon +C_{1}}}}+C_{2}}

المعادلات الخطية حتى الرتبة n

معادلة تفاضليةطريقة الحلحل عام
معاملات الدالة من الدرجة الأولى، الخطية، غير المتجانسة [ 30 ]

دyدx+P(x)y=سؤال(x){\displaystyle {\frac {dy}{dx}}+P(x)y=Q(x)}

عامل التكامل:هـxP(λ)دλ.{\displaystyle e^{\int ^{x}P(\lambda )\,d\lambda }.}y=هـ-xP(λ)دλ[xهـλP(ε)دεسؤال(λ)دλ+ج]{\displaystyle y=e^{-\int ^{x}P(\lambda )\,d\lambda }\left[\int ^{x}e^{\int ^{\lambda }P(\varepsilon )\,d\varepsilon }Q(\lambda )\,d\lambda +C\right]}
معاملات دالة خطية غير متجانسة من الدرجة الثانية

د2yدx2+2ص(x)دyدx+(ص(x)2+ص(x))y=q(x){\displaystyle {\frac {d^{2}y}{dx^{2}}}+2p(x){\frac {dy}{dx}}+\left(p(x)^{2}+p'(x)\right)y=q(x)}

عامل التكامل:هـxP(λ)دλ{\displaystyle e^{\int ^{x}P(\lambda )\,d\lambda }}y=هـ-xP(λ)دλ[x(ξهـλP(ε)دεسؤال(λ)دλ)دξ+ج1x+ج2]{\displaystyle y=e^{-\int ^{x}P(\lambda )\,d\lambda }\left[\int ^{x}\left(\int ^{\xi }e^{\int ^{\lambda }P(\varepsilon )\,d\varepsilon }Q(\lambda )\,d\lambda \right)d\xi +C_{1}x+C_{2}\right]}
معاملات ثابتة من الدرجة الثانية، خطية، غير متجانسة [ 34 ]

د2yدx2+بدyدx+جy=ر(x){\displaystyle {\frac {d^{2}y}{dx^{2}}}+b{\frac {dy}{dx}}+cy=r(x)}

الوظيفة التكميليةyج{\displaystyle y_{c}}: يفترضyج(x)=هـαx{\displaystyle y_{c}(x)=e^{\alpha x}}، استبدل وحل متعددة الحدود فيα{\displaystyle \alpha }لإيجاد الدوال المستقلة خطيًاهـαجx{\displaystyle e^{\alpha _{j}x}}.

التكامل الخاصyص{\displaystyle y_{p}}بشكل عام، طريقة تغيير المعاملات ، وإن كانت بسيطة للغاية.ر(x){\displaystyle r(x)}قد يكون التفتيش مفيدًا. [ 30 ]

y=yج+yص{\displaystyle y=y_{c}+y_{p}}

لوب2>4ج{\displaystyle b^{2}>4c}، ثم

yج=ج1هـ-x2(ب+ب2-4ج)+ج2هـ-x2(ب-ب2-4ج){\displaystyle y_{c}=C_{1}e^{-{\frac {x}{2}}\,\left(b+{\sqrt {b^{2}-4c}}\right)}+C_{2}e^{-{\frac {x}{2}}\,\left(b-{\sqrt {b^{2}-4c}}\right)}}

لوب2=4ج{\displaystyle b^{2}=4c}، ثم

yج=(ج1x+ج2)هـ-بx2{\displaystyle y_{c}=(C_{1}x+C_{2})e^{-{\frac {bx}{2}}}}

لوب2<4ج{\displaystyle b^{2}<4c}، ثم

yج=هـ-بx2[ج1الخطيئة(x4ج-ب22)+ج2كوس(x4ج-ب22)]{\displaystyle y_{c}=e^{-{\frac {bx}{2}}}\left[C_{1}\sin \left(x\,{\frac {\sqrt {4c-b^{2}}}{2}}\right)+C_{2}\cos \left(x\,{\frac {\sqrt {4c-b^{2}}}{2}}\right)\right]}

ن{\displaystyle n}معاملات ثابتة من الرتبة n، خطية، غير متجانسة [ 34 ]

ج=0نبجدجyدxج=ر(x){\displaystyle \sum _{j=0}^{n}b_{j}{\frac {d^{j}y}{dx^{j}}}=r(x)}

الوظيفة التكميليةyج{\displaystyle y_{c}}: يفترضyج(x)=هـαx{\displaystyle y_{c}(x)=e^{\alpha x}}، استبدل وحل متعددة الحدود فيα{\displaystyle \alpha }لإيجاد الدوال المستقلة خطيًاهـαجx{\displaystyle e^{\alpha _{j}x}}.

التكامل الخاصyص{\displaystyle y_{p}}بشكل عام، طريقة تغيير المعاملات ، وإن كانت بسيطة للغاية.ر(x){\displaystyle r(x)}قد يكون التفتيش مفيدًا. [ 30 ]

y=yج+yص{\displaystyle y=y_{c}+y_{p}}

منذαج{\displaystyle \alpha _{j}}هي حلول متعددة الحدود من الدرجةن{\displaystyle n}:ج=1ن(α-αج)=0{\textstyle \prod _{j=1}^{n}(\alpha -\alpha _{j})=0}ثم: لـαج{\displaystyle \alpha _{j}}جميعها مختلفة، yج=ج=1نججهـαجx{\displaystyle y_{c}=\sum _{j=1}^{n}C_{j}e^{\alpha _{j}x}} لكل جذرαج{\displaystyle \alpha _{j}}كرركج{\displaystyle k_{j}}مرات، yج=ج=1ن(=1كججج،x-1)هـαجx{\displaystyle y_{c}=\sum _{j=1}^{n}\left(\sum _{\ell =1}^{k_{j}}C_{j,\ell }x^{\ell -1}\right)e^{\alpha _{j}x}} بالنسبة للبعضαج{\displaystyle \alpha _{j}}معقد، ثم وضعαج=χج+أناγج{\displaystyle \alpha _{j}=\chi _{j}+i\gamma _{j}}وباستخدام صيغة أويلر ، يمكن كتابة بعض الحدود في النتائج السابقة بالشكل التالي: ججهـαجx=ججهـχجxكوس(γجx+φج){\displaystyle C_{j}e^{\alpha _{j}x}=C_{j}e^{\chi _{j}x}\cos(\gamma _{j}x+\varphi _{j})} أينφج{\displaystyle \varphi _{j}}هو ثابت اعتباطي (إزاحة الطور).

تخمين الحلول

عندما تفشل جميع الطرق الأخرى لحل معادلة تفاضلية عادية، أو في الحالات التي يكون لدينا فيها حدسٌ ما حول شكل حل المعادلة التفاضلية، يُمكن أحيانًا حل المعادلة التفاضلية ببساطة عن طريق تخمين الحل والتحقق من صحته. [ 35 ] لاستخدام هذه الطريقة، نخمن ببساطة حلاً للمعادلة التفاضلية، ثم نعوض بهذا الحل في المعادلة التفاضلية للتحقق مما إذا كان يحقق المعادلة. إذا كان كذلك، فلدينا حلٌ خاص للمعادلة التفاضلية، وإلا نبدأ من جديد ونجرب تخمينًا آخر. على سبيل المثال، يُمكننا تخمين أن حل المعادلة التفاضلية له الشكل التالي:y=أهـأناαت{\displaystyle y=Ae^{i\alpha t}}لأن هذا حل شائع جداً يتصرف فيزيائياً بطريقة جيبية.

في حالة معادلة تفاضلية عادية من الرتبة الأولى غير متجانسة، نحتاج أولاً إلى إيجاد حل للجزء المتجانس من المعادلة التفاضلية، والمعروف أيضاً بالمعادلة المتجانسة المرتبطة بها، ثم إيجاد حل للمعادلة غير المتجانسة بأكملها عن طريق التخمين. وأخيراً، نجمع هذين الحلين معاً للحصول على الحل العام للمعادلة التفاضلية العادية، وهو:

الحل العام=الحل العام للمعادلة المتجانسة المرتبطة+حل خاص{\displaystyle {\text{general solution}}={\text{general solution of the associated homogeneous equation}}+{\text{particular solution}}}

برنامج لحل المعادلات التفاضلية العادية

  • ماكسيما ، نظام جبر حاسوبي مفتوح المصدر .
  • COPASI ، حزمة برامج مجانية ( رخصة فنية 2.0 ) لتكامل وتحليل المعادلات التفاضلية العادية.
  • MATLAB ، تطبيق حاسوبي تقني (مختبر المصفوفات)
  • GNU Octave ، لغة برمجة عالية المستوى، مخصصة في المقام الأول للحسابات العددية.
  • Scilab ، تطبيق مفتوح المصدر للحسابات العددية.
  • Maple ، تطبيق خاص لإجراء العمليات الحسابية الرمزية.
  • Mathematica ، تطبيق احتكاري مخصص في المقام الأول للحسابات الرمزية.
  • SymPy ، حزمة بايثون يمكنها حل المعادلات التفاضلية العادية بشكل رمزي
  • جوليا (لغة برمجة) ، وهي لغة عالية المستوى مخصصة في المقام الأول للحسابات العددية.
  • SageMath ، وهو تطبيق مفتوح المصدر يستخدم بنية تشبه لغة بايثون مع مجموعة واسعة من القدرات التي تغطي عدة فروع من الرياضيات.
  • SciPy ، حزمة بايثون تتضمن وحدة تكامل المعادلات التفاضلية العادية.
  • Chebfun ، حزمة مفتوحة المصدر، مكتوبة بلغة MATLAB ، للحساب باستخدام الدوال بدقة تصل إلى 15 رقمًا.
  • GNU R ، بيئة حسابية مفتوحة المصدر مخصصة في المقام الأول للإحصاء، والتي تتضمن حزمًا لحل المعادلات التفاضلية العادية.

انظر أيضاً

ملحوظات

  1. دينيس ج. زيل (15 مارس 2012). مدخل إلى المعادلات التفاضلية مع تطبيقات النمذجة . سينجايج ليرنينج. ISBN 978-1-285-40110-2أُرشف من المصدر الأصلي بتاريخ 17 يناير 2020. تم الاطلاع عليه بتاريخ 11 يوليو 2019 .
  2. ما هو أصل مصطلح "المعادلات التفاضلية العادية"؟ . hsm.stackexchange.com . Stack Exchange . تم الاطلاع عليه بتاريخ 28-07-2016 .
  3. كاراس، تيرو؛ أيتالا، ميكا؛ أيلا، تيمو؛ لاين، سامولي (2022). "توضيح فضاء تصميم النماذج التوليدية القائمة على الانتشار". arXiv : 2206.00364 [ cs.CV ].
  4. بوتشر، ج. س. (15 ديسمبر 2000). "الأساليب العددية للمعادلات التفاضلية العادية في القرن العشرين" . مجلة الرياضيات الحسابية والتطبيقية . التحليل العددي 2000. المجلد السادس: المعادلات التفاضلية العادية والمعادلات التكاملية. 125 (1): 1-29 . رمز Bibcode : 2000JCoAM.125....1B . doi : 10.1016/S0377-0427(00)00455-6 . ISSN 0377-0427 . 
  5. غرينبيرغ، مايكل د. (2012). المعادلات التفاضلية العادية . هوبوكين، نيوجيرسي: وايلي. ISBN 978-1-118-23002-2.
  6. أكتون، فورمان س. (1990). الطرق العددية الفعالة . سبكتروم. واشنطن العاصمة: الجمعية الرياضية الأمريكية. ISBN 978-1-4704-5727-3.
  7. دينيس، بياكاتوندا (2020-12-10). "نظرة عامة على الطرق العددية والتحليلية لحل المعادلات التفاضلية العادية". arXiv : 2012.07558 [ math.HO ].
  8. الرياضيات للكيميائيين، دي إم هيرست، دار ماكميلان للنشر ، 1976، (بدون رقم ISBN) رقم SBN: 333-18172-7
  9. كريزيج (1972 ، ص 64) 
  10. سيمونز (1972 ، ص 1، 2) 
  11. هاليداي وريزنيك (1977 ، ص 78) 
  12. ^ تيبلر (1991 ، ص 78-83) 
  13. 1 2 هاربر (1976 ، ص 127) 
  14. كريزيج (1972 ، ص 2) 
  15. سيمونز (1972 ، ص 3) 
  16. 1 2 كريزيج (1972 ، ص 24)
  17. سيمونز (1972 ، ص 47)
  18. هاربر (1976 ، ص 128)
  19. كريزيج (1972 ، ص 12) 
  20. ^ آشر وبيتزولد (1998 ، ص 12) 
  21. آخيم إيلشمان؛ تيمو رايس (2014). دراسات في المعادلات التفاضلية الجبرية II . سبرينغر. ص 104-105 . ISBN  978-3-319-11050-9.
  22. ^ آشر وبيتزولد (1998 ، ص 5) 
  23. كريزيج (1972 ، ص 78) 
  24. كريزيج (1972 ، ص 4) 
  25. فارديا ت. هايمو (1985). "المعادلات التفاضلية ذات الزمن المحدود". المؤتمر الرابع والعشرون لمعهد مهندسي الكهرباء والإلكترونيات حول التحكم واتخاذ القرارات، 1985. الصفحات 1729-1733. doi : 10.1109 /CDC.1985.268832 . S2CID 45426376 .  
  26. كريل ، 1866، 1868
  27. دريسنر (1999 ، ص 9) 
  28. لوجان، ج. (2013). الرياضيات التطبيقية (الطبعة الرابعة).
  29. ^ آشر وبيتزولد (1998 ، ص 13) 
  30. 1 2 3 4 5 6 7 8 9 10 المعادلات التفاضلية الأولية ومسائل القيم الحدية (الطبعة الرابعة)، دبليو إي بويس، آر سي ديبريما، وايلي إنترناشونال، جون وايلي وأولاده، 1986، رقم ISBN 0-471-83824-1
  31. بوسكاين؛ شيتور 2011، ص 21
  32. 1 2 دليل الرياضيات للصيغ والجداول (الطبعة الثالثة)، إس. ليبشوتز، إم آر شبيغل، جيه. ليو، سلسلة شوم أوتلاين، 2009، ISC_2N 978-0-07-154855-7
  33. التحليل الأولي الإضافي، آر. بورتر، جي. بيل وأولاده (لندن)، 1978، رقم ISBN 0-7135-1594-5
  34. 1 2 الأساليب الرياضية للفيزياء والهندسة، كيه إف رايلي، إم بي هوبسون، إس جيه بينس، مطبعة جامعة كامبريدج، 2010، ISC_2N 978-0-521-86153-3
  35. فيلدر، غاري؛ فيلدر، كيني (2016). الأساليب الرياضية في الهندسة والفيزياء: مواضيع تمهيدية . هوبوكين، نيوجيرسي: وايلي. ص 39. ISBN  978-1-118-44960-8.

مراجع

فهرس