الترشيح: تشكيل الطيف الترددي
كل إشارة تصادفها في العالم الحقيقي تحتوي على مزيج من الترددات — بعضها مرغوب وبعضها غير مرغوب. يلتقط الميكروفون ضوضاء الغرفة إلى جانب صوت المتحدث. تحمل قراءة المستشعر كلاً من القياس الفيزيائي وتشويش 60 هرتز الكهربائي. يستقبل جهاز الراديو عشرات القنوات دفعةً واحدة. الترشيح هو العملية التي تُمرِّر ترددات معينة بشكل انتقائي بينما تُضعف غيرها، وهو من أكثر العمليات الأساسية في معالجة الإشارات الرقمية.
هناك عائلتان رئيسيتان من المرشحات الرقمية: FIR (الاستجابة النبضية المحدودة) وIIR (الاستجابة النبضية اللانهائية). تركّز هذه الوحدة على مرشحات FIR — الأبسط والأكثر قابلية للتنبؤ، والأقوى في كثير من النواحي. إن فهم آلية عمل مرشح FIR سيفتح أمامك كل ما يليه: طرق التصميم، وأنواع المرشحات، وخاصية الطور الخطي المذهلة.
الفكرة الجوهرية: لا تغذية راجعة
يحسب مرشح FIR ناتجه بوصفه مجموعاً موزوناً للعينة الحالية وعدد محدود من عينات الإدخال السابقة. هذا كل شيء. لا توجد تغذية راجعة — لا يُعاد إدخال الناتج في الحسابات المستقبلية. هذا الاختيار الهيكلي الواحد يمنح مرشحات FIR كل خواصها المرغوبة.
رياضياً، الناتج عند الزمن n هو:
المعاملات b_0, b_1, …, b_M هي مقابض التصميم. اختيارها بحكمة يُشكّل الاستجابة الترددية للمرشح إلى أي نطاق تمرير ونطاق إيقاف وعرض انتقالي تحتاجه.
الاستجابة النبضية هي المرشح
تذكّر من الوحدة 3 أن نظام LTI يتميز بالكامل باستجابته النبضية h[n]. بالنسبة لمرشح FIR، تساوي الاستجابة النبضية تماماً مجموعة المعاملات:
بما أن h[n] ذات طول محدود، فإن الناتج هو ببساطة الالتواء بين الإدخال وهذا التسلسل المحدود: y[n] = x[n] * h[n]. لا يوجد شيء آخر. ناتج المرشح في أي لحظة هو حاصل ضرب نقطي متحرك بين آخر M+1 عينات إدخال وشعاع المعاملات.
الاستجابة الترددية
الاستجابة الترددية لمرشح FIR هي تحويل DTFT لاستجابته النبضية، محسوبةً على دائرة الوحدة (z = e^{jω}):
على خلاف المرشح التناظري الذي يتقيد ترتيبه بتعقيد الدائرة، يمكن جعل مرشح FIR طويلاً بشكل عشوائي — مزيد من المعاملات يعني انتقالات أحدّ ونطاقات تمرير أكثر استواءً، بتكلفة حوسبة أعلى.
الاستقرار غير المشروط
من أهم المزايا العملية لمرشحات FIR أنها مستقرة دائماً. تذكّر معيار استقرار BIBO: النظام مستقر إذا وفقط إذا كانت جميع أقطابه داخل دائرة الوحدة تماماً. بما أن مرشح FIR لا يحتوي على تغذية راجعة، فإن دالة نقله متعددة حدود في z⁻¹ — لها أصفار لكن لا أقطاب (أو بمعنى آخر، كل أقطابها عند الأصل z = 0، داخل دائرة الوحدة تلقائياً).
يُنتج المرشح غير المستقر مخرجات تنمو بلا حدود حتى لإدخال محدود — كارثة في نظام صوتي أو تحكم آني. مرشحات IIR، التي تحتوي تغذية راجعة، يمكن أن تصبح غير مستقرة إذا لم تُصَمَّم معاملاتها بعناية أو إذا دفعت تأثيرات الطول المحدود للكلمة قطباً خارج دائرة الوحدة. مع مرشحات FIR، ببساطة لن تقلق بشأن هذا أبداً. غيّر أي معامل لأي قيمة ويبقى المرشح مستقراً.
لا أقطاب ← لا عدم استقرار. مرشحات FIR مستقرة BIBO بشكل غير مشروط.خاصية الطور الخطي
ربما تكون أشهر خاصية لمرشحات FIR هي القدرة على تحقيق الطور الخطي. المرشح يمتلك طوراً خطياً إذا كانت استجابة طوره خطاً مستقيماً: ∠H(e^{jω}) = −αω لثابت α ما. الطور الخطي يعني أن كل مركبة ترددية تتأخر بنفس المقدار بالضبط (α عينة). شكل موجة الإشارة المارة عبر مرشح ذي طور خطي يُحفظ — فقط التوقيت يتحول.
يتحقق الطور الخطي حين تكون معاملات الاستجابة النبضية متناظرة (أو مضادة للتناظر):
شرط التناظر هذا قيد قوي: يُنصّف عدد الضربات المطلوبة (لأن b_k = b_{M−k}، ونصف المعاملات فقط فريدة) ويضمن أن طور الاستجابة الترددية خطي — وهو ما لا تستطيع مرشحات IIR تحقيقه بالضبط أبداً.
التكلفة الحوسبية
ثمن مرشحات FIR هو الحوسبة: إنتاج كل عينة ناتج يتطلب M+1 عملية ضرب وM جمعاً. مرشح تمريري منخفض من 101 معاملاً (M = 100) يتطلب 101 عملية ضرب-تراكم لكل عينة. عند صوت 48 كيلوهرتز، هذا 4.85 مليون عملية MAC في الثانية — تافهة لمعالج DSP أو متحكم حديث، لكنها اعتبار مهم في التطبيقات ذات معدلات الأخذ العالية جداً أو الطاقة المنخفضة جداً.
| طول المرشح | عدد المعاملات | MAC / عينة | الاستخدام النموذجي |
|---|---|---|---|
| قصير | 11 – 31 | 11 – 31 | تنعيم بسيط، تقوية مسبقة |
| متوسط | 51 – 127 | 51 – 127 | معادِل صوت، مكافحة التضمين |
| طويل | 255 – 1023 | 255 – 1023 | تشعّبات حادة، معادلة القناة |
| طويل جداً | 4096+ | 4096+ | تصحيح الغرفة، صدى الالتواء |
بالنسبة لمرشحات FIR الطويلة جداً (مئات إلى آلاف من المعاملات)، يصبح الالتواء المباشر مكلفاً. تستغل طريقتا الجمع المتداخل والحفظ المتداخل تحويل FFT لتطبيق الالتواء في نطاق التردد بتكلفة O(N log N) لكل كتلة — تسريع هائل عندما يكون المرشح أطول بكثير من كتلة الإدخال.
مرشحات FIR في التطبيق العملي
تظهر مرشحات FIR في كل مكان في معالجة الإشارات الرقمية:
- يحسب مرشح FIR ناتجه كمجموع موزون للعينة الحالية وM عينة إدخال سابقة — بلا تغذية راجعة ولا تكرار.
- معادلة الفروق هي y[n] = Σ b_k · x[n−k]، حيث المعاملات M+1 b_k هي معاملات التصميم.
- الاستجابة النبضية تساوي تسلسل المعاملات: h[k] = b_k. لها طول محدود — وهذا سبب التسمية "الاستجابة النبضية المحدودة".
- مرشحات FIR مستقرة BIBO بشكل غير مشروط — بلا أقطاب، لا شيء يُعرّضها لعدم الاستقرار، بصرف النظر عن قيم المعاملات.
- الاستجابة النبضية المتناظرة (h[n] = h[M−n]) تضمن الطور الخطي: كل تردد يتأخر بنفس المقدار، فيُحفظ شكل الموجة.
- التكلفة الحوسبية هي M+1 عملية MAC لكل عينة ناتج؛ تُنفَّذ مرشحات FIR الطويلة عبر طريقة الجمع المتداخل بتحويل FFT.
- تُستخدم مرشحات FIR في كل مكان: معادلة الصوت، وتشكيل النبضات في الاتصالات، وترشيح تعدد المعدلات، ومعالجة الإشارات الطبية الحيوية.