Quantum algorithms for optimal control: large-scale finite-horizon linear-quadratic regulation
تقدم هذه الورقة خوارزميات كمومية تحل مسائل التنظيم الخطي-التربيعي ذات الأفق المحدود واسعة النطاق عن طريق إعادة صياغتها كأنظمة إطلاق متعدد (multiple-shooting) خطية ذات أعداد شرطية مستقلة عن الأفق، محققةً تدرجاً لوغاريتمياً متعدد الحدود في حجم النظام ومدة التحكم دون الحاجة إلى البناء الصريح لمؤثر حل ريكاتي.
تخيل عالماً يجب فيه على الآلات اتخاذ قرارات في أجزاء من الثانية للحفاظ على توازنها، أو حيث يجب على شبكات الطاقة أن تتكيف فوراً مع الطلب المتغير، أو حيث يجب على المركبات الفضائية أن تسلك مسارات معقدة بدقة متناهية. في كل هذه الحالات، يعتمد المهندسون على إطار رياضي يُعرف باسم التحكم الأمثل. الهدف بسيط: توجيه نظام من نقطة البداية إلى وجهة مرغوبة مع استخدام أقل قدر ممكن من الطاقة أو الوقت. لعقود من الزمن، كانت الأداة الأكثر موثوقية لحل هذه المشكلات هي طريقة تُسمى التنظيم الخطي التربيعي. وهي تعمل بشكل رائع عندما يسلك النظام سلوكاً في خط مستقيم يمكن التنبؤ به، وتنمو فيها تكلفة ارتكاب الأخطاء بطريقة سلسة ومتوقعة. ومع ذلك، عندما تكون النافذة الزمنية للمهمة ثابتة ومحدودة - كأن يكون هناك عدد محدد من الثواني لهبوط طائرة بدون طيار أو لاستقرار محطة طاقة - تصبح الرياضيات صعبة للغاية. فالطرق المطلوبة لإيجاد المسار المثالي تتضمن لغزاً معقداً وغير خطي يزداد صعوبة بشكل أسسي كلما كبر حجم النظام أو طالت النافذة الزمنية. وغالباً ما تصطدم الحواسيب الكلاسيكية، التي نستخدمها اليوم، بجدار هنا، حيث تجد صعوبة في معالجة الأرقام بالسرعة الكافية لتكون مفيدة للتطبيقات واسعة النطاق في الوقت الفعلي.
لقد طور فريق من الباحثين الآن نهجاً جديداً يستخدم مبادئ الحوسبة الكمومية لتجاوز هذه الاختناقات. فبدلاً من محاولة حل اللغز غير الخطي الصعب مباشرة، وجدوا طريقة لإعادة صياغة المشكلة بأكملها. لقد أخذوا المعادلات القياسية التي تصف كيفية حركة النظام وكيفية تراكم تكلفة الأخطاء، ودمجوها في نظام واحد ضخم من المعادلات الخطية. فكر في الأمر كأخذ عقدة متشابكة من الخيوط وفردها على طاولة؛ فبينما يصعب فك العقدة، يكون الترتيب المسطح سهلاً للقياس. ومن خلال تنظيم المشكلة بهذه الطريقة، حولوا تحدياً زمنياً فوضوياً إلى تحدٍ هيكلي ثابت. هذا الصياغة الجديدة تسمح لهم باستخدام خوارزمية كمومية متخصصة مصممة لحل الأنظمة الخطية، وهي مهمة يتفوق فيها الحاسوب الكمومي بشكل كبير على نظائره الكلاسيكية.
أثبت الباحثون أن هذه الطريقة الجديدة تعمل بكفاءة حتى للفترات الزمنية الطويلة جداً والأنظمة الكبيرة جداً، بشرط أن يكون النظام مستقراً بما يكفي ليتم التحكم فيه. وقد أظهروا أن صعوبة الحساب لا تنفجر مع طول النافذة الزمنية، وهي مشكلة واجهت المحاولات السابقة. وفي الواقع، أظهروا أنه بالنسبة لعائلة محددة من الأنظمة - مثل سلسلة طويلة من الكتل المتصلة بنوابض، وهو نموذج شائع لكل شيء من الجسور إلى السلاسل الجزيئية - فإن صعوبة الحساب تنمو ببطء شديد لدرجة أنها تصبح شبه مهملة مع كبر حجم النظام. وهذا يعني أن الحاسوب الكمومي يمكنه، من الناحية النظرية، إيجاد مدخلات التحكم المثلى لنظام ضخم في جزء ضئيل من الوقت الذي قد يستغرقه حاسوب كلاسيكي خارق.
ولضمان سلامة نظريتهم، أجرى الفريق تجارب عددية مكثفة على حاسوب كلاسيكي لمحاكاة كيفية سلوك طريقتهم. وقد اختبروا ثلاثة أنواع مختلفة من الأنظمة: سلسلة الكتل والنوابض، وخط من المركبات التي تتبع قائداً، وتفاعلاً كيميائياً ينتشر عبر مجال ما. وفي كل حالة، أكدت عمليات المحاكاة أن إعادة الصياغة الرياضية لطريقتهم كانت دقيقة وأن النظام ظل يعمل بشكل جيد، أو "جيد التكيّف"، بغض النظر عن طول مدة المحاكاة أو عدد أجزاء النظام. ووجدوا أن الطريقة يمكنها تقدير أفضل تكلفة ممكنة وإجراءات التحكم الدقيقة المطلوبة في أي لحظة معينة بدقة عالية. وتشير النتائج إلى أنه بالنسبة للأنظمة ذات عدد ثابت من مدخلات التحكم، فإن الوقت المطلوب لإيجاد الحل ينمو ببطء شديد مع حجم النظام لدرجة أنه يصبح ثابتاً فعلياً للأغراض العملية.
لا يدعي هذا العمل أنه قد حل كل مشكلات التحكم الموجودة، كما لا يدعي أن الحواسيب الكمومية جاهزة لتشغيل هذه الخوارزميات على الأجهزة المتاحة اليوم. ويؤكد الباحثون بوضوح أن نتائجهم تستند إلى إثباتات نظرية ومحاكاة. ومع ذلك، فقد وضعوا مساراً واضحاً للمضي قدماً. فمن خلال إثبات أن الهيكل الرياضي الأساسي لهذه المشكلات يمكن ترويضه وجعله فعالاً للمعالجة الكمومية، فقد أزالوا عائقاً نظرياً كبيراً. وتشير نتائجهم إلى أنه بمجرد نضوج الأجهزة الكمومية، سيكون من الممكن حل مشكلات التحكم واسعة النطلة والحساسة للوقت والتي تقع حالياً خارج نطاق الوصول، مما قد يحدث ثورة في كيفية إدارة البنية التحتية المعقدة، وتوجيه المركبات ذاتية القيادة، وتثبيت أنظمة الطاقة الحيوية. إن الرؤية الجوهرية هي أنه بتغيير المنظور تجاه المشكلة، يصبح المستحيل مجرد أمر صعب، ويصبح الصعب أمراً يمكن التعامل معه.
صياغة المشكلة يتناول هذا العمل التحدي الحسابي لحل مشكلة المنظم الخطي التربيعي (LQR) المستمر ذات المدى المحدود واسع النطاق. تتضمن المشكلة تقليل دالة التكلفة التربيعية عبر فترة زمنية محددة [0,T] خاضعة لديناميكيات خطية x˙(t)=Ax(t)+Bu(t). وخلافاً لـ LQR ذي المدى اللانهائي، الذي يختزل في معادلة ريكاتي جبرية (ARE)، تتطلب حالة المدى المحدود حل معادلة ريكاتي تفاضلية (DRE) تعتمد على الزمن. وتعد الطرق العددية الكلاسيكية لمعادلة DRE مكلفة حسابياً بالنسبة لأبعاد الحالة الكبيرة (n) والآفاق الزمنية الطويلة (T)، حيث غالباً ما تتدرج حدود زمنية متعددة أو أسوأ مع هذه المعلمات. علاوة على ذلك، أظهرت النهج الكمومية الحالية للمعادلات التفاضلية المصفوفية غير الخطية، مثل تلك القائمة على تحويلات موبيوس، أنها تعاني من أرقام شرطية (condition numbers) تزداد أسياً مع الأفق الزمني T، مما يجعلها غير عملية لمهام التحكم طويلة الأمد.
المنهجية يقترح المؤلفون إطار عمل لخوارزمية كمومية يتجنب الحل المباشر لمعادلة ريكاتي غير الخطية. بدلاً من ذلك، تعيد المنهجية صياغة مشكلة التحكم الأمثل باستخدام مبدأ بونترياغين الأقصى (PMP)، والذي ينتج عنه مسألة قيم حدية نقطتين خطية (TPBVP) تتضمن الحالة x(t) والمتغير المترفق (costate) λ(t).
إعادة صياغة الرماية المتعددة (Multiple-Shooting Reformulation): يتم تمثيل مسألة TPBVP المستمرة بشكل متقطع باستخدام مخطط "غاوس-ليجندر رونج-كوتا". تُجمع متغيرات الحالة والمتغير المترفق عند نقاط الشبكة الزمنية في متجه q، وتُجمع معادلات التطور مع الشروط الحدية في نظام خطي متفرق كبير Lexq=b. هذا النهج "للرماية المتعددة" يحول مشكلة التحكم التفاضلية إلى مسألة جبر خطي.
تحليل التكييف عبر فك الارتباط بين ريكاتي وليابونوف: المساهمة النظرية الحاسمة هي الإثبات بأن الرقم الشرطي للنظام الخطي الناتج، κ(Lex)، محدود بشكل موحد بالنسبة للأفق الزمني T، بشرما كان النظام يستوفي فرضيات القابلية للاستقرار والقدرة على الكشف. يحقق المؤلفون ذلك من خلال تطبيق تغيير متغيرات "ريكاتي-ليابونوف" الذي يفك ارتباط النظام إلى أنظمة فرعية مستقرة أمامياً ومستقرة خلفياً. يعمل هذا التحويل على فصل الأنماط المتنامية، مما يضمن اضمحلال تأثير نقاط الشبكة أسياً مع تباعدها، وبالتالي منع الانفجار الأسي للرقم الشرطي مع T.
خوارزمية الأنظمة الخطية الكمومية (QLSA): يتم حل النظام الخطي باستخدام خوارزمية الأنظمة الخطية الكمومية (QLSA) القائمة على تحويل القيمة المفردة الكمومي (QSVT) وتقنيات الترميز الكتلي (block-encoding). تقوم الخوارزمية ببناء ترميزات كتلية لمصفوفات النظام ومعكوس المؤثر المتفرق.
استخراج المخرجات: لا تخرج الخوارزمية المسار الكلاسيكي الكامل. بدلاً من ذلك، توفر تقديرات كلاسيكية للتكلفة المثلى J∗ ومتجه التحكم الأمثل u∗(tk) عند نقطة شبكية زمنية واحدة مختارة tk. تُستخرج هذه الكميات عبر اختبارات هادامارد وتقدير السعة المطبق على الحالة الكمومية التي تمثل الحل.
المساهمات والنتائج الرئيسية
التكييف الموحد مع الأفق: تثبت الورقة أن الرقم الشرطي لنظام الرماية المتعددة هو رقم مستقل بالأفق الزمني T. يعتمد هذا النتيجة على الاضمحلال الأسي للمنتشر المغلق (closed-loop propagator) تحت فرضيات القابلية للاستقرار والقدرة على الكشف، وهي خاصية ثبت أنها تتحقق بشكل موحد لعائلة من سلاسل الكتلة والزنبرك المقترنة.
تدرج التعقيد: تحت افتراض الترميز الكتلي للمصفوفات المتفرقة ومعلمات الاستقرار المحدودة، تتدرج تعقيدات الاستعلام والبوابات للخوارزميات الكمومية المقترحة لوغاريتمياً مع بُعد الحالة n والأفق الزمني T. وتحديداً:
تقدير التكلفة المثلى: O~(ε−1polylog(n,d,T)).
تقدير التحكم الأمثل: O~(d3/2ε−1polylog(n,d,T)). حيث ε هو الخطأ المستهدف، و d هو بُعد التحكم، و O~ تخفي العوامل اللوغاريتمية في الخطأ واحتمالية الفشل.
مثال موحد البعد: يوضح المؤلفون أنه بالنسبة لعائلة محددة من سلاسل الكتلة والزنبرك المقترنة مع تحكمات طرفية ثابتة، تظل معلمات الاستقرار ومعايير المصفوفات محدودة بشكل مستقل عن عدد الكتل. وهذا يؤكد أن التدرج اللوغاريتمي ينطبق على الأنظمة واسعة النطاق ذات الصلة فيزيائياً ذات الطاقة الأولية الإجمالية المحدودة.
التحقق العددي: تؤكد التجارب العددية الكلاسيكية على سلاسل الكتلة والزنبرك، وأفواج المركبات، وأنظمة التفاعل والانتشار، حدود التكييف النظرية ودقة تمثيل "غاوس" المتقطع. وتؤكد التجارب أن الرقم الشرطي يظل مستقراً مع زيادة الأفق، وأن خطأ التقطيع يتناقص بسرعة مع عدد مراحل "رونج-كوتا".
الأهمية والادعاءات تدعي الورقة أنها توفر خطوة تأسيسية نحو التسريع الكمومي للتحكم الأمثل واسع النطاق. من خلال إعادة صياغة مشكلة LQR ذات المدى المحدود كمسألة نظام خطي ذات رقم شرطي مستقل عن الأفق، يتغلب المؤلفون على عائق رئيسي حد سابقاً من النهج الكمومية في هذا المجال (تحديداً الاعتماد الأسي على T الموجود في طرق خطية معادلة ريكاتي التفاضلية السابقة).
تكمن الأهمية في القدرة على تقدير التكاليف والتحكمات المثلى للأنظمة عالية الأبعاد عبر فترات زمنية طويلة بموارد تتدرج بكفاءة (لوغاريتمياً) بدلاً من التدرج الأسي. ويؤكد العمل أن حلول ريكاتي وليابونوف تُستخدم فقط للتحليل النظري للتكييف وليس لحسابها بواسطة الخوارزمية الكمومية، مما يتجنب العبء الحسابي لحل المعادلات غير الخطية على حاسوب كمومي. ويحافظ المؤلفون على نطاق متواضع، بالتركيز على تقدير كميات محددة عند نقاط زمنية مختارة بدلاً من إعادة بناء مسار التحكم الكلاسيكي الكامل، وهو ما يتطلب القراءة عند كل عقدة. تُقدم النتائج كإطار نظري مدعوم بالأدلة العددية، مع اقتراح أعمال مستقبلية لتهيئة الأنظمة ضعيفة التخميد وتحسين حدود الموارد لوظائف تحكم محددة.