In this paper, we investigate the potential of reconfigurable intelligent surfaces (RISs) in facilitating passive/device-free three-dimensional (3D) drone localization within existing cellular infrastructure operating at millimeter-wave (mmWave) frequencies and employing multiple antennas at the transceivers. The developed localization system operates in the bi-static mode without requiring direct communication between the drone and the base station. We analyze the theoretical performance limits via Fisher information analysis and Cramér Rao lower bounds (CRLBs). Furthermore, we develop a low-complexity yet effective drone localization algorithm based on coordinate gradient descent and examine the impact of factors such as radar cross section (RCS) of the drone and training overhead on system performance. It is demonstrated that integrating RIS yields significant benefits over its RIS-free counterpart, as evidenced by both theoretical analyses and numerical simulations.
معرّف الورقة : 2404.14879العنوان : Device-Free 3D Drone Localization in RIS-Assisted mmWave MIMO Networksالمؤلفون : Jiguang He, Charles Vanwynsberghe, Hui Chen, Chongwen Huang, Aymen Fakhreddineالتصنيف : eess.SP (معالجة الإشارات)تاريخ النشر : أبريل 2024 (نسخة arXiv مسبقة، محدثة في 2 يناير 2025)رابط الورقة : https://arxiv.org/abs/2404.14879 تبحث هذه الورقة في إمكانية استخدام الأسطح الذكية القابلة لإعادة التشكيل (RIS) في نطاق موجات ميلليمترية (mmWave) مع البنية التحتية الخلوية متعددة الهوائيات لتحقيق تحديد موقع الطائرات بدون أجهزة/سلبي ثلاثي الأبعاد (3D). يعمل نظام التحديد الموضع المطور في الوضع ثنائي الأساس دون الحاجة إلى اتصال مباشر بين الطائرة والمحطة الأساسية. تم تحليل حدود الأداء النظرية من خلال تحليل معلومات فيشر وحد كرامير-راو السفلي (CRLB). بالإضافة إلى ذلك، تم تطوير خوارزمية تحديد موقع طائرة منخفضة التعقيد لكن فعالة بناءً على الانحدار الإحداثي، وتم دراسة تأثير مقطع الرادار العاكس (RCS) للطائرة والنفقات التدريبية على أداء النظام. تُظهر التحليلات النظرية والمحاكاة الرقمية مزايا كبيرة للنظام المدمج مع RIS مقارنة بالأنظمة بدون RIS.
المشكلة الأساسية التي تعالجها هذه الورقة هي: تحقيق تحديد موقع سلبي ثلاثي الأبعاد للطائرات الحرة في البنية التحتية الخلوية الموجودة، أي تحديد موقعها المكاني بدقة دون الحاجة إلى اتصال نشط بين الطائرة والمحطة الأساسية.
متطلبات الأمان : مع انتشار تطبيقات الطائرات بدون طيار، يزداد الطلب على الكشف والتحديد الموضع للطائرات غير المصرح بها أو الخبيثةإعادة استخدام البنية التحتية : الاستفادة الكاملة من البنية التحتية لشبكات الجيل الخامس/السادس الخلوية الموجودة، مما يتجنب نشر أنظمة رادار متخصصة إضافيةدمج التكنولوجيا : تطبيق مهم لنموذج دمج الاتصالات والاستشعار (ISAC)أنظمة الرادار التقليدية : تتطلب معدات رادار متخصصة، وتكاليف عالية ونطاق تغطية محدودأنظمة التحديد الموضع النشطة : تتطلب اتصالاً بين الهدف والمحطة الأساسية، غير مناسبة للأهداف غير التعاونيةالأنظمة السلبية الموجودة : الأداء محدودة في نطاق موجات ميلليمترية، دقة التحديد الموضع غير كافيةالاستفادة من RIS لتعزيز مسارات انتشار الإشارات وتحسين أداء التحديد الموضع السلبي تحقيق تحديد موقع ثلاثي الأبعاد عالي الدقة في أنظمة MIMO موجات ميلليمترية تطوير خوارزميات منخفضة التعقيد عملية تقترب من حدود الأداء النظرية تصميم النظام : اقتراح معمارية نظام تحديد موقع طائرة سلبي بموجات ميلليمترية مدعومة بـ RISالتحليل النظري : تحليل حدود الأداء النظرية للنظام من خلال مصفوفة معلومات فيشر و CRLBتطوير الخوارزمية : تطوير خوارزمية تحديد موقع ثلاثية المراحل منخفضة التعقيد بناءً على الانحدار الإحداثي (CGD)تقييم الأداء : تحليل شامل لتأثير المعاملات الرئيسية مثل RCS والنفقات التدريبية على أداء النظامالتحقق من المزايا : إثبات من خلال النظرية والمحاكاة المزايا الكبيرة لنظام RIS المدعوم مقارنة بالطرق التقليديةالإدخال : إشارات استقبال متعددة الفترات الزمنية في جهاز المستخدم المستقبل (UE)
الإخراج : إحداثيات الموقع ثلاثية الأبعاد للطائرة ( x D , y D , z D ) (x_D, y_D, z_D) ( x D , y D , z D ) القيود : لا يوجد اتصال مباشر بين الطائرة والمحطة الأساسية، يعمل النظام في الوضع ثنائي الأساس
يتكون النظام من ثلاثة مكونات رئيسية:
محطة موجات ميلليمترية أساسية (BS) : مجهزة بمصفوفة مستوية موحدة من M B = M B , y × M B , z M_B = M_{B,y} \times M_{B,z} M B = M B , y × M B , z هوائيسطح ذكي قابل لإعادة التشكيل (RIS) : يحتوي على M R = M R , x × M R , y M_R = M_{R,x} \times M_{R,y} M R = M R , x × M R , y وحدة انعكاسجهاز مستخدم استقبال (UE) : مجهز بـ M U M_U M U هوائييتم استخدام نموذج قناة Saleh-Valenzuela، مع الأخذ في الاعتبار خصائص التشتت الضئيلة في نطاق موجات ميلليمترية. القنوات الرئيسية تشمل:
قناة BS-RIS H 1 ∈ C M R × M B \mathbf{H}_1 \in \mathbb{C}^{M_R \times M_B} H 1 ∈ C M R × M B :
H 1 = e − j 2 π d 1 / λ ρ 1 α x ( θ r , 1 , ϕ r , 1 ) ⊗ α y ( θ r , 1 , ϕ r , 1 ) ( α y ( θ t , 1 , ϕ t , 1 ) ⊗ α z ( ϕ t , 1 ) ) H \mathbf{H}_1 = \frac{e^{-j2\pi d_1/\lambda}}{\sqrt{\rho_1}} \boldsymbol{\alpha}_x(\theta_{r,1}, \phi_{r,1}) \otimes \boldsymbol{\alpha}_y(\theta_{r,1}, \phi_{r,1}) \left(\boldsymbol{\alpha}_y(\theta_{t,1}, \phi_{t,1}) \otimes \boldsymbol{\alpha}_z(\phi_{t,1})\right)^H H 1 = ρ 1 e − j 2 π d 1 / λ α x ( θ r , 1 , ϕ r , 1 ) ⊗ α y ( θ r , 1 , ϕ r , 1 ) ( α y ( θ t , 1 , ϕ t , 1 ) ⊗ α z ( ϕ t , 1 ) ) H قنوات أخرى : قناة BS-الطائرة h 2 \mathbf{h}_2 h 2 ، قناة RIS-الطائرة h 3 \mathbf{h}_3 h 3 ، قناة الطائرة-UE h 4 \mathbf{h}_4 h 4 ، وغيرهاتستخدم المحطة الأساسية طريقة متعددة الحزم:
حزمة ثابتة f 0 \mathbf{f}_0 f 0 موجهة نحو RIS حزم متغيرة بالزمن f k \mathbf{f}_k f k موجهة نحو منطقة السماء حزم RIS متغيرة بالزمن ω k \boldsymbol{\omega}_k ω k تغطي منطقة الهدف نموذج الإشارة المستقبلة:
y k = ζ h 4 ( h 3 diag ( ω k ) H 1 F k s k + h 2 F k s k ) + H 5 F k s k + n k \mathbf{y}_k = \zeta\mathbf{h}_4(\mathbf{h}_3\text{diag}(\boldsymbol{\omega}_k)\mathbf{H}_1\mathbf{F}_k\mathbf{s}_k + \mathbf{h}_2\mathbf{F}_k\mathbf{s}_k) + \mathbf{H}_5\mathbf{F}_k\mathbf{s}_k + \mathbf{n}_k y k = ζ h 4 ( h 3 diag ( ω k ) H 1 F k s k + h 2 F k s k ) + H 5 F k s k + n k
حيث ζ \zeta ζ هو مقطع الرادار العاكس (RCS) للطائرة.
استخدام طريقة الانحدار الإحداثي (CGD) لتقدير القناة المتسلسلة:
دالة الهدف:
min f ( h 2 , h 3 , h 4 ) = ∥ Y − P / 2 H ~ Ω ˉ − P / 2 H ^ F ˉ ∥ F 2 \min f(\mathbf{h}_2,\mathbf{h}_3,\mathbf{h}_4) = \left\|\mathbf{Y} - \sqrt{P/2}\tilde{\mathbf{H}}\bar{\boldsymbol{\Omega}} - \sqrt{P/2}\hat{\mathbf{H}}\bar{\mathbf{F}}\right\|_F^2 min f ( h 2 , h 3 , h 4 ) = Y − P /2 H ~ Ω ˉ − P /2 H ^ F ˉ F 2
التحديث التكراري:
h 4 ( k ) = h 4 ( k − 1 ) − η ∂ f ∂ h 4 ∣ h 4 = h 4 ( k − 1 ) \mathbf{h}_4^{(k)} = \mathbf{h}_4^{(k-1)} - \eta \frac{\partial f}{\partial \mathbf{h}_4}\bigg|_{\mathbf{h}_4=\mathbf{h}_4^{(k-1)}} h 4 ( k ) = h 4 ( k − 1 ) − η ∂ h 4 ∂ f h 4 = h 4 ( k − 1 )
استخراج معاملات الزاوية من خلال البحث ثنائي الأبعاد:
θ ^ t , 2 , ϕ ^ t , 2 = arg max θ t , 2 , ϕ t , 2 ∣ h ^ 2 ( α y ( θ t , 2 , ϕ t , 2 ) ⊗ α z ( ϕ t , 2 ) ) ∣ \hat{\theta}_{t,2}, \hat{\phi}_{t,2} = \arg\max_{\theta_{t,2},\phi_{t,2}} \left|\hat{\mathbf{h}}_2(\boldsymbol{\alpha}_y(\theta_{t,2}, \phi_{t,2}) \otimes \boldsymbol{\alpha}_z(\phi_{t,2}))\right| θ ^ t , 2 , ϕ ^ t , 2 = arg max θ t , 2 , ϕ t , 2 h ^ 2 ( α y ( θ t , 2 , ϕ t , 2 ) ⊗ α z ( ϕ t , 2 ))
تعيين تقديرات الزاوية إلى الموقع ثلاثي الأبعاد بناءً على مبدأ المربعات الصغرى:
p ^ D = ( B B + B R + B U ) − 1 ( B B p B + B R p R + B U p U ) \hat{\mathbf{p}}_D = (\mathbf{B}_B + \mathbf{B}_R + \mathbf{B}_U)^{-1}(\mathbf{B}_B\mathbf{p}_B + \mathbf{B}_R\mathbf{p}_R + \mathbf{B}_U\mathbf{p}_U) p ^ D = ( B B + B R + B U ) − 1 ( B B p B + B R p R + B U p U )
الاستفادة من المسارات المتعددة : الاستفادة المتزامنة من مسارات BS-الطائرة-UE و BS-RIS-الطائرة-UEتصميم الحزم : ضمان تعامد حزم المحطة الأساسية، تقليل تداخل المسار المباشرخوارزمية منخفضة التعقيد : طريقة CGD لها تعقيد حسابي أقل مقارنة بتقليل الحد الأدنى للنورم النوويحدود نظرية : توفير تحليل معلومات فيشر و CRLB كمعايير أداءتكوين الهوائيات : M B = 8 × 8 M_B = 8 \times 8 M B = 8 × 8 , M U = 4 × 4 M_U = 4 \times 4 M U = 4 × 4 , M R = 6 × 6 M_R = 6 \times 6 M R = 6 × 6 مواقع العقد :
BS: (0, 0, 26) RIS: (0, 0.5, 25.5) الطائرة: (3, 3, 30) UE: (2, 2, 24) معاملات القناة : مؤشر فقدان المسار Γ = 2 \Gamma = 2 Γ = 2 ، عرض النطاق B = 20 B = 20 B = 20 MHzتباين الضوضاء : σ 2 = − 101 \sigma^2 = -101 σ 2 = − 101 dBmاستخدام جذر متوسط الخطأ التربيعي (RMSE) كمؤشر تقييم رئيسي:
RMSE = E { ( p D − p ^ D ) T ( p D − p ^ D ) } \text{RMSE} = \sqrt{E\{(\mathbf{p}_D - \hat{\mathbf{p}}_D)^T(\mathbf{p}_D - \hat{\mathbf{p}}_D)\}} RMSE = E {( p D − p ^ D ) T ( p D − p ^ D )}
مع RIS مقابل بدون RIS : التحقق من مكاسب أداء RISCRLB النظري : بمثابة معيار الحد الأعلى للأداءتُظهر نتائج المحاكاة في نطاق SNR من -10dB إلى 10dB:
نظام مع RIS : يُظهر مزايا كبيرة في منطقة SNR المنخفضةفجوة الأداء : عند SNR -10dB، يقلل نظام RIS RMSE بحوالي رتبة واحدة من حيث الحجم مقارنة بنظام بدون RISالقرب من النظرية : أداء الخوارزمية المقترحة قريبة من حد CRLB النظري السفليتأثير عدد فترات التدريب المختلفة K:
كلما زادت قيمة K، زادت دقة التحديد الموضع توجد نقطة تشبع في الأداء، والعائد الهامشي لزيادة K يتناقص التأثير الكبير لمقطع الرادار العاكس على الأداء:
ζ = 0.5 \zeta = 0.5 ζ = 0.5 : أداء نسبية سيئةζ = 1.0 \zeta = 1.0 ζ = 1.0 : معيار الأداء القياسيζ = 2.0 \zeta = 2.0 ζ = 2.0 : تحسن كبير في الأداءميزة SNR المنخفضة : مكاسب أداء RIS الأكثر أهمية في ظروف نسبة الإشارة إلى الضوضاء المنخفضةحساسية RCS : مقطع الرادار العاكس للطائرة هو عامل رئيسي يؤثر على دقة التحديد الموضعفعالية الخوارزمية : تحقق خوارزمية CGD أداءً قريبة من الحد النظري بتعقيد منخفضتقنيات التحديد الموضع 5G/6G : طرق التحديد الموضع بناءً على موجات ميلليمترية و MIMO واسع النطاقدمج الرادار والاستشعار : الكشف والتحديد الموضع للأهداف في أنظمة ISACالاتصالات المدعومة بـ RIS : تطبيقات الأسطح الذكية القابلة لإعادة التشكيل في الاتصالات اللاسلكيةمقارنة بالرادار التقليدي : الاستفادة من البنية التحتية الخلوية الموجودة، تقليل تكاليف النشرمقارنة بالتحديد الموضع النشط : قابل للتطبيق على الأهداف غير التعاونية، توسيع نطاق التطبيقاتمقارنة بأبحاث RIS الموجودة : أول تطبيق لـ RIS في سيناريو تحديد موقع الطائرة السلبييحسن RIS بشكل كبير أداء تحديد موقع الطائرة السلبي في نظام موجات ميلليمترية تحقق خوارزمية CGD المقترحة أداءً ممتازة بتعقيد منخفض معاملات النظام (RCS، النفقات التدريبية) لها تأثير مهم على الأداء افتراضات القناة : تأخذ في الاعتبار فقط مسارات الرؤية المباشرة (LoS)، قد توجد انتشارات متعددة المسارات في البيئات الفعليةالسيناريو الثابت : لم تأخذ في الاعتبار تأثير حركة الطائرة على دقة التحديد الموضعالمزامنة المثالية : افتراض المزامنة الزمنية والترددية المثالية للنظامالمواقع المعروفة : افتراض أن مواقع RIS و UE معروفة بدقةتحسين توزيع الطاقة : دراسة استراتيجية توزيع الطاقة المثلى بين حزم المحطة الأساسية المتعددةالتتبع الديناميكي : توسيع نطاق تتبع مسار الطائرات بدون طيار المتحركةمعالجة المسارات المتعددة : تحديد موقع قوي في بيئات الانتشار غير المباشرالنشر الفعلي : التحقق من أداء الخوارزمية في البيئات الحقيقيةابتكار قوي : أول تطبيق لتقنية RIS في تحديد موقع الطائرة السلبي، بابتكار قوياكتمال نظري : توفير تحليل كامل لمعلومات فيشر واشتقاق CRLBخوارزمية عملية : خوارزمية CGD منخفضة التعقيد، سهلة التطبيق العمليتحليل شامل : تحليل منهجي لتأثير المعاملات الرئيسية على الأداءالتحقق الكافي : النتائج النظرية والمحاكاة تؤكد بعضها البعضافتراضات مثالية : قد لا تنطبق افتراضات قناة LoS على البيئات الفعليةسيناريوهات محدودة : تنطبق بشكل أساسي على تحديد موقع الطائرات العالية في المساحات المفتوحةقيود التجربة : عدم وجود التحقق من الاختبارات الفعلية، فقط نتائج المحاكاةمعالجة التداخل : اعتبار غير كافٍ لقمع التداخل في بيئات كهرومغناطيسية معقدةالقيمة الأكاديمية : توفير أفكار جديدة لأبحاث RIS في تطبيقات الاستشعارالآفاق العملية : إمكانيات تطبيق في المدن الذكية والمراقبة الأمنية وغيرهادفع التكنولوجيا : تعزيز تطوير وتوحيد تقنية ISACأنظمة مراقبة الطائرات بدون طيار على أسطح المباني العالية في المدن شبكات الكشف عن الطائرات بدون طيار حول المطارات حماية المجال الجوي للمرافق المهمة تطبيقات الخدمات المضافة لشبكات 5G/6G تستشهد الورقة بأعمال مهمة متعددة في المجالات ذات الصلة، بما في ذلك أبحاث متقدمة في تقنيات تحديد موقع 5G، أنظمة ISAC، اتصالات RIS وغيرها، مما يوفر أساساً نظرياً قوياً لهذا العمل.