עקרונות הסימולציה הממוחשבת

פרק 1

מהי סימולציה ממוחשבת של מעגל?

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

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

שלוש אבני הבניין של כל סימולציה

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

מטרת הפרק

נראה כיצד שלושת הרכיבים הללו משתלבים יחד. נצא מבחירת צומת הייחוס, נגיע לבניית מטריצת האדמיטנס $$Y$$ ולווקטור הזרמים $$I$$ , ונסיים בפתרון $V=Y^{-1}I$ ובחילוץ זרם בודד באמצעות כלל קרמר. בסוף נראה כיצד לממש את כל התהליך בכמה שורות של MATLAB.

פרק 2

בחירת צומת ייחוס

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

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

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

דוגמה: מעגל עם ארבעה צמתים

נניח מעגל בעל ארבעה צמתים $$a,b,c,d$$ שאליו מחוברים מספר מקורות זרם, כל אחד בעוצמה $1\,\mathrm{A}$ . אם נבחר את $$b$$ כצומת ייחוס, נגדיר $$V_b=0$$ ונאמוד שלושה מתחי צמתים בלבד – $$V_a, V_c, V_d$$ . הזרם $$i_2$$ בענף כלשהו ייגזר ישירות מהפרשי המתחים שנקבל. בחירת צומת ייחוס שונה תשנה את הערכים המוחלטים של מתחי הצמתים, אך לא תשנה את הזרם בענף – זה אינווריאנט פיזיקלי של המעגל.

שימו לב

בחירת צומת ייחוס שגוי – למשל בחירת צומת מבודד או צומת המחובר בלבד דרך קבל אידיאלי ב־DC – יכולה להפוך את מטריצת האדמיטנס לסינגולרית (לא הפיכה). הפותר יתריע על Floating Node, וכל הסימולציה תיכשל. וודאו תמיד שצומת הייחוס מחובר חשמלית לשאר המעגל.

פרק 3

ניתוח צמתים – הבסיס המתמטי

לאחר בחירת צומת הייחוס, מתבצע צעד אלגברי מרכזי: ניתוח צמתים (Nodal Analysis). על כל אחד מהצמתים שאינם צומת הייחוס מופעל חוק הזרמים של קירכהוף (KCL) – סכום הזרמים היוצאים מהצומת שווה לסכום הזרמים הנכנסים אליו.

\sum_{k}i_{k}(t)=0\qquad\text{at every non-reference node}

זרם דרך רכיב פסיבי מבוטא באמצעות הפרש המתחים שעל פניו והאדמיטנס שלו. עבור נגד בעל התנגדות $$R$$ והאדמיטנס המתאים $$G=1/R$$ , הזרם בין צומת $$i$$ לצומת $$j$$ הוא:

i_{ij}=G\bigl(V_{i}-V_{j}\bigr)

הצבת ביטויים אלה במשוואות ה־KCL לכל צומת מובילה למערכת ליניארית של משוואות במשתני מתחי הצמתים. הצורה המסודרת של המערכת – לאחר איסוף איברים – היא בדיוק המבנה המטריציוני שמהווה את לב הסימולציה.

למה דווקא ניתוח צמתים?

קיימת שיטה אלטרנטיבית – ניתוח לולאות (Mesh Analysis) – העובדת על זרמי לולאה במקום מתחי צמתים. עם זאת, סימולטורים מודרניים מעדיפים כמעט תמיד את ניתוח הצמתים בגרסה המורחבת (MNA – Modified Nodal Analysis). הסיבה היא שמספר משתני המתח קטן בדרך כלל ממספר משתני הלולאה, וקל יותר לתת לכל מקור או רכיב מודל אדמיטנס.

פרק 4

בניית מטריצת האדמיטנס

לב הסימולציה הוא מטריצה ריבועית $$Y$$ בגודל $n\times n$ , כאשר $$n$$ הוא מספר הצמתים שאינם צומת הייחוס. כל איבר במטריצה מקיים כלל הרכבה פשוט ושיטתי:

Y_{ii}=\sum_{k}G_{ik},\qquad Y_{ij}=-G_{ij}\quad(i\neq j)

המשמעות בעברית: על האלכסון הראשי של המטריצה ניצב סכום האדמיטנסים המחוברים לצומת $$i$$ . מחוץ לאלכסון, באיבר $Y_{ij}$ , ניצב מינוס האדמיטנס המקשר את הצומת $$i$$ לצומת $$j$$ . אם בין שני הצמתים אין חיבור ישיר, האיבר המתאים שווה לאפס.

תכונות מבניות של מטריצת האדמיטנס

סימטריות

במעגלים פסיביים, האדמיטנס בין צמתים $$i$$ ו־ $$j$$ אינו תלוי בכיוון, ולכן $Y_{ij}=Y_{ji}$ . המטריצה סימטרית – תכונה המנוצלת על ידי הפותר לשיפור ביצועים.

דומיננטיות אלכסונית

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

דלילות (Sparsity)

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

דוגמה אינטואיטיבית – ארבעה צמתים

נניח שלאחר בחירת צומת הייחוס נותרו שלושה צמתים פעילים, ובינם מחוברים נגדים עם אדמיטנסים $G_{1},G_{2},G_{3},G_{4}$ . מטריצת האדמיטנס נראית מבחינה סכמטית כך:

Y=\begin{bmatrix} G_{1}+G_{2} & -G_{2} & 0\\ -G_{2} & G_{2}+G_{3}+G_{4} & -G_{4}\\ 0 & -G_{4} & G_{4}+G_{5} \end{bmatrix}

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

פרק 5

וקטור הזרמים

במקביל לבניית מטריצת האדמיטנס, מורכב וקטור הזרמים העצמאיים, המסומן $$I$$ . כל איבר $I_{i}$ מייצג את הזרם הנטו שמוזרק לתוך צומת $$i$$ על ידי מקורות זרם חיצוניים.

I_{i}=\sum I_{\text{in},i}-\sum I_{\text{out},i}

המוסכמה הנפוצה: זרם שנכנס לצומת נחשב חיובי, וזרם היוצא ממנו נחשב שלילי. אם לצומת אין מקור זרם חיצוני, האיבר המתאים בוקטור $$I$$ שווה לאפס.

טיפול במקורות מתח

סימולטורים אמיתיים אינם מטפלים במקורות מתח כפי שהם, אלא ממירים אותם למקורות זרם שקולים באמצעות טרנספורמציית נורטון, או משתמשים בגרסה מורחבת בשם MNA (Modified Nodal Analysis). בגרסה זו, זרם המקור עצמו נכנס כמשתנה לא ידוע נוסף למערכת – ולכן המטריצה גדלה בשורה ובעמודה לכל מקור מתח.

תובנה חשובה

הוקטור $$I$$ הוא הצד הימני של מערכת המשוואות, וקובע את "התרגיל" שעל המחשב לפתור. מטריצת $$Y$$ מתארת את הטופולוגיה והערכים של הרכיבים הפסיביים – היא תכונה של המעגל. הוקטור $$I$$ לעומת זאת מתאר את ההפעלה החיצונית של המעגל – הוא תכונה של הניסוי.

פרק 6

פתרון מערכת המשוואות

לאחר שבנינו את המטריצה $$Y$$ ואת הוקטור $$I$$ , אנו מקבלים את המשוואה הליניארית המרכזית של ניתוח המעגל:

Y\,V=I

כאשר $$V$$ הוא וקטור מתחי הצמתים הלא ידועים. הפתרון הפורמלי הוא:

V=Y^{-1}\,I

בכל סימולטור, הפותר אינו מחשב את $Y^{-1}$ במפורש. במקום זאת הוא משתמש בפירוקי מטריצה יעילים – LU עבור מטריצות צפופות, ופירוקים מיוחדים למטריצות דלילות. הפתרון בפועל מבוצע באופן הבא:

פירוק LU

המטריצה $$Y$$ מתפרקת למכפלת $$L$$ משולשת תחתונה ו־ $$U$$ משולשת עליונה: $$Y=LU$$ . הפירוק הוא נומרית יציב כאשר המטריצה דומיננטית אלכסונית.

פתרון קדמי

פתרון $$LZ=I$$ עבור וקטור הביניים $$Z$$ . מערכת זו פתירה בקלות בגלל המבנה המשולש של $$L$$ , שורה אחר שורה.

פתרון אחורי

פתרון $$UV=Z$$ עבור וקטור מתחי הצמתים. שוב, המבנה המשולש מאפשר פתרון מהיר, מהשורה האחרונה כלפי מעלה.

חישוב זרמי הענפים

לאחר ש־ $$V$$ ידוע, זרמי הענפים נגזרים ישירות: $i_{ij}=G_{ij}(V_{i}-V_{j})$ . צעד זה הוא לולאה פשוטה על הענפים, ללא חישוב מטריציוני נוסף.

פרק 7

כלל קרמר לחילוץ מתח או זרם בודד

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

V_{k}=\frac{\det\bigl(Y_{k}\bigr)}{\det(Y)}

המטריצה $Y_{k}$ נוצרת מן המטריצה $$Y$$ על ידי החלפת העמודה ה־ $$k$$ שלה בוקטור הזרמים $$I$$ . הצורך לחשב שתי דטרמיננטות בלבד מאפשר עיתים פתרון אנליטי קומפקטי, במיוחד למעגלים קטנים שבהם מתבצע ניתוח ידני.

מתי כלל קרמר עדיף?

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

זרם בענף באמצעות מתחי הצמתים

לאחר חישוב $V_{i}$ ו־ $V_{j}$ באמצעות כלל קרמר, הזרם בענף ביניהם – למשל הזרם $i_{2}$ בענף בעל אדמיטנס $G_{2}$ – ניתן ישירות על ידי:

i_{2}=G_{2}\bigl(V_{i}-V_{j}\bigr)

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

פרק 8

דוגמת חישוב מלאה

כדי לחבר את כל הפרקים הקודמים, נשתמש בדוגמה מייצגת המבוססת על מבנה טיפוסי של מעגל DC: מעגל בעל ארבעה צמתים $$a,b,c,d$$ , אליו מחוברים מספר מקורות זרם בעוצמה $1\,\mathrm{A}$ , ובינם נגדים בעלי אדמיטנסים שונים. המטרה: לחשב את הזרם $i_{2}$ בענף מסוים.

Example circuit for nodal analysis and admittance matrix formulation

איור 1: מעגל לדוגמה עבור ניתוח צמתים. הצומת $b$ נבחר כצומת הייחוס, ולכן $V_b=0$. לאחר בחירת הצומת, הסימולטור בונה את מטריצת האדמיטנס $Y$, את וקטור הזרמים $I$, פותר את המערכת $YV=I$, ולבסוף מחשב את הזרם $i_2$ מתוך הפרש המתחים בין הצמתים $a$ ו־$c$.

שלב 1 – בחירת צומת ייחוס

נבחר את הצומת $$b$$ כצומת ייחוס, ולכן $V_{b}=0$ . נותרו שלושה מתחי צמתים לא ידועים: $V_{a}, V_{c}, V_{d}$ .

שלב 2 – בניית מטריצת האדמיטנס

לפי כללי ההרכבה מפרק 4, נבנה מטריצה סימטרית $3\times 3$ . נסמן באופן כללי את האדמיטנסים הסכמטיים הקרובים לצמתים:

Y=\begin{bmatrix} G_{a} & -G_{ac} & -G_{ad}\\ -G_{ac} & G_{c} & -G_{cd}\\ -G_{ad} & -G_{cd} & G_{d} \end{bmatrix}

איברי האלכסון $G_{a},G_{c},G_{d}$ הם סכום האדמיטנסים המחוברים לכל צומת בהתאמה. שאר האיברים הם מינוס האדמיטנסים המקשרים בין שני הצמתים המתאימים.

שלב 3 – וקטור הזרמים

לפי כיוון מקורות הזרם בדוגמה, נניח שלצומת $$a$$ נכנס זרם נטו של $+1\,\mathrm{A}$ , לצומת $$c$$ זרם של $-1\,\mathrm{A}$ , ולצומת $$d$$ זרם של $+1\,\mathrm{A}$ :

I=\begin{bmatrix} +1\\ -1\\ +1\end{bmatrix}\quad[\mathrm{A}]

שלב 4 – פתרון מתחי הצמתים

פותרים את המערכת $Y\,V=I$ ומקבלים את וקטור המתחים:

V=Y^{-1}\,I=\begin{bmatrix} V_{a}\\ V_{c}\\ V_{d}\end{bmatrix}

שלב 5 – חילוץ הזרם $i_{2}$

נניח שהזרם $i_{2}$ זורם בענף שבין הצמתים $$a$$ ו־ $$c$$ , וכי האדמיטנס של הענף הוא $G_{ac}$ . אז:

i_{2}=G_{ac}\bigl(V_{a}-V_{c}\bigr)

לחילופין, באמצעות כלל קרמר ניתן לחשב את $V_{a}$ ואת $V_{c}$ ישירות מבלי לפתור את כל המערכת – כל אחד מהם הוא יחס של דטרמיננטות $3\times 3$ , וההפרש ביניהם מוכפל ב־ $G_{ac}$ נותן את התשובה הסופית.

אינווריאנט פיזיקלי

אם נחזור על אותה דוגמה כאשר צומת הייחוס הוא $$c$$ במקום $$b$$ , נקבל ערכים שונים עבור $V_{a},V_{b},V_{d}$ , אך ההפרש $V_{a}-V_{c}$ (וכמובן הזרם $i_{2}$ ) יישאר זהה לחלוטין. זוהי תכונה יסודית: הזרם הוא גודל פיזיקלי, ובחירת צומת הייחוס היא רק קונבנציה חישובית.

פרק 9

יישום מטריציוני ב־MATLAB

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

% ===== Step 1: Define node admittances =====

% Conductance between nodes (Siemens)

G1 = 1; G2 = 0.5; G3 = 0.25;

G4 = 1; G5 = 0.5;

% ===== Step 2: Build the admittance matrix Y =====

% Reference node is chosen separately; here 3 active nodes remain

Y = [ G1+G2 -G2 0 ;

-G2 G2+G3+G4 -G4 ;

0 -G4 G4+G5 ];

% ===== Step 3: Build the current source vector I =====

I = [ +1; -1; +1 ]; % Net injected current at each node [A]

% ===== Step 4: Solve the linear system =====

V = Y \ I; % Backslash uses LU with partial pivoting

disp('Node voltages [V]:'); disp(V);

% ===== Step 5: Cramer's rule for a single node voltage =====

k = 1; % Index of node of interest

Yk = Y; Yk(:, k) = I;

Vk = det(Yk) / det(Y);

fprintf('V_%d (via Cramer) = %.4f V\n', k, Vk);

% ===== Step 6: Branch current between nodes i and j =====

i_node = 1; j_node = 2; G_branch = G2;

i2 = G_branch * (V(i_node) - V(j_node));

fprintf('Branch current i2 = %.4f A\n', i2);

% ===== Step 7: Sanity check via KCL at each node =====

residual = Y*V - I;

fprintf('Max KCL residual = %.2e\n', max(abs(residual)));

שיקול ביצועים

ב־MATLAB, אופרטור \ (backslash) הוא הדרך המומלצת לפתור מערכות ליניאריות. הוא בוחר באופן אוטומטי את האלגוריתם המתאים על בסיס מבנה המטריצה – משולש, סימטרי, דליל, או כללי. שימוש בו לרוב מהיר ויציב יותר משימוש ב־inv(Y)*I.

הימנעו מהיפוך מפורש

כתיבת V = inv(Y) * I נראית טבעית, אך היא איטית ופחות יציבה נומרית מהאופרטור \. ההבדל הופך משמעותי במעגלים בני מאות צמתים, או כאשר מטריצת האדמיטנס מותנית גרוע.

מעבר ממעגלי DC למעגלי AC

השינוי היחיד הדרוש במעבר ל־AC הוא להחליף את הגדלים הריאליים בגדלים מרוכבים: התנגדויות הופכות לעכבות, ואדמיטנסים הופכים למספרים מרוכבים. מטריצת $$Y$$ תהיה מטריצה מרוכבת, וגם $$V$$ ו־ $$I$$ יהיו מרוכבים. MATLAB מטפל במספרים מרוכבים באופן שקוף לחלוטין – אותו הקוד יעבוד ללא שינוי.

פרק 10

הרחבות: AC, רכיבים אקטיביים ותלות בזמן

הניסוח $Y\,V=I$ הוא נקודת המוצא. המעבר למצבים מתקדמים יותר מתבצע על ידי שימור המבנה הכללי תוך הרחבת המודלים של איברי המטריצה והוקטור.

ניתוח AC הרמוני

עוברים לפזה (Phasor Domain): כל אדמיטנס הופך למרוכב, $Y_{C}=j\omega C$ עבור קבל ו־ $Y_{L}=1/(j\omega L)$ עבור סליל. הפתרון מתקבל בכל תדר $\omega$ בנפרד.

סריקת תדר

לולאה על תדרים $\omega_{1},\ldots,\omega_{m}$ , פתרון $Y(\omega)V=I$ בכל אחד מהם, ובניית תגובת תדר מלאה. כך מתקבלים גרפי Bode במישרין מהמודל המתמטי.

תחום הזמן (Transient)

דיסקרטיזציה של הקבלים והסלילים בשיטות אינטגרציה (Backward Euler, Trapezoidal). בכל צעד זמן נפתרת מערכת $$YV=I$$ ריאלית, כאשר $$I$$ מעודכן לפי המצב הקודם.

רכיבים תלויים

VCVS, CCVS, VCCS, CCCS – כל אחד מהם תורם איברים מחוץ לאלכסון של $$Y$$ , באופן השובר את הסימטריה. הניתוח הופך מורכב יותר, אך המבנה הכללי נשמר.

רכיבים לא ליניאריים

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

ניתוח רעש

כל רכיב מקבל מקור רעש שקול (תרמי, Shot, Flicker). הניתוח מתבצע באותה תשתית של $$Y$$ , ומתקבלת צפיפות ספקטרלית של מתח או זרם בכל נקודה במעגל.

פרק 11

שיקולים נומריים

האלגוריתם הכללי נשמע פשוט, אך מעבר מתורה לפועל דורש מודעות לכמה מלכודות נומריות.

מטריצת אדמיטנס סינגולרית

נגרמת בדרך כלל מצומת מבודד או מבחירת צומת ייחוס לא חוקית (Floating Node). הפותר יזרוק שגיאת Conditioning. הפתרון: ודאו רציפות חשמלית של כל המעגל ובחירה תקינה של אדמה.

טווח דינמי רחב של ערכים

עירוב של נגדים בני $1\,\Omega$ ו־ $10^{12}\,\Omega$ באותה מטריצה מוביל ל־Condition Number עצום ולאיבוד דיוק. ב־MATLAB ניתן לבדוק זאת באמצעות cond(Y).

שימוש ב־inv במקום ב־backslash

היפוך מפורש של מטריצה הוא איטי וצובר שגיאות עיגול. תמיד העדיפו את אופרטור \ או פירוק LU מפורש על inv().

איבוד הסימטריה והדלילות

המרת המטריצה לפורמט צפוף כאשר ניתן לשמור עליה דלילה גורמת לאיבוד מהירות ולשימוש מיותר בזיכרון. ב־MATLAB השתמשו בפונקציה sparse() במעגלים גדולים.

השוואת תוצאות הפותר לחישוב ידני

בדיקה תמיד מומלצת – ניתן לחשב KCL שיורי ( $$YV-I$$ ) ולוודא שהוא קרוב לאפס. אם הוא לא, כנראה שיש טעות בהרכבת $$Y$$ או ב־ $$I$$ .

פרק 12

סיכום

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

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

שבעה כללי האצבע של ניתוח מטריציוני

1. בחרו צומת ייחוס עם הכי הרבה ענפים.

2. בנו את

\(Y\)

לפי כלל אלכסון/לא־אלכסון.

3. וקטור

\(I\)

מורכב רק ממקורות חיצוניים.

4. פתרו עם backslash, לא עם inv.

5. זרם בענף = אדמיטנס × הפרש מתחים.

6. בדקו את

YV-I\approx 0

בסוף.

7. כלל קרמר טוב למעגלים קטנים בלבד.

השלב הבא

לאחר הבנת התהליך עבור DC, מומלץ להתקדם לניתוח AC הרמוני ולסריקות תדר. שם תראו כיצד אותה תשתית מטריציונית מובילה ישירות לדיאגרמות Bode, ולתכנון מסננים אנלוגיים.

מעבר למעגלי זרם חילופין