العملية التي تُعرّف أنظمة LTI
في الدرس السابق أثبتنا أن نظام الزمن المنفصل الخطي الثابت (LTI) يُوصَّف كلياً باستجابته النبضية h[n]. لكن كيف نحسب فعلياً الخرج y[n] لأي إدخال x[n]؟ الجواب هو الالتواء — ولعله أهم عملية رياضية في معالجة الإشارات بأسرها.
الالتواء ليس مجرد معادلة تُحفظ. إنه ينبثق طبيعياً من حقيقتين نعرفهما مسبقاً: يمكن تحليل أي إشارة إلى مجموع موزون من النبضات المُزاحة (خاصية الغربلة)، وأن استجابة نظام LTI لنبضة مُزاحة ومُعامَلة هي نسخة مُزاحة ومُعامَلة من h[n] (بفضل الخطية وثبات الزمن). اجمع الحقيقتين وسيتراءى لك مجموع الالتواء مباشرة.
يستخدم الرمز y[n] = (x * h)[n] النجمة للدلالة على الالتواء. انتبه: في معالجة الإشارات المستمرة تعني النجمة أيضاً المرافق المركب — لكن السياق يُزيل الالتباس دائماً. في DSP للزمن المنفصل، x[n] * h[n] تعني حصراً مجموع الالتواء أعلاه.
من أين يأتي مجموع الالتواء؟
نبدأ من خاصية الغربلة للنبضة الوحدوية δ[n]: يمكن كتابة أي إشارة x[n] كتراكب من نبضات موزونة ومُزاحة:
نطبّق الآن عملية نظام LTI T{·} على طرفي المعادلة. بفضل الخطية، تمر T خلال المجموع ومعاملات التحجيم. بفضل ثبات الزمن، استجابة δ[n − k] هي h[n − k] (نسخة مُزاحة بـ k من h[n]). لذا:
يُوضّح هذا الاشتقاق لماذا الالتواء ليس مجرد خوارزمية ذكية: إنه النتيجة الحتمية لما يعنيه أن يكون النظام LTI. أي نظام خطي ثابت الزمن لا بد أن يُنتج خرجه عبر الالتواء. وإذا لم تستطع كتابة نظام بهذه الصورة، فهو إمّا غير خطي أو متغيّر مع الزمن.
التفسير البياني: اعكس، ازحز، اضرب، اجمع
الشعار الرباعي "اعكس، ازحز، اضرب، اجمع" يصف كيفية تقييم مجموع الالتواء عند أي مؤشر خرج n. تصوّر أنك تُمرّر نسخة معكوسة من h فوق x، وتضرب العناصر المتقابلة عند كل موضع، ثم تجمع النواتج.
كرّر ذلك لكل قيمة n لبناء متتالية الخرج الكاملة. قد يبدو الأمر آلياً، لكن التفسير البصري قوي: التواء إشارة مع مرشّح صندوقي (مستطيل من الآحاد) يحسب متوسطاً محلياً — كل عينة خرج هي متوسط نافذة منزلقة. الالتواء مع نواة مشتقة يُبرز التغيرات المفاجئة. والالتواء مع نواة صدى يُضيف نسخة مؤخَّرة. كل مرشّح محدود الطول هو "نواة التوائه" نفسها.
ليكن x[n] = {1, 2, 3} (غير صفري عند n = 0, 1, 2) و h[n] = {1, 1} (غير صفري عند n = 0, 1). طول الخرج = 3 + 2 − 1 = 4 عينات.
y[0] = x[0]h[0] = 1·1 = 1
y[1] = x[0]h[1] + x[1]h[0] = 1·1 + 2·1 = 3
y[2] = x[1]h[1] + x[2]h[0] = 2·1 + 3·1 = 5
y[3] = x[2]h[1] = 3·1 = 3
y[n] = {1, 3, 5, 3} — المجموع المتحرك للنافذة ذات العنصرين.الخصائص الجبرية للالتواء
يخضع الالتواء لثلاثة قوانين جبرية تُشابه قوانين الضرب العادي وتُبسّط تحليل الأنظمة المتسلسلة والمتوازية:
الخاصية الإبدالية مفاجئة بشكل خاص: تقول إن تطبيق مرشّح h على إشارة x يعطي النتيجة ذاتها تماماً كتطبيق الإشارة x على المرشّح h. التفسيران يبدوان مختلفين تماماً، لكن الرياضيات واحدة. وعملياً، يختار المهندسون أيّ المتغيّرين يعكسونه ويُمرّرونه بحسب أيّهما أرخص حساباً للإشارات التي بين أيديهم.
الخاصية الترابطية هي أساس تصميم المرشّحات: نظامان LTI متسلسلان (خرج الأول يغذّي الثاني) يتصرفان كنظام LTI واحد استجابته النبضية h₁ * h₂. ومعنى ذلك أن ترتيب المرشّحات المتسلسلة لا يؤثر على الخرج النهائي، وأن بإمكانك حساب سلسلة كاملة من المرشّحات مسبقاً ودمجها في نواة واحدة قبل التشغيل.
الالتواء مع النبضة: العملية المحايدة
النبضة الوحدوية δ[n] هي العنصر المحايد للالتواء، تماماً كما أن الضرب في 1 يترك العدد دون تغيير:
بشكل أعم، الالتواء مع نبضة مُزاحة δ[n − n₀] يُنتج تأخيراً نقياً: x[n] * δ[n − n₀] = x[n − n₀]. وهذا يخبرنا أن خط التأخير المثالي نظام LTI استجابته النبضية عينة واحدة عند الزمن n₀. كل ذاكرة تأخير رقمية وكل وحدة صدى وكل مُعالج دوران صوتي مبني على هذا المبدأ.
طول الخرج والتعقيد الحسابي
حين تحتوي x[n] على N عينة غير صفرية وh[n] على M عينة غير صفرية، يحتوي الالتواء y[n] = x * h على N + M − 1 عينة غير صفرية بالضبط. ويسهل أن ترى ذلك: أول عينة خرج غير صفرية تظهر لحظة تتقاطع فيها أول عينة غير صفرية من x مع أول عينة غير صفرية من h، وآخر عينة تظهر لحظة تعبر فيها آخر عينة غير صفرية في كل منهما الأخرى.
يتطلب الحساب المباشر لمجموع الالتواء ما يقارب N · M عملية ضرب-تراكم — O(N·M) إجمالاً. ومع الإشارات الطويلة والمرشّحات الطويلة يصبح هذا مكلفاً. يتيح تحويل فورييه السريع (FFT) حسابه في O((N+M) log(N+M)) عملية عبر مبرهنة الالتواء التي تنص على أن الالتواء في المجال الزمني يساوي الضرب نقطةً بنقطة في المجال الترددي. وهذا الالتواء المبني على FFT هو حصان العمل في معالجة الصوت الحديثة، وتصفية الصور، وأنظمة الاتصالات.
لمقطع صوتي مدته ثانية واحدة عند 44,100 هرتز (N = 44,100 عينة) مُصفَّى باستجابة نبضية غرفة مدتها 0.5 ثانية (M = 22,050 عينة)، يحتاج الالتواء المباشر إلى نحو 44,100 × 22,050 ≈ 972 مليون ضربة. يخفّضها الالتواء بـ FFT إلى نحو (66,150) × log₂(66,150) ≈ 1.1 مليون عملية — أسرع بما يقرب من 1000 مرة. وعملياً، تُطبَّق طريقة الجمع بالتراكب (overlap-add) أو الحفظ بالتراكب (overlap-save) التواءَ FFT على كتل متتابعة من التدفّق.
مباشر: O(N·M) | FFT: O((N+M) log(N+M))FIR وIIR: عائلتان من أنظمة LTI
يُقسّم الالتواء طبيعياً أنظمة LTI إلى عائلتين بحسب مدة الاستجابة النبضية:
أنظمة FIR (الاستجابة النبضية المنتهية) لها استجابة نبضية h[n] غير صفرية لعدد منتهٍ فقط من العينات. مجموع التواويها منتهٍ — سهل الحساب المباشر ومضمون الاستقرار. مرشّحات المتوسط المتحرك والمرشّحات المتطابقة ومعظم المرشّحات المُصمَّمة بنوافذ هي FIR. والخرج عند أي لحظة n يعتمد على نافذة منتهية من الإدخالات السابقة.
أنظمة IIR (الاستجابة النبضية اللانهائية) لها استجابة نبضية تمتد إلى ما لا نهاية. تُنفَّذ عادةً عبر معادلات الفروق لا الالتواء المباشر، لأن جمع عدد لانهائي من الحدود غير عملي. خُذ مرشّحاً عودياً من الرتبة الأولى: y[n] = αy[n−1] + x[n]. استجابته النبضية h[n] = αnu[n] — تتضاءل أُسّياً لكنها لا تنتهي فعلياً أبداً. يمكن لمرشّحات IIR تحقيق انتقائية ترددية حادة بعدد معاملات أقل بكثير من مرشّحات FIR المكافئة، لكنها تستوجب تحليل استقرار دقيقاً (|α| < 1 في المثال أعلاه).
- الالتواء y[n] = Σ x[k]·h[n−k] هو معادلة الخرج الأساسية لأي نظام LTI؛ ويُشتق مباشرة من الخطية وثبات الزمن وخاصية الغربلة.
- الطريقة البيانية الرباعية — اعكس h، ازحزه بمقدار n، اضربه مع x نقطةً بنقطة، اجمع — تُعطي تفسيراً حدسياً لتقييم الالتواء عند كل مؤشر خرج.
- الالتواء إبدالي (x*h = h*x)، ترابطي ((x*h₁)*h₂ = x*(h₁*h₂))، وتوزيعي على الجمع.
- النبضة الوحدوية δ[n] هي المحايد: x[n] * δ[n] = x[n]؛ والنبضة المُزاحة δ[n−n₀] تمثّل تأخيراً نقياً.
- التواء إشارة طولها N مع نواة طولها M يُنتج خرجاً طوله N+M−1؛ التعقيد المباشر O(N·M)، وبـ FFT ينخفض إلى O((N+M)log(N+M)).
- أنظمة FIR لها استجابة نبضية منتهية ومضمونة الاستقرار؛ وأنظمة IIR لها استجابة نبضية لانهائية تُنفَّذ عبر العودية وتستلزم التحقق من الاستقرار.