التخفيض (Decimation) | DSP-101

تقليل معدل البيانات

لكل إشارة رقمية معدل أخذ عينات — عدد العينات المُجمَّعة في الثانية. كثيرًا ما يكون هذا المعدل أعلى مما يتطلبه التطبيق. جهاز الاستقبال اللاسلكي عريض النطاق يرقمن بـ 100 ميغاهرتز رغم أن الإشارة المطلوبة لا تحتل سوى 200 كيلوهرتز. محطة معالجة الصوت تسجّل بـ 96 كيلوهرتز وإن كان الناتج النهائي سيُشغَّل بـ 44.1 كيلوهرتز. في كلا الحالين، حمل العينات الزائدة يُضيّع الذاكرة والتخزين وطاقة المعالجة. التخفيض (Decimation) هو العملية المنهجية لتقليل معدل أخذ العينات بحذف العينات بطريقة رياضية مضبوطة.

التخفيض بعامل M يعني الاحتفاظ بكل عينة M-th وحذف الباقي؛ والناتج إشارة بمعدل 1/M من المعدل الأصلي. لكن ثمة محذور حاسم: لا يمكن ببساطة التخلص من العينات دون إزالة الترددات التي ستُولّد الالتباس الترددي (Aliasing) في الإشارة ذات المعدل المنخفض. الحصول على هذه الخطوة بصورة صحيحة هو جوهر التخفيض.

الفكرة الأساسية

التخفيض بعامل M = مرشح تمرير منخفض لمنع الالتباس يليه الاحتفاظ بكل عينة M-th. المرشح إلزامي — إغفاله يُدخل التباسًا يُفسد الإشارة بشكل دائم.

التخفيض الساذج: ما الذي يحدث؟

افترض أن لديك إشارة بمعدل f_s = 8 كيلوهرتز وتريد تخفيضها إلى 2 كيلوهرتز (M = 4). النهج الساذج: احتفظ بالعينة 0، وأسقط العينات 1 و2 و3، ثم احتفظ بالعينة 4، وهكذا. تسلسل الإخراج هو y[n] = x[4n].

المشكلة هي نظرية نايكويست المطبَّقة على المعدل الجديد. عند 2 كيلوهرتز، لا يمكن تمثيل سوى الترددات الأقل من 1 كيلوهرتز دون التباس. لكن الإشارة الأصلية بـ 8 كيلوهرتز قد تحتوي طاقة حتى 4 كيلوهرتز. عندما تحتفظ بكل عينة رابعة، يضغط محور التردد بمقدار 4: مكوّن عند 3 كيلوهرتز يُصبح 3000/4 = 750 هرتز في الإخراج — فوق مكوّن 750 هرتز الأصيل. يختلطان بشكل لا يمكن فصله، والتلوث لا رجعة منه؛ لا معالجة لاحقة تستطيع إلغاء الالتباس متى وقع.

التخفيض

y[n] = x[Mn]

التسلسل المُخفَّض y[n] يُؤخذ بالاحتفاظ بكل عينة M-th من x[n]. في نطاق التردد يُضغط الطيف بمقدار M وتُجمع M نسخة — مما يُسبب التباسًا إن لم تكن الإشارة محدودة النطاق عند π/M.

مرشح مكافحة الالتباس

الحل واضح من حيث المبدأ: قبل التخفيض بعامل M، طبّق مرشح تمرير منخفض بتردد قطع f_c = f_s/(2M). يُزيل هذا كل طاقة عند حد نايكويست الجديد أو فوقه، فلا يبقى شيء يُسبب التباسًا. يتكون المُخفِّض الكامل من مرحلتين متتاليتين: أولًا المرشح التمريري المنخفض، ثم منتقي العينات.

تردد القطع

f_c = \frac{f_s}{2M} \quad\Leftrightarrow\quad \omega_c = \frac{\pi}{M}

يجب أن يقطع مرشح مكافحة الالتباس عند تردد نايكويست الخاص بمعدل الإخراج أو دونه — أي f_c = f_s/(2M)، وهو نصف معدل أخذ العينات للإخراج f_s/M. (وافصل بين المصطلحين: معدَّل نايكويست هو أدنى معدل لأخذ العينات وقيمته 2f_max؛ أما تردد نايكويست فهو أعلى تردد قابل للتمثيل، وقيمته نصف معدل أخذ العينات. ولو نصّفنا تردد نايكويست هنا لهبط تردد القطع بمقدار أوكتاف كامل.) وبالتردد الرقمي المُعيَّر (راديان/عينة على المعدل الأصلي)، تردد القطع هو π/M.

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

قاعدة تصميم المرشح

للتخفيض بعامل M، صمّم مرشح FIR تمريريًا منخفضًا بتردد قطع ω_c = π/M (راديان/عينة على معدل الإدخال). رتبة المرشح تحدد المفاضلة بين حدة انتقالية النطاق وتكلفة الحساب: مزيد من النقرات ← تسيُّل أحدّ لكن ضرب أكثر لكل عينة إدخال.

ما الذي يحدث في نطاق التردد

منظور نطاق التردد يجعل التخفيض حدسيًا. الإشارة التي يحتل طيفها النطاق [0, π/M] راديان/عينة (بعد إزالة كل شيء فوقه بالمرشح) — عندما نحتفظ بكل عينة M-th، يكون DTFT للإخراج هو مجموع M نسخة منزاحة من طيف الإدخال، كل منها مقسومة على M:

طيف الإشارة المُخفَّضة

Y\!\left(e^{j\omega}\right) = \frac{1}{M}\sum_{k=0}^{M-1} X\!\left(e^{j(\omega - 2\pi k)/M}\right)

طيف الإخراج Y(e^jω) هو مجموع M نسخة منزاحة ترددياً من طيف الإدخال المُرشَّح، تفصل بينها مسافة 2π/M. إن كان الإدخال محدود النطاق عند π/M، لا تتداخل هذه النسخ — لا التباس.

لأن مرشح مكافحة الالتباس يضمن أن طيف الإدخال صفر فوق π/M، فإن النسخ M في المجموع لا تتداخل. طيف الإخراج هو نسخة مُمتدة من نطاق العبور للإشارة — الترددات التي كانت في [0, π/M] تمتد الآن لتشغل [0, π] عند المعدل الجديد (الأقل). الإشارة تُضغط طيفيًا: نفس المعلومات تشغل عرض النطاق الكامل للإشارة الناتجة.

1/M

معدل الإخراج نسبةً إلى الإدخال

π/M

الحد الأقصى لتردد القطع الآمن (راديان/عينة)

M×

عامل الضغط الطيفي في الإخراج

التنفيذ الكفء: التحليل متعدد الأوجه

مُخفِّض ساذج يحسب جميع عينات إخراج المرشح ثم يتخلى عن M−1 من كل M — هدر واضح. التحليل متعدد الأوجه (Polyphase Decomposition) يستغل بنية التخفيض لتفادي حساب العينات التي ستُحذف.

مرشح FIR طوله N يُقسَّم إلى M مرشحات فرعية (فروع متعددة الأوجه)، كل منها طوله تقريبًا N/M. الفرع k يحتوي معاملات المرشح عند المؤشرات k، k+M، k+2M، ... . للمُخفِّض، يُلزم فرع واحد فقط لكل عينة إخراج — يتقلص الحساب بعامل M. تنخفض التكلفة الإجمالية من N عملية ضرب لكل عينة إدخال إلى تقريبًا N/M لكل عينة إخراج.

التحليل متعدد الأوجه

E_k[n] = h[nM + k], \quad k = 0, 1, \ldots, M-1

يُحلَّل المرشح h[n] إلى M فرع متعدد الأوجه E_k[n] = h[nM + k]. كل فرع يعمل على المعدل المنخفض للإخراج، مما يُقلل الحساب بعامل M.

اعتبارات عملية

التخفيض الصحيح مقابل النسبي: يجب أن يكون العامل M عددًا صحيحًا موجبًا. لتخفيض بنسبة غير صحيحة (مثلًا 44.1 كيلوهرتز إلى 32 كيلوهرتز = 320:441)، تجمع رفع العينات والتخفيض: ارفع أولًا بعامل L، ثم خفِّض بعامل M، مختارًا L وM بحيث تكون L/M هي النسبة المطلوبة. هذا ما يُعرف بـ التحويل النسبي لمعدل العينات.

التخفيض متعدد المراحل: لقيم M الكبيرة (مثلًا M = 64)، يستلزم مرشح FIR أحادي المرحلة معاملات كثيرة لتلبية متطلبات نطاق الانتقال الضيقة. النهج متعدد المراحل يُسلسل مُخفِّضات أصغر (مثلًا تخفيض بـ 4 ثم 4 ثم 4) بمرشحات أبسط في كل مرحلة، مما يُقلل التكلفة الإجمالية بشكل كبير.

مرشحات CIC للتخفيض الشديد: مرشحات CIC (Cascaded Integrator-Comb) هي بنية حسابية خالية من الضرب، مثالية لنسب تخفيض عالية جدًا (M = 256 أو أكثر). تُستخدم شيوعًا في مُحوِّلات ΣΔ وواجهات SDR الأمامية.

الدرس القادم: الاستيفاء (Interpolation/Upsampling) — العملية العكسية: إدراج عينات لرفع المعدل، مع مرشح مكافحة الصور الطيفية الذي يمنع ظهور نسخ طيفية في الإخراج.

أبرز النقاط

التخفيض بعامل M يُقلل معدل أخذ العينات بعامل M بالاحتفاظ بكل عينة M-th؛ وهو الأداة الأساسية لتقليل معدل البيانات في معالجة الإشارات متعددة المعدلات.
التخفيض الساذج يُسبب التباسًا: تتداخل النسخ الطيفية وتُفسد الإشارة بشكل دائم. مرشح التمرير المنخفض لمكافحة الالتباس (قطع π/M) إلزامي قبل التخفيض.
في نطاق التردد، التخفيض بعامل M يضغط الطيف بمقدار M ويجمع M نسخة منزاحة؛ إن كانت الإشارة محدودة النطاق عند π/M، لا تتداخل النسخ.
التحليل متعدد الأوجه يُلغي الحساب المُهدَر بتشغيل المرشح على المعدل المنخفض للإخراج، مما يُخفض التكلفة بعامل M تقريبًا.
نسب التخفيض الكبيرة تُنفَّذ بكفاءة على مراحل متعاقبة؛ مرشحات CIC تتعامل مع التخفيض الشديد (M في المئات) بدون ضرب.
التخفيض هو النصف الأول من تحويل معدل العينات؛ النسب غير الصحيحة تجمع التخفيض مع الاستيفاء للتحويل بين معدلات عينات اعتباطية.