An Exponentially stable Extended Kalman Filter with Estimate dependent Process noise Covariance for Chemical Reaction Networks

ملخص تقني: مرشح كالمان الموسع مستقر أسياً مع تغاير ضوضاء عملية يعتمد على التقدير لشبكات التفاعل الكيميائي

1. بيان المشكلة

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

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

2. المنهجية

يقترح المؤلفون إطار عمل جديد لـ EKF حيث يتم تحديث تغاير ضوضاء العملية عند كل خطوة زمنية بناءً على تقدير الحالة الحالي، والمستمد من معادلة لانجفان الكيميائية (CLE).

• إطار النموذج: يتم نمذجة النظام كعملية عشوائية منفصلة زمنياً حيث يكون حد الانتشار (الضوضاء) دالة لمتجه الميل $A(x_k)$. ويُعرف التغاير المقترح كما يلي:
  $$Q_k = \delta V \tilde{A}(\hat{x}^+_k) \tilde{A}(\hat{x}^+_k)^T V^T$$
  حيث $\delta$ هو فترة أخذ العينات، و $V$ هي مصفوفة القياس المتسق (stoichiometric matrix)، و $\tilde{A}$ هو الجذر التربيعي لمتجه الميل.
• تحليل الاستقرار: لضمان موثوقية المرشح، يقوم المؤلفون بإجراء تحليل استقرار عشوائي صارم. يهدفون إلى إثبات أن خطأ التقدير $e_k$ هو مقيد أسياً في المعنى التربيعي المتوسط، مما يعني:
  $$E(\|e_k\|^2) \leq C_1 \gamma^k E(\|e_0\|^2) + C_0$$
  حيث $0 < \gamma < 1$.
• المنهج الرياضي: يستخدم الإثبات دوال تشبه دالة ليابونوف (Lyapunov-like functions) ومنطق السوبر-مارنتينج (supermartingale). يستنتج المؤلفون متباينة حدودية بدلالة الفاصل الزمني لأخذ العينات $\delta$. ومن خلال تطبيق قاعدة ديكارت لعلامات الإشارة، يثبتون وجود عتبة موجبة فريدة ($\delta_{max}$) لفترة أخذ العينات. إذا تم الإبقاء على الفاصل الزمني لأخذ العينات $\delta$ تحت هذه العتبة، فإن المرشح يضمن الاستقرار العشوائي.

3. المساهمات الرئيسية

1. تصميم مرشح قائم على المبادئ الأولى: بدلاً من الضبط التجريبي، يقدم المؤلفون طريقة لتصميم تغاير ضوضاء العملية بناءً على الخصائص الفيزيائية للتفاعلات الكيميائية (CLE).
2. إثبات الاستقرار: يقدمون إثباتاً رياضياً رسمياً بأن EKF مع تغاير ضوضاء يعتمد على الحالة يكون مستقراً أسياً من الناحية العشوائية.
3. حد فترة أخذ العينات: يستنتجون حداً أعلى محدداً للفاصل الزمني لأخذ العينات ($\delta < \delta_{max}$) مطلوب للحفاظ على استقرار المرشح، مما يجسّر الفجوة بين الديناميكيات البيولوجية في الزمن المستمر والتقدير في الزمن المنفصل.

4. النتائج والتحقق

تم التحقق من الإطار النظري باستخدام نموذج تعبير جيني معقد وغير خطي (يتضمن بوليميراز الحمض النووي الريبي، والحمض النووي الريبي المرسال، والريبوسومات، والبروتينات).

• العتبة العددية: بالنسبة لمعلمات التعبير الجيني المحددة المستخدمة، حسب المؤلفون الحد الأقصى المسموح به للفاصل الزمني لأخذ العينات $\delta_{max} \approx 0.00054$.
• أداء المحاكاة: باستخدام فاصل زمني لأخذ العينات قدره $\delta = 0.0005$ (والذي يستوفي شرط الاستقرار)، أظهرت عمليات المحاكاة:
  • التتبع الدقيق: تتبعت أعداد جميع الجزيئات الحيوية المقدرة (RNA، mRNA، إلخ) القيم الحقيقية بدقة رغم الضوضاء الذاتية العالية.
  • خطأ مقيد: أظهرت مصفوفة تغاير الخطأ اللاحق ومتوسط مربع الخطأ (MSE) سلوكاً تذبذبياً ولكن مقيداً، مما يؤكد أن الخطأ لا يتشتت.
  • التقارب: أظهر متوسط مربع الخطأ (MSE) اتجاهاً تناقصياً أسياً نحو حالة مستقرة مقيدة، وهو ما يتوافق مع التنبؤ النظري بالاستقرار الأسي.

5. الأهمية

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