קבוצת מנדלברוט

**קבוצת מַנְדֶלְבְּרוֹט** היא קבוצה של מספרים מרוכבים שמוגדרת באמצעות כללים פשוטים, אך הקבוצה כולה היא בעלת מורכבות גדולה. תנאי השייכות לקבוצה עדין ביותר, ובקרבת שפתה מופיעה התנהגות פרקטלית בה התמונה נשארת מורכבת בכל קנה מידה. ויזואליזציה של הקבוצה הופיע לראשונה במאמר של ברוקס ומטלסקי שנכתב ב 1978. זמן קצר לאחר מכן ובאופן בלתי תלוי, החל בנואה מנדלברוט לחקור אותה. עבודותיו של מנדלברוט זכו לפרסום רב והקבוצה נקראה על שמו. קבוצת מנדלברוט מהווה אובייקט מחקר בסיסי בתחום הפרקטלים, ובעזרתה חוקרים תופעות במערכות דינמיות מרוכבות.

## היסטוריה
קבוצת מנדלברוט היא אחד הנושאים הנחקרים בתחום הפרקטלים, שהחל לתפוס תאוצה בתחילת המאה ה-20, כשבין אבותיו נמנים פייר פאטו וגסטון ז'וליה. איור  של קבוצת מנדלברוט, ברזולוציה נמוכה מאוד, הופיע לראשונה במאמר של מתמטיקאים רוברט ברוקס ופטר מטלסקי  בשנת 1978. בשנת 1980 גילה בנואה מנדלברוט את הקבוצה תוך ביצוע אנימציות של קבוצת ז'וליה במחשב, וכתב מאמר ראשוני בו היא נכללה. בשנת 1985 פרסמו שני מתמטיקאים נוספים מאמר נרחב על קבוצה זו, וקראו לה על שם מנדלברוט. בהמשך, הפכה קבוצת מנדלברוט לאובייקט מחקר בסיסי בתורת המערכות הדינמיות המרוכבות.

החל משנת 1985 הפך עניין שרטוט קבוצת מנדלברוט למושא עניינם של רבים, וגם פורסמו אלגוריתמים לשרטוט הקבוצה ותמונות צבעוניות שלה.

## הגדרה
קבוצת מנדלברוט מוגדרת באופן הבא:
לכל מספר מרוכב c, אפשר להגדיר באופן רקורסיבי סדרה, שאיברה הראשון a_0(c)=0, והמשכה מחושב על-פי הכלל a_{n+1}(c)=a_n(c)^2+c. הסדרה עשויה להיות חסומה או בלתי-חסומה, תלוי בערכו של c. קבוצת מנדלברוט מורכבת מן המספרים c \in \mathbb{C} שעבורם הסדרה a_n(c) חסומה.

מנקודת מבט שונה במקצת, לכל c, מוגדרת קבוצת ז'וליה J_c כשפה של קבוצת הנקודות z \in \mathbb{C} שעבורן הסדרה f_c^n(z) חסומה (כאשר f_c(z) = z^2+c). ז'וליה ופאטו הוכיחו שאם הסדרה f_c^n(0) חסומה (כלומר c שייכת לקבוצת מנדלברוט) אז J_c קשירה, ואם היא אינה חסומה, אז J_c בלתי קשירה לחלוטין. משפט זה מוביל להגדרה שקולה נוספת של הקבוצה: זהו המקום הגאומטרי הקשיר של משפחת הפולינומים f_c(z)=z^2+c, כלומר אוסף הנקודות עבורן קבוצת ז'וליה המתאימה היא קשירה.

## המערכת הדינמית <math>z\mapsto f_c(z)=z^2+c</math>

השאלה האם מספר c נמצא בקבוצת מנדלברוט ומקומו בקבוצה זאת מספקת מידע רב על המערכת הדינמית המוגדרת על ידי הפונקציה  f_c(z):=z^2+c. זאת אחת המערכות הדינמיות הפשוטות ביותר המוגדרות על ידי פונקציה לא הפיכה: מדובר בפולינום ממעלה שנייה במשתנה מְרוּכַּב אחד. מערכות דינמיות המוגדרות על ידי פונקציות ליניאריות הן פשוטת מאוד. כל מערכת דינמית המוגדרת על ידי פולינום ממעלה שנייה ממשתנה מרוכב אחד שקולה למערכת המוגדרת על ידי הפונקציה f_c. כך שניתן לחשוב על דוגמה זאת כעל "הדוגמה הלא טריוויאלית הראשונה". קבוצת מנדלברוט מאפשרת טיפול סימולטני בכל המערכות הדינמיות האלה לכל c.

למרות פשטותה היחסית של מערכת זו היא מדגימה התנהגויות מורכבות מאוד כפי שניתן לראות מהמורכבות של קבוצת מנדלברוט. ישנן שאלות פתוחות רבות לגבי התנהגויות אלה אך גם תוצאות רבות שכבר הוכחו.

### שאלות כלליות בחקר מערכות דינמיות
באופן כללי מערכת דינמית בדידה היא זוג, המורכב מקבוצה X ופונקציה f מ-X לעצמה. במקרה שלנו X היא המישור המרוכב \C ו-f היא הפונקציה f_c עבור פרמטרc\in \C.

השאלה הבסיסית בחקר מערכות דינמיות היא זו: בהינתן נקודה x_0\in X, איך מתנהגת הסדרה f^{n}(x_0) שמתקבלת מהפעלת הפונקציה על x_0 ואז על תוצאת ההפעלה וכך הלאה עד אינסוף. סדרות כאלה נקראות מסלולים.

#### יציבות של נקודות שבת ומסלולים מחזוריים
התנהגות אפשרית אחת היא התכנסות לנקודת שֶבת. נקודת שֶבֶת היא נקודה x\in X המקיימת f(x)=x. נקודת שבת x\in X נקראת יציבה אם בהינתן נקודה x_0\in X הקרובה מספיק ל x, אז הסדרה f^n(x_0) מתכנסת לx.

בהינתן נקודת שבת יציבה יש מסלולים רבים שמתכנסים אליה. משימת איתור נקודות שבת שקולה לפתרון משוואה, ולכן קלה יחסית. בדרך כלל ניתן לבדוק האם נקודת שבת היא יציבה על ידי ניתוח הנגזרת של הפונקציה בנקודת השבת, ולכן גם זו משימה לא קשה.

מיון נקודות השבת של מערכת דינמית היא השלב הראשון בחקירתה.

התנהגות אפשרית אחת היא היתכנות למסלול מחזורי. מסלול מחזורי הוא מסלול f^n(x) המקיים f^{n+k}(x)=f^{n}(x) עבור k מסוים וכל n. מסלול שבו במקום מסוים מופיעה הנקודה המקורית. אומרים שמסלול מתכנס למסלול מחזורי אם החל ממקום מסוים הוא קרוב כרצוננו למסלול מחזורי. מסלול מחזורי נקרא יציב אם בהינתן נקודה x_0 קרובה מספיק לאחת מנקודות המסלול, הסדרה f^n(x_0) מתכנסת למסלול המחזורי.

ניתוח המסלולים המחזוריים ממחזור נתון k ויציבותם שקול לניתוח נקודות השבת של f^k. מכאן שעבור k ספציפי מדובר במשימה פשוטה למדי. אולם ניתוח כל המסלולים המחזוריים ויצובותם הוא משימה מסובכת בהרבה.

#### קבוצות ז'וליה ופאטו

קל יותר לנתח מערכות דינמיות כאשר הקבוצה X היא קומפקטית. המישור המרוכב אינו קומפקטי. ניתן להחליף אותו בספירה של רימן (כלומר להוסיף לו נקודה באינסוף) ואז נקבל מרחב קומפקטי. קל לראות שהפונקציה f_c מתרחבת ברציפות לספירה של רימן (על ידי כך שקובעים את האינסוף להיות נקודת שבת שלה). מכאן ניתוח המערכת הדינמית שלנו שקול לניתוח המערכת הדינמית שמתקבלת על הספירה של רימן.

בהינתן מערכת דינמית על קבוצה קומפקטית X ניתן לחלק את X לשתי קבוצות:
* קבוצת פאטו - אוסף הנקודות x_0 כך ששינוי קטן בהן מוביל לשינוי קטן בסדרה f^n(x_0)
* קבוצת ז'וליה - אוסף הנקודות x_0 כך ששינוי קטן בהן מוביל לשינוי גדול בסדרה f^n(x_0)

### תכונות ייחודיות למערכת הדינמית <math>z\mapsto f_c(z)=z^2+c</math>
עבור המערכת הדינמית המוגדרת על ידי הפונקציה f_c(z) מתקימות התכונות הבאות:
# עבור כל x_0 המקיים |x_0|>2 הסדרה f_c^n(x_0) שואפת לאינסוף. זה מאפשר לקבוע בקלות יחסית שנקודה מסוימת לא נמצאת בקבוצת מנדלברוט.
# למערכת יש לא יותר ממסלול מחזורי אחד.
# למערכת יש מסלול מחזורי אם ורק אם המסלול המתחיל ב-0 מתכנס למסלול מחזורי.
# קבוצת ז'וליה של המערכת קשירה אם ורק אם c נמצאת בקבוצת מנדלברוט.
# אם c לא נמצאת בקבוצת מנדלברוט אז קבוצת ז'וליה של המערכת הומיאומורפית לקבוצת קנטור ובפרט בלתי קשירה לחלוטין.
# אם למערכת יש נקודת שבת יציבה אז קבוצת ז'וליה של המערכת היא עקום פשוט (רציף ובדרך כלל לא חלק) סגור. בהתאם לפי משפט העקום של ז'ורדן, לקבוצת פאטו של המערכת יש בדיוק שני רכיבי קשירות (האחד תחום על ידי קבוצת ז'וליה והאחר מחוצה לה).
# אם למערכת יש מסלול מחזורי יציב (שאינו נקודת שבת) אז לקבוצת פאטו של המערכת יש אינסוף רכיבי קשירות.

## המבנה של קבוצת מנדלברוט
במבט חטוף על קבוצת מנדלברוט ניתן לראות שהיא מורכבת מ"חלקים" המהווים צורת גאומטריות פשוטת יחסית (לא פרקטלים). ניתן להגדיר באופן ריגורוזי חלקים אלה בתור רכיבי קשירות של הפנים של קבוצת מנדלברוט. משערים שפנים זה צפוף בקבוצת מנדלברוט.

### הקרדיואידה הראשית
הגדול והבולט בין רכיבים אלו נקרא הקרדיואידה הראשית.
צורתה הגאומטרית היא קרדיואידה. ניתן להגדיר אותה כאוסף כל הנקודות c עבורם למערכת z\mapsto f_c(z):=z^2+c יש נקודת שבת יציבה. קל לראות שקבוצה זאת מתוארת על ידי התנאי \exist z \text{ s.t. } f(z)=z \text{ and }  |f'(z)|
קל לראות שתנאי זה מתקיים אם ורק אם קיים \mu בעיגול היחידה כך ש:
c= \frac\mu2\left(1-\frac\mu2\right)
מכאן שהשפה של הקבוצה הזאת מתוארת על ידי המשוואה הפרמטרית:
c(t)= \frac{e^{it}}2\left(1-\frac{e^{it}}2\right)

### רכיבים היפרבוליים
רכיב נקרא רכיב היפרבולי אם עבור כל נקודה בו למערכת יש מסלול מחזורי יציב. מחזור המסלול במקרה זה הוא קבוע ונקרא המחזור של הרכיב. משערים שכל רכיבי הקשירות של הפנים של קבוצת מנדלברוט הם היפרבוליים.

עבור n נתון, תיאור כל הרכיבים ההיפרבוליים ממחזור n היא משימה הניתנת לביצוע באמצעות חישוב: התנאי לכך שנקודה נמצאת ברכיב היפרבולי ממחזור n הוא \exist z \text{ s.t. } f^{n}(z)=z \text{ and }  |(f^{n})'(z)| זהו תנאי אלגברי-למחצה ולכן על פי משפט זינדנברג-טרסקי (שיש לו הוכחות קנסטרוקטיביות) ניתן לשכתב תנאי זה ללא כמתים אלא רק עם אי-שוויונים פולינומיים על החלק הממשי והמרוכב של c. לאחר מכן ניתן לחשב את מספר רכיבי הקשירות ואת האי-שוויון שמתאר כל אחד מהם. אולם כאשר n גדל חישוב ישיר נהיה לא מעשי מהר מאוד.

עבור נקודה c הנמצאת ברכיב היפרבולי ממחזור n לקבוצת פאטו של המערכת הנתונה יש אינסוף רכיבי קשירות. n הנקודות של המסלול המחזורי נמצאים ב-n רכיבים שונים של קבוצת פאטו. הסגורים של רכיבים אלה נוגעים בנקודה אחת. ניתן להגדיר סדר ציקלי על רכיבים אלה לפי הסדר שבו הם נוגעים זה בזה (נגד כיוון השעון). סדר זה נותן סדר ציקלי על נקודות המסלול המחזורי. ניתן להראות שהפונקציה f_c מבצעת הזזה בקבוע של נקודות המסלול המחזורי לפי הסדר הציקלי הנ"ל. גודל הזזה k זר ל-
n. מכאן שאנו מקבלים שבר מצומצם \frac{k}{n}. מספר אינו תלוי בנקודה c ברכיב ההיפרבולי ובכך מהווה אינבריאנט של הרכיב ההיפרבולי.

### רכיבים היפרבוליים ראשיים
ניתן למיין ולתאר במפורש את כל הרכיבים המשיקים לקרדאוידה הראשית. רכיבים אלה נקראים רכיבים היפרבוליים ראשיים. הם כולם בצורת עיגול. לכל מספר רציונלי r בקטע [0,1] קיים ויחיד רכיב היפרבולי ראשי שהאינבריאנט הרציונלי שהוגדר מעלה עבורו הוא r. רכיב זה משיק לקרדיאוד הראשי בנקודה c(t)= \frac{e^{ir}}2\left(1-\frac{e^{ir}}2\right)

### אנטנות
בסמוך לכל רכיב היפרבולי ראשי, ניתן לזהות חלק של קבוצת מנדלברוט הנראה כמעין קו (מאוד לא חלק ומכיל לאורכו רכיבים נוספים) שמתפצל למספר קווים הנראים גם הם כך. חלקים אלה של קבוצת מנדלברוט נקראים לעיתים "אנטנות". ניתן להראות שהמספר הכולל של האנטנות בהתפצלות הראשונה בסמוך לרכיב הפרבולי ראשי שווה למחזור של רכיב זה.

## הצגת הקבוצה כפרקטל
את קבוצת מנדלברוט מתארים במישור המרוכב שבו הצירים מייצגים את החלק הממשי והמדומה של כל מספר. הנקודות השייכות לקבוצה (כלומר, הערכים של c שעבורם הסדרה a_n(c) חסומה) נצבעות בשחור, וכל נקודה אחרת מקבלת צבע התלוי במספר האיברים בסדרה שערכם המוחלט קטן ממספר מסוים – 2 בדרך כלל.

## הדגמה
בהדגמה שלהלן מתבצע שינוי קנה המידה. כל תמונה היא הגדלה של אזור מסוים בתמונה הקודמת. ההגדלה הכוללת מהתמונה הראשונה לאחרונה מגיעה עד כדי 60 מיליארד.

## אמנות וקבוצת מנדלברוט
יש אנשים שנהנים למצוא צורות גאומטריות מעניינות בקבוצת מנדלברוט. להלן דוגמאות לאזורים יפים בקבוצת מנדלברוט:

## הכללה
ממוזער|שמאל|קבוצת מנדלברוט המוכללת עבור <math>0 \le d \le 5</math>
את קבוצת מנדלברוט ניתן להכליל באופן הבא: במקום להביט בפונקציה f_c(z) = z^2+c ולהציב בה ערכים, מביטים בפונקציה הכללית f_{c.d}(z) = z^d+c, עבור d ממשי כלשהו. כאשר d מספר שלם, המשפט אודות הקשר בין הגדרה זו להגדרה בעזרת קבוצת ז'וליה נשאר תקף.

כאשר d משתנה, ניתן לשרטט את הקבוצה ולקבל תמונה שונה בכל פעם. כך למשל, כאשר d נמצא בין 0 ל-5 (כמו בתמונה המצורפת), ניתן לראות את קבוצת מנדלברוט "נולדת" מכלום ומתפתחת.

## קירוב לפאי בעזרת קבוצת מנדלברוט
שימוש נוסף של קבוצת מנדלברוט הוא האפשרות לקרב את פאי בעזרת מונחים וכלים מנושא זה. כידוע, המספר c=\frac{1}{4} הוא המספר הממשי הגדול ביותר בקבוצת מנדלברוט; עוד ידוע כי כל פעם שהסדרה a_n(c) חסומה, אז היא קטנה מ-2 בהכרח. כעת, לכל מספר ממשי c> \frac{1}{4}, מסמנים ב-N_c את מספר הפעמים שיש לבצע את האיטרציה עד שעוברים את 2, כלומר המספר הראשון כך ש-a_{N_c}(c) > 2.

מסתבר שהסדרה n \cdot N_{\frac{1}{4}+1/n^2} מתכנסת ל-\pi.

נקודות נוספות בקבוצת מנדלברוט אשר מאפשרות לקרב את פאי כוללות נקודות הנמצאות על התפר בין הצורות של הקבוצה. עבור נקודות כאלו, המסלול אשר בעזרתו שואפים אל הנקודה משתנה ולרוב איננו כולל מספרים ממשיים.

### הסבר אינטואיטיבי אודות הקשר בין קבוצת מנדלברוט לפאי
עבור הפונקציה הממשית f(x) = x^2 + \frac{1}{4}, הישר y=x משיק לפונקציה; לכן, בעת הפעלת התהליך של סדרת מנדלברוט, נוצרת מעין פונקציית מדרגה תחומה בין הפרבולה לבין הישר, אשר שואפת בגובהה אל \frac{1}{2} אך לא מגיעה. לעומת זאת, בפונקציה f(x) = x^2 + \frac{1}{4} + \epsilon, הישר y=x לא משיק לפרבולה, ולכן נוצרת ביניהם פונקציית מדרגה. במובן מסוים, פאי הוא ה"זמן" שלוקח לאותה פונקציית מדרגה לעבור את \frac{1}{2}.

הסבר נוסף כלהלן: כאמור לעיל, הפונקציה הממשית f(x) = x^2 + \frac{1}{4} + \epsilon עוברת את \frac{1}{2}. צורת הגרף סביב ערך זה (בטווח) מתקרבת לצורה של פונקציית הטנגנס ככל שאפסילון קטן יותר. הקשר של פונקציה זו לפאי מעולם הטריגונומטריה ברור.

## ראו גם
*דינמיקה הולומורפית
** קבוצת מולטיברוט
** קבוצת ז'וליה
* פרקטל

## קישורים חיצוניים

* 
* "Mandelbrot Set - Online Generator"
* "עשר דקות של תקריב לתוך הפרקטל"
* קורס על דינמקה מרוכבת הכלל את ההגדרה של קבוצות מנדלברוט.
* וידאו של 3Blue1Brown על דינמיקה הולומורפית.
* וידאו של Matheologer על דינמיקה הולומורפית.
* וידאו של numberphile על דינמיקה הולומורפית.
* סדרת סרטוני יוטיוב של numberphile מאת הולי קריגר. חלק גדול מהסרטונים עוסקים בקבוצת מנדלברוט.
* Robert L. Devaney, The Complex Geometry of the Mandelbrot Set, Boston University

## הערות שוליים

קטגוריה:ערכים שבהם תבנית בריטניקה אינה מתאימה
קטגוריה:פרקטלים
קטגוריה:מערכות דינמיות