تصميم الفلتر

تصميم المرشحات هو عملية تصميم مرشح لمعالجة الإشارات يلبي مجموعة من المتطلبات، قد يكون بعضها متعارضًا. والهدف هو إيجاد تصميم للمرشح يلبي كل متطلب بدرجة مقبولة.

يمكن وصف عملية تصميم المرشح بأنها مسألة تحسين. يمكن أتمتة بعض أجزاء عملية التصميم، ولكن قد يكون من الضروري الاستعانة بمصمم خبير للحصول على نتيجة جيدة.

يُعد تصميم المرشحات الرقمية موضوعًا معقدًا. [ 1 ] على الرغم من سهولة فهم المرشحات وحسابها، إلا أن التحديات العملية لتصميمها وتنفيذها كبيرة وتخضع لبحوث متقدمة.

متطلبات التصميم النموذجية

المتطلبات النموذجية التي يمكن أخذها في الاعتبار في عملية التصميم هي:

دالة التردد

يُعدّ استجابة التردد المطلوبة معيارًا مهمًا . ويحدد انحدار منحنى الاستجابة وتعقيده رتبة المرشح وجدواه.

يحتوي المرشح التكراري من الدرجة الأولى على مُكوّن واحد فقط يعتمد على التردد. هذا يعني أن ميل استجابة التردد يقتصر على 6 ديسيبل لكل أوكتاف . وهذا غير كافٍ للعديد من الأغراض. ولتحقيق ميول أكثر حدة، يلزم استخدام مرشحات من درجات أعلى.

فيما يتعلق بدالة التردد المطلوبة، قد تكون هناك أيضًا دالة ترجيح مصاحبة ، تصف، لكل تردد، مدى أهمية أن تقترب دالة التردد الناتجة من الدالة المطلوبة.

من الأمثلة النموذجية لدالة التردد ما يلي:

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

تأخير المرحلة والمجموعة

  • يُمرر مرشح تمرير الكل جميع الترددات دون تغيير، ولكنه يُغير طور الإشارة. يمكن استخدام مرشحات من هذا النوع لمعادلة تأخير المجموعة في المرشحات التكرارية. كما يُستخدم هذا المرشح في تأثيرات الفايزر .
  • محول هيلبرت هو مرشح تمرير كلي محدد يمرر الموجات الجيبية بسعة غير متغيرة ولكنه يحرك طور كل موجة جيبية بمقدار ±90 درجة.
  • مرشح التأخير الجزئي هو مرشح تمرير كامل له تأخير مجموعة أو طور محدد وثابت لجميع الترددات.

استجابة النبضة

توجد علاقة مباشرة بين دالة التردد للمرشح واستجابته النبضية: فالأولى هي تحويل فورييه للثانية. وهذا يعني أن أي شرط على دالة التردد هو شرط على الاستجابة النبضية، والعكس صحيح.

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

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

السببية

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

يمكن أن تكون المرشحات التي لا تعمل في الوقت الفعلي (مثل مرشحات معالجة الصور) غير سببية. وقد تُصمم هذه المرشحات بحيث يكون تأخيرها معدوماً.

استقرار

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

المنطقة

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

وفقًا لعلاقة عدم اليقين في تحويل فورييه، يجب أن يتجاوز حاصل ضرب عرض دالة استجابة النبضة للمرشح وعرض دالة التردد الخاصة به قيمة ثابتة معينة. هذا يعني أن أي شرط يتعلق بموقعية المرشح يستلزم أيضًا حدًا أقصى لعرض دالة التردد الخاصة به. وبالتالي، قد لا يكون من الممكن تلبية متطلبات موقعية دالة استجابة النبضة للمرشح ومتطلبات موقعية دالة التردد الخاصة به في آنٍ واحد. وهذا مثال نموذجي على المتطلبات المتناقضة.

التعقيد الحسابي

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

توجد عدة طرق يمكن من خلالها أن تختلف التعقيدات الحسابية للمرشح. على سبيل المثال، تتناسب رتبة المرشح تقريبًا مع عدد العمليات. وهذا يعني أنه باختيار مرشح ذي رتبة منخفضة، يمكن تقليل وقت الحساب.

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

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

اعتبارات أخرى

يجب أيضاً تحديد كيفية تطبيق المرشح:

المرشحات التناظرية

يتم تغطية تصميم المرشحات التناظرية الخطية في الغالب في قسم المرشحات الخطية .

المرشحات الرقمية

تُصنف المرشحات الرقمية إلى أحد شكلين أساسيين، وفقًا لكيفية استجابتها لنبضة وحدة واحدة :

معدل أخذ العينات

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

مضاد التعرج

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

الأساس النظري

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

  • مبدأ عدم اليقين بين مجالي الزمن والتردد
  • نظرية امتداد التباين
  • السلوك التقاربي لمجال واحد مقابل الانقطاعات في المجال الآخر

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

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

نظرية امتداد التباين

يتركσs2{\displaystyle \sigma _{s}^{2}}ليكن تباين إشارة الإدخال، وليكنσو2{\displaystyle \sigma _{f}^{2}}ليكن تباين المرشح. تباين استجابة المرشح،σر2{\displaystyle \sigma _{r}^{2}}، ثم يتم تحديده بواسطة

σر2{\displaystyle \sigma _{r}^{2}}=σs2{\displaystyle \sigma _{s}^{2}}+σو2{\displaystyle \sigma _{f}^{2}}

هذا يعني أنσر>σو{\displaystyle \sigma _{r}>\sigma _{f}}وهذا يعني أن تحديد موقع خصائص مختلفة، مثل النبضات أو الخطوات في استجابة المرشح، محدود بعرض المرشح في مجال الإشارة. إذا طُلب تحديد موقع دقيق، فنحن بحاجة إلى مرشح ذي عرض صغير في مجال الإشارة، وبناءً على مبدأ عدم اليقين، لا يمكن أن يكون عرضه في مجال التردد صغيرًا بشكل تعسفي.

الانقطاعات مقابل السلوك التقاربي

لتكن f(t) دالة، ولتكنF(ω){\displaystyle F(\omega )}ليكن تحويل فورييه الخاص بها. هناك نظرية تنص على أنه إذا كانت المشتقة الأولى للدالة F غير المتصلة من الرتبةن0{\displaystyle n\geq 0}إذن، فإن الدالة f لها انحلال تقاربي مثلت-ن-1{\displaystyle t^{-n-1}}.

ومن نتائج هذه النظرية أن دالة التردد للمرشح يجب أن تكون سلسة قدر الإمكان للسماح لاستجابته النبضية بأن يكون لها اضمحلال سريع، وبالتالي عرض قصير.

المنهجية

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

هناك طريقة أخرى لإيجاد مرشح FIR منفصل، وهي تحسين المرشح الموصوف في بحث كنوتسون وآخرون، والذي يقلل من تكامل مربع الخطأ بدلاً من قيمته القصوى. في شكله الأساسي، يتطلب هذا النهج وجود دالة تردد مثالية للمرشح.Fأنا(ω){\displaystyle F_{I}(\أوميغا )}يتم تحديدها مع دالة ترجيح التردددبليو(ω){\displaystyle W(\omega )}ومجموعة من الإحداثياتxك{\displaystyle x_{k}}في مجال الإشارة حيث توجد معاملات المرشح.

دالة الخطأε{\displaystyle \varepsilon }يُعرَّف بأنه

ε=دبليو(Fأنا-F{و})2{\displaystyle \varepsilon =\|W\cdot (F_{I}-{\mathcal {F}}\{f\})\|^{2}}

أينو(x){\displaystyle f(x)}هو المرشح المنفصل وF{\displaystyle {\mathcal {F}}}هي تحويل فورييه المنفصل زمنيًا المعرف على مجموعة الإحداثيات المحددة. المعيار المستخدم هنا، رسميًا، هو المعيار المعتاد علىل2{\displaystyle L^{2}}المسافات. هذا يعني أنε{\displaystyle \varepsilon }يقيس الانحراف بين دالة التردد المطلوبة للمرشح،Fأنا{\displaystyle F_{I}}، ووظيفة التردد الفعلية للمرشح المُنفذ،F{و}{\displaystyle {\mathcal {F}}\{f\}}ومع ذلك، يخضع الانحراف أيضًا لدالة الترجيح.دبليو{\displaystyle W}قبل حساب دالة الخطأ.

بمجرد تحديد دالة الخطأ، يتم تحديد المرشح الأمثل بواسطة المعاملاتو(x){\displaystyle f(x)}مما يقللε{\displaystyle \varepsilon }يمكن تحقيق ذلك عن طريق حل مسألة المربعات الصغرى المناظرة. عمليًا،ل2{\displaystyle L^{2}}يجب تقريب المعيار عن طريق جمع مناسب على نقاط منفصلة في مجال التردد. ومع ذلك، بشكل عام، ينبغي أن يكون عدد هذه النقاط أكبر بكثير من عدد المعاملات في مجال الإشارة للحصول على تقريب مفيد.

التحسين المتزامن في كلا المجالين

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

انظر أيضاً

مراجع

  1. فالديز، م. إي. "المرشحات الرقمية" . شبكات جي آر إم . تم الاسترجاع في 13 يوليو 2020 .
  2. ستوري، مايك (سبتمبر 1997). "تفسير مقترح لبعض الاختلافات المسموعة بين المواد الصوتية ذات معدل العينة العالي والمواد الصوتية ذات معدل العينة التقليدي" (ملف PDF) . شركة dCS المحدودة. مؤرشف (ملف PDF) من الأصل في 28 نوفمبر 2009.
  3. روبجونز، هيو (أغسطس 2016). "دقة MQA في المجال الزمني وجودة الصوت الرقمي" . soundonsound.com . ساوند أون ساوند. مؤرشف من الأصل في 10 مارس 2023.
  4. بروكيس، جون ج.؛ مانولاكيس، ديميتريس ج. (2007). معالجة الإشارات الرقمية ( الطبعة الرابعة). أبر سادل ريفر، نيوجيرسي: بيرسون برنتيس هول. الصفحات 678-699 . ISBN   978-0-13-187374-2.
  5. رابينر، لورانس ر.، وغولد، برنارد، 1975: نظرية وتطبيق معالجة الإشارات الرقمية (إنجلوود كليفس، نيو جيرسي: برنتيس هول، إنك.) ISBN 0-13-914101-4

فهرس

  • أ. أنطونيو (1993). المرشحات الرقمية: التحليل والتصميم والتطبيقات (  الطبعة الثانية). ماكجرو هيل، نيويورك، نيويورك. ISBN 978-0-07-002117-4.
  • أ. أنطونيو (2006). معالجة الإشارات الرقمية: الإشارات والأنظمة والمرشحات . ماكجرو هيل، نيويورك، نيويورك. ISBN 978-0-07-145424-7.
  • SWA Bergen; A. Antoniou (2005). "تصميم مرشحات رقمية غير تكرارية باستخدام دالة النافذة فوق الكروية" . مجلة EURASIP لمعالجة الإشارات التطبيقية . 2005 (12): 1910. Bibcode : 2005EJASP2005...44B . doi : 10.1155/ASP.2005.1910 .
  • أ. ج. ديكزكي (أكتوبر 1972). "توليف المرشحات الرقمية المتكررة باستخدام معيار الحد الأدنى للخطأ p". معاملات IEEE في الصوتيات والإلكترونيات الصوتية . 20 (4): 257-263 . Bibcode : 1972ITAuE..20..257D . doi : 10.1109/TAU.1972.1162392 .
  • جيه كيه كايزر (1974). "تصميم مرشح رقمي غير تكراري باستخدام دالة النافذة I 0 -sinh". وقائع ندوة IEEE الدولية لنظرية الدوائر (ISCAS74) لعام 1974. سان فرانسيسكو، كاليفورنيا. الصفحات 20-23 . 
  • هـ. كنوتسون؛ م. أندرسون؛ جيه ويكلوند (يونيو 1999). “تصميم الفلتر المتقدم”. بروك. الندوة الاسكندنافية حول تحليل الصور، كانجيرلوسواك، جرينلاند .
  • إس كيه ميترا (1998). معالجة الإشارات الرقمية: منهج قائم على الحاسوب . ماكجرو هيل، نيويورك، نيويورك. رقم ISBN 978-0-07-286546-2.
  • A.V. Oppenheim; R.W. Schafer; J.R. Buck (1999). Discrete-Time Signal Processing. Prentice-Hall, Upper Saddle River, NJ. ISBN 978-0-13-754920-7.
  • T.W. Parks; J.H. McClellan (March 1972). "Chebyshev Approximation for Nonrecursive Digital Filters with Linear Phase". IEEE Transactions on Circuit Theory. 19 (2): 189–194. Bibcode:1972ITCT...19..189P. doi:10.1109/TCT.1972.1083419.
  • L.R. Rabiner; J.H. McClellan; T.W. Parks (April 1975). "FIR Digital Filter Design Techniques Using Weighted Chebyshev Approximation". Proceedings of the IEEE. 63 (4): 595–610. Bibcode:1975IEEEP..63..595R. doi:10.1109/PROC.1975.9794. S2CID 12579115.