معادله دیفرانسیل خطی — از صفر تا صد (+ دانلود فیلم آموزش رایگان)

۱۷۵۵۰ بازدید
آخرین به‌روزرسانی: ۱۳ اردیبهشت ۱۴۰۲
زمان مطالعه: ۴۴ دقیقه
دانلود PDF مقاله
معادله دیفرانسیل خطی — از صفر تا صد (+ دانلود فیلم آموزش رایگان)

در راستای ارائه مفاهیم مرتبط با معادلات دیفرانسیل، روش‌های حل این معادلات اعم از مرتبه اول، دوم و مراتب بالاتر توضیح داده ‌شدند. با این حال در این مطلب قصد داریم تا معادله دیفرانسیل خطی مرتبه اول را مورد بررسی قرار دهیم. البته مطالعه مطالب معادلات دیفرانسیل، معادلات مرتبه دوم، معادلات دیفرانسیل مرتبه بالاتر و معادلات ناهمگن نیز به منظور درک هرچه بهتر این مطلب، خالی از لطف نخواهد بود.

997696
محتوای این مطلب جهت یادگیری بهتر و سریع‌تر آن، در انتهای متن به صورت ویدیویی نیز ارائه شده است.

پاسخ معادله دیفرانسیل خطی

یک معادله دیفرانسیل مرتبه اول می‌تواند خطی یا غیر خطی باشد. همان‌طور که پیش‌تر نیز اشاره شده، اگر در یک معادله دیفرانسیل توابعِ y یا مشتقاتشان در هم ضرب شده باشند، آن معادله غیر خطی خواهد بود. برای نمونه رابطه زیر، یک معادله دیفرانسیل غیرخطی را نشان می‌دهد.

y2+y+x2+1=0 \large {y ^ { \prime }} ^ 2 + y + x ^ 2 + 1 = 0

دلیل غیر خطی بودن رابطه فوق وجود ترمِ y2 \large {y ^ { \prime } } ^ 2 است. در این مطلب قصد داریم تا معادلات دیفرانسیل خطی مرتبه اول را حل کنیم. شکل کلی این معادلات به صورت زیر است.

dydt+p(t)y=g(t) \large \begin {equation} \frac { { d y } } { { d t } } + p \left ( t \right ) y = g \left ( t \right ) \end {equation}
رابطه ۱

توجه داشته باشید که در رابطه فوق، هر دو تابع (p(t و (g(t پیوسته هستند.

پاسخ معادله دیفرانسیل خطی

به منظور بدست آوردن پاسخ معادله ۱، طرفین آن را در تابعی جادویی - که فعلا آن را نمی‌شناسیم - ضرب می‌کنیم. با فرض این که اسم این تابع μ(t) \mu ( t ) باشد، به معادله زیر می‌رسیم. البته در ادامه این تابع را عامل انتگرال‌ساز صدا خواهیم زد.

μ(t)dydt+μ(t)p(t)y=μ(t)g(t) \large \begin {equation} \mu \left ( t \right ) \frac { { d y } } { { d t } } + \mu \left ( t \right ) p \left ( t \right ) y = \mu \left ( t \right ) g \left ( t \right ) \end {equation}

فرض کنید تابع μ(t) \mu ( t ) به شکلی است که رابطه زیر برای آن برقرار است.

μ(t)p(t)=μ(t) \large \begin {equation} \mu \left ( t \right ) p \left ( t \right ) = \mu ^ { \prime } \left ( t \right ) \end {equation}

در ادامه نحوه بدست آوردن μ(t) \mu ( t ) را توضیح خواهیم داد. با در نظر گرفتن فرض فوق، رابطه زیر بدست می‌آید.

μ(t)dydt+μ(t)y=μ(t)g(t) \large \begin {equation} \mu \left ( t \right ) \frac { { d y } } { { d t } } + \mu ^ { \prime } \left ( t \right ) y = \mu \left ( t \right ) g \left ( t \right ) \end {equation}
رابطه ۲

با استفاده از قوانین مشتق‌گیری می‌دانید که سمت چپ معادله بالا، برابر با مشتق حاصل‌ضرب μ(t)y(t) \mu ( t ) y ( t ) است. بنابراین سمت چپ معادله برابر است با:

μ(t)dydt+μ(t)y=(μ(t)y(t)) \large \mu \left ( t \right ) \frac { { d y } } { { d t } } + \mu ^ { \prime } \left ( t \right ) y = { \left ( { \mu \left ( t \right ) y \left ( t \right ) } \right ) ^ \prime }

با استفاده از رابطه فوق و رابطه ۲، عبارت زیر بدست می‌آید.

(μ(t)y(t))=μ(t)g(t) \large \begin {equation} { \left ( { \mu \left ( t \right ) y \left ( t \right ) } \right ) ^ \prime } = \mu \left ( t \right ) g \left ( t \right ) \end {equation}

حال به منظور بدست آوردن تابع y از طرفین رابطه فوق انتگرال می‌گیریم. با انجام این کار داریم:

(μ(t)y(t))dt=μ(t)g(t)dt \large \int { { { { \left ( { \mu \left ( t \right ) y \left ( t \right ) } \right ) } ^ \prime } \, d t } } = \int { { \mu \left ( t \right ) g \left ( t \right ) \, d t } }

μ(t)y(t)+c=μ(t)g(t)dt \large \begin {equation} \mu \left ( t \right ) y \left ( t \right ) + c = \int { { \mu \left ( t \right ) g \left ( t \right ) \, d t } } \end {equation}

توجه داشته باشید که در هنگام انتگرال‌گیری ثابت c باید در نظر گرفته شود. در غیر این‌صورت به پاسخی اشتباه خواهیم رسید. با توجه به عبارت بدست آمده در بالا، تابع y به صورت زیر بدست می‌آید.

μ(t)y(t)=μ(t)g(t)dtcy(t)=μ(t)g(t)dtcμ(t) \large \begin {align*} \mu \left ( t \right ) y \left ( t \right ) & = \int { { \mu \left ( t \right ) g \left ( t \right ) \, d t } } - c \\ y \left ( t \right ) & = \frac { { \int { { \mu \left ( t \right ) g \left ( t \right ) \, d t } } - c } } { { \mu \left ( t \right ) } } \end {align*}

مقدار c- ناشناخته است، بنابراین می‌توان آن را به صورت ثابتی نامعلوم در نظر گرفت. نهایتا پاسخ یک معادله دیفرانسیل خطی مرتبه اول، برابر با عبارت زیر بدست می‌آید.

y(t)=μ(t)g(t)dt+cμ(t) \large \boxed { \begin {equation} y \left ( t \right ) = \frac { { \int { { \mu \left ( t \right ) g \left ( t \right ) \, d t } } + c } } { { \mu \left ( t \right ) } } \end {equation} }
رابطه ۳

البته همان‌طور که می‌دانید مثبت یا منفی فرض کردن c، تاثیری در پاسخ نهایی نخواهد داشت.

بدست آوردن عامل انتگرال‌ساز (μ(t) \large \mu ( t ) )

در بالا متوجه شدیم که با ضرب یک معادله دیفرانسیل خطی مرتبه اول در تابعی جادویی تحت عنوان μ(t) \mu ( t ) می‌توان پاسخ y را با استفاده از فرمولی کلی بر حسب آن بدست آورد.

حال سوال این است تابع μ(t) \large \mu ( t ) به چه شکل بدست می‌آید؟ بدین منظور از فرض بیان شده در بالا استفاده می‌کنیم.

μ(t)p(t)=μ(t) \large \mu \left ( t \right ) p \left ( t \right ) = \mu ^ { \prime } \left ( t \right )

با تقسیم کردن طرفین رابطه فوق به μ(t) \mu ( t ) به رابطه زیر می‌رسیم.

μ(t)μ(t)=p(t) \large \frac { { \mu ^ { \prime } \left ( t \right ) } } { { \mu \left ( t \right ) } } = p \left ( t \right )

سمت چپ رابطه بالا معادل با (ln μ(t)) ( \ln { \ {\mu (t)} } )' است. بنابراین می‌توان رابطه فوق را به صورت زیر بازنویسی کرد.

(lnμ(t))=p(t) \large { \left ( { \ln \mu \left ( t \right ) } \right ) ^ \prime } = p \left ( t \right )

با انتگرال‌گیری، رابطه زیر بدست می‌آید.

lnμ(t)+k=p(t)dtlnμ(t)=p(t)dt+k \large \begin {align*} \ln \mu \left ( t \right ) + k &= \int { { p \left ( t \right ) \, d t } } \\ \ln \mu \left ( t \right ) & = \int { { p \left ( t \right ) \, d t } } + k \end {align*}

با نوشتن طرفین رابطه فوق به صورت نمایی، تابع μ(t) \mu ( t ) بدست خواهد آمد.

μ(t)=ep(t)dt+k \large \mu \left ( t \right ) = { { \bf{e} } ^ { \int { { p \left ( t \right ) \, d t } } \, + \, k } }

با استفاده از قوانین لگاریتمی نهایتا تابع μ(t) \mu ( t ) به صورت زیر بدست می‌آید.

μ(t)=kep(t)dt \large \begin {equation} \mu \left ( t \right ) = k \, { { \bf { e } } ^ { \int { { p \left ( t \right ) \, d t } } } } \end {equation}

پاسخ نهایی معادله دیفرانسیل خطی

با بدست آمدن تابع μ(t) \mu ( t ) ، پاسخ معادله نیز به صورت زیر بدست می‌آید.

y(t)=kep(t)dtg(t)dt+ckep(t)dt =kep(t)dtg(t)dt+ckep(t)dt =ep(t)dtg(t)dt+ckep(t)dt \large \large \begin {align*} y \left ( t \right ) & = \frac { { \int { { k \, { { \bf{ e } } ^ { \int { { p \left ( t \right ) \, d t } } } } g \left ( t \right ) \, d t } } + c } } { { k \, { { \bf{e} } ^ { \int { { p \left( t \right ) \, d t } } } } } } \\~\\ & = \frac { { k \int { { \, { { \bf{e}} ^ { \int { { p \left ( t \right ) \, d t } } } } g \left ( t \right ) \, d t } } + c } } { { k \, { { \bf{e}} ^ { \int { { p \left ( t \right ) \, d t } } } } } } \\~\\ & = \frac { { \int { { \,{ { \bf{e}} ^ { \int { { p \left ( t \right ) \, d t } } } } g \left ( t \right ) \, d t } } + \frac { c } { k } } } { { \,{ { \bf{e}} ^ { \int { { p \left ( t \right)\, d t } } } } } } \end {align*}

با فرض کردن c/k به عنوان یک ثابت، نهایتا پاسخ معادله دیفرانسیل خطی مرتبه اول، به صورت زیر قابل بیان است.

y(t)=μ(t)g(t)dt+cμ(t) \large \boxed { \begin {equation} y \left ( t \right ) = \frac { { \int { { \mu \left ( t \right ) g \left ( t \right ) \, d t } } + c } } { { \mu \left( t \right ) } } \end {equation} }
رابطه ۴

در رابطه فوق تابع μ(t) \mu ( t ) ، عامل انتگرال‌ساز نامیده شده و با استفاده از رابطه زیر محاسبه می‌شود.

μ(t)=ep(t)dt \large \begin {equation} \mu \left ( t \right ) = \, { { \bf{e } } ^ { \int { { p\left ( t \right ) \, d t } } } } \end {equation}

مراحل حل معادله دیفرانسیل خطی مرتبه اول

به منظور بدست آوردن پاسخ یک معادله دیفرانسیل خطی مرتبه اول، به ترتیب زیر عمل کنید:

  1. معادله را به شکل استاندارد (در قالب رابطه ۱) بنویسید.
  2. عامل انتگرال سازِ μ(t) \mu ( t ) را بدست آورید.
  3. کل معادله را در عامل انتگرال‌ساز ضرب کرده و نشان دهید که سمت چپ از قانون مشتقِ (μ(t)y(t)) \left ( { \mu \left ( t \right ) y \left ( t \right ) } \right ) ^ { \prime } پیروی می‌کند.
  4. از طرفین عبارت بدست آمده، انتگرال بگیرید. در این مرحله در مورد ثابت‌ها دقت کنید.
  5. پس از انتگرال‌گیری y را بدست آورید.

در ادامه مثال‌هایی ذکر شده که مطالعه آن‌ها را توصیه می‌کنیم.

مثال ۱

پاسخ معادله دیفرانسیل خطی زیر را بدست آورید.

dvdt=9.80.196v \large \frac { { d v } } { { d t } } = 9 .8 - 0 . 1 9 6 v

در ابتدا، معادله را به شکل استاندارد و مطابق با رابطه زیر بیان می‌کنیم:

dvdt+0.196v=9.8 \large \frac { { d v } } { { d t } } + 0.196 v = 9.8

با توجه به رابطه فوق p(t)=0.196 p(t)=0.196 بوده، در نتیجه عامل انتگرال‌ساز برابر است با:

μ(t)=e0.196dt=e0.196t\large \mu \left( t \right) = {{\bf{e}} ^ { \int { { 0.196dt} } } } = { { \bf{e} } ^ { 0 .196 t } }

با ضرب کردن طرفین معادله در عامل انتگرال‌ساز، می‌توان دید که سمت چپ عبارت بدست آمده برابر با مشتق تابع e0.196tv e ^ { 0.196t } v است.

e0.196tdvdt+0.196e0.196tv=9.8e0.196t(e0.196tv)=9.8e0.196t \large \begin {align*} { { \bf{e} } ^ { 0.196 t } } \frac { { d v } } { { d t } } + 0.196 { {\bf{e} } ^ { 0.196t } } v & = 9.8 { { \bf{e } } ^ { 0.196 t } } \\ { \left( { { { \bf{e } } ^ { 0.196t } } v } \right ) ^ \prime } & = 9.8 { {\bf{e} } ^ {0.196t } } \end {align*}

با انتگرال‌گیری از طرفین عبارت فوق، داریم:

(e0.196tv)dt=9.8e0.196tdte0.196tv+k=50e0.196t+c \large \begin {align*} \int { { { { \left ( { { { \bf{e} } ^ { 0.196 t } } v } \right ) } ^ \prime } \, d t } } & = \int { { 9.8 { { \bf{e} } ^ { 0.196 t } } \, d t } } \\ { { \bf{e } } ^ {0.196t} } v + k & = 50 { { \bf{e}} ^ { 0.196 t } } + c \end{align*}

توجه داشته باشید که برآیند دو مقدار c و k را می‌توان برابر با ثابت c در نظر گرفت. بنابراین نهایتا تابع (v(t برابر با عبارت زیر بدست می‌آید.

v(t)=50+ce0.196t \large v \left ( t \right ) = 5 0 + c { { \bf{e} } ^ { - 0.196 t } }

با توجه به رابطه فوق متوجه می‌شویم که چرا ضریب ثابت c دارای اهمیت است. در شکل زیر نمودار تابع به ازای چند c ترسیم شده است.

معادله دیفرانسیل خطی

در حقیقت ثابت c با استفاده از شرایط اولیه یا شرایط مرزی بدست می‌آیند. در ادامه مثالی ذکر شده که در آن ضریب ثابت نیز بدست می‌آید.

مثال ۲

پاسخ معادله دیفرانسیل خطی زیر را بیابید.

dvdt=9.80.196v,  v(0)=48 \large \frac { { d v } } { { d t } } = 9.8 - 0.196 v \hspace {0.25in} , \ \ v \left ( 0 \right) = 48

مسئله دقیقا مشابه با مثال ۱ بوده و تنها تفاوت آن‌ها در ارائه شرایط اولیه است. بنابراین پاسخِ معادله برابر است با:

v(t)=50+ce0.196t \large v \left ( t \right ) = 5 0 + c { { \bf{e} } ^ { - 0.196 t } }

بنابراین ضریب ثابت در t=0 به صورت زیر بدست می‌آید.

48=v(0)=50+cc=2 \large 48 = v \left( 0 \right ) = 50 + c \hspace {0.25in} \Rightarrow \hspace {0.25in} c = - 2

بنابراین پاسخ نهایی برابر است با:

v=502e0.196t \large v = 50 - 2{{\bf{e}}^{ - 0.196t}}

مثال ۳

پاسخ معادله دیفرانسیل خطی زیر را با توجه به شرط اولیه ارائه شده بیابید.

cos(x)y+sin(x)y=2cos3(x)sin(x)1   :y(π4)=32,0x<π2 \large \cos \left ( x \right ) y ^ { \prime } + \sin \left( x \right)y = 2{\cos ^3} \left( x \right)\sin \left( x \right) - 1 \ \ \ : \hspace {0.25in} y \left ( {\frac { \pi } { 4 } } \right ) = 3 \sqrt 2 ,\hspace{0.25in}\,\,\,\,\,\,0 \le x < \frac {\pi } { 2 }

در ابتدا معادله به فرم استاندارد نوشته می‌شود.

y+sin(x)cos(x)y=2cos2(x)sin(x)1cos(x)y+tan(x)y=2cos2(x)sin(x)sec(x) \large \begin{align*}y' + \frac{{\sin \left( x \right)}}{{\cos \left( x \right)}}y & = 2{\cos ^2}\left( x \right)\sin \left( x \right) - \frac{1}{{\cos \left( x \right)}}\\ y' + \tan \left( x \right)y &= 2{\cos ^2}\left( x \right)\sin \left( x \right) - \sec \left( x \right)\end{align*}

در مرحله بعد، عامل انتگرال‌ساز به صورت زیر بدست می‌آید.

μ(t)=etan(x)dx=elnsec(x)=elnsec(x)=sec(x) \large \mu \left( t \right) = {{\bf{e}}^{\int{{\tan \left( x \right)\,dx}}}} = {{\bf{e}}^{\ln \left| {\sec \left( x \right)} \right|}} = {{\bf{e}}^{\ln \,\,\sec \left( x \right)}} = \sec \left( x \right)

با ضرب عامل انتگرال‌ساز در طرفین معادله‌ی اصلی، به عبارت زیر می‌رسیم.

sec(x)y+sec(x)tan(x)y=2sec(x)cos2(x)sin(x)sec2(x)(sec(x)y)=2cos(x)sin(x)sec2(x) \large \begin{align*}\sec \left( x \right)y' + \sec \left( x \right)\tan \left( x \right)y & = 2\sec \left( x \right){\cos ^2}\left( x \right)\sin \left( x \right) - {\sec ^2}\left( x \right)\\ {\left( {\sec \left( x \right)y} \right)^\prime } & = 2\cos \left( x \right)\sin \left( x \right) - {\sec ^2}\left( x \right)\end{align*}

پس از انتگرال‌گیری از طرفین رابطه بالا، عبارت زیر بدست می‌آید.

(sec(x)y(x))dx=2cos(x)sin(x)sec2(x)dxsec(x)y(x)=sin(2x)sec2(x)dxsec(x)y(x)=12cos(2x)tan(x)+c \large \begin{align*}\int{{{{\left( {\sec \left( x \right)y\left( x \right)} \right)}^\prime }\,dx}} & = \int{{2\cos \left( x \right)\sin \left( x \right) - {{\sec }^2}\left( x \right)\,dx}}\\ \sec \left( x \right)y\left( x \right) & = \int{{\sin \left( {2x} \right) - {{\sec }^2}\left( x \right)\,dx}}\\ \sec \left( x \right)y\left( x \right) & = - \frac{1}{2}\cos \left( {2x} \right) - \tan \left( x \right) + c\end{align*}

توجه داشته باشید که در عبارت فوق، از رابطه مثلثاتی sin(2θ)=2sinθcosθ \sin \left( {2\theta } \right) = 2\sin \theta \cos \theta استفاده شده است. با توجه به عبارت بالا، تابع y برابر می‌شود با:

y(x)=12cos(x)cos(2x)cos(x)tan(x)+ccos(x)=12cos(x)cos(2x)sin(x)+ccos(x) \large \begin{align*}y\left( x \right) & = - \frac{1}{2}\cos \left( x \right)\cos \left( {2x} \right) - \cos \left( x \right)\tan \left( x \right) + c\cos \left( x \right)\\ & = - \frac{1}{2}\cos \left( x \right)\cos \left( {2x} \right) - \sin \left( x \right) + c\cos \left( x \right)\end{align*}

در مرحله بعد، با توجه به شرایط اولیه، ثابتِ c نیز به صورت زیر بدست می‌آید.

32=y(π4)=12cos(π4)cos(π2)sin(π4)+ccos(π4)32=22+c22c=7 \large \begin{align*}3\sqrt 2 = y\left( {\frac{\pi }{4}} \right) & = - \frac{1}{2}\cos \left( {\frac{\pi }{4}} \right)\cos \left( {\frac{\pi }{2}} \right) - \sin \left( {\frac{\pi }{4}} \right) + c\cos \left( {\frac{\pi }{4}} \right)\\ 3\sqrt 2 & = - \frac{{\sqrt 2 }}{2} + c\frac{{\sqrt 2 }}{2}\\ c & = 7\end{align*}

با بدست آمدن c، پاسخ نهایی نیز برابر با تابع زیر بدست می‌آید.

y(x)=12cos(x)cos(2x)sin(x)+7cos(x) \large y\left( x \right) = - \frac{1}{2}\cos \left( x \right)\cos \left( {2x} \right) - \sin \left( x \right) + 7\cos \left( x \right)

نمودار تابعِ بدست آمده، در ادامه ترسیم شده است.

معادله دیفرانسیل خطی

در این مطلب نحوه حل معادله دیفرانسیل خطی مرتبه اول توضیح داده شد. با این حال در آینده روش‌های حل معادلات دیفرانسیل مرتبه بالاتر را نیز توضیح خواهیم داد.

در صورت علاقه‌مندی به مباحث مرتبط در زمینه ریاضی، آموزش‌های زیر نیز به شما پیشنهاد می‌شوند:

^^

فیلم‌ های آموزش معادله دیفرانسیل خطی — از صفر تا صد (+ دانلود فیلم آموزش رایگان)

فیلم آموزشی پاسخ معادله دیفرانسیل خطی

دانلود ویدیو

فیلم آموزشی تعیین عامل انتگرال‌ساز

دانلود ویدیو

فیلم آموزشی حل مثال معادله دیفرانسیل خطی

دانلود ویدیو
بر اساس رای ۵۴ نفر
آیا این مطلب برای شما مفید بود؟
اگر بازخوردی درباره این مطلب دارید یا پرسشی دارید که بدون پاسخ مانده است، آن را از طریق بخش نظرات مطرح کنید.
منابع:
Pauls Online Notes
۳ دیدگاه برای «معادله دیفرانسیل خطی — از صفر تا صد (+ دانلود فیلم آموزش رایگان)»

سلام. وقت بخیر. ممنون بابت سایت خوبتون. در فیلم آموزشی
اول (پاسخ معادله دیفرانسیل خطی) در زمان 6:36 یه خطای کوچک رخ داده:
پاسخ انتگرال (Ln(t برابر میشه با:
(t(Ln(t)-1
که در ویدئو پاسخ به این صورت بیان شده:
(t-1)Ln(t)
که به طبع ادامه پاسخ سوال هم اشتباه شده.

سلامت باشین خیلی اسون و علمی توضیح دادین . بخصوص توضیح مختصر انجام کار تو این بخش. ممنون از مدرسین خوب سایت و محتوای خوبش ?

آقا خدا خیرتون بده من چه در دوره کاردانی و چه در دوره کارشناسی با مطالب شما کارم راه افتاد. متاسفانه تو بقیه سایت های ایرانی (خیلی عذر میخوام) چیزی جز خزعبل نیست مثل آدمم درست و کامل توضیح نمیدن که آدم یه چیزی بفهمه، ولی اینجا راست حسینی همه چی درج میکنین و مثال درست حسابیم میزنین به طوریکه وقتی تو جزوم به مشکل میخورم مطلب شما من یاری میکنه به شدت سپاسگزار و قدردانم بابت راه اندازی این وبلاگ و زحمات شما آقای عوض زاده

نظر شما چیست؟

نشانی ایمیل شما منتشر نخواهد شد. بخش‌های موردنیاز علامت‌گذاری شده‌اند *