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

في الرياضيات ، وخاصة الجبر الخطي والتحليل العددي ، تعد عملية غرام-شميدت أو خوارزمية غرام-شميدت طريقة لإيجاد مجموعة من متجهين أو أكثر متعامدين مع بعضهم البعض.
بحسب التعريف التقني، هي طريقة لإنشاء أساس متعامد من مجموعة من المتجهات في فضاء الضرب الداخلي ، وغالبًا ما يكون الفضاء الإقليديمزودة بالضرب الداخلي القياسي . تأخذ عملية غرام-شميدت مجموعة محدودة ومستقلة خطيًا من المتجهاتلـ k ≤ n وتولد مجموعة متعامدةوهذا يشمل نفسفضاء فرعي ذو أبعاد منمثل.
سُميت هذه الطريقة نسبةً إلى يورغن بيدرسن غرام وإرهارد شميدت ، لكن بيير سيمون لابلاس كان على دراية بها قبل غرام وشميدت. [ 1 ] في نظرية تفكيكات زمرة لي ، يتم تعميمها بواسطة تفكيك إيواساوا .
يؤدي تطبيق عملية غرام-شميدت على متجهات الأعمدة لمصفوفة ذات رتبة عمودية كاملة إلى تحليل QR (يتم تحليلها إلى مصفوفة متعامدة ومصفوفة مثلثية ).
وصف

الإسقاط المتجهي لمتجهعلى متجه غير صفرييُعرَّف على النحو التالي [ ملاحظة 1 ] أينيرمز إلى الضرب النقطي للمتجهاتووهذا يعني أنهو الإسقاط المتعامد لـعلى الخط الذي يمتد بواسطة. لوإذا كان متجه الصفر،يُعرَّف بأنه المتجه الصفري.
منحمتجهات مستقلة خطيًا غير صفريةتحدد عملية غرام-شميدت المتجهاتعلى النحو التالي:
التسلسلهو النظام المطلوب من المتجهات المتعامدة، والمتجهات المعياريةتشكيل مجموعة متعامدة . حساب المتتاليةتُعرف هذه العملية باسم تعامد غرام-شميدت ، وحساب المتتاليةيُعرف باسم التعامد غرام-شميدت .
للتحقق من أن هذه الصيغ تُنتج متتالية متعامدة، احسب أولاًعن طريق استبدال الصيغة أعلاه بـنحصل على صفر. ثم نستخدم هذا لحسابمرة أخرى عن طريق استبدال الصيغة بـنحصل على صفر. لأي قيمة عشوائيةيتم إثبات ذلك عن طريق الاستقراء الرياضي .
هندسياً، تتم هذه الطريقة على النحو التالي: لحساب، إنها مشاريعبشكل متعامد على الفضاء الفرعيتم إنشاؤه بواسطةوهو نفس الفضاء الجزئي الناتج عنالمتجهويُعرَّف ذلك بأنه الفرق بينوهذا الإسقاط، المضمون أن يكون متعامدًا مع جميع المتجهات في الفضاء الفرعي.
تنطبق عملية غرام-شميدت أيضًا على متتالية خطية مستقلة قابلة للعد لا نهائية { vᵢ } ᵢ . والنتيجة هي متتالية متعامدة (أو متعامدة معيارية) {uᵢ} ᵢ بحيث يكون ، بالنسبة للعدد الطبيعي n : الفضاء الجبري لـهو نفسه مثل ذلك الخاص بـ.
إذا طُبقت عملية غرام-شميدت على متتالية مرتبطة خطيًا، فإنها تُخرج المتجه 0 علىالخطوة رقم 1، بافتراض أنهو مزيج خطي منإذا كان المطلوب إنتاج أساس متعامد، فيجب على الخوارزمية اختبار وجود متجهات صفرية في المخرجات وتجاهلها لأنه لا يمكن أن يكون لأي مضاعف للمتجه الصفري طول يساوي 1. وسيكون عدد المتجهات التي تنتجها الخوارزمية هو بُعد الفضاء الذي تغطيه المدخلات الأصلية.
نوع من أنواع عملية غرام-شميدت باستخدام التكرار المتسامي المطبق على سلسلة لا نهائية (ربما غير قابلة للعد) من المتجهاتينتج عنه مجموعة من المتجهات المتعامدة.معبحيث يكون لأي، إتمام فترةهو نفسه مثل ذلك الخاص بـعلى وجه الخصوص ، عند تطبيقها على أساس (جبري) لفضاء هيلبرت (أو، بشكل أعم، أساس لأي فضاء جزئي كثيف)، فإنها تُنتج أساسًا متعامدًا (تحليليًا وظيفيًا). لاحظ أنه في الحالة العامة، غالبًا ما تكون المتباينة الصارمةويبقى هذا صحيحًا، حتى لو كانت المجموعة الأولية مستقلة خطيًا، ومدىليس بالضرورة أن يكون فضاءً جزئياً من امتداد(بل هو فضاء جزئي من اكتماله).
مثال
الفضاء الإقليدي
ضع في اعتبارك مجموعة المتجهات التالية في(مع الضرب الداخلي التقليدي )
الآن، قم بإجراء اختبار غرام-شميدت للحصول على مجموعة متعامدة من المتجهات:
نتحقق من أن المتجهاتوهي متعامدة بالفعل: مع ملاحظة أنه إذا كان حاصل الضرب النقطي لمتجهين يساوي صفرًا فإنهما متعامدان.
بالنسبة للمتجهات غير الصفرية، يمكننا بعد ذلك تطبيع المتجهات عن طريق قسمة أحجامها كما هو موضح أعلاه:
ملكيات
يرمز بـنتيجة تطبيق عملية غرام-شميدت على مجموعة من المتجهاتينتج عن ذلك خريطة.
له الخصائص التالية:
- إنه مستمر
- وهو يحافظ على التوجيه بمعنى أن.
- يتوافق مع الخرائط المتعامدة:
يتركأن تكون متعامدة (بالنسبة للجداء الداخلي المعطى). عندئذٍ لدينا
علاوة على ذلك، ينتج عن نسخة مُعَلمة من عملية غرام-شميدت انكماش تشوه (قوي) للمجموعة الخطية العامةعلى المجموعة المتعامدة.
الاستقرار العددي
عند تنفيذ هذه العملية على جهاز كمبيوتر، فإن المتجهاتغالباً ما تكون هذه العمليات غير متعامدة تماماً، بسبب أخطاء التقريب . بالنسبة لعملية غرام-شميدت كما هو موضح أعلاه (والتي يشار إليها أحياناً باسم "غرام-شميدت الكلاسيكية")، فإن فقدان التعامد هذا سيء بشكل خاص؛ لذلك، يقال إن عملية غرام-شميدت (الكلاسيكية) غير مستقرة عددياً .
يمكن تثبيت عملية غرام-شميدت بتعديل بسيط؛ ويُشار إلى هذه النسخة أحيانًا باسم غرام-شميدت المعدلة أو MGS. يُعطي هذا النهج نفس نتيجة الصيغة الأصلية في الحساب الدقيق، ويُقلل الأخطاء في الحساب ذي الدقة المحدودة.
بدلاً من حساب المتجه u k كـ يتم حسابه على النحو التالي
تُستخدم هذه الطريقة في الرسوم المتحركة السابقة، عندما يكون الوسيطيُستخدم المتجه عند تعامد المتجه الأزرق.
إليك وصف آخر للخوارزمية المعدلة. بالنظر إلى المتجهاتفي خطوتنا الأولى نقوم بإنتاج المتجهاتعن طريق إزالة المكونات على طول اتجاهفي الصيغ،بعد هذه الخطوة، أصبح لدينا بالفعل متجهان متعامدان من المتجهات المطلوبة.، أيلكننا صنعنا أيضًامتعامد بالفعل معبعد ذلك، نقوم بتعامد تلك المتجهات المتبقية معوهذا يعني أننا نحسبعن طريق الطرحلقد قمنا الآن بتخزين المتجهاتحيث تكون المتجهات الثلاثة الأولى بالفعلوالمتجهات المتبقية متعامدة بالفعل معوكما هو واضح الآن، فإن الخطوة التالية هي التعامدضدوباتباع هذه الطريقة نجد المجموعة الكاملة من المتجهات المتعامدةإذا كانت المتجهات المتعامدة مطلوبة، فإننا نقوم بالتطبيع أثناء تقدمنا، بحيث تتحول المقامات في صيغ الطرح إلى واحد.
الخوارزمية
تُنفذ خوارزمية MATLAB التالية عملية التعامد الكلاسيكية لغرام-شميدت. المتجهات v1 ، ...، vk (أعمدة المصفوفة) ،V بحيث V(:,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 } كمصفوفةثم تطبيق طريقة الحذف الغاوسي على المصفوفة الموسعةسينتج ذلك المتجهات المتعامدة بدلاً منومع ذلك، فإن المصفوفةيجب تحويلها إلى شكل صفوف متدرجة ، باستخدام عملية الصف فقط المتمثلة في إضافة مضاعف عددي لصف إلى صف آخر. [ 3 ] على سبيل المثال، بأخذكما ذكرنا سابقاً، لدينا
واختزال هذا إلى شكل صفوف متدرجة ينتج
ثم تصبح المتجهات المعيارية كما في المثال أعلاه.
صيغة المحدد
يمكن التعبير عن نتيجة عملية غرام-شميدت في صيغة غير تكرارية باستخدام المحددات .
أينو، لـ،هو المحدد غرام
لاحظ أن التعبير لـهو محدد "رسمي"، أي أن المصفوفة تحتوي على كل من الكميات العددية والمتجهات؛ يتم تعريف معنى هذا التعبير على أنه نتيجة لتوسيع العامل المرافق على طول صف المتجهات.
إن صيغة المحدد لـ Gram-Schmidt أبطأ حسابيًا (بشكل أسي) من الخوارزميات المتكررة الموصوفة أعلاه؛ وهي ذات أهمية نظرية في المقام الأول.
تم التعبير عنها باستخدام الجبر الهندسي
باستخدام الرموز المستخدمة في الجبر الهندسي ، يمكن التعبير عن النتائج غير المعيارية لعملية غرام-شميدت على النحو التالي: وهو ما يعادل التعبير باستخدامالعامل المعرّف أعلاه. ويمكن التعبير عن النتائج بشكل مكافئ على النحو التالي [ 4 ] وهو ما يرتبط ارتباطًا وثيقًا بالتعبير باستخدام المحددات أعلاه.
البدائل
تستخدم خوارزميات التعامد الأخرى تحويلات هاوسهولدر أو دورانات جيفنز . وتتميز الخوارزميات التي تستخدم تحويلات هاوسهولدر باستقرار أكبر من عملية غرام-شميدت المستقرة. من ناحية أخرى، تنتج عملية غرام-شميدتالمتجه المتعامد بعدفي التكرار رقم 1، بينما ينتج عن التعامد باستخدام انعكاسات هاوسهولدر جميع المتجهات في النهاية فقط. وهذا يجعل عملية غرام-شميدت هي الطريقة الوحيدة القابلة للتطبيق على الطرق التكرارية مثل تكرار أرنولدي .
ثمة بديل آخر مدفوع باستخدام تحليل تشوليسكي لعكس مصفوفة المعادلات العادية في المربعات الصغرى الخطية .لنفترض أن لدينا مصفوفة ذات رتبة عمودية كاملة ، ويجب جعل أعمدتها متعامدة.هيرميتية وموجبة تمامًا ، لذا يمكن كتابتها على النحو التاليباستخدام تحليل تشوليسكي . المصفوفة المثلثية السفليةالمصفوفة التي تحتوي على عناصر قطرية موجبة تمامًا قابلة للعكس . عندئذٍ، تكون أعمدة المصفوفةتكون متعامدة وتغطي نفس الفضاء الجزئي الذي تغطيه أعمدة المصفوفة الأصليةالاستخدام الصريح للمنتجيجعل هذا الخوارزمية غير مستقرة، خاصةً إذا كان رقم حالة المنتج كبيرًا. ومع ذلك، تُستخدم هذه الخوارزمية عمليًا وتُطبّق في بعض حزم البرامج نظرًا لكفاءتها العالية وبساطتها.
في ميكانيكا الكم، توجد عدة مخططات للتعامد ذات خصائص أفضل من مخطط غرام-شميدت الأصلي لتطبيقات معينة. ومع ذلك، لا يزال هذا المخطط خوارزمية شائعة وفعالة حتى في حسابات البنية الإلكترونية الأكثر تعقيدًا. [ 5 ]
تعقيد وقت التشغيل
يمكن إجراء عملية التعامد باستخدام خوارزمية غرام-شميدت في وقت متعدد الحدود قوي . ويشابه تحليل وقت التشغيل تحليل عملية الحذف الغاوسي . [ 6 ] : 40
انظر أيضاً
مراجع
- ↑ تشيني الابن، إليوت وارد ؛ كينكيد، ديفيد (2009). الجبر الخطي: النظرية والتطبيقات . سودبري، ماساتشوستس: جونز وبارتليت. الصفحات 544، 558. ISBN 978-0-7637-5020-6.
- ^ جولوب وفان لون 1996 ، §5.2.8.
- ↑ بورسيل، لايل؛ تريمبل، إس واي (1 يناير 1991). "التعامد باستخدام طريقة غرام-شميدت عن طريق حذف غاوس". المجلة الرياضية الأمريكية الشهرية . 98 (6): 544-549 . doi : 10.2307/2324877 . JSTOR 2324877 .
- ↑ دورن، كريس جيه إل ؛ لاسنبي، أنتوني (2007). الجبر الهندسي للفيزيائيين . مطبعة جامعة كامبريدج. ص 124. ISBN 978-0-521-71595-9.
- ↑ بورسيل، يوكيهيرو؛ وآخرون (2011). "حسابات مبدئية لحالات الإلكترون في سلك نانوي من السيليكون يحتوي على 100,000 ذرة على حاسوب K". وقائع المؤتمر الدولي لعام 2011 للحوسبة عالية الأداء والشبكات والتخزين والتحليل . الصفحات 1:1–1:11. doi : 10.1145/2063384.2063386 . ISBN 9781450307710. S2CID 14316074 .
- ^ غروتشل، مارتن ؛ الأماكن القريبة : شريفر ، ألكسندر (1993)، الخوارزميات الهندسية والتحسين التوافقي ، الخوارزميات والتوافقيات، المجلد. 2 ( الطبعة الثانية)، Springer-Verlag، برلين، دوى : 10.1007 / 978-3-642-78240-4 ، ISBN 978-3-642-78242-8MR 1261419
ملحوظات
- ↑ في الحالة المركبة، يفترض هذا أن الضرب الداخلي خطي في المتغير الأول ومترافق خطي في المتغير الثاني. في الفيزياء، يُعدّ الاصطلاح الأكثر شيوعًا هو الخطية في المتغير الثاني، وفي هذه الحالة نُعرّف
مصادر
- باو الثالث، ديفيد؛ تريفثين، لويد ن. (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.
روابط خارجية
- "التعامد" ، موسوعة الرياضيات ، دار نشر EMS ، 2001 [1994]
- درس رياضيات من كلية هارفي مود حول خوارزمية غرام-شميدت
- أقدم الاستخدامات المعروفة لبعض كلمات الرياضيات: G يحتوي المدخل "Gram-Schmidt orthogonalization" على بعض المعلومات والمراجع حول أصول الطريقة.
- عروض توضيحية: عملية غرام شميدت في المستوى وعملية غرام شميدت في الفضاء
- تطبيق غرام-شميدت للتعامد
- روتين NAG Gram–Schmidt للتعامد لـ n متجهات من الرتبة m
- البرهان: ريموند بوزيو، كينان كيدويل. "برهان خوارزمية غرام-شميدت للتعامد" (الإصدار 8). PlanetMath.org.
- الجبر الخطي
- التحليل الوظيفي
