صلة متسلسلات فورييه

الأساس المثلثي

في الدرس الثاني أثبتنا أن L²([−π, π]) هو فضاء هيلبرت — فضاء حاصل ضرب داخلي مكتمل. الآن نتساءل: هل يمتلك أساساً متعامداً طبيعياً؟ الجواب نعم، وهذا الأساس هو المنظومة المثلثية.

نعرّف الدوال: φ₀(t) = 1/√(2π) وللـ n ≥ 1: φₙᶜ(t) = cos(nt)/√π وφₙˢ(t) = sin(nt)/√π. هذه الدوال متعامدة اثنتين اثنتين بموجب حاصل الضرب الداخلي في L²، ولكل منها نظيرة وحدة. تُشكّل معاً مجموعة متعامدة طبيعية لا تنتهي وهي أيضاً مكتملة — لا توجد دالة غير صفرية تعامد جميعها.

تعامد المنظومة المثلثية

\int_{-\pi}^{\pi} \cos(mt)\cos(nt)\,dt = \pi\,\delta_{mn}, \quad \int_{-\pi}^{\pi}\sin(mt)\sin(nt)\,dt = \pi\,\delta_{mn}, \quad \int_{-\pi}^{\pi}\cos(mt)\sin(nt)\,dt = 0

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

معاملات فورييه كحواصل ضرب داخلية

لأن المنظومة المثلثية مجموعة متعامدة طبيعية مكتملة في L²([−π, π])، يمكن تطوير أي f ∈ L² على شكل متسلسلة لا نهائية. معاملات التطوير هي بالضبط حواصل الضرب الداخلية لـ f مع كل دالة أساس:

معاملات فورييه

a_n = \frac{1}{\pi}\int_{-\pi}^{\pi} f(t)\cos(nt)\,dt, \quad b_n = \frac{1}{\pi}\int_{-\pi}^{\pi} f(t)\sin(nt)\,dt

معاملات فورييه a₀ وaₙ وbₙ هي حواصل ضرب داخلية في L². تقيس مدى تمثيل كل تردد أساسي في f. هذا مشابه تماماً لـ v_i = ⟨v, eᵢ⟩ في ℝⁿ — استخراج الإحداثيات بالإسقاط على متجهات الأساس.

متسلسلة فورييه هي التطوير اللانهائي الناتج:

متسلسلة فورييه

f(t) = \frac{a_0}{2} + \sum_{n=1}^{\infty}\bigl[a_n\cos(nt) + b_n\sin(nt)\bigr]

تتقارب المتسلسلة بمعنى نظيرة L²: ‖f − Sₙf‖ → 0 عندما N → ∞. هذا لا يضمن التقارب النقطي في كل مكان، لكن للدوال المتعددة الملاسة تتقارب المتسلسلة نقطياً عند كل نقطة استمرار.

المجاميع الجزئية كإسقاطات متعامدة

المجموع الجزئي من الدرجة N، وهو Sₙf، ليس مجرد اقتطاع — بل هو الإسقاط المتعامد لـ f على الفضاء الجزئي المنتهي الأبعاد Vₙ = span{1, cos t, sin t, …, cos(Nt), sin(Nt)}.

بموجب مبرهنة الإسقاط، Sₙf هو العنصر الوحيد في Vₙ الأقرب إلى f بمعنى نظيرة L². لا يمكن لأي تركيبة خطية أخرى من التوافقيات الـ N+1 الأولى أن تحقق نتيجة أفضل. هذا التفسير الهندسي يشرح لماذا تُعدّ متسلسلات فورييه الأداة الطبيعية لتقريب الدوال ذات البنية الدورية.

خاصية أفضل تقريب

S_N f = \underset{g\,\in\,V_N}{\arg\min}\,\|f - g\|_{L^2}, \qquad \langle f - S_N f,\, g\rangle = 0 \;\forall\, g \in V_N

المجموع الجزئي من الدرجة N هو المصغّر للخطأ في L² على جميع التركيبات الخطية من دوال الأساس الـ N+1 الأولى. الخطأ f − Sₙf متعامد مع Vₙ — هذا هو شرط التعامد الذي يُميّز الإسقاط.

البقية متعامدة مع الفضاء الجزئي

خاصية جوهرية للإسقاط المتعامد: الخطأ f − Sₙf يحقق ⟨f − Sₙf, g⟩ = 0 لكل g ∈ Vₙ. هذا يعني أن الخطأ لا يحمل أي مركّب في الاتجاهات التي قاربناها مسبقاً. كل توافقية إضافية ندرجها تلتقط قدراً من الخطأ المتبقي، وذلك القدر الملتقط يكون دوماً متعامداً مع كل ما أُدرج من قبل.

متسلسلة فورييه المركّبة

العمل بالأسس الأسية المركّبة eⁱⁿᵗ يُبسّط الجبر بشكل كبير. نعرّف ψₙ(t) = eⁱⁿᵗ/√(2π) لـ n ∈ ℤ. تُشكّل هذه مجموعة متعامدة طبيعية مكتملة في L² المركّب ([−π, π]) المزوّد بحاصل الضرب الداخلي ⟨f, g⟩ = ∫₋π^π f(t)g̅(t) dt.

معاملات فورييه المركّبة

c_n = \frac{1}{2\pi}\int_{-\pi}^{\pi} f(t)\,e^{-int}\,dt, \qquad f(t) = \sum_{n=-\infty}^{\infty} c_n\,e^{int}

المعامل المركّب cₙ هو حاصل الضرب الداخلي لـ f مع دالة الأساس الأسية المركّبة رقم n. تربطهما بالشكل الحقيقي العلاقتان: aₙ = 2Re(cₙ) وbₙ = −2Im(cₙ). يدمج الترميز المركّب أوضاع الجيب وجيب التمام في صيغة واحدة مدرّجة على جميع الأعداد الصحيحة.

التقارب وظاهرة غيبس

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

التقارب النقطي

تنص مبرهنة ديريخليه: إذا كانت f متعددة الملاسة على [−π, π]، فتتقارب متسلسلة فورييه عند كل نقطة استمرار إلى f(t)، وعند كل انقطاع قفزي x₀ إلى المتوسط (f(x₀⁺) + f(x₀⁻))/2. هذه إفادة نقطية، أقوى من تقارب L² لكنها تتطلب ملاسة أكبر.

ظاهرة غيبس

قرب الانقطاع القفزي، تتجاوز المجاميع الجزئية القيمة الحقيقية بنحو 9% من ارتفاع القفزة، بصرف النظر عن عدد الحدود المضمّنة. هذا التجاوز — ظاهرة غيبس — قيد جوهري: Sₙf هو أفضل تقريب بمعنى L² من Vₙ، لكن "الأفضل في L²" لا يعني "صغير في كل مكان".

الأثر الهندسي

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

تحويل فورييه المنفصل

عملياً، تُؤخذ عينات الإشارات عند N نقطة زمنية منفصلة t_k = 2πk/N للـ k = 0, 1, …, N−1. يمتلك فضاء المتجهات المركّبة ℂᴺ تناظره المنفصل الخاص للمنظومة المثلثية: متجهات أساس DFT w_n = (1, ωⁿ, ω²ⁿ, …, ω^((N−1)n))/√N حيث ω = e^(i2π/N).

DFT كضرب مصفوفة بمتجه

X[k] = \sum_{n=0}^{N-1} x[n]\,\omega^{kn}, \quad \omega = e^{-i2\pi/N}, \qquad \mathbf{X} = F\mathbf{x},\; F F^* = I

مصفوفة DFT F يكون عنصرها F_{kn} = ω^(kn)/√N. هي مصفوفة موحّدة: F*F = I. حساب DFT إذن مكافئ للتحول إلى أساس متعامد طبيعي في ℂᴺ. تحسب خوارزمية FFT هذا الضرب في O(N log N) بدلاً من O(N²) باستغلال البنية العودية لمصفوفة DFT.

تطبيقات مُعاد رؤيتها عبر الجبر الخطي

التحليل الطيفي

حساب كثافة الطيف الطاقوي لإشارة هو حساب ‖Px‖² لكل إسقاط Px على الفضاء الجزئي المُمتدّ بمتجه أساس DFT واحد. الطيف الطاقوي هو مربّع سعات الإحداثيات في أساس DFT. الجبر الخطي يوضّح لماذا تكون أطياف الطاقة غير سالبة (لأنها مربّعات نظائر) ولماذا تتجمع عبر الفضاءات الجزئية المتعامدة (بارسيفال).

الالتفاف الدائري كضرب قطري

في أساس DFT، يصبح الالتفاف الدائري ضرباً نقطياً: F(x ⊛ h) = √N · (Fx) ⊙ (Fh)، حيث ⊙ هو الضرب مكوّناً بمكوّن. والجذر √N هو ثمن التطبيع الأحادي بـ1/√N الذي اعتمدناه للمصفوفة F؛ فمع DFT غير المطبَّع يختفي هذا العامل وتصبح المتطابقة F(x ⊛ h) = (Fx) ⊙ (Fh)، وهي الصيغة التي ستراها مذكورة عادةً. وفي الحالتين الحقيقة البنيوية واحدة، وهي تصحّ لأن متجهات أساس DFT هي متجهات ذاتية لكل عملية التفاف دائري. تطبيق مرشّح إذن يتضمن: (١) الإسقاط على أساس DFT، (٢) تحجيم كل مركّب بالاستجابة الترددية للمرشّح، (٣) الإسقاط العكسي. هذا هو التقطير القطري — أهم عملية في الجبر الخطي، مُطبَّقة الآن على معالجة الإشارات.

مبرهنة الالتفاف

F(\mathbf{x} \circledast \mathbf{h}) = \sqrt{N}\,(F\mathbf{x}) \odot (F\mathbf{h}) \quad\Longleftrightarrow\quad Y[k] = H[k]\,X[k]

الالتفاف الدائري في المجال الزمني يساوي الضرب النقطي (هاداماد) في المجال الترددي. ويظهر الجذر √N لأننا طبّعنا F بـ1/√N لتصبح أحادية؛ أما الصيغة المكافئة Y[k] = H[k]X[k] فمكتوبة بدلالة تحويلات DFT غير المطبَّعة، وفيها يذوب هذا العامل. هذه هي علاقة القيمة الذاتية–المتجه الذاتي: الأسس الأسية المركّبة eⁱⁿᵗ هي دوال ذاتية لأي عملية خطية ثابتة زمنياً، والاستجابة الترددية H(n) للمرشّح هي القيمة الذاتية المقابلة. الترشيح هو التقطير القطري.

OFDM في معادلة واحدة

يُعدّل OFDM رموز البيانات المستقلة X[0], X[1], …, X[N−1] على متجهات أساس DFT، ثم يُرسل إشارة المجال الزمني x = F⁻¹X. في جهة الاستقبال، يُطبّق القناة التفافاً دائرياً مع h. الإشارة المستقبَلة y = h * x + ضوضاء تحقق Fy = (Fh) ⊙ (Fx) + F(ضوضاء). بعد إزالة البادئة الدورية، يُطبّق المستقبِل F للحصول على Y[k] = H[k] X[k] + W[k] — N قناة سلمية مستقلة، لكل منها كسب مركّب H[k]. هذه المعادلة أحادية الصنبور ممكنة بالضبط لأن DFT يُقطّر الالتفاف الدائري قطرياً.

الترشيح المطابق كحاصل ضرب داخلي

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

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

متسلسلة فورييه هي إسقاط متعامد: كل معامل aₙ أو bₙ هو حاصل ضرب داخلي في L²، والمجموع الجزئي من الدرجة N هو أقرب عنصر من الفضاء الجزئي ذي N توافقية إلى f. الأسس الأسية المركّبة تُشكّل أساساً متعامداً طبيعياً مكتملاً؛ DFT هو النسخة المنتهية الأبعاد مع مصفوفة تغيير أساس موحّدة. الالتفاف الدائري يُقطَّر بواسطة DFT — الترشيح هو تغيير الأساس والتحجيم بالقيم الذاتية ثم العودة. ظاهرة غيبس تُظهر أن الأمثلية بمعنى L² لا تعني الدقة النقطية عند الانقطاعات. هذا التوحيد — تحليل فورييه كجبر خطي — هو الجسر المفاهيمي بين نظرية حاصل الضرب الداخلي المجردة ومعالجة الإشارات العملية.