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

في الرياضيات ، المعادلة التفاضلية الخطية هي معادلة تفاضلية خطية في الدالة المجهولة ومشتقاتها، لذا يمكن كتابتها على الصورة التالية :أ0(x)y+أ1(x)y+أ2(x)y"+أن(x)y(ن)=ب(x){\displaystyle a_{0}(x)y+a_{1}(x)y'+a_{2}(x)y''\cdots +a_{n}(x)y^{(n)}=b(x)} حيث أن a 0 ( x ) , ..., a n ( x ) و b ( x ) هي دوال قابلة للتفاضل بشكل عشوائي لا تحتاج إلى أن تكون خطية، و y ′, ..., y ( n ) هي المشتقات المتتالية لدالة غير معروفة y للمتغير x .

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

أنواع الحلول

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

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

المصطلحات الأساسية

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

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

المؤثر التفاضلي الخطي

المؤثر التفاضلي الأساسي من الرتبة i هو دالة تربط أي دالة قابلة للتفاضل بمشتقتها من الرتبة i ، أو في حالة وجود عدة متغيرات، بإحدى مشتقاتها الجزئية من الرتبة i . ويُشار إليه عادةً بـ دأنادxأنا{\displaystyle {\frac {d^{i}}{dx^{i}}}} في حالة الدوال أحادية المتغير ، و أنا1++أنانx1أنا1xنأنان{\displaystyle {\frac {\partial ^{i_{1}+\cdots +i_{n}}}{\partial x_{1}^{i_{1}}\cdots \partial x_{n}^{i_{n}}}}} في حالة الدوال ذات n متغير. تشمل المؤثرات التفاضلية الأساسية المشتقة من الرتبة 0، وهي دالة التطابق.

المؤثر التفاضلي الخطي (يُختصر في هذه المقالة بالمؤثر الخطي أو ببساطة المؤثر ) هو توليفة خطية من المؤثرات التفاضلية الأساسية، حيث تكون الدوال القابلة للتفاضل هي المعاملات. في حالة المتغير الواحد، يكون للمؤثر الخطي الشكل [ 1 ].أ0(x)+أ1(x)ددx++أن(x)دندxن،{\displaystyle a_{0}(x)+a_{1}(x){\frac {d}{dx}}+\cdots +a_{n}(x){\frac {d^{n}}{dx^{n}}},} حيث أن a 0 ( x ), ..., a n ( x ) هي دوال قابلة للتفاضل، والعدد الصحيح غير السالب n هو رتبة المؤثر (إذا لم تكن a n ( x ) هي الدالة الصفرية ).

ليكن L مؤثرًا تفاضليًا خطيًا. يُرمز عادةً لتطبيق L على دالة f بـ Lf أو Lf ( X ) عند الحاجة إلى تحديد المتغير (يجب عدم الخلط بين هذا وبين الضرب). المؤثر التفاضلي الخطي هو مؤثر خطي ، لأنه يحول المجاميع إلى مجاميع، والضرب في عدد قياسي إلى الضرب في نفس العدد القياسي.

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

تتيح لغة المؤثرات كتابة مختصرة للمعادلات التفاضلية: إذا ل=أ0(x)+أ1(x)ددx++أن(x)دندxن،{\displaystyle L=a_{0}(x)+a_{1}(x){\frac {d}{dx}}+\cdots +a_{n}(x){\frac {d^{n}}{dx^{n}}},} إذا كان مؤثرًا تفاضليًا خطيًا، فإن المعادلة أ0(x)y+أ1(x)y+أ2(x)y"++أن(x)y(ن)=ب(x){\displaystyle a_{0}(x)y+a_{1}(x)y'+a_{2}(x)y''+\cdots +a_{n}(x)y^{(n)}=b(x)} قد يُعاد كتابتها لy=ب(x).{\displaystyle Ly=b(x).}

قد يكون هناك عدة صيغ مختلفة لهذه الصيغة؛ على وجه الخصوص، قد يظهر متغير التفاضل بشكل صريح أو لا يظهر في y والطرف الأيمن من المعادلة، مثل Ly ( x ) = b ( x ) أو Ly = b .

نواة المؤثر التفاضلي الخطي هي نواته كتحويل خطي، أي الفضاء المتجهي لحلول المعادلة التفاضلية (المتجانسة) Ly = 0 .

في حالة المؤثر التفاضلي العادي من الرتبة n ، تنص نظرية وجود كاراثيودوري على أنه في ظل شروط بسيطة للغاية، تكون نواة L فضاءً متجهيًا ذا بُعد n ، وأن حلول المعادلة Ly ( x ) = b ( x ) لها الشكل التالي: S0(x)+ج1S1(x)++جنSن(x)،{\displaystyle S_{0}(x)+c_{1}S_{1}(x)+\cdots +c_{n}S_{n}(x),} حيث c₁ , ..., cₙ أعداد اختيارية. عادةً، تتحقق فرضيات نظرية كاراثيودوري في فترة I ، إذا كانت الدوال b , a₀ , ... , aₙ متصلة في I ، وكان هناك عدد حقيقي موجب k بحيث يكون | aₙ ( x ) | > k لكل x في I.

معادلة متجانسة ذات معاملات ثابتة

تكون المعادلة التفاضلية الخطية المتجانسة ذات معاملات ثابتة إذا كانت على الصورة التالية: أ0y+أ1y+أ2y"++أنy(ن)=0{\displaystyle a_{0}y+a_{1}y'+a_{2}y''+\cdots +a_{n}y^{(n)}=0} حيث تمثل a1 ، ...، an أعدادًا (حقيقية أو مركبة). بعبارة أخرى، يكون لها معاملات ثابتة إذا تم تعريفها بواسطة مؤثر خطي ذي معاملات ثابتة.

يعود تاريخ دراسة هذه المعادلات التفاضلية ذات المعاملات الثابتة إلى ليونارد أويلر ، الذي قدم الدالة الأسية.هـx{\displaystyle e^{x}}وهو الحل الوحيد للمعادلةو=و{\displaystyle f'=f}بحيثو(0)=1{\displaystyle f(0)=1}ويترتب على ذلك أن المشتقة النونية لـهـجx{\displaystyle e^{cx}}يكونجنهـجx{\displaystyle c^{n}e^{cx}}وهذا يسمح بحل المعادلات التفاضلية الخطية المتجانسة بسهولة تامة.

يترك أ0y+أ1y+أ2y"++أنy(ن)=0{\displaystyle a_{0}y+a_{1}y'+a_{2}y''+\cdots +a_{n}y^{(n)}=0} لتكن معادلة تفاضلية خطية متجانسة ذات معاملات ثابتة (أي أن a 0 ، ...، a n هي أعداد حقيقية أو مركبة).

إن البحث عن حلول لهذه المعادلة التي تأخذ الشكل e αx يكافئ البحث عن الثوابت α بحيث أ0هـαx+أ1αهـαx+أ2α2هـαx++أنαنهـαx=0.{\displaystyle a_{0}e^{\alpha x}+a_{1}\alpha e^{\alpha x}+a_{2}\alpha ^{2}e^{\alpha x}+\cdots +a_{n}\alpha ^{n}e^{\alpha x}=0.} يُظهر تحليل العامل المشترك e αx (الذي لا يساوي الصفر أبدًا) أن α يجب أن يكون جذرًا لكثير الحدود المميزأ0+أ1ت+أ2ت2++أنتن{\displaystyle a_{0}+a_{1}t+a_{2}t^{2}+\cdots +a_{n}t^{n}} من المعادلة التفاضلية، وهي الطرف الأيسر من المعادلة المميزةأ0+أ1ت+أ2ت2++أنتن=0.{\displaystyle a_{0}+a_{1}t+a_{2}t^{2}+\cdots +a_{n}t^{n}=0.}

عندما تكون هذه الجذور جميعها متميزة ، يكون لدينا n حلاً متميزاً، ليست بالضرورة حقيقية، حتى لو كانت معاملات المعادلة حقيقية. يمكن إثبات استقلال هذه الحلول خطياً ، من خلال حساب محدد فاندرموند لقيم هذه الحلول عند x = 0، ...، n – 1. تشكل هذه الحلول مجتمعةً أساساً للفضاء المتجهي لحلول المعادلة التفاضلية (أي نواة المؤثر التفاضلي).

مثال

y-2y+2y"-2y+y=0{\displaystyle y''''-2y'''+2y''-2y'+y=0} لها المعادلة المميزة z4-2z3+2z2-2z+1=0.{\displaystyle z^{4}-2z^{3}+2z^{2}-2z+1=0.} يحتوي هذا على أصفار، i و -i و 1 ( تعددية 2). وبالتالي، فإن أساس الحل هو هـأناx،هـ-أناx،هـx،xهـx.{\displaystyle e^{ix},\;e^{-ix},\;e^{x},\;xe^{x}.} وبالتالي فإن الأساس الحقيقي للحل هو كوسx،الخطيئةx،هـx،xهـx.{\displaystyle \cos x,\;\sin x,\;e^{x},\;xe^{x}.}

في حالة وجود جذور بسيطة فقط لكثير الحدود المميز ، فإن ما سبق يوفر أساسًا كاملاً لفضاء متجهات الحلول. أما في حالة وجود جذور متعددة ، فيلزم وجود حلول أكثر استقلالًا خطيًا لتكوين أساس. وتأخذ هذه الحلول الشكل التالي: xكهـαx،{\displaystyle x^{k}e^{\alpha x},} حيث k عدد صحيح غير سالب، وα جذر لكثير الحدود المميز ذي التعددية m ، و k < m . لإثبات أن هذه الدوال حلول، تجدر الإشارة إلى أنه إذا كان α جذرًا لكثير الحدود المميز ذي التعددية m ، فيمكن تحليل كثير الحدود المميز إلى P ( t )( tα ) m . وبالتالي، فإن تطبيق المؤثر التفاضلي على المعادلة يكافئ تطبيق المؤثر m مرة.ددx-α{\textstyle {\frac {d}{dx}}-\alpha }ثم المؤثر الذي تكون فيه P متعددة الحدود المميزة. وفقًا لنظرية الإزاحة الأسية ، (ددx-α)(xكهـαx)=كxك-1هـαx،{\displaystyle \left({\frac {d}{dx}}-\alpha \right)\left(x^{k}e^{\alpha x}\right)=kx^{k-1}e^{\alpha x},}

وبالتالي نحصل على الصفر بعد تطبيق k + 1 منددx-α{\textstyle {\frac {d}{dx}}-\alpha }.

وبما أن مجموع تعدد جذور كثير الحدود يساوي درجة كثير الحدود، فإن عدد الحلول المذكورة أعلاه يساوي رتبة المعادلة التفاضلية، وتشكل هذه الحلول أساسًا للفضاء المتجه للحلول .

في الحالة الشائعة حيث تكون معاملات المعادلة حقيقية، يكون من الأنسب عمومًا وجود أساس للحلول يتكون من دوال حقيقية . يمكن الحصول على هذا الأساس من الأساس السابق بملاحظة أنه إذا كان a + ib جذرًا لكثير الحدود المميز، فإن aib هو أيضًا جذر من نفس الرتبة. وبالتالي، يتم الحصول على أساس حقيقي باستخدام صيغة أويلر ، واستبدالxكهـ(أ+أناب)x{\displaystyle x^{k}e^{(a+ib)x}}وxكهـ(أ-أناب)x{\displaystyle x^{k}e^{(a-ib)x}}بواسطةxكهـأxكوس(بx){\displaystyle x^{k}e^{ax}\cos(bx)}وxكهـأxالخطيئة(بx){\displaystyle x^{k}e^{ax}\sin(bx)}.

حالة من الدرجة الثانية

يمكن كتابة معادلة تفاضلية خطية متجانسة من الرتبة الثانية y"+أy+بy=0،{\displaystyle y''+ay'+by=0,} ومتعددة الحدود المميزة لها هي ر2+أر+ب.{\displaystyle r^{2}+ar+b.}

إذا كان a و b عددين حقيقيين ، فهناك ثلاث حالات للحلول، وذلك بحسب المميز D = - 4b . في جميع الحالات الثلاث ، يعتمد الحل العام على ثابتين اختياريين c₁ و c₂ .

  • إذا كانت D > 0 ، فإن لكثير الحدود المميز جذرين حقيقيين مختلفين α و β . في هذه الحالة، يكون الحل العام هوج1هـαx+ج2هـβx.{\displaystyle c_{1}e^{\alpha x}+c_{2}e^{\beta x}.}
  • إذا كانت D = 0 ، فإن متعددة الحدود المميزة لها جذر مضاعف a /2 ، والحل العام هو(ج1+ج2x)هـ-أx/2.{\displaystyle (c_{1}+c_{2}x)e^{-ax/2}.}
  • إذا كانت D < 0 ، فإن متعددة الحدود المميزة لها جذران مترافقان مركبان α ± βi ، والحل العام هوج1هـ(α+βأنا)x+ج2هـ(α-βأنا)x،{\displaystyle c_{1}e^{(\alpha +\beta i)x}+c_{2}e^{(\alpha -\beta i)x},}والتي يمكن إعادة كتابتها بصيغة حقيقية، باستخدام صيغة أويلر كما يلي:هـαx(ج1كوس(βx)+ج2الخطيئة(βx)).{\displaystyle e^{\alpha x}(c_{1}\cos(\beta x)+c_{2}\sin(\beta x)).}

لإيجاد الحل y ( x ) الذي يحقق y (0) = d1 و y ′(0) = d2 ، نساوي قيم الحل العام المذكور أعلاه عند الصفر ومشتقته عند تلك النقطة بـ d1 و d2 على التوالي . ينتج عن ذلك نظام خطي من معادلتين خطيتين في المجهولين c1 و c2 . حل هذا النظام يعطي حلًا لما يُسمى بمسألة كوشي ، حيث تُحدد قيم الحل التفاضلي ومشتقته عند الصفر .

معادلة غير متجانسة ذات معاملات ثابتة

يمكن كتابة معادلة غير متجانسة من الرتبة n ذات معاملات ثابتة على النحو التاليy(ن)(x)+أ1y(ن-1)(x)++أن-1y(x)+أنy(x)=و(x)،{\displaystyle y^{(n)}(x)+a_{1}y^{(n-1)}(x)+\cdots +a_{n-1}y'(x)+a_{n}y(x)=f(x),} حيث أن a 1 ، ...، a n هي أعداد حقيقية أو مركبة، و f هي دالة معطاة لـ x ، و y هي الدالة المجهولة (من أجل التبسيط، سيتم حذف " ( x ) " فيما يلي).

توجد عدة طرق لحل هذه المعادلة. وتعتمد أفضل طريقة على طبيعة الدالة f التي تجعل المعادلة غير متجانسة. إذا كانت f عبارة عن توليفة خطية من دوال أسية وجيبية، فيمكن استخدام صيغة الاستجابة الأسية . أما إذا كانت f ، بشكل أعم، عبارة عن توليفة خطية من دوال على الصورة xⁿe⁻ᵏ⁺ ، و xⁿcos ( ax ) ، و xⁿsin ( ax ) ، حيث n عدد صحيح غير سالب، و a ثابت (لا يشترط أن يكون متساويًا في جميع الحدود)، فيمكن استخدام طريقة المعاملات غير المحددة . وبشكل أعم، تُطبق طريقة الإفناء عندما تحقق f معادلة تفاضلية خطية متجانسة ، وعادةً ما تكون دالة هولونومية .

الطريقة الأكثر عمومية هي تغيير الثوابت ، والتي يتم عرضها هنا.

الحل العام للمعادلة المتجانسة المرتبطة y(ن)+أ1y(ن-1)++أن-1y+أنy=0{\displaystyle y^{(n)}+a_{1}y^{(n-1)}+\cdots +a_{n-1}y'+a_{n}y=0} يكون y=u1y1++uنyن،{\displaystyle y=u_{1}y_{1}+\cdots +u_{n}y_{n},} حيث ( y₁ , ..., yn ) أساسٌ للفضاء المتجهي للحلول، و u₁ , ... , yn ثوابت اختيارية . تستمد طريقة تغيير الثوابت اسمها من الفكرة التالية: بدلاً من اعتبار u₁ , ..., yn ثوابت، يمكن اعتبارها دوال مجهولة يجب تحديدها لجعل y حلاً للمعادلة غير المتجانسة. ولتحقيق هذا الغرض، تُضاف القيود التالية :0=u1y1+u2y2++uنyن0=u1y1+u2y2++uنyن0=u1y1(ن-2)+u2y2(ن-2)++uنyن(ن-2)،{\displaystyle {\begin{aligned}0&=u'_{1}y_{1}+u'_{2}y_{2}+\cdots +u'_{n}y_{n}\\0&=u'_{1}y'_{1}+u'_{2}y'_{2}+\cdots +u'_{n}y'_{n}\\&\;\;\vdots \\0&=u'_{1}y_{1}^{(n-2)}+u'_{2}y_{2}^{(n-2)}+\cdots +u'_{n}y_{n}^{(n-2)},\end{aligned}}} مما يستلزم (بحسب قاعدة الضرب والاستقراء ) y(أنا)=u1y1(أنا)++uنyن(أنا){\displaystyle y^{(i)}=u_{1}y_{1}^{(i)}+\cdots +u_{n}y_{n}^{(i)}} لـ i = 1، ...، n – 1 ، و y(ن)=u1y1(ن)++uنyن(ن)+u1y1(ن-1)+u2y2(ن-1)++uنyن(ن-1).{\displaystyle y^{(n)}=u_{1}y_{1}^{(n)}+\cdots +u_{n}y_{n}^{(n)}+u'_{1}y_{1}^{(n-1)}+u'_{2}y_{2}^{(n-1)}+\cdots +u'_{n}y_{n}^{(n-1)}.}

باستبدال y ومشتقاتها في المعادلة الأصلية بهذه التعبيرات، وباستخدام حقيقة أن y1 ، ...، yn هي حلول للمعادلة المتجانسة الأصلية، نحصل علىو=u1y1(ن-1)++uنyن(ن-1).{\displaystyle f=u'_{1}y_{1}^{(n-1)}+\cdots +u'_{n}y_{n}^{(n-1)}.}

تشكل هذه المعادلة والمعادلات السابقة التي يكون طرفها الأيسر صفرًا نظامًا من n معادلة خطية بدلالة u′1 ، ...، u′n ، حيث معاملات هذه المعادلات دوال معلومة ( f ، و yᵢ ، ومشتقاتها). يمكن حل هذا النظام بأي طريقة من طرق الجبر الخطي . بحساب الدوال الأصلية، نحصل على u₁ ، ... ، uₙ ، ومن ثم y = u₁y₁ + ... + uₙyₙ .

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

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

الصيغة العامة لمعادلة تفاضلية خطية عادية من الرتبة 1، بعد قسمة معامل y ′( x ) ، هي: y(x)=و(x)y(x)+ز(x).{\displaystyle y'(x)=f(x)y(x)+g(x).}

إذا كانت المعادلة متجانسة، أي g ( x ) = 0 ، فيمكن إعادة كتابتها وتكاملها: yy=و،سجلy=ك+F،{\displaystyle {\frac {y'}{y}}=f,\qquad \log y=k+F,} حيث k ثابت تكامل اختياري وF=ودx{\displaystyle F=\textstyle \int f\,dx}أي دالة أصلية للدالة f . وبالتالي، فإن الحل العام للمعادلة المتجانسة هو y=جهـF،{\displaystyle y=ce^{F},} حيث c = e k ثابت اختياري.

بالنسبة للمعادلة غير المتجانسة العامة، من المفيد ضرب طرفي المعادلة بمقلوب e F لحل المعادلة المتجانسة. [ 2 ] وهذا يعطي yهـ-F-yوهـ-F=زهـ-F.{\displaystyle y'e^{-F}-yfe^{-F}=ge^{-F}.} كما-وهـ-F=ددx(هـ-F)،{\displaystyle -fe^{-F}={\tfrac {d}{dx}}\left(e^{-F}\right),}تسمح قاعدة الضرب بإعادة كتابة المعادلة على النحو التالي :ددx(yهـ-F)=زهـ-F.{\displaystyle {\frac {d}{dx}}\left(ye^{-F}\right)=ge^{-F}.} وبالتالي، فإن الحل العام هو y=جهـF+هـFزهـ-Fدx،{\displaystyle y=ce^{F}+e^{F}\int ge^{-F}dx,} حيث c هو ثابت التكامل، و F هو أي دالة أصلية لـ f (تغيير الدالة الأصلية يؤدي إلى تغيير ثابت التكامل).

نهج بديل

لنضرب طرفي المعادلة y-و(x)y=ز(x).{\displaystyle y'-f(x)y=g(x).}بواسطة دالةح(x){\displaystyle h(x)}ح(x)y-ح(x)و(x)y=ح(x)ز(x).{\displaystyle h(x)y'-h(x)f(x)y=h(x)g(x).}

نريد أن يكون الجانب الأيسر مشتقًا كليًا:د(ح(x)y(x))/دx{\displaystyle d(h(x)y(x))/dx}باستخدام قاعدة الضرب،د(ح(x)y)=ح(x)دy+yدح(x){\displaystyle d(h(x)y)=h(x)dy+ydh(x)}بالمقارنة مع المعادلة الأصلية،-ح(x)و(x)y=yح(x){\displaystyle -h(x)f(x)y=yh'(x)}أودح(x)/ح(x)=-و(x)دx{\displaystyle dh(x)/h(x)=-f(x)dx}أو بالتكامل نحصل على ح(x)=هـ-ودx{\displaystyle h(x)=e^{-\textstyle \int f\,dx}} والتكاملد(ح(x)y(x))/دx=ح(x)ز(x){\displaystyle d(h(x)y(x))/dx=h(x)g(x)} نحصل y(x)=ج+ح(x)ز(x)دxح(x){\displaystyle y(x)={\frac {C+\textstyle \int h(x)g(x)\,dx}{h(x)}}}

مثال

حل المعادلة y(x)+y(x)x=3x.{\displaystyle y'(x)+{\frac {y(x)}{x}}=3x.} المعادلة المتجانسة المرتبطةy(x)+y(x)x=0{\displaystyle y'(x)+{\frac {y(x)}{x}}=0}أعطِ yy=-1x،{\displaystyle {\frac {y'}{y}}=-{\frac {1}{x}},} إنه y=جx.{\displaystyle y={\frac {c}{x}}.}

قسمة المعادلة الأصلية على أحد هذه الحلول تعطي xy+y=3x2.{\displaystyle xy'+y=3x^{2}.} إنه (xy)=3x2،{\displaystyle (xy)'=3x^{2},}xy=x3+ج،{\displaystyle xy=x^{3}+c,} و y(x)=x2+ج/x.{\displaystyle y(x)=x^{2}+c/x.} بالنسبة للحالة الأولية y(1)=α،{\displaystyle y(1)=\alpha ,} يحصل المرء على الحل المحدد y(x)=x2+α-1x.{\displaystyle y(x)=x^{2}+{\frac {\alpha -1}{x}}.}

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

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

يمكن تحويل أي معادلة تفاضلية خطية عادية ونظام من هذه المعادلات إلى نظام من المعادلات التفاضلية الخطية من الدرجة الأولى عن طريق إضافة متغيرات لجميع المشتقات باستثناء أعلى رتبة. أي، إذاy،y"،...،y(ك){\displaystyle y',y'',\ldots ,y^{(k)}}إذا ظهرت في معادلة، فيمكن استبدالها بدوال مجهولة جديدة .y1،...،yك{\displaystyle y_{1},\ldots ,y_{k}}التي يجب أن تحقق المعادلاتy=y1{\displaystyle y'=y_{1}}وyأنا=yأنا+1،{\displaystyle y_{i}'=y_{i+1},}for i = 1, ..., k 1 .

يمكن عادةً حل نظام خطي من الدرجة الأولى، يحتوي على n دالة مجهولة و n معادلة تفاضلية، لإيجاد مشتقات الدوال المجهولة. إذا لم يكن الأمر كذلك، فهذا نظام تفاضلي جبري ، وهذه نظرية مختلفة. لذلك، فإن الأنظمة التي نتناولها هنا تأخذ الشكل التالي: y1(x)=ب1(x)+أ1،1(x)y1++أ1،ن(x)yنyن(x)=بن(x)+أن،1(x)y1++أن،ن(x)yن،{\displaystyle {\begin{aligned}y_{1}'(x)&=b_{1}(x)+a_{1,1}(x)y_{1}+\cdots +a_{1,n}(x)y_{n}\\[1ex]&\;\;\vdots \\[1ex]y_{n}'(x)&=b_{n}(x)+a_{n,1}(x)y_{1}+\cdots +a_{n,n}(x)y_{n},\end{aligned}}} أينبن{\displaystyle b_{n}}وأأنا،ج{\displaystyle a_{i,j}}هي دوال لـ x . في تدوين المصفوفات، يمكن كتابة هذا النظام (مع حذف " ( x ) "). y=أy+ب.{\displaystyle \mathbf {y} '=A\mathbf {y} +\mathbf {b} .}

طريقة الحل مشابهة لطريقة حل المعادلات التفاضلية الخطية من الدرجة الأولى، ولكن مع وجود تعقيدات ناتجة عن عدم تبادلية ضرب المصفوفات.

يترك u=أu.{\displaystyle \mathbf {u} '=A\mathbf {u} .} لتكن المعادلة المتجانسة المرتبطة بمعادلة المصفوفة أعلاه. تشكل حلولها فضاءً متجهيًا ذا بُعد n ، وبالتالي فهي أعمدة مصفوفة مربعة من الدوال .يو(x){\displaystyle U(x)}، والتي لا يكون محددها هو الدالة الصفرية. إذا كان n = 1 ، أو كانت A مصفوفة من الثوابت، أو بشكل أعم، إذا كانت A تتبادل مع مشتقتها الأصلية .ب=أدx{\displaystyle \textstyle B=\int Adx}عندئذٍ ، يمكن اختيار U مساويًا للدالة الأسية لـ B. في الواقع، في هذه الحالات، يكون لدى المرء ددxخبرة(ب)=أخبرة(ب).{\displaystyle {\frac {d}{dx}}\exp(B)=A\exp(B).} في الحالة العامة، لا يوجد حل مغلق للمعادلة المتجانسة، ويتعين على المرء استخدام إما طريقة عددية ، أو طريقة تقريبية مثل توسيع ماغنوس .

بمعرفة المصفوفة U ، يكون الحل العام للمعادلة غير المتجانسة هو y(x)=يو(x)y0+يو(x)يو-1(x)ب(x)دx،{\displaystyle \mathbf {y} (x)=U(x)\mathbf {y_{0}} +U(x)\int U^{-1}(x)\mathbf {b} (x)\,dx,} حيث مصفوفة الأعمدةy0{\displaystyle \mathbf {y_{0}} }هو ثابت تكامل اختياري .

إذا تم إعطاء الشروط الابتدائية على النحو التالي y(x0)=y0،{\displaystyle \mathbf {y} (x_{0})=\mathbf {y} _{0},} الحل الذي يحقق هذه الشروط الأولية هو y(x)=يو(x)يو-1(x0)y0+يو(x)x0xيو-1(ت)ب(ت)دت.{\displaystyle \mathbf {y} (x)=U(x)U^{-1}(x_{0})\mathbf {y_{0}} +U(x)\int _{x_{0}}^{x}U^{-1}(t)\mathbf {b} (t)\,dt.}

رتبة أعلى بمعاملات متغيرة

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

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

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

ومع ذلك، فقد تم حل حالة الرتبة الثانية ذات المعاملات النسبية بشكل كامل بواسطة خوارزمية كوفاسيتش .

معادلة كوشي-أويلر

تُعد معادلات كوشي-أويلر أمثلة على المعادلات من أي رتبة، ذات المعاملات المتغيرة، والتي يمكن حلها بشكل صريح. وهي معادلات من الشكل التالي: xنy(ن)(x)+أن-1xن-1y(ن-1)(x)++أ0y(x)=0،{\displaystyle x^{n}y^{(n)}(x)+a_{n-1}x^{n-1}y^{(n-1)}(x)+\cdots +a_{0}y(x)=0,} أينأ0،...،أن-1{\displaystyle a_{0},\ldots ,a_{n-1}}هي معاملات ثابتة.

الوظائف الشاملة

الدالة الهولونومية ، وتسمى أيضًا الدالة المحدودة D ، هي دالة تمثل حلاً لمعادلة تفاضلية خطية متجانسة ذات معاملات متعددة الحدود.

معظم الدوال الشائعة في الرياضيات هي دوال هولونومية أو نواتج قسمة دوال هولونومية. في الواقع، تشمل الدوال الهولونومية كثيرات الحدود ، والدوال الجبرية ، واللوغاريتمات ، والدوال الأسية ، ودوال الجيب ، وجيب التمام ، وجيب التمام الزائدي ، وجيب التمام الزائدي ، والدوال المثلثية العكسية ، والدوال الزائدية العكسية ، والعديد من الدوال الخاصة مثل دوال بيسل والدوال فوق الهندسية .

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

تتجلى فائدة مفهوم الدوال الهولونومية في نظرية زيلبرغر، والتي تلي ذلك. [ 3 ]

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

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

انظر أيضاً

مراجع

  1. جيرشنفيلد 1999، ص 9
  2. الدافع : قياسًا على أسلوب إكمال المربع، نكتب المعادلة على الصورة y ′ − fy = g ، ونحاول تعديل الطرف الأيسر ليصبح مشتقة. تحديدًا، نبحث عن "عامل تكامل" h = h ( x ) بحيث يؤدي الضرب به إلى جعل الطرف الأيسر مساويًا لمشتقة hy ، أي hy ′ − hfy = ( hy )′ . هذا يعني أن h ′ = −hf ، وبالتالي h = e −∫ f dx = e F ، كما هو موضح في النص.
  3. 1 2 3 زيلبرغر، دورون. منهج الأنظمة الهولونومية لهويات الدوال الخاصة . مجلة الرياضيات الحسابية والتطبيقية. 32.3 (1990): 321-368
  4. بينوا، أ.، تشيزاك، ف.، داراس، أ.، جيرهولد، س.، ميزاروبا، م.، وسالفي، ب. (سبتمبر 2010). القاموس الديناميكي للدوال الرياضية (DDMF) . في المؤتمر الدولي للبرمجيات الرياضية (ص 35-41). سبرينغر، برلين، هايدلبرغ.
  • بيركوف، غاريت وروتا، جيان كارلو (1978)، المعادلات التفاضلية العادية ، نيويورك: جون وايلي وأولاده، ISBN 0-471-07411-X
  • غيرشنفيلد، نيل (1999)، طبيعة النمذجة الرياضية ، كامبريدج، المملكة المتحدة: مطبعة جامعة كامبريدج، رقم ISBN 978-0-521-57095-4
  • روبنسون، جيمس سي. (2004)، مقدمة في المعادلات التفاضلية العادية ، كامبريدج، المملكة المتحدة: مطبعة جامعة كامبريدج، ISBN 0-521-82650-0