الحذف محفوظًا
يحل حذف غاوس المعادلة Ax=b بطرح مضاعفات صف من الصفوف التي تحته حتى تصبح المصفوفة مثلثية عليا U. والمضاعِفات التي استعملتها تستحق الحفظ: ضعها تحت قطر مصفوفة مثلثية سفلى L قطرها آحاد فيكون A=LU. يسجّل هذا التحليل عملية الحذف كلها، فيُنجَز العمل مرة واحدة ويُعاد استخدامه.
A=LU,Ly=b,Ux=y(1) للحل نحل أولًا Ly=b بالتعويض الأمامي (من الأعلى إلى الأسفل) ثم Ux=y بالتعويض الخلفي (من الأسفل إلى الأعلى). وكلٌّ منهما يكلّف نحو n2 عملية، أما التحليل فيكلّف نحو 32n3. ومع كثرة الطرف الأيمن (سيناريوهات كثيرة على المصفوفة نفسها) تحلّل مرة واحدة وتدفع n2 فقط لكل سيناريو.
تحليل LU يدويًالتكن
A=(2415). المضاعِف هو
4/2=2: اطرح ضعف الصف 1 من الصف 2 فتحصل على
U=(2013) و
L=(1201). تحقّق:
LU=(2412+3)=A. لنحل
Ax=(3,9):
Ly=b يعطي
y1=3 و
y2=9−2⋅3=3؛ ثم
Ux=y يعطي
x2=1 و
2x1+1=3 فيكون
x=(1,1).
- L: مثلثية سفلى، آحاد على القطر، ومضاعِفات الحذف تحته.
- U: مثلثية عليا ناتج الحذف؛ قطرها يضم المحاور.
- التمحور: بدّل الصفوف ليكون المحور كبيرًا فيُنتج PA=LU بمصفوفة تبديل P.
المحور الصفري أو الضئيل يحتاج إلى تبديل صفوففي (0111) المحور الأول 0 فيقسم الحذف البسيط على صفر مع أن المصفوفة قابلة للقلب تمامًا. وتبديل الصفين يحل المشكلة. أما المحور الضئيل (لا الصفري) فلا يسبب انهيارًا لكنه ينتج مضاعِفات ضخمة ويفسد الدقة، لذلك تستعمل المكتبات الحقيقية التمحور دائمًا. المختبر: ما الذي يُقيَّم
- المصفوفة L15 نقطة
- المصفوفة U15 نقطة
- محدّد من U15 نقطة
- الحل بـ L وU15 نقطة
- لماذا التمحور15 نقطة
- خطوات حل LU10 نقطة
- لماذا التحليل مرة واحدة15 نقطة