مصفوفة ثلاثية الأقطار

في الجبر الخطي ، المصفوفة ثلاثية الأقطار هي مصفوفة شريطية تحتوي على عناصر غير صفرية فقط على القطر الرئيسي ، والقطر الفرعي/القطر السفلي (أول قطر أسفله)، والقطر العلوي/القطر العلوي (أول قطر أعلى القطر الرئيسي). على سبيل المثال، المصفوفة التالية ثلاثية الأقطار :

(1400341002340013).{\displaystyle {\begin{pmatrix}1&4&0&0\\3&4&1&0\\0&2&3&4\\0&0&1&3\\\end{pmatrix}}.}

محدد المصفوفة ثلاثية الأقطار يُعطى بواسطة امتداد عناصرها. [ 1 ]

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

ملكيات

المصفوفة ثلاثية الأقطار هي مصفوفة هيسنبرغ علوية وسفلية في آن واحد . [ 2 ] وبالتحديد، المصفوفة ثلاثية الأقطار هي مجموع مباشر لـ p مصفوفة من الرتبة 1×1 و q مصفوفة من الرتبة 2×2، بحيث يكون p + q / 2 = n ، وهو بُعد المصفوفة ثلاثية الأقطار. مع أن المصفوفة ثلاثية الأقطار العامة ليست بالضرورة متناظرة أو هيرميتية ، فإن العديد من المصفوفات التي تظهر عند حل مسائل الجبر الخطي تتمتع بإحدى هاتين الخاصيتين. علاوة على ذلك، إذا كانت المصفوفة ثلاثية الأقطار الحقيقية A تحقق الشرط a <sub> k , k +1</sub> حيث a <sub> k +1, k </sub> > 0 لجميع قيم k ، بحيث تكون إشارات عناصرها متناظرة، فإنها تُشابه المصفوفة الهيرميتية، وذلك بتغيير قطري لمصفوفة الأساس. وبالتالي، فإن قيمها الذاتية حقيقية. إذا استبدلنا المتباينة الصارمة بـ a k , k +1 a k +1, k 0، فبسبب الاستمرارية، تظل القيم الذاتية مضمونة أن تكون حقيقية، ولكن المصفوفة لم تعد بالضرورة مشابهة لمصفوفة هيرميتية. [ 3 ]

تشكل مجموعة جميع المصفوفات ثلاثية الأقطار من الرتبة n × n فضاءً متجهيًا ذا أبعاد 3n -2 .

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

المحدد

يمكن حساب محدد المصفوفة ثلاثية الأقطار A من الرتبة n من علاقة تكرارية ثلاثية الحدود . [ 4 ] اكتب f1 = | a1 | = a1 ( أي أن f1 هو محدد المصفوفة 1×1 التي تتكون فقط من a1 ) ، ولتكن    

ون=|أ1ب1ج1أ2ب2ج2بن-1جن-1أن|.{\displaystyle f_{n}={\begin{vmatrix}a_{1}&b_{1}\\c_{1}&a_{2}&b_{2}\\&c_{2}&\ddots &\ddots \\&&\ddots &\ddots &b_{n-1}\\&&&c_{n-1}&a_{n}\end{vmatrix}}.}

تُسمى المتتالية ( f i ) بالمتتالية المستمرة ، وهي تحقق علاقة التكرار.

ون=أنون-1-جن-1بن-1ون-2{\displaystyle f_{n}=a_{n}f_{n-1}-c_{n-1}b_{n-1}f_{n-2}}

بقيم ابتدائية f 0  =  1 و f −1  =  0. تكلفة حساب محدد المصفوفة ثلاثية الأقطار باستخدام هذه الصيغة خطية في n ، بينما التكلفة تكعيبية بالنسبة للمصفوفة العامة.

الانقلاب

معكوس مصفوفة ثلاثية الأقطار غير منفردة T

تي=(أ1ب1ج1أ2ب2ج2بن-1جن-1أن){\displaystyle T={\begin{pmatrix}a_{1}&b_{1}\\c_{1}&a_{2}&b_{2}\\&c_{2}&\ddots &\ddots \\&&\ddots &\ddots &b_{n-1}\\&&&c_{n-1}&a_{n}\end{pmatrix}}}

يُعطى بواسطة

(تي-1)أناج={(-1)أنا+جبأنابج-1θأنا-1ϕج+1/θن لو أنا<جθأنا-1ϕج+1/θن لو أنا=ج(-1)أنا+ججججأنا-1θج-1ϕأنا+1/θن لو أنا>ج{\displaystyle (T^{-1})_{ij}={\begin{cases}(-1)^{i+j}b_{i}\cdots b_{j-1}\theta _{i-1}\phi _{j+1}/\theta _{n}&{\text{ if }}i<j\\\theta _{i-1}\phi _{j+1}/\theta _{n}&{\text{ if }}i=j\\(-1)^{i+j}c_{j}\cdots c_{i-1}\theta _{j-1}\phi _{i+1}/\theta _{n}&{\text{ if }}i>j\\\end{cases}}}

حيث تحقق θ i علاقة التكرار

θأنا=أأناθأنا-1-بأنا-1جأنا-1θأنا-2أنا=2،3،...،ن{\displaystyle \theta _{i}=a_{i}\theta _{i-1}-b_{i-1}c_{i-1}\theta _{i-2}\qquad i=2,3,\ldots ,n}

مع الشروط الابتدائية θ₀ = 1 ، θ₁ = a₁ ، و ϕᵢ تحقق    

ϕأنا=أأناϕأنا+1-بأناجأناϕأنا+2أنا=ن-1،...،1{\displaystyle \phi _{i}=a_{i}\phi _{i+1}-b_{i}c_{i}\phi _{i+2}\qquad i=n-1,\ldots ,1}

مع الشروط الأولية ϕ n +1  =  1 و ϕ n  = a n . [ 5 ] [ 6 ] 

يمكن حساب الحلول المغلقة لحالات خاصة مثل المصفوفات المتناظرة التي تتساوى فيها جميع العناصر القطرية وغير القطرية [ 7 ] أو مصفوفات توبليتز [ 8 ] ، وكذلك للحالة العامة. [ 9 ] [ 10 ]

بشكل عام، معكوس المصفوفة ثلاثية الأقطار هو مصفوفة شبه قابلة للفصل ، والعكس صحيح. [ 11 ] يمكن كتابة معكوس المصفوفة ثلاثية الأقطار المتناظرة كمصفوفة أحادية الزوج (تُعرف أيضًا باسم المصفوفة شبه القابلة للفصل والممثلة بالمولد ) على الصورة [ 12 ] [ 13 ].

(α1-β1-β1α2-β2-βن-1-βن-1αن)-1=(أ1ب1أ1ب2أ1بنأ1ب2أ2ب2أ2بنأ1بنأ2بنأنبن)=(أمين(أنا،ج)بالأعلى(أنا،ج)){\displaystyle {\begin{pmatrix}\alpha _{1}&-\beta _{1}\\-\beta _{1}&\alpha _{2}&-\beta _{2}\\&\ddots &\ddots &\ddots &\\&&\ddots &\ddots &-\beta _{n-1}\\&&&-\beta _{n-1}&\alpha _{n}\end{pmatrix}}^{-1}={\begin{pmatrix}a_{1}b_{1}&a_{1}b_{2}&\cdots &a_{1}b_{n}\\a_{1}b_{2}&a_{2}b_{2}&\cdots &a_{2}b_{n}\\\vdots &\vdots &\ddots &\vdots \\a_{1}b_{n}&a_{2}b_{n}&\cdots &a_{n}b_{n}\end{pmatrix}}=\left(a_{\min(i,j)}b_{\max(i,j)}\right)}

أين{أأنا=βأناβن-1دلتاأنادلتانبنبأنا=β1βأنا-1د1دأنا{\displaystyle {\begin{cases}\displaystyle a_{i}={\frac {\beta _{i}\cdots \beta _{n-1}}{\delta _{i}\cdots \delta _{n}\,b_{n}}}\\\displaystyle b_{i}={\frac {\beta _{1}\cdots \beta _{i-1}}{d_{1}\cdots d_{i}}}\end{cases}}} مع{دن=αن،دأنا-1=αأنا-1-βأنا-12دأنا،أنا=ن،ن-1،،2،دلتا1=α1،دلتاأنا+1=αأنا+1-βأنا2دلتاأنا،أنا=1،2،،ن-1.{\displaystyle {\begin{cases}d_{n}=\alpha _{n},\quad d_{i-1}=\alpha _{i-1}-{\frac {\beta _{i-1}^{2}}{d_{i}}},&i=n,n-1,\cdots ,2,\\\delta _{1}=\alpha _{1},\quad \delta _{i+1}=\alpha _{i+1}-{\frac {\beta _{i}^{2}}{\delta _{i}}},&i=1,2,\cdots ,n-1.\end{cases}}}

حل النظام الخطي

نظام المعادلات Ax  = b لـ  بRن{\displaystyle b\in \mathbb {R} ^{n}}يمكن حلها بواسطة شكل فعال من أشكال الحذف الغاوسي عندما تكون المصفوفة A ثلاثية الأقطار تسمى خوارزمية المصفوفة ثلاثية الأقطار ، وتتطلب O ( n ) عملية. [ 14 ]

القيم الذاتية

عندما تكون المصفوفة ثلاثية الأقطار أيضًا من نوع Toeplitz ، يوجد حل بسيط مغلق الشكل لقيمها الذاتية، وهو: [ 15 ] [ 16 ]

أ-2بجكوس(كπن+1)،ك=1،...،ن.{\displaystyle a-2{\sqrt {bc}}\cos \left({\frac {k\pi }{n+1}}\right),\qquad k=1,\ldots ,n.}

المصفوفة الثلاثية الأقطار المتناظرة الحقيقية لها قيم ذاتية حقيقية، وتكون جميع القيم الذاتية متميزة (بسيطة) إذا كانت جميع العناصر غير القطرية غير صفرية. [ 17 ] توجد طرق عديدة للحساب العددي للقيم الذاتية لمصفوفة ثلاثية الأقطار متناظرة حقيقية بدقة محدودة اختيارية، وعادةً ما تتطلبيا(ن2){\displaystyle O(n^{2})}عمليات لمصفوفة بحجمن×ن{\displaystyle n\times n}على الرغم من وجود خوارزميات سريعة لا تتطلب سوى (بدون حساب متوازٍ)يا(نسجلن){\displaystyle O(n\log n)}[ 18 ]

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

التشابه مع المصفوفة ثلاثية الأقطار المتناظرة

بالنسبة للمصفوفات ثلاثية الأقطار غير المتناظرة ، يمكن حساب تحليل القيم الذاتية باستخدام تحويل التشابه. بالنظر إلى مصفوفة ثلاثية الأقطار حقيقية غير متناظرة

تي=(أ1ب1ج1أ2ب2ج2بن-1جن-1أن){\displaystyle T={\begin{pmatrix}a_{1}&b_{1}\\c_{1}&a_{2}&b_{2}\\&c_{2}&\ddots &\ddots \\&&\ddots &\ddots &b_{n-1}\\&&&c_{n-1}&a_{n}\end{pmatrix}}}

أينبأناجأنا{\displaystyle b_{i}\neq c_{i}}افترض أن كل حاصل ضرب للعناصر غير القطرية موجب تمامًابأناجأنا>0{\displaystyle b_{i}c_{i}>0}وحدد مصفوفة التحويلد{\displaystyle D}بواسطة [ 20 ]

د:=التشخيص(دلتا1،...،دلتان)لدلتاأنا:={1،أنا=1جأنا-1...ج1بأنا-1...ب1،أنا=2،...،ن.{\displaystyle D:=\operatorname {diag} (\delta _{1},\dots ,\delta _{n})\quad {\text{for}}\quad \delta _{i}:={\begin{cases}1&,\,i=1\\{\sqrt {\frac {c_{i-1}\dots c_{1}}{b_{i-1}\dots b_{1}}}}&,\,i=2,\dots ,n\,.\end{cases}}}

تحويل التشابهد-1تيد{\displaystyle D^{-1}TD}ينتج عنه مصفوفة ثلاثية الأقطار متناظرةج{\displaystyle J}بواسطة: [ 21 ] [ 20 ]

ج:=د-1تيد=(أ1علامةب1ب1ج1علامةب1ب1ج1أ2علامةب2ب2ج2علامةب2ب2ج2علامةبن-1بن-1جن-1علامةبن-1بن-1جن-1أن).{\displaystyle J:=D^{-1}TD={\begin{pmatrix}a_{1}&\operatorname {sgn} b_{1}\,{\sqrt {b_{1}c_{1}}}\\\operatorname {sgn} b_{1}\,{\sqrt {b_{1}c_{1}}}&a_{2}&\operatorname {sgn} b_{2}\,{\sqrt {b_{2}c_{2}}}\\&\operatorname {sgn} b_{2}\,{\sqrt {b_{2}c_{2}}}&\ddots &\ddots \\&&\ddots &\ddots &\operatorname {sgn} b_{n-1}\,{\sqrt {b_{n-1}c_{n-1}}}\\&&&\operatorname {sgn} b_{n-1}\,{\sqrt {b_{n-1}c_{n-1}}}&a_{n}\end{pmatrix}}\,.}

لاحظ أنتي{\displaystyle T}وج{\displaystyle J}لها نفس القيم الذاتية.

برمجة الحاسوب

التحويل الذي يُختزل مصفوفة عامة إلى شكل هيسنبرغ يُختزل مصفوفة هيرميتية إلى شكل ثلاثي الأقطار. لذا، فإن العديد من خوارزميات القيم الذاتية ، عند تطبيقها على مصفوفة هيرميتية، تُختزل مصفوفة هيرميت المُدخلة إلى شكل ثلاثي الأقطار (حقيقي متناظر) كخطوة أولى. [ 22 ]

يمكن تخزين المصفوفة ثلاثية الأقطار بكفاءة أكبر من المصفوفة العامة باستخدام نظام تخزين خاص . على سبيل المثال، تقوم حزمة LAPACK Fortran بتخزين مصفوفة ثلاثية الأقطار غير متناظرة من الرتبة n في ثلاث مصفوفات أحادية البعد، إحداها بطول n تحتوي على العناصر القطرية، واثنتان بطول n 1 تحتويان على العناصر تحت القطرية والعناصر فوق القطرية .

التطبيقات

التجزئة المكانية لمعادلة الانتشار أو الحرارة أحادية البعد

u(ت،x)ت=α2u(ت،x)x2{\displaystyle {\frac {\partial u(t,x)}{\partial t}}=\alpha {\frac {\partial ^{2}u(t,x)}{\partial x^{2}}}}

ينتج عن استخدام الفروق المركزية المحدودة من الرتبة الثانية ما يلي:

(u1(ت)تu2(ت)تuشمال(ت)ت)=αΔx2(-210...01-21001-210...01-2)(u1(ت)u2(ت)uشمال(ت)){\displaystyle {\begin{pmatrix}{\frac {\partial u_{1}(t)}{\partial t}}\\{\frac {\partial u_{2}(t)}{\partial t}}\\\vdots \\{\frac {\partial u_{N}(t)}{\partial t}}\end{pmatrix}}={\frac {\alpha }{\Delta x^{2}}}{\begin{pmatrix}-2&1&0&\ldots &0\\1&-2&1&\ddots &\vdots \\0&\ddots &\ddots &\ddots &0\\\vdots &&1&-2&1\\0&\ldots &0&1&-2\end{pmatrix}}{\begin{pmatrix}u_{1}(t)\\u_{2}(t)\\\vdots \\u_{N}(t)\\\end{pmatrix}}}

مع ثابت التقطيع Δx{\displaystyle \Delta x}المصفوفة ثلاثية الأقطارأأنا=-2{\displaystyle a_{i}=-2}وبأنا=جأنا=1{\displaystyle b_{i}=c_{i}=1}ملاحظة: لم يتم تحديد أي شروط حدودية بشكل صريح، ولكن هذه المصفوفة تتوافق مع شروط حدود نيومان (تدرج صفري).

انظر أيضاً

ملحوظات

  1. توماس موير (1960). رسالة في نظرية المحددات . منشورات دوفر . الصفحات 516-525 . 
  2. هورن، روجر أ.؛ جونسون، تشارلز ر. (1985). تحليل المصفوفات . مطبعة جامعة كامبريدج. ص 28. ISBN  0521386322.
  3. هورن وجونسون، صفحة 174
  4. الميكاوي، MEA (2004). "حول معكوس المصفوفة ثلاثية الأقطار العامة". الرياضيات التطبيقية والحساب . 150 (3): 669-679 . doi : 10.1016/S0096-3003(03)00298-4 .
  5. دا فونسيكا، سي إم (2007). "حول القيم الذاتية لبعض المصفوفات ثلاثية الأقطار" . مجلة الرياضيات الحسابية والتطبيقية . 200 : 283-286 . doi : 10.1016/j.cam.2005.08.047 .
  6. عثماني، ر. أ. (1994). "عكس مصفوفة جاكوبي ثلاثية الأقطار" . الجبر الخطي وتطبيقاته . 212-213 : 413-414 . doi : 10.1016/0024-3795(94)90414-6 .
  7. هو، جي واي؛ أوكونيل، آر إف (1996). "التحليل العكسي للمصفوفات ثلاثية الأقطار المتناظرة". مجلة الفيزياء أ: الرياضية والعامة . 29 (7): 1511. رمز Bibcode : 1996JPhA...29.1511H . doi : 10.1088/0305-4470/29/7/020 .
  8. هوانغ، واي؛ ماكول، دبليو إف (1997). "التحليل العكسي للمصفوفات ثلاثية الأقطار العامة". مجلة الفيزياء أ: الرياضية والعامة . 30 (22): 7919. Bibcode : 1997JPhA...30.7919H . doi : 10.1088/0305-4470/30/22/026 .
  9. مالك، ر.ك. (2001). "معكوس المصفوفة ثلاثية الأقطار" . الجبر الخطي وتطبيقاته . 325 ( 1-3 ): 109-139 . doi : 10.1016/S0024-3795(00)00262-7 .
  10. كيليتش، إي. (2008). "صيغة صريحة لإيجاد معكوس مصفوفة ثلاثية الأقطار باستخدام الكسور المستمرة العكسية". الرياضيات التطبيقية والحساب . 197 : 345-357 . doi : 10.1016/j.amc.2007.07.046 .
  11. راف فاندبريل؛ مارك فان باريل؛ نيكولا ماستروناردي (2008). حسابات المصفوفات والمصفوفات شبه المنفصلة. المجلد الأول: الأنظمة الخطية . مطبعة جامعة جونز هوبكنز. النظرية 1.38، ص 41. ISBN 978-0-8018-8714-7.
  12. ميوران، جيرارد (1992). "مراجعة حول معكوس المصفوفات المتناظرة ثلاثية الأقطار والمصفوفات ثلاثية الأقطار الكتلية" . مجلة SIAM لتحليل المصفوفات وتطبيقاتها . 13 (3): 707-728 . doi : 10.1137/0613045 .
  13. بوسو، سيباستيان (2024). "المصفوفات ثلاثية الأقطار والمصفوفات أحادية الزوج والمجموع العكسي لمصفوفتين أحاديتي الزوج" . الجبر الخطي وتطبيقاته . 699 : 129-158 . arXiv : 2304.06100 . doi : 10.1016/j.laa.2024.06.018 .
  14. غولوب، جين هـفان لون، تشارلز ف. (1996). حسابات المصفوفات ( الطبعة الثالثة). مطبعة جامعة جونز هوبكنز. ISBN  0-8018-5414-8.
  15. نوشيز، س.؛ باسكويني، ل.؛ رايشل، ل. (2013). "مصفوفات توبليتز ثلاثية الأقطار: خصائصها وتطبيقاتها الجديدة". الجبر الخطي العددي مع التطبيقات . 20 (2): 302. doi : 10.1002/nla.1811 .
  16. يمكن كتابة هذا أيضًا على النحو التاليأ+2بجكوس(كπ/(ن+1)){\displaystyle a+2{\sqrt {bc}}\cos(k\pi /{(n+1)})}لأنكوس(x)=-كوس(π-x){\displaystyle \cos(x)=-\cos(\pi -x)}كما هو موضح في: كولكارني، د.؛ شميدت، د.؛ تسوي، س.ك. (1999). "القيم الذاتية لمصفوفات شبه توبليتز ثلاثية الأقطار" (ملف PDF) . الجبر الخطي وتطبيقاته . 297 ( 1-3 ): 63-80 . doi : 10.1016/S0024-3795(99)00114-7 .
  17. بارليت، بي إن (1997) [1980]. مسألة القيم الذاتية المتناظرة . كلاسيكيات في الرياضيات التطبيقية. المجلد 20. سيام. ISBN  0-89871-402-8. OCLC 228147822 . 
  18. كوكلي، إي إس؛ روخلين، في. (2012). "خوارزمية سريعة للتجزئة والتغلب لحساب أطياف المصفوفات الثلاثية القطرية المتناظرة الحقيقية" . التحليل التوافقي التطبيقي والحسابي . 34 (3): 379-414 . doi : 10.1016/j.acha.2012.06.003 .
  19. ديلون، إندرجيت سينغ (1997). خوارزمية جديدة من رتبة O(n² ) لمسألة القيم الذاتية/المتجهات الذاتية المتناظرة ثلاثية الأقطار (ملف PDF) (أطروحة دكتوراه). جامعة كاليفورنيا، بيركلي. ص 8. CSD-97-971، ADA637073. 
  20. 1 2 كريير، م. (1994). "عمليات الولادة والوفاة التحليلية: منهج فضاء هيلبرت". العمليات العشوائية وتطبيقاتها . 49 (1): 65-74 . doi : 10.1016/0304-4149(94)90112-0 .
  21. ميوران، جيرار (أكتوبر 2008). "المصفوفات ثلاثية الأقطار" (PDF) - عبر معهد الرياضيات الحاسوبية، جامعة هونغ كونغ المعمدانية.
  22. إيدلمان، يولي؛ جوهبرغ، إسرائيل؛ جيمينياني، لوكا (1 يناير 2007). "حول الاختزال السريع لمصفوفة شبه قابلة للفصل إلى صيغ هيسنبرغ وثلاثية الأقطار" . الجبر الخطي وتطبيقاته . 420 (1): 86-101 . doi : 10.1016/j.laa.2006.06.028 . ISSN 0024-3795 .