פתיחת התפריט הראשי

אלגוריתם לנצוש

אלגוריתם לָנְצוֹשׁ (Lanczos) הוא אלגוריתם שפיתח קורנליוס לנצוש ב-1950[1]. האלגוריתם מוצא חלק (או את כל) הערכים עצמיים והווקטורים העצמיים של מטריצה הרמיטית (או סימטרית) בגודל במספר מוגבל של פעולות. כמות הפעולות הנדרשת למציאת כל הערכים העצמיים והווקטורים העצמיים של מטריצה היא , אבל בעזרת אלגוריתם לנצוש ניתן למצוא חלק מהערכים העצמיים בפחות פעולות.[2]

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

האלגוריתםעריכה

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

פלט: מטריצה   בגודל   עם עמודות אורתונורמליות ומטריצה תלת-אלכסונית   בגודל  . הערכים העצמיים של המטריצה   "קרובים" לחלק מהערכים העצמיים של המטריצה המקורית  . בשביל לקבל את הערכים העצמיים של   אפשר להשתמש למשל באלגוריתם QR. מכיוון שהמטריצה   היא מטריצה תלת-אלכסונית, הערכים העצמיים שלה מחושבים ב   פעולות.

אלגוריתם:

  1. יהי   וקטור שרירותי עם נורמה אוקלידית  .
  2. צעד התחלתי:
    1. יהי  .
    2. יהי  .
    3. יהי  .
  3. עבור   בצע:
    1. יהי   (נורמה אוקלידית).
    2. אם   אז  ,
      אחרת בחר את   להיות וקטור שרירותי עם נורמה אוקלידית   שיהיה אורתוגונלי לכל  .
    3. יהי  .
    4. יהי  .
    5. יהי  .
  4. תהי   המטריצה עם עמודות  . תהי  .
הערה   עבור  .

מימושעריכה

אחד המימושים המקובלים לאלגוריתם נמצא בחבילה ARPACK[3]. לדוגמה, המימוש באוקטבה (Octave) מתבסס על מימוש זה[4].

לקריאה נוספתעריכה

  • Trefethen, Lloyd N., and David Bau III. Numerical linear algebra. Vol. 50. Siam, 1997.
  • Stewart, Gilbert W. Matrix Algorithms: Volume II: Eigensystems. Society for Industrial and Applied Mathematics, 2001.

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

  1. ^ Lanczos, C. "An iteration method for the solution of the eigenvalue problem of linear differential and integral operators", J. Res. Nat’l Bur. Std. 45, 255-282 (1950).
  2. ^ קיים אנלוג לאלגוריתם עבור מטריצות שאינן סימטריות בשם אלגוריתם ארנולדי
  3. ^ http://www.caam.rice.edu/software/ARPACK/
  4. ^ https://www.gnu.org/software/octave/doc/interpreter/Sparse-Linear-Algebra.html#doc_002deigs