جسر بين التناظري والرقمي
قدّم الدرس السابق ثلاث عائلات كلاسيكية من المرشحات التناظرية — باترورث، وتشيبيشيف، والإهليليجي — كل منها يتميّز بتوزيع أقطاب وأصفار محددة في مستوى s. والسؤال الآن: كيف نحوّل المرشح التناظري H_a(s) إلى مرشح رقمي H(z)؟ والإجابة التي تُستخدم في كل تصميم عملي تقريباً هي التحويل الثنائي الخطي.
التحويل الثنائي الخطي هو تعويض يُربط كل نقطة في مستوى s بنقطة فريدة على دائرة الوحدة في مستوى z — والعكس صحيح. يضمن هذا التحويل أن المرشح التناظري المستقر (أقطاب في النصف الأيسر من مستوى s) يُحوَّل دائماً إلى مرشح رقمي مستقر (أقطاب داخل دائرة الوحدة). لا يوجد تراكب ترددي — وهذه ميزة حاسمة على طريقة الاستجابة النبضية الثابتة الأقدم منها.
صيغة التعويض
يُعوّض التحويل الثنائي الخطي المتغيّر الترددي المركّب s بدالة في z:
هذا التعويض جبري ودقيق — لا توجد تقريب في التحويل نفسه. العيب الوحيد الذي يُدخله هو تشويه غير خطي لمحور التردد، يجب تعويضه.
تشويه التردد
تقييم التعويض الثنائي الخطي على دائرة الوحدة (z = e^{jω}) يكشف كيف تتحوّل الترددات التناظرية Ω إلى ترددات رقمية ω:
هذا الضغط هو مصدر تشويه التردد. بالقرب من التيار المستمر (ω ≈ 0) يكون التحويل شبه خطي: Ω ≈ ω/T. لكن كلما اقترب ω من π، تُضغط الترددات التناظرية بشدة. ترددات القطع الحادة عند Ω_c في المجال التناظري تقابل تردداً رقمياً مختلفاً ω_c — وهما ليسا متساويين.
تشويه التردد يُشوّه محور التردد — لكن بالنسبة للمرشحات التي تُحدَّد متطلباتها بحواف نطاق التمرير ونطاق الإيقاف لا بالشكل الدقيق، فهذا مقبول. طابع الاستجابة المتموج أو الرتيب داخل كل نطاق يُحفظ. فقط مواضع حواف الانتقال تتحوّل. يُصحّح التحويل المسبق هذه المواضع بالضبط.
طابع النطاق محفوظ — تتغيّر ترددات الحواف فقطالتحويل المسبق للترددات الحرجة
الحل لتشويه التردد هو التحويل المسبق للنموذج التناظري الأولي قبل تطبيق التحويل الثنائي الخطي. إذا أردت أن يكون نقطة −3 ديسيبل في المرشح الرقمي (أو حافة النطاق ذو التموّج المتساوي) عند التردد الرقمي ω_c، فصمّم النموذج التناظري بتردد قطع عند التردد المُحوَّل مسبقاً:
يضمن التحويل المسبق أن تردداً حرجاً واحداً أو اثنين — حافة نطاق التمرير، حافة نطاق الإيقاف، أو نقطة −3 ديسيبل — تقعان في المكان الصحيح بعد أن يضغط التحويل الثنائي الخطي محور التردد. ستظل باقي الترددات مُشوَّهة، لكن المتطلبات ستُلبَّى في المواضع التي تهمّ.
سير عمل التصميم الكامل
يجمع تصميم النموذج التناظري الأولي مع التحويل الثنائي الخطي إجراءً عملياً كاملاً لتصميم المرشح IIR:
النتيجة مرشح IIR رقمي تلبّي استجابته السعوية المتطلبات الأصلية بالضبط عند الترددات الحرجة المُحوَّلة مسبقاً. عملياً، تتولى أدوات التصميم مثل bilinear() في MATLAB وscipy.signal.bilinear() في Python هذا الجبر تلقائياً.
التحويل الثنائي الخطي مقابل الاستجابة النبضية الثابتة
قبل أن يُصبح التحويل الثنائي الخطي معياراً، كان المهندسون يستخدمون طريقة الاستجابة النبضية الثابتة — أخذ عينات من الاستجابة النبضية التناظرية لإنشاء المرشح الرقمي. يُفضَّل التحويل الثنائي الخطي في معظم التطبيقات لأن:
| الخاصية | الاستجابة النبضية الثابتة | التحويل الثنائي الخطي |
|---|---|---|
| الإيهام الترددي | نعم — نطاق الإيقاف يُراكَب | لا — تحويل كامل من s إلى z |
| محور التردد | خطي (لكن مع تراكب) | مُشوَّه (لكن بلا تراكب) |
| أنواع المرشحات المناسبة | تمرير منخفض فقط (يجب أن يكون نطاق الإيقاف مهملاً قرب نايكويست) | أي نوع — LP, HP, BP, BS |
| الاستجابة النبضية | محفوظة (أخذ عينات) | تتغيّر بسبب التشويه |
| الاستخدام الصناعي | نادر (تاريخي) | معياري |
تصميم المرشحات العالية والنطاقية والحاجزة
سير عمل التحويل الثنائي الخطي لا يقتصر على المرشحات ذات التمرير المنخفض. لتصميم مرشحات التمرير العالي، والنطاقي، والحاجب، تُضاف خطوة تحويل تردد إضافية بين النموذج التناظري ذي التمرير المنخفض والتعويض الثنائي الخطي:
butter(N,'high')، cheby1(N,Rp,'bandpass')، وscipy.signal.iirdesign() تُعالج جميع التحويلات والتحويل الثنائي الخطي تلقائياً.الاستقرار ودائرة الوحدة
من أجمل خصائص التحويل الثنائي الخطي ضمان الاستقرار. النصف الأيسر من مستوى s (Re{s} < 0) يتحوّل بالضبط إلى داخل دائرة الوحدة في مستوى z (|z| < 1). المحور التخيلي (jΩ) يتحوّل إلى دائرة الوحدة ذاتها (|z| = 1). النصف الأيمن يتحوّل إلى خارج دائرة الوحدة.
لأن التحويل الثنائي الخطي تحويل أحادي الاتجاه بين مستوى s ومستوى z (باستثناء z = −1 ↔ s = ∞)، فإن النموذج التناظري السببي المستقر — بجميع أقطابه في النصف الأيسر — يُنتج دائماً مرشحاً رقمياً سببياً مستقراً بجميع أقطابه داخل دائرة الوحدة. لا حاجة لفحص الاستقرار بعد التحويل.
النصف الأيسر ← داخل دائرة الوحدة: الاستقرار محفوظ بالبنية- يُعوّض التحويل الثنائي الخطي s = (2/T)(z−1)/(z+1) لتحويل النموذج التناظري H_a(s) إلى مرشح IIR رقمي H(z) بلا تراكب ترددي.
- يضغط التحويل محور التردد التناظري كله (من −∞ إلى +∞) في النطاق الرقمي (−π إلى π) — مُنشئاً تشويه التردد.
- يُصحّح التحويل المسبق هذا التشويه: صمّم النموذج التناظري عند Ω_c = (2/T)·tan(ω_c·T/2) ليقع تردد القطع الرقمي بالضبط حيث تريد.
- الاستقرار مضمون: أقطاب التناظري في النصف الأيسر تتحوّل دائماً إلى أقطاب رقمية داخل دائرة الوحدة.
- تفوّق التحويل الثنائي الخطي على طريقة الاستجابة النبضية الثابتة لأنه يُزيل الإيهام ويعمل مع جميع أنواع المرشحات.
- تحويلات التردد الإضافية (LP→HP, LP→BP, LP→BS) تُؤلَّف مع التحويل الثنائي الخطي لتصميم أي نوع مرشح من نموذج تمرير منخفض.
- عملياً، يُؤتمَت سير العمل كاملاً بدوال
butter()،cheby1()،ellip()في MATLAB وSciPy.