# خوارزمية قيمة ذاتية

في الجبر الخطي، من أهم المسائل تصميم خوارزميات ذات كفاءة عالية ومستقرة لإيجاد القيم الممتلكة لمصفوفة. يمكن بواسطة خوارزمية القيمة الذاتية إيجاد المتجهات الممتلكة أيضًا.

## كثيرة حدود مميزة

إذا علمت مصفوفة مربعة A، فإن قيمة ممتلكة λ ومتجهها الممتلك المصاحب v هما من التعريف، زوجان خاضعان للعلاقة
    $A{\mathbf {v} }=\lambda {\mathbf {v} }\,\!,$
حيث v غير صفرية. بالمثل، $(A-\lambda I)\mathbf {v} =0$ (حيث ${\mathcal}I$ مصفوفة الوحدة) مقتضية أن $\det(A-\lambda I)=0$ . يمكن توسيع هذا المحدد إلى كثيرة حدود في λ، المعروفة بكثيرة حدود مميزة للمصفوفة A. أحد الطرق المعروفة لإيجاد القيم الممتلكة لمصفوفة صغيرة هي بإيجاد جذور كثيرة الحدود المميزة.
للأسف فإن هذه الطريقة محدودة. لا يمكن حل كثيرة حدود من الرتبة n > 4 بواسطة عدد محدد من العمليات الحسابية المتعاقبة.هناك طرق أكثر كفاءة وهي خوارزميات إيجاد الجذور لكثيرات حدود برتب أعلى.

## معاودة القوى

الفكرة الأساسية لهذه الطريقة هي اختيار اعتباطي لمتجه أولي b ومن ثم ضربه المتكرر بالمصفوفة، وبشكل معاود احتساب Ab, A²b, A³b,…. بفرض أن القيم الممتلكة مرتبة حسب قيمتها بحيث تكون λ$_1$ الأكبر، وبمتجه ذاتي مصاحب v$_1$.بعدها كل تكرار يقايس المركبة b باتجاه v$_1$ بمقدار λ$_1$، وكل اتجاه آخر بقيمة أقل (بفرض |λ$_2$| < |λ$_1$|).باستثناء مجموعة قياس صفر، لأي متجه أولي ستتقارب النتيجة إلى متجه ممتلك متوافق مع القيمة الممتلكة الغالبة.

## القيم الممتلكة لمصفوفة

في الرياضيات، خصوصًا في الجبر الخطي، تعد كثيرة الحدود المميزة أداة هامة عند وصف القيم الممتلكة لمصفوفة مربعة.

### الأنواع
#### القيم الممتلكة ذات المصفوفات 2×2

يمكن إيجاد حل تحليلي لقيم ذاتية لمصفوفات 2×2 مباشرة من الصيغة الرباعية:
    $A={\begin{bmatrix}a&b\\c&d\end{bmatrix}}$
فتكون كثيرة الحدود المميزة هي
    ${\rm {det{\begin{bmatrix}a-\lambda &b\\c&d-\lambda \end{bmatrix}}=(a-\lambda )(d-\lambda )-bc=\lambda ^{2}-(a+d)\lambda +(ad-bc)}}$
وتكون الحلول
    $\lambda ={\frac {a+d}{2}}\pm {\sqrt {{\frac {(a+d)^{2}}{4}}+bc-ad}}={\frac {a+d}{2}}\pm {\frac {\sqrt {4bc+(a-d)^{2}}}{2}}$
لاحظ أن كثيرة الحدود المميزة لمصفوفة 2×2 يمكن كتابته بدلالة الأثر ${\mathcal}{\rm {tr}}(A)=a+d$ والمحدد ${\mathcal}{\rm {det}}(A)=ad-bc$ بالشكل
    ${\rm {det}}{\begin{bmatrix}a-\lambda &b\\c&d-\lambda \end{bmatrix}}={\rm {det}}\left[A-\lambda I_{2}\right]=\lambda ^{2}-\lambda {\rm {tr}}(A)+{\rm {det}}(A)$
حيث $I_{2}$ هي مصفوفة الوحدة 2×2. لذا يمكن كتابة حلول القيم الذاتية لمصفوفة 2×2 بالصورة
    $\lambda ={\frac {1}{2}}\left({\rm {tr}}(A)\pm {\sqrt {{\rm {tr}}^{2}(A)-4{\rm {det}}(A)}}\right)$

#### قيم ممتلكة ذات مصفوفات 3×3

إذا كانت
    $A={\begin{bmatrix}a&b&c\\d&e&f\\g&h&i\end{bmatrix}}$
فإن كثيرة الحدود المميزة لـA هو
    $\det {\begin{bmatrix}a-\lambda &b&c\\d&e-\lambda &f\\g&h&i-\lambda \end{bmatrix}}=-\lambda ^{3}+\lambda ^{2}(a+e+i)+\lambda (db+gc+fh-ae-ai-ei)+(aei-afh-dbi+dch+gbf-gce).$
بالمثل يمكن كتابة كثيرة الحدود المميزة لمصفوفة 3×3 بدلالة الأثر ${\mathcal}{\rm {tr}}(A)$ و المحدد ${\mathcal}{\rm {det}}(A)$ على الصورة
    ${\rm {det}}\left[A-\lambda I_{3}\right]=-\lambda ^{3}+\lambda ^{2}{\rm {tr}}(A)+\lambda {\frac {1}{2}}\left[{\rm {tr}}(A^{2})-{\rm {tr}}^{2}(A)\right]+{\rm {det}}(A)$
حيث $I_{3}$ مصفوفة الوحدة 3×3.

#### طرق تقريبية برمجية لإيجاد القيم الممتلكة

في حال لم يكن لديك برامجيات علمية ذات دالة قيمة ممتلكة، فيما يلي مجموعة من الخطوات التي يمكنك استعمالها كشفرة عامة لمساعدتك في إنشاء برنامج لحل المسائل.
للنوع 2×2:
للنوع 3×3:

#### القيم الممتلكة لمصفوفة متماثلة 3x3

ملاحظة: لاتعمل هذه الطريقة في حال المصفوفات المفردة (ذات قيمة أو أكثر ممتلكة صفرية).

#### القيم الممتلكة والمتجهات الممتلكة لمصفوفات خاصة

بالنسبة للمصفوفات التي تحقق الشرط $A^{2}=\alpha$ يمكن كتابة صيغ خاصة للقيم الممتلكة الممكنة.
    $P_{+}={\frac {1}{2}}\left(I+{\frac {A}{\sqrt {\alpha }}}\right)$
    $P_{-}={\frac {1}{2}}\left(I-{\frac {A}{\sqrt {\alpha }}}\right)$
مع
    $AP_{+}={\sqrt {\alpha }}P_{+}\quad AP_{-}=-{\sqrt {\alpha }}P_{-}$
و
    $P_{+}P_{+}=P_{+}\quad P_{-}P_{-}=P_{-}\quad P_{+}P_{-}=P_{-}P_{+}=0$
وهذا يعطي الحل مرة أخرى للمتطابقة
    $I=P_{+}+P_{-}={\frac {1}{2}}\left(I+{\frac {A}{\sqrt {\alpha }}}\right)+{\frac {1}{2}}\left(I-{\frac {A}{\sqrt {\alpha }}}\right)$

## مراجع

* ↑ Press، William H.؛ Teukolsky، Saul A.؛ Vetterling، William T.؛ Flannery، Brian P. (1992). Numerical Recipes in C (ط. 2nd). Cambridge University Press. ISBN:0-521-43108-5.
* بوابة جبر
* بوابة رياضيات