عملية غرام-شميدت

الخطوتان الأوليان من عملية غرام-شميدت

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

بحسب التعريف التقني، هي طريقة لإنشاء أساس متعامد من مجموعة من المتجهات في فضاء الضرب الداخلي ، وغالبًا ما يكون الفضاء الإقليديRن{\displaystyle \mathbb {R} ^{n}}مزودة بالضرب الداخلي القياسي . تأخذ عملية غرام-شميدت مجموعة محدودة ومستقلة خطيًا من المتجهاتS={v1،...،vك}{\displaystyle S=\{\mathbf {v} _{1},\ldots ,\mathbf {v} _{k}\}}لـ kn وتولد مجموعة متعامدةS={u1،...،uك}{\displaystyle S'=\{\mathbf {u} _{1},\ldots ,\mathbf {u} _{k}\}}وهذا يشمل نفسك{\displaystyle k}فضاء فرعي ذو أبعاد منRن{\displaystyle \mathbb {R} ^{n}}مثلS{\displaystyle S}.

سُميت هذه الطريقة نسبةً إلى يورغن بيدرسن غرام وإرهارد شميدت ، لكن بيير سيمون لابلاس كان على دراية بها قبل غرام وشميدت. [ 1 ] في نظرية تفكيكات زمرة لي ، يتم تعميمها بواسطة تفكيك إيواساوا .

يؤدي تطبيق عملية غرام-شميدت على متجهات الأعمدة لمصفوفة ذات رتبة عمودية كاملة إلى تحليل QR (يتم تحليلها إلى مصفوفة متعامدة ومصفوفة مثلثية ).

وصف

يتم تنفيذ عملية غرام-شميدت المعدلة على ثلاثة متجهات غير متعامدة ومستقلة خطيًا لقاعدة لـR3{\displaystyle \mathbb {R} ^{3}}انقر على الصورة للاطلاع على التفاصيل. تم شرح التعديل في قسم الاستقرار العددي من هذه المقالة.

الإسقاط المتجهي لمتجهv{\displaystyle \mathbf {v} }على متجه غير صفريu{\displaystyle \mathbf {u} }يُعرَّف على النحو التالي [ ملاحظة 1 ]مشروعu(v)=v،uu،uu،{\displaystyle \operatorname {proj} _{\mathbf {u} }(\mathbf {v} )={\frac {\langle \mathbf {v} ,\mathbf {u} \rangle }{\langle \mathbf {u} ,\mathbf {u} \rangle }}\,\mathbf {u} ,} أينv،u{\displaystyle \langle \mathbf {v} ,\mathbf {u} \rangle }يرمز إلى الضرب النقطي للمتجهاتu{\displaystyle \mathbf {u} }وv{\displaystyle \mathbf {v} }وهذا يعني أنمشروعu(v){\displaystyle \operatorname {proj} _{\mathbf {u} }(\mathbf {v} )}هو الإسقاط المتعامد لـv{\displaystyle \mathbf {v} }على الخط الذي يمتد بواسطةu{\displaystyle \mathbf {u} }. لوu{\displaystyle \mathbf {u} }إذا كان متجه الصفر،مشروعu(v){\displaystyle \operatorname {proj} _{\mathbf {u} }(\mathbf {v} )}يُعرَّف بأنه المتجه الصفري.

منحك{\displaystyle k}متجهات مستقلة خطيًا غير صفريةv1،...،vك{\displaystyle \mathbf {v} _{1},\ldots ,\mathbf {v} _{k}}تحدد عملية غرام-شميدت المتجهاتu1،...،uك{\displaystyle \mathbf {u} _{1},\ldots ,\mathbf {u} _{k}}على النحو التالي: u1=v1،هـ1=u1u1u2=v2-مشروعu1(v2)،هـ2=u2u2u3=v3-مشروعu1(v3)-مشروعu2(v3)،هـ3=u3u3u4=v4-مشروعu1(v4)-مشروعu2(v4)-مشروعu3(v4)،هـ4=u4u4    uك=vك-ج=1ك-1مشروعuج(vك)،هـك=uكuك.// _ {2}-\operatorname {proj} _{\mathbf {u} _{1}}(\mathbf {v} _{2}),&\!\mathbf {e} _{2}&={\frac {\mathbf {u} _{2}}{\|\mathbf {u} _{2}\|}}\\\mathbf {u} _{3}&=\mathbf {v} _ {3}-\اسم المشغل {proj} _{\mathbf {u} _{1}}(\mathbf {v} _{3})-\operatorname {proj} _{\mathbf {u} _{2}}(\mathbf {v} _{3}),&\!\mathbf {e} _{3}&={\frac {\mathbf {u} _{3}}{\|\mathbf {u} _ {3}\|}}\\\mathbf {u} _{4}&=\mathbf {v} _ {4}-\operatorname {proj} _{\mathbf {u} _{1}}(\mathbf {v} _{4})-\operatorname {proj} _{\mathbf {u} _{2}}(\mathbf {v} _{4})-\operatorname {مشروع} _ {\mathbf {u} _ {3}}(\mathbf {v} _{4}),&\!\mathbf {e} _{4}&={\mathbf {u} _{4} \over \|\mathbf {u} _{4}\|}\\&{}\ \ \vdots &&{}\ \ \vdots \\\mathbf {u} _{k}&=\mathbf {v} _ {k}-\sum _{j=1}^{k-1}\operatorname {proj} _{\mathbf {u} _{j}}(\mathbf {v} _{k}),&\!\mathbf {e} _{k}&={\frac {\mathbf {u} _{k}}{\|\mathbf {u} _{k}\|}}.\end{محاذاة}}}

التسلسلu1،...،uك{\displaystyle \mathbf {u} _{1},\ldots ,\mathbf {u} _{k}}هو النظام المطلوب من المتجهات المتعامدة، والمتجهات المعياريةهـ1،...،هـك{\displaystyle \mathbf {e} _{1},\ldots ,\mathbf {e} _{k}}تشكيل مجموعة متعامدة . حساب المتتاليةu1،...،uك{\displaystyle \mathbf {u} _{1},\ldots ,\mathbf {u} _{k}}تُعرف هذه العملية باسم تعامد غرام-شميدت ، وحساب المتتاليةهـ1،...،هـك{\displaystyle \mathbf {e} _{1},\ldots ,\mathbf {e} _{k}}يُعرف باسم التعامد غرام-شميدت .

للتحقق من أن هذه الصيغ تُنتج متتالية متعامدة، احسب أولاًu1،u2{\displaystyle \langle \mathbf {u} _{1},\mathbf {u} _{2}\rangle }عن طريق استبدال الصيغة أعلاه بـu2{\displaystyle \mathbf {u} _{2}}نحصل على صفر. ثم نستخدم هذا لحسابu1،u3{\displaystyle \langle \mathbf {u} _{1},\mathbf {u} _{3}\rangle }مرة أخرى عن طريق استبدال الصيغة بـu3{\displaystyle \mathbf {u} _{3}}نحصل على صفر. لأي قيمة عشوائيةك{\displaystyle k}يتم إثبات ذلك عن طريق الاستقراء الرياضي .

هندسياً، تتم هذه الطريقة على النحو التالي: لحسابuأنا{\displaystyle \mathbf {u} _{i}}، إنها مشاريعvأنا{\displaystyle \mathbf {v} _{i}}بشكل متعامد على الفضاء الفرعييو{\displaystyle U}تم إنشاؤه بواسطةu1،...،uأنا-1{\displaystyle \mathbf {u} _{1},\ldots ,\mathbf {u} _{i-1}}وهو نفس الفضاء الجزئي الناتج عنv1،...،vأنا-1{\displaystyle \mathbf {v} _{1},\ldots ,\mathbf {v} _{i-1}}المتجهuأنا{\displaystyle \mathbf {u} _{i}}ويُعرَّف ذلك بأنه الفرق بينvأنا{\displaystyle \mathbf {v} _{i}}وهذا الإسقاط، المضمون أن يكون متعامدًا مع جميع المتجهات في الفضاء الفرعييو{\displaystyle U}.

تنطبق عملية غرام-شميدت أيضًا على متتالية خطية مستقلة قابلة للعد لا نهائية { vᵢ } . والنتيجة هي متتالية متعامدة (أو متعامدة معيارية) {uᵢ} بحيث يكون ، بالنسبة للعدد الطبيعي n : الفضاء الجبري لـv1،...،vن{\displaystyle \mathbf {v} _{1},\ldots ,\mathbf {v} _{n}}هو نفسه مثل ذلك الخاص بـu1،...،uن{\displaystyle \mathbf {u} _{1},\ldots ,\mathbf {u} _{n}}.

إذا طُبقت عملية غرام-شميدت على متتالية مرتبطة خطيًا، فإنها تُخرج المتجه 0 علىأنا{\displaystyle i}الخطوة رقم 1، بافتراض أنvأنا{\displaystyle \mathbf {v} _{i}}هو مزيج خطي منv1،...،vأنا-1{\displaystyle \mathbf {v} _{1},\ldots ,\mathbf {v} _{i-1}}إذا كان المطلوب إنتاج أساس متعامد، فيجب على الخوارزمية اختبار وجود متجهات صفرية في المخرجات وتجاهلها لأنه لا يمكن أن يكون لأي مضاعف للمتجه الصفري طول يساوي 1. وسيكون عدد المتجهات التي تنتجها الخوارزمية هو بُعد الفضاء الذي تغطيه المدخلات الأصلية.

نوع من أنواع عملية غرام-شميدت باستخدام التكرار المتسامي المطبق على سلسلة لا نهائية (ربما غير قابلة للعد) من المتجهات(vα)α<λ{\displaystyle (v_{\alpha })_{\alpha <\lambda }}ينتج عنه مجموعة من المتجهات المتعامدة.(uα)α<κ{\displaystyle (u_{\alpha })_{\alpha <\kappa }}معκλ{\displaystyle \kappa \leq \lambda }بحيث يكون لأيαλ{\displaystyle \alpha \leq \lambda }، إتمام فترة{uβ:β<مين(α،κ)}{\displaystyle \{u_{\beta }:\beta <\min(\alpha ,\kappa )\}}هو نفسه مثل ذلك الخاص بـ{vβ:β<α}{\displaystyle \{v_{\beta }:\beta <\alpha \}}على وجه الخصوص ، عند تطبيقها على أساس (جبري) لفضاء هيلبرت (أو، بشكل أعم، أساس لأي فضاء جزئي كثيف)، فإنها تُنتج أساسًا متعامدًا (تحليليًا وظيفيًا). لاحظ أنه في الحالة العامة، غالبًا ما تكون المتباينة الصارمةκ<λ{\displaystyle \kappa <\lambda }ويبقى هذا صحيحًا، حتى لو كانت المجموعة الأولية مستقلة خطيًا، ومدى(uα)α<κ{\displaystyle (u_{\alpha })_{\alpha <\kappa }}ليس بالضرورة أن يكون فضاءً جزئياً من امتداد(vα)α<λ{\displaystyle (v_{\alpha })_{\alpha <\lambda }}(بل هو فضاء جزئي من اكتماله).

مثال

الفضاء الإقليدي

ضع في اعتبارك مجموعة المتجهات التالية فيR2{\displaystyle \mathbb {R} ^{2}}(مع الضرب الداخلي التقليدي ) S={v1=[31]،v2=[22]}.{\displaystyle S=\left\{\mathbf {v} _{1}={\begin{bmatrix}3\\1\end{bmatrix}},\mathbf {v} _{2}={\begin{bmatrix}2\\2\end{bmatrix}}\right\}.}

الآن، قم بإجراء اختبار غرام-شميدت للحصول على مجموعة متعامدة من المتجهات: u1=v1=[31]{\displaystyle \mathbf {u} _{1}=\mathbf {v} _{1}={\begin{bmatrix}3\\1\end{bmatrix}}}u2=v2-مشروعu1(v2)=[22]-مشروع[31][22]=[22]-810[31]=[-2/56/5].{\displaystyle \mathbf {u} _{2}=\mathbf {v} _{2}-\operatorname {proj} _{\mathbf {u} _{1}}(\mathbf {v} _{2})={\begin{bmatrix}2\\2\end{bmatrix}}-\operatorname {proj} _{\left[{\begin{smallmatrix}3\\1\end{smallmatrix}}\right]}{\begin{bmatrix}2\\2\end{bmatrix}}={\begin{bmatrix}2\\2\end{bmatrix}}-{\frac {8}{10}}{\begin{bmatrix}3\\1\end{bmatrix}}={\begin{bmatrix}-2/5\\6/5\end{bmatrix}}.}

نتحقق من أن المتجهاتu1{\displaystyle \mathbf {u} _{1}}وu2{\displaystyle \mathbf {u} _{2}}هي متعامدة بالفعل: u1،u2=[31]،[-2/56/5]=-65+65=0،{\displaystyle \langle \mathbf {u} _{1},\mathbf {u} _{2}\rangle =\left\langle {\begin{bmatrix}3\\1\end{bmatrix}},{\begin{bmatrix}-2/5\\6/5\end{bmatrix}}\right\rangle =-{\frac {6}{5}}+{\frac {6}{5}}=0,} مع ملاحظة أنه إذا كان حاصل الضرب النقطي لمتجهين يساوي صفرًا فإنهما متعامدان.

بالنسبة للمتجهات غير الصفرية، يمكننا بعد ذلك تطبيع المتجهات عن طريق قسمة أحجامها كما هو موضح أعلاه: هـ1=110[31]{\displaystyle \mathbf {e} _{1}={\frac {1}{\sqrt {10}}}{\begin{bmatrix}3\\1\end{bmatrix}}}هـ2=14025[-2/56/5]=110[-13].{\displaystyle \mathbf {e} _{2}={\frac {1}{\sqrt {40 \over 25}}}{\begin{bmatrix}-2/5\\6/5\end{bmatrix}}={\frac {1}{\sqrt {10}}}{\begin{bmatrix}-1\\3\end{bmatrix}}.}

ملكيات

يرمز بـجي إس(v1،...،vك){\displaystyle \operatorname {GS} (\mathbf {v} _{1},\dots ,\mathbf {v} _{k})}نتيجة تطبيق عملية غرام-شميدت على مجموعة من المتجهاتv1،...،vك{\displaystyle \mathbf {v} _{1},\dots ,\mathbf {v} _{k}}ينتج عن ذلك خريطةجي إس:(Rن)ك(Rن)ك{\displaystyle \operatorname {GS} \colon (\mathbb {R} ^{n})^{k}\to (\mathbb {R} ^{n})^{k}}.

له الخصائص التالية:

  • إنه مستمر
  • وهو يحافظ على التوجيه بمعنى أنأو(v1،...،vك)=أو(جي إس(v1،...،vك)){\displaystyle \operatorname {or} (\mathbf {v} _{1},\dots ,\mathbf {v} _{k})=\operatorname {or} (\operatorname {GS} (\mathbf {v} _{1},\dots ,\mathbf {v} _{k}))}.
  • يتوافق مع الخرائط المتعامدة:

يتركز:RنRن{\displaystyle g\colon \mathbb {R} ^{n}\to \mathbb {R} ^{n}}أن تكون متعامدة (بالنسبة للجداء الداخلي المعطى). عندئذٍ لدينا جي إس(ز(v1)،...،ز(vك))=(ز(جي إس(v1،...،vك)1)،...،ز(جي إس(v1،...،vك)ك)){\displaystyle \operatorname {GS} (g(\mathbf {v} _{1}),\dots ,g(\mathbf {v} _{k}))=\left(g(\operatorname {GS} (\mathbf {v} _{1},\dots ,\mathbf {v} _{k})_{1}),\dots ,g(\operatorname {GS} (\mathbf {v} _{1},\dots ,\mathbf {v} _{k})_{k})\right)}

علاوة على ذلك، ينتج عن نسخة مُعَلمة من عملية غرام-شميدت انكماش تشوه (قوي) للمجموعة الخطية العامةجيل(Rن){\displaystyle \mathrm {GL} (\mathbb {R} ^{n})}على المجموعة المتعامدةيا(Rن){\displaystyle O(\mathbb {R} ^{n})}.

الاستقرار العددي

عند تنفيذ هذه العملية على جهاز كمبيوتر، فإن المتجهاتuك{\displaystyle \mathbf {u} _{k}}غالباً ما تكون هذه العمليات غير متعامدة تماماً، بسبب أخطاء التقريب . بالنسبة لعملية غرام-شميدت كما هو موضح أعلاه (والتي يشار إليها أحياناً باسم "غرام-شميدت الكلاسيكية")، فإن فقدان التعامد هذا سيء بشكل خاص؛ لذلك، يقال إن عملية غرام-شميدت (الكلاسيكية) غير مستقرة عددياً .

يمكن تثبيت عملية غرام-شميدت بتعديل بسيط؛ ويُشار إلى هذه النسخة أحيانًا باسم غرام-شميدت المعدلة أو MGS. يُعطي هذا النهج نفس نتيجة الصيغة الأصلية في الحساب الدقيق، ويُقلل الأخطاء في الحساب ذي الدقة المحدودة.

بدلاً من حساب المتجه u k كـ uك=vك-مشروعu1(vك)-مشروعu2(vك)--مشروعuك-1(vك)،{\displaystyle \mathbf {u} _{k}=\mathbf {v} _{k}-\operatorname {proj} _{\mathbf {u} _{1}}(\mathbf {v} _{k})-\operatorname {proj} _{\mathbf {u} _{2}}(\mathbf {v} _{k})-\cdots -\operatorname {proj} _{\mathbf {u} _{k-1}}(\mathbf {v} _{k}),} يتم حسابه على النحو التالي uك(1)=vك-مشروعu1(vك)،uك(2)=uك(1)-مشروعu2(uك(1))،uك(ك-2)=uك(ك-3)-مشروعuك-2(uك(ك-3))،uك(ك-1)=uك(ك-2)-مشروعuك-1(uك(ك-2))،هـك=uك(ك-1)uك(ك-1){\displaystyle {\begin{aligned}\mathbf {u} _{k}^{(1)}&=\mathbf {v} _{k}-\operatorname {proj} _{\mathbf {u} _{1}}(\mathbf {v} _{k}),\\\mathbf {u} _{k}^{(2)}&=\mathbf {u} _{k}^{(1)}-\operatorname {proj} _{\mathbf {u} _{2}}\left(\mathbf {u} _{k}^{(1)}\right),\\&\;\;\vdots \\\mathbf {u} _{k}^{(k-2)}&=\mathbf {u} _{k}^{(k-3)}-\operatorname {proj} _{\mathbf {u} _{k-2}}\left(\mathbf {u} _{k}^{(k-3)}\right),\\\mathbf {u} _{k}^{(k-1)}&=\mathbf {u} _{k}^{(k-2)}-\operatorname {proj} _{\mathbf {u} _{k-1}}\left(\mathbf {u} _{k}^{(k-2)}\right),\\\mathbf {e} _{k}&={\frac {\mathbf {u} _{k}^{(k-1)}}{\left\|\mathbf {u} _{k}^{(k-1)}\right\|}}\end{aligned}}}

تُستخدم هذه الطريقة في الرسوم المتحركة السابقة، عندما يكون الوسيطv3{\displaystyle \mathbf {v} '_{3}}يُستخدم المتجه عند تعامد المتجه الأزرقv3{\displaystyle \mathbf {v} _{3}}.

إليك وصف آخر للخوارزمية المعدلة. بالنظر إلى المتجهاتv1،v2،...،vن{\displaystyle \mathbf {v} _{1},\mathbf {v} _{2},\dots ,\mathbf {v} _{n}}في خطوتنا الأولى نقوم بإنتاج المتجهاتv1،v2(1)،...،vن(1){\displaystyle \mathbf {v} _{1},\mathbf {v} _{2}^{(1)},\dots ,\mathbf {v} _{n}^{(1)}}عن طريق إزالة المكونات على طول اتجاهv1{\displaystyle \mathbf {v} _{1}}في الصيغ،vك(1):=vك-vك،v1v1،v1v1{\displaystyle \mathbf {v} _{k}^{(1)}:=\mathbf {v} _{k}-{\frac {\langle \mathbf {v} _{k},\mathbf {v} _{1}\rangle }{\langle \mathbf {v} _{1},\mathbf {v} _{1}\rangle }}\mathbf {v} _{1}}بعد هذه الخطوة، أصبح لدينا بالفعل متجهان متعامدان من المتجهات المطلوبة.u1،...،uن{\displaystyle \mathbf {u} _{1},\dots ,\mathbf {u} _{n}}، أيu1=v1،u2=v2(1){\displaystyle \mathbf {u} _{1}=\mathbf {v} _{1},\mathbf {u} _{2}=\mathbf {v} _{2}^{(1)}}لكننا صنعنا أيضًاv3(1)،...،vن(1){\displaystyle \mathbf {v} _{3}^{(1)},\dots ,\mathbf {v} _{n}^{(1)}}متعامد بالفعل معu1{\displaystyle \mathbf {u} _{1}}بعد ذلك، نقوم بتعامد تلك المتجهات المتبقية معu2=v2(1){\displaystyle \mathbf {u} _{2}=\mathbf {v} _{2}^{(1)}}وهذا يعني أننا نحسبv3(2)،v4(2)،...،vن(2){\displaystyle \mathbf {v} _{3}^{(2)},\mathbf {v} _{4}^{(2)},\dots ,\mathbf {v} _{n}^{(2)}}عن طريق الطرحvك(2):=vك(1)-vك(1)،u2u2،u2u2{\displaystyle \mathbf {v} _{k}^{(2)}:=\mathbf {v} _{k}^{(1)}-{\frac {\langle \mathbf {v} _{k}^{(1)},\mathbf {u} _{2}\rangle }{\langle \mathbf {u} _{2},\mathbf {u} _{2}\rangle }}\mathbf {u} _{2}}لقد قمنا الآن بتخزين المتجهاتv1،v2(1)،v3(2)،v4(2)،...،vن(2){\displaystyle \mathbf {v} _{1},\mathbf {v} _{2}^{(1)},\mathbf {v} _{3}^{(2)},\mathbf {v} _{4}^{(2)},\dots ,\mathbf {v} _{n}^{(2)}}حيث تكون المتجهات الثلاثة الأولى بالفعلu1،u2،u3{\displaystyle \mathbf {u} _{1},\mathbf {u} _{2},\mathbf {u} _{3}}والمتجهات المتبقية متعامدة بالفعل معu1،u2{\displaystyle \mathbf {u} _{1},\mathbf {u} _{2}}وكما هو واضح الآن، فإن الخطوة التالية هي التعامدv4(2)،...،vن(2){\displaystyle \mathbf {v} _{4}^{(2)},\dots ,\mathbf {v} _{n}^{(2)}}ضدu3=v3(2){\displaystyle \mathbf {u} _{3}=\mathbf {v} _{3}^{(2)}}وباتباع هذه الطريقة نجد المجموعة الكاملة من المتجهات المتعامدةu1،...،uن{\displaystyle \mathbf {u} _{1},\dots ,\mathbf {u} _{n}}إذا كانت المتجهات المتعامدة مطلوبة، فإننا نقوم بالتطبيع أثناء تقدمنا، بحيث تتحول المقامات في صيغ الطرح إلى واحد.

الخوارزمية

تُنفذ خوارزمية MATLAB التالية عملية التعامد الكلاسيكية لغرام-شميدت. المتجهات v1 ، ...، vk (أعمدة المصفوفة) ،V بحيث V(:,j)تكونج{\displaystyle j}يتم استبدال المتجهات (المتجه رقم 1) بمتجهات متعامدة (أعمدة من U) التي تمتد على نفس الفضاء الفرعي.

دالة U = غرامشميدت ( V )[ n , k ] = size ( V );U = zeros ( n , k );U (:، 1 ) = V (:، 1 ) / القاعدة ( V (:، 1 ))؛لكل i = 2 : kU (:, i ) = V (:, i );لكل j = 1 : i - 1U (:, i ) = U (:, i ) - ( U (:, j ) '* U (:, i )) * U (:, j );نهايةU (:, i ) = U (:, i ) / norm ( U (:, i ));نهايةنهاية

تبلغ تكلفة هذه الخوارزمية تقاربياً O( nk² ) من عمليات الفاصلة العائمة، حيث n هو بُعد المتجهات. [ 2 ]

عن طريق الحذف الغاوسي

إذا تم كتابة الصفوف { v 1 , ..., v k } كمصفوفةأ{\displaystyle A}ثم تطبيق طريقة الحذف الغاوسي على المصفوفة الموسعة[أأتي|أ]{\displaystyle \left[AA^{\mathsf {T}}|A\right]}سينتج ذلك المتجهات المتعامدة بدلاً منأ{\displaystyle A}ومع ذلك، فإن المصفوفةأأتي{\displaystyle AA^{\mathsf {T}}}يجب تحويلها إلى شكل صفوف متدرجة ، باستخدام عملية الصف فقط المتمثلة في إضافة مضاعف عددي لصف إلى صف آخر. [ 3 ] على سبيل المثال، بأخذv1=[31]،v2=[22]{\displaystyle \mathbf {v} _{1}={\begin{bmatrix}3&1\end{bmatrix}},\mathbf {v} _{2}={\begin{bmatrix}2&2\end{bmatrix}}}كما ذكرنا سابقاً، لدينا [أأتي|أ]=[108318822]{\displaystyle \left[AA^{\mathsf {T}}|A\right]=\left[{\begin{array}{rr|rr}10&8&3&1\\8&8&2&2\end{array}}\right]}

واختزال هذا إلى شكل صفوف متدرجة ينتج [10.80.3101-0.250.75]{\displaystyle \left[{\begin{array}{rr|rr}1&.8&.3&.1\\0&1&-.25&.75\end{array}}\right]}

ثم تصبح المتجهات المعيارية هـ1=10.32+12[0.31]=110[31]{\displaystyle \mathbf {e} _{1}={\frac {1}{\sqrt {.3^{2}+.1^{2}}}}{\begin{bmatrix}.3&.1\end{bmatrix}}={\frac {1}{\sqrt {10}}}{\begin{bmatrix}3&1\end{bmatrix}}}هـ2=10.252+0.752[-0.250.75]=110[-13]،{\displaystyle \mathbf {e} _{2}={\frac {1}{\sqrt {.25^{2}+.75^{2}}}}{\begin{bmatrix}-.25&.75\end{bmatrix}}={\frac {1}{\sqrt {10}}}{\begin{bmatrix}-1&3\end{bmatrix}},} كما في المثال أعلاه.

صيغة المحدد

يمكن التعبير عن نتيجة عملية غرام-شميدت في صيغة غير تكرارية باستخدام المحددات .

هـج=1دج-1دج|v1،v1v2،v1vج،v1v1،v2v2،v2vج،v2v1،vج-1v2،vج-1vج،vج-1v1v2vج|{\displaystyle \mathbf {e} _{j}={\frac {1}{\sqrt {D_{j-1}D_{j}}}}{\begin{vmatrix}\langle \mathbf {v} _{1},\mathbf {v} _{1}\rangle &\langle \mathbf {v} _{2},\mathbf {v} _{1}\rangle &\cdots &\langle \mathbf {v} _{j},\mathbf {v} _{1}\rangle \\\langle \mathbf {v} _{1},\mathbf {v} _{2}\rangle &\langle \mathbf {v} _{2},\mathbf {v} _{2}\rangle &\cdots &\langle \mathbf {v} _{j},\mathbf {v} _{2}\rangle \\\vdots &\vdots &\ddots &\vdots \\\langle \mathbf {v} _{1},\mathbf {v} _{j-1}\rangle &\langle \mathbf {v} _{2},\mathbf {v} _{j-1}\rangle &\cdots &\langle \mathbf {v} _{j},\mathbf {v} _{j-1}\rangle \\\mathbf {v} _{1}&\mathbf {v} _{2}&\cdots &\mathbf {v} _{j}\end{vmatrix}}}

uج=1دج-1|v1،v1v2،v1vج،v1v1،v2v2،v2vج،v2v1،vج-1v2،vج-1vج،vج-1v1v2vج|{\displaystyle \mathbf {u} _{j}={\frac {1}{D_{j-1}}}{\begin{vmatrix}\langle \mathbf {v} _{1},\mathbf {v} _{1}\rangle &\langle \mathbf {v} _{2},\mathbf {v} _{1}\rangle &\cdots &\langle \mathbf {v} _{j},\mathbf {v} _{1}\rangle \\\langle \mathbf {v} _{1},\mathbf {v} _{2}\rangle &\langle \mathbf {v} _{2},\mathbf {v} _{2}\rangle &\cdots &\langle \mathbf {v} _{j},\mathbf {v} _{2}\rangle \\\vdots &\vdots &\ddots &\vdots \\\langle \mathbf {v} _{1},\mathbf {v} _{j-1}\rangle &\langle \mathbf {v} _{2},\mathbf {v} _{j-1}\rangle &\cdots &\langle \mathbf {v} _{j},\mathbf {v} _{j-1}\rangle \\\mathbf {v} _{1}&\mathbf {v} _{2}&\cdots &\mathbf {v} _{j}\end{vmatrix}}}

أيند0=1{\displaystyle D_{0}=1}و، لـج1{\displaystyle j\geq 1}،دج{\displaystyle D_{j}}هو المحدد غرام

دج=|v1،v1v2،v1vج،v1v1،v2v2،v2vج،v2v1،vجv2،vجvج،vج|.{\displaystyle D_{j}={\begin{vmatrix}\langle \mathbf {v} _{1},\mathbf {v} _{1}\rangle &\langle \mathbf {v} _{2},\mathbf {v} _{1}\rangle &\cdots &\langle \mathbf {v} _{j},\mathbf {v} _{1}\rangle \\\langle \mathbf {v} _{1},\mathbf {v} _{2}\rangle &\langle \mathbf {v} _{2},\mathbf {v} _{2}\rangle &\cdots &\langle \mathbf {v} _{j},\mathbf {v} _{2}\rangle \\\vdots &\vdots &\ddots &\vdots \\\langle \mathbf {v} _{1},\mathbf {v} _{j}\rangle &\langle \mathbf {v} _{2},\mathbf {v} _{j}\rangle &\cdots &\langle \mathbf {v} _{j},\mathbf {v} _{j}\rangle \end{vmatrix}}.}

لاحظ أن التعبير لـuك{\displaystyle \mathbf {u} _{k}}هو محدد "رسمي"، أي أن المصفوفة تحتوي على كل من الكميات العددية والمتجهات؛ يتم تعريف معنى هذا التعبير على أنه نتيجة لتوسيع العامل المرافق على طول صف المتجهات.

إن صيغة المحدد لـ Gram-Schmidt أبطأ حسابيًا (بشكل أسي) من الخوارزميات المتكررة الموصوفة أعلاه؛ وهي ذات أهمية نظرية في المقام الأول.

تم التعبير عنها باستخدام الجبر الهندسي

باستخدام الرموز المستخدمة في الجبر الهندسي ، يمكن التعبير عن النتائج غير المعيارية لعملية غرام-شميدت على النحو التالي: uك=vك-ج=1ك-1(vكuج)uج-1 ،{\displaystyle \mathbf {u} _{k}=\mathbf {v} _{k}-\sum _{j=1}^{k-1}(\mathbf {v} _{k}\cdot \mathbf {u} _{j})\mathbf {u} _{j}^{-1}\ ,} وهو ما يعادل التعبير باستخداممشروع{\displaystyle \operatorname {proj} }العامل المعرّف أعلاه. ويمكن التعبير عن النتائج بشكل مكافئ على النحو التالي [ 4 ]uك=vكvك-1v1(vك-1v1)-1،{\displaystyle \mathbf {u} _{k}=\mathbf {v} _{k}\wedge \mathbf {v} _{k-1}\wedge \cdot \cdot \cdot \wedge \mathbf {v} _{1}(\mathbf {v} _{k-1}\wedge \cdot \cdot \cdot \wedge \mathbf {v} _{1})^{-1},} وهو ما يرتبط ارتباطًا وثيقًا بالتعبير باستخدام المحددات أعلاه.

البدائل

تستخدم خوارزميات التعامد الأخرى تحويلات هاوسهولدر أو دورانات جيفنز . وتتميز الخوارزميات التي تستخدم تحويلات هاوسهولدر باستقرار أكبر من عملية غرام-شميدت المستقرة. من ناحية أخرى، تنتج عملية غرام-شميدتج{\displaystyle j}المتجه المتعامد بعدج{\displaystyle j}في التكرار رقم 1، بينما ينتج عن التعامد باستخدام انعكاسات هاوسهولدر جميع المتجهات في النهاية فقط. وهذا يجعل عملية غرام-شميدت هي الطريقة الوحيدة القابلة للتطبيق على الطرق التكرارية مثل تكرار أرنولدي .

ثمة بديل آخر مدفوع باستخدام تحليل تشوليسكي لعكس مصفوفة المعادلات العادية في المربعات الصغرى الخطية .V{\displaystyle V}لنفترض أن لدينا مصفوفة ذات رتبة عمودية كاملة ، ويجب جعل أعمدتها متعامدة.V*V{\displaystyle V^{*}V}هيرميتية وموجبة تمامًا ، لذا يمكن كتابتها على النحو التاليV*V=لل*،{\displaystyle V^{*}V=LL^{*},}باستخدام تحليل تشوليسكي . المصفوفة المثلثية السفليةل{\displaystyle L}المصفوفة التي تحتوي على عناصر قطرية موجبة تمامًا قابلة للعكس . عندئذٍ، تكون أعمدة المصفوفةيو=V(ل-1)*{\displaystyle U=V\left(L^{-1}\right)^{*}}تكون متعامدة وتغطي نفس الفضاء الجزئي الذي تغطيه أعمدة المصفوفة الأصليةV{\displaystyle V}الاستخدام الصريح للمنتجV*V{\displaystyle V^{*}V}يجعل هذا الخوارزمية غير مستقرة، خاصةً إذا كان رقم حالة المنتج كبيرًا. ومع ذلك، تُستخدم هذه الخوارزمية عمليًا وتُطبّق في بعض حزم البرامج نظرًا لكفاءتها العالية وبساطتها.

في ميكانيكا الكم، توجد عدة مخططات للتعامد ذات خصائص أفضل من مخطط غرام-شميدت الأصلي لتطبيقات معينة. ومع ذلك، لا يزال هذا المخطط خوارزمية شائعة وفعالة حتى في حسابات البنية الإلكترونية الأكثر تعقيدًا. [ 5 ]

تعقيد وقت التشغيل

يمكن إجراء عملية التعامد باستخدام خوارزمية غرام-شميدت في وقت متعدد الحدود قوي . ويشابه تحليل وقت التشغيل تحليل عملية الحذف الغاوسي . [ 6 ] : 40

انظر أيضاً

مراجع

  1. تشيني الابن، إليوت وارد ؛ كينكيد، ديفيد (2009). الجبر الخطي: النظرية والتطبيقات . سودبري، ماساتشوستس: جونز وبارتليت. الصفحات  544، 558. ISBN 978-0-7637-5020-6.
  2. ^ جولوب وفان لون 1996 ، §5.2.8.
  3. بورسيل، لايل؛ تريمبل، إس واي (1 يناير 1991). "التعامد باستخدام طريقة غرام-شميدت عن طريق حذف غاوس". المجلة الرياضية الأمريكية الشهرية . 98 (6): 544-549 . doi : 10.2307/2324877 . JSTOR 2324877 . 
  4. دورن، كريس جيه إل ؛ لاسنبي، أنتوني (2007). الجبر الهندسي للفيزيائيين . مطبعة جامعة كامبريدج. ص 124. ISBN  978-0-521-71595-9.
  5. بورسيل، يوكيهيرو؛ وآخرون (2011). "حسابات مبدئية لحالات الإلكترون في سلك نانوي من السيليكون يحتوي على 100,000 ذرة على حاسوب K". وقائع المؤتمر الدولي لعام 2011 للحوسبة عالية الأداء والشبكات والتخزين والتحليل . الصفحات 1:1–1:11. doi : 10.1145/2063384.2063386 . ISBN   9781450307710. S2CID 14316074 . 
  6. ^ غروتشل، مارتن ؛ الأماكن القريبة : شريفر ، ألكسندر (1993)، الخوارزميات الهندسية والتحسين التوافقي ، الخوارزميات والتوافقيات، المجلد. 2 ( الطبعة الثانية)، Springer-Verlag، برلين، دوى : 10.1007 / 978-3-642-78240-4 ، ISBN   978-3-642-78242-8MR 1261419 

ملحوظات

  1. في الحالة المركبة، يفترض هذا أن الضرب الداخلي خطي في المتغير الأول ومترافق خطي في المتغير الثاني. في الفيزياء، يُعدّ الاصطلاح الأكثر شيوعًا هو الخطية في المتغير الثاني، وفي هذه الحالة نُعرّفمشروعu(v)=u،vu،uu.{\displaystyle \operatorname {proj} _{\mathbf {u} }(\mathbf {v} )={\frac {\langle \mathbf {u} ,\mathbf {v} \rangle }{\langle \mathbf {u} ,\mathbf {u} \rangle }}\,\mathbf {u} .}

مصادر

  • باو الثالث، ديفيد؛ تريفثين، لويد ن. (1997)، الجبر الخطي العددي ، فيلادلفيا: جمعية الرياضيات الصناعية والتطبيقية، ISBN 978-0-89871-361-9.
  • جولوب، جين هـفان لون، تشارلز ف. (1996)، حسابات المصفوفات (  الطبعة الثالثة)، جونز هوبكنز، ISBN 978-0-8018-5414-9.
  • غرويب، فيرنر (1975)، الجبر الخطي (  الطبعة الرابعة)، سبرينغر.
  • سوليفيريز، سي إي؛ غاليانو، إي. (1985)، "التعامد على المستوى: منهج هندسي" (ملف PDF) ، مجلة الفيزياء المكسيكية ، 31 (4): 743-758 ، مؤرشف من النسخة الأصلية (ملف PDF) بتاريخ 2014-03-07 ، تم استرجاعه بتاريخ 2013-06-22.