من نطاق الزمن إلى نطاق التردد | DSP-101

كل إشارة هي مجموع من الجيبيات

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

التفكير بالترددات يؤتي ثماره في معالجة الإشارات الرقمية. المرشحات تُصمَّم بشكل طبيعي بتحديد الترددات التي تمر وتلك التي تُحجب. الضوضاء التي تتركز في نطاق تردد معين يسهل قمعها حين تظهر كقمة طيفية. الإشارات الصوتية والموسيقية وإشارات الراديو لها بصمات تردد مميزة تكشف بنية خفية لا تراها في نطاق الزمن. نطاق التردد ليس مجرد وسيلة مريحة — بل هو في الغالب الموطن الأكثر ملاءمة لحل مشكلاتنا.

تحويل فورييه للزمن المتقطع (DTFT)

يُحوّل DTFT متتالية الزمن المتقطع x[n] إلى دالة مستمرة في التردد X(e^jω). على عكس تحويل Z الذي يستخدم متغيرًا مركبًا عامًا z، يُقيّم DTFT تحويل Z على دائرة الوحدة z = e^jω. النتيجة دالة مركبة القيمة لتردد زاوي معيّر ω (راديان لكل عينة)، وهي دورية بدورة 2π.

تعريف DTFT

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

يجمع DTFT المتتالية x[n] المُرجَّحة بأسس مركبة e^{−jωn}. النتيجة X(e^{jω}) دالة مستمرة دورية في ω بدورة 2π. يتطلب التقارب أن تكون المتتالية قابلة للجمع المطلق.

يستعيد DTFT العكسي x[n] بتكامل X(e^jω) على فترة واحدة:

DTFT العكسي

x[n] = \frac{1}{2\pi}\int_{-\pi}^{\pi} X(e^{j\omega})\,e^{j\omega n}\,d\omega

تكامل X(e^{jω}) e^{jωn} على [−π, π] مقسومًا على 2π يسترجع المتتالية الأصلية x[n] بدقة. هذا يؤكد أن DTFT تمثيل كامل وقابل للعكس.

التردد المعيَّر

في الزمن المتقطع، يُقاس التردد ω بـراديان لكل عينة، ويتراوح بين −π و π. للتحويل إلى تردد فيزيائي بالهرتز، اضرب في f_s/(2π). تردد نايكويست f_s/2 يقابل ω = π.

ω = 2π f / f_s — فـ ω = π يقابل f = f_s/2 (نايكويست)

خصائص DTFT الأساسية

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

الخطية

مبدأ التراكب

DTFT{ax[n] + by[n]} = aX(e^jω) + bY(e^jω). فكّك الإشارات المعقدة إلى مكونات أبسط، حوّل كلًا منها، ثم اجمع النتائج.

إزاحة زمنية

إزاحة طور

DTFT{x[n − n₀]} = e^−jωn₀ X(e^jω). تأخير الإشارة بمقدار n₀ عينة يضرب طيفها بعامل طور خطي — المقدار لا يتغير.

الالتواء

الضرب

DTFT{x[n] ∗ h[n]} = X(e^jω) H(e^jω). الالتواء في الزمن يساوي الضرب النقطي في التردد — أساس الترشيح الخطي.

التناظر

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

للإشارة الحقيقية x[n]، |X(e^jω)| زوجي و ∠X(e^jω) فردي. هذا التناظر المترافق يعني أن الطيف الكامل محدد بقيمه على [0, π].

ثنائية الالتواء-الضرب هي على الأرجح أهم خاصية. تعني أن الترشيح — وهو التواء في نطاق الزمن — يصبح ضربًا نقطيًا في نطاق التردد. لتصميم مرشح تمرير منخفض، تحدد ببساطة H(e^jω) = 1 لـ |ω| < ω_c وH(e^jω) = 0 غير ذلك، ثم تجد h[n] المقابل عبر DTFT العكسي.

أزواج DTFT الشائعة

عدد من أزواج DTFT يتكرر باستمرار في معالجة الإشارات الرقمية. حفظ هذه الأزواج — أو القدرة على اشتقاقها بسرعة — يوفر قدرًا كبيرًا من الحسابات.

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

\delta[n] \xrightarrow{\text{DTFT}} 1 \quad (\text{flat spectrum for all }\omega)

النبضة الواحدة δ[n] لها طيف مسطح (أبيض) — طاقة متساوية عند جميع الترددات. لهذا تُعدّ الإشارة الاختبارية المثالية لأنظمة LTI.

الأسية السببية

a^n u[n] \xrightarrow{\text{DTFT}} \frac{1}{1 - a\,e^{-j\omega}}, \quad |a| < 1

الأسية السببية a^n u[n] (مع |a| < 1 للتقارب) تتحول إلى استجابة تردد نسبية. يبلغ المقدار ذروته قرب ω = 0 لقيم a حقيقية موجبة وقرب ω = π لقيم سالبة.

النافذة المستطيلة (المرشح المثالي)

H(e^{j\omega}) = \begin{cases} 1 & |\omega| \leq \omega_c \\ 0 & \omega_c < |\omega| \leq \pi \end{cases} \xrightarrow{\text{IDTFT}} h[n] = \frac{\omega_c}{\pi}\operatorname{sinc}\!\left(\frac{\omega_c n}{\pi}\right)

المرشح المثالي لتمرير المنخفض له طيف مستطيل (1 لـ |ω| ≤ ωc، و0 غير ذلك). DTFT العكسي له هو متتالية sinc المستمرة لا نهائيًا، ولهذا يجب اقتطاعها وتطبيق نافذة عليها في التطبيقات العملية.

القيد: التردد المستمر

DTFT رائع رياضيًا لكنه غير مريح حسابيًا. مخرجه X(e^jω) دالة مستمرة في ω — يأخذ قيمًا لا نهاية لها. الحاسوب الرقمي لا يستطيع تخزين أو معالجة سوى مجموعة منتهية من الأرقام. لا يمكن تمثيل أو حساب DTFT الكامل على الحاسوب دون تقريب أو أخذ عينات.

القضية الثانية هي الكفاءة. حتى تقييم X(e^jω) عند تردد واحد ω₀ يتطلب جمعًا على جميع عينات x[n] البالغة N، وهو O(N) لكل نقطة تردد. التقييم عند K نقطة تردد يكلف O(NK) — كبير جدًا للإشارات الطويلة.

لماذا لا يكفي DTFT وحده

يحوّل DTFT متتالية زمن متقطع إلى دالة تردد مستمرة. تحتاج الأنظمة الرقمية إلى تمثيل تردد متقطع — مجموعة منتهية من الأرقام — ليكون قابلًا للحساب. الحل هو أخذ عينات من DTFT عند N نقطة متساوية التباعد.

أخذ عينات من DTFT: تحويل DFT

إذا قيّمنا DTFT عند N تردد متساوي التباعد ω_k = 2πk/N لـ k = 0, 1, …, N−1، نحصل على تحويل فورييه المتقطع (DFT):

تعريف DFT

X[k] = \sum_{n=0}^{N-1} x[n]\,e^{-j2\pi kn/N}, \quad k = 0, 1, \ldots, N-1

X[k] هو DFT للمتتالية x[n] ذات N نقطة. كل bin k يقابل التردد ωk = 2πk/N راديان/عينة، أو fk = k·fs/N هرتز. DFT متقطع في الزمن والتردد معًا — قابل للحساب الكامل.

DFT هو DTFT المأخوذة عيناته عند N نقطة متساوية التباعد على دائرة الوحدة. الأهم من ذلك أنه تحويل منتهٍ وقابل للعكس: N عدد مركب يدخل وN عدد مركب يخرج، والمتتالية الأصلية تُستعاد بدقة عبر DFT العكسي. هذا ما جعله العمود الفقري لتحليل الطيف الحسابي.

غير أن أخذ عينات من DTFT في التردد له ثمن. تمامًا كما أن أخذ العينات في الزمن يُنشئ دورية في التردد (الترادف)، فأخذ العينات في التردد يُنشئ دورية في الزمن — DFT يعامل x[n] ضمنيًا كدورة واحدة من إشارة دورية. هذه الطبيعة الدائرية (الدورية) لـ DFT لها تداعيات مهمة للنوافذ والحشو بالأصفار والالتواء الدائري، وكلها تُستكشف في الدروس التالية.

الدرس التالي يتعمق في صيغة DFT، معنى كل bin تردد، دقة التردد، وDFT باعتباره عملية مصفوفية — بناء الأدوات الحسابية لتحليل الطيف.

أبرز ما تعلمته

كل إشارة هي مجموع موزون من الجيبيات؛ نطاق التردد يكشف هذا التحليل مما يجعل تصميم المرشحات وتحليل الضوضاء طبيعيين وبديهيين.
يحوّل DTFT متتالية الزمن المتقطع x[n] إلى دالة مستمرة دورية بدورة 2π هي X(e^jω) للتردد الزاوي المعيّر ω (راديان/عينة).
التردد المعيّر ω = π يقابل تردد نايكويست f_s/2؛ العلاقة هي ω = 2πf/f_s.
الالتواء في الزمن يساوي الضرب في التردد — أساس ترشيح أنظمة LTI في نطاق التردد.
DTFT مستمر في التردد وبالتالي لا يُحسب مباشرة؛ تحتاج الآلة الرقمية إلى تمثيل تردد منتهٍ ومتقطع.
أخذ عينات من DTFT عند N نقطة متساوية التباعد ينتج DFT: X[k] = Σ x[n] e^−j2πkn/N، وهو قابل للحساب والعكس الكامل.
يعامل DFT x[n] ضمنيًا كإشارة دورية بدورة N — خاصية محورية تقود اعتبارات النوافذ والحشو بالأصفار.