כיצד להשתמש Scipy עבור שורשים נומריים למצוא בהנדסה סימלציה

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

שם מקור: SciPy Root Finding

SciPy מציע מספר פונקציות למציאת שורש, בעיקר בתוך מודול:0. . פונקציות אלה יכולות להתמודד עם משוואות בודדות או מערכות של משוואות.הם מתאימים לסוגים שונים של בעיות, כולל משוואות לא לינאריות נפוצות בסימולציות הנדסיות.

שימוש ב-'root' Function

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

דוגמה:

⁇ python

ייבוא מספרד כמו np

מקור: scipy.optimizeיבו שורש

משוואה (x):

חזרה x*2 - 4

ראשי התיבות של: [1]

פתרון = שורש (הפצה, ראשית gues)

הדפסה (רזולוציה.x)

⁇

מערכות של משוואות

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

דוגמה:

⁇ python

מערכת הגנה (x):

[x]2 + x[1] - 1 x [x]

ראשי התיבות של: (0.5, 0.5)

פתרון = שורש (מערכת, ראשית gues)

הדפסה (רזולוציה.x)

בחירת השיטה הנכונה

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

ניתן להגדיר שיטות באמצעות פרמטר פרמטר (FLT 3: 3):

⁇ python

פתרון = שורש (הפצה, ראשונית gues, Method='hybr')

⁇