دستگاه معادلات دیفرانسیل خطی با ضرایب متغیر — از صفر تا صد
در مطالب قبلی از مجموعه آموزشهای ریاضیات مجله فرادرس، درباره دستگاه معادلات دیفرانسیل خطی با ضرایب ثابت بحث کردیم. در این آموزش، روش حل دستگاه معادلات دیفرانسیل خطی با ضرایب متغیر را بررسی میکنیم.
دستگاه معادلات دیفرانسیل خطی با ضرایب متغیر
یک دستگاه معادلات دیفرانسیل معمولی با ضرایب متغیر، بهصورت زیر نوشته میشود:
که در آن، توابع مجهولی هستند که در بازه پیوسته و مشتقپذیرند. ضرایب و جملات آزاد ، توابعی پیوسته در بازه هستند.
دستگاه معادلات را بهصورت ماتریسی-برداری زیر مینویسیم:
که در آن:
$$ \large { { \mathbf{X} } \left( t \right) = \left[ { \begin{array}{*{20}{c}}<br /> { { x _ 1 } \left( t \right)}\\<br /> { { x _ 2 } \left( t \right)}\\<br /> \vdots \\<br /> { { x _ n }\left( t \right)}<br /> \end{array}} \right],\;\;}\kern-0.3pt<br /> { { A \left( t \right) \text{ = }}\kern0pt{\left[ {\begin{array}{*{20}{c}}<br /> { { a _ { 1 1 } } \left( t \right)}&{ { a _ { 1 2 } } \left( t \right)}& \vdots &{{a_{1n}}\left( t \right)}\\<br /> { { a _ { 2 1 } } \left( t \right)}& { { a _{ 2 2 } } \left( t \right)}& \vdots & { { a _ { 2 n } } \left( t \right) }\\<br /> \cdots & \cdots & \cdots & \cdots \\<br /> { { a _ { n 1 } } \left( t \right)}&{ { a _ { n 2 } } \left( t \right) }& \vdots & { { a _ { n n } } \left( t \right)}<br /> \end{array}} \right],\;\;}}\kern-0.3pt<br /> { { \mathbf{f}}\left( t \right) = \left[ {\begin{array}{* { 2 0 } { c } }<br /> { { f _ 1 }\left( t \right)}\\<br /> { { f _ 2 }\left( t \right)}\\<br /> \vdots \\<br /> { { f _ n } \left( t \right)}<br /> \end{array}} \right].} $$
در حالت کلی، ماتریس و بردار توابع و دارای هر دو مقدار حقیقی و مختلط هستند.
دستگاه همگن متناظر با ضرایب متغیر بهفرم برداری زیر است:
دستگاه اساسی جوابها و ماتریس اساسی
توابع برداری ، ، و در بازه وابسته خطی هستند، اگر اعداد ، ، و همگی صفر نباشند و رابطه زیر برقرار باشد:
اگر این معادله فقط در شرایطِ
برقرار باشد، توابع برداری ، در بازه دادهشده مستقل خطی نامیده میشوند.
هر دستگاه با جواب مستقل خطی ، ، و یک دستگاه اساسی یا پایه از جوابها نامیده میشود.
ماتریس مربعی که ستونهای آن از جوابهای مستقل خطی ، ، و تشکیل شده، ماتریس اساسی دستگاه معادلات نامیده میشود و بهفرم زیر است:
$$ \large \require{cancel}<br /> \mathbf{ X ’ }\left( t \right) = A \left( t \right)\mathbf{X}\left( t \right) + \mathbf{f}\left( t \right),\;\;\Rightarrow<br /> {\cancel{ \Phi ’ \left( t \right)\mathbf{C}\left( t \right)} + \Phi \left( t \right)\mathbf{C’}\left( t \right) } \\ \large<br /> = {\cancel{A\left( t \right)\Phi \left( t \right)\mathbf{C}\left( t \right)} + \mathbf{ f } \left( t \right),\;\;}\Rightarrow<br /> {\Phi \left( t \right)\mathbf{C’}\left( t \right) = \mathbf{f}\left( t \right).} $$
که در آن، مؤلفههای بردار جوابهای مستقل خطی ، ، و هستند.
لازم به ذکر است که ماتریس اساسی غیرمنفرد است، یعنی ماتریس معکوس وجود دارد. از آنجایی که ماتریس اساسی، جواب مستقل خطی دارد، از جایگذاری آن در دستگاه همگن، رابطه زیر بهدست میآید:
معادله اخیر را از سمت راست، در ماتریس معکوس ضرب میکنیم:
رابطه منتجه، بهطور منحصربهفرد، یک دستگاه معادلات همگن را ارائه میکند که ماتریس اساسی را نتیجه میدهد.
جواب عمومی دستگاه همگن برحسب ماتریس اساسی، بهصورت زیر است:
که در آن، یک بردار بُعدی شامل اعداد دلخواه است.
در اینجا، حالت خاص دستگاه معادلات همگن را بررسی میکنیم. اگر ضرب ماتریس در انتگرال این ماتریس، جابهجاییپذیر باشد، یعنی:
ماتریس اساسی این دستگاه معادلات، با رابطه زیر بیان میشود:
این ویژگی، در ماتریسهای متقارن و بهطور خاص در ماتریسهای قطری برقرار است.
فرمول رونسکین و لیوویل
دترمینان ماتریس اساسی ، «رونسکین» (Wronskian) دستگاه جوابهای ، ، و نامیده میشود:
$$ \large {W\left( t \right) }={ W\left[ {{\mathbf{x}_1},{\mathbf{x}_2}, \ldots ,{\mathbf{x}_n}} \right] \text{ = }}\kern0pt<br /> {\left| {\begin{array}{*{20}{c}}<br /> {{x_{11}}\left( t \right)}&{{x_{12}}\left( t \right)}& \vdots &{{x_{1n}}\left( t \right)}\\<br /> {{x_{21}}\left( t \right)}&{{x_{22}}\left( t \right)}& \vdots &{{x_{2n}}\left( t \right)}\\<br /> \cdots & \cdots & \cdots & \cdots \\<br /> {{x_{n1}}\left( t \right)}&{{x_{n2}}\left( t \right)}& \vdots &{{x_{nn}}\left( t \right)}<br /> \end{array}} \right|.} $$
رونسکین، برای بررسی استقلال خطی جوابها مفید است. با استفاده از رونسکین میتوان موارد زیر را بیان کرد:
- جوابهای ، ، و دستگاه همگن، یک دستگاه اساسی تشکیل میدهند اگر و فقط اگر رونسکین متناظر، در هیچ نقطهای از بازه صفر نباشد.
- جوابهای ، ، و در بازه وابسته خطی هستند اگر و فقط اگر رونسکین در این بازه بهصورت تحلیلی صفر باشد.
رونسکین جوابهای ، ، و با «فرمول لیوویل» (Liouville’s Formula) محاسبه میشود:
که در آن، اثر ماتریس (یعنی مجموع درایههای قطر اصلی) است:
در حالتی که جواب خصوصی معلوم باشد، میتوان از فرمول لیوویل برای تشکیل جواب عمومی دستگاه همگن استفاده کرد.
روش تغییر ثابتها (روش لاگرانژ)
اکنون درباره دستگاههای ناهمگن بحث میکنیم که میتوان آنها را بهفرم برداری-ماتریسی زیر نوشت:
جواب عمومی چنین دستگاهی، مجموع جواب عمومی دستگاه همگن متناظر و یک جواب خصوصی دستگاه ناهمگن است. یعنی:
که در آن، یک ماتریس اساسی و یک بردار دلخواه است.
متداولترین روش برای حل دستگاههای ناهمگن، روش تغییر ثابتها (روش لاگرانژ) است. در این روش، بهجای بردار ثابت ، بردار در نظر گرفته میشود که مؤلفههای آن، توابع پیوسته مشتقپذیری از متغیر مستقل هستند، یعنی فرض میکنیم:
با جایگذاری رابطه اخیر در دستگاه ناهمگن، بردار مجهول بهدست میآید:
$$ \large \require{cancel}<br /> {\mathbf{X’}\left( t \right) = A\left( t \right)\mathbf{X}\left( t \right) + \mathbf{f}\left( t \right),\;\;}\Rightarrow<br /> {\cancel{\Phi’\left( t \right)\mathbf{C}\left( t \right)} + \Phi \left( t \right)\mathbf{C’}\left( t \right) } \\ \large<br /> = {\cancel{A\left( t \right)\Phi \left( t \right)\mathbf{C}\left( t \right)} + \mathbf{f}\left( t \right),\;\;}\Rightarrow<br /> {\Phi \left( t \right)\mathbf{C’}\left( t \right) = \mathbf{f}\left( t \right).} $$
از آنجایی که ماتریس نامنفرد است، معادله اخیر را از چپ در ضرب میکنیم:
بعد از انتگرالگیری، ماتریس بهدست میآید.
مثالها
در ادامه، چند مثال را بررسی میکنیم.
مثال ۱
دستگاه خطی معادلاتی با جوابهای زیر را بنویسید:
$$ \large { { \mathbf{ x } _ 1 }\left( t \right) = \left[ {\begin{array}{*{ 2 0 } { c } }<br /> 2\\<br /> t<br /> \end{array}} \right],\;\;}\kern-0.3pt<br /> {{\mathbf{ x } _ 2 }\left( t \right) = \left[ {\begin{array}{*{ 2 0 } { c } }<br /> t\\<br /> { { t ^ 2 } }<br /> \end{array}} \right],\;\;}\kern-0.3pt<br /> {t \ne 0.} $$
حل: در این مسئله، ماتریس اساسی دستگاه بهصورت زیر است:
$$ \large \Phi \left( t \right) = \left[ {\begin{array}{ * { 2 0 } { c } }<br /> 2&t\\<br /> t&{ { t ^ 2 } }<br /> \end{array}} \right]. $$
اکنون ماتریس معکوس را محاسبه میکنیم:
$$ \large { { \Delta \left( \Phi \right) = \left| { \begin{array}{* { 2 0 } { c } }<br /> 2&t\\<br /> t& { { t ^ 2 } }<br /> \end{array}} \right| }={ 2{t^2} – {t^2} = {t^2},\;\;}}\\ \large \Rightarrow<br /> {{\Phi ^ { – 1 } } \left( t \right) = \frac{ 1 }{ { \Delta \left( \Phi \right) } } C _ { i j } ^ T }<br /> = {\frac{ 1 }{ { { t ^ 2 } } } { \left[ {\begin{array} { * { 2 0 } { c } }<br /> { {t ^ 2 } } &{ – t}\\<br /> { – t }&2<br /> \end{array}} \right]^T} }\\ \large<br /> = {\frac{1}{ { {t ^ 2 } } } \left[ {\begin{array}{ * { 2 0 } { c } }<br /> { { t ^ 2 } } & { – t } \\<br /> { – t } & 2<br /> \end{array}} \right] }<br /> = {\left[ {\begin{array}{*{20}{c}}<br /> 1&{ – \frac{1}{t}}\\<br /> { – \frac{ 1 } { t } } &{\frac{2}{ { { t ^ 2 } } } }<br /> \end{array}} \right].} $$
در معادلات بالا، همسازههای متناظر با درایههای ماتریس اساسی هستند.
ماتریس ضرایب دستگاه معادلات بهصورت زیر است:
مشتق ماتریس اساسی برابر است با (درایه به درایه حساب میشود):
$$ \large { \Phi ’ } \left( t \right) = \left[ {\begin{array} { * { 2 0 } { c } }<br /> 0&1\\<br /> 1& { { 2 t } }<br /> \end{array}} \right]. $$
بنابراین، داریم:
$$ \large {A\left( t \right) }={ \left[ {\begin{array} { * { 2 0 } { c } }<br /> 0&1\\<br /> 1& { 2 t }<br /> \end{array}} \right]\left[ {\begin{array} { * { 2 0 } { c } }<br /> 1 & { – \frac{ 1 } { t } } \\<br /> { – \frac{1}{t } } & {\frac{2} { { { t ^ 2 } } } }<br /> \end{array}} \right] } \\ \large<br /> = {\left[ {\begin{array}{*{20}{c}}<br /> {0 – \frac{1}{t}}&{0 + \frac{2}{{{t^2}}}}\\<br /> {1 – 2}&{ – \frac{ 1 } { t } + \frac{ 4 } { t } }<br /> \end{array}} \right] }<br /> = {\left[ {\begin{array}{ * { 2 0 } { c } }<br /> { – \frac{ 1 } { t } } & { \frac{ 2 } { { { t ^ 2 } } } } \\<br /> { – 1 }&{\frac{ 3 } { t } }<br /> \end{array}} \right].} $$
در نتیجه، دستگاه معادلات که جوابهای آن، و است، بهصورت زیر خواهد بود:
مثال ۲
جواب عمومی دستگاه معادلاتِ
را با دانستن جواب زیر بیابید:
$$ \large { \mathbf{ X } _ 1 }\left( t \right) = \left[ {\begin{array}{ * { 2 0 } { c } }<br /> { { x _ 1 } \left( t \right)}\\<br /> { { y _ 1 } \left( t \right) }<br /> \end{array}} \right] = \left[ {\begin{array}{ * { 2 0 } { c } }<br /> 1\\<br /> t<br /> \end{array}} \right]. $$
حل: جواب مستقل خطی دوم را با تابع برداری زیر نمایش میدهیم:
$$ \large { \mathbf{ X } _ 2 } \left( t \right) = \left[ {\begin{array}{ * { 2 0 } { c } }<br /> { { x _ 2 } \left( t \right)}\\<br /> { { y _ 2 }\left( t \right)}<br /> \end{array}} \right] = \left[ {\begin{array}{ * { 2 0 } { c } }<br /> u\\<br /> v<br /> \end{array}} \right] $$
که شرایط اولیه، و است.
میتوانیم از فرمول لیوویل استفاده کنیم:
$$ \large {W\left( t \right) = \left| {\begin{array}{*{20}{c}}<br /> 1&u\\<br /> t&v<br /> \end{array}} \right| }<br /> = {\left| {\begin{array}{*{20}{c}}<br /> 1&0\\<br /> t&1<br /> \end{array}} \right|{e^{\,\int\limits_0^t {\left( {A\left( \tau \right)} \right)d\tau } }} } \\ \large<br /> = {1 \cdot {e^{\,\int\limits_0^t {\left( { – \tau + \tau } \right)d\tau } }} }<br /> = {{e^{\,\int\limits_0^t {0d\tau } }} }<br /> = {{e^0} = 1.} $$
در نتیجه، رابطه بین توابع مجهول و بهدست میآید:
معادله دوم دستگاه اصلی را در نظر بگیرید. با جایگذاری جواب ، داریم:
از معادله قبل جمله را میتوان بهصورت زیر نوشت:
با جایگذاری معادله اخیر در معادله دیفرانسیل تابع ، داریم:
$$ \large \require{cancel}<br /> {\frac{{dv}}{{dt}} = \left( {1 – {t^2}} \right)u + tv,\;\;}\Rightarrow<br /> {\frac{{dv}}{{dt}} = \left( {1 – {t^2}} \right)u + {t^2}u + t,\;\;} \\ \large\Rightarrow<br /> {\frac{{dv}}{{dt}} = u – \cancel{{t^2}u} + \cancel{{t^2}u} + t,\;\;}\Rightarrow<br /> {\frac{{dv}}{{dt}} = u + t,\;\;}\Rightarrow<br /> {t\frac{{dv}}{{dt}} = tu + {t^2}.} $$
با توجه به اینکه ، یک معادله دیفرانسیل مرتبه اول برای تابع بهدست میآید:
یا
ابتدا جواب متناظر با معادله همگن را پیدا میکنیم:
که در آن، یک عدد دلخواه است.
اکنون، جواب معادله ناهمگن را با استفاده از روش تغییر پارامتر تعیین میکنیم:
بعد از جایگذاری، توصیف مشتق را بهدست میآوریم:
با انتگرالگیری از معادله فوق، مقدار بهدست میآید:
تابع با فرمول زیر بیان میشود:
اکنون یافتن تابع ساده است:
بنابراین، جواب دوم دستگاه بهصورت زیر است:
$$ \large {{\mathbf{X}_2}\left( t \right) = \left[ {\begin{array}{*{20}{c}}<br /> {{x_2}\left( t \right)}\\<br /> {{y_2}\left( t \right)}<br /> \end{array}} \right] }<br /> = {\left[ {\begin{array}{*{20}{c}}<br /> u\\<br /> v<br /> \end{array}} \right] }<br /> = {\left[ {\begin{array}{*{20}{c}}<br /> t\\<br /> {{t^2} + 1}<br /> \end{array}} \right].} $$
جواب کلی را نیز میتوان بهصورت زیر نوشت:
$$ \large {\mathbf{X}\left( t \right) }={ {C_1}{\mathbf{X}_1}\left( t \right) + {C_2}{\mathbf{X}_2}\left( t \right) }<br /> = {{{C_1}\left[ {\begin{array}{*{20}{c}}<br /> 1\\<br /> t<br /> \end{array}} \right] }+{ {C_2}\left[ {\begin{array}{*{20}{c}}<br /> t\\<br /> {{t^2} + 1}<br /> \end{array}} \right],}} $$
که در آن، و ثابتهایی دلخواه هستند.
با سلام و وقت بخیر، یک دستگاه معادلات دیفرانسیل خطی مشابه بالا با ضرایب متغیر داریم، بدون هیچ اطلاعاتی، چگونه میتوانیم جواب عمومی را پیدا کنیم؟ با تشکر از جواب گویی