نظرية التحكم1960متوسط12 دقيقة قراءة

نهج جديد لمسائل الترشيح الخطي والتنبُّؤ

A New Approach to Linear Filtering and Prediction Problems

Kalman, R. E. — ASME Journal of Basic Engineering

المشكلة

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

الإسهام

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

الأثر

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

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

لو اعتمد على البوصلة وحدها فسينحرف تدريجياً. ولو اعتمد على السونار وحده فستقفز قراءاته عشوائياً. لكنه حين يُرجّح كل مصدر بمقدار ثقته فيه ويمزجهما معاً، يحصل على تقدير لموقعه أدقّ ممّا يستطيع أيٌّ منهما تقديمه بمفرده.

هذا بالضبط ما يفعله مرشّح كالمن: يمزج تنبُّؤاً ناقصاً مع قياس مُشوَّش بطريقة منهجية، ويُعدّل درجة الثقة بكلٍّ منهما تلقائياً في كل خطوة.

السقف: مرشّح وينر وحدوده

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

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

ما كان يحتاجه المهندسون هو مرشّح يعمل في المجال الزمني، ويتعامل مع أنظمة تتغيّر خصائصها الإحصائية، ويُحدَّث تكرارياً قياساً بقياس، ويتوسّع بسهولة إلى أبعاد متعددة.

الفكرة المحورية: التفكير في فضاء الحالة

الفكرة الجوهرية عند كالمن كانت تغيير طريقة التفكير بالمسألة بالكامل. بدلاً من تحليل الترددات كما فعل وينر، صاغ كالمن المسألة في فضاء الحالة: السؤال لم يَعُد «ما المُركّبات الترددية لهذه الإشارة؟» بل أصبح «ما للنظام الذي أنتج هذا القياس المُشوَّش؟»

الفكرة في جوهرها بسيطة: العالَم الحقيقي له حالة داخلية (موضع، سرعة، درجة حرارة…) لا نستطيع قياسها مباشرةً، وإنما نرصدها من خلال مستشعرات مُشوَّشة. بين كل قياسين تتطوّر هذه الحالة وفق قوانين فيزيائية معروفة أو وفق نموذج رياضي. وكلٌّ من التطوّر والقياس تعتريه ضوضاء عشوائية.

معادلتان فقط تصفان هذا العالَم:

xk=Fxk1+Buk1+wk1x_{k} = F\,x_{k-1} + B\,u_{k-1} + w_{k-1}
معادلة انتقال الحالة — كيف تتطوّر الحالة المخفيةF = مصفوفة انتقال الحالة (الفيزياء) · B u = مُدخل تحكم اختياري · w = ضوضاء العملية (عدم اليقين في النموذج ذاته، مسحوبة من 𝒩(0, Q))
zk=Hxk+vkz_{k} = H\,x_{k} + v_{k}
معادلة القياس — ما يراه المستشعرH = مصفوفة القياس (تربط الحالة المخفية بما يمكن ملاحظته) · v = ضوضاء القياس (عدم دقة المستشعر، مسحوبة من 𝒩(0, R))

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

افتح في المختبر
العالَم ذو الطبقتين: الحالة المخفية تتطوّر عبر F وتُرصد عبر H. بدِّل لترى تأثير الضوضاء.
تستيقظ التجربة عند وصولك…

الخوارزمية: تنبّأ ثم حدِّث

مرشّح كالمن يعمل كحلقة بسيطة تتكرر مع كل قياس جديد، وفيها خطوتان فقط:

الخطوة 1 — التنبُّؤ. خذ النموذج الفيزيائي واستعمله لتُسقط تقدير الحالة الحالي وعدم اليقين المرافق له خطوة إلى الأمام في الزمن. ما تحصل عليه الآن هو تقدير مُسبق — أفضل تخمين لديك قبل أن تصلك أي بيانات جديدة. تخيّل بحّاراً يحسب موقعه من البوصلة فقط: «بناءً على اتجاهي وسرعتي، يُفترض أنني هنا

الخطوة 2 — التحديث. عندما يصل القياس الجديد، احسب كسب كالمن — وهو معامل يُحدّد كم تثق بالقياس مقارنةً بالتنبُّؤ. ادمج الاثنين بالنسبة المناسبة، ثم قلِّص عدم اليقين: لأنك الآن تملك معلومات أكثر ممّا كنت تملكه قبل القياس.

النقطة المفتاحية هنا أن التنبؤ والقياس كلاهما يحمل قدراً من عدم اليقين، والمرشّح يوازن بينهما بالشكل الأمثل. حين يكون النموذج واثقاً والمستشعر مُشوَّشاً، يميل المرشّح نحو النموذج. وحين يكون المستشعر دقيقاً والنموذج غير متيقّن، يميل نحو المستشعر. كسب كالمن هو الآلية التي تجد نقطة التوازن المثلى تلقائياً.

افتح في المختبر
تابع دورة التنبؤ-التحديث على هدف متحرك. القطع الأزرق هو عدم يقين التنبؤ، والبرتقالي عدم يقين القياس، والأخضر التقدير المدموج.
تستيقظ التجربة عند وصولك…

الصياغة الرياضية: خمس معادلات تُشغّل العالَم

كل ما في مرشّح كالمن يتلخّص في خمس معادلات، موزّعة بين مرحلتي التنبؤ والتحديث. كل معادلة تؤدي دوراً محدداً — لا شيء هنا للزينة.

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

كسب كالمن: قلب المرشّح

كسب كالمن KK هو مصفوفة تجيب عن سؤال واحد: «بأيّ مقدار أُصحّح تنبُّؤي؟» تخيّله قرصاً يدور بين حالتين متطرفتين:

  • K0K \to 0: تجاهل القياس كلياً — «نموذجي ممتاز والمستشعر لا يُعتمد عليه.»
  • KIK \to I: تجاهل التنبؤ كلياً — «نموذجي لا يعرف شيئاً، المستشعر هو المرجع الوحيد.»

في الواقع لا يبلغ KK أياً من الحدّين، بل يقع في منطقة وسطى ويتعدّل تلقائياً في كل خطوة زمنية وفقاً لعدم اليقين النسبي. في المراحل الأولى حين يكون تقديرك خشناً، يكون KK كبيراً لأن كل قياس يحمل معلومات كثيرة. ومع تراكم القياسات وتحسُّن التقدير، يتقلّص KK تدريجياً — فأنت تعرف تقريباً أين أنت، وقراءة مُشوَّشة واحدة لن تُحرّك تقديرك كثيراً.

Kk=Pkk1HTHPkk1HT+RK_k = \frac{P_{k|k-1}\,H^T}{H\,P_{k|k-1}\,H^T + R}
كسب كالمن — وزن المزج الأمثل (الصيغة العددية)البسط = مقدار عدم يقين التنبؤ (في فضاء القياس) · المقام = عدم اليقين الكلي (تنبؤ + قياس) · حين R → 0 يؤول K → 1 (ثق بالمستشعر). حين P → 0 يؤول K → 0 (ثق بالنموذج).

لماذا التوزيع الغاوسي؟ الأمثلية وجمال الحل المغلق

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

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

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

تتبُّع جسم متحرك: المرشّح أثناء العمل

لنأخذ مثالاً بسيطاً: كرة تتحرك في بُعد واحد. الحالة هنا هي [الموضع,السرعة]T[\text{الموضع},\, \text{السرعة}]^T. بين كل قياسين تتحرك الكرة وفق الموضعالجديد=الموضع+السرعة×Δt\text{الموضع}_{\text{الجديد}} = \text{الموضع} + \text{السرعة} \times \Delta t. والمستشعر المتاح يقيس الموضع فقط — لا السرعة — مع قدر من الضوضاء.

في العرض التفاعلي أدناه: النقاط الحمراء هي قراءات المستشعر المُشوَّشة، والخط الأزرق هو تقدير مرشّح كالمن. لاحظ أن الخط الأزرق أنعم بكثير من النقاط الحمراء وأقرب إلى المسار الحقيقي (الأخضر) — وهذا الفرق يظهر بوضوح أكبر كلما زادت ضوضاء المستشعر.

افتح في المختبر
اضبط ضوضاء المستشعر (R) وضوضاء العملية (Q) لترى كيف يتكيّف المرشّح. R عالية → خط ناعم يثق بالنموذج. Q عالية → خط متعرّج يثق بالقياسات.
تستيقظ التجربة عند وصولك…

الفكرة ذاتها في شيفرة برمجية

مرشّح كالمن لتتبُّع موضع أحادي البُعد، كاملاًpython

مبسَّط لإظهار الفكرة — ليس التنفيذ الحقيقي.

import numpy as np

def kalman_filter(measurements, F, H, Q, R, x0, P0):
    """
    F: مصفوفة انتقال الحالة          (n x n)
    H: مصفوفة القياس                 (m x n)
    Q: مصفوفة تغاير ضوضاء العملية     (n x n)
    R: مصفوفة تغاير ضوضاء القياس      (m x m)
    x0: تقدير الحالة الابتدائي         (n x 1)
    P0: مصفوفة التغاير الابتدائية      (n x n)
    """
    x = x0.copy()
    P = P0.copy()
    estimates = []

    for z in measurements:
        # ── التنبؤ ─────────────────────────────────
        x_pred = F @ x                    # تنبؤ الحالة
        P_pred = F @ P @ F.T + Q          # تنبؤ التغاير

        # ── التحديث ──────────────────────────────────
        y = z - H @ x_pred                # الابتكار (المفاجأة)
        S = H @ P_pred @ H.T + R          # تغاير الابتكار
        K = P_pred @ H.T @ np.linalg.inv(S)  # كسب كالمن

        x = x_pred + K @ y               # الحالة المُصحَّحة
        P = (np.eye(len(x)) - K @ H) @ P_pred  # التغاير المُصحَّح

        estimates.append(x.copy())

    return np.array(estimates)

# مثال: تتبُّع الموضع والسرعة من مستشعر لا يقيس إلا الموضع
dt = 1.0
F = np.array([[1, dt],    # الموضع = الموضع + السرعة × dt
               [0,  1]])   # السرعة = السرعة (نموذج سرعة ثابتة)
H = np.array([[1, 0]])    # المستشعر يقيس الموضع فقط
Q = np.array([[0.1, 0], [0, 0.1]])  # ضوضاء عملية صغيرة
R = np.array([[1.0]])     # تباين ضوضاء المستشعر
x0 = np.array([[0], [1]])  # ابدأ عند 0 بسرعة 1
P0 = np.eye(2) * 10       # عدم يقين كبير في البداية

كالمن مقابل وينر: ماذا تغيّر

ما فعله كالمن لم يكن تحسيناً تدريجياً لمرشّح وينر — بل كان تغييراً جذرياً في طريقة التفكير بالمسألة:

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

التغاير: المرشّح يعرف ما لا يعرفه

من أهمّ ما يُميّز مرشّح كالمن أنه لا يُعطيك تقديراً فقط، بل يُرفقه بـمؤشر لمدى ثقتك فيه. مصفوفة التغاير PP تُخبرك تحديداً بمقدار عدم اليقين في كل بُعد من أبعاد الحالة.

في خطوة التنبؤ يكبر PP — لأنك تبتعد عن آخر قياس فعلي وتعتمد على نموذج غير مثالي. وفي خطوة التحديث يتقلّص PP — لأن القياس الجديد أضاف معلومات. مع مرور الوقت يتقارب PP نحو قيمة مستقرة تعكس التوازن الأساسي بين ضوضاء العملية وضوضاء القياس.

هذا «الوعي الذاتي» هو ما يجعل مرشّح كالمن أقوى بكثير من أساليب بسيطة كالمتوسط المتحرك أو التنعيم الأُسِّي: تلك تُعطيك تقديراً لكن دون أي مقياس لمدى موثوقيته.

الأثر: من القمر إلى جيبك

  1. 1960

    كالمن ينشر الورقة

    رودولف كالمن يقدّم المُقدِّر التكراري للحالة في مجلة ASME. الورقة تعيد صياغة الترشيح من مجال التردد إلى فضاء الحالة.

  2. 1961

    امتداد كالمن-بوسي للزمن المتصل

    كالمن وبوسي يُوسّعان المرشّح ليشمل أنظمة الزمن المتصل، ممّا يفتح الباب لتطبيقه على الأنظمة التناظرية والمعادلات التفاضلية.

  3. 1962

    نظام ملاحة أبولو

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

  4. 1973

    تصميم نظام GPS

    بنية نظام GPS تدمج مرشّح كالمن في جوهرها لمزج إشارات الأقمار الصناعية مع القياسات القصورية. وصف أحدهم نظام GPS بأكمله بأنه «مرشّح كالمن واحد هائل».

  5. 1970s

    نماذج ARIMA والسلاسل الزمنية

    نماذج بوكس-جينكنز ARIMA ترتبط ارتباطاً وثيقاً بمرشّح كالمن. أي نموذج ARIMA يمكن صياغته في فضاء الحالة وتقديره عبر مرشّح كالمن، ممّا يربط نظرية التحكم بالاقتصاد القياسي.

  6. 2000s

    دمج مستشعرات الهواتف الذكية

    كل هاتف ذكي يدمج بيانات مقياس التسارع والجيروسكوب ومقياس المغناطيسية وGPS عبر مرشّحات كالمن لتوفير توجّه سلس وعدّ خطوات وتتبُّع موقع دقيق.

  7. 2010s

    المركبات ذاتية القيادة

    السيارات ذاتية القيادة تستخدم مرشّحات كالمن الممتدة وغير المُقادة لدمج الكاميرا والليدار والرادار وGPS في تقدير متماسك لحالة المركبة والعالم المحيط بها.

  8. 2021

    S4 ونماذج فضاء الحالة العميقة

    نموذج فضاء الحالة المُهيكَل للتسلسلات (S4) يُعيد أفكار فضاء الحالة إلى التعلم العميق محققاً نتائج قوية في نمذجة التسلسلات بعيدة المدى — سليلٌ فكري مباشر لصياغة كالمن في فضاء الحالة.

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

ما وراء الخطية: الامتدادات والقيود

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

  • مرشّح كالمن الممتد (EKF): يُقرّب النظام اللاخطي بتمديد تايلور من الدرجة الأولى حول التقدير الحالي. هو الأكثر شيوعاً ويُستخدم في معظم أجهزة GPS وملاحة المركبات الفضائية. بسيط لكنه قد يفشل حين تكون اللاخطية حادة.
  • مرشّح كالمن غير المُقاد (UKF): بدلاً من التقريب الخطي، يختار مجموعة صغيرة من «نقاط سيغما» تُمثّل المتوسط والتغاير، ويُمرّرها عبر الدالة اللاخطية الفعلية، ثم يحسب الإحصاءات من النتائج. أكثر دقة من EKF في الأنظمة شديدة اللاخطية.
  • مرشّحات الجسيمات: تتخلّى عن فرضية الغاوسية تماماً. تُمثّل التوزيع الاحتمالي بسحابة من العيّنات المُوزَّنة (الجسيمات). قادرة على التعامل مع أي شكل من اللاخطية وأي نوع من الضوضاء، لكن بتكلفة حسابية أعلى بكثير.

كل امتداد يأتي بمقايضة بين الدقة والكلفة الحسابية. ومرشّح كالمن الخطي الأصلي يظل الخيار الأفضل كلما صحّت فرضيتا الخطية والغاوسية — وهذا يحصل في التطبيقات العملية أكثر ممّا قد تتوقع.

المرجعKalman, R. E.. A New Approach to Linear Filtering and Prediction Problems. ASME Journal of Basic Engineering, 1960.

مصطلحات هذه الورقة