مصفوفة توبليتز

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

[أبجدهـوأبجدزوأبجحزوأبأناحزوأ].{\displaystyle \qquad {\begin{bmatrix}a&b&c&d&e\\f&a&b&c&d\\g&f&a&b&c\\h&g&f&a&b\\i&h&g&f&a\end{bmatrix}}.}

أين×ن{\displaystyle n\times n}مصفوفةأ{\displaystyle A}من الشكل

أ=[أ0أ-1أ-2أ-(ن-1)أ1أ0أ-1أ2أ1أ-1أ-2أ1أ0أ-1أن-1أ2أ1أ0]{\displaystyle A={\begin{bmatrix}a_{0}&a_{-1}&a_{-2}&\cdots &\cdots &a_{-(n-1)}\\a_{1}&a_{0}&a_{-1}&\ddots &&\vdots \\a_{2}&a_{1}&\ddots &\ddots &\ddots &\vdots \\\vdots &\ddots &\ddots &\ddots &a_{-1}&a_{-2}\\\vdots &&\ddots &a_{1}&a_{0}&a_{-1}\\a_{n-1}&\cdots &\cdots &a_{2}&a_{1}&a_{0}\end{bmatrix}}}

هي مصفوفة توبليتز . إذاأنا،ج{\displaystyle i,j}عنصر منأ{\displaystyle A}يُشار إليه بـأأنا،ج{\displaystyle A_{i,j}}ثم لدينا

أأنا،ج=أأنا+1،ج+1=أأنا-ج.{\displaystyle A_{i,j}=A_{i+1,j+1}=a_{ij}.}

مصفوفة توبليتز ليست بالضرورة مربعة .

حل نظام توبليتز

معادلة مصفوفية على الشكل

أx=ب{\displaystyle Ax=b}

يُطلق عليه نظام توبليتز إذاأ{\displaystyle A}هي مصفوفة توبليتز. إذاأ{\displaystyle A}هون×ن{\displaystyle n\times n}إذا كانت مصفوفة توبليتز، فإن النظام يحتوي على الأكثر على 1000 مصفوفة توبليتز فقط. 2ن-1{\displaystyle 2n-1}قيم فريدة، بدلاً منن2{\displaystyle n^{2}}لذلك قد نتوقع أن يكون حل نظام توبليتز أسهل، وهذا هو الحال بالفعل.

يمكن حل أنظمة توبليتز بواسطة خوارزميات مثل خوارزمية شور أو خوارزمية ليفينسون فييا(ن2){\displaystyle O(n^{2})}[ 1 ] [ 2 ] وقد ثبت أن متغيرات هذه الطريقة الأخيرة مستقرة بشكل ضعيف (أي أنها تُظهر استقرارًا عدديًا للأنظمة الخطية جيدة التكييف ). [ 3 ] ويمكن أيضًا استخدام الخوارزميات لإيجاد محدد مصفوفة توبليتز فييا(ن2){\displaystyle O(n^{2})}الوقت. [ 4 ]

يمكن أيضًا تحليل مصفوفة توبليتز (أي تحليلها إلى عوامل)يا(ن2){\displaystyle O(n^{2})}[ 5 ] خوارزمية Bareiss لتحليل LU مستقرة. [ 6 ] يوفر تحليل LU طريقة سريعة لحل نظام Toeplitz، وكذلك لحساب المحدد. باستخدام رتبة الإزاحة، نحصل على طريقة تتطلبيا~(αω-1ن){\displaystyle {\tilde {O}}({\alpha ^{\omega -1}}n)}العمليات باستخدام خوارزميات ضرب المصفوفات السريعة ، حيثα{\displaystyle \alpha }هي الرتبة و2.37ω<3{\displaystyle ^{\sim}2.37\leq \omega <3}[ 7 ]

ملكيات

  • أنن×ن{\displaystyle n\times n}يمكن تعريف مصفوفة توبليتز على أنها مصفوفةأ{\displaystyle A}أينأأنا،ج=جأنا-ج{\displaystyle A_{i,j}=c_{ij}}، بالنسبة للثوابتج1-ن،...،جن-1{\displaystyle c_{1-n},\ldots ,c_{n-1}}مجموعةن×ن{\displaystyle n\times n}مصفوفات توبليتز هي فضاء جزئي من الفضاء المتجهي لـن×ن{\displaystyle n\times n}المصفوفات (في ظل جمع المصفوفات والضرب القياسي).
  • يمكن إضافة مصفوفتين من نوع توبليتز فييا(ن){\displaystyle O(n)}الوقت (عن طريق تخزين قيمة واحدة فقط لكل قطر) وضربه فييا(ن2){\displaystyle O(n^{2})}وقت.
  • مصفوفات توبليتز متناظرة حول المحور . مصفوفات توبليتز المتناظرة متناظرة مركزياً وثنائية التناظر .
  • ترتبط مصفوفات توبليتز ارتباطًا وثيقًا بمتسلسلات فورييه ، إذ يمكن تمثيل عملية الضرب بمتعددة حدود مثلثية ، بعد ضغطها إلى فضاء محدود الأبعاد، بواسطة هذه المصفوفة. وبالمثل، يمكن تمثيل الالتفاف الخطي كضرب بمصفوفة توبليتز.
  • تُعتبر مصفوفات توبليتز مكافئة تقاربياً للمصفوفات الدائرية مع ازدياد الأبعاد، وهي نتيجة تُعرف باسم نظرية غريناندر-سيغو . [ 8 ] هذه الخاصية الدائرية التقاربية هي السبب في أن تحويل فورييه المنفصل يُقطّر مصفوفات توبليتز الكبيرة تقريبًا، وهي أساس فعالية تقدير الكثافة الطيفية القائم على تحويل فورييه المنفصل للعمليات المستقرة .
  • تتبادل مصفوفات توبليتز تقاربياً . وهذا يعني أنها تقطر في نفس الأساس عندما يؤول بُعد الصف والعمود إلى اللانهاية.
  • بالنسبة لمصفوفات توبليتز المتناظرة، يوجد التفكيك
1أ0أ=جيجيتي-(جي-أنا)(جي-أنا)تي{\displaystyle {\frac {1}{a_{0}}}A=GG^{\operatorname {T} }-(GI)(GI)^{\operatorname {T} }}
أينجي{\displaystyle G}هو الجزء المثلث السفلي من1أ0أ{\displaystyle {\frac {1}{a_{0}}}A}.
أ-1=1α0(ببتي-ججتي){\displaystyle A^{-1}={\frac {1}{\alpha _{0}}}(BB^{\operatorname {T} }-CC^{\operatorname {T} })}
أينب{\displaystyle B}وج{\displaystyle C}هي مصفوفات توبليتز المثلثية السفلية وج{\displaystyle C}هي مصفوفة مثلثية سفلية تمامًا. [ 9 ]

الالتفاف المنفصل

يمكن بناء عملية الالتفاف كضرب مصفوفات، حيث يتم تحويل أحد المدخلات إلى مصفوفة توبليتز. على سبيل المثال، التفافح{\displaystyle h}وx{\displaystyle x}يمكن صياغتها على النحو التالي:

y=ح*x=[ح1000ح2ح1ح3ح200ح3ح10حم-1ح2ح1حمحم-1ح20حمحم-200حم-1حم-2حمحم-1000حم][x1x2x3xن]{\displaystyle y=h\ast x={\begin{bmatrix}h_{1}&0&\cdots &0&0\\h_{2}&h_{1}&&\vdots &\vdots \\h_{3}&h_{2}&\cdots &0&0\\\vdots &h_{3}&\cdots &h_{1}&0\\h_{m-1}&\vdots &\ddots &h_{2}&h_{1}\\h_{m}&h_{m-1}&&\vdots &h_{2}\\0&h_{m}&\ddots &h_{m-2}&\vdots \\0&0&\cdots &h_{m-1}&h_{m-2}\\\vdots &\vdots &&h_{m}&h_{m-1}\\0&0&0&\cdots &h_{m}\end{bmatrix}}{\begin{bmatrix}x_{1}\\x_{2}\\x_{3}\\\vdots \\x_{n}\end{bmatrix}}}
yتي=[ح1ح2ح3حم-1حم][x1x2x3xن00000x1x2x3xن00000x1x2x3...xن00000x1xن-2xن-1xن00000x1xن-2xن-1xن].{\displaystyle y^{T}={\begin{bmatrix}h_{1}&h_{2}&h_{3}&\cdots &h_{m-1}&h_{m}\end{bmatrix}}{\begin{bmatrix}x_{1}&x_{2}&x_{3}&\cdots &x_{n}&0&0&0&\cdots &0\\0&x_{1}&x_{2}&x_{3}&\cdots &x_{n}&0&0&\cdots &0\\0&0&x_{1}&x_{2}&x_{3}&\ldots &x_{n}&0&\cdots &0\\\vdots &&\vdots &\vdots &\vdots &&\vdots &\vdots &&\vdots \\0&\cdots &0&0&x_{1}&\cdots &x_{n-2}&x_{n-1}&x_{n}&0\\0&\cdots &0&0&0&x_{1}&\cdots &x_{n-2}&x_{n-1}&x_{n}\end{bmatrix}}.}

يمكن توسيع هذا النهج لحساب الارتباط الذاتي ، والارتباط المتبادل ، والمتوسط ​​المتحرك، وما إلى ذلك.

مصفوفة توبليتز اللانهائية

مصفوفة توبليتز ثنائية اللانهاية (أي المدخلات مفهرسة بواسطةZ×Z{\displaystyle \mathbb {Z} \times \mathbb {Z} })أ{\displaystyle A}يُحدث مؤثرًا خطيًا على2{\displaystyle \ell ^{2}}.

أ=[أ0أ-1أ-2أ-3أ1أ0أ-1أ-2أ2أ1أ0أ-1أ3أ2أ1أ0].{\displaystyle A={\begin{bmatrix}&\vdots &\vdots &\vdots &\vdots \\\cdots &a_{0}&a_{-1}&a_{-2}&a_{-3}&\cdots \\\cdots &a_{1}&a_{0}&a_{-1}&a_{-2}&\cdots \\\cdots &a_{2}&a_{1}&a_{0}&a_{-1}&\cdots \\\cdots &a_{3}&a_{2}&a_{1}&a_{0}&\cdots \\&\vdots &\vdots &\vdots &\vdots \end{bmatrix}}.}

يكون المؤثر المستحث محدودًا إذا وفقط إذا كانت معاملات مصفوفة توبليتزأ{\displaystyle A}هي معاملات فورييه لدالة محدودة أساسًاو{\displaystyle f}.

في مثل هذه الحالات،و{\displaystyle f}يُطلق عليه رمز مصفوفة توبليتزأ{\displaystyle A}والمعيار الطيفي لمصفوفة توبليتزأ{\displaystyle A}يتزامن معل{\displaystyle L^{\infty }}معيار رمزه. يمكن إيجاد البرهان في النظرية 1.1 من بوتشر وغرودسكي. [ 10 ]

انظر أيضاً

ملحوظات

  1. بريس وآخرون، 2007 ، §2.8.2 مصفوفات توبليتز
  2. هايز 1996 ، الفصل 5.2.6
  3. كريشنا ووانغ 1993
  4. موناهان 2011 ، §4.5 أنظمة توبليتز
  5. برنت 1999
  6. بويانكزيك وآخرون 1995
  7. بوستان، أ.؛ جينرود، س.-ب.؛ شوست، إ. (2008). "حل الأنظمة الخطية المهيكلة ذات رتبة الإزاحة الكبيرة". علوم الحاسوب النظرية . 407 ( 1-3 ): 155-181 . doi : 10.1016/j.tcs.2008.05.014 .
  8. غريناندر، أولف؛ سيغو، غابور (1958). أشكال توبليتز وتطبيقاتها . بيركلي، كاليفورنيا: مطبعة جامعة كاليفورنيا.
  9. موخرجي ومايتي 1988
  10. بوتشر وغرودسكي 2012

مراجع

للمزيد من القراءة