قراءة
وضع العرض التقديمي

ما هو مرشح FIR؟

~13 دقيقة قراءة الدرس 1 من الوحدة 7

الترشيح: تشكيل الطيف الترددي

كل إشارة تصادفها في العالم الحقيقي تحتوي على مزيج من الترددات — بعضها مرغوب وبعضها غير مرغوب. يلتقط الميكروفون ضوضاء الغرفة إلى جانب صوت المتحدث. تحمل قراءة المستشعر كلاً من القياس الفيزيائي وتشويش 60 هرتز الكهربائي. يستقبل جهاز الراديو عشرات القنوات دفعةً واحدة. الترشيح هو العملية التي تُمرِّر ترددات معينة بشكل انتقائي بينما تُضعف غيرها، وهو من أكثر العمليات الأساسية في معالجة الإشارات الرقمية.

هناك عائلتان رئيسيتان من المرشحات الرقمية: FIR (الاستجابة النبضية المحدودة) وIIR (الاستجابة النبضية اللانهائية). تركّز هذه الوحدة على مرشحات FIR — الأبسط والأكثر قابلية للتنبؤ، والأقوى في كثير من النواحي. إن فهم آلية عمل مرشح FIR سيفتح أمامك كل ما يليه: طرق التصميم، وأنواع المرشحات، وخاصية الطور الخطي المذهلة.

الفكرة الجوهرية: لا تغذية راجعة

يحسب مرشح FIR ناتجه بوصفه مجموعاً موزوناً للعينة الحالية وعدد محدود من عينات الإدخال السابقة. هذا كل شيء. لا توجد تغذية راجعة — لا يُعاد إدخال الناتج في الحسابات المستقبلية. هذا الاختيار الهيكلي الواحد يمنح مرشحات FIR كل خواصها المرغوبة.

رياضياً، الناتج عند الزمن n هو:

معادلة الفروق لمرشح FIR
y[n] = \sum_{k=0}^{M} b_k\, x[n-k] = b_0 x[n] + b_1 x[n-1] + \cdots + b_M x[n-M]
y[n] هو الناتج، x[n−k] هي عينات الإدخال الحالية والـ M السابقة، وb_k هي معاملات المرشح (الأوزان). يمتلك المرشح M+1 معاملاً وترتيبه M.

المعاملات b_0, b_1, …, b_M هي مقابض التصميم. اختيارها بحكمة يُشكّل الاستجابة الترددية للمرشح إلى أي نطاق تمرير ونطاق إيقاف وعرض انتقالي تحتاجه.

الاستجابة النبضية هي المرشح

تذكّر من الوحدة 3 أن نظام LTI يتميز بالكامل باستجابته النبضية h[n]. بالنسبة لمرشح FIR، تساوي الاستجابة النبضية تماماً مجموعة المعاملات:

الاستجابة النبضية
h[n] = \begin{cases} b_n & 0 \le n \le M \\ 0 & \text{otherwise} \end{cases}
h[k] = b_k لـ k = 0, 1, …, M؛ وصفر في كل مكان آخر. للاستجابة النبضية طول محدود M+1 — لهذا السبب يُسمى المرشح "الاستجابة النبضية المحدودة".

بما أن h[n] ذات طول محدود، فإن الناتج هو ببساطة الالتواء بين الإدخال وهذا التسلسل المحدود: y[n] = x[n] * h[n]. لا يوجد شيء آخر. ناتج المرشح في أي لحظة هو حاصل ضرب نقطي متحرك بين آخر M+1 عينات إدخال وشعاع المعاملات.

الاستجابة الترددية

الاستجابة الترددية لمرشح FIR هي تحويل DTFT لاستجابته النبضية، محسوبةً على دائرة الوحدة (z = e^{jω}):

الاستجابة الترددية
H(e^{j\omega}) = \sum_{k=0}^{M} b_k\, e^{-j\omega k}
H(e^{jω}) متعددة حدود في e^{−jω}. يصف المقدار |H(e^{jω})| كيف يتضخم كل تردد أو يتضاءل؛ وتصف زاوية الطور ∠H(e^{jω}) تأخير الزمن المفروض على كل تردد.

على خلاف المرشح التناظري الذي يتقيد ترتيبه بتعقيد الدائرة، يمكن جعل مرشح FIR طويلاً بشكل عشوائي — مزيد من المعاملات يعني انتقالات أحدّ ونطاقات تمرير أكثر استواءً، بتكلفة حوسبة أعلى.

الاستقرار غير المشروط

من أهم المزايا العملية لمرشحات FIR أنها مستقرة دائماً. تذكّر معيار استقرار BIBO: النظام مستقر إذا وفقط إذا كانت جميع أقطابه داخل دائرة الوحدة تماماً. بما أن مرشح FIR لا يحتوي على تغذية راجعة، فإن دالة نقله متعددة حدود في z⁻¹ — لها أصفار لكن لا أقطاب (أو بمعنى آخر، كل أقطابها عند الأصل z = 0، داخل دائرة الوحدة تلقائياً).

لماذا يهمّ الاستقرار

يُنتج المرشح غير المستقر مخرجات تنمو بلا حدود حتى لإدخال محدود — كارثة في نظام صوتي أو تحكم آني. مرشحات IIR، التي تحتوي تغذية راجعة، يمكن أن تصبح غير مستقرة إذا لم تُصَمَّم معاملاتها بعناية أو إذا دفعت تأثيرات الطول المحدود للكلمة قطباً خارج دائرة الوحدة. مع مرشحات FIR، ببساطة لن تقلق بشأن هذا أبداً. غيّر أي معامل لأي قيمة ويبقى المرشح مستقراً.

لا أقطاب ← لا عدم استقرار. مرشحات FIR مستقرة BIBO بشكل غير مشروط.

خاصية الطور الخطي

ربما تكون أشهر خاصية لمرشحات FIR هي القدرة على تحقيق الطور الخطي. المرشح يمتلك طوراً خطياً إذا كانت استجابة طوره خطاً مستقيماً: ∠H(e^{jω}) = −αω لثابت α ما. الطور الخطي يعني أن كل مركبة ترددية تتأخر بنفس المقدار بالضبط (α عينة). شكل موجة الإشارة المارة عبر مرشح ذي طور خطي يُحفظ — فقط التوقيت يتحول.

يتحقق الطور الخطي حين تكون معاملات الاستجابة النبضية متناظرة (أو مضادة للتناظر):

متناظرة
h[n] = h[M − n]
تناظر مرآة حول المعامل المركزي. يُعطي تشويه طور صفرياً. الحالة الأكثر شيوعاً — تُستخدم في مرشحات FIR التمريرية المنخفضة والعالية والنطاقية.
مضادة للتناظر
h[n] = −h[M − n]
تناظر فردي حول المركز. يُضيف إزاحة طور 90° إلى الطور الخطي. تُستخدم للمفاضلات ومحولات هيلبرت.

شرط التناظر هذا قيد قوي: يُنصّف عدد الضربات المطلوبة (لأن 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 ذات الطور الخطي لتشكيل الاستجابة الترددية دون إدخال تشويه طور مسموع.
الاتصالات
تشكيل النبضات
مرشحات الجذر التربيعي لجيب تمام الارتفاع (RRC) هي مرشحات FIR تُستخدم في المودم الرقمي والراديو المحمول للحد من النطاق الترددي والتحكم في التداخل بين الرموز.
تعدد المعدلات
مكافحة التضمين والتصوير
مرشحات التخفيض والاستيفاء من الوحدة 11 هي دائماً تقريباً FIR: استقرارها غير المشروط وطورها الخطي يجعلانها مثالية لسلاسل تحويل معدل الأخذ.
الطبي الحيوي
معالجة ECG / EEG
تحتاج مرشحات ECG وEEG إلى طور خطي للحفاظ على شكل الموجة — تشويه طوري بسيط قد يجعل مركب QRS الطبيعي يبدو مرضياً. مرشحات 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 في كل مكان: معادلة الصوت، وتشكيل النبضات في الاتصالات، وترشيح تعدد المعدلات، ومعالجة الإشارات الطبية الحيوية.
التالي طريقة Windowed-Sinc نظرة عامة على الوحدة السابق تحليل الطيف الآني