متعدد الأطياف

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

تعريف

في منشوراته الأصلية، يعتبر بريلينجر عملية ثابتة (الإشارة).z(ت){\displaystyle z(t)}ومعاملها التراكمي متعدد الأوقات . [ 1 ]و(τ1،،τن-1)=جن(z(ت)،z(ت+τ1)،،z(ت+τن-1)).{\displaystyle f(\tau _{1},\cdots ,\tau _{n-1})=C_{n}(z(t),z(t+\tau _{1}),\cdots ,z(t+\tau _{n-1})).} لا يعتمد التراكم علىت{\displaystyle t}منذz(ت){\displaystyle z(t)}ثابت. وهذا يسمح بتعريفن{\displaystyle n}طيف متعدد الأطياف من الرتبة -th Sz(ن){\displaystyle S_{z}^{(n)}}إشارة ثابتةz(ت){\displaystyle z(t)}بواسطة جن(z(ω1)،،z(ωن))=2πدلتا(ω1++ωن)Sz(ن)(ω1،،ωن-1){\displaystyle C_{n}{\big (}z(\omega _{1}),\cdots ,z(\omega _{n}){\big )}=2\pi \delta (\omega _{1}+\cdots +\omega _{n})S_{z}^{(n)}(\omega _{1},\cdots ,\omega _{n-1})} باستخدام تحويل فورييهz(ω)=-+هـأناωτz(τ)دτ،{\displaystyle z(\omega )=\int _{-\infty }^{+\infty }e^{i\omega \tau }z(\tau )d\tau ,} ودالة ديراك دلتادلتا{\displaystyle \delta }والذي يبدو أنه ناتج عن استقلالو{\displaystyle f}علىت{\displaystyle t}[ 2 ]

يمكن إعادة كتابة طيف الرتبة الثانية على النحو التالي: Sz(2)(ω)=-+هـأناωτج2(z(ت+τ)،z(ت))دτ،{\displaystyle S_{z}^{(2)}(\omega )=\int _{-\infty }^{+\infty }e^{i\omega \tau }C_{2}(z(t+\tau ),z(t))d\tau ,} أينج2(x،y)=xy-xy{\displaystyle C_{2}(x,y)=\langle xy\rangle -\langle x\rangle \langle y\rangle }هي التباين المشترك . تمثيل لـSz(2){\displaystyle S_{z}^{(2)}}وهو ما يعادل تعريف الكثافة الطيفية لـz{\displaystyle z}(طيف القدرة) كتحويل فورييه لدالة الارتباط الذاتي لـz(ت){\displaystyle z(t)}.

تاريخ الأطياف المتعددة

تمّ وضع أسس نظرية الأطياف ذات الرتب العليا في أعمال بريلينجر بدءًا من خمسينيات القرن العشرين. [ 1 ] وقدّم مندل ونيكياس أولى المراجعات الرئيسية حول هذا الموضوع في أوائل تسعينيات القرن العشرين، [ 3 ] [ 4 ] مما أدى إلى ظهور "مجموعة أدوات تحليل الأطياف ذات الرتب العليا" (HOSA)، وهي مكتبة برمجية تضمنت الطيف الثنائي وبعض إحصاءات الرتبة الرابعة. [ 5 ] ومع ذلك، صرّح بيركلوند في عام 2003: "نظريًا، يمكن تطبيق الأطياف المتعددة لحل العديد من المشكلات المهمة في معالجة الإشارات وتحليل البيانات. ولكن عمليًا، يُثبط المرء بسبب ضعف الخصائص الإحصائية لمعظم مُقدِّرات الأطياف المتعددة." [ 6 ] في عام 2026، تم اقتراح وتحليل مُقدِّرات غير متحيزة قائمة على العزوم للأطياف ذات الرتب العليا. [ 7 ] وهي مُطبَّقة في مكتبات برمجية مُسرَّعة بواسطة وحدة معالجة الرسومات (GPU) تعتمد على ArrayFire [ 8 ] ، أو بديلًا عن ذلك، تعتمد على PyTorch. [ 9 ]

الأرقام الخاصة بالأطياف المتعددة

يفترض تعريف الأطياف المتعددة عملية ثابتةz(ت){\displaystyle z(t)}بيانات من العالم الحقيقي معروفة بدقة زمنية لا نهائية في فترة زمنية لا نهائية.z(ت){\displaystyle z(t)}لذا، لا يُسمح إلا بحساب تقديرات الأطياف المتعددة. وتتميز هذه التقديرات بدقة طيفية محدودة ودقة محدودة.

المقدرات الكلاسيكية

تبدأ طرق فورييه المباشرة عادةً بتحويل فورييه للإشارة بأكملها، مما ينتج عنه دقة تردد عالية مع تباين كبير (خطأ). بعد ذلك، يمكن تطبيق دالة تنعيم لموازنة الدقة مع التباين. بدلاً من ذلك، يتم حساب متوسط ​​أطياف أجزاء من الإشارة (طريقة ويلش). [ 10 ]

تعتمد طريقة نافذة التأخير على تقدير التراكمات في المجال الزمني والتطبيق اللاحق لتحويلات فورييه.

تعتمد طرق التحليل متعدد المراحل على تقسيم الإشارة بأكملها إلى مراحل متعامدة. ثم يتم حساب متوسط ​​الأطياف لكل مرحلة على حدة. [ 6 ]

تعتمد الطرق الكلاسيكية على مقدرات التراكم التي تكون غير متحيزة تقاربياً فقط.

المقدرات غير المتحيزة القائمة على التراكمات

توفر المعادلات التالية تقديرًا غير متحيز ( مُقدِّرًا ) للأطياف المتعددة حتى الرتبة الرابعة. [ 7 ] معاملات فورييه لـz(ت){\displaystyle z(t)}أك=تيشمالج=0شمال-1زجzجهـ2πأناكج/شمال،{\displaystyle a_{k}={\frac {T}{N}}\sum _{j=0}^{N-1}g_{j}z_{j}e^{2\pi ikj/N}\,,} يتم حسابها من خلال العملية zج=z(جتي/شمال-ت0){\displaystyle z_{j}=z(jT/N-t_{0})}والتي تنقسم إلى نوافذ بطولتي{\displaystyle T}معشمال{\displaystyle N}عدد النقاط لكل نافذة. دالة النافذةزج=ز(جتي/شمال){\displaystyle g_{j}=g(jT/N)}يقلل من التسرب الطيفي.

يتم الحصول على الأطياف المتعددة التقريبية من تراكمات معاملات فورييه: متوسط: Sz(1)شمالج1(أ0)تيج=0شمال-1زجطيف القدرة: Sz(2)(ωك)شمالج2(أك،أك*)تيج=0شمال-1زجزج*الطيف الثنائي: Sz(3)(ωك،ωل)شمالج3(أك،أل،أك+ل*)تيج=0شمال-1زج2زج*ثلاثي الأطياف: Sz(4)(ωك،ωل،ωص)شمالج4(أك،أل،أص،أك+ل+ص*)تيج=0شمال-1زج3زج*.{\displaystyle {\begin{aligned}&{\text{المتوسط: }}&S_{z}^{(1)}&\approx {\frac {NC_{1}(a_{0})}{T\sum _{j=0}^{N-1}g_{j}}}\\\\[2pt]&{\text{طيف القدرة: }}&S_{z}^{(2)}(\omega _{k})&\approx {\frac {NC_{2}(a_{k},a_{k}^{*})}{T\sum _{j=0}^{N-1}g_{j}g_{j}^{*}}}\\\\[2pt]&{\text{طيف ثنائي: }}&S_{z}^{(3)}(\omega _{k},\omega _{l})&\approx {\frac { {NC_{3}(a_{k},a_{l},a_{k+l}^{*})}{T\sum _{j=0}^{N-1}g_{j}^{2}g_{j}^{*}}}\\\\[2pt]&{\text{Trispectrum: }}&S_{z}^{(4)}(\omega _{k},\omega _{l},\omega _{p})&\approx {\frac {NC_{4}(a_{k},a_{l},a_{p},a_{k+l+p}^{*})}{T\sum _{j=0}^{N-1}g_{j}^{3}g_{j}^{*}}}.\\\end{aligned}}} مسافات المواضع الطيفيةωك=2πك/تي{\displaystyle \omega _{k}=2\pi k/T}يتناقص مع زيادة طول النافذةتي{\displaystyle T}مما يؤدي إلى دقة طيفية أعلى. التراكماتجن{\displaystyle C_{n}}يمكن تقديرها بواسطة إحصائية k التي توفر مقدرات غير متحيزةجن{\displaystyle c_{n}} ( المُقدِّر ). [ 11 ] [ 12 ] [ 13 ] يتم إعطاء المُقدِّرات بواسطة [ 14 ]ج2(x،y)=مم-1(xy¯-x¯y¯)ج3(x،y،z)=م2(م-1)(م-2)(xyz¯-xy¯z¯-xz¯y¯-yz¯x¯+2x¯y¯z¯)ج4(x،y،z،w)=م2(م-1)(م-2)(م-3)[(م+1)xyzw¯-(م+1)(xyz¯w¯+xyw¯z¯+xzw¯y¯+yzw¯x¯)-(م-1)(xy¯zw¯+xz¯yw¯+xw¯yz¯)+2م(xy¯z¯w¯+xz¯y¯w¯+xw¯y¯z¯+yz¯x¯w¯+yw¯x¯z¯+zw¯x¯y¯)-6مx¯y¯z¯w¯]،{\displaystyle {\begin{aligned}&c_{2}(x,y)={\frac {m}{m-1}}({\overline {xy}}-{\overline {x}}\,{\overline {y}})\\\\[2pt]&c_{3}(x,y,z)={\frac {m^{2}}{(m-1)(m-2)}}({\overline {xyz}}-{\overline {xy}}\,{\overline {z}}-{\overline {xz}}\,{\overline {y}}-{\overline {yz}}\,{\overline {x}}+2{\overline {x}}\,{\overline {y}}\,{\overline {z}})\\\\[2pt]&c_{4}(x,y,z,w)={\frac {m^{2}}{(m-1)(m-2)(m-3)}}{\big [}(m+1){\overline {xyzw}}\\&\qquad \qquad \qquad \qquad \qquad -(m+1)({\overline {xyz}}\,{\overline {w}}+{\overline {xyw}}\,{\overline {z}}+{\overline {xzw}}\,{\overline {y}}+{\overline {yzw}}\,{\overline {x}})\\&\qquad \qquad \qquad \qquad \qquad -(m-1)({\overline {xy}}\,{\overline {zw}}+{\overline {xz}}\,\,{\overline {yw}}+{\overline {xw}}\,{\overline {yz}})\\&\qquad \qquad \qquad \qquad \qquad +2m({\overline {xy}}\,{\overline {z}}\,{\overline {w}}+{\overline {xz}}\,{\overline {y}}\,{\overline {w}}+{\overline {xw}}\,{\overline {y}}\,{\overline {z}}+{\overline {yz}}\,{\overline {x}}\,{\overline {w}}+{\overline {yw}}\,{\overline {x}}\,{\overline {z}}+{\overline {zw}}\,{\overline {x}}\,{\overline {y}})\\&\qquad \qquad \qquad \qquad \qquad -6m{\overline {x}}\,{\overline {y}}\,{\overline {z}}\,{\overline {w}}{\big ]},\end{aligned}}} أين(..)¯{\displaystyle {\overline {(..)}}}يشير إلى متوسطم{\displaystyle m}عينات.

التطبيقات

تُستخدم الأطياف الثنائية في التطبيقات الطبية، لا سيما في تحليل تخطيطات الدماغ الكهربائية [ 15 ] ، وفي الهندسة للكشف عن الأعطال في مكونات الآلات الدوارة [ 16 ] . كما أصبحت الأطياف المتعددة أداةً في الفيزياء الكمية لتحليل سجلات القياسات الكمية [ 17 ] . وتتراوح تطبيقات الأطياف المتعددة الكمية من النقل الكمي، مرورًا بالكشف عن الفوتون المفرد، وصولًا إلى قياسات ضوضاء الدوران [ 18 ] [ 19 ] [ 20 ] [ 21 ] .

مراجع

  1. 1 2 بريلينجر، ديفيد ر. (1965). "مقدمة في الأطياف المتعددة" . حوليات الإحصاء الرياضي . 36 (5): 1351-1374 . doi : 10.1214/aoms/1177699896 .
  2. أوبيلوده، نيلز؛ فريكه، كريستيان؛ فليندت، كريستيان؛ هولز، فرانك؛ هاوغ، رولف ج. (2012). "قياس إحصائيات التيار ذي التردد المحدود في ترانزستور أحادي الإلكترون" . نيتشر كوميونيكيشنز . 3 : 612. arXiv : 1201.2163 . doi : 10.1038/ncomms1620 .
  3. مندل، جيري م. (1991). "مقدمة في الإحصاءات ذات الرتبة العليا (الأطياف) في معالجة الإشارات ونظرية الأنظمة: نتائج نظرية وبعض التطبيقات" . وقائع معهد مهندسي الكهرباء والإلكترونيات . 79 (3): 278-305 . doi : 10.1109/5.75086 .
  4. نيكاس، سي إل (1993). "معالجة الإشارات باستخدام أطياف من الرتبة العليا" . مجلة معالجة الإشارات التابعة لمعهد مهندسي الكهرباء والإلكترونيات . 10 (3): 10-37 . doi : 10.1109/79.221324 .
  5. سوامي، أنانثرام؛ مندل، جيري م.؛ نيكاس، كريسوستوموس ل. (1993). مجموعة أدوات تحليل الطيف من الرتبة العليا: دليل المستخدم (ملف PDF) . كومبتون، كاليفورنيا: شركة يونايتد سيجنالز آند سيستمز.
  6. 1 2 بيركيلوند، ينجفي؛ هانسن، ألفريد. القوى، إدوارد ج. (2003). “مقدرات متعددة الأطياف”. معالجة الإشارات . 83 (3): 545-559 . دوى : 10.1016 / S0165-1684 (02)00484-X .
  7. 1 2 سيف، ماركوس؛ قربانيتيماد، أرمين؛ فاغنر، فابيان؛ هاغيله، دانيال (2026). "التقدير الصحيح للأطياف ذات الرتبة الأعلى: من التحديات النظرية إلى التنفيذ العملي متعدد القنوات في SignalSnap" . معالجة الإشارات الرقمية . 173 105893. doi : 10.1016/j.dsp.2026.105893 .
  8. "SignalSnap: Spectra at your fingertips" . GitHub . MarkusSifft . تم الاطلاع عليه في 26 فبراير 2026. SignalSnap عبارة عن مجموعة أدوات مفتوحة المصدر مكتوبة بلغة بايثون لتحليل الطيف عالي الرتبة للسلاسل الزمنية ، حيث تحسب بكفاءة أطياف القدرة، والأطياف الثنائية، والأطياف الثلاثية، وتدعم الحوسبة المعجلة بالأجهزة.
  9. "SignalSnap-PyTorch: إعادة كتابة SignalSnap باستخدام PyTorch، مع إضافة خيارات متعددة القنوات وتحسين الأداء" . GitHub . ArminGEtemad . تاريخ الاسترجاع: 26 فبراير 2026. تمت إعادة كتابة SignalSnap باستخدام PyTorch، مع خيارات متعددة القنوات وتحسين الأداء؛ ينقل هذا الإصدار الواجهة الخلفية من ArrayFire إلى PyTorch لتحسين سهولة الاستخدام والأداء.
  10. هينيتش، ميلفين جيه؛ كلاي، كريستيان إس. (1968). "تطبيق تحويل فورييه المنفصل في تقدير أطياف القدرة، والتماسك، والأطياف الثنائية للبيانات الجيوفيزيائية". مراجعات الجيوفيزياء . 6 (3): 347-363 . doi : 10.1029/RG006i003p00347 .
  11. فيشر، رونالد أ. (1928). "العزوم وعزوم الضرب لتوزيعات المعاينة". وقائع الجمعية الرياضية بلندن . السلسلة الثانية. 30 (1): 199-238 . doi : 10.1112/plms/s2-30.1.199 .
  12. كيندال، موريس ج. (1943). النظرية المتقدمة للإحصاء . المجلد 1 ( الطبعة الأولى). لندن: تشارلز جريفين وشركاه.  
  13. كوك، إم بي (1951). "إحصائيات k ثنائية المتغيرات وتراكمات توزيع المعاينة المشترك لها". Biometrika . 38 (1/2): 179–195 . doi : 10.1093/biomet/38.1-2.179 .
  14. شيفزيك، فابيان؛ هاجيل، دانيال (2019). "مُقدِّرات غير متحيزة جاهزة للاستخدام للعزوم التراكمية متعددة المتغيرات، بما في ذلك مُقدِّر يتفوق على غيره"x3{\displaystyle x^{3}}". arXiv : 1904.12154 [ math.ST ].
  15. سيجل، جيه سي؛ تشامبلين، إن بي (1994). "التحليل الطيفي الثنائي لتخطيط كهربية الدماغ: مراجعة". مجلة المراقبة السريرية . 10 (6): 392-404 . doi : 10.1007/BF01618421 . PMID 7836975 . 
  16. باركر، بي إي؛ كاسادا، تي إيه؛ فوربس، جي إل (2000). "تشخيص الأعطال باستخدام الكشف عن التغيرات الإحصائية في المجال الطيفي الثنائي". أنظمة ميكانيكية ومعالجة الإشارات . 14 (4): 561-570 . doi : 10.1006/mssp.2000.1299 .
  17. نوريس، لي م.؛ باز-سيلفا، جيراردو أ.؛ فيولا، لورينزا. "مطيافية ضوضاء الكيوبت لبيئات إزالة التماسك غير الغاوسية" . رسائل المراجعة الفيزيائية . 116 (15) 150503. arXiv : 1512.01575 . doi : 10.1103/PhysRevLett.116.150503 .
  18. إيماري، سي.؛ ماركوس، دي.؛ أغوادو، آر.؛ براندس، تي. "إحصائيات العد المعتمدة على التردد في الموصلات النانوية المتفاعلة". مجلة Physical Review B. 76 ( 16): 161404(R). arXiv : cond-mat/0703781 . doi : 10.1103/PhysRevB.76.161404 .
  19. هاجيل، دانيال؛ شيفزيك، فابيان (2018). "العزوم التراكمية، والتراكمات، وأطياف قياسات الضوضاء الكمومية المستمرة من الرتبة العليا". مجلة Physical Review B. 98 ( 20) 205143. arXiv : 1611.02077 . doi : 10.1103/PhysRevB.98.205143 .
  20. سيف، ماركوس؛ كورتزمان، أنيكا؛ كيرسكي، ينس؛ شوت، روديجر؛ ويك، أندرياس د.؛ جيلر، مارتن؛ لوركه، أكسل؛ هاجيله، دانيال (2021). "الأطياف الكمومية المتعددة لنمذجة وتقييم قياسات النقل الكمومي: نهج موحد لنظام القياس القوي والضعيف" . مجلة Physical Review Research . 3 (3) 033123. arXiv : 2011.07992 . doi : 10.1103/PhysRevResearch.3.033123 .
  21. سيف، ماركوس؛ كورتزمان، أ.؛ كيرسكي، ج.؛ شوت، ر.؛ لودفيج، أ.؛ ويك، أ.د.؛ لوركه، أ.؛ جيلر، م.؛ هاجيل، دانيال (2024). "نهج الأطياف المتعددة الكمومية لديناميكيات الباعثات الكمومية الوامضة عند معدلات الفوتونات المنخفضة دون تجميع: جعل كل فوتون مهمًا". مجلة Physical Review A. 109 ( 6) 062210. doi : 10.1103/PhysRevA.109.062210 .