مكمل شور

مكمل شور هو أداة رئيسية في مجالات الجبر الخطي ، ونظرية المصفوفات ، والتحليل العددي ، والإحصاء .

يُعرَّف هذا المفهوم لمصفوفة كتلية . لنفترض أن p و q عددان صحيحان غير سالبين بحيث يكون p + q > 0 ، ولنفترض أن A و B و C و D هي على التوالي مصفوفات من الأعداد المركبة بأبعاد p × p و p × q و q × p و q × q . م=[أبجد]{\displaystyle M={\begin{bmatrix}A&B\\C&D\end{bmatrix}}} بحيث تكون M عبارة عن مصفوفة ( p + q ) × ( p + q ).

إذا كانت D قابلة للعكس، فإن مكمل شور للكتلة D من المصفوفة M هو المصفوفة p × p المعرفة بـ م/د:=أ-بد-1ج.{\displaystyle M/D:=A-BD^{-1}C.} إذا كانت المصفوفة A قابلة للعكس، فإن مكمل شور للكتلة A من المصفوفة M هو المصفوفة q × q المعرفة بـ م/أ:=د-جأ-1ب.{\displaystyle M/A:=D-CA^{-1}B.} في حالة كون A أو D مفردة ، فإن استبدال المعكوس المعمم للمعاكوس على M/A و M/D ينتج عنه مكمل شور المعمم .

سُمّي مُكمِّل شور نسبةً إلى إساي شور [ 1 ] الذي استخدمه لإثبات مُبرهنة شور ، على الرغم من استخدامه سابقًا. [ 2 ] وكانت إميلي فيرجينيا هاينسورث أول من أطلق عليه اسم مُكمِّل شور . [ 3 ] ويُشار إلى مُكمِّل شور أحيانًا باسم خريطة فيشباخ نسبةً إلى الفيزيائي هيرمان فيشباخ . [ 4 ]

خلفية

ينشأ مكمل شور عند إجراء عملية حذف غاوسي جزئي على المصفوفة M. ولحذف العناصر الموجودة أسفل القطر الرئيسي، يتم ضرب المصفوفة M بمصفوفة مثلثية سفلية جزئية على اليمين كما يلي: م=[أبجد][أبجد][أناص0-د-1جأناq]=[أ-بد-1جب0د]،{\displaystyle {\begin{aligned}&M={\begin{bmatrix}A&B\\C&D\end{bmatrix}}\quad \to \quad {\begin{bmatrix}A&B\\C&D\end{bmatrix}}{\begin{bmatrix}I_{p}&0\\-D^{-1}C&I_{q}\end{bmatrix}}={\begin{bmatrix}A-BD^{-1}C&B\\0&D\end{bmatrix}},\end{aligned}}} حيث I p تُمثل مصفوفة الوحدة من الرتبة p × p . ونتيجة لذلك، فإن مكمل شورم/د=أ-بد-1ج{\displaystyle M/D=A-BD^{-1}C}يظهر في الجزء العلوي الأيسر من المربع p × p .

مواصلة عملية الاستبعاد بعد هذه النقطة (أي إجراء عملية استبعاد جاوس-جوردان الكتلية )، [أ-بد-1جب0د][أناص-بد-10أناq][أ-بد-1جب0د]=[أ-بد-1ج00د]،\displaystyle \begin{aligned}&{\begin{bmatrix}A-BD^{-1}C&B\\0&D\end{bmatrix}}\quad \to \quad {\begin{bmatrix}I_{p}&-BD^{-1}\\0&I_{q}\end{bmatrix}}{\begin{bmatrix}A-BD^{-1}C&B\\0&D\end{bmatrix}}={\begin{bmatrix}A-BD^{-1}C&0\\0&D\end{bmatrix}},\end{aligned}}} يؤدي ذلك إلى تحليل LDU لـ M ، والذي يُقرأ كالتالي: م=[أبجد]=[أناصبد-10أناq][أ-بد-1ج00د][أناص0د-1جأناq].{\displaystyle {\begin{aligned}M&={\begin{bmatrix}A&B\\C&D\end{bmatrix}}={\begin{bmatrix}I_{p}&BD^{-1}\\0&I_{q}\end{bmatrix}}{\begin{bmatrix}A-BD^{-1}C&0\\0&D\end{bmatrix}}{\begin{bmatrix}I_{p}&0\\D^{-1}C&I_{q}\end{bmatrix}}.\end{aligned}}}وبالتالي، يمكن التعبير عن معكوس M باستخدام D 1 ومعكوس مكمل شور، بافتراض وجوده، على النحو التالي: م-1=[أبجد]-1=([أناصبد-10أناq][أ-بد-1ج00د][أناص0د-1جأناq])-1=[أناص0-د-1جأناq][(أ-بد-1ج)-100د-1][أناص-بد-10أناq]=[(أ-بد-1ج)-1-(أ-بد-1ج)-1بد-1-د-1ج(أ-بد-1ج)-1د-1+د-1ج(أ-بد-1ج)-1بد-1]=[(م/د)-1-(م/د)-1بد-1-د-1ج(م/د)-1د-1+د-1ج(م/د)-1بد-1].{\displaystyle {\begin{aligned}M^{-1}={\begin{bmatrix}A&B\\C&D\end{bmatrix}}^{-1}={}&\left({\begin{bmatrix}I_{p}&BD^{-1}\\0&I_{q}\end{bmatrix}}{\begin{bmatrix}A-BD^{-1}C&0\\0&D\end{bmatrix}}{\begin{bmatrix}I_{p}&0\ \D^{-1}C&I_{q}\end{bmatrix}}\right)^{-1}\\={}&{\begin{bmatrix}I_{p}&0\\-D^{-1}C&I_{q}\end{bmatrix}}{\begin{bmatrix}\left(A-BD^{-1}C\right)^{-1}&0\\0&D^{-1}\end{bmatrix}}{\begin{bmatrix}I_{p}&-BD^{-1} \\0&I_{q}\end{bmatrix}}\\[4pt]={}&{\begin{bmatrix}\left(A-BD^{-1}C\right)^{-1}&-\left(A-BD^{-1}C\right)^{-1}BD^{-1}\\-D^{-1}C\left(A-BD^{-1}C\right)^{-1}&D^{-1}+D^{-1}C\left(A-BD^{-1}C\right)^{-1}BD^ {-1}\end{bmatrix}}\\[4pt]={}&{\begin{bmatrix}\left(M/D\right)^{-1}&-\left(M/D\right)^{-1}BD^{-1}\\-D^{-1}C\left(M/D\right)^{-1}&D^{-1}+D^{-1}C\left(M/D\right)^{-1}BD^{-1}\end{bmatrix}}.\end{aligned}}} تنشأ العلاقة المذكورة أعلاه من عمليات الحذف التي تتضمن D 1 و M/D . ويمكن اشتقاق مكافئ مع تبديل أدوار A و D. وبمساواة تعبيرات M 1 التي تم الحصول عليها بهاتين الطريقتين المختلفتين، يمكن إثبات مبرهنة معكوس المصفوفة ، التي تربط بين مكملي شور للمصفوفة M : M/D و M/A (انظر "الاشتقاق من تحليل LDU" في متطابقة مصفوفة وودبري §  البراهين البديلة ).

ملكيات

  • إذا كان كل من p و q يساوي 1 (أي أن A و B و C و D كلها أعداد قياسية)، فإننا نحصل على الصيغة المألوفة لمعكوس مصفوفة 2×2:
م-1=1أد-بج[د-ب-جأ]{\displaystyle M^{-1}={\frac {1}{AD-BC}}\left[{\begin{matrix}D&-B\\-C&A\end{matrix}}\right]}
بشرط أن يكون الفرق بين AD و BC غير صفري.  
  • بشكل عام، إذا كانت المصفوفة A قابلة للعكس، فإن
م=[أبجد]=[أناص0جأ-1أناq][أ00د-جأ-1ب][أناصأ-1ب0أناq]،م-1=[أ-1+أ-1ب(م/أ)-1جأ-1-أ-1ب(م/أ)-1-(م/أ)-1جأ-1(م/أ)-1]{\displaystyle {\begin{aligned}M&={\begin{bmatrix}A&B\\C&D\end{bmatrix}}={\begin{bmatrix}I_{p}&0\\CA^{-1}&I_{q}\end{bmatrix}}{\begin{bmatrix}A&0\\0&D-CA^{-1}B\end{bmatrix}}{\begin{bmatrix}I_ {p}&A^{-1}B\\0&I_{q}\end{bmatrix}},\\[4pt]M^{-1}&={\begin{bmatrix}A^{-1}+A^{-1}B(M/A)^{-1}CA^{-1}&-A^{-1}B(M/A)^{-1}\\-(M/A)^{-1}CA^{-1}&(M/A)^{-1}\end{bmatrix}}\end{aligned}}}
كلما وُجد هذا العكس.
  • (صيغة شور) عندما تكون المصفوفة A ، أو D ، قابلة للعكس، يتضح أيضًا أن محدد المصفوفة M يُعطى بالصيغة التالية:
المحقق(م)=المحقق(أ)المحقق(د-جأ-1ب){\displaystyle \det(M)=\det(A)\det \left(D-CA^{-1}B\right)}، على التوالى
المحقق(م)=المحقق(د)المحقق(أ-بد-1ج){\displaystyle \det(M)=\det(D)\det \left(A-BD^{-1}C\right)}،
والتي تعمم صيغة المحدد للمصفوفات 2 × 2.
  • ( صيغة غوتمان لجمع الرتب) إذا كانت D قابلة للعكس، فإن رتبة M تُعطى بواسطة
رتبة(م)=رتبة(د)+رتبة(أ-بد-1ج){\displaystyle \operatorname {rank} (M)=\operatorname {rank} (D)+\operatorname {rank} \left(A-BD^{-1}C\right)}

تطبيق لحل المعادلات الخطية

ينشأ مكمل شور بشكل طبيعي عند حل نظام من المعادلات الخطية مثل [ 7 ]

[أبجد][xy]=[uv]{\displaystyle {\begin{bmatrix}A&B\\C&D\end{bmatrix}}{\begin{bmatrix}x\\y\end{bmatrix}}={\begin{bmatrix}u\\v\end{bmatrix}}}.

بافتراض أن المصفوفة الفرعيةأ{\displaystyle A}إذا كانت قابلة للعكس، فيمكننا حذفهاx{\displaystyle x}من المعادلات، كما يلي.

x=أ-1(u-بy).{\displaystyle x=A^{-1}(u-By).}

بإدخال هذا التعبير في المعادلة الثانية نحصل على

(د-جأ-1ب)y=v-جأ-1u.{\displaystyle \left(D-CA^{-1}B\right)y=v-CA^{-1}u.}

نشير إلى هذا باسم المعادلة المختزلة التي تم الحصول عليها عن طريق حذفx{\displaystyle x}من المعادلة الأصلية. تُسمى المصفوفة التي تظهر في المعادلة المختزلة مكمل شور للكتلة الأولىأ{\displaystyle A}فيم{\displaystyle M}:

S =دهـو د-جأ-1ب{\displaystyle S\ {\overset {\underset {\mathrm {def} }{}}{=}}\ D-CA^{-1}B}.

وبحل المعادلة المختزلة، نحصل على

y=S-1(v-جأ-1u).{\displaystyle y=S^{-1}\left(v-CA^{-1}u\right).}

وبتعويض هذا في المعادلة الأولى نحصل على

x=(أ-1+أ-1بS-1جأ-1)u-أ-1بS-1v.{\displaystyle x=\left(A^{-1}+A^{-1}BS^{-1}CA^{-1}\right)u-A^{-1}BS^{-1}v.}

يمكننا التعبير عن المعادلتين أعلاه على النحو التالي:

[xy]=[أ-1+أ-1بS-1جأ-1-أ-1بS-1-S-1جأ-1S-1][uv].{\displaystyle {\begin{bmatrix}x\\y\end{bmatrix}}={\begin{bmatrix}A^{-1}+A^{-1}BS^{-1}CA^{-1}&-A^{-1}BS^{-1}\\-S^{-1}CA^{-1}&S^{-1}\end{bmatrix}}{\begin{bmatrix}u\\v\end{bmatrix}}.}

وبالتالي، فإن صيغة معكوس مصفوفة الكتل هي:

[أبجد]-1=[أ-1+أ-1بS-1جأ-1-أ-1بS-1-S-1جأ-1S-1]=[أناص-أ-1بأناq][أ-1S-1][أناص-جأ-1أناq].{\displaystyle {\begin{bmatrix}A&B\\C&D\end{bmatrix}}^{-1}={\begin{bmatrix}A^{-1}+A^{-1}BS^{-1}CA^{-1}&-A^{-1}BS^{-1}\\-S^{-1}CA^{-1}&S^{-1}\end{bmatrix}}={\begin{bmatrix}I_{p}&-A^{-1}B\\&I_{q}\end{bmatrix}}{\begin{bmatrix}A^{-1}&\\&S^{-1}\end{bmatrix}}{\begin{bmatrix}I_{p}&\\-CA^{-1}&I_{q}\end{bmatrix}}.}

على وجه الخصوص، نرى أن مكمل شور هو معكوس مكمل شور.2،2{\displaystyle 2,2}إدخال كتلة معكوسم{\displaystyle M}.

من الناحية العملية، يحتاج المرءأ{\displaystyle A}يجب أن تكون هذه الخوارزمية دقيقة عددياً وأن تكون في حالة جيدة .

تُعدّ هذه الطريقة مفيدة في الهندسة الكهربائية لتقليل أبعاد معادلات الشبكة. وهي مفيدة بشكل خاص عندما تكون عناصر متجه الخرج صفرًا. على سبيل المثال، عندماu{\displaystyle u}أوv{\displaystyle v}إذا كانت قيمة المتغير صفرًا، فيمكننا حذف الصفوف المرتبطة بها من مصفوفة المعاملات دون أي تغييرات على بقية متجه الإخراج.v{\displaystyle v}إذا كانت القيمة صفرًا، فإن المعادلة أعلاه لـx{\displaystyle x}يتقلص إلى x=(أ-1+أ-1بS-1جأ-1)u{\displaystyle x=\left(A^{-1}+A^{-1}BS^{-1}CA^{-1}\right)u}وبالتالي تقليل أبعاد مصفوفة المعاملات مع تركu{\displaystyle u}دون تعديل. يُستخدم هذا بشكل مفيد في الهندسة الكهربائية حيث يُشار إليه باسم إزالة العقدة أو تقليل كرون .

تطبيقات في نظرية الاحتمالات والإحصاء

لنفترض أن متجهي العمود العشوائيين X و Y ينتميان إلى Rⁿ و R⁼ على التوالي ، وأن المتجه ( X , Y ) في Rⁿ⁺⁼ له توزيع طبيعي متعدد المتغيرات، ومصفوفة تباينه هي مصفوفة متماثلة موجبة التحديد .

Σ=[أببتيج]،{\displaystyle \Sigma =\left[{\begin{matrix}A&B\\B^{\mathrm {T} }&C\end{matrix}}\right],}

أين أRن×ن{\textstyle A\in \mathbb {R} ^{n\times n}}هي مصفوفة التغاير لـ X ،جRم×م{\textstyle C\in \mathbb {R} ^{m\times m}}هي مصفوفة التغاير لـ Y وبRن×م{\textstyle B\in \mathbb {R} ^{n\times m}}هي مصفوفة التغاير بين X و Y.

إذن، فإن التغاير الشرطي لـ X بشرط Y هو مكمل شور لـ C فيΣ{\textstyle \Sigma }[ 8 ]

كوف(X|Y)=أ-بج-1بتيهـ(X|Y)=هـ(X)+بج-1(Y-هـ(Y)){\displaystyle {\begin{aligned}\operatorname {Cov} (X\mid Y)&=A-BC^{-1}B^{\mathrm {T} }\\\operatorname {E} (X\mid Y)&=\operatorname {E} (X)+BC^{-1}(Y-\operatorname {E} (Y))\end{aligned}}}

إذا أخذنا المصفوفةΣ{\displaystyle \Sigma }إذا كان ما سبق ليس تباينًا لمتجه عشوائي ، بل تباينًا لعينة ، فقد يكون له توزيع ويشارت . في هذه الحالة، يكون مكمل شور لـ C فيΣ{\displaystyle \Sigma }كما أن لديها شركة توزيع ويشارت.

شروط التحديد الإيجابي والتحديد شبه الإيجابي

لتكن X مصفوفة متناظرة من الأعداد الحقيقية معطاة بـ X=[أببتيج].{\displaystyle X=\left[{\begin{matrix}A&B\\B^{\mathrm {T} }&C\end{matrix}}\right].} ثم باستخدام صيغة هاينسورث لجمع العطالة ، نجد

  • إذا كانت A قابلة للعكس، فإن X تكون موجبة التحديد إذا وفقط إذا كانت A ومتممتها X/A موجبتي التحديد: [ 2 ] : 34
X0أ0،X/أ=ج-بتيأ-1ب0.{\displaystyle X\succ 0\Leftrightarrow A\succ 0,X/A=C-B^{\mathrm {T} }A^{-1}B\succ 0.}
  • إذا كانت C قابلة للعكس، فإن X تكون موجبة التحديد إذا وفقط إذا كانت C ومكملتها X/C موجبتي التحديد:
X0ج0،X/ج=أ-بج-1بتي0.{\displaystyle X\succ 0\Leftrightarrow C\succ 0,X/C=A-BC^{-1}B^{\mathrm {T} }\succ 0.}
  • إذا كانت A موجبة تمامًا، فإن X تكون موجبة شبه تمامًا إذا وفقط إذا كانت متممتها X/A موجبة شبه تمامًا: [ 2 ] : 34
لو أ0، ثم X0X/أ=ج-بتيأ-1ب0.{\displaystyle {\text{If }}A\succ 0,{\text{ then }}X\succeq 0\Leftrightarrow X/A=C-B^{\mathrm {T} }A^{-1}B\succeq 0.}
  • إذا كانت C موجبة التحديد، فإن X تكون موجبة شبه التحديد إذا وفقط إذا كانت المتممة X/C موجبة شبه التحديد:
لو ج0، ثم X0X/ج=أ-بج-1بتي0.{\displaystyle {\text{If }}C\succ 0,{\text{ then }}X\succeq 0\Leftrightarrow X/C=A-BC^{-1}B^{\mathrm {T} }\succeq 0.}

يمكن أيضًا استنتاج العبارتين الأولى والثالثة [ 7 ] من خلال النظر في القيمة التي تُقلل الكمية uتيأu+2vتيبتيu+vتيجv،{\displaystyle u^{\mathrm {T} }Au+2v^{\mathrm {T} }B^{\mathrm {T} }u+v^{\mathrm {T} }Cv,\,} كدالة لـ v (لـ u ثابتة ).

علاوة على ذلك، بما أن [أببتيج]0[جبتيبأ]0{\displaystyle \left[{\begin{matrix}A&B\\B^{\mathrm {T} }&C\end{matrix}}\right]\succ 0\Longleftrightarrow \left[{\begin{matrix}C&B^{\mathrm {T} }\\B&A\end{matrix}}\right]\succ 0} وبالمثل بالنسبة للمصفوفات شبه المحددة الموجبة، فإن العبارة الثانية (والرابعة على التوالي) تأتي مباشرة من العبارة الأولى (والثالثة على التوالي).

يوجد أيضًا شرط كافٍ وضروري لشبه التحديد الموجب لـ X بدلالة مكمل شور المعمم. [ 2 ] تحديدًا،

  • X0أ0،ج-بتيأزب0،(أنا-أأز)ب=0{\displaystyle X\succeq 0\Leftrightarrow A\succeq 0,C-B^{\mathrm {T} }A^{g}B\succeq 0,\left(I-AA^{g}\right)B=0\,}و
  • X0ج0،أ-بجزبتي0،(أنا-ججز)بتي=0،{\displaystyle X\succeq 0\Leftrightarrow C\succeq 0,A-BC^{g}B^{\mathrm {T} }\succeq 0,\left(I-CC^{g}\right)B^{\mathrm {T} }=0,}

أينأز{\displaystyle A^{g}}يشير إلى معكوس معمّم لـأ{\displaystyle A}.

انظر أيضاً

مراجع

  1. ^ شور، ج. (1917). "Über Potenzreihen die im Inneren des Einheitskreises beschränkt sind" . J. رين يو أنجيواندتي الرياضيات . 147 : 205– 232. دوى : 10.1515/crll.1917.147.205 .
  2. 1 2 3 4 تشانغ، فوزين (2005). تشانغ، فوزين (محرر). مكمل شور وتطبيقاته . الطرق العددية والخوارزميات. المجلد 4. سبرينغر. doi : 10.1007/b105056 . ISBN  0-387-24271-6.
  3. Haynsworth, EV, "حول مكمل شور"، ملاحظات بازل الرياضية ، #BNB 20، 17 صفحة، يونيو 1968.
  4. فيشباخ، هيرمان (1958). "النظرية الموحدة للتفاعلات النووية". حوليات الفيزياء . 5 (4): 357-390 . doi : 10.1016/0003-4916(58)90007-1 .
  5. كرابتري، دوغلاس إي.؛ هاينسورث، إميلي في. (1969). "متطابقة لمكمل شور لمصفوفة" . وقائع الجمعية الرياضية الأمريكية . 22 (2): 364-366 . doi : 10.1090/S0002-9939-1969-0255573-1 . ISSN 0002-9939 . S2CID 122868483 .  
  6. ديفرييندت، كاريل (2022). "المقاومة الفعّالة تتجاوز المسافة: لابلاس، سيمبليس، ومُكمِّل شور" . الجبر الخطي وتطبيقاته . 639 : 24-49 . arXiv : 2010.04521 . doi : 10.1016/j.laa.2022.01.002 . S2CID 222272289 . 
  7. 1 2 بويد، إس. وفاندنبيرغ، إل. (2004)، "التحسين المحدب"، مطبعة جامعة كامبريدج (الملحق أ.5.5)
  8. فون ميزس، ريتشارد (1964). "الفصل الثامن.9.3". النظرية الرياضية للاحتمالات والإحصاء . دار النشر الأكاديمية. ISBN 978-1483255385.{{cite book}}عدم توافق رقم ISBN / التاريخ ( مساعدة )