نظریه انشعاب | به زبان ساده
در آموزشهای پیشین مجله فرادرس، با نظریه آشوب و فراکتالها آشنا شدیم. در این آموزش نظریه انشعاب را معرفی میکنیم. «نظریه انشعاب» یا «نظریه دوشاخگی» (Bifurcation Theory) آن دسته از سیستمهای دینامیکی را بررسی میکند که تغییرات کوچک در پارامترهای آنها، موجب تغییر ناگهانی در رفتار کیفی یا توپولوژیکیشان میشود.
انشعاب چیست؟
بسیاری از معادلات دیفرانسیل - اگر نگوییم همه آنها - به پارامترها وابسته هستند. بسته به مقادیر این پارامترها، رفتار کیفی جوابهای یک سیستم میتواند کاملاً متفاوت باشد. به عبارت دیگر، میگوییم دستگاه معادله دیفرانسیلِ
$$ \large {\mathbf { \dot x } } ( t) = \mathbf f ( \mathbf x (t) , \mu ) $$
یک انشعاب در مقدار $$ \mu = \mu _ c $$ دارد، اگر در هنگام عبور پارامتر $$ \mu $$ از $$ \mu _ c $$ تغییری ناگهانی در ساختار مسیر موجود به وجود آید. به عبارت دیگر، در انشعاب با تغییر تعداد و یا پایداری نقطه تعادل سیستم موجه خواهیم بود.
چند مثال از نظریه انشعاب
در این بخش، برای درک مفهوم نظریه انشعاب، چند مثال را بررسی میکنیم.
مثال اول نظریه انشعاب
معادله زیر را در نظر بگیرید:
$$ \large \dot { x } = \mu - x ^ 2 . $$
اگر $$ \mu > 0 $$ باشد، دو نقطه تعادل $$ x ^* = \pm \sqrt { \mu } $$ خواهیم داشت. مشتق سمت راست معادله $$ D f ( x , \mu ) = - 2 x $$ است. مقدار این مشتق را در نقاط ثابت به دست میآوریم: $$ D f (\sqrt { \mu} , \mu ) = - 2 \sqrt { \mu } < 0 $$، که نشان میدهد نقطه تعادل $$ x ^ * = - \sqrt { u } $$ پایدار است؛ و $$ D f ( - \sqrt { \mu } , \mu ) = 2 \sqrt { \mu} > 0 $$ که مبین ناپایدار بودن نقطه تعادل $$ x ^ * = - \sqrt { \mu } $$ است. اگر $$ \mu < 0 $$ باشد، تعادلی وجود نخواهد داشت. وقتی $$ \mu = 0 $$ باشد، سیستم تنها یک نقطه تعادل، یعنی $$ x ^ * = 0 $$، دارد. در این حالت، نقطه تعادل «غیرهذلولوی» (Nonhyperbolic) است، زیرا $$ D f ( 0 , 0 ) = 0 $$، و نمیتوانیم از خطیسازی برای تحلیل پایداری آن استفاده کنیم. در این مورد، «رخساره فاز» (Phase Portrait) میتواند راهگشا باشد.
نتیجه میگیریم که نقطه تعادل $$ x ^ * = 0 $$ یک نقطه یا گره زینی ناپایدار است. این سیستم دارای یک انشعاب نقطه زینی در $$ \mu = 0 $$ است (شکل ۲ را ببینید).
مثال دوم نظریه انشعاب
معادله زیر داده شده است:
$$ \large \dot { x } = \mu x - x ^ 2 . $$
این معادله دو نقطه تعادل دارد: $$ x ^ * = 0 $$ و $$ x ^ * = \mu $$. مشتق فضای برداری $$ D f ( x , \mu ) = \mu - 2 x $$ است. با محاسبه این مشتق در $$ x ^ * = 0 $$، مقدار $$ D f ( 0 , \mu ) = \mu $$ را به دست خواهیم آورد که به معنی پایداری نقطه تعادل در صورت $$ \mu < 0 $$ و ناپایداری آن برای $$ \mu > 0 $$ است. همچنین، مقدار $$ D f ( \mu , \mu ) = - \mu $$ را داریم که بیان میکند $$ x ^ * = \mu $$ پایدار است اگر $$ \mu > 0 $$ باشد و ناپایدار است اگر $$ \mu < 0 $$.
اگر $$ \mu = 0 $$، سیستم نقطه تعادل $$ x ^ * = 0 $$ را خواهد داشت. از آنجا که $$ D f ( 0 , 0 ) = 0 $$، نقطه تعادل $$ x ^ * = 0 $$ غیرهذلولوی است. مشابه قبل، از رخساره فاز برای بررسی پایداری استفاده میکنیم. به نظر میرسد که رخساره فاز در این حالت مشابه همانی است که در شکل ۱ دیدیم. بنابراین، وقتی $$ \mu = 0 $$ باشد، نقطه تعادل $$ x ^ * = 0 $$ یک نقطه زینی و در نتیجه ناپایدار است. در این حالت، میگوییم سیستم در مقدار پارامتر $$ \mu = 0 $$ یک «انشعاب تَرابحرانی» (Transcritical Bifurcation) دارد. وقتی پارامتر $$ \mu $$ از مقدار $$ \mu = 0 $$ عبور کند، پایداری تعادل سیستم تغییر میکند.
مثال سوم نظریه انشعاب
معادله زیر را در نظر بگیرید:
$$ \large \dot { x } = \mu x - x ^ 3 . $$
اگر $$ \mu > 0 $$ باشد، این معادله سه نقطه تعادل خواهد داشت: $$ x ^ * = 0$$، $$ x ^ * = \pm \sqrt { \mu} $$. مشتق تابع برابر با $$ D f (x , \mu ) = \mu - 3 x ^ 2 $$ است. بنابراین، $$ D f ( 0 , \mu ) = \mu > 0 $$ است که نتیجه میدهد نقطه ثابت $$ x ^ * = 0 $$ ناپایدار است. همچنین، $$ D f ( \pm \sqrt {\mu } , \mu ) = - 2 \mu < 0 $$، که در نتیجه آن میتوان گفت که هر دو نقطه تعادل دیگر پایدار هستند. اگر $$ \mu < 0 $$ باشد، تنها نقطه تعادل $$ x ^ * = 0 $$ و $$ D f ( 0 , \mu ) = \mu < 0 $$ است، بنابراین، نقطه تعادل پایدار خواهد بود. اگر $$ \mu = 0 $$ باشد، مجدداً، تنها نقطه تعادل $$ x ^ * = 0 $$ است. از آنجا که $$ D f ( 0 , 0 ) = 0 $$، نقطه تعادل غیرهذلولوی است. مشابه مثالهای قبل، رخساره فاز شکل ۴ را خواهیم داشت.
مشاهده میکنیم که نقطه تعادل در این حالت پایدار است. در پارامتر $$ \mu = 0 $$، سیستم در وضعیتی است که «انشعاب چنگالی» (Pitchfork Bifurcation) نامیده میشود (شکل ۵ را ببینید).
انشعاب هوپف
یک «انشعاب هوپف» (Hopf Bifurcation) زمانی رخ میدهد که با تغییر پارامتر $$ \mu $$ یک جواب دورهای (متناوب) یا چرخه حدی در اطراف نقطه تعادل به وجود آید یا از آن دور شود. وقتی یک چرخه حدی پایدار حول نقطه تعادل ناپایدار قرار داشته باشد، انشعاب یک «انشعاب هوپف اَبَربحرانی» (Supercritical Hopf Bifurcation) نامیده میشود. اگر چرخه حدی ناپایدار و حول یک نقطه تعادل پایدار باشد، آنگاه انشعاب یک «انشعاب هوپف زیربحرانی» (Subcritical Hopf Bifurcation) نامیده میشود.
قبل از بیان قضیه، مثالی از یک انشعاب هوپف را در یک سیستم دوبعدی در مختصات قطبی بررسی میکنیم.
مثال اول انشعاب هوپف
سیستم زیر را در نظر بگیرید:
$$ \large \begin {align*}
\dot r & = r ( \mu - r ^ 2 ) , \\
\dot \theta & = -1 , \\
r & \ge 0 .
\end {align*} $$
تنها نقطه بحرانی این سیستم $$ r ^ * = 0 $$، یعنی همان مبدأ است. از آنجا که $$ \dot \theta < 0 $$ است، مسیرها در جهت عقریههای ساعت حول مبدأ حرکت میکنند. اگر $$ \mu = 0 $$ باشد، آنگاه $$ \dot r = - r ^ 3 $$ است. برای $$ r $$ غیرصفر، نامساوی $$ \dot r < 0 $$ را داریم. بنابراین، مسیر بستهای وجود ندارد و وقتی $$ t \to \infty $$ همه مسیرها به مبدأ متمایل میشوند.
مبدأ یک «کانون» (Focus) پایدار است (شکل ۶ را مشاهده کنید). اگر $$ \mu < 0 $$ باشد، آنگاه $$ \mu - r ^ 2 <0 $$ است (برای هر $$r$$).
مشابه مورد قبلی، برای مقادیر غیرصفر $$ r $$ نامساوی $$ \dot r < 0 $$ برقرار است. همچنین، مسیر بستهای وجود ندارد و مبدأ یک کانون پایدار است (شکل ۷).
حال اگر $$ \mu > 0 $$ باشد، آنگاه برای $$ r \in ( \sqrt {\mu} , \infty)$$ نامساوی $$ \dot r < 0 $$ و برای همه مقادیر $$ r $$ نامساوی $$ \mu - r ^ 2 < 0 $$ را داریم. مبدأ یک کانون ناپایدار است و یک مدار یا مسیر بسته پایدار $$ r = \sqrt { \mu} $$ وجود دارد (شکل ۸ را مشاهده کنید).
در این حالت، یک انشعاب هوپف اَبَربحرانی در مقدار پارامتر $$ \mu = 0 $$ رخ داده است (شکل ۹ را ببینید).
قضیه انشعاب هوپف
سیستم «مسطح» (Planar) زیر را در نظر بگیرید:
$$ \large \begin {align*}
\dot x & = f _ \mu ( x , y ) , \\
\dot y & = g _ \mu ( x , y ) ,
\end {align*} $$
که در آن، $$\mu$$ یک پارامتر است. فرض کنید این سیستم یک نقطه ثابت دارد که بدون از دست دادن کلیت مسئله، میتوانیم فرض کنیم در $$ ( x , y ) = ( 0 , 0 )$$ واقع شده است. فرض کنید مقادیر ویژه سیستم خطیشده حول نقطه ثابت $$ \lambda (\mu) , \overline \lambda (\mu) =\alpha (\mu) \pm i \beta (\mu) $$ باشند. همچنین فرض کنید برای مقدار مشخصی از $$ \mu $$ (که فرض میکنیم $$0$$ باشد)، شرایط زیر برقرارند:
- $$ \alpha (0) = 0 $$ و $$ \beta (0) = \omega \neq 0 $$ که $$\text{sgn} (\omega ) = \text{sgn} [\partial g _ \mu / \partial x ) \vert _{\mu = 0 } ( 0 , 0 ) ] $$ (شرایط غیرهذلولوی بودن: مقادیر ویژه موهومی).
- $$ \frac { d \alpha (\mu)} { d \mu} \vert _{\mu = 0 } = d \neq 0 $$ (شرط تقاطعپذیری: مقادیر ویژه محور موهومی را با سرعت غیرصفر قطع میکنند).
- $$ a \neq 0 $$، که $$ a = \frac { 1 } { 1 6 } ( f _ { xxx} + f _ { x yy} + g _ { xxy } + g _{ yyy} ) + \frac { 1 } { 16 \omega } ( f _ { x y } ( f _ { xx} + f _ {yy} ) - g _ { x y } ( g _ { yx } ( g _ {xx} + g _ { yy}) - f _ { xx} g _ { xx} + f _ { yy} g _ { yy} $$ با $$ f _ { x y } = ( \partial ^ 2 f _ \mu / \partial _ x \partial _ y ) \vert _ { \mu = 0} ( 0 , 0 ) $$ و غیره (شرط عمومیت).
در نتیجه، اگر $$ a d < 0 $$ یا $$ a d > 0 $$ باشد، به ترتیب، یک منحنی یکتا از جوابهاب دورهای از مبدأ به ناحیه $$ \mu > 0 $$ یا $$ \mu < 0 $$ انشعاب ایجاد میکنند. به ازای $$ \mu > 0 $$، مبدأ یک نقطه ثابت پایدار است و برای $$ \mu < 0 $$ یک نقطه ثابت ناپایدار خواهد بود، در حالی که اگر مبدأ در $$ \mu = 0 $$ (که جوابهای دورهای وجود دارند) ناپایدار باشد، جوابهای دورهای پایدار هستند.
دامنه مسیرهای بسته دورهای تا $$ \sqrt { |\mu|} $$ افزایش مییابد، در حالی که دوره تناوب آنها وقتی $$ |\mu | $$ به صفر میل میکند، تا $$ 2 \pi / |\omega |$$ هم میرسد.
مثال دوم انشعاب هوپف
«معادله لینارد» (Liénard Equation) را در نظر بگیرید:
$$ \large \ddot x - ( \mu - x ^ 2 ) \dot x + x = 0 . $$
اگر $$ u = x $$ و $$ v = \dot x $$ را در نظر بگیریم، معادله را میتوان به صورت یک دستگاه معادلات مرتبه اول دوبعدی نوشت:
$$ \large \begin {align*}
\dot u & = v , \\
\dot v & = - u + ( \mu - u ^ 2 ) v .
\end {align*} $$
تنها نقطه تعادل مبدأ است. ماتریس ژاکوبی سیستم خطیشده حول مبدأ به صورت زیر است:
$$ \large \left ( \begin {array} { c } 0 & 1 \\ -1 & \mu \end {array} \right ). $$
مقادیر ویژه ماتریس ژاکوبی $$ \alpha ( \mu ) \pm \beta ( \mu) = \mu / 2 \pm i \sqrt { 4 - \mu ^ 2 } / 2 $$ هستند. توجه کنید که $$ \alpha ( 0 ) = 0 $$ و $$ \omega = \beta ( 0 ) = - 1 \neq 0 $$ است. همچنین، $$ d = \frac { d \alpha ( \mu)} { d \mu } \vert _ { \mu = 0 } = 1 / 2 \neq 0 $$ است. در نهایت، $$ a = - \frac 18 \neq 0 $$. بنابراین، همه شرایط قضیه انشعاب هوپف برقرار است. در نتیجه، $$ a d = - 1 / 16 < 0 $$ و مبدأ برای $$ \mu < 0 $$ پایدار است (شکل ۱۰ را ببینید) و برای $$ \mu > 0 $$ ناپایدار است، در حالی که یک مسیر بسته دورهای پایدار وجود دارد (شکل ۱۱). سیستم دارای یک انشعاب هوپف ابربحرانی در $$ \mu = 0 $$ است.
مثال سوم انشعاب هوپف
«مدل شلکوف» (Sel’kov Model) گلیکولیز یا قندکافت را که فرایند از بین بردن قند با سلولهای زنده برای به دست آوردن انرژی است، در نظر بگیرید:
$$ \large \begin {align*}
\dot x & = - x + a y + x ^ 2 \\
\dot y & = b - a y - x ^ 2 y ,
\end {align*} $$
که $$ x $$ و $$ y $$، به ترتیب، غلظت ADP و F6P هستند و $$ a , b > 0 $$. نقطه تعادل سیستم به صورت زیر است:
$$ \large x ^ * = b , \; y ^ * = \frac { b } { a + b ^ 2 } . $$
اگر مختصات را با قرار دادن $$ \widetilde { x } = x - b $$ و $$ \widetilde { y } = y - \frac { b } { a + b ^ 2 } $$ تغییر دهیم، آنگاه $$ x = \widetilde { x } + b $$، $$ y = \widetilde { y } + \frac { b } { a + b ^ 2 } $$، $$ \dot{ \widetilde{x}} = \dot x $$ و $$ \dot {\widetilde { {y } }} = \dot y$$ را خواهیم داشت. آنگاه اگر $$ \widetilde {x} = x $$ و $$ \widetilde {y} = y $$ را تغییر نام دهیم، سیستم به صورت زیر درخواهد آمد:
$$ \large \begin {align*}
\dot { x } & = - ( x + b ) + a \left ( y + \frac { b } { a + b ^ { 2 } } \right ) + ( x + b ) ^ { 2 } \left ( y + \frac { b } { a + b ^ { 2 } } \right ) \\
\dot { y } & = b - a \left ( y + \frac { b } { a + b ^ { 2 } } \right ) - ( x + b ) ^ { 2 } \left ( y + \frac { b } { a + b ^ { 2 } } \right )
\end {align*} $$
که نقطه تعادل آن $$ ( 0 , 0 ) $$ است.
ماتریس ژاکوبی این سیستم در نقطه تعادل به صورت زیر است:
$$ \large \left ( \begin {array} { c c }
- 1 + \frac { 2 b ^ { 2 } } { a + b ^ { 2 } } & a + b ^ { 2 } \\
- \frac { 2 b ^ { 2 } } { a + b ^ { 2 } } & - a - b ^ { 2}
\end {array} \right) $$
مقادیر ویژه ماتریس ژاکوبی به صورت زیر هستند:
$$ \large \lambda _ { \pm } ( a , b ) = \frac { a + a ^ { 2 } -b ^ { 2 } + 2 a b ^ { 2 } + b ^ { 4 } \pm i \sqrt { 4 \left ( a + b ^ { 2 } \right ) ^ { 3 } - \left ( a ( 1 +a ) + ( - 1 + 2 a ) b ^ { 2 } + b ^ { 4 } \right ) ^ { 2 } }} {- 2 \left ( a + b ^ { 2 } \right ) } $$
داریم:
$$ \large \alpha ( a , b ) = R e \left ( \lambda _ { \pm } ( a , b ) \right ) = \frac { a + a ^ { 2 } - b ^ {2 } + 2 a b ^ { 2 } + b ^ { 4 } } { - 2 \left ( a + b ^ { 2 } \right ) } $$
و
$$ \large \beta ( a , b ) = \operatorname {Im} \left ( \lambda _ { \pm } ( a , b ) \right ) = \frac { \sqrt { 4 \left ( a + b ^ { 2 } \right ) ^ { 3 } - \left ( a ( 1 + a ) + ( - 1 + 2 a ) b ^ { 2 } + b ^ { 4 } \right ) ^ { 2 } } } { - 2 \left ( a + b ^ { 2 } \right ) } $$
برای برقراری شرط ۱ قضیه انشعاب هوپف، باید داشته باشیم:
$$ \large b = b _ { 1 } ( a ) = \sqrt { ( 1 - \sqrt { 1 - 8 a } - 2 a ) / 2 } $$
یا
$$ \large b = b _ { 2 } ( a ) = \sqrt { ( 1 + \sqrt { 1 - 8 a } - 2 a ) / 2 } $$
مشتق بخش حقیقی مقادیر ویژه نسبت به پارامتر $$ b $$ در هریک از این مقادیر، به صورت زیر است:
$$ \large \left . \frac { d \alpha ( a , b ) } { d b } \right | _ { b = b_ { 1 } ( a ) } = \frac { \sqrt { 2 - 1 6 a } \sqrt { 1 -\sqrt { 1 - 8 a } - 2 a } } { 1 - \sqrt { 1 - 8 a } } , $$
و
$$ \large \left . \frac { d \alpha ( a , b ) } { d b } \right | _ { b =b _ { 2 } ( a ) } = \frac { \sqrt { 2 - 1 6 a } \sqrt { 1 + \sqrt { 1 - 8 a } - 2 a } } { - 1 - \sqrt { 1 - 8 a } } $$
عبارت کمّی شرط ۳ قضیه انشعاب هوپف به صورت زیر است:
$$ \large \left. \left ( - \frac { 1 } { 8 } -\frac { ( a + b ) ^ { 2 } \left ( \frac { b } { a + b ^ { 2 } } + y \right ) ^ { 2 } } { 2 \sqrt { 4 \left ( a + b ^ { 2 } \right ) ^ { 3 } - \left ( a ( 1 + a) + ( - 1 + 2 a ) b ^ { 2 } + b ^ {4 } \right ) ^ { 2 } } } \right ) \right | _ { b = b _ { 1 } , b _ { 2 } } ( 0 , 0 ) . $$
اکنون مقدار پارامتر $$ \alpha = 0.1 $$ را تصحیح میکنیم. در نتیجه، $$ b_{1}(a)=b_{1}(0.1) \approx 0.419992 $$ و $$ b_{2}(a)=b_{2}(0.1) \approx 0.789688 $$ را خواهیم داشت. همچنین، $$\omega_{1}=\beta\left(0.1, b_{1}(0.1)\right) \approx-0.525731$$، $$\omega_{2}=\beta\left(0.1, b_{2}(0.1)\right) \approx-0.850651$$، $$d_{1}=\left.\frac{d \alpha(0.1, b)}{d b}\right|_{b=b_{1}(0.1)} \approx 0.679561$$ و $$ d_{2}=\left.\frac{d \alpha(0.1, b)}{d b}\right|_{b=b_{2}(0.1)} \approx-0.488054 $$ را داریم. در نهایت، مقادیر عبارت شرط ۳ قضیه $$a_{1} \approx-1.223 $$ و $$ a_{2} \approx-0.475021 $$ هستند.
بنابراین، برای هر دو مقدار پارامتر $$ b $$، قضیه انشعاب هوپف اعمال میشود. در هریک از دو حالت، $$ b = b _ 1 \approx 0.419992 $$ و $$ b = b _ 2 \approx 0.789688 $$، یک انشعاب هوپف وجود دارد. دو نامساوی $$ d _ 1 > 0 $$ و $$ a _ 1 d _ 1 < 0 $$ را داریم. بنابراین، مبدأ برای $$ b < b _ 1 $$ پایدار است و یک چرخه حدی پایدار برای $$ b > b _ 1 $$ وجود دارد که $$ b $$ به اندازه کافی کوچک است. علاوه بر این، $$ d _ 2 < 0 $$ و $$ a _ 2 d _ 2 > 0 $$ است. در نتیجه، مبدأ در ناحیه $$ b _ 1 < b < b _ 2 $$ ناپایدار است. در این ناحیه یک چرخه حدی یکتا و پایدار وجود دارد (شکلهای ۱۲، ۱۳ و ۱۴ را ببینید).
جمعبندی
نظریه انشعاب یک زمینه پژوهشی گسترده و فعال است. در این آموزش، با چند مثال ساده، برخی از انواع مختلف انشعابها را در دستگاه معادلات دیفرانسیل معمولی یکبعدی و دوبعدی بررسی کردیم. همچنین نمونههای انشعاب را از طریق رخساره فاز، نمودارهای انشعاب و قضیه انشعاب هوپف مورد تجزیه و تحلیل قرار دادیم.
سلام، خیلی ممنون که مطالب ریاضی هم قرار دادید. متن خیلی خوبی داره و واضح بیان شده. امیدوارم بیشتر از این مطالب قرار بدید. واقعا متشکرم.