نجمّع كل حاجة مع بعض
طول الوحدة العاشرة، بنينا مجموعة أدوات نظرية متكاملة: الـ autocorrelation والـ cross-correlation لقياس التشابه بين الـ signals، وإطار الـ Neyman-Pearson detection لأخد قرارات binary في وجود الـ noise، والـ matched filter (المرشح المتوافق) لاستخلاص أقصى SNR. ودلوقتي هنستخدم الأدوات دي في تلاتة من أهم مجالات التطبيق في الهندسة الكهربائية: الرادار، والسونار، والاتصالات الرقمية.
كل مجال بيطرح نفس المشكلة الأساسية — كشف signal شكلها معروف مدفونة في الـ noise وتوصيفها — بس بقنوات فيزيائية مختلفة، ومقاييس زمن وfrequency مختلفة، ومقاييس أداء مختلفة. وإنك تشوف إزاي نفس الإطار الرياضي بيشتغل في المجالات دي كلها بيعمّق الحدس وبيبيّن ليه النظرية اللي طوّرناها في الوحدة دي مش تجريد أكاديمي، دي أساس هندسي ثابت.
الأنظمة التلاتة كلها بتبعت أو بتستقبل waveform معروفة، وبتطبّق matched filter (أو correlator)، وبتقارن الخرج بـthreshold، وبعدين بتقرا التأخير الزمني أو الـ frequency shift أو هوية الـ symbol. الخوارزمية الأساسية واحدة؛ الطبقة الفيزيائية بس هي اللي بتختلف.
الرادار: المسافة والسرعة والـ resolution
الرادار (Radio Detection And Ranging) بيبعت نبضة راديو وبيسمع الصدى اللي راجع من الهدف. وزمن الرحلة رايح جاي τ بيحدّد المسافة R = cτ/2 حيث c سرعة الضوء. وكشف الصدى هو بالظبط binary hypothesis test: H0 (مفيش هدف، noise بس) قصاد H1 (فيه هدف، signal + noise). والـ statistic المثالية لاختبار Neyman-Pearson هي الـ cross-correlation بين الـ waveform المستقبَلة والمبعوتة — يعني خرج الـ matched filter.
الـ matched filter بيعمل وظيفتين في نفس الوقت في الرادار. الأولى، بيكبّر الـ peak SNR بتاع الصدى، فبيخلّي الـ detection ممكن على أبعد مسافة ممكنة. التانية، لما الـ waveform المبعوتة تبقى نبضة chirp (linear frequency modulation)، الـ matched filter بيضغط الصدى لنبضة ضيقة عرضها تقريبًا 1/B حيث B هو الـ bandwidth اللي بيتكنس. والـ pulse compression ده بيفصل الـ range resolution (اللي B بيتحكّم فيها) عن الطاقة المبعوتة (اللي مدة النبضة T بتتحكّم فيها)، فبيخلّي الرادار الحديث يجمع بين مدى detection بعيد وrange resolution عالية في نفس الوقت.
الهدف المتحرّك بيعمل Doppler shift fd = 2v/λ، حيث v السرعة الـ radial وλ الـ wavelength. والـ shift ده بيبعّد spectrum الصدى عن spectrum الـ signal المبعوتة، فبيقلّل خرج الـ matched filter. وعشان تقيس السرعة، الرادارات الحديثة بتستخدم bank من الـ matched filters، كل واحد متظبط على Doppler shift مختلف، أو اللي بيكافئه: correlation ثنائي الأبعاد على التأخير والـ Doppler — الـ ambiguity function، وهي الـ cross-correlation بين الـ signal المبعوتة ونسخة منها متزحزحة بالـ Doppler.
السونار: الصوتيات تحت المية
السونار (Sound Navigation And Ranging) هو النسخة الصوتية من الرادار. بدل الموجات الكهرومغناطيسية بيستخدم موجات صوت بتمشي في المية بسرعة حوالي 1500 م/ث — أبطأ من الضوء بحوالي 200,000 مرة، بس الـ attenuation بتاعها أقل بكتير عند الـ frequencies المستخدمة. والرياضيات هي هي: ابعت waveform معروفة، طبّق matched filter على الصدى المستقبَل، وقارن الخرج بـ threshold.
السونار الـ active بيبعت نبضة وبيسمع الرجوع. أما السونار الـ passive فبيسمع بس، وبيستخدم الـ cross-correlation بين sensorين عشان يحدّد اتجاه مصدر الصوت ومسافته. وفي السونار الـ passive، "الـ signal المعروفة" بيحل محلها الـ cross-correlation بين hydrophoneين متباعدين في المكان: قمة الـ correlation عند التأخير τ بتشاور إن نفس الـ wavefront وصلت للـ sensorين بفرق زمن τ، ومن الفرق ده بنطلّع زاوية الوصول.
أكبر تحدي في السونار هو الـ multipath propagation. الصوت بيرتد من سطح المحيط، ومن قاع البحر، ومن الطبقات الحرارية، فبيوصل للـ receiver كذا نسخة من الـ signal المبعوتة، كل واحدة متأخّرة ومتوهّنة بدرجة مختلفة. وخرج الـ matched filter بيبيّن قمم كتير — قمة لكل مسار — والـ receiver لازم يميّز رجوع المسار المباشر عشان ميغلطش في تقدير المسافة. وده شبه جدًا مسألة الـ channel equalization في الاتصالات الرقمية.
hydrophoneين بينهم مسافة d بيستقبلوا نفس الـ wavefront في أوقات مختلفة شوية. وقمة الـ cross-correlation عند التأخير τ بتدّي الـ time difference of arrival (TDOA): τ = d sinθ/c، حيث θ زاوية الـ azimuth وc سرعة الصوت في المية.
τTDOA = d sinθ / c — الـ azimuth من فرق التأخيرالاتصالات الرقمية: الـ demodulation وتقدير القناة
في نظام الاتصالات الرقمية، الـ transmitter بيبعت سلسلة symbols، كل واحد بيمثّله pulse shape g[n]. والقناة بتضيف noise، ويمكن انعكاسات multipath. والـ receiver لازم يحدّد أنهي symbol اتبعت في كل لحظة، ويخلّي احتمال القرار الغلط أصغر ما يمكن. وده كمان hypothesis test: H0 (الـ bit 0، النبضة -g[n]) قصاد H1 (الـ bit 1، النبضة +g[n]).
الـ matched filter للمسألة دي الـ impulse response بتاعه h[n] = g*[N−1−n]، مراية الـ pulse shape في الزمن مع الـ conjugate. وإنك تاخد sample من الخرج في آخر كل فترة symbol بيدّي الـ sufficient statistic للقرار المثالي. ولـ signaling antipodal binary في AWGN، احتمال الغلط بيبقى:
لما القناة تدخّل multipath (يعني intersymbol interference، ISI)، الـ matched filter لوحده مبيكفّيش. الـ signal المستقبَلة هي convolution السلسلة المبعوتة مع الـ impulse response بتاعة القناة، والنبضات المتراكبة بتخلّي كل symbol يتداخل مع اللي جنبه. وعشان كده الـ receivers الحديثة بتجمع بين الـ matched filtering وequalizer — زي الـ minimum mean square error (MMSE) equalizer أو Viterbi decoder على trellis — عشان تشيل الـ ISI قبل ما تاخد القرارات.
تقدير القناة بيستخدم الـ cross-correlation على طول. الـ transmitter بيبعت pilot sequence p[n] الـ autocorrelation بتاعها شبه impulse (زي PN sequence). والـ receiver بيحسب الـ cross-correlation بين الـ signal المستقبَلة والـ pilot sequence المعروفة عشان يقدّر الـ impulse response بتاعة القناة h[n]: قمة الـ correlation عند التأخير k بتدّي معامل القناة عند الـ tap k. ودي نفس طريقة الـ matched filter، بس متطبّقة على مسألة الـ channel identification بدل مسألة الـ detection.
الـ GPS: correlation في بُعدين
الـ acquisition بتاع GPS حالة بارزة لازم فيها تعمل correlation على مجهولين في نفس الوقت: الـ code phase (أنهي إزاحة chip نستخدم) والـ Doppler shift بتاع الـ carrier (لأن القمر بيتحرّك بالنسبة للـ receiver). والـ receiver بيحتفظ بشبكة بحث ثنائية الأبعاد وبيحسب الـ correlation بين الـ signal الجاية وكل تركيبة من الـ code phase وفرضية الـ Doppler.
الـ receivers الحديثة بتاعة GPS بتستخدم الـ FFT عشان تسرّع البحث ده. لفرضية Doppler ثابتة، الـ correlation على كل الـ code phases هو circular cross-correlation وممكن يتحسب في O(N log N) عملية بالـ FFT، بدل O(N²) للبحث المباشر. وأول ما الخانة الصح (الـ code phase، الـ Doppler) تتحدّد بقمة بتعدّي الـ detection threshold، الـ receiver بيدخل وضع الـ tracking، وبيستعين بـ phase-locked loops (PLL) وdelay-locked loops (DLL) عشان يحافظ على الـ synchronization.
الـ acquisition بتاع GPS بيحسب correlation الـ PRN code اللي فيه 1023 chip على كل الـ code phases الممكنة. البحث المباشر: O(1023²) عملية ضرب لكل Doppler bin. بالـ FFT: O(1023 log 1023) لكل Doppler bin — تسريع حوالي 100×. والـ chips الحديثة بتشغّل مئات الـ Doppler bins بالتوازي، فبتلقط الـ signal في أقل من ثانية.
الـ ultrasound: التصوير بالـ correlation
الـ ultrasound الطبي قريب جدًا من السونار. الـ transducer بيطلّع نبضة صوت قصيرة (عادةً 2–15 ميجاهرتز) وبيسجّل الأصداء اللي راجعة من حدود الأنسجة. والـ time of flight بيدّي عمق كل عاكس، والـ amplitude بتاع خرج الـ matched filter بيدّي معامل الانعكاس. ومسح الـ beam إلكترونيًا بيطلّع صورة مقطعية ثنائية الأبعاد (B-mode) في الـ real time.
الـ Doppler ultrasound بيقيس سرعة سريان الدم بتتبّع الـ Doppler shift بتاع الأصداء من كرات الدم الحمرا المتحرّكة. ودي بالظبط مسألة Doppler الرادار، بس متصغّرة لـ frequencies صوتية ومسافات بالسنتيمترات. ونفس مقايضة الـ ambiguity function بتنطبق هنا: النبضة الطويلة بتحدّد السرعة بدقة بس بتدّي axial resolution (في العمق) ضعيفة؛ والنبضة القصيرة بتدّي resolution حادة في العمق بس resolution ضعيفة في السرعة.
محرّك الـ correlation: هيكل مشترك
لو بصّينا ورا، كل تطبيقات الدرس ده بتشترك في نفس سلسلة المعالجة بالظبط:
1. ابعت (أو راقب) waveform معروفة s[n].
2. استقبل r[n] = s[n−τ] * h[n] + w[n] (متأخّرة، ومعمول لها convolution مع القناة، مع noise).
3. احسب الـ correlation: Rrs[k] = Σ r[n] s*[n−k].
4. الـ threshold: أعلن detection لو max|Rrs[k]| > γ.
5. قدّر: اقرا التأخير أو الـ Doppler أو الـ symbol من مكان القمة.
الهيكل ده هو الـ receiver بتاع الـ matched filter/الـ correlator، وهو مثالي بمعنى Neyman-Pearson وبمعنى أقل BER. ودرجات الحرية الوحيدة اللي فاضلة هي اختيار الـ waveform s[n] — وهي اللي بترسم الـ ambiguity function والـ side lobes بتاعة الـ autocorrelation — والـ threshold γ اللي بيوازن بين Pfa وPd.
- الرادار والسونار والـ GPS والاتصالات والـ ultrasound الطبي كلهم بيحلّوا نفس مسألة الـ detection/التقدير باستخدام الـ matched filtering والـ cross-correlation.
- مسافة الرادار بتتحدّد بتأخير الرحلة رايح جاي R = cτ/2؛ والـ range resolution بعد الـ pulse compression هي ΔR = c/(2B).
- الأهداف المتحرّكة بتعمل Doppler shift fd = 2v/λ. والـ ambiguity function بتوصف الـ resolution المشتركة في التأخير والـ Doppler.
- السونار الـ passive بيستخدم الـ cross-correlation بين sensorين عشان يلاقي الـ TDOA، وده بيدّي زاوية الـ azimuth لمصدر الصوت.
- في الاتصالات الرقمية، الـ matched filter بيحقق أقل BER = Q(√(2Eb/N0)) للـ signaling الـ antipodal الـ binary في AWGN.
- تقدير القناة بيستخدم pilot sequences الـ autocorrelation بتاعها شبه impulse. والـ correlation بين الـ signal المستقبَلة والـ pilot sequence بيرجّع الـ impulse response بتاعة القناة على طول.
- الـ acquisition بتاع GPS هو بحث correlation ثنائي الأبعاد على الـ code phase والـ Doppler. وطرق الـ FFT بتنزّل التعقيد من O(N²) لـ O(N log N) لكل Doppler bin.
- كل الأنظمة دي بتشترك في نفس سلسلة المعالجة: تبعت waveform معروفة ← تعمل correlation بينها وبين الـ signal المستقبَلة ← threshold ← تقدّر التأخير/الـ Doppler/الـ symbol.