معادلات دیفرانسیل مرتبه بالا — به زبان ساده (+ دانلود فیلم آموزش رایگان)
۵۰۸۵ بازدید
آخرین بهروزرسانی: ۲۱ تیر ۱۴۰۲
زمان مطالعه: ۹۰ دقیقه
در این آموزش، معادلات دیفرانسیل خطی مرتبه nاُم را بررسی میکنیم. از آنجایی که حل این معادلات، تعمیمی از حل معادلات دیفرانسیل مرتبه اول و مرتبه دوم است، پیشنهاد میکنیم قبل از خواندن این مطلب، آموزشهای «معادلات دیفرانسیل مرتبه اول» و «معادلات دیفرانسیل مرتبه دوم» را مطالعه کنید. در این آموزش، روشهای اصلی حل معادلات مرتبه بالا را بیان میکنم. توجه کنید که منظور از مرتبه بالا، مراتب بالاتر از ۲ است.
قضیه زیر بیان میکند که میتوان یک حل منحصربهفرد برای یک مسئله مقدار اولیه، دو رابطه اخیر ارائه کرد.
قضیه ۱
فرض کنید توابع p0,p1,…,pn−1 و g(t) همه در بازه باز I که شامل t0 است، پیوسته باشند. آنگاه یک حل منحصربهفرد برای مسئله مقدار اولیه (۲) و (۳) در همه tهای بازه I وجود دارد.
این قضیه، تعمیمی ساده از همان قضیهای است که برای معادلات مرتبه اول بیان کردیم.
ابتدا درباره معادله دیفرانسیل مرتبه n همگن بحث میکنیم. بنابراین، معادله همگن زیر را در نظر بگیرید:
y(n)+pn−1(t)y(n−1)+⋯+p1(t)y′+p0(t)y=0(4)
فرض کنید همه n پاسخ y1(t),y2(t),…,yn(t) معادله (۴) را میدانیم؛ آنگاه با تعمیمی اصل برهمنهی داریم:
y(t)=c1y1(t)+c2y2(t)+⋯+cnyn(t)
که حل معادله (۴) است. حال این پرسش مطرح میشود که آیا جواب بالا حل عمومی معادله دیفرانسیل مرتبه n است؟
برای آنکه این جواب، عمومی باشد، باید بتوانیم ثابتهای c1,c2,…,cn را برای هر t0 (که در بازه I برای قضیه ۱ بیان شد) و هر y1,y2,…,yn بهدست آوریم. بهعبارت دیگر، باید ضرایب c1,c2,…,cn را پیدا کنیم که در معادله زیر صدق کنند:
فرض کنید توابع p0,p1,…,pn−1 در بازه باز I پیوسته هستند. همچنین فرض کنید توابع y1(t),y2(t),…yn(t) حل معادله (۴) باشند. اگر W(y1,y2,…yn)(t)=0 آنگاه توابع y1(t),y2(t),…yn(t) یک مجموعه پایه برای پاسخ معادله تشکیل میدهند و جواب عمومی معادله (۴) بهصورت زیر خواهد بود:
y(t)=c1y1(t)+c2y2(t)+⋯+cnyn(t)
لازم به یادآوری است اگر مجموعهای از جوابها، یک مجموعه پایه جواب را تشکیل دهند، مجموعهای از توابع مستقل خطی خواهند بود. طبق قضیه زیر میتوان جواب معادلات دیفرانسیل ناهمگن مرتبه n را ارائه کرد.
قضیه ۳
فرض کنید Y1(t) و Y2(t) دو جواب معادله (۲)، و توابع y1(t),y2(t),…yn(t) یک مجموعه پایه از جوابهای معادله دیفرانسیل ناهمگن (۴) باشند. آنگاه،
Y1(t)−Y2(t)
یک جواب برای معادله (۴) بوده و میتوان آن را بهصورت زیر نوشت:
Y1(t)−Y2(t)=c1y1(t)+c2y2(t)+⋯+cnyn(t)
اکنون اگر Y(t) را بهعنوان جواب عمومی معادله (۲) و YP(t) را هر جواب دیگری برای آن در نظر بگیریم، با استفاده از قضیه فوق میتوان گفت:
فرض میکنیم حل معادله بالا بهفرم y(t)=ert باشد و آن را در معادله دیفرانسیل جایگذاری میکنیم. با کمی سادهسازی داریم:
ert(anrn+an−1rn−1+⋯+a1r+a0)=0
برای آنکه معادله بالا برقرار باشد، باید تساوی زیر را داشته باشیم:
anrn+an−1rn−1+⋯+a1r+a0=0
معادله بالا، چندجملهای یا معادله مشخصه نامیده میشود و با یافتن ریشههای آن میتوان پاسخ معادله دیفرانسیل را بهدست آورد. میدانیم که یک معادله دیفرانسیل مرتبه nاُم، تعداد n ریشه (با احتساب ریشههای تکراری) دارد.
در معادله مشخصه بالا، با سه دسته پاسخ مواجه هستیم: اول ریشههای مجزا و حقیقی، دوم ریشههای تکراری و سوم، ریشههای مختلط.
فرض کنید یک معادله دیفرانسیل مرتبه نهم داریم. ممکن است ترکیب ریشهها اینگونه باشد: ۳ ریشه حقیقی مجزا، یک ریشه که چهار بار تکرار شده و ۲ ریشه مزدوج مختلط. لازم به ذکر است که چون ضرایب معادله حقیقی هستند، اگر ریشه مختلط وجود داشته باشد، همراه با مزدوج آن خواهد بود.
بنابراین، برای هر معادله دیفرانسیل مرتبه nاُم باید n تابع مستقل خطی (یعنی یک مجموعه پایه از جوابها) تشکیل دهیم. البته قبل از این کار باید ماتریس رونسکین را بررسی کنیم که صفر نباشد. در ادامه، جوابهای متناظر با هر دسته از ریشههای معادله مشخصه را بیان میکنیم.
(۱) ریشههای حقیقی و مجزا: فرض کنید ریشههای معادله مشخصه بهصورت r1,r2,…,rk و مجزا از هم باشند. جواب متناظر با هریک از این ریشهها بهصورت زیر است:
er1t,er2t,⋯,erkt
(۲) ریشههای حقیقی تکراری: اکنون فرض میکنیم ریشه تکراری داشته باشیم. این مورد نیز مشابه موردی است که در مورد معادلات دیفرانسیل مرتبه دوم گفتیم. فرض کنید r ریشهای است که k بار تکرار شده است. بنابراین، تعداد k پاسخ بهصورت زیر داریم:
ert,tert,⋯,tk−1ert
(۳) ریشههای مختلط: حال پاسخ مربوط به ریشههای مزدوج مختلط را بررسی میکنیم. فرض کنید ریشههای مزدوج مختلط r=λ±μi را داریم که تکرار هم نشدهاند. در این حالت، دو پاسخ استاندارد داریم:
eλtcos(μt)eλtsin(μt)
(۴) ریشههای مختلط تکراری: مسئلهای که در این مورد وجود دارد، این است که ریشههای تکراری مختلط در معادلات دیفرانسیل مرتبه چهارم به بالا وجود دارند. اکنون فرض کنید r=λ±μi به تعداد k بار تکرار شده باشد. در این حالت، مشابه ریشههای تکراری حقیقی عمل میکنیم و 2k پاسخ مختلط خواهیم داشت:
توجه کنید که ما پاسخهای حقیقی را میخواهیم. بنابراین از فرمول اویلر استفاده میکنیم و قسمتهای حقیقی و موهومی عبارات بالا را بهصورت جداگانه مینویسیم:
قبل از اینکه چند مثال را بررسی کنیم، باید نکتهای را یادآور شویم و آن این است که یافتن ریشه چندجملهایهایی با درجه بالاتر از ۳، بهصورت دستی، کار دشواری است. بنابراین برای حل آنها باید از نرمافزاهای مناسب کمک بگیریم.
با توجه به معادله بالا، دو ریشه r1=21 و r2=−2 داریم که دومی ۳ بار تکرار شده است. طبق مواردی که در بالا گفتیم، وقتی ریشه تکراری داریم توانهایی از t در پاسخ ظاهر میشود. بنابراین، حل عمومی معادله بهصورت زیر خواهد بود:
همانطور که از معادله بالا مشخص است، یک ریشه حقیقی r=1، دو ریشه تکراری r=5 و دو ریشه مزدوج مختلط r=2±i داریم. بنابراین، حل معادله بالا بهصورت زیر است:
حل معادله بالا، ساده است. برای یافتن ریشههای آن باید ریشههای چهارم −16 را محاسبه کنیم. برای آشنایی بیشتر با روش یافتن ریشهها میتوانید آموزش «توان و ریشه اعداد مختلط — از صفر تا صد» را مطالعه کنید. ریشههای چهارم عدد −16 بهصورت زیر قابل محاسبه است:
اگر g(t) یک تابع نمایی، چندجملهای، سینوسی، کسینوسی یا جمع و ضریبی از این توابع باشد، باید شکل آن را مشابه معادله مرتبه دوم، حدس بزنیم. پس از آن میتوانیم ضرایب را محاسبه کنیم.
مثال ۵
معادله دیفرانسیل زیر را حل کنید.
y(3)−12y′′+48y′−64y=12−32e−8t+2e4t
حل: ابتدا جواب عمومی را محاسبه میکنیم. بنابراین، معادله مشخصه بهصورت زیر خواهد بود:
r3−12r2+48r−64=(r−4)3=0⇒r=4
در نتیجه، یک ریشه حقیقی داریم که سه بار تکرار شده است. بنابراین، جواب مکمل بهصورت زیر خواهد بود:
yc(t)=c1e4t+c2te4t+c3t2e4t
اکنون جواب خصوصی را حدس میزنیم:
YP=A+Be−8t+Ce4t
توجه کنید که جمله آخر معادله بالا، در پاسخ عمومی نیز وجود دارد. بنابراین میتوان چنین تعبیر کرد که چهار ریشه r=4 داریم. در مورد ریشههای تکراری گفتیم که ضریبی از توان t را باید در آنها ضرب کرد. بنابراین، جواب خصوصی را بهصورت زیر تصحیح میکنیم:
YP=A+Be−8t+Ct3e4t
برای یافتن ضرایب، سه بار از معادله فوق مشتق میگیریم و آن را برابر با سمت راست معادله دیفرانسیل اصلی قرار میدهیم:
یکی دیگر از راههای حل معادلات، استفاده از تبدیل لاپلاس است. این روش را در قالب یک مثال توضیح میدهیم.
مثال ۶
مسئله مقدار اولیه زیر را حل کنید.
y′′′−4y′′=4t+3u6(t)e30−5t,y(0)=−3y′(0)=1y′′(0)=4
حل: ابتدا باید بررسی کنیم جابهجایی زمانی تابع پله u(t) صحیح باشد:
y′′′−4y′′=4t+3u6(t)e−5(t−6)
میبینیم که جابهجایی e−5t بهدرستی انجام شده است. اکنون باید از معادله تبدیل لاپلاس بگیریم. همانطور که میدانیم، تبدیل لاپلاس معادله مرتبه سوم به صورت زیر است:
اکنون تبدیل لاپلاس جواب مورد نظر را در اختیار داریم. برای بهدست آوردن جواب، باید F(s) و G(s) را به کسرهای جزئی تجزیه کرده و عکس تبدیل لاپلاس آنها را محاسبه کنیم:
سید سراج حمیدی دانشآموخته مهندسی برق است و به ریاضیات و زبان و ادبیات فارسی علاقه دارد. او آموزشهای مهندسی برق، ریاضیات و ادبیات مجله فرادرس را مینویسد.