נוסחת קרמר

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

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

==נוסחת קרמר==
כידוע מאלגברה ליניארית, למערכת משוואות ריבועית (כלומר, מספר המשתנים שווה למספר המשוואות) המיוצגת על ידי <math> \ A x=b</math>, כאשר <math>\ A</math> היא מטריצה ריבועית, ו-<math>\ b</math> הוא וקטור עמודה, קיים פתרון יחיד אם ורק אם <math> \ \det A \ne 0 </math>. 

על פי נוסחת קרמר, הרכיב ה-<math>\ k</math> של וקטור הפתרון <math>\ x</math> נתון על ידי
<math display="block">\ x_k=\frac{\det A_k}{\det A}</math>
כאשר <math> \ A_k</math> היא המטריצה המתקבלת על ידי החלפת העמודה ה-<math>\ k</math> שבמטריצה <math>\ A</math> בווקטור <math>\ b</math>.

==דוגמה==
נתונה מערכת המשוואות

<math> 
\begin{align}
\begin{cases}
 \ x \ + 2y + 3z = 2 \\
 4x +5y +6z =2 \\
 7x + 8y + 8z =4 
\end{cases}
\end{align}
</math>

כלומר המערכת מיוצגת על ידי המטריצה
<math>A=\begin{pmatrix} 1 & 2 & 3 \\ 4 & 5 & 6 \\ 7 & 8 & 8 \end{pmatrix}</math>, והווקטור <math> b=\begin{pmatrix} {\color{red}2} \\ {\color{red}2}\\ {\color{red}4} \end{pmatrix} </math>.

נחשב את הדטרמיננטות:

<math>\det A =\det\begin{pmatrix}1 & 2 & 3 \\ 4 & 5 & 6 \\ 7 & 8 & 8 \end{pmatrix}=3 </math> 

<math> \ \det A_1=\det\begin{pmatrix}{\color{red}2} & 2 & 3 \\ {\color{red}2} & 5 & 6 \\ {\color{red}4} & 8 & 8 \end{pmatrix}=-12</math>

<math>\det A_2 =\det\begin{pmatrix}1 & {\color{red}2} & 3 \\ 4 & {\color{red}2} & 6 \\ 7 & {\color{red}4} & 8 \end{pmatrix}=18 </math>

<math>\det A_3 =\det\begin{pmatrix}1 & 2 & {\color{red}2} \\ 4 & 5 & {\color{red}2} \\ 7 & 8 & {\color{red}4} \end{pmatrix}=-6 </math> 

והפתרון נתון על ידי
<math> \begin{array}{lclcl} x & = & \frac{-12}{3} & = & -4 \\

y & = & \frac{18}{3} & = & 6 \\

z & = & \frac{-6}{3} & = & -2 \end{array}
</math>

==הוכחה==
===הוכחה בעזרת התכונות של פונקציית נפח===
נניח כי נתונה המערכת <math>\ Ax=b </math>. נסמן את עמודות המטריצה ב <math>\ \alpha_1, \dots, \alpha_n</math>. הטענה כי הווקטור <math> \ (x_1,\dots,x_n)</math> פותר את המערכת היא בעצם הטענה כי
<math display="block">\ x_1 \alpha_1+\dots +x_n \alpha_n=\sum_{i=1}^{n} x_i \alpha_i =b</math>
נחשוב על הדטרמיננטה כעל פונקציית נפח, המקבלת כארגומנטים את עמודות המטריצה:
<math>\ \det A = \det(\alpha_1,\dots,\alpha_n) </math>

הדטרמיננטה <math> \ \det A_k </math> מתקבלת מהחלפת העמודה ה-k בעמודה b. כלומר:
<math display="block">\ \det A_k=\det (\alpha_1,\dots,\alpha_{k-1},b,\alpha_{k+1}\dots,\alpha_n) = \det (\alpha_1,\dots,\alpha_{k-1},\sum x_i \alpha_i,\alpha_{k+1},\dots,\alpha_n)</math>

מכיוון שהדטרמיננטה, כפונקציית נפח, היא ליניארית בכל רכיב, מתקבל

<math display="block">\ \det A_k=\sum_1^n x_i \det(\alpha_1,\dots,\alpha_{k-1},\alpha_i,\alpha_{k+1},\dots,\alpha_n) </math>

ומתכונת פונקציית הנפח, לכל <math> i\ne k</math> מתקיים כי <math> \ \det(\alpha_1,\dots,\alpha_{k-1},\alpha_i,\alpha_{k+1},\dots,\alpha_n)=0</math> ולכן נותרנו עם

<math display="block">\ \det A_k=x_k \det (\alpha_1,\dots,\alpha_{k-1},\alpha_k,\alpha_{k+1},\dots,\alpha_n)=x_k\det A</math>

ולכן <math>\ x_k =\frac{\det A_k}{\det A}</math>.

===הוכחה בעזרת הרחבת המטריצה===

נניח כי נתונה מערכת לא הומוגנית:

<div style="text-align: center;">
<math>
\begin{pmatrix}
a_{11} & \dots & a_{1n} \\
\vdots & \ddots & \vdots \\
a_{n1} & \dots & a_{nn}
\end{pmatrix}
\begin{pmatrix}
x_1 \\
\vdots\\
x_n \end{pmatrix} = \begin{pmatrix}
b_1 \\
\vdots\\
b_n \end{pmatrix} </math>
</div>
כאשר <math>\ (x_1, \dots, x_n)</math> הוא וקטור הפתרון (היחיד) של המערכת. לחישוב האיבר ה-<math>\ k</math> של <math>\ x</math>, נוסיף למערכת משוואה אחת, כך:

<div style="text-align: center;">
<math>
\begin{pmatrix}
a_{11} & \dots & \dots & \dots & a_{1n} & b_1 \\
\vdots & \ddots & \ddots & \ddots & \vdots & \vdots \\
a_{n1} & \dots & \dots & \dots & a_{nn} & b_n \\
0 & \dots & 1 & \dots & 0 & x_k 
\end{pmatrix}
\begin{pmatrix}
x_1 \\
\vdots\\
x_n\\ 
-1 \end{pmatrix} = \begin{pmatrix}
0 \\
\vdots \\
\vdots\\
0 \end{pmatrix} </math>
</div>
כאשר בשורה התחתונה במטריצה כל האיברים הם אפס, פרט לעמודה ה-<math>\ k</math> ולעמודה ה-<math>\ n+1</math>. מכיוון שהווקטור <math>\ x</math> פותר את המערכת, קל לראות על ידי העברת אגפים כי הווקטור המורחב <math> \ (x_1, x_2, \dots , x_n, -1) </math> פותר את המערכת המורחבת, שהיא כעת מערכת הומוגנית. טענה מאלגברה ליניארית אומרת כי למערכת משוואות הומוגנית יש פתרון לא טריוויאלי אם ורק אם הדטרמיננטה מתאפסת. הווקטור המורחב אינו וקטור האפס והוא פותר את המערכת, ולכן הדטרמיננטה של המטריצה המורחבת חייבת להתאפס. נפתח את הדטרמיננטה לפי השורה התחתונה, ונשים לב כי המינור ה-<math>\ k</math> הוא פשוט הדטרמיננטה של המטריצה <math>\ A_k </math> מוכפל בגורם <math>\ (-1) ^{n-k}</math>, מכיוון שכדי להביא את העמודה האחרונה למקום ה-<math>\ k</math> יש לבצע על העמודות תמורה שהיא מחזור באורך <math>\ n-k+1</math>. בנוסף, מהנוסחה לפיתוח המינור יש להכפיל בגורם <math>\ (-1)^{n+1+k}</math> ולכן מתקבל:

<math>\ (-1)^{n+1+k+n-k}\det A_k + (-1)^{n+1+n+1}x_k\det A = x_k\det A-\det A_k =0</math> 

כלומר, <math>\ x_k =\frac{\det A_k}{\det A}</math>

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

* 
* 
*

קרמר
קטגוריה:דטרמיננטות