تقسيم مصفوفة

في القواعد الرياضية للجبر الخطي، تقسيم المصفوفة هو تعبير يستخدم لتمثيل مصفوفة ما في صورة جمع أو طرح عدة مصفوفات. وتعتمد العديد من الطرق التكرارية (مثل نظام المعادلات التفاضلية) على حل معادلات المصفوفة مباشرة.

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

نحن نسعي لحل المعادلة:

$\mathbf {A} \mathbf {x} =\mathbf {k} ,$

(1)

حيث A هي مصفوفة n × n، وk هو متجه عمودي معطى يحتوي على n عنصر، فنقسم المصفوفة A إلى:

$\mathbf {A} =\mathbf {B} -\mathbf {C} ,$

(2)

حيث B وC هما مصفوفتان n × n. نقول أن A = B − C هي مصفوفة مقسمة، إذا كانت B⁻¹ ≥ 0 وC ≥ 0. ونفرض أن:

$\mathbf {B} \mathbf {x} =\mathbf {g} ,$

(3)

حيث g هو متجه عمودي معطى يمكن حله مباشرة للمتجه x. وبذلك يكون:

$\mathbf {B} \mathbf {x} ^{(m+1)}=\mathbf {C} \mathbf {x} ^{(m)}+\mathbf {k} ,\quad m=0,1,2,\ldots ,$

(4)

$\mathbf {x} ^{(m+1)}=\mathbf {B} ^{-1}\mathbf {C} \mathbf {x} ^{(m)}+\mathbf {B} ^{-1}\mathbf {k} ,\quad m=0,1,2,\ldots$

(5)

والمصفوفة D = B⁻¹C لا تحتوي على عناصر سالبة طالما (معادلة 2) تمثل تقسيم اعتيادي لـA.[1]

الطرق التكرارية

يمكن وصف العديد من الطرق التكرارية على أنها تقسيم للمصفوفة. إذا كانت عناصر القطر الرئيسي للمصفوفة A لا تساوي صفر، يتم التعبير عن A على أنها:

$\mathbf {A} =\mathbf {D} -\mathbf {U} -\mathbf {L} ,$

(6)

حيث D هي الجزء القطري لـA، وU وL يمثلان على الترتيب مصفوفتان مثلثيتان علوية وسفلية بأبعاد n × n، ويصبح لدينا التالي.

يمكن التعبير عن طريقة جاكوبي في صورة مصفوفة كتقسيم

$\mathbf {x} ^{(m+1)}=\mathbf {D} ^{-1}(\mathbf {U} +\mathbf {L} )\mathbf {x} ^{(m)}+\mathbf {D} ^{-1}\mathbf {k} .$ [2][3]

(7)

كذلك طريقة جاوس-سيدل

$\mathbf {x} ^{(m+1)}=(\mathbf {D} -\mathbf {L} )^{-1}\mathbf {U} \mathbf {x} ^{(m)}+(\mathbf {D} -\mathbf {L} )^{-1}\mathbf {k} .$ [3][4]

(8)

كذلك طريقة (successive over-relaxation)

$\mathbf {x} ^{(m+1)}=(\mathbf {D} -\omega \mathbf {L} )^{-1}[(1-\omega )\mathbf {D} +\omega \mathbf {U} ]\mathbf {x} ^{(m)}+\omega (\mathbf {D} -\omega \mathbf {L} )^{-1}\mathbf {k} .$ [3][5]

(9)

أمثلة

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

$\mathbf {A} ={\begin{pmatrix}6&-2&-3\\-1&4&-2\\-3&-1&5\end{pmatrix}},\quad \mathbf {k} ={\begin{pmatrix}5\\-12\\10\end{pmatrix}}$

(10)

${\begin{aligned}&\mathbf {B} ={\begin{pmatrix}6&0&0\\0&4&0\\0&0&5\end{pmatrix}},\quad \mathbf {C} ={\begin{pmatrix}0&2&3\\1&0&2\\3&1&0\end{pmatrix}},\end{aligned}}$

(11)

{\begin{aligned}&\mathbf {A^{-1}} ={\frac {1}{47}}{\begin{pmatrix}18&13&16\\11&21&15\\13&12&22\end{pmatrix}},\quad \mathbf {B^{-1}} ={\begin{pmatrix}{\frac {1}{6}}&0&0\\[4pt]0&{\frac {1}{4}}&0\\[4pt]0&0&{\frac {1}{5}}\end{pmatrix}},\end{aligned}}

{\begin{aligned}\mathbf {D} =\mathbf {B^{-1}C} ={\begin{pmatrix}0&{\frac {1}{3}}&{\frac {1}{2}}\\[4pt]{\frac {1}{4}}&0&{\frac {1}{2}}\\[4pt]{\frac {3}{5}}&{\frac {1}{5}}&0\end{pmatrix}},\quad \mathbf {B^{-1}k} ={\begin{pmatrix}{\frac {5}{6}}\\[4pt]-3\\[4pt]2\end{pmatrix}}.\end{aligned}}

وحيث أن B⁻¹ ≥ 0 وC ≥ 0، فإن التقسيم (معادلة 11) هو تقسيم اعتيادي.

$\mathbf {x} ^{(m+1)}={\begin{pmatrix}0&{\frac {1}{3}}&{\frac {1}{2}}\\[4pt]{\frac {1}{4}}&0&{\frac {1}{2}}\\[4pt]{\frac {3}{5}}&{\frac {1}{5}}&0\end{pmatrix}}\mathbf {x} ^{(m)}+{\begin{pmatrix}{\frac {5}{6}}\\[4pt]-3\\[4pt]2\end{pmatrix}},\quad m=0,1,2,\ldots$

(12)

ويصبح بذلك الحل الأمثل للمعادلة:

$\mathbf {x} ={\begin{pmatrix}2\\-1\\3\end{pmatrix}}.$

(13)

يظهر الجدول التالي أول تكرارات للمعادلة (12)، بدءً بـ $x (0) = (0.0, 0.0, 0.0) T$ . ومن الجدول نستطيع ملاحظة أن الطريقة تتجه نحو الحل (معادلة 13).

$x_{1}^{(m)}$	$x_{2}^{(m)}$	$x_{3}^{(m)}$
0.0	0.0	0.0
0.83333	-3.0000	2.0000
0.83333	-1.7917	1.9000
1.1861	-1.8417	2.1417
1.2903	-1.6326	2.3433
1.4608	-1.5058	2.4477
1.5553	-1.4110	2.5753
1.6507	-1.3235	2.6510
1.7177	-1.2618	2.7257
1.7756	-1.2077	2.7783
1.8199	-1.1670	2.8238

طريقة جاوس-سيدل

نظرًا لأن القطر الرئيسي للمصفوفة A لا يحتوي على عناصر صفرية، فيمكن كتابة المصفوفة بالصيغة التالية:

$\mathbf {D} ={\begin{pmatrix}6&0&0\\0&4&0\\0&0&5\end{pmatrix}},\quad \mathbf {U} ={\begin{pmatrix}0&2&3\\0&0&2\\0&0&0\end{pmatrix}},\quad \mathbf {L} ={\begin{pmatrix}0&0&0\\1&0&0\\3&1&0\end{pmatrix}}.$

(14)

{\begin{aligned}&\mathbf {(D-L)^{-1}} ={\frac {1}{120}}{\begin{pmatrix}20&0&0\\5&30&0\\13&6&24\end{pmatrix}},\end{aligned}}

{\begin{aligned}&\mathbf {(D-L)^{-1}U} ={\frac {1}{120}}{\begin{pmatrix}0&40&60\\0&10&75\\0&26&51\end{pmatrix}},\quad \mathbf {(D-L)^{-1}k} ={\frac {1}{120}}{\begin{pmatrix}100\\-335\\233\end{pmatrix}}.\end{aligned}}

$\mathbf {x} ^{(m+1)}={\frac {1}{120}}{\begin{pmatrix}0&40&60\\0&10&75\\0&26&51\end{pmatrix}}\mathbf {x} ^{(m)}+{\frac {1}{120}}{\begin{pmatrix}100\\-335\\233\end{pmatrix}},\quad m=0,1,2,\ldots$

(15)

يظهر الجدول التالي أول تكرارات للمعادلة (15)، بدءً بـ $x (0) = (0.0, 0.0, 0.0) T$ . ومن الجدول نستطيع ملاحظة أن الطريقة تتجه نحول الحل (معادلة 13)، لكن أسرع من طريقة جاكوبي الموضحة أعلاه.

$x_{1}^{(m)}$	$x_{2}^{(m)}$	$x_{3}^{(m)}$
0.0	0.0	0.0
0.8333	-2.7917	1.9417
0.8736	-1.8107	2.1620
1.3108	-1.5913	2.4682
1.5370	-1.3817	2.6459
1.6957	-1.2531	2.7668
1.7990	-1.1668	2.8461
1.8675	-1.1101	2.8985
1.9126	-1.0726	2.9330
1.9423	-1.0479	2.9558
1.9619	-1.0316	2.9708

ملاحظات

Varga (1960, pp. 121–122)
Burden & Faires (1993, p. 408)
Varga (1962, p. 88)
Burden & Faires (1993, p. 411)
Burden & Faires (1993, p. 416)

مراجع

Burden, Richard L.; Faires, J. Douglas (1993), Numerical Analysis (الطبعة 5th), Boston: Prindle, Weber and Schmidt, ISBN 0-534-93219-3 الوسيط |CitationClass= تم تجاهله (مساعدة); الوسيط |separator= تم تجاهله (مساعدة)CS1 maint: ref=harv (link).
Varga, Richard S. (1960). "Factorization and Normalized Iterative Methods". In Langer, Rudolph E. (المحرر). Boundary Problems in Differential Equations. Madison: University of Wisconsin Press. صفحات 121–142. LCCN 60-60003. الوسيط |CitationClass= تم تجاهله (مساعدة)
Varga, Richard S. (1962), Matrix Iterative Analysis, New Jersey: Prentice-Hall, LCCN 62-21277 الوسيط |CitationClass= تم تجاهله (مساعدة); الوسيط |separator= تم تجاهله (مساعدة)CS1 maint: ref=harv (link).

بوابة رياضيات

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.

[1] Varga (1960, pp. 121–122)

[2] Burden & Faires (1993, p. 408)

[Varga_1962_88-3] Varga (1962, p. 88)

[4] Burden & Faires (1993, p. 411)

[5] Burden & Faires (1993, p. 416)