قراءة
وضع القصص

التسرّب الطيفي عملياً

~١٥ دقيقة قراءة الدرس ٤ من الوحدة ٥

مشكلة القياس

نظرياً، تحويلة DFT لجيبٍ نقيّ هي خانة (bin) واحدة غير صفرية. عملياً نادراً ما تكون كذلك: فطاقة النغمة الواحدة تنتشر عبر حزم كثيرة مجاورة، وهو ما يُسمّى التسرّب الطيفي. أرجع الدرس السابق التسرّب إلى سببه — نافذة الرصد المحدودة تضرب الإشارة ببوابة مستطيلة، والضرب في الزمن هو التواءٌ في التردد، فيُستبدل بكل خطّ طيفي حقيقي نسخةٌ منزاحة من طيف النافذة نفسها. وهذا الدرس عمّا يفعله ذلك بـالقياس: كيف يُفسد التسرّب السعة التي تقرؤها من خانة، وكيف تجعله يختفي حين تستطيع، وكيف تُحدّه حين لا تستطيع.

والرهان عمليّ. فمن يقرأ محلّل طيف يحتاج أن يعرف هل تقع القمة عند −40 dB لأن النغمة ضعيفة فعلاً، أم لأن تسرّباً من جارٍ قويّ رفع أرضية الضوضاء، أم لأن النغمة وقعت بين خانتين فقيست أقلّ بنحو 4 dB. الثلاثة كلها ظواهر تسرّب، ولكلٍّ منها علاج مختلف.

أخذ العيّنات المتوافق: متى يختفي التسرّب

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

شرط أخذ العيّنات المتوافق
f_0 = k\,\frac{f_s}{N}, \quad k \in \{0, 1, \dots, N-1\}
تقع النغمة عند f₀ على الخانة k تماماً — ولا تتسرّب إلى أي خانة أخرى — حين تكون f₀ مضاعفاً صحيحاً لتباعد الخانات f_s/N، أي حين يتّسع سجلّ الـ N عيّنة لعددٍ صحيح من الدورات الكاملة.

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

خسارة التموّج: تأثير سياج الأوتاد

لا تُبلّغ DFT عن الطيف إلا عند مراكز خاناتها — فهي تأخذ عيّنات من الطيف المستمر الكامن عند N نقطة منفصلة، كمن يرى مشهداً من خلف سياج أوتاد. والنغمة التي تقع بين خانتين تُقاس على حاشية الفصّ الرئيسي للنافذة لا على قمّته، فتُقرأ سعتها منخفضة. وأسوأ الحالات نغمةٌ تقع في منتصف المسافة بين خانتين تماماً.

للنافذة المستطيلة يبلغ هذا النقص في أسوأ الحالات 3.92 dB — فتُقرأ السعة عند 64% فقط من قيمتها الحقيقية — وهو خطأ كبير بما يكفي ليهمّ في أي قياس سعة. هذا النقص الحتمي المعتمد على التردد هو خسارة التموّج (وتُسمّى أيضاً خسارة سياج الأوتاد). وهو ليس ضوضاء، بل انحيازٌ يمكنك توقّعه من موقع النغمة بالنسبة للحزم. والنوافذ ذات الفصّ الرئيسي الأعرض والأكثر استواءً تتموّج أقل: تُنزل نافذة هان أسوأ الحالات إلى 1.42 dB، وتُصمَّم نوافذ القمّة المسطّحة (flat-top) خصّيصاً لإبقائها تحت 0.1 dB — ولهذا هي الخيار الافتراضي للقياس الدقيق للسعة رغم أن فصّها الرئيسي العريض يُفسد دقّة التمييز الترددي.

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

يحكم التسرّب أيضاً كيف تقيس DFT المنفّذة بنافذة الضوضاء والقدرة، وهذا هو الجسر إلى كثافة الطيف الاستطاعية في الوحدة 9، الدرس 1. فالنافذة مرشّح، وكل خانة من DFT تمرّر لا ترددها وحده بل حزمة ضوضاء حوله يشكّلها الفصّ الرئيسي والفصوص الجانبية. وعرض المرشّح المستطيل المثالي الذي يمرّر القدر نفسه من قدرة الضوضاء هو عرض النطاق النبضي المكافئ (ENBW):

عرض النطاق النبضي المكافئ
\text{ENBW} = N\,\frac{\sum_{n=0}^{N-1} w[n]^2}{\left(\sum_{n=0}^{N-1} w[n]\right)^2}
يعتمد ENBW، بوحدات خانات DFT، على معاملات النافذة w[n] وحدها. النافذة المستطيلة لها ENBW = خانة واحدة (وهو الأضيق الممكن)؛ وكل نافذة متناقصة أعرض منها.

وهو ما يتيح تحويل طيف مقدارٍ مُنفَّذ بنافذة إلى كثافة طيف استطاعية معايَرة بوحدة الواط لكل هرتز — الوحدة التي تستعملها الوحدة 9 دون أن تشتقّها. فمربّع مقدار خطٍّ واحد قدرة؛ وقسمته على معدّل أخذ العيّنات وطاقة النافذة يحوّله إلى كثافة:

كثافة الطيف من DFT مُنفَّذة بنافذة
\hat{S}_{xx}(f_k) = \frac{2\,\lvert X_w[k]\rvert^2}{f_s \sum_{n=0}^{N-1} w[n]^2}
الدوريّة المُعدّلة (modified periodogram): |X_w[k]|² مقسومةً على f_s وطاقة النافذة Σw[n]² تعطي كثافة الطيف أحادية الجانب بوحدة قدرة لكل هرتز (للحزم الداخلية). ودون حدّ Σw[n]² يكون التقدير منحازاً بحسب النافذة التي اخترتها.

وذات الـ ENBW التي تعاير أرضية الضوضاء هي سبب أن النافذة قليلة التسرّب تكلّفك حساسيةً للنغمات الضعيفة: فنشرُ ضوضاء خانةٍ على عرض خانتين يرفع أرضية الضوضاء بمقدار 10 log₁₀(2) ≈ 3 dB بالنسبة للنافذة المستطيلة. لذا فاختيار النافذة دائماً مقايضةٌ ثلاثية بين دقّة التمييز الترددي، ودقّة السعة، والقدرة على رؤية إشارة صغيرة بجوار كبيرة.

اختيار النافذة

المرجع الكلاسيكي هو مسح فريدريك هاريس عام 1978، الذي يجدول هذه المقاييس لعشرات النوافذ. والخمس أدناه تغطّي المدى المفيد من «بلا نافذة» إلى «كبح الفصوص الجانبية المتطرّف»:

مقاييس أداء النوافذ (هاريس، 1978)
النافذةالفصّ الرئيسي −3 dB (حزم)ذروة الفصّ الجانبي (dB)خسارة التموّج (dB)ENBW (حزم)
مستطيلة0.89−133.921.00
هان (Hann)1.44−321.421.50
هامنغ (Hamming)1.30−431.751.36
بلاكمان (Blackman)1.68−581.101.73
بلاكمان-هاريس (رباعية الحدود)1.90−920.832.00

اقرأ الجدول مقايضةً واحدة تُمسح من الأعلى إلى الأسفل: كلما تناقصت النافذة أشدّ، هبطت ذروة الفصّ الجانبي (−13 dB إلى −92 dB، فتصير النغمة الضعيفة مرئيّة بجوار القويّة) وتقلّصت خسارة التموّج — لكن يتّسع الفصّ الرئيسي وينمو ENBW، فتسوء دقّة التمييز الترددي والحساسية للضوضاء. وصفّ بلاكمان-هاريس هو الطرف قليل التسرّب: فصوصه الجانبية عند −92 dB قادرة على كشف نغمة تقلّ 90 dB عن جارتها، بثمن فصٍّ رئيسي يزيد على ضعف عرض المستطيلة.

استراتيجيات عملية

إن تحكّمت بالمنبع
خذ العيّنات متوافقاً
اربط المولّد والمرقِّم بمرجع زمني واحد واختر N بحيث تكون النغمة محاذية لخانة. يختفي التسرّب بلا ثمنٍ من النافذة — أسلوب مختبر المعايرة.
لدقّة السعة
استعمل نافذة قمّة مسطّحة
خسارة تموّج تحت 0.1 dB أينما وقعت النغمة. فصّها الرئيسي عريض فتسوء لتمييز النغمات المتقاربة — لكنها الخيار الصحيح حين يجب أن يكون المستوى دقيقاً.
نغمة ضعيفة بجوار قويّة
استعمل فصوصاً جانبية منخفضة
بلاكمان-هاريس، أو نافذة كايزر بـ β عالية. فصوص عند −90 dB تمنع حاملةً قويّة من دفن جارتها الضعيفة، بمقايضة الدقّة الترددية بالمدى الديناميكي.
للمستويات المطلقة
صحّح لأثر النافذة
اقسم على كسب النافذة المتوافق (Σw[n]) وعلى ENBW قبل أن تذكر قدرةً أو كثافة — وإلا كانت القراءة منحازة بالنافذة التي صادف أن اخترتها.

افترض هذا كلّه أن حساب DFT نفسها ميسور. والوحدة 6 ترفع هذا الافتراض: تحويلة فورييه السريعة تُنزل كلفة DFT من O(N²) إلى O(N log N)، وهو ما يجعل التحليل الطيفي الآنيّ — وكل ما في هذا الدرس — عملياً.

أبرز النقاط
  • التسرّب الطيفي بصمة النافذة المحدودة: يُستبدل بالخطّ الحقيقي نسخةٌ من طيف النافذة، فتنتشر طاقة نغمةٍ واحدة عبر حزم كثيرة.
  • أخذ العيّنات المتوافق — نغمة عند f₀ = k·f_s/N بالضبط — يُنتج تسرّباً صفرياً، لأن كل خانة أخرى تقع على صفرٍ من طيف النافذة. ولا ينجح إلا حين تتحكّم بالإشارة.
  • خسارة التموّج نقصٌ متوقّع في السعة للنغمات الواقعة بين الحزم: 3.92 dB في أسوأ حالات المستطيلة، و1.42 dB لهان، وتحت 0.1 dB للقمّة المسطّحة.
  • عرض النطاق النبضي المكافئ (ENBW = N·Σw²/(Σw)²) عرضُ المرشّح المستطيل الذي يمرّر قدرة الضوضاء نفسها؛ وهو يعاير DFT المنفّذة بنافذة إلى كثافة طيف بوحدة W/Hz.
  • كل نافذة تقايض دقّة التمييز الترددي بمستوى الفصوص الجانبية ودقّة السعة والحساسية للضوضاء؛ وجدول هاريس 1978 يقيس المقايضة لكلٍّ منها.
  • القاعدة العملية: خذ العيّنات متوافقاً إن استطعت؛ والقمّة المسطّحة للسعة؛ وبلاكمان-هاريس للمدى الديناميكي؛ وصحّح دوماً للكسب المتوافق ولـ ENBW قبل ذكر قدرة مطلقة.
التالي لماذا تحويل DFT بطيء — مشكلة O(N²) نظرة عامة على الوحدة السابق دقة التردد والتنويف