من معادلات الفروق إلى تحويل Z | DSP-101

لماذا نحتاج إلى أداة جبرية؟

في الوحدة الثالثة، درسنا أنظمة LTI من خلال الاستجابة النبضية h[n] والالتواء. لكن حساب y[n] = x[n] * h[n] يتطلب جمع متسلسلة لا نهائية لمرشحات IIR، وتحليل الاستقرار يستدعي فحص كل عينة من h[n] على حدة. أما تصميم المرشح — اختيار المعاملات لتحقيق استجابة ترددية مطلوبة — فهو أصعب بكثير في نطاق الزمن وحده.

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

تعريف تحويل Z

X(z) = \sum_{n=-\infty}^{\infty} x[n]\,z^{-n}

يحوّل تحويل Z المتتالية x[n] إلى دالة X(z) في المتغير المركب z. كل عينة x[n] توزّن بقوة z المرفوعة إلى −n، مما يضغط المتتالية كلها في كيان جبري واحد.

المتغير z مركّب — z = re^jω بالصيغة القطبية. عند r = 1 تُرسم دائرة الوحدة، وX(z) المحسوبة على دائرة الوحدة تساوي تحويل DTFT: X(e^jω). إذن يعمم تحويل Z تحويل DTFT إلى المستوى المركب كله، مما يتيح التحليل خارج محور التردد.

أزواج تحويل Z الشائعة

عدد محدود من المتتاليات يظهر باستمرار لدرجة أن تحويلاتها تُحفظ عن ظهر قلب. اشتقاقها من التعريف تمرين مفيد: عوّض تعريف المتتالية في المجموع وطبّق صيغ المتسلسلات الهندسية.

النبضة الوحدة

\mathcal{Z}\{\delta[n]\} = 1

δ[n] تتحول إلى 1 — صالحة لجميع قيم z. أبسط زوج، يؤكد أن تحويل Z للنبضة الوحدة هو الواحد.

الخطوة الوحدة

\mathcal{Z}\{u[n]\} = \frac{z}{z-1}, \quad |z|>1

u[n] تتحول إلى z/(z−1)، صالحة لـ |z| > 1. القطب عند z = 1 يعكس الطابع غير المتلاشي للخطوة.

المتتالية الأسية

\mathcal{Z}\{a^n u[n]\} = \frac{z}{z-a}, \quad |z|>|a|

المتتالية الأسية aⁿu[n] أحادية الجانب تتحول إلى z/(z−a)، صالحة لـ |z| > |a|. مستقرة عند |a| < 1.

متتالية جيب التمام

\mathcal{Z}\{\cos(\omega_0 n)\,u[n]\} = \frac{z^2 - z\cos\omega_0}{z^2 - 2z\cos\omega_0 + 1}

متتالية جيب التمام تُنتج قطبين مركبين على دائرة الوحدة عند ±ω₀. تُستخدم لتحليل الاستجابة المستقرة الجيبية.

هذه الأزواج، مع خصائص تحويل Z أدناه، تشكل مجموعة أدوات متكاملة. بدلاً من تقييم المجموع اللانهائي في كل مرة، يبني المهندسون X(z) بتفكيك x[n] إلى نسخ مقياسة ومزاحة من متتاليات معروفة وتطبيق الخطية.

الخصائص الرئيسية

تكمن قوة تحويل Z في خصائصه. ثلاث خصائص بالغة الأهمية في تحليل المرشحات.

الخاصية 1

الخطية

Z{ax[n] + by[n]} = aX(z) + bY(z). مبدأ التراكب ساريٌ — المجموع الموزون في الزمن يتحول إلى مجموع موزون من التحويلات.

الخاصية 2

إزاحة الزمن

Z{x[n−k]} = z^−kX(z). تأخير k عينة يضرب في z^−k في نطاق z — أساس معامل التأخير.

الخاصية 3

الالتواء ← ضرب

Z{x[n] * h[n]} = X(z)·H(z). يتحول الالتواء إلى ضرب — السبب الجوهري لأهمية تحويل Z في تحليل أنظمة LTI.

الخاصية 4

التقييس في z

Z{aⁿx[n]} = X(z/a). الضرب بـ aⁿ في الزمن يقيّس مستوى z. يُستخدم في إزاحة التردد وتحويلات المرشح النموذجي.

خاصية إزاحة الزمن

\mathcal{Z}\{x[n-k]\} = z^{-k}\,X(z)

وحدة التأخير z⁻¹ في نطاق z تقابل تأخير عينة واحدة في نطاق الزمن. هذه الخاصية تحول مباشرةً معادلات الفروق إلى معادلات جبرية في z.

تحويل معادلات الفروق

خاصية إزاحة الزمن هي الجسر بين معادلات الفروق ونطاق z. خذ LCCDE عامة بمعاملات {b_k} و{a_k}. طبّق تحويل Z على الطرفين واستخدم الخطية وخاصية إزاحة الزمن على كل حد:

الصورة الجبرية في نطاق Z

Y(z)\sum_{k=0}^{M}a_k z^{-k} = X(z)\sum_{k=0}^{N}b_k z^{-k}

تطبيق تحويل Z على LCCDE يحولها من معادلة زمنية تكرارية إلى معادلة متعددة حدود نسبية — قابلة للحل جبرياً.

نسبة Y(z)/X(z) هي دالة النقل H(z) — تحويل Z للاستجابة النبضية. هي نسبة متعددتَي حدود في z⁻¹ (أو ما يعادله في z)، تحددها معاملات المرشح بالكامل. كل نظام LTI يمكن تمثيله بـ LCCDE له دالة نقل نسبية.

IIR من الدرجة الأولى: نطاق الزمن ← نطاق Z

خذ مرشح IIR البسيط: y[n] = ay[n−1] + x[n]. طبّق تحويل Z على الطرفين:

Y(z) = az⁻¹Y(z) + X(z)

اجمع Y(z): Y(z)(1 − az⁻¹) = X(z)

H(z) = Y(z)/X(z) = 1/(1 − az⁻¹) = z/(z − a)

تحويل Z أحادي الجانب مقابل ثنائي الجانب

التعريف أعلاه يجمع على كل n من −∞ إلى +∞ — هذا هو تحويل Z الثنائي. يعالج المتتاليات السببية وغير السببية معاً لكنه يستلزم تحديد منطقة التقارب (ROC) بعناية لتحديد المعكوس بصورة فريدة.

في الممارسة الهندسية، معظم الإشارات سببية (صفرية لـ n < 0)، لذا يُستخدم تحويل Z أحادي الجانب الذي يجمع من n = 0 إلى +∞. الصيغة أحادية الجانب مفيدة بشكل خاص لحل معادلات الفروق بشروط ابتدائية.

تحويل Z أحادي الجانب

X(z) = \sum_{n=0}^{\infty} x[n]\,z^{-n}

يجمع تحويل Z أحادي الجانب على n ≥ 0 فقط. للمتتاليات السببية النتيجة مطابقة لتحويل Z الثنائي. للمتتاليات غير السببية تختلف الصيغتان.

العلاقة بتحويل لابلاس

تحويل Z هو النظير المتقطع لتحويل لابلاس. تتجلى الصلة حين نتذكر عملية الأخذ بالعينات: الإشارة المستمرة المأخوذة بالعينات x(t) تكتب كـ x[n] = x(nT)، وتحويل لابلاس للإشارة المأخوذة بالعينات يتضمن حدوداً بصيغة e^sT. بتعريف z = e^sT تتضح العلاقة بين المستوى s والمستوى z.

التعيين من المستوى s إلى المستوى z

z = e^{sT}

التعويض z = e^(sT) يعيّن المستوى s على المستوى z. المحور التخيلي jω في المستوى s يتعيّن على دائرة الوحدة في المستوى z. النصف الأيسر من المستوى يتعيّن على داخل دائرة الوحدة.

لهذا التعيين تداعيات عميقة. في المستوى s، يستلزم استقرار النظام المستمر أن تكون أقطابه في النصف الأيسر (الجزء الحقيقي سالب). بعد التعيين، تحل الأقطاب ذات النصف الأيسر داخل دائرة الوحدة في المستوى z. إذن استقرار الزمن المتقطع — جميع الأقطاب داخل دائرة الوحدة — هو النظير الرقمي المباشر لاستقرار الزمن المستمر.

المستوى s مقابل المستوى z في لمحة

منطقة الاستقرار: النصف الأيسر (s) ↔ داخل دائرة الوحدة (z)

محور التردد: المحور التخيلي jω (s) ↔ دائرة الوحدة e^jω (z)

منطقة عدم الاستقرار: النصف الأيمن (s) ↔ خارج دائرة الوحدة (z)

لابلاس ↔ تحويل Z، تماماً كما المستمر ↔ المتقطع

الدرس التالي يستكشف منطقة التقارب (ROC) — مجموعة قيم z التي يتقارب عندها مجموع تحويل Z. تحدد ROC ما إذا كانت الإشارة سببية أم مضاد-سببية أم ثنائية الجانب، وهي ضرورية لعكس تحويل Z بشكل صحيح.

النقاط الرئيسية

يعيّن تحويل Z المتتالية المتقطعة x[n] إلى دالة X(z) في المتغير المركب z، محولاً معادلات الفروق التكرارية إلى معادلات متعددة حدود جبرية.
التعريف: X(z) = Σ x[n]·z⁻ⁿ. عند تقييمه على دائرة الوحدة |z| = 1 يساوي تحويل DTFT.
الأزواج الأساسية: δ[n]↔1، u[n]↔z/(z−1)، aⁿu[n]↔z/(z−a).
الخصائص الأساسية: الخطية، إزاحة الزمن (تأخير k → ضرب في z⁻ᵏ)، الالتواء ← ضرب.
أي LCCDE تصبح دالة نقل نسبية H(z) = Y(z)/X(z) — نسبة متعددتَي حدود في z.
تحويل Z أحادي الجانب (الجمع من n=0) يُستخدم للأنظمة السببية ومسائل الشروط الابتدائية.
التعيين z = e^(sT) يربط تحويل لابلاس بتحويل Z؛ النصف الأيسر من المستوى s يتعيّن داخل دائرة الوحدة في المستوى z.