خريطة لوجستية

يُعرض سلوك الخريطة اللوجستية في شكل مخطط شبكة العنكبوت . يُظهر الرسم المتحرك تغير السلوك مع زيادة قيمة المعامل ( r في الشكل) من 1 إلى 4، بدءًا من قيمة ابتدائية للدالة تبلغ 0.2.

الخريطة اللوجستية هي نظام ديناميكي منفصل يُعرَّف بواسطة معادلة الفرق التربيعية

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

استُخدمت الخريطة في البداية من قِبل إدوارد لورنز في ستينيات القرن العشرين لعرض خصائص الحلول غير المنتظمة في الأنظمة المناخية. [ 1 ] وشاع استخدامها في ورقة بحثية نُشرت عام 1976 من قِبل عالم الأحياء روبرت ماي ، [ ماي، روبرت م. (1976) 1 ] ، جزئيًا كنموذج ديموغرافي منفصل زمنيًا يُشابه المعادلة اللوجستية التي وضعها بيير فرانسوا فيرهولست . [ 2 ] ومن بين الباحثين الآخرين الذين ساهموا في دراسة الخريطة اللوجستية: ستانيسواف أولام ، وجون فون نيومان ، وبيكا ميربيرغ ، وأولكسندر شاركوفسكي ، ونيكولاس متروبوليس ، وميتشل فيجنباوم . [ 3 ]

مثالان تمهيديان

مثال على الأنظمة الديناميكية

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

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

بفرض قيمة مناسبة للمعامل r وإجراء الحسابات بدءًا من حالة ابتدائيةx0{\displaystyle x_{0}}، فنحصل على المتتاليةx0{\displaystyle x_{0}}،x1{\displaystyle x_{1}}،x2{\displaystyle x_{2}}...، والتي يمكن تفسيرها على أنها سلسلة من الخطوات الزمنية في تطور النظام.

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

هناك طريقة أخرى لفهم هذا التسلسل وهي تكرار الخريطة اللوجستية (الممثلة هنا بواسطةو(x){\displaystyle f(x)}) إلى الحالة الأوليةx0{\displaystyle x_{0}}: [ ديفاني 1989 2 ]x1=و(x0)،x2=و(x1)=و(و(x0))،x3=و(x2)=و(و(و(x0)))،x4=...{\displaystyle {\begin{aligned}x_{1}&=f(x_{0}),\\x_{2}&=f(x_{1})=f(f(x_{0})),\\x_{3}&=f(x_{2})=f(f(f(x_{0}))),\\x_{4}&=\dots \\\end{aligned}}}

كان هذا هو النهج الأولي لهنري بوانكاريه لدراسة الأنظمة الديناميكية ، وفي نهاية المطاف الفوضى، انطلاقًا من دراسة النقاط الثابتة ، أو بعبارة أخرى، الحالات التي لا تتغير بمرور الوقت (أي عندماxن=...=x1=x0=و(x0){\displaystyle x_{n}=...=x_{1}=x_{0}=f(x_{0})}تنشأ العديد من الأنظمة الفوضوية، مثل مجموعة ماندلبروت، من تكرار دوال تربيعية غير خطية بسيطة للغاية مثل الخريطة اللوجستية. [ 5 ]

مثال على نموذج ديموغرافي

بأخذ نموذج السكان البيولوجي كمثال، فإن x<sub> n </sub> هو عدد بين الصفر والواحد، ويمثل نسبة السكان الحاليين إلى الحد الأقصى الممكن للسكان . [ ماي، روبرت م. (1976) 2 ] تهدف معادلة الفرق غير الخطية هذه إلى رصد تأثيرين:

  • التكاثر ، حيث يزداد عدد السكان بمعدل يتناسب مع عدد السكان الحالي عندما يكون حجم السكان صغيراً.
  • المجاعة (الوفيات المرتبطة بالكثافة)، حيث سينخفض ​​معدل النمو بمعدل يتناسب مع القيمة التي تم الحصول عليها من خلال طرح عدد السكان الحالي من "القدرة الاستيعابية" النظرية للبيئة.

القيم المعتادة للمعامل r هي تلك الموجودة في الفترة [0،  4] ، بحيث تبقى x<sub> n</sub> محصورة في الفترة [0،  1] . حالة r = 4 للخريطة اللوجستية هي تحويل غير خطي لكل من خريطة إزاحة البت وحالة μ = 2 لخريطة الخيمة . إذا كانت r > 4 ، فإن هذا يؤدي إلى أحجام سكانية سالبة. (لا تظهر هذه المشكلة في نموذج ريكر القديم ، الذي يُظهر أيضًا ديناميكيات فوضوية). يمكن أيضًا النظر في قيم r في الفترة [−2،  0] ، بحيث تبقى x<sub> n</sub> محصورة في الفترة [−0.5،  1.5] . [ 6 ]

توصيف الخريطة اللوجستية

بالنسبة لسلوك التسلسل من r=0.02 إلى r=4، يمكن تصور الإحداثي الأفقي على أنه الزمن، والإحداثي الرأسي إما كموقع في الفضاء عند الزمن t أو كحجم السكان عند الزمن t

يُظهر الرسم المتحرك سلوك التسلسلxن{\displaystyle x_{n}}على قيم مختلفة للمعامل r. الملاحظة الأولى هي أن المتتالية لا تتباعد وتبقى محدودة لقيم r بين 0 و4. من الممكن ملاحظة الظواهر النوعية التالية بترتيب زمني:

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

كما ذُكر سابقًا، فإن الخريطة اللوجستية هي دالة تربيعية عادية. ومن الأسئلة المهمة في مجال الأنظمة الديناميكية: كيف يتغير سلوك المسار بتغير قيمة المعامل r ؟ وبحسب قيمة r ، قد يكون سلوك مسار الخريطة اللوجستية بسيطًا أو معقدًا. [ Thompson & Stewart 1 ] فيما يلي، سنشرح كيف يتغير سلوك الخريطة اللوجستية مع ازدياد قيمة r .

المجال، والرسوم البيانية، والنقاط الثابتة

رسم بياني للخريطة اللوجستية (العلاقة بينxن+1{\displaystyle x_{n+1}}وxن{\displaystyle x_{n}}). يكون شكل الرسم البياني على شكل قطع مكافئ، ويتغير رأس القطع المكافئ بتغير المعامل r.

كما ذُكر سابقًا، يمكن استخدام الخريطة اللوجستية كنموذج لدراسة تقلبات حجم الجماعة. في هذه الحالة، يُمثل المتغير x في الخريطة اللوجستية عدد أفراد الكائن الحي مقسومًا على الحد الأقصى لحجم الجماعة، وبالتالي فإن القيم الممكنة لـ x محصورة بين 0 و1. لهذا السبب، غالبًا ما يُناقش سلوك الخريطة اللوجستية بحصر نطاق المتغير في الفترة [0، 1]. [ Hirsch, Smale & Devaney 1 ]

إذا قصرنا المتغيرات على النطاق 0 ≤ x ≤ 1، فإن نطاق المعامل r سيقتصر بالضرورة على النطاق من 0 إلى 4 (0 ≤ r ≤ 4). وذلك لأنه إذاxن{\displaystyle x_{n}}إذا كانت القيمة ضمن النطاق [0، 1]، فإن القيمة القصوى لـxن+1{\displaystyle x_{n+1}}r/4. وبالتالي، عندما r > 4، فإن قيمةxن+1{\displaystyle x_{n+1}}يمكن أن تتجاوز القيمة 1. من ناحية أخرى، عندما تكون قيمة r سالبة، يمكن أن تأخذ قيمة x قيمًا سالبة. [ Hirsch, Smale & Devaney 1 ]

يمكن أيضًا استخدام رسم بياني للخريطة لمعرفة الكثير عن سلوكها. الرسم البياني للخريطة اللوجستيةxن+1=رxن(1-xن){\displaystyle x_{n+1}=rx_{n}(1-x_{n})}هو المنحنى المستوي الذي يرسم العلاقة بينxن{\displaystyle x_{n}}وxن+1{\displaystyle x_{n+1}}، معxن{\displaystyle x_{n}}(أو x) على المحور الأفقي وxن+1{\displaystyle x_{n+1}}(أو f (x)) على المحور الرأسي. يبدو الرسم البياني للخريطة اللوجستية كما يلي، باستثناء الحالة r = 0:

له شكل قطع مكافئ برأس عند [ جوليك 1 ]

عند تغيير قيمة r، يتحرك رأس القطع المكافئ لأعلى أو لأسفل، ويتغير شكل القطع المكافئ . بالإضافة إلى ذلك، يتقاطع القطع المكافئ للخريطة اللوجستية مع المحور الأفقي (الخط الذي عنده r = 0).xن+1=0{\displaystyle x_{n+1}=0}) عند نقطتين . نقطتا التقاطع هما(xن،xن+1)=(0،0){\displaystyle (x_{n},x_{n+1})=(0,0)}و(xن،xن+1)=(1،0){\displaystyle (x_{n},x_{n+1})=(1,0)}، ومواقع نقاط التقاطع هذه ثابتة ولا تعتمد على قيمة r.

مثال على إسقاط شبكة العنكبوت لمسار على الرسم البياني للخريطة اللوجستية، ومواقع النقاط الثابتةxو1{\displaystyle x_{f1}}وxو2{\displaystyle x_{f2}}على الرسم البياني.

تُعدّ الرسوم البيانية للخرائط، وخاصةً تلك التي تمثل متغيرًا واحدًا مثل الخريطة اللوجستية، أساسيةً لفهم سلوك الخريطة. ومن استخدامات الرسوم البيانية توضيح النقاط الثابتة. ارسم خطًا مستقيمًا y = x (خطًا بزاوية 45 درجة) على الرسم البياني للخريطة. إذا وُجدت نقطة يتقاطع فيها هذا الخط مع الرسم البياني، فإن تلك النقطة تُعتبر نقطة ثابتة. رياضيًا، النقطة الثابتة هي

يعني ذلك نقطة لا تتغير عند تطبيق الخريطة. سنرمز إلى النقطة الثابتة بـxو{\displaystyle x_{f}}في حالة الخريطة اللوجستية، يتم الحصول على النقطة الثابتة التي تحقق المعادلة (2-2) عن طريق حلرx(1-x)=x{\displaystyle rx(1-x)=x}.

(باستثناء r = 0). يُعد مفهوم النقاط الثابتة ذا أهمية أساسية في الأنظمة الديناميكية المنفصلة.

هناك تقنية بيانية أخرى يمكن استخدامها لرسم الخرائط ذات المتغير الواحد وهي إسقاط شبكة العنكبوت . بعد تحديد قيمة أوليةx0{\displaystyle x_{0}}على المحور الأفقي، ارسم خطًا رأسيًا من القيمة الأوليةx0{\displaystyle x_{0}}ارسم خطًا أفقيًا من نقطة تقاطع منحنى f(x) مع الخط y = x بزاوية 45 درجة، ثم ارسم خطًا رأسيًا من نقطة التقاطع إلى منحنى f(x). بتكرار هذه العملية، يتشكل مخطط يشبه شبكة العنكبوت أو الدرج على المستوى. هذا الرسم البياني يُكافئ حساب المسار بيانيًا، ويمثل مخطط شبكة العنكبوت المسار بدءًا منx0{\displaystyle x_{0}}. يسمح هذا الإسقاط برؤية السلوك العام للمسار بنظرة سريعة.

يعتمد السلوك على r

تُظهر الصورة أدناه سعة ومحتوى التردد لخريطة لوجستية تُكرر نفسها لقيم معلمات تتراوح من 2  إلى  4. وكما هو الحال، يمكن ملاحظة سلوك خطي أولي، ثم سلوك فوضوي، ليس فقط في المجال الزمني (يسارًا)، بل وخاصة في مجال التردد أو الطيف (يمينًا)، أي أن الفوضى موجودة على جميع المقاييس كما هو الحال في شلال طاقة كولموغوروف ، بل وتنتشر من مقياس إلى آخر. [ Thompson & Stewart 2 ]

بتغيير المعامل r ، يُلاحظ السلوك التالي:

الحالة عندما يكون 0 ≤ r < 1

أولاً، عندما تكون قيمة المعامل r تساوي صفرًا،x1=0{\displaystyle x_{1}=0}بغض النظر عن القيمة الأوليةx0{\displaystyle x_{0}}بمعنى آخر، فإن مسار الخريطة اللوجستية عندما تكون a = 0 هو مسار تكون فيه جميع القيم بعد القيمة الأولية تساوي 0، لذلك لا يوجد الكثير للبحث فيه في هذه الحالة.

بعد ذلك، عندما تكون قيمة المعامل r في النطاق 0 < r < 1،xن{\displaystyle x_{n}}يتناقص بشكل رتيب لأي قيمة منx0{\displaystyle x_{0}}بين 0 و 1. أي،xن{\displaystyle x_{n}}يتقارب إلى 0 في حالة n → ∞. [ جوليك 2 ] النقطة التيxن{\displaystyle x_{n}}النقطة الثابتة هي نقطة التقاربxو1{\displaystyle x_{f1}}كما هو موضح في المعادلة (2-3). تُسمى النقاط الثابتة من هذا النوع، حيث تتقارب المدارات حولها، مستقرة تقاربياً، أو مستقرة، أو جاذبة. وعلى العكس من ذلك، إذا كانت المدارات حولxو{\displaystyle x_{f}}الابتعاد عنxو{\displaystyle x_{f}}مع ازدياد الزمن n، النقطة الثابتةxو{\displaystyle x_{f}}يُطلق عليه اسم غير مستقر أو تنافري. [ جوليك 3 ]

مخطط العنكبوت (يسار) والسلسلة الزمنية (n مقابل xn) (يمين) للمعامل r = 0.9. يتقارب المسار بشكل رتيب إلى 0.

إحدى الطرق الشائعة والبسيطة لمعرفة ما إذا كانت نقطة ثابتة مستقرة تقاربياً هي حساب مشتقة الدالة f. [ جوليك 4 ] تُعبّر هذه المشتقة على النحو التالي:و(x){\displaystyle f'(x)}، xو{\displaystyle x_{f}}تكون مستقرة تقاربياً إذا تحقق الشرط التالي.

ميول المماس لنقطة ثابتة مستقرة تقاربياً (يسار) ونقطة ثابتة غير مستقرة (يمين) وحالة المدارات المحيطة بها

يمكن ملاحظة ذلك من خلال رسم الخريطة: إذا كان ميل المماس للمنحنى عندxو{\displaystyle x_{f}}إذا كانت القيمة بين -1 و 1، فإنxو{\displaystyle x_{f}}مستقر، والمدار حوله ينجذب إلىxو{\displaystyle x_{f}}مشتق الخريطة اللوجستية هو

لذلك، بالنسبة لـ x = 0 و 0 < r < 1، فإن 0 < f '(0) < 1، وبالتالي النقطة الثابتةxو1{\displaystyle x_{f1}}= 0 يحقق المعادلة (3-1).

ومع ذلك، فإن طريقة التمييز باستخدام المعادلة (3-1) لا تعرف نطاق المدارات منxو{\displaystyle x_{f}}الذين ينجذبون إلىxو{\displaystyle x_{f}}يضمن ذلك فقط أن x ضمن نطاق معين منxو{\displaystyle x_{f}}سوف تتقارب. في هذه الحالة، يكون نطاق القيم الأولية التي تتقارب إلى 0 هو النطاق الكامل [0، 1]، ولكن للتأكد من ذلك، يلزم إجراء دراسة منفصلة.

يمكن إيجاد طريقة لتحديد ما إذا كانت النقطة الثابتة غير مستقرة عن طريق اشتقاق الدالة بنفس الطريقة. بالنسبة لـ r<1، إذا كانت النقطة الثابتةxو{\displaystyle x_{f}}غير مستقر إذا

إذا كانت قيمة المعامل تقع في النطاق 0 < r < 1، فإن النقطة الثابتة الأخرىxو2=1-1/أ{\displaystyle x_{f2}=1-1/a} إنها سالبة وبالتالي لا تقع في النطاق [0، 1]، ولكنها موجودة كنقطة ثابتة غير مستقرة.

الحالة عندما يكون 1 ≤ r ≤ 2

في الحالة العامة عندما تكون قيمة r بين 1 و 2، سيقترب عدد السكان بسرعة من القيمة r − 1 / r ، بغض النظر عن عدد السكان الأولي.

يحدث التشعب العابر للخريطة اللوجستية عند r = 1. أما بالنسبة لـ r < 1،xو2{\displaystyle x_{f2}}توجد خارج [0، 1] كنقطة ثابتة غير مستقرة، ولكن بالنسبة لـ r = 1، تتصادم النقطتان الثابتتان، وبالنسبة لـ r > 1،xو2{\displaystyle x_{f2}}يظهر بين [0، 1] كنقطة ثابتة مستقرة.

عندما تكون قيمة المعامل r = 1، يتقارب مسار الخريطة اللوجستية إلى الصفر كما في السابق، ولكن سرعة التقارب أبطأ عند r = 1. النقطة الثابتة 0 عند r = 1 مستقرة تقاربياً، لكنها لا تحقق المعادلة (3-1). في الواقع، تعمل طريقة التمييز القائمة على المعادلة (3-1) بتقريب الخريطة إلى الرتبة الأولى بالقرب من النقطة الثابتة. عندما تكون r = 1، لا يصح هذا التقريب، ويُحدد الاستقرار أو عدم الاستقرار بواسطة الحدود التربيعية للخريطة، أو بعبارة أخرى، بواسطة الاضطراب من الرتبة الثانية.

عند رسم المنحنى عند r = 1، يكون المنحنى مماسًا للقطر المائل بزاوية 45 درجة عند x = 0. في هذه الحالة، تكون النقطة الثابتةxو2=1-1/ر{\displaystyle x_{f2}=1-1/r}، والتي توجد في النطاق السالب لـ0<ر<1{\displaystyle 0<r<1}، يكون xو2=0{\displaystyle x_{f2}=0}. لxو2=0{\displaystyle x_{f2}=0}أي أنه مع ازدياد قيمة r، تزداد قيمةxو2{\displaystyle x_{f2}} يقترب من الصفر، وعند r = 1 تحديداً،xو2{\displaystyle x_{f2}}يصطدم بـ xو1=0{\displaystyle x_{f1}=0}يؤدي هذا التصادم إلى ظاهرة تُعرف باسم التشعب العابر . يُستخدم مصطلح التشعب لوصف تغير نوعي في سلوك نظام ديناميكي. في هذه الحالة، يحدث التشعب العابر عندما تتناوب استقرارية النقاط الثابتة فيما بينها . أي عندما تكون قيمة r أقل من 1 .xو1{\displaystyle x_{f1}}مستقر وxو2{\displaystyle x_{f2}}غير مستقر، ولكن عندما تكون قيمة r أكبر من 1،xو1{\displaystyle x_{f1}}غير مستقر وxو2{\displaystyle x_{f2}}مستقر. تُسمى قيم المعاملات التي يحدث عندها التشعب بنقاط التشعب. في هذه الحالة، r = 1 هي نقطة التشعب.

نقطة ثابتةxو2=1-1/ر{\displaystyle x_{f2}=1-1/r}مثال على التقارب المتناقص بشكل رتيب إلى (r = 1.2، x 0 = 0.6)
نقطة ثابتةxو2=1-1/أ{\displaystyle x_{f2}=1-1/a}مثال على التقارب المتزايد بشكل رتيب إلى (r = 1.8، x 0 = 0.2)

نتيجةً للتفرع، يتقارب مدار الخريطة اللوجستية إلى نقطة النهاية.xو2=1-1/ر{\displaystyle x_{f2}=1-1/r}بدلاً منxو1=0{\displaystyle x_{f1}=0}. على وجه الخصوص، إذا كانت المعلمة1<ر2{\displaystyle 1<r\leq 2}ثم المسار الذي يبدأ من قيمةx0{\displaystyle x_{0}}في الفترة (0، 1)، باستثناء 0 و1، تتقارب إلىxو2{\displaystyle x_{f2}}عن طريق الزيادة أو النقصان بشكل رتيب. ويعتمد اختلاف نمط التقارب على نطاق القيمة الأولية. 0<x0<1-1/ر{\displaystyle 0<x_{0}<1-1/r}

في حالة1-1/ر<x0<1/ر{\displaystyle 1-1/r<x_{0}<1/r} ثم يتقارب بشكل رتيب، 1/ر<x0<1{\displaystyle 1/r<x_{0}<1}، تتقارب الدالة بشكل رتيب باستثناء الخطوة الأولى.

علاوة على ذلك، النقطة الثابتةxو1=0{\displaystyle x_{f1}=0}يصبح غير مستقر بسبب التشعب، ولكنه يستمر في الوجود كنقطة ثابتة حتى بعد أن تصبح قيمة r أكبر من 1. هذا لا يعني أنه لا توجد قيمة ابتدائية أخرى غير r.xو1{\displaystyle x_{f1}}التي يمكنها الوصول إلى هذه النقطة الثابتة غير المستقرة xو1{\displaystyle x_{f1}}. هذا هوx0=1{\displaystyle x_{0}=1}وبما أن الخريطة اللوجستية تحقق الشرط f(1) = 0 بغض النظر عن قيمة r، فإن تطبيق الخريطة مرة واحدة علىx0=1{\displaystyle x_{0}=1}يرسمها إلىxو1=0{\displaystyle x_{f1}=0}. النقطة مثل x = 1 التي يمكن الوصول إليها مباشرة كنقطة ثابتة من خلال عدد محدود من تكرارات الخريطة تسمى نقطة ثابتة نهائية.

الحالة عندما يكون 2 ≤ r ≤ 3

عندما تكون قيمة r بين 2 و3، سيقترب عدد السكان في النهاية من القيمة نفسها r − 1 / r ، ولكنه سيتذبذب حول هذه القيمة لفترة من الزمن. معدل التقارب خطي، باستثناء r = 3 ، حيث يكون بطيئًا للغاية، وأقل من الخطي (انظر ذاكرة التشعب ).

عندما تكون قيمة المعامل 2 < r < 3، باستثناء القيمتين الابتدائيتين 0 و1، تكون النقطة الثابتةxو2=1-1/ر{\displaystyle x_{f2}=1-1/r}يكون الوضع مماثلاً لما هو عليه عندما يكون 1 < r ≤ 2. ومع ذلك، في هذه الحالة، لا يكون التقارب رتيبًا. فعندما يقترب المتغير منxو2{\displaystyle x_{f2}}يصبح أكبر وأصغر منxو2{\displaystyle x_{f2}}ويتكرر ذلك بشكل متكرر، ويتبع مسارًا متقاربًا يتذبذب حولxو2{\displaystyle x_{f2}}.

القيمة التي يتم تعيينها إلىxو2{\displaystyle x_{f2}}بتطبيق عملية الربط مرة واحدةو(x~و2)=xو2{\displaystyle f({\tilde {x}}_{f2})=x_{f2}} -->

رسم متحرك لإسقاط العنكبوت عند a = 2.8، يتقارب حول نقطة ثابتة.

بشكل عام، تُعدّ مخططات التشعب مفيدةً لفهم التشعبات. هذه المخططات عبارة عن رسوم بيانية للنقاط الثابتة (أو النقاط الدورية، كما هو موضح أدناه) x كدالة للمعامل a، حيث يُمثّل a المحور الأفقي و x المحور الرأسي. وللتمييز بين النقاط الثابتة المستقرة وغير المستقرة، تُرسَم المنحنيات الأولى أحيانًا كخطوط متصلة، بينما تُرسَم المنحنيات الثانية كخطوط متقطعة. عند رسم مخطط تشعب للخريطة اللوجستية، يكون لدينا خط مستقيم يُمثّل النقطة الثابتة.xو1=0{\displaystyle x_{f1}=0} وخط مستقيم يمثل النقطة الثابتةxو2=1-1/أ{\displaystyle x_{f2}=1-1/a} يمكن ملاحظة أن المنحنيات التي تمثل a و b تتقاطع عند r = 1، وأن الاستقرار يتم تبديله بين الاثنين.

مخطط تشعب الخريطة اللوجستية للمعاملات من 0 إلى 3. يمثل الخط الأزرق النقطة الثابتةxو1=0{\displaystyle x_{f1}=0}والخط الأحمر يمثل النقطة الثابتةxو2=1-1/ر{\displaystyle x_{f2}=1-1/r}يمثل.

الحالة عندما يكون 3 ≤ r ≤ 3.44949

في الحالة العامة، عندما تكون قيمة r بين 3 و 1  + √6 ≈ 3.44949، فإن عدد السكان سيقترب من تذبذبات دائمة بين قيمتين. هاتان القيمتان تعتمدان على r ويتم تحديدهما بالمعادلة [ 6 ]. x±=12ر(ر+1±(ر-3)(ر+1)){\displaystyle x_{\pm }={\frac {1}{2r}}\left(r+1\pm {\sqrt {(r-3)(r+1)}}\right)}.

عندما تكون قيمة المعامل r = 3 بالضبط، فإن المدار يحتوي أيضًا على نقطة ثابتةxو2=1-1/ر{\displaystyle x_{f2}=1-1/r}ومع ذلك ، تتقارب المتغيرات بشكل أبطأ مما كانت عليه عندما2<ر<3{\displaystyle 2<r<3}. متىر=3{\displaystyle r=3}المشتقو(xو2){\displaystyle f'(x_{f2})}تصل قيمة r إلى -1 ولا تعود تحقق المعادلة (3-1). عندما تتجاوز r القيمة 3،و(xو2)<-1{\displaystyle f'(x_{f2})<-1}، وxو2{\displaystyle x_{f2}}تصبح نقطة ثابتة غير مستقرة. أي أن تشعبًا آخر يحدث عندر=3{\displaystyle r=3}.

لر=3{\displaystyle r=3}يحدث نوع من التشعب يُعرف باسم تشعب مضاعفة الفترة .ر>3{\displaystyle r>3}لم يعد المدار يتقارب إلى نقطة واحدة، بل يتناوب بين قيم كبيرة وصغيرة حتى بعد مرور فترة زمنية كافية. على سبيل المثال، بالنسبة لـر=3.3{\displaystyle r=3.3}، يتناوب المتغير بين القيمتين 0.4794... و 0.8236....

مخطط العنكبوت والسلسلة الزمنية لـ a = 3.3. ينجذب المدار إلى نقطة مستقرة ذات دورية 2.

يُطلق على المدار الذي يتكرر دوريًا عبر نفس القيم اسم المدار الدوري. في هذه الحالة، يكون السلوك النهائي للمتغير عندما n → ∞ مدارًا دوريًا بفترتين. تُسمى كل قيمة (نقطة) تُكوّن مدارًا دوريًا نقطة دورية. في المثال حيث a = 3.3، فإن 0.4794... و0.8236... هي نقاط دورية. إذا كانت قيمة x معينة نقطة دورية، ففي حالة وجود نقطتين دوريتين، فإن تطبيق الدالة مرتين على x سيعيدها إلى حالتها الأصلية.

إذا طبقنا معادلة الخريطة اللوجستية (1-2) على هذه المعادلة، فسنحصل على

وهذا يعطينا المعادلة التالية من الدرجة الرابعة. حلول هذه المعادلة هي النقاط الدورية. في الواقع، النقطتان الثابتتانxو1=0{\displaystyle x_{f1}=0}وxو2=1-1/ر{\displaystyle x_{f2}=1-1/r} كما أنها تحقق المعادلة (3-4). لذلك، من بين حلول المعادلة (3-5)، اثنان منها يتوافقان معxو1{\displaystyle x_{f1}}وxو2{\displaystyle x_{f2}}والحلان المتبقيان هما نقطتان دوريتان بفترة 2. لنرمز إلى النقطتين الدوريتين بفترة 2 بالرمز التالي:xو1(2){\displaystyle x_{f1}^{(2)}}وxو2(2){\displaystyle x_{f2}^{(2)}}على التوالي. وبحل المعادلة (3-5)، يمكننا الحصول عليها كما يلي

يمكن تطبيق نظرية مماثلة حول استقرار النقاط الثابتة على النقاط الدورية أيضًا. أي أن النقطة الدورية التي تجذب المدارات المحيطة بها تُسمى نقطة دورية مستقرة تقاربياً، بينما تُسمى النقطة الدورية التي تبتعد عنها المدارات المحيطة بها نقطة دورية غير مستقرة. من الممكن تحديد استقرار النقاط الدورية بنفس طريقة تحديد استقرار النقاط الثابتة. في الحالة العامة، لنفترضوك(x){\displaystyle f^{k}(x)}بعد k تكرارًا للخريطة. ليكن(وك)(x){\displaystyle (f^{k})'(x)}ليكن المشتق دوك(x)/دx{\displaystyle df^{k}(x)/dx}النقطة الدورية kxو(ك){\displaystyle x_{f}^{(k)}}. لوxو(ك){\displaystyle x_{f}^{(k)}}يرضي :

ثمxو(ك){\displaystyle x_{f}^{(k)}}مستقر تقاربياً.

ثمxو(ك){\displaystyle x_{f}^{(k)}}غير مستقر.

يمكن فهم المناقشة السابقة حول استقرار النقاط الدورية بسهولة من خلال رسم بياني، تمامًا كما هو الحال مع النقاط الثابتة. في هذا الرسم البياني، يمثل المحور الأفقي xn والمحور الرأسيxن+2{\displaystyle x_{n+2}}ويتم رسم منحنى يوضح العلاقة بينxن+2{\displaystyle x_{n+2}}وxن{\displaystyle x_{n}}تمثل نقاط تقاطع هذا المنحنى مع الخط ذي الزاوية 45° نقاطًا تحقق المعادلة (3-4)، لذا فإن هذه التقاطعات تمثل نقاطًا ثابتة ونقاطًا دورية ثنائية. إذا رسمنا مخططًا بيانيًا للخريطة اللوجستيةو2(x){\displaystyle f^{2}(x)}يمكننا أن نلاحظ أن ميل المماس عند النقطة الثابتةxو2{\displaystyle x_{f2}}يتجاوز 1 عند الحدودر=3{\displaystyle r=3}ويصبح غير مستقر. في الوقت نفسه، يظهر تقاطعان جديدان، وهما النقطتان الدوريتان.xو1(2){\displaystyle x_{f1}^{(2)}}وxو2(2){\displaystyle x_{f2}^{(2)}}.

العلاقة بينxن+2{\displaystyle x_{n+2}}وxن{\displaystyle x_{n}}عندما تكون قيمة r تساوي 2.7، قبل حدوث تشعب مضاعفة الدورة. يتقارب المدار إلى نقطة ثابتةxو2{\displaystyle x_{f2}}.
العلاقة بينxن+2{\displaystyle x_{n+2}}وxن{\displaystyle x_{n}}عندما r = 3. ميل المماس عند النقطة الثابتةxو2{\displaystyle x_{f2}}يساوي 1 بالضبط، ويحدث تشعب مضاعفة الفترة.
العلاقة بينxن+2{\displaystyle x_{n+2}}وxن{\displaystyle x_{n}}عندما تكون قيمة r = 3.3.xو2{\displaystyle x_{f2}}يصبح غير مستقر ويتقارب المدار إلى النقاط الدوريةxو1(2){\displaystyle x_{f1}^{(2)}}وxو2(2){\displaystyle x_{f2}^{(2)}}.

عندما نحسب فعليًا معاملات التفاضل لنقطتين دوريتين للخريطة اللوجستية، نحصل على

عند تطبيق ذلك على المعادلة (3-7)، تصبح قيمة المعامل a كما يلي:

يمكن ملاحظة أن النقاط الدورية ذات البعد 2 تكون مستقرة تقاربياً عندما يكون هذا النطاق3<ر<1+6{\displaystyle 3<r<1+{\sqrt {6}}}أي عندما تتجاوز قيمة r1+6=3.44949...{\displaystyle 1+{\sqrt {6}}=3.44949...}، لم تعد النقاط الدورية ذات الفترة 2 مستقرة تقاربياً ويتغير سلوكها.

تنجذب جميع القيم الأولية تقريبًا في الفترة [0، 1] إلى النقاط الدورية ذات الفترة 2، ولكنxو1=0{\displaystyle x_{f1}=0}وxو2=1-1/أ{\displaystyle x_{f2}=1-1/a} تبقى هذه النقطة ثابتة غير مستقرة في الفترة [0,1]. وتستمر هذه النقاط الثابتة غير المستقرة في البقاء ضمن الفترة [0,1] حتى مع زيادة قيمة r. لذلك، عندما تكون القيمة الابتدائية بالضبطxو1{\displaystyle x_{f1}}أوxو2{\displaystyle x_{f2}}لا ينجذب المدار إلى نقطة دورية ثنائية. علاوة على ذلك، عندما تكون القيمة الأولية هي النقطة الثابتة النهائية لـxو1{\displaystyle x_{f1}}أو النقطة الثابتة النهائية لـxو2{\displaystyle x_{f2}}لا ينجذب المدار إلى نقطة دورية ثنائية. يوجد عدد لا نهائي من هذه النقاط الثابتة النهائية في الفترة [0، 1]. ومع ذلك، فإن عدد هذه النقاط ضئيل للغاية مقارنةً بمجموعة الأعداد الحقيقية [0، 1].

الحالة عندما يكون 3.44949 ≤ r ≤ 3.56995

عندما تتراوح قيمة r بين 3.44949 و3.54409 (تقريبًا)، فإنّ التوزيع السكاني، في جميع الظروف الأولية تقريبًا، سيقترب من تذبذبات دائمة بين أربع قيم. وهذه القيمة الأخيرة هي جذر لكثير حدود من الدرجة الثانية عشرة (المتتالية A086181 في OEIS ) .

مع ازدياد قيمة r لتتجاوز 3.54409، تقترب أعداد السكان، في معظم الحالات الابتدائية، من التذبذب بين 8 قيم، ثم 16، ثم 32، وهكذا. تتناقص أطوال فترات المعلمات التي تُنتج تذبذبات بطول مُحدد بسرعة؛ وتقترب النسبة بين أطوال فترتي تشعب متتاليتين من ثابت فيجنباوم δ ≈ 4.66920 . يُعد هذا السلوك مثالًا على سلسلة مضاعفة الدورة .

عندما تتجاوز قيمة المعامل r1+6=3.44949...{\displaystyle 1+{\sqrt {6}}=3.44949...}تصبح النقاط الدورية الثنائية المستقرة سابقًا غير مستقرة، وتتولد نقاط دورية رباعية مستقرة، وينجذب المدار نحو تذبذب دوري رباعي. أي أن تشعب مضاعفة الدورة يحدث مرة أخرى عندر=3.44949...{\displaystyle r=3.44949...}قيمة x عند النقطة الدورية ذات 4 هي أيضًا

بما أن المعادلة (3-11) تحقق الشرط، فإن حلها يسمح بإيجاد قيم x عند النقاط الدورية ذات الأربع دورات. مع ذلك، فإن المعادلة (3-11) هي معادلة من الدرجة السادسة عشرة، وحتى لو استخرجنا الحلول الأربعة للنقاط الثابتة والنقاط الدورية ذات الدورتين، فإنها تبقى معادلة من الدرجة الثانية عشرة. لذا، لم يعد من الممكن حل هذه المعادلة للحصول على دالة صريحة لـ a تمثل قيم النقاط الدورية ذات الأربع دورات بنفس طريقة تمثيلها للنقاط الدورية ذات الدورتين.

أمثلة على نقاط التفرع حتى 256 دورة
الفرع رقم kالفترة 2 كيلونقطة التفرع أ ك
123.0000000
243.4494896
383.5440903
4163.5644073
5323.5687594
6643.5696916
71283.5698913
82563.5699340

مع ازدياد قيمة a، تخضع النقطة المستقرة ذات الدورة الرباعية لمضاعفة دورية أخرى، مما ينتج عنه نقطة مستقرة ذات دورة ثمانية. ومع ازدياد قيمة a، تحدث تفرعات مضاعفة الدورة بشكل لانهائي: 16، 32، 64، ... وهكذا، حتى الوصول إلى دورة لانهائية، أي مدار لا يعود أبدًا إلى قيمته الأصلية. تُسمى هذه السلسلة اللانهائية من تفرعات مضاعفة الدورة بالتتابع. وبينما تحدث هذه التفرعات بشكل لانهائي، تتناقص الفترات بين قيم a التي تحدث عندها بشكل هندسي. وبالتالي، يحدث عدد لا نهائي من تفرعات مضاعفة الدورة قبل أن تصل قيمة المعامل a إلى قيمة محدودة. لنفترض أن التفرع من الدورة 1 إلى الدورة 2 الذي يحدث عند r = 3 يُعتبر أول تفرع لمضاعفة الدورة. عندئذٍ، في هذا التتابع من تفرعات مضاعفة الدورة، تظهر نقطة مستقرة ذات دورة 2k عند نقطة التفرع رقم k. لنفترض أن نقطة التشعب رقم k هي a، ونرمز لها بـ a k. في هذه الحالة، من المعروف أنرك{\displaystyle r_{k}}يتقارب إلى القيمة التالية عندما k → ∞. (المتتالية A098587 في OEIS )

علاوة على ذلك، من المعروف أن معدل انخفاض ak يصل إلى قيمة ثابتة في النهاية، كما هو موضح في المعادلة التالية.

تُسمى هذه القيمة لـ δ ثابت فيجنباوم لأنه تم اكتشافه بواسطة الفيزيائي الرياضي ميتشل فيجنباوم. ر{\displaystyle r_{\infty }}تُسمى هذه النقطة بنقطة فيجنباوم. في سلسلة مضاعفة الفترة،وم{\displaystyle f^{m}} وو2م{\displaystyle f^{2m}}تتميز هذه الدوال بأنها متطابقة محليًا بعد إجراء تحويل قياسي مناسب. ويمكن إيجاد ثابت فيجنباوم باستخدام تقنية تُسمى إعادة التطبيع، والتي تستغل هذا التشابه الذاتي. كما أن الخصائص التي تُظهرها الدالة اللوجستية في سلسلة مضاعفة الدورة هي خصائص عامة في فئة أوسع من الدوال، كما سيتم توضيحه لاحقًا.

للحصول على نظرة عامة على السلوك النهائي للمدار عند قيمة معينة لمعامل محدد، يُعدّ مخطط التشعب التقريبي، أو ما يُعرف بمخطط المدار، مفيدًا. في هذا المخطط، يُمثّل المحور الأفقي المعامل r، بينما يُمثّل المحور الرأسي المتغير x، كما هو الحال في مخطط التشعب. باستخدام الحاسوب، تُحدّد المعاملات، وتُجرى، على سبيل المثال، 500 تكرار. بعد ذلك، تُهمل أول 100 نتيجة، وتُرسم نتائج الـ 400 المتبقية فقط. يسمح هذا بتجاهل السلوك الانتقالي الأولي، مع الحفاظ على السلوك التقاربي للمدار. على سبيل المثال، عند رسم نقطة واحدة لـ r، تُمثّل نقطة ثابتة، وعند رسم m نقطة لـ r، فإنها تُشير إلى مدار دوري ذي m دورة. عند رسم مخطط المدار للخريطة اللوجستية، يُمكن ملاحظة كيفية انقسام الفرع الذي يُمثّل المدار الدوري المستقر، والذي يُمثّل سلسلة من تشعبات مضاعفة الدورة.

مخطط تشعب سلسلة تشعب مضاعفة الفترة التي تحدث بين المعلماتر1=3{\displaystyle r_{1}=3}وأ=3.56994...{\displaystyle a_{\infty }=3.56994...}بعد 64 فترة (أ5{\displaystyle a_{5}}تصبح المسافة ضيقة للغاية وتكاد تنهار.

عندما يكون المعاملر=ر{\displaystyle r=r_{\infty }}إنها بالضبط نقطة تراكم سلسلة مضاعفة الفترة، المتغيرxن{\displaystyle x_{n}}ينجذب إلى مدارات غير دورية لا تغلق أبدًا. بعبارة أخرى، توجد نقطة دورية ذات دورة لانهائية عندر{\displaystyle r_{\infty }}يُطلق على هذا المدار غير الدوري اسم جاذب فيجنباوم .2{\displaystyle 2^{\infty }}الجاذب. الجاذب مصطلح يُستخدم للإشارة إلى منطقة تتميز بخاصية جذب المدارات المحيطة بها، وهو المدار الذي ينجذب إليها في النهاية ويستمر فيه. النقاط الثابتة الجاذبة والنقاط الدورية المذكورة أعلاه هي أيضاً من عائلة الجاذبات.

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

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

الحالة عندما يكون 3.56995 < r < 4

ملخص نوعي
تطور الشروط الأولية المختلفة كدالة لـ r (المعامل k من الشكل يتوافق مع المعامل r من التعريف الوارد في المقالة.)
تطور الشروط الأولية المختلفة كدالة لـ r مع الانحياز (المعامل k من الشكل يتوافق مع المعامل r من التعريف الوارد في المقالة.)
  • عند قيمة r ≈ 3.56995 (التسلسل A098587 في OEIS ) ، تبدأ الفوضى، في نهاية سلسلة مضاعفة الدورة. انطلاقًا من جميع الشروط الابتدائية تقريبًا، لم نعد نرى تذبذبات ذات دورة محدودة. تؤدي الاختلافات الطفيفة في عدد السكان الأولي إلى نتائج مختلفة جذريًا بمرور الوقت، وهي سمة أساسية للفوضى.
  • يُقارن هذا الرقم ويُفهم على أنه مكافئ لرقم رينولدز لبداية الظواهر الفوضوية الأخرى، مثل الاضطراب، ومماثل لدرجة الحرارة الحرجة للانتقال الطوري . في جوهره، يحتوي فضاء الطور على فضاء فرعي كامل من الحالات ذات متغيرات ديناميكية إضافية لتوصيف الحالة المجهرية للنظام، ويمكن فهم هذه المتغيرات على أنها دوامات في حالة الاضطراب، ومعاملات ترتيب في حالة الانتقالات الطورية .
  • تُظهر معظم قيم r التي تتجاوز 3.56995 سلوكًا فوضويًا، ولكن لا تزال هناك نطاقات معزولة معينة من r تُظهر سلوكًا غير فوضوي؛ تُسمى هذه النطاقات أحيانًا بجزر الاستقرار . على سبيل المثال، بدءًا من 1  + √8 [ 8 ] (حوالي 3.82843)، يوجد نطاق من قيم r يُظهر تذبذبًا بين ثلاث قيم، وبالنسبة لقيم r الأعلى قليلًا ، يُظهر تذبذبًا بين ست قيم، ثم اثنتي عشرة قيمة، وهكذا. 
  • فير=1+8=3.8284...{\displaystyle r=1+{\sqrt {8}}=3.8284...}[ 9 ] تظهر الدورة المستقرة ذات الفترة 3.
  • يُطلق على تطور السلوك الفوضوي للمتتالية اللوجستية مع تغير قيمة المعامل r من حوالي 3.56995 إلى حوالي 3.82843 اسم سيناريو بوميو-مانفيل ، والذي يتميز بمرحلة دورية (صفائحية) تتخللها فترات من السلوك غير الدوري. يُستخدم هذا السيناريو في أجهزة أشباه الموصلات. [ 10 ] توجد نطاقات أخرى تُنتج تذبذبًا بين 5 قيم، إلخ؛ وتحدث جميع فترات التذبذب لبعض قيم r . نافذة مضاعفة الدورة ذات المعامل c هي نطاق من قيم r يتكون من سلسلة من النطاقات الفرعية. يحتوي النطاق الفرعي k على قيم r التي تُشكل دورة مستقرة (دورة تجذب مجموعة من النقاط الأولية ذات قياس الوحدة) دورتها 2kc . تُسمى هذه السلسلة من النطاقات الفرعية سلسلة التوافقيات . [ ماي، روبرت م. (1976) 1 ] في نطاق فرعي ذي دورة مستقرة دورية مقدارها 2k * c ، توجد دورات غير مستقرة دورية مقدارها 2kc لجميع قيم k < k * . تُسمى قيمة r في نهاية التسلسل اللانهائي للنطاقات الفرعية بنقطة تراكم سلسلة التوافقيات. مع ازدياد قيمة r ، تظهر سلسلة من النوافذ الجديدة بقيم c مختلفة . النافذة الأولى هي عند c = 1 ؛ وتظهر جميع النوافذ اللاحقة التي تتضمن قيم c فردية بترتيب تنازلي لـ c بدءًا من قيمة c كبيرة كيفما كانت . [ ماي، روبرت م. (1976) 1 ] [ 11 ]
  • فير=3.678...،x=0.728...{\displaystyle r=3.678...,x=0.728...}يتقاطع نطاقان فوضويان من مخطط التشعب عند نقطة ميسيوريفيتش الأولى للخريطة اللوجستية. ويحقق ذلك المعادلات التالية:ر3-2ر2-4ر-8=0،x=1-1/ر{\displaystyle r^{3}-2r^{2}-4r-8=0,x=1-1/r}[ 12 ]
  • بعد تجاوز قيمة r = 4 ، تخرج جميع القيم الابتدائية تقريبًا من الفترة [0,1] وتتباعد. تشكل مجموعة الشروط الابتدائية التي تبقى ضمن الفترة [0,1] مجموعة كانتور ، وتكون الديناميكيات المقيدة بهذه المجموعة فوضوية. [ 13 ]

لكل قيمة لـ توجد دورة مستقرة واحدة على الأكثر. إذا وُجدت دورة مستقرة، فهي مستقرة عالميًا، وتجذب جميع النقاط تقريبًا. [ 14 ] : 13 بعض قيم r التي لها دورة مستقرة ذات فترة معينة، لها عدد لا نهائي من الدورات غير المستقرة ذات فترات مختلفة.

مخطط التشعب للخريطة اللوجستية. يظهر الجاذب لأي قيمة للمعامل r على الخط الرأسي عند تلك القيمة r .

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

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

تكبير المنطقة الفوضوية من الخريطة
مناطق مستقرة داخل المنطقة الفوضوية، حيث يحدث تشعب مماس عند الحدود بين الجاذب الفوضوي والدوري، مما يعطي مسارات متقطعة كما هو موضح في سيناريو بوميو-مانفيل

يمكننا أيضًا النظر في القيم السالبة لـ r :

  • بالنسبة لقيم r بين -2 و -1، فإن المتتالية اللوجستية تتميز أيضاً بسلوك فوضوي. [ 6 ]
  • عندما تكون قيمة r بين -1 و 1  - √6 ، وعندما تكون قيمة x0 بين 1/ r و 1 - 1/ r ، فإن عدد السكان سيقترب من تذبذبات دائمة بين قيمتين، كما هو الحال عندما تكون قيمة r بين 3 و 1 + √6 ، ويتم حسابها بنفس الصيغة. [ 6 ]   
ظهور الفوضى
مدارات فوضوية للخريطة اللوجستية عندما تكون قيمة r = 3.82. المربعات البرتقالية هي مدارات تبدأ منx0=0.1234{\displaystyle x_{0}=0.1234}والدوائر الزرقاء والخضراء هي مدارات تبدأ منx^0=0.1234+10-9{\displaystyle {\hat {x}}_{0}=0.1234+10^{-9}}.
المسار بدءًا من x₀ = 0.1234 و ˆx، والفرق في المدارات بدءًا منx0=0.1234+10-9{\displaystyle x_{0}=0.1234+10^{-9}}ينمو بشكل أُسّي. المحور الرأسي هوΔxن=|xن-x^ن|{\displaystyle \Delta x_{n}=|x_{n}-{\hat {x}}_{n}|}، معروضة على مقياس لوغاريتمي .

عندما تتجاوز قيمة المعامل rر=3.56994...{\displaystyle r_{\infty }=3.56994...}تُظهر الخريطة اللوجستية سلوكًا فوضويًا. وبشكل عام، الفوضى هي سلوك معقد وغير منتظم يحدث على الرغم من أن معادلة الفرق التي تصف الخريطة اللوجستية لا تحتوي على أي غموض احتمالي، وأن الحالة التالية محددة بشكل كامل وفريد. نطاقر>ر{\displaystyle r>r_{\infty }}يُطلق على جزء من الخريطة اللوجستية اسم المنطقة الفوضوية.

من خصائص الفوضى عدم القدرة على التنبؤ بها، وهو ما يرمز إليه مصطلح " تأثير الفراشة" . ويعود ذلك إلى خاصية الفوضى التي تجعل اختلافًا طفيفًا في الحالة الابتدائية يؤدي إلى اختلاف كبير في الحالة اللاحقة. ففي سياق نظام ديناميكي منفصل، إذا كان لدينا قيمتان ابتدائيتانx0{\displaystyle x_{0}}وx^0{\displaystyle {\hat {x}}_{0}} بغض النظر عن مدى قربهم، بمجرد أن يتقدم الزمن إلى حد معين، فإن كل وجهةxن{\displaystyle x_{n}}وx^ن{\displaystyle {\hat {x}}_{n}}قد تختلف بشكل كبير. على سبيل المثال، استخدمر=3.95،x0=0.1،x^0=x0+10-9{\displaystyle r=3.95,x_{0}=0.1,{\hat {x}}_{0}=x_{0}+10^{-9}} إذا تم حساب المدارات باستخدام قيمتين أوليتين متشابهتين للغاية، 0 = 0.1000000001، فإن الفرق ينمو إلى قيم كلية يمكن رؤيتها بوضوح على الرسم البياني بعد حوالي 29 تكرارًا.

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

هنا، تعني كلمة log اللوغاريتم الطبيعي . λ هي المسافة بين المدارين (xن{\displaystyle x_{n}}و x^ن{\displaystyle {\hat {x}}_{n}}تشير القيمة الموجبة لـ λ إلى أن النظام حساس للشروط الابتدائية، بينما تشير القيمة الصفرية أو السالبة إلى أن النظام غير حساس لها. عند حساب λ عدديًا، يمكن التأكد من أن λ تبقى ضمن نطاق القيم الصفرية أو السالبة.ر<ر{\displaystyle r<r_{\infty }}وأن قيمة λ يمكن أن تأخذ قيمًا موجبة في النطاقر>ر{\displaystyle r>r_{\infty }}.

نافذة، متقطعة

حتى أبعد من ذلكر{\displaystyle r_{\infty }}لا يعتمد السلوك ببساطة على المعامل r. تكمن العديد من البنى الرياضية المعقدة في المنطقة الفوضوية لـر>ر{\displaystyle r>r_{\infty }}في هذه المنطقة ، لا تدوم الفوضى إلى الأبد؛ إذ تعود المدارات الدورية المستقرة للظهور. سلوك لـر<أ4{\displaystyle r_{\infty }<a\leq 4}يمكن تقسيمها بشكل عام إلى نوعين:

  • نقطة دورية مستقرة: في هذه الحالة، يكون أس ليابونوف سالبًا.
  • المدارات غير الدورية: في هذه الحالة، يكون أس ليابونوف موجبًا.

منطقة النقاط الدورية المستقرة الموجودة لـ rر<ر4{\displaystyle r_{\infty }<r\leq 4}تُسمى هذه المنطقة نافذة دورية، أو ببساطة نافذة. إذا نظرنا إلى منطقة فوضوية في مخطط مداري، فإن منطقة المدارات غير الدورية تبدو كسحابة من نقاط لا حصر لها، حيث تمثل النوافذ الفراغات المتناثرة المحاطة بهذه السحابة.

مخطط مداري للخريطة اللوجستية من r = 3.55 إلى r = 4 (يشار إلى المعلمة بـ r في المخطط)

في كل نافذة، سلسلة من التشعبات التي تضاعفت فيها الفترة والتي حدثت من قبلر=3.56994...{\displaystyle r_{\infty }=3.56994...}يتكرر الأمر. ولكن بدلاً من المدارات الدورية المستقرة السابقة ذات قيمة 2^k، تتولد مدارات دورية مستقرة جديدة مثل 3×2^k و5×2^k. النافذة الأولى لها دورة مقدارها p، وتُسمى النوافذ التي يحدث عندها تضاعف الدورة بنوافذ الدورة p، وهكذا. على سبيل المثال، توجد نافذة دورتها 3 في المنطقة المحيطة بـ 3.8284 < a < 3.8415، وضمن هذه المنطقة تكون مضاعفات الدورة كالتالي: 3، 6، 12، 24، ...، 3×2^k، ....

فوضى عابرة عند a = 3.8285. يتصرف النظام بشكل فوضوي حتى ينجذب إلى مدار دوري 3.

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

يوجد عدد لا نهائي من النوافذ في النطاق a∞ < a < 4. تتميز هذه النوافذ بفترات زمنية مختلفة، حيث توجد نافذة ذات فترة زمنية لكل عدد طبيعي أكبر من أو يساوي ثلاثة. مع ذلك، لا تظهر كل نافذة مرة واحدة فقط. فكلما زادت قيمة p، زاد تكرار ظهور نافذة بتلك الفترة. تظهر نافذة ذات فترة زمنية 3 مرة واحدة فقط، بينما تظهر نافذة ذات فترة زمنية 13 مرة 315 مرة. عند ظهور مدار دوري بقيمة 3 في النافذة ذات الفترة الزمنية 3، يكتمل ترتيب سزاركوفسكي، وبذلك تُشاهد جميع المدارات بجميع الفترات الزمنية.

إذا اقتصرنا على الحالة التي يكون فيها p عددًا أوليًا، فإن عدد النوافذ ذات الفترة p هو

تم اشتقاق هذه الصيغة لـ p ليكون عددًا أوليًا، ولكن في الواقع من الممكن حساب عدد النقاط الدورية المستقرة p بدقة جيدة لـ p غير الأولي أيضًا.

يكون عرض النافذة (الفرق بين نقطة بداية النافذة ونقطة نهايتها) أوسع ما يمكن للنوافذ ذات الدورة 3، ويضيق مع زيادة الدورة. على سبيل المثال، يبلغ عرض النافذة ذات الدورة 13 حوالي 3.13 × 10⁻⁶. وتشير التقديرات التقريبية إلى أن حوالي 10% من[ر،4]{\displaystyle [r_{\infty },4]}تقع في منطقة النافذة، بينما تهيمن المدارات الفوضوية على الباقي.

يحدث التحول من الفوضى إلى نافذة مع زيادة قيمة r نتيجةً لتشعب مماس، حيث يكون منحنى الخريطة مماسًا لقطر الخط y = x لحظة التشعب، وتؤدي تغييرات المعلمات اللاحقة إلى نقطتين ثابتتين يتقاطع فيهما المنحنى مع الخط. بالنسبة لنافذة دورية p، تكون الخريطة المتكررةوص(x){\displaystyle f^{p}(x)}يُظهر هذا النظام تشعبًا مماسًا، مما ينتج عنه مدارات مستقرة دورية من الرتبة p. القيمة الدقيقة لنقطة التشعب لنافذة دورية من الرتبة 3 معروفة، وإذا كانت قيمة نقطة التشعب هذه r هير3{\displaystyle r_{3}}، ثمر3=1+8=3.828427...{\displaystyle r_{3}=1+{\sqrt {8}}=3.828427...}يمكن فهم الخطوط العريضة لهذا التشعب من خلال النظر إلى الرسم البياني لـو3(x){\displaystyle f^{3}(x)}(المحور الرأسي)xن+3{\displaystyle x_{n+3}}المحور الأفقيxن{\displaystyle x_{n}}).

رسم بياني لـو3(x){\displaystyle f^{3}(x)}عندما تكون قيمة r أقل بقليل من 3. لا يكون الرسم البياني مماسًا إلا عند النقاط الثابتة، ولا توجد نقاط دورية مدتها 3.
عندما تكون قيمة a هي 3 بالضبط، فإن الرسم البياني يمس القطر عند ثلاث نقاط بالضبط، مما ينتج عنه ثلاث نقاط دورية.
عندما تكون قيمة a أكبر بقليل من 3، يمر الرسم البياني بالقطر وينقسم إلى نقاط دورية مستقرة وغير مستقرة ذات 3 دورات.

عندما ننظر إلى سلوكxن{\displaystyle x_{n}}عندما تكون قيمة r = 3.8282، وهي أصغر قليلاً من نقطة التفرعر3{\displaystyle r_{3}}نلاحظ أنه بالإضافة إلى التغيرات غير المنتظمة، يوجد سلوك يتغير دوريًا بثلاث دورات تقريبًا، وتحدث هذه الدورات بالتناوب. يُطلق على هذا النوع من السلوك الدوري اسم "التدفق الصفائحي"، بينما يُطلق على السلوك غير المنتظم اسم "الانفجار"، قياسًا على الموائع. لا يوجد انتظام في طول الدورات الزمنية للانفجارات والتدفقات الصفائحية، وتتغير بشكل غير منتظم. ومع ذلك، عندما نلاحظ السلوك عند r = 3.828327، وهو أقرب إلىر3{\displaystyle r_{3}}يكون متوسط ​​طول التدفقات الصفائحية أطول، بينما يكون متوسط ​​طول الانفجارات أقصر مما هو عليه عندما تكون قيمة r تساوي 3.8282. وإذا زدنا قيمة r أكثر، فإن طول التدفقات الصفائحية يزداد أكثر فأكثر، وعندر3{\displaystyle r_{3}}يتحول إلى ثلاث فترات مثالية.

سلسلة زمنية عندما r = 3.8282
سلسلة زمنية عندما r = 3.828327

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

أنماط القنوات التي تظهر في الرسم البياني للدالة f3 (x)

يمكن فهم آلية التقطع أيضًا من خلال الرسم البياني للخريطة. عندمار{\displaystyle r}أصغر قليلاً منر3{\displaystyle r_{3}}هناك فجوة صغيرة جدًا بين الرسم البياني لـو3(x){\displaystyle f^{3}(x)}والقطري . تُسمى هذه الفجوة قناة، وتحدث العديد من تكرارات الخريطة أثناء مرور المدار عبر القناة الضيقة. أثناء المرور عبر هذه القناة،xن{\displaystyle x_{n}}وxن+3{\displaystyle x_{n+3}}تتقارب المسارات بشدة، وتتغير المتغيرات بشكل أشبه بدوران مداري ثلاثي دوري. وهذا يُشابه التدفق الصفائحي. يغادر المسار في النهاية القناة الضيقة، لكنه يعود إليها مجددًا نتيجةً للبنية العامة للخريطة. وأثناء مغادرته القناة، يتصرف بشكل فوضوي. وهذا يُشابه الانفجار.

شريط، تشطيب نافذة

بالنظر إلى المجال الفوضوي بأكمله، سواء كان فوضويًا أو مُقسَّمًا إلى نافذة، فإن القيمتين القصوى والدنيا على المحور الرأسي لمخطط المدار (الحدود العليا والسفلى للمُستقطِب) محصورتان ضمن نطاق مُحدد. وكما هو موضح في المعادلة (2-1)، تُعطى القيمة القصوى للخريطة اللوجستية بالقيمة r/4، وهي الحد الأعلى للمُستقطِب. أما الحد الأدنى للمُستقطِب فيُعطى بالنقطة f(r/4) حيث تُرسم القيمة r/4. في النهاية، تعتمد القيمتان القصوى والدنيا اللتان يتحرك عندهما xn على مخطط المدار على المعامل r.

وأخيرًا، بالنسبة لـ r = 4، يمتد المدار على النطاق الكامل [0، 1].

عند رصد خريطة مدارية، يتسم توزيع النقاط بتظليل مميز. تشير المناطق الداكنة إلى أن المتغير يأخذ قيمًا في جوار هذه المناطق، بينما تشير المناطق الفاتحة إلى أن المتغير يأخذ قيمًا في جوار المناطق الداكنة. تعود هذه الاختلافات في تردد النقاط إلى شكل الرسم البياني للخريطة اللوجستية. يجذب الجزء العلوي من الرسم البياني، بالقرب من r/4، مدارات ذات تردد عالٍ، وتصبح المنطقة القريبة من f(r/4) التي تُرسم من هناك ذات تردد عالٍ أيضًا، والمنطقة القريبة منو2(ر/4){\displaystyle f^{2}(r/4)}تصبح النقاط التي يتم رسمها من هناك متكررة للغاية، وهكذا. يتميز توزيع كثافة النقاط الناتجة عن الخريطة بكمية تسمى المقياس الثابت أو دالة التوزيع، ويكون المقياس الثابت للمستقطب قابلاً للتكرار بغض النظر عن القيمة الأولية.

بالنظر إلى بداية المنطقة الفوضوية في مخطط المدار، مباشرةً بعد نقطة التراكمر=3.56994{\displaystyle r_{\infty }=3.56994}في الدورة الأولى - دورة التضاعف المتتالية - يُلاحظ أن المدار مُقسّم إلى عدة مناطق فرعية تُسمى نطاقات . عند وجود نطاقات متعددة، يتحرك المدار عبر كل نطاق بترتيب منتظم، لكن القيم داخل كل نطاق غير منتظمة. تُسمى هذه المدارات الفوضوية بفوضى النطاقات أو الفوضى الدورية، وتُسمى الفوضى التي تحتوي على k نطاقًا بفوضى النطاقات k. تقع فوضى النطاقين تقريبًا في النطاق 3.590 < r < 3.675.

بنية النطاق. لأنهـص{\displaystyle e_{p}}مع تناقص المسافة بين الخطوط بسرعة، لا يمكن عرض أكثر من ثمانية نطاقات. تقع الخطوط العلوية والسفلية، التي تحتوي على المدارات، ضمن نطاق المعادلة (3-16).

مع انخفاض قيمة r أكثر من الطرف الأيسر للفوضى ثنائية النطاق، r = 3.590، يتضاعف عدد النطاقات، تمامًا كما هو الحال في تشعب مضاعفة الفترة. لنفترضهـص{\displaystyle e_{p}}(لـ p = 1، 2، 4، ...، 2k، ...) تُشير إلى نقاط التشعب حيث ينقسم فوضى النطاق p − 1 إلى فوضى النطاق p، أو حيث تندمج فوضى النطاق p في فوضى النطاق p − 1. بعد ذلك، وكما هو الحال في تشعب مضاعفة الفترة، يتراكم ep إلى قيمة معينة عندما p → ∞. عند نقطة التراكم هذههـ{\displaystyle e_{\infty }}يصبح عدد النطاقات لانهائيًا، وقيمةهـ{\displaystyle e_{\infty }}يساوي قيمةر{\displaystyle r_{\infty }}.

بنية هرمية ذاتية التشابه لخريطة المسار الكاملة للخريطة اللوجستية

وبالمثل، بالنسبة لنقاط التشعب في سلسلة تشعب مضاعفة الدورة التي ظهرت قبل a∞، فلنرمز بـ ap (حيث p = 1، 2، 4، ...، 2k، ...) إلى نقاط التشعب التي تتفرع عندها p مدارات دورية مستقرة إلى p + 1 مدارات دورية مستقرة. في هذه الحالة، إذا نظرنا إلى مخطط المدارات منر2{\displaystyle r_{2}}لهـ2{\displaystyle e_{2}}يوجد نسختان مبسطتان من مخطط المدار العالمي منر1{\displaystyle r_{1}}لهـ1{\displaystyle e_{1}}في مخطط المدارات منر2{\displaystyle r_{2}}لهـ2{\displaystyle e_{2}}وبالمثل، إذا نظرنا إلى مخطط المدارات منر4{\displaystyle r_{4}}لهـ4{\displaystyle e_{4}}يوجد أربعة إصدارات مُختزلة من مخطط المدارات العالمي من a1 إلى e1 في مخطط المدارات منر4{\displaystyle r_{4}}لهـ4{\displaystyle e_{4}}وبالمثل ، هناك نسخ مخفضة من مخطط المدار العالمي في مخطط المدار من ap إلى ep، ويحتوي هيكل التفرع للخريطة اللوجستية على تسلسل هرمي متشابه ذاتيًا لانهائي.

بنية هرمية ذاتية التشابه لنوافذ الخريطة اللوجستية

توجد أيضًا سلسلة هرمية ذاتية التشابه من هياكل التشعب داخل النوافذ. تتبع سلاسل تشعب مضاعفة الدورة داخل النافذة نفس مسار سلاسل تشعب الدورة 2k. أي أن هناك عددًا لا نهائيًا من تشعبات مضاعفة الدورة داخل النافذة، وبعد ذلك يصبح السلوك فوضويًا مرة أخرى. على سبيل المثال، في نافذة ذات دورة 3، تنتهي سلسلة المدارات الدورية المستقرة عندأ3{\displaystyle a_{3\infty }}≈ 3.8495. بعدأ3{\displaystyle a_{3\infty }}عندما تكون قيمة a ≈ 3.8495، يصبح السلوك فوضى نطاقية من مضاعفات العدد ثلاثة. ومع زيادة قيمة a منأ3{\displaystyle a_{3\infty }}تتداخل هذه الفوضى النطاقية أيضًا في أزواج، حتى يصبح لدينا في نهاية النافذة ثلاثة نطاقات. وضمن هذه النطاقات داخل النافذة، يوجد عدد لا نهائي من النوافذ. في النهاية، تحتوي النافذة على نسخة مصغرة من مخطط المدار الكامل لـ 1 ≤ a ≤ 4، وضمن النافذة يوجد تسلسل هرمي متشابه ذاتيًا من التفرعات.

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

عندما r = 4

مخطط العنكبوت للخريطة اللوجستية بمعامل r = 4 (يسار) وسلسلة زمنية تصل إلى n = 500 (يمين) للقيمة الأوليةx0{\displaystyle x_{0}}= 0.3.

عندما تكون قيمة المعامل r تساوي 4، يصبح السلوك فوضويًا على كامل النطاق [0، 1]. في هذه الحالة، يصل مُعامل ليابونوف λ إلى أقصى قيمة له، وتكون الحالة في أقصى درجات الفوضى. يمكن حساب قيمة λ للخريطة اللوجستية عند r = 4 بدقة، وهي λ = log 2. على الرغم من عدم وجود تعريف رياضي دقيق للفوضى حتى الآن، إلا أنه يمكن إثبات أن الخريطة اللوجستية عند r = 4 فوضوية على النطاق [0، 1] وفقًا لأحد التعريفات المعروفة للفوضى.

رسم بياني للمقياس الثابت ρ(x) لـ r = 4. يوضح الرسم البياني النقطي التردد الفعلي للنقاط التي تم الحصول عليها على مدار 10000 تكرار (مع قياس الارتفاع وفقًا لـ ρ (x)).

يمكن أيضًا إعطاء المقياس الثابت لكثافة النقاط، ρ(x)، بواسطة الدالة الدقيقة ρ(x) لـ r = 4 :

هنا، تعني ρ(x) أن نسبة النقاط xn التي تقع في الفترة المتناهية الصغر [x,x+dx] عند تكرار الخريطة تُعطى بالعلاقة ρ(x) dx. يتميز التوزيع التكراري للخريطة اللوجستية مع r = 4 بكثافة عالية بالقرب من جانبي [0, 1]، وبأقل كثافة عند x = 0.5.

عندما تكون قيمة r تساوي 4، فإنه بالإضافة إلى المدارات الفوضوية، توجد أيضًا مدارات دورية ذات فترات زمنية مختلفة. بالنسبة لعدد طبيعي n، يكون الرسم البياني لـور=4ن(x){\displaystyle f_{r=4}^{n}(x)}هو منحنى مع2ن-1{\displaystyle 2^{n-1}}القمم و2ن-1-1{\displaystyle 2^{n-1}-1}الوديان، وكلها مماسة للنقطتين 0 و1. وبالتالي، فإن عدد نقاط التقاطع بين القطر والرسم البياني هو2ن{\displaystyle 2^{n}}وهناك2ن{\displaystyle 2^{n}}نقاط ثابتة منون(x){\displaystyle f^{n}(x)}تُدرج النقاط الدورية ذات الفترة n دائمًا في هذه النقاط .2ن{\displaystyle 2^{n}}النقاط الثابتة، لذا فإن أي مدار دوري n موجود لـور=4ن(x){\displaystyle f_{r=4}^{n}(x)}وبالتالي ، عندما تكون قيمة r تساوي 4، يوجد عدد لا نهائي من النقاط الدورية على الفترة [0، 1]، لكن جميع هذه النقاط الدورية غير مستقرة. علاوة على ذلك، في المجموعة اللانهائية غير القابلة للعد في الفترة [0، 1]، يكون عدد النقاط الدورية لانهائيًا قابلًا للعد، ولذلك فإن جميع المدارات تقريبًا التي تبدأ من قيم ابتدائية ليست دورية بل غير دورية.

إذا قمنا بتحويل مدار الخريطة اللوجستيةور=4{\displaystyle f_{r=4}}بتحويلها إلى سلسلة من الأصفار والآحاد، يمكننا إعادة إنتاج أي سلسلة من الرموز.

من أهم جوانب الفوضى طبيعتها المزدوجة: الحتمية والعشوائية. تُعدّ الأنظمة الديناميكية عمليات حتمية، ولكن عند تقريب نطاق المتغيرات بشكل مناسب، تصبح غير قابلة للتمييز عن العمليات العشوائية. في حالة الخريطة اللوجستية ذات r = 4، يمكن وصف نتيجة كل رمية عملة معدنية بمسار الخريطة اللوجستية. ويمكن توضيح ذلك على النحو التالي.

لنفترض أن قطعة نقدية تُرمى باحتمالية 1/2 للظهور على صورة أو كتابة، وتُرمى القطعة بشكل متكرر. إذا كانت الصورة تساوي 0 والكتابة تساوي 1، فإن نتيجة رميات متكررة (صورة، كتابة، صورة، كتابة، إلخ) ستكون سلسلة من الرموز مثل 01001.... من ناحية أخرى، بالنسبة للمسارx0،x1،x2،...{\displaystyle x_{0},x_{1},x_{2},...}في الخريطة اللوجستية، تُحوّل القيم الأقل من x = 0.5 إلى 0، والقيم الأكبر من x = 0.5 إلى 1، ويُستبدل المسار بسلسلة رموز تتكون من أصفار وواحدات. على سبيل المثال، إذا كانت القيمة الأولية هيx0=0.2{\displaystyle x_{0}=0.2}، ثمx1=0.64{\displaystyle x_{1}=0.64}،x2=0.9216{\displaystyle x_{2}=0.9216}،x3=0.28901{\displaystyle x_{3}=0.28901}...، لذا سيكون المسار هو سلسلة الرموز 0110.... لنفترضSج{\displaystyle S_{C}}ليكن سلسلة الرموز الناتجة عن رمية العملة السابقة، وSل{\displaystyle S_{L}}لتكن سلسلة الرموز الناتجة عن الخريطة اللوجستية الأخيرة. الرموز الموجودة في سلسلة الرموزSج{\displaystyle S_{C}}تم تحديدها عن طريق رمي عملة معدنية عشوائيًا، لذا فإن أي نمط من أنماط تسلسل الأرقام ممكن. لذلك، مهما كانت السلسلةSل{\displaystyle S_{L}}توجد خريطة لوجستية مماثلة فيSج{\displaystyle S_{C}}والأمر "اللافت" هو أن العكس صحيح أيضاً: مهما كانت سلسلةSج{\displaystyle S_{C}}ويمكن تحقيق ذلك من خلال مسار الخريطة اللوجستيةSل{\displaystyle S_{L}}باختيار القيم الأولية المناسبة. أي، لأيSج{\displaystyle S_{C}}توجد نقطة فريدةx0{\displaystyle x_{0}}في الفترة [0، 1] بحيثSج=Sل{\displaystyle S_{C}=S_{L}}.

عندما تكون قيمة r أكبر من 4

بالنسبة للخريطة اللوجستية مع r = 4.5، فإن المسارات التي تبدأ من أي نقطة تقريبًا في [0، 1] تتجه نحو اللانهاية السالبة.

عندما تتجاوز المعلمة r القيمة 4، فإن الرأس r /4 من الرسم البياني للخريطة اللوجستية يتجاوز 1. وبقدر ما يخترق الرسم البياني 1، يمكن للمسارات أن تهرب من [0، 1] .

يُعدّ التشعب عند r = 4 نوعًا من الأزمات، وتحديدًا أزمة حدودية. في هذه الحالة، يصبح الجاذب عند [0, 1] غير مستقر وينهار، وبما أنه لا يوجد جاذب خارجه، فإن المسار يتباعد إلى ما لا نهاية.

من جهة أخرى، توجد مدارات تبقى ضمن النطاق [0، 1] حتى لو كانت قيمة r أكبر من 4. ومن الأمثلة الواضحة على ذلك النقاط الثابتة والنقاط الدورية في النطاق [0، 1]، والتي تبقى ضمن هذا النطاق. مع ذلك، توجد أيضًا مدارات أخرى تبقى ضمن النطاق [0، 1] غير النقاط الثابتة والنقاط الدورية.

يتركأ0{\displaystyle A_{0}}لتكن الفترة التي يكون فيها x بحيث f(x) > 1. كما ذكرنا سابقًا، بمجرد أن يكون المتغيرxن{\displaystyle x_{n}}يدخلأ0{\displaystyle A_{0}}، يتباعد إلى سالب ما لا نهاية. يوجد أيضًارن{\displaystyle r_{n}}x في [0, 1] الذي يُقابلأ0{\displaystyle A_{0}}بعد تطبيق الخريطة مرة واحدة. يتم تقسيم هذه الفترة x إلى قسمين، يُطلق عليهما مجتمعين اسمأ1{\displaystyle A_{1}}وبالمثل، هناك أربع فترات زمنية تتوافق معأ1{\displaystyle A_{1}}بعد تطبيق واحد للخريطة، والتي تسمى مجتمعةأ2{\displaystyle A_{2}}وبالمثل، توجد فترات 2nأن{\displaystyle A_{n}}ذلك المدىأ0{\displaystyle A_{0}}بعد n تكرارًا. لذلك، فإن الفترةΛ{\displaystyle \Lambda }تم الحصول عليها عن طريق إزالةأن{\displaystyle A_{n}}من [0، 1] عدد لا نهائي من المرات كما يلي هي مجموعة من المدارات التي تبقى في I.

عملية الإزالةأن{\displaystyle A_{n}}من [0، 1] يشبه بناء مجموعة كانتور المذكورة أعلاه، وفي الواقع توجد Λ في [0، 1] كمجموعة كانتور (مجموعة فرعية مغلقة ومنفصلة تمامًا وكاملة من [0، 1]). علاوة على ذلك، علىΛ{\displaystyle \Lambda }الخريطة اللوجستيةور>4{\displaystyle f_{r>4}}فوضوي.

عندما تكون قيمة r أقل من الصفر

نظرًا لأن الخريطة اللوجستية دُرست كثيرًا كنموذج بيئي، نادرًا ما نوقشت حالة كون المعامل r سالبًا. عندما تتناقص قيمة a من الصفر، حيث −1 < r < 0، تقترب الخريطة تقاربًا مقاربًا من نقطة ثابتة مستقرة xf = 0، ولكن عندما تتجاوز a قيمة −1، تتفرع إلى نقطتين دوريتين، وكما هو الحال في القيم الموجبة، تمر بتشعب مضاعفة الدورة وتصل إلى حالة الفوضى. أخيرًا، عندما تنخفض قيمة a إلى أقل من −2، تتباعد الخريطة إلى ما لا نهاية.

مخطط مداري للمعامل r من -2 إلى 4. يتباعد المدار عندما يتجاوز المعامل a هذا النطاق، سواء على الجانبين السالب أو الموجب.

حلول دقيقة للحالات الخاصة

بالنسبة لخريطة لوجستية ذات معلمات محددةأ{\displaystyle a}، حل دقيق يتضمن الوقت بشكل صريحن{\displaystyle n}والقيمة الأوليةx0{\displaystyle x_{0}}تم الحصول عليها على النحو التالي.

عندما r = 4

عندما r = 2

عندما r = −2

بالنظر إلى الحلول الثلاثة الدقيقة المذكورة أعلاه، فإن جميعها

الفوضى والخريطة اللوجستية

مخطط شبكي للخريطة اللوجستية، يوضح السلوك الفوضوي لمعظم قيم r > 3.57
الدالة اللوجستية f ( باللون الأزرق ) ونسخها المتكررة f₂ و f₃ و f₄ و f₅ عندما r = 3.5 . على سبيل المثال ، لأي قيمة ابتدائية على المحور الأفقي، تعطي f₄ قيمة التكرار بعد أربع دورات.

إن بساطة الخريطة اللوجستية تجعلها نقطة انطلاق شائعة الاستخدام لدراسة مفهوم الفوضى. ويمكن وصف الفوضى بشكل عام بأنها سمة من سمات الأنظمة الفوضوية: [ ديفاني 1989 3 ] (انظر: الديناميكيات الفوضوية ).

هذه خصائص الخريطة اللوجستية لمعظم قيم r بين 3.57 و4 تقريبًا (كما ذُكر أعلاه). [ ماي، روبرت م. (1976) 1 ] أحد الأسباب الشائعة لهذه الحساسية للشروط الابتدائية هو أن الخريطة تمثل عملية طي وتمديد متكررة للفضاء الذي تُعرَّف عليه. في حالة الخريطة اللوجستية، يمكن اعتبار معادلة الفرق التربيعية التي تصفها بمثابة عملية تمديد وطيّ على الفترة (0،1) . [ 4 ]

يوضح الشكل التالي عملية التمديد والطي عبر سلسلة من تكرارات الخريطة. يُظهر الشكل (أ)، على اليسار، مخطط بوانكاريه ثنائي الأبعاد لمساحة حالة الخريطة اللوجستية عند r = 4 ، ويُظهر بوضوح المنحنى التربيعي لمعادلة الفرق ( 1 ). مع ذلك، يُمكننا تضمين التسلسل نفسه في مساحة حالة ثلاثية الأبعاد، لدراسة البنية الأعمق للخريطة. يُوضح الشكل (ب) ذلك، مُبينًا كيف تبدأ النقاط المتقاربة في البداية بالتباعد، لا سيما في مناطق x t التي تُقابل الأجزاء الأكثر انحدارًا من المخطط.

تُظهر مخططات بوانكاريه ثنائية وثلاثية الأبعاد بنية التمدد والطي للخريطة اللوجستية
تُظهر مخططات بوانكاريه ثنائية وثلاثية الأبعاد بنية التمدد والطي للخريطة اللوجستية

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

عند r = 2، الدالةرx(1-x){\displaystyle rx(1-x)}يتقاطعy=x{\displaystyle y=x}عند النقطة القصوى تحديدًا، يكون التقارب إلى نقطة التوازن من رتبةدلتا2ن{\displaystyle \delta ^{2^{n}}}وبالتالي، تُسمى نقطة التوازن "نقطة فائقة الاستقرار". ومعامل ليابونوف الخاص بها هو-{\displaystyle -\infty }. وتُظهر حجة مماثلة وجود حالة فائقة الاستقرارر{\displaystyle r}القيمة ضمن كل فترة زمنية يكون فيها للنظام الديناميكي دورة مستقرة. ويمكن ملاحظة ذلك في مخطط أس ليابونوف على شكل انخفاضات حادة. [ 16 ]

بما أن الخريطة محصورة في فترة على خط الأعداد الحقيقية، فإن بُعدها أقل من أو يساوي واحدًا. تُعطي التقديرات العددية بُعد ارتباط قدره0.500 ± 0.005 ( جراسبرجر ، 1983)، وبُعد هاوسدورف حوالي 0.538 ( جراسبرجر ، 1981)، وبُعد معلومات يُقارب 0.5170976 ( جراسبرجر، 1983) لـ r ≈ 3.5699456 (بداية الفوضى). ملاحظة: يمكن إثبات أن بُعد الارتباط يقع بالتأكيد بين 0.4926 و0.5024.

مع ذلك، غالبًا ما يكون من الممكن تقديم بيانات دقيقة وموثوقة حول احتمالية حالة مستقبلية في نظام فوضوي. إذا كان لنظام ديناميكي (ربما فوضوي) جاذب ، فإنه يوجد مقياس احتمالي يُعطي النسبة المئوية للوقت الذي يقضيه النظام على المدى الطويل في مختلف مناطق الجاذب. في حالة الخريطة اللوجستية ذات المعامل r = 4 وحالة ابتدائية في الفترة (0,1) ، يكون الجاذب أيضًا هو الفترة (0,1) ، ويتوافق المقياس الاحتمالي مع توزيع بيتا ذي المعاملين a = 0.5 و b = 0.5 . تحديدًا، [ 17 ] يكون المقياس الثابت هو

1πx(1-x).{\displaystyle {\frac {1}{\pi {\sqrt {x(1-x)}}}}.}

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

خريطة لوجستية باستخدام دالة أس ليابونوف

التمثيل البياني

يمكن تصور مخطط التشعب للخريطة اللوجستية باستخدام كود بايثون التالي :

استيراد numpy كـ npاستورد مكتبة matplotlib.pyplot باسم pltالفاصل الزمني = ( 2.8 ، 4 ) للبداية، النهايةالدقة = 0.0001عدد التكرارات = 600 عدد التكراراتعدد المخططات = 200lims = np.zeros ( reps )fig , biax = plt.subplots ( )fig.set_size_inches ( 16 , 9 )lims [ 0 ] = np.random.rand ( )for r in np.arange ( interval [ 0 ] , interval [ 1 ], accuracy ) :for i in range ( reps - 1 ):lims [ i + 1 ] = r * lims [ i ] * ( 1 - lims [ i ])biax.plot ([ r ] * numtoplot , lims [ reps - numtoplot : ], "b. " , markersize = 0.02 )biax.set ( xlabel = "r" , ylabel = " x" , title = "logistic map " )plt.show ( )

حالات خاصة من الخريطة

الحد الأعلى عندما يكون 0 ≤ r ≤ 1

على الرغم من أن الحلول الدقيقة للعلاقة التكرارية لا تتوفر إلا في عدد قليل من الحالات، إلا أنه من المعروف وجود حد أعلى مغلق الصيغة للدالة اللوجستية عندما يكون 0 r ≤ 1. [ 18 ] هناك جانبان لسلوك الدالة اللوجستية ينبغي أن يغطيهما الحد الأعلى في هذا النطاق: التضاؤل ​​الهندسي التقاربي مع ثبات r، والتضاؤل ​​الأولي السريع عندما تكون x₀ قريبة من 1، والناتج عن الحد (1 − xₙ) في العلاقة التكرارية . ويغطي الحد التالي كلا هذين التأثيرين:

ن{0،1،...}وx0،ر[0،1]،xنx0ر-ن+x0ن.{\displaystyle \forall n\in \{0,1,\ldots \}\quad {\text{and}}\quad x_{0},r\in [0,1],\quad x_{n}\leq {\frac {x_{0}}{r^{-n}+x_{0}n}}.}

الحل عندما r = 4

يمكن في الواقع حلّ الحالة الخاصة عندما r = 4 بدقة، وكذلك الحالة عندما r = 2 ؛ [ 19 ] أما الحالة العامة فلا يمكن التنبؤ بها إلا إحصائيًا. [ 20 ] الحل عندما r = 4 هو: [ 19 ] [ 21 ]

xن=الخطيئة2(2نθπ)،{\displaystyle x_{n}=\sin ^{2}\left(2^{n}\theta \pi \right),}

حيث تُعطى معلمة الشرط الابتدائي θ بالعلاقة التالية

θ=1πالخطيئة-1(x0).{\displaystyle \theta ={\tfrac {1}{\pi }}\sin ^{-1}\left({\sqrt {x_{0}}}\right).}

بالنسبة لقيم θ النسبية ، بعد عدد محدود من التكرارات، تتحول x<sub> n </sub> إلى متتالية دورية. لكن معظم قيم θ غير نسبية، وبالنسبة لقيم θ غير النسبية ، لا تتكرر x<sub> n </sub> أبدًا - فهي غير دورية. توضح معادلة الحل هذه بوضوح السمتين الرئيسيتين للفوضى - التمدد والطي: يُظهر العامل 2<sup> n</sup> النمو الأسي للتمدد، مما ينتج عنه اعتماد حساس على الشروط الابتدائية ، بينما تحافظ دالة الجيب التربيعي على x<sub> n </sub> مطوية ضمن النطاق [0,1] .

بالنسبة لـ r = 4، يكون الحل المكافئ بدلالة الأعداد المركبة بدلاً من الدوال المثلثية هو [ 19 ]

xن=-α2ن-α-2ن+24{\displaystyle x_{n}={\frac {-\alpha ^{2^{n}}-\alpha ^{-2^{n}}+2}{4}}}

حيث α هو أحد الأعداد المركبة

α=1-2x0±(1-2x0)2-1{\displaystyle \alpha =1-2x_{0}\pm {\sqrt {\left(1-2x_{0}\right)^{2}-1}}}

بمعامل يساوي 1. وكما أن دالة الجيب التربيعي في الحل المثلثي لا تؤدي إلى انكماش أو تمدد مجموعة النقاط التي تمت زيارتها، فإن هذا التأثير في الحل الأخير يتم تحقيقه بواسطة المعامل الواحد لـ α .

وعلى النقيض من ذلك، فإن الحل عندما يكون r = 2 هو [ 19 ]

xن=12-12(1-2x0)2ن{\displaystyle x_{n}={\tfrac {1}{2}}-{\tfrac {1}{2}}\left(1-2x_{0}\right)^{2^{n}}}

لـ x 0 ∈ [0,1) . بما أن (1 − 2 x 0 ) ∈ (−1,1) لأي قيمة لـ x 0 بخلاف النقطة الثابتة غير المستقرة 0، فإن الحد (1 − 2 x 0 ) 2 n يؤول إلى 0 عندما يؤول n إلى اللانهاية، لذا فإن x n يؤول إلى النقطة الثابتة المستقرة 1 / 2 .

إيجاد دورات بأي طول عندما r = 4

في حالة r = 4 ، تكون متتالية التكرارات فوضوية في معظم الحالات الابتدائية. ومع ذلك، يوجد عدد لا نهائي من الحالات الابتدائية التي تؤدي إلى دورات، وبالفعل توجد دورات بطول k لجميع الأعداد الصحيحة k > 0. يمكننا استغلال العلاقة بين الخريطة اللوجستية والتحويل الثنائي (المعروف أيضًا باسم خريطة إزاحة البتات ) لإيجاد دورات بأي طول. إذا كان x يتبع الخريطة اللوجستية x <sub>n +1</sub> = 4x <sub> n</sub> (1 − x<sub> n</sub> ) و y يتبع التحويل الثنائي

yن+1={2yن0yن<122yن-112yن<1،{\displaystyle y_{n+1}={\begin{cases}2y_{n}&0\leq y_{n}<{\tfrac {1}{2}}\\2y_{n}-1&{\tfrac {1}{2}}\leq y_{n}<1,\end{cases}}}

ثم يرتبط الاثنان بتماثل شكلي

xن=الخطيئة2(2πyن).{\displaystyle x_{n}=\sin ^{2}\left(2\pi y_{n}\right).}

يُطلق على التحويل الثنائي أيضًا اسم خريطة إزاحة البتات، لأنه عند كتابة y بالترميز الثنائي، تُحرك الخريطة النقطة الثنائية خانة واحدة إلى اليمين (وإذا أصبحت البتة الموجودة على يسار النقطة الثنائية "1"، فإن هذه "1" تُغير إلى "0"). على سبيل المثال، تحدث دورة بطول 3 إذا كان للتكرار تسلسل متكرر من 3 بتات في تمثيله الثنائي (وهو ليس أيضًا تسلسلًا متكررًا من بت واحد): 001، 010، 100، 110، 101، أو 011. يُحوّل التكرار 001001001... إلى 010010010...، والذي يُحوّل بدوره إلى 100100100...، والذي يُحوّل بدوره إلى 001001001... الأصلي؛ لذا فهذه دورة من 3 لخريطة إزاحة البتات. أما التسلسلات المتكررة الثلاثة الأخرى في التوسع الثنائي فتعطي الدورة الثلاثية 110110110... → 101101101... → 011011011... → 110110110.... يمكن تحويل أي من هذه الدورات الثلاثية إلى صيغة كسرية: على سبيل المثال، يمكن كتابة الدورة الثلاثية الأولى المعطاة على النحو التالي : 1/7 2/7 4/7 1/7 . باستخدام التحويل أعلاه من خريطة إزاحة البت إلىر=4{\displaystyle r=4}تُعطي الخريطة اللوجستية الدورة اللوجستية المقابلة 0.611260467... → 0.950484434... → 0.188255099... → 0.611260467.... وبالمثل، يُمكننا ترجمة دورة الإزاحة الثنائية الثلاثية الأخرى إلى دورتها اللوجستية المقابلة. كما يُمكن إيجاد دورات بأي طول k في خريطة الإزاحة الثنائية، ثم ترجمتها إلى الدورات اللوجستية المقابلة.

مع ذلك، ولأن معظم الأعداد في الفترة [0,1) غير نسبية، فإن معظم الشروط الابتدائية لخريطة إزاحة البتات تؤدي إلى عدم دورية الفوضى. هذه إحدى الطرق لإثبات أن خريطة لوجستية r = 4 فوضوية في معظم الشروط الابتدائية.

عدد الدورات ذات الطول الأدنى k = 1، 2، 3، ... للخريطة اللوجستية ذات r = 4 ( خريطة الخيمة ذات μ = 2 ) هو متتالية أعداد صحيحة معروفة (المتتالية A001037 في OEIS ) : 2، 1، 2، 3، 6، 9، 18، 30، 56، 99، 186، 335، 630، 1161 .... هذا يعني أن الخريطة اللوجستية ذات r = 4 تحتوي على نقطتين ثابتتين، ودورة واحدة بطول 2، ودورتين بطول 3، وهكذا . تأخذ هذه المتتالية شكلاً بسيطاً للغاية للعدد الأولي k : 2 ⋅ 2k − 11 / k . على سبيل المثال: 2  2 13 − 1 − 1 / 13 = 630 هو عدد الدورات التي طولها 13. وبما أن هذه الحالة من الخريطة اللوجستية فوضوية بالنسبة لجميع الشروط الأولية تقريبًا، فإن جميع هذه الدورات ذات الطول المحدود غير مستقرة. 

عالمية

فئة من التطبيقات التي تُظهر سلوكًا متجانسًا

رسم بياني لخريطة الجيب (4-1)
مخطط مداري لخريطة الجيب (4-1)

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

هنا، المجال هو 0 ≤ b ≤ 1 و 0 ≤ x ≤ 1. تُظهر خريطة الجيب (4-1) سلوكًا مطابقًا نوعيًا لخريطة اللوجستي (1-2) : مثل خريطة اللوجستي، تصبح فوضوية أيضًا عبر مسار مضاعفة الدورة مع زيادة المعامل b، وعلاوة على ذلك، مثل خريطة اللوجستي، تُظهر أيضًا نافذة في المنطقة الفوضوية .

كل من الخريطة اللوجستية وخريطة الجيب عبارة عن خرائط أحادية البعد تقوم بتحويل الفترة [0، 1] إلى [0، 1] وتفي بالخاصية التالية، والتي تسمى أحادية النمط .

و(0)=و(1)=0{\displaystyle f(0)=f(1)=0}الدالة قابلة للتفاضل، وتوجد نقطة حرجة وحيدة c في الفترة [0, 1] بحيثو(ج)=0{\displaystyle f'(c)=0}بشكل عام، إذا كانت الخريطة أحادية البعد ذات المعلمة الواحدة والمتغير الواحد أحادية النمط ويمكن تقريب الرأس بواسطة متعدد حدود من الدرجة الثانية، فعندئذٍ، بغض النظر عن الشكل المحدد للخريطة، سيحدث تسلسل لانهائي من التشعبات يضاعف الفترة لنطاق المعلمة 3 ≤ r ≤ 3.56994...، وتكون النسبة δ المعرفة بالمعادلة (3-13) مساوية لثابت فيجنباوم، 4.669 ....

يُعد نمط المدارات الدورية المستقرة التي تظهر من الخريطة اللوجستية نمطًا عالميًا أيضًا . بالنسبة للخريطة أحادية النمط،xن+1=جو(xن){\displaystyle x_{n+1}=cf(x_{n})}مع المعامل c، تستمر المدارات الدورية المستقرة ذات الفترات المختلفة في الظهور ضمن مجال معامل تكون فيه النقطتان الثابتتان غير مستقرتين، ومن المعروف أن نمط ظهورها (عدد المدارات الدورية المستقرة ذات فترة معينة وترتيب ظهورها) نمط مشترك . بعبارة أخرى، بالنسبة لهذا النوع من الخرائط، يكون تسلسل المدارات الدورية المستقرة هو نفسه بغض النظر عن الشكل المحدد للخريطة . بالنسبة للخريطة اللوجستية، يكون مجال المعامل 3 < a < 4، أما بالنسبة لخريطة الجيب (4-1)، فإن مجال المعامل للتسلسل المشترك للمدارات الدورية المستقرة هو 0.71... < b < 1. يُطلق على هذا التسلسل العالمي للمدارات الدورية المستقرة اسم تسلسل U.

بالإضافة إلى ذلك، تتميز الدالة اللوجستية بخاصية أن مشتقتها الشوارزية سالبة دائمًا على الفترة [0، 1]. المشتقة الشوارزية للدالة f (من الفئة C3) هي

في الواقع، عندما نحسب مشتقة شوارز للدالة اللوجستية، نحصل على

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

رسم الخرائط المترافقة الطوبولوجية

لنرمز بالرمز ∘ إلى تركيب التطبيقات. بشكل عام، بالنسبة للفضاء الطوبولوجي X و Y، فإن التطبيقين f : X → X و g : Y → Y يتكونان من تشاكل طوبولوجي h  : X → Y.

يُقال إن الدالتين f و g مترافقتان طوريًا إذا كانتا تحققان العلاقة التالية : f = g. يلعب مفهوم الترافق الطوري دورًا هامًا في دراسة الأنظمة الديناميكية . تُظهر الدالتان المترافقتان طوريًا f و g سلوكًا متطابقًا جوهريًا، وإذا كان سلوك f دوريًا، فإن g يكون دوريًا أيضًا، وإذا كان سلوك f فوضويًا، فإن g يكون فوضويًا أيضًا .

على وجه الخصوص، إذا كان التشاكل h خطيًا، فإن الدالتين f و g تُقالان مترافقتان خطيًا. كل دالة تربيعية مترافقة خطيًا مع كل دالة تربيعية أخرى. ومن ثم،

هي مترافقة خطية للخريطة اللوجستية لأي قيمة للمعامل a . تُسمى المعادلتان (4-6) و(4-7) أيضًا بالخرائط اللوجستية . وعلى وجه الخصوص، يُعد الشكل (4-7) مناسبًا للحسابات العددية التي تستغرق وقتًا طويلاً، لأنه يتطلب جهدًا حسابيًا أقل .

منظر مداري لخريطة الخيمة (4-8). تربطها علاقة طوبولوجية مترافقة مع الخريطة اللوجستية a = 4.

علاوة على ذلك، الخريطة اللوجستيةوأ=4{\displaystyle f_{a=4}}لر=4{\displaystyle r=4}وهي مترافقة طوبولوجيًا مع خريطة الخيمة التالية T ( x ) وخريطة إزاحة برنولي B ( x ) .

يمكن استخدام علاقات الاقتران الطوري هذه لإثبات أن الخريطة اللوجستيةوأ=4{\displaystyle f_{a=4}}هي فوضوية تمامًا، ولإيجاد الحل الدقيق (3-19) لـور=4{\displaystyle f_{r=4}} .

بدلاً من ذلك، وبتقديم مفهوم الأنظمة الديناميكية الرمزية، ضع في اعتبارك خريطة الإزاحة التالية σ المعرفة على فضاء السلسلة الرمزية المكون من سلاسل من 0 و1 كما هو موضح أعلاه :

هنا،sأنا{\displaystyle s_{i}}تكون قيمته 0 أو 1. على المجموعةΛ{\displaystyle \Lambda }الخريطة اللوجستية التي تم تقديمها في المعادلة (3-18)ور>4{\displaystyle f_{r>4}}وهي مترافقة طوبولوجيًا مع خريطة الإزاحة، لذا يمكننا استخدام هذا لاستنتاج ذلك.ور>4{\displaystyle f_{r>4}}علىΛ{\displaystyle \Lambda }فوضوي .

مسار مضاعفة الدورة نحو الفوضى

في الخريطة اللوجستية، لدينا دالةور(x)=رx(1-x){\displaystyle f_{r}(x)=rx(1-x)}ونريد دراسة ما يحدث عند تكرار عملية رسم الخريطة عدة مرات. قد تستقر الخريطة عند نقطة ثابتة، أو دورة ثابتة، أو حالة من الفوضى. عندما تستقر الخريطة في دورة ثابتة مستقرة ذات طول معين، فإنها تستقر عند نقطة ثابتة مستقرة ذات طول معين.ن{\displaystyle n}، سنجد أن الرسم البياني لـورن{\displaystyle f_{r}^{n}}والرسم البياني لـxx{\displaystyle x\mapsto x}يتقاطع عندن{\displaystyle n}النقاط، وميل الرسم البياني لـورن{\displaystyle f_{r}^{n}}محصور في(-1،+1){\displaystyle (-1,+1)}عند تلك التقاطعات.

على سبيل المثال، عندمار=3.0{\displaystyle r=3.0}لدينا تقاطع واحد، بميل محدود في(-1،+1){\displaystyle (-1,+1)}مما يشير إلى أنها نقطة ثابتة مفردة مستقرة.

مثلر{\displaystyle r}زيادة إلى ما هو أبعدر=3.0{\displaystyle r=3.0}، تنقسم نقطة التقاطع إلى نقطتين، وهو ما يُعرف بمضاعفة الدورة. على سبيل المثال، عندمار=3.4{\displaystyle r=3.4}توجد ثلاث نقاط تقاطع، النقطة الوسطى غير مستقرة، والنقطتان الأخريان مستقرتان.

مثلر{\displaystyle r}الأساليبر=3.45{\displaystyle r=3.45}ثم يحدث تضاعف آخر للدورة بنفس الطريقة. وتتكرر مضاعفات الدورة بشكل متزايد، حتى تصل إلى حد معين.ر3.56994567{\displaystyle r\approx 3.56994567}تصبح مضاعفات الفترة لانهائية، وتصبح الخريطة فوضوية. هذا هو مسار مضاعفة الفترة إلى الفوضى .

العلاقة بينxن+2{\displaystyle x_{n+2}}وxن{\displaystyle x_{n}}متىأ=2.7{\displaystyle a=2.7}قبل حدوث تشعب مضاعفة الفترة، يتقارب المدار نحو النقطة الثابتة.xو2{\displaystyle x_{f2}}.
العلاقة بينxن+2{\displaystyle x_{n+2}}وxن{\displaystyle x_{n}}متىأ=3{\displaystyle a=3}ميل المماس عند النقطة الثابتةxو2{\displaystyle x_{f2}}. يساوي 1 بالضبط، ويحدث تشعب مضاعفة الفترة.
العلاقة بينxن+2{\displaystyle x_{n+2}}وxن{\displaystyle x_{n}}متىأ=3.3{\displaystyle a=3.3}النقطة الثابتةxو2{\displaystyle x_{f2}}يصبح غير مستقر، وينقسم إلى دورة مستقرة دورية-2.
متىر=3.0{\displaystyle r=3.0}لدينا تقاطع واحد، بميل يساوي تمامًا+1{\displaystyle +1}مما يشير إلى أنها على وشك أن تخضع لمضاعفة دورية.
متىر=3.4{\displaystyle r=3.4}توجد ثلاث نقاط تقاطع، النقطة الوسطى غير مستقرة، والنقطتان الأخريان مستقرتان.
متىر=3.45{\displaystyle r=3.45}توجد ثلاث نقاط تقاطع، النقطة الوسطى منها غير مستقرة، أما النقطتان الأخريان فلهما ميل ثابت+1{\displaystyle +1}مما يشير إلى أنها على وشك أن تخضع لمضاعفة دورية أخرى.
متىر3.56994567{\displaystyle r\approx 3.56994567}، هناك عدد لا نهائي من التقاطعات، وقد وصلنا إلى الفوضى عبر مسار مضاعفة الفترة .

حد التوسع

الاقتراب من حد التوسع كـر{\displaystyle r}الأساليبر*=3.5699{\displaystyle r^{*}=3.5699\cdots }من الأسفل.
عند نقطة الفوضىر*=3.5699{\displaystyle r^{*}=3.5699\cdots }بينما نكرر مضاعفة الفترةور*1،ور*2،ور*4،ور*8،ور*16،...{\displaystyle f_{r^{*}}^{1},f_{r^{*}}^{2},f_{r^{*}}^{4},f_{r^{*}}^{8},f_{r^{*}}^{16},\dots }، تبدو الرسوم البيانية متشابهة، باستثناء أنها تتقلص نحو المنتصف، وتدور بمقدار 180 درجة، وتتقارب إلى شكل كسري.

بالنظر إلى الصور، يمكن للمرء أن يلاحظ أنه عند نقطة الفوضىر*=3.5699{\displaystyle r^{*}=3.5699\cdots }منحنىور*{\displaystyle f_{r^{*}}^{\infty }}يبدو كشكل كسري. علاوة على ذلك، مع تكرار مضاعفة الدوراتور*1،ور*2،ور*4،ور*8،ور*16،...{\displaystyle f_{r^{*}}^{1},f_{r^{*}}^{2},f_{r^{*}}^{4},f_{r^{*}}^{8},f_{r^{*}}^{16},\dots }، تبدو الرسوم البيانية متشابهة، باستثناء أنها تتقلص باتجاه المنتصف، وتدور بمقدار 180 درجة.

يشير هذا إلى وجود حدٍّ للتوسع: إذا ضاعفنا الدالة بشكل متكرر، ثم قمنا بتوسيعها بمقدارα{\displaystyle \alpha }لثابت معينα{\displaystyle \alpha }:و(x)-αو(و(-x/α)){\displaystyle f(x)\mapsto -\alpha f(f(-x/\alpha ))}ثم عند الحد، سنحصل على دالةز{\displaystyle g}ذلك يرضيز(x)=-αز(ز(-x/α)){\displaystyle g(x)=-\alpha g(g(-x/\alpha ))}هذه دالة فيجنباوم ، التي تظهر في معظم مسارات مضاعفة الدورة المؤدية إلى الفوضى (وبالتالي فهي مثال على الشمولية ). علاوة على ذلك، مع تناقص فترات مضاعفة الدورة، تتقارب النسبة بين فترتي مضاعفة دورة متتاليتين إلى حدٍّ معين، وهو ثابت فيجنباوم الأول.دلتا=4.6692016{\displaystyle \delta =4.6692016\cdots }.

بالنسبة للقيم الخاطئة لعامل القياسα{\displaystyle \alpha }لا تتقارب الخريطة إلى حد معين، ولكن عندماα=2.5029...{\displaystyle \alpha =2.5029\dots }، إنها تتقارب.
عند نقطة الفوضىر*=3.5699{\displaystyle r^{*}=3.5699\cdots }، بينما نكرر تكرار المعادلة الوظيفيةو(x)-αو(و(-x/α)){\displaystyle f(x)\mapsto -\alpha f(f(-x/\alpha ))}معα=2.5029...{\displaystyle \alpha =2.5029\dots }، نجد أن الخريطة تتقارب إلى حد معين.

الثابتα{\displaystyle \alpha }يمكن إيجادها عدديًا بتجربة العديد من القيم الممكنة. بالنسبة للقيم الخاطئة، لا تتقارب الدالة إلى حد معين، ولكن عندما تكونα=2.5029...{\displaystyle \alpha =2.5029\dots }، يتقارب. هذا هو ثابت فيجنباوم الثاني.

نظام فوضوي

في ظل النظام الفوضوي،ور{\displaystyle f_{r}^{\infty }}، يصبح حد تكرارات الخريطة عبارة عن نطاقات داكنة فوضوية تتخللها نطاقات ساطعة غير فوضوية.

في ظل النظام الفوضوي،ور{\displaystyle f_{r}^{\infty }}، يصبح حد تكرارات الخريطة عبارة عن نطاقات داكنة فوضوية تتخللها نطاقات ساطعة غير فوضوية.

حدود أخرى للتوسع

متىر{\displaystyle r}الأساليبر3.8494344{\displaystyle r\approx 3.8494344}لدينا نهج آخر لمضاعفة الدورات في التعامل مع الفوضى، ولكن هذه المرة مع دورات 3، 6، 12، ... وهذا النهج له نفس ثوابت فيجنباوم.دلتا،α{\displaystyle \delta ,\alpha }حدو(x)-αو(و(-x/α)){\textstyle f(x)\mapsto -\alpha f(f(-x/\alpha ))}وهي أيضاً نفس دالة فيجنباوم . هذا مثال على العمومية .

خريطة لوجستية تقترب من حد التوسع الفوضوي المضاعف للدورةر*=3.84943...{\displaystyle r^{*}=3.84943\dots }من الأسفل. عند الحد، يكون لهذا الشكل نفس شكل شكلر*=3.5699{\displaystyle r^{*}=3.5699\cdots }، لأن جميع مسارات مضاعفة الفترة المؤدية إلى الفوضى هي نفسها (الشمولية).

يمكننا أيضًا النظر في مسار مضاعفة الفترة ثلاث مرات نحو الفوضى عن طريق اختيار سلسلة منر1،ر2،...{\displaystyle r_{1},r_{2},\dots }بحيثرن{\displaystyle r_{n}}وهي أدنى قيمة في الفترة-3ن{\displaystyle 3^{n}}نافذة مخطط التشعب. على سبيل المثال، لدينار1=3.8284،ر2=3.85361،...{\displaystyle r_{1}=3.8284,r_{2}=3.85361,\dots }، مع الحدر=3.854077963...{\displaystyle r_{\infty }=3.854077963\dots }. هذا له زوج مختلف من ثوابت فيجنباومدلتا=55.26...،α=9.277...{\displaystyle \delta =55.26\dots ,\alpha =9.277\dots }[ 22 ] وور{\displaystyle f_{r}^{\infty }}يتقارب إلى النقطة الثابتةو(x)-αو(و(و(-x/α))){\displaystyle f(x)\mapsto -\alpha f(f(f(-x/\alpha )))}كمثال آخر، يمتلك التكرار الدوري من الدرجة الرابعة زوجًا من ثوابت فيجنباوم يختلف عن زوج ثوابت التكرار الدوري من الدرجة الثانية، على الرغم من أن التكرار الدوري من الدرجة الرابعة يتحقق من خلال تكرارين دوريين. عرّف بالتفصيلر1،ر2،...{\displaystyle r_{1},r_{2},\dots }بحيثرن{\displaystyle r_{n}}وهي أدنى قيمة في الفترة-4ن{\displaystyle 4^{n}}نافذة مخطط التشعب. ثم لدينار1=3.960102،ر2=3.9615554،...{\displaystyle r_{1}=3.960102,r_{2}=3.9615554,\dots }، مع الحدر=3.96155658717...{\displaystyle r_{\infty }=3.96155658717\dots }. هذا له زوج مختلف من ثوابت فيجنباومدلتا=981.6...،α=38.82...{\displaystyle \delta =981.6\dots ,\alpha =38.82\dots }.

بشكل عام، لكل مسار مضاعف للدورة يؤدي إلى الفوضى زوج خاص به من ثوابت فيجنباوم. في الواقع، يوجد عادةً أكثر من زوج واحد. على سبيل المثال، بالنسبة لدورة 7-pling، يوجد ما لا يقل عن 9 أزواج مختلفة من ثوابت فيجنباوم. [ 22 ]

عمومًا،3دلتا2α2{\textstyle 3\delta \approx 2\alpha ^{2}}وتصبح العلاقة دقيقة عندما يتزايد كلا العددين إلى ما لا نهاية:ليمدلتا/α2=2/3{\displaystyle \lim \delta /\alpha ^{2}=2/3}.

شمولية فيجنباوم للخرائط أحادية البعد

عمومية الخرائط أحادية البعد ذات القيم القصوى المكافئة وثوابت فيجنباومدلتا=4.669201...{\displaystyle \delta =4.669201...}،α=2.502907...{\displaystyle \alpha =2.502907...}[ 23 ] [ 24 ]

الزيادة التدريجية لـجي{\displaystyle G}على فترات[0،){\displaystyle [0,\infty )}تغير الديناميكيات من منتظمة إلى فوضوية [ 25 ] بنفس مخطط التشعب النوعي مثل تلك الخاصة بالخريطة اللوجستية.

تقدير إعادة التطبيع

يمكن تقدير ثوابت فيجنباوم عن طريق حجة إعادة التطبيع. (القسم 10.7، [ 16 ] ).

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

عرّف العائلةور(x)=-(1+ر)x+x2{\displaystyle f_{r}(x)=-(1+r)x+x^{2}}للعائلة نقطة توازن عند الصفر، ور{\displaystyle r}عند ازديادها، تخضع لتشعب مضاعفة الفترة عندر=ر0،ر1،ر2،...{\displaystyle r=r_{0},r_{1},r_{2},...}.

يحدث التشعب الأول عندر=ر0=0{\displaystyle r=r_{0}=0}بعد حدوث تشعب مضاعفة الدورة، يمكننا إيجاد المدار المستقر ذي الدورة 2 عن طريقور(ص)=q،ور(q)=ص{\displaystyle f_{r}(p)=q,f_{r}(q)=p}، مما ينتج عنه{ص=12(ر+ر(ر+4))q=12(ر-ر(ر+4)){\displaystyle {\begin{cases}p={\frac {1}{2}}(r+{\sqrt {r(r+4)}})\\q={\frac {1}{2}}(r-{\sqrt {r(r+4)}})\end{cases}}}في مرحلة مار=ر1{\displaystyle r=r_{1}}يخضع المدار المستقر ذو الفترة 2 لتشعب مضاعف للفترة مرة أخرى، مما ينتج عنه مدار مستقر ذو فترة 4. لمعرفة شكل المدار المستقر، نقوم بالتكبير حول منطقةx=ص{\displaystyle x=p}باستخدام التحويل الأفينيتي(x)=x/ج+ص{\displaystyle T(x)=x/c+p}والآن، باستخدام الجبر الروتيني، لدينا(تي-1ور2تي)(x)=-(1+S(ر))x+x2+يا(x3){\displaystyle (T^{-1}\circ f_{r}^{2}\circ T)(x)=-(1+S(r))x+x^{2}+O(x^{3})}أينS(ر)=ر2+4ر-2،ج=ر2+4ر-3ر(ر+4){\displaystyle S(r)=r^{2}+4r-2,c=r^{2}+4r-3{\sqrt {r(r+4)}}}في حواليS(ر)=0{\displaystyle S(r)=0}ثم يحدث التشعب الثاني، وبالتاليS(ر1)0{\displaystyle S(r_{1})\approx 0}.

عن طريق التشابه الذاتي، يحدث التشعب الثالث عندماS(ر)ر1{\displaystyle S(r)\approx r_{1}}وهكذا دواليك.رنS(رن+1){\displaystyle r_{n}\approx S(r_{n+1})}، أورن+1رن+6-2{\displaystyle r_{n+1}\approx {\sqrt {r_{n}+6}}-2}بتكرار هذه الخريطة، نجدر=ليمنرنليمنS-ن(0)=12(17-3){\displaystyle r_{\infty }=\lim _{n}r_{n}\approx \lim _{n}S^{-n}(0)={\frac {1}{2}}({\sqrt {17}}-3)}، وليمنر-رنر-رن+1S(ر)1+17{\displaystyle \lim _{n}{\frac {r_{\infty }-r_{n}}{r_{\infty }-r_{n+1}}}\approx S'(r_{\infty })\approx 1+{\sqrt {17}}}.

وبالتالي، لدينا التقديراتدلتا1+17=5.12...{\displaystyle \delta \approx 1+{\sqrt {17}}=5.12...}، وαر2+4ر-3ر2+4ر-2.24...{\displaystyle \alpha \approx r_{\infty }^{2}+4r_{\infty }-3{\sqrt {r_{\infty }^{2}+4r_{\infty }}}\approx -2.24...}هذه القيم تقع ضمن نطاق 10% من القيم الحقيقية.

العلاقة بالمعادلة التفاضلية العادية اللوجستية

تُظهر الخريطة اللوجستية العديد من خصائص الحلول الدورية والفوضوية، بينما تُظهر المعادلة التفاضلية العادية اللوجستية حلولًا منتظمة، تُعرف عادةً بدالة سيجمويد على شكل حرف S. يمكن اعتبار الخريطة اللوجستية النظير المتقطع للمعادلة التفاضلية العادية اللوجستية، وقد نُوقشت علاقتهما باستفاضة في الأدبيات. [ 26 ]

الخريطة اللوجستية كنموذج للسكان البيولوجيين

نموذج السكان المنفصل

على الرغم من أن لورنز استخدم الخريطة اللوجستية في عام 1964، [ 21 ] إلا أنها اكتسبت شعبية واسعة النطاق بفضل أبحاث عالم الأحياء الرياضي البريطاني روبرت ماي، وأصبحت معروفة على نطاق واسع كصيغة لدراسة التغيرات في أعداد الكائنات الحية. في مثل هذه الخريطة اللوجستية لأعداد الكائنات الحية، يكون المتغيرxن{\displaystyle x_{n}}يمثل هذا العدد عدد الكائنات الحية التي تعيش في بيئة معينة (أو بتعبير أدق، حجم الجماعة). علاوة على ذلك، يُفترض عدم مغادرة أي كائنات حية للبيئة وعدم دخول أي كائنات حية خارجية إليها (أو عدم وجود تأثير يُذكر حتى في حالة الهجرة)، والنموذج الرياضي المستخدم لدراسة الزيادة أو النقصان في عدد السكان في مثل هذه الحالة هو الخريطة اللوجستية في علم الأحياء الرياضي.

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

لنفترض أن N يمثل عدد أفراد نوع واحد في بيئة معينة. أبسط نموذج لنمو السكان هو النموذج الذي يستمر فيه نمو السكان بمعدل ثابت بالنسبة لعدد الأفراد. يُسمى هذا النوع من نماذج نمو السكان بالنموذج المالتوسي، ويمكن التعبير عنه كما يلي :

هنا، يُمثل N<sub>n</sub> عدد الأفراد في الجيل n، وα معدل نمو السكان، وهو ثابت موجب . مع ذلك، في النموذج (5-1)، يستمر عدد السكان في النمو بلا حدود، مما يجعله نموذجًا غير واقعي لمعظم الظواهر في العالم الحقيقي . نظرًا لوجود حد أقصى لعدد الأفراد الذين يمكن للبيئة استيعابهم، فمن الطبيعي أن ينخفض ​​معدل النمو α مع ازدياد عدد السكان N<sub>n</sub> . يُطلق على هذا التغير في معدل النمو نتيجة لتغيرات كثافة السكان اسم تأثير الكثافة . معادلة الفرق التالية هي أبسط نموذج مُحسَّن يعكس تأثير الكثافة في النموذج (5-1) .

هنا، يُمثل a أقصى معدل نمو ممكن في البيئة، بينما يُمثل b قوة تأثير كثافة الأفراد. يفترض النموذج (5-2) أن معدل النمو يتناقص ببساطة بما يتناسب مع عدد الأفراد . لنفترض أن N<sub>n</sub> في المعادلة (5-2) هو

بعد إجراء تحويل المتغيرات، يتم اشتقاق الخريطة اللوجستية التالية :

عند استخدام المعادلة (5-2) أو المعادلة (5-4) كحجم لجماعة كائن حي، إذا أصبحت قيمة Nn أو xn سالبة، فإنها تفقد معناها كحجم للجماعة. ولمنع ذلك، يُشترط أن تكون قيمة x0 الابتدائية 0 ≤ x0 ≤ 1، وأن يكون المعامل a 0 ≤ r ≤ 4.

بدلاً من ذلك، يمكننا افتراض حجم سكان أقصى K يمكن للبيئة أن تدعمه، واستخدام هذا لـ

يمكن اشتقاق الخريطة اللوجستية من خلال النظر في معادلة فرقية تتضمن تأثيرات الكثافة بالشكل التالي:xن=شمالن/ك{\displaystyle x_{n}=N_{n}/K} ، حيث المتغيرxن{\displaystyle x_{n}}يمثل نسبة عدد الأفرادشمالن{\displaystyle N_{n}}إلى الحد الأقصى لعدد الأفراد K.

تجزئة المعادلة اللوجستية

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

هنا، يُمثل N عدد الكائنات الحية أو كثافتها السكانية، و t الزمن المتصل، و K و r معلمات. K هي القدرة الاستيعابية، و r هو معدل الزيادة الطبيعية الجوهرية، والذي يكون عادةً موجبًا. يُشير الطرف الأيسر من هذه المعادلة dN/dt إلى معدل تغير حجم السكان عند الزمن t .

مثال على حل المعادلة اللوجستية. بعد مرور الوقت t، يتقارب حجم السكان N إلى القدرة الاستيعابية K بغض النظر عن القيمة الأولية.

تبدو المعادلة اللوجستية (5-6) مشابهة للخريطة اللوجستية (5-4)، لكن سلوك الحلول يختلف تمامًا عن سلوك الخريطة اللوجستية . طالما أن القيمة الابتدائية N₀ موجبة، فإن حجم المجموعة N للمعادلة اللوجستية يتقارب دائمًا بشكل رتيب إلى K.

يمكن اشتقاق الخريطة اللوجستية بتطبيق طريقة أويلر، وهي طريقة لحل المعادلات التفاضلية العادية من الرتبة الأولى عدديًا، على هذه المعادلة اللوجستية. [ملاحظة 2] تستخدم طريقة أويلر فاصلًا زمنيًا (حجم خطوة زمنية) Δt لتقريب معدل النمو dN/dt، والذي يُقارب كما يلي :

يؤدي هذا التقريب إلى الخريطة اللوجستية التالية :

أينxن{\displaystyle x_{n}}و a في هذه المعادلة مرتبطة بالمعلمات والمتغيرات الأصلية وحجم الخطوة الزمنية كما يلي :

إذا كانت قيمة Δt صغيرة بما يكفي، فإن المعادلة (5-8) تُعدّ تقريبًا صحيحًا للمعادلة الأصلية (5-6)، وتتطابق مع حلّها عندما تقترب Δt من الصفر . من جهة أخرى، عندما تكبر قيمة Δt، ينحرف الحلّ عن الحلّ الأصلي . علاوة على ذلك، وبسبب العلاقة في المعادلة (5-10)، فإن زيادة Δt تُكافئ زيادة قيمة المعامل a . وبالتالي، فإن زيادة Δt لا تزيد الخطأ عن المعادلة الأصلية فحسب، بل تُنتج أيضًا سلوكًا فوضويًا في الحلّ .

تحديد الموقع

As described above, in biological population dynamics, the logistic map is one of the models of discrete growth processes. However, unlike the laws of physics, the logistic map as a model of biological population size is not derived from direct experimental results or universally valid principles . Although there is some rationality in the way it is derived, it is essentially a "model" thought up in one's mind . May, who made the logistic map famous, did not claim that the model he was discussing accurately represented the increase and decrease in population size . Historically, continuous-time models based on differential equations have been widely used in the study of biological population dynamics, and the application of these continuous-time models has deepened our understanding of biological population dynamics . As a discrete-time population model that takes into account density effects, the Ricker model, in which the population size is not negative, is more realistic .

Generally speaking, mathematical models can provide important qualitative information about population dynamics, but their results should not be taken too seriously without experimental support. Even if the conclusions of mathematical models deviate from those of biological studies, mathematical modeling is still useful because it can provide a useful control. Biological issues may be raised by reviewing the model construction process and settings, or the biological knowledge and assumptions that the model is based on. Although the logistic map is too simple to be realistic as a population model, its results suggest that a variety of population fluctuations may occur due to the dynamics inherent in the population itself, regardless of random influences from the environment.

Applications

Coupled map system

The degree of freedom or dimension of a one-variable logistic map as a system is one . On the other hand, in the real natural world, it is thought that there are many chaotic systems with many degrees of freedom, not only in time but also in space . Alternatively, the synchronization phenomenon of oscillators performing chaotic motion is also a research subject . To investigate such things, there is a method of coupled maps that couples many difference equations (maps) . The logistic map is often used as a subject of coupled map model research . The reason for this is that the logistic map itself has already been well investigated as a typical model of chaos, and there is an accumulation of research on it .

There are various methods for the specific coupling in the coupled map model . Suppose a total of N maps are coupled, and the state of the i-th map at time n is represented by xn(i){\displaystyle x_{n}(i)} . In a method called globally coupled maps, xn+1(i){\displaystyle x_{n+1}(i)} is formulated as follows :

In the current field of coupled oscillators, the simplest model is the following, in which two oscillators, x and y, are coupled by a difference in variables :

In these equations, f( x ) is the specific map to incorporate into the coupled map model, and applies here if the logistic map is used .

Changes in two variables (top) and their difference (bottom) in the coupled map model ( 6-2 ) with a = 3.8 and D = 0.43 . The two variables suddenly become out of sync after synchronization, and then return to the sync state.

In equations ( 6-1 ) and ( 6-2 ), ε and D are parameters called coupling coefficients, which indicate the strength of the coupling between the maps . On the other hand, when the logistic map is incorporated into a coupled map model, the parameter a of the logistic map indicates the strength of the nonlinearity of the model . By changing the value of a and the value of ε or D, various phenomena appear in the coupled map system of logistic maps. For example, in model ( 6-2 ), when D is increased to a value Dc or more, x and y oscillate chaotically while synchronously . Even below Dc, not only do chaotic oscillations occur in a continuous manner . When D is in a certain range, x and y oscillate with two periods even though r = 4 . When a = 3.8, behavior in which synchronous and asynchronous states alternate continuously is also observed .

In a study of the application of the logistic map to a globally coupled map with a large degree of freedom ( 6-1 ), a phenomenon called chaotic itinerancy was found . This is a phenomenon in which the orbit traverses a region in phase space that is said to be the remains of an attractor, repeating the cycle from an orderly state in which several clusters oscillate together to a disordered state, then to another cluster state, then back to the disordered state again, and so on .

Pseudorandom number generator

In the fields of computer simulation and information security, the creation of pseudorandom numbers using a computer is an important technique, and one of the methods for generating pseudorandom numbers is the use of chaos. Although a pseudorandom number generator based on chaos with sufficient performance has not yet been realized, several methods have been proposed. Several researchers have also investigated the possibility of creating a pseudorandom number generator based on chaos for the logistic map.

Parameter r = 4 is often used for pseudorandom number generation using the logistic map. Historically, as described below, in 1947, shortly after the birth of electronic computers, Stanisław Ulam and John von Neumann also pointed out the possibility of a pseudorandom number generator using the logistic map with r = 4. However, the distribution of points for the logistic map fr=4{\displaystyle f_{r=4}} is as shown in equation ( 3-17 ), and the numbers that are generated are biased toward 0 and 1. Therefore, some processing is required to obtain unbiased uniform random numbers. Methods for doing so include:

A method for converting the obtained values to a uniform distribution using the tent map ( 4-8 ). The resulting number is converted to either 0 or 1 using a threshold, as in the coin tossing analogy above, and this process is repeated to obtain a uniformly random bit string. In addition, the sequences xn{\displaystyle x_{n}} and xn+1{\displaystyle x_{n}+1} obtained by the logistic map are strongly correlated, which makes it problematic for pseudorandom sequences. One way to solve this is to generate the sequence x0,x1,x2,...{\displaystyle x_{0},x_{1},x_{2},...} for each iteration of the map, rather than generating the sequence x0,xτ,x2τ,...{\displaystyle x_{0},x_{\tau },x_{2\tau },...} for some number of iterations τ > 1. For example, it is said that good pseudorandom numbers can be obtained for method 1 with τ > 10 or τ > 13, and for method 2 with τ > 16.

A common problem with digitally calculating chaos using a computer is that, because a computer has a finite calculation precision, it is in principle impossible to obtain a truly aperiodic sequence, which is the nature of chaos, and instead outputs a finite periodic sequence. Even if aperiodic sequences cannot be obtained in principle, sequences with as long a period as possible are desirable for generating pseudorandom numbers. However, when the periodicity of the sequence actually output by the logistic map fr=4{\displaystyle f_{r=4}} in single-precision floating-point calculations was investigated, it was reported that the period of the sequence actually output is much smaller than the maximum period possible from the number of bits allocated, and from this point of view, it has been pointed out that pseudorandom number generation by the logistic map is inferior to existing pseudorandom number generators such as the Mersenne Twister. In addition, with the logistic map, fr=4{\displaystyle f_{r=4}} there is a risk that the value will fall to the fixed point 0 during the calculation and remain constant. On the other hand, the logistic map always takes values in the open interval (0, 1), so it can be calculated without problems not only with floating point but also with fixed point, and can enjoy the advantages of fixed point arithmetic. It has been pointed out that fixed point has a longer period than floating point for the same number of bits, and that unintended convergence to 0 can be eliminated.

Extension to complex numbers

Correspondence between the orbit diagram of a variation of the logistic map (top) and the Mandelbrot set (bottom)

Dynamical systems defined by complex analytic functions are also of interest. An example is the dynamical system defined by the quadratic function:

where the parameter c and the variable z are complex numbers. This map is essentially the same as the logistic map (1–2). As mentioned above, the map (6–3) is topologically conjugate to the logistic map (1–2) through a linear function.

When the iteration of the map (6–3) is calculated with a fixed parameter c and varying the initial value z0{\displaystyle z_{0}}, a set of z0{\displaystyle z_{0}} such that zn{\displaystyle z_{n}} does not diverge to infinity as n → ∞ is called a filled Julia set. Furthermore, the boundary of a filled Julia set is called a Julia set. When the iteration of the map (6–3) is calculated with a fixed initial value z0=0{\displaystyle z_{0}=0} and varying the parameter c, a set of c such that z does not diverge to infinity is called a Mandelbrot set. The Julia sets and Mandelbrot sets of the map (6–3) generate fractal figures that are described as "mystical looking" and "extremely mysterious".

In particular, in the Mandelbrot set, each disk in the diagram corresponds to a region of asymptotically stable periodic orbits of a certain period. By juxtaposing the logistic map orbit diagram with the Mandelbrot set diagram, it is possible to see that the asymptotically stable fixed points, period doubling bifurcations, and period-three windows of the logistic map orbit diagram correspond on the real axis to the Mandelbrot set diagram.

When there is a time delay

The trajectory of the delayed logistic map. The initial values (x0,y0){\displaystyle (x_{0},y_{0})} are the same in both figures, but at the bifurcation point r = 2, the trajectory is attracted to a closed curve (left) and a point (right).

If we interpret the logistic map as a model of the population of each generation of organisms, it is possible that the population of the next generation will affect not only the population of the current generation, but also the population of the generation before that. An example of such a case is

where the number of individuals in the previous generation, xn1{\displaystyle x_{n-1}}, is included in the equation as a negative density effect . If xn+1=yn{\displaystyle x_{n+1}=y_{n}}, then equation ( 6-4 ) can be replaced by the following two-variable difference equation .

This dynamical system is used to study bifurcation of quasi-periodic attractors and is called the delayed logistic map . The delayed logistic map exhibits a Neimark–Sacker bifurcation at r = 2, where the asymptotically stable fixed point becomes unstable and an asymptotically stable invariant curve forms around the unstable fixed point .

Research history

Before Chaos was named

Before the iteration of maps became relevant to dynamical systems, mathematicians Gaston Julia and Pierre Fatou studied the iteration of complex functions. Julia and Fatou's work was broad, focusing on analytic functions. In particular, they studied the behavior of the following complex quadratic function, also shown in equation (6–3), in the 1920s.

Julia and Fatu also recognized chaotic behavior in Julia sets, but because there was no computer graphics at the time, no one followed suit and their research stalled. Research on complex dynamical systems then declined until the late 1970s, and it was not until the appearance of Benoit Mandelbrot and others that the rich dynamical behavior exhibited by maps on the complex plane was noticed.

Stanislaw Ulam
John von Neumann

In 1947, mathematicians Stanislaw Ulam and John von Neumann wrote a short paper entitled "On combination of stochastic and deterministic processes" in which they

They pointed out that pseudorandom numbers can be generated by the repeated composition of quadratic functions such as. In modern terms, this equation corresponds to the logistic map with r = 4. At that time, the word "chaos" had not yet been used, but Ulam and von Neumann were already paying attention to the generation of complex sequences using nonlinear functions. In their report, Ulam and von Neumann also clarified that the map (7–2) and the tent map are topologically conjugate, and that the invariant measure of the sequence of this map is given by equation (3–17).

There have since been some detailed investigations of quadratic maps of the form with arbitrary parameter a. Between 1958 and 1963, Finnish mathematician Pekka Mylberg developed the

This line of research is essential for dynamical systems, and Mühlberg has also investigated the period-doubling branching cascades of this map, showing the existence of an accumulation point λ = 1.401155189.... Others, such as the work of the Soviet Oleksandr Sharkovsky in 1964, the French Igor Gumowski and Christian Mila in 1969, and Nicholas Metropolis in 1973, have revealed anomalous behavior of simple one-variable difference equations such as the logistic map.

Robert May's research

Robert May (photographed in 2009)

Later, in the early 1970s, mathematical biologist Robert May encountered the model of equation (1–2) while working on an ecological problem. May introduced equation (1–2), i.e., the logistic map, by discretizing the logistic equation in time. He mathematically analyzed the behavior of the logistic map, and published his results in 1973 and 1974. Numerical experiments were performed on the logistic map to investigate the change in its behavior depending on the parameter r.

In 1976, he published a paper in Nature entitled "Simple mathematical models with very complicated dynamics". This paper was a review paper that focused on the logistic map and emphasized and drew attention to the fact that even simple nonlinear functions can produce extremely complex behaviors such as period-doubling bifurcation cascades and chaos. This paper in particular caused a great stir and was accepted by the scientific community due to May's status as a mathematical biologist, the clarity of his research results, and above all, the shocking content that a simple parabolic equation can produce surprisingly complex behavior. Through May's research, the logistic map attracted many researchers to chaos research and became such a famous mathematical model that it is said to have restarted the flow of chaos research.

After May's research

May also drew attention to the paper by using the term "chaos", which was used by Tien-Yen Li and James York in their paper "Period three implies chaos". Although some disagree, Li and York's paper is considered the first to use the word "chaos" as a mathematical term, and is credited with coining the term "chaos" to refer to deterministic, chaotic behavior. Li and York completed the paper in 1973, but when they submitted it to The American Mathematical Monthly, they were told that it was too technical and that it should be significantly rewritten to make it easier to understand, and it was rejected. The paper was then left unrevised. However, the following year, in 1974, May came to give a special guest lecture at the University of Maryland where Lee and York were working, and talked about the logistic map. At the time, May did not yet understand what was happening in the chaotic domain of the logistic map, but Lee and York were also unaware of the period-doubling cascade of the logistic map. Excited by May's talk, Lee and York caught up with May after the lecture and told him about their results, which surprised May. Lee and York quickly rewrote the rejected paper, and the resubmitted paper was published in 1975.

Mitchell Feigenbaum (photographed in 2006)

Also, around 1975, mathematical physicist Mitchell Feigenbaum noticed a scaling law in which the branching values converged in a geometric progression when he looked at the period-doubling cascade of the logistic map, and discovered the existence of a constant, now called the Feigenbaum constant, through numerical experiments. May and George Oster had also noticed the scaling law, but they were unable to follow it in depth. Feigenbaum discovered that the same constant also appeared in the sine map shown in equation (4–1), and realized that this scaling law had a universality that went beyond the logistic map. In 1980, a rigorous proof of this result was given by Pierre Collé, Jean-Pierre Eckman, Oscar Rumford, and others. At about the same time as Feigenbaum, or later, physicists discovered the same period doubling cascade and the Feigenbaum constant in real life, and chaos, which had previously been seen as a strictly mathematical phenomenon, had a major impact on the field of physics as well.

However, there is criticism of the tendency to downplay research results from before the popularity of chaos, and to attribute many of those results to rediscoverers who used the logistic map, etc. May himself respects the existence of previous research, but positions his own achievement as not being "the first to independently discover the strange mathematical behavior of quadratic maps", but as being one of the "last researchers to emphasize their broad implications in science". Mathematician Robert Devaney states the following before explaining the logistic map in his book:

This means that by simply iterating the quadratic function fλ(x)=λx(1x){\displaystyle f_{\lambda }(x)=\lambda x(1-x)} (also known as the logistic map), we can predict the fate of the initial population x0{\displaystyle x_{0}}. This sounds simple, but I dare to point out that it was only in the late 1990s, after the efforts of hundreds of mathematicians, that the iteration of this simple quadratic function was fully understood.

Occurrences and similar systems

  • In a toy model for discrete laser dynamics: xGx(1tanh(x)){\displaystyle x\rightarrow Gx(1-\tanh(x))}, where x{\displaystyle x} stands for electric field amplitude, G{\displaystyle G}[27] is laser gain as bifurcation parameter.
  • Hofstadter sequences are an example of one dimensional quasi-random, aperiodic, chaotic sequences again defined by recursion, a very special case is the logistic map

See also

Notes

  1. Lorenz, Edward N. (1964-02-01). "The problem of deducing the climate from the governing equations". Tellus. 16 (1): 1–11. Bibcode:1964Tell...16....1L. doi:10.3402/tellusa.v16i1.8893. ISSN 0040-2826.
  2. Weisstein, Eric W."Logistic Equation". MathWorld.
  3. See § Research history.
  4. 123Gleick, James (1987). Chaos: Making a New Science. London: Penguin Books. ISBN 978-0-14-009250-9.
  5. Mandelbrot, Benoit B. (2004). Fractals and Chaos, The Mandelbrot Set and Beyond. pp. 259–267. doi:10.1007/978-1-4757-4017-2. ISBN 978-1-4419-1897-0.
  6. 1234Tsuchiya, Takashi; Yamagishi, Daisuke (February 11, 1997). "The Complete Bifurcation Diagram for the Logistic Map". Z. Naturforsch. 52a (6–7): 513–516. Bibcode:1997ZNatA..52..513T. doi:10.1515/zna-1997-6-708. S2CID 101491730.
  7. Bonet, J.; Martínez-Giménez, F.; Peris, A. (2001). "A Banach space which admits no chaotic operator". Bulletin of the London Mathematical Society. 33 (2): 196–198. doi:10.1112/blms/33.2.196. S2CID 121429354.
  8. Zhang, Cheng (October 2010). "Period three begins". Mathematics Magazine. 83 (4): 295–297. doi:10.4169/002557010x521859. S2CID 123124113.
  9. Bechhoefer, John (1996-04-01). "The Birth of Period 3, Revisited". Mathematics Magazine. 69 (2): 115–118. doi:10.1080/0025570X.1996.11996402. ISSN 0025-570X.
  10. Jeffries, Carson; Pérez, José (1982). "Observation of a Pomeau–Manneville intermittent route to chaos in a nonlinear oscillator". Physical Review A. 26 (4): 2117–2122. Bibcode:1982PhRvA..26.2117J. doi:10.1103/PhysRevA.26.2117. S2CID 119466337.
  11. Baumol, William J.; Benhabib, Jess (February 1989). "Chaos: Significance, Mechanism, and Economic Applications". Journal of Economic Perspectives. 3 (1): 77–105. doi:10.1257/jep.3.1.77.
  12. "Misiurewicz Point of the Logistic Map". sprott.physics.wisc.edu. Retrieved 2023-05-08.
  13. Teschl, Gerald (2012). Ordinary Differential Equations and Dynamical Systems. Amer. Math Soc. ISBN 978-0-8218-8328-0.
  14. Collet, Pierre; Eckmann, Jean-Pierre (1980). Iterated Maps on the Interval as Dynamical Systems. Birkhauser. ISBN 978-3-7643-3026-2.
  15. "Topological transitivity". Encyclopedia of Mathematics.
  16. 12Strogatz, Steven (2019). "10.1: Fixed Points and Cobwebs". Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering (2nd ed.). Boca Raton. ISBN 978-0-367-09206-1. OCLC 1112373147.{{cite book}}: CS1 maint: location missing publisher (link)
  17. Jakobson, M. (1981). "Absolutely continuous invariant measures for one-parameter families of one-dimensional maps". Communications in Mathematical Physics. 81 (1): 39–88. Bibcode:1981CMaPh..81...39J. doi:10.1007/BF01941800. S2CID 119956479.
  18. Campbell, Trevor; Broderick, Tamara (2017). "Automated scalable Bayesian inference via Hilbert coresets". arXiv:1710.05053 [stat.ML].
  19. 1234Schröder, Ernst (1870). "Ueber iterirte Functionen". Mathematische Annalen. 3 (2): 296–322. doi:10.1007/BF01443992. S2CID 116998358.
  20. Little, M.; Heesch, D. (2004). "Chaotic root-finding for a small class of polynomials"(PDF). Journal of Difference Equations and Applications. 10 (11): 949–953. arXiv:nlin/0407042. doi:10.1080/10236190412331285351. S2CID 122705492.
  21. 12Lorenz, Edward (1964). "The problem of deducing the climate from the governing equations". Tellus. 16 (February): 1–11. Bibcode:1964Tell...16....1L. doi:10.3402/tellusa.v16i1.8893.
  22. 12Delbourgo, R.; Hart, W.; Kenny, B. G. (1985-01-01). "Dependence of universal constants upon multiplication period in nonlinear maps". Physical Review A. 31 (1): 514–516. Bibcode:1985PhRvA..31..514D. doi:10.1103/PhysRevA.31.514. ISSN 0556-2791. PMID 9895509.
  23. Feigenbaum, M. J. (1976) "Universality in complex discrete dynamics", Los Alamos Theoretical Division Annual Report 1975-1976
  24. Feigenbaum, Mitchell (1978). "Quantitative universality for a class of nonlinear transformations". Journal of Statistical Physics. 19 (1): 25–52. Bibcode:1978JSP....19...25F. CiteSeerX 10.1.1.418.9339. doi:10.1007/BF01020332. S2CID 124498882.
  25. Okulov, A Yu; Oraevskiĭ, A N (1984). "Regular and stochastic self-modulation in a ring laser with nonlinear element". Soviet Journal of Quantum Electronics. 14 (2): 1235–1237. Bibcode:1984QuEle..14.1235O. doi:10.1070/QE1984v014n09ABEH006171.
  26. Shen, Bo-Wen; Pielke, Roger A.; Zeng, Xubin (2023-08-12). "The 50th Anniversary of the Metaphorical Butterfly Effect since Lorenz (1972): Multistability, Multiscale Predictability, and Sensitivity in Numerical Models". Atmosphere. 14 (8): 1279. Bibcode:2023Atmos..14.1279S. doi:10.3390/atmos14081279. ISSN 2073-4433.
  27. Okulov, A Yu; Oraevskiĭ, A N (1986). "Space–temporal behavior of a light pulse propagating in a nonlinear nondispersive medium". J. Opt. Soc. Am. B. 3 (5): 741–746. Bibcode:1986JOSAB...3..741O. doi:10.1364/JOSAB.3.000741. S2CID 124347430.

Citations

  1. Gulick 1995, p. 16, Example 3
  2. Gulick 1995, p. 36
  3. Gulick 1995, p. 9
  4. Gulick 1995, p. 10
  1. Devaney 1989, p. 27.
  2. Devaney 1989, p. 2.
  3. Devaney 1989, p. 50
  1. 1234May 1976.
  2. May 1976, formula 2 and 3

References