تضاؤل ​​الدقة

إن تخفيف الدقة ( DOP )، أو التخفيف الهندسي للدقة ( GDOP )، هو مصطلح يستخدم في الملاحة عبر الأقمار الصناعية وهندسة الجيوماتكس لتحديد انتشار الخطأ كتأثير رياضي لهندسة قمر الملاحة على دقة قياس الموقع.

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

مقدمة

نشأ مفهوم تدهور الدقة (DOP) مع مستخدمي نظام الملاحة لوران-سي . [ 1 ] وتتلخص فكرة تدهور الدقة الهندسي في توضيح كيفية تأثير أخطاء القياس على تقدير الحالة النهائية. ويمكن تعريف ذلك على النحو التالي: [ 2 ]

GDOP=Δ(موقع الإخراج)Δ(البيانات المقاسة){\displaystyle \operatorname {GDOP} ={\frac {\Delta ({\text{موقع الإخراج}})}{\Delta ({\text{البيانات المقاسة}})}}}

من الناحية النظرية، يمكنك أن تتخيل هندسياً الأخطاء التي تحدث في القياس والتي تؤدي إلىΔ(البيانات المقاسة){\displaystyle \Delta ({\text{البيانات المقاسة}})}تغيير المصطلح. من الناحية المثالية، لن تؤدي التغييرات الطفيفة في البيانات المقاسة إلى تغييرات كبيرة في موقع الإخراج. وعلى النقيض من ذلك، يكون الحل شديد الحساسية لأخطاء القياس. يوضح الشكل على اليمين تفسير هذه الصيغة، مبينًا سيناريوهين محتملين: أحدهما ذو دقة مقبولة والآخر ذو دقة ضعيفة.

مع الانتشار الواسع لأنظمة الملاحة عبر الأقمار الصناعية ، أصبح هذا المصطلح أكثر استخدامًا. وبإهمال تأثيرات الغلاف الأيوني [ 3 ] والغلاف الجوي السفلي [ 4 ] ، فإن إشارة الأقمار الصناعية للملاحة تتمتع بدقة ثابتة. ولذلك، تلعب الهندسة النسبية بين القمر الصناعي وجهاز الاستقبال دورًا رئيسيًا في تحديد دقة المواقع والأوقات المُقدَّرة. ونظرًا للهندسة النسبية لأي قمر صناعي بالنسبة لجهاز الاستقبال، فإن دقة المدى الزائف للقمر الصناعي تُترجم إلى مُكوِّن مُقابل في كل بُعد من أبعاد الموقع الأربعة التي يقيسها جهاز الاستقبال (أي،x{\displaystyle x}،y{\displaystyle y}،z{\displaystyle z}، وت{\displaystyle t}تعتمد دقة تحديد الموقع في كل بُعد من أبعاد قياسات جهاز الاستقبال على دقة إشارات عدة أقمار صناعية، وذلك وفقًا لموقعها النسبي. فعندما تكون أقمار الملاحة المرئية متقاربة في السماء، يُقال إن التداخل ضعيف وقيمة DOP عالية؛ أما عندما تكون متباعدة، فيُقال إن التداخل قوي وقيمة DOP منخفضة. لنفترض وجود حلقتين متداخلتين، أو ما يُعرف بالحلقات الدائرية ، لهما مركزان مختلفان. إذا تداخلتا بزاوية قائمة، فإن أقصى مدى للتداخل يكون أصغر بكثير مما لو تداخلتا بشكل شبه متوازٍ. وبالتالي، فإن انخفاض قيمة DOP يدل على دقة تحديد موقع أفضل نظرًا لزيادة المسافة الزاوية بين الأقمار الصناعية المستخدمة لحساب موقع الوحدة. ومن العوامل الأخرى التي قد تزيد من قيمة DOP الفعالة وجود عوائق مثل الجبال أو المباني القريبة.

يمكن التعبير عن DOP بعدد من القياسات المنفصلة:

HDOP
التخفيف الأفقي للدقة
VDOP
التخفيف الرأسي للدقة
PDOP
تخفيف دقة الموضع (ثلاثي الأبعاد)
TDOP
تدهور الدقة مع مرور الوقت
GDOP
التخفيف الهندسي للدقة

تُستنتج هذه القيم رياضياً من مواقع الأقمار الصناعية القابلة للاستخدام. وتتيح أجهزة استقبال الإشارة عرض هذه المواقع ( سكايبلوت ) بالإضافة إلى قيم DOP.

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

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

تفسير

قيمة DOPالتقييم [ 5 ]وصف
< 1مثاليأعلى مستوى ممكن من الثقة يُستخدم للتطبيقات التي تتطلب أعلى دقة ممكنة في جميع الأوقات.
1-2ممتازعند هذا المستوى من الثقة، تعتبر القياسات الموضعية دقيقة بما يكفي لتلبية جميع التطبيقات باستثناء التطبيقات الأكثر حساسية.
2-5جيديمثل هذا المستوى الحد الأدنى المناسب لاتخاذ قرارات دقيقة. ويمكن استخدام قياسات الموقع لتقديم اقتراحات موثوقة للملاحة أثناء الرحلة للمستخدم.
5-10معتدليمكن استخدام قياسات الموقع لإجراء الحسابات، ولكن لا يزال بالإمكان تحسين جودة التثبيت. يُنصح برؤية السماء بشكل أوضح.
10-20عدليمثل هذا مستوى ثقة منخفضًا. ينبغي تجاهل قياسات الموقع أو استخدامها فقط للإشارة إلى تقدير تقريبي للغاية للموقع الحالي.
> 20فقيرفي هذا المستوى، ينبغي تجاهل القياسات.

عوامل DOP هي دوال للعناصر القطرية لمصفوفة التغاير للمعلمات، معبر عنها إما في إطار جيوديسي عالمي أو محلي.

حساب

كخطوة أولى في حساب درجة الاحتمالية، [ 5 ] ضع في اعتبارك متجهات الوحدة من جهاز الاستقبال إلى القمر الصناعيأنا{\displaystyle i}: (xأنا-xRأنا،yأنا-yRأنا،zأنا-zRأنا)،Rأنا=(xأنا-x)2+(yأنا-y)2+(zأنا-z)2\begin{aligned}\left(\frac{x_{i}-x}{R_{i}},\frac{y_{i}-y}{R_{i}},\frac{z_{i}-z}{R_{i}}\right),\R_{i}=\sqrt{(x_{i}-x)^2+(y_{i}-y)^2+(z_{i}-z)^2}}\end{aligned}}}

أينx،y،z{\displaystyle x,y,z}يشير إلى موضع جهاز الاستقبال وxأنا،yأنا،zأنا{\displaystyle x_{i},y_{i},z_{i}}لنرمز إلى موقع القمر الصناعي i. صغ المصفوفة A، والتي (بالنسبة لأربع معادلات متبقية لقياس المدى الزائف) هي: أ=[x1-xR1y1-yR1z1-zR11x2-xR2y2-yR2z2-zR21x3-xR3y3-yR3z3-zR31x4-xR4y4-yR4z4-zR41]{\displaystyle A={\begin{bmatrix}{\frac {x_{1}-x}{R_{1}}}&{\frac {y_{1}-y}{R_{1}}}&{\frac {z_{1}-z}{R_{1}}}&1\\{\frac {x_{2}-x}{R_{2}}}&{\frac {y_{2}-y}{R_{2}}}&{\frac {z_{2}-z}{R_{2}}}&1\\{\frac {x_{3}-x}{R_{3}}}&{\frac {y_{3}-y}{R_{3}}}&{\frac {z_{3}-z}{R_{3}}}&1\\{\frac {x_{4}-x}{R_{4}}}&{\frac {y_{4}-y}{R_{4}}}&{\frac {z_{4}-z}{R_{4}}}&1\end{bmatrix}}}

العناصر الثلاثة الأولى من كل صف في المصفوفة A هي مركبات متجه وحدة من جهاز الاستقبال إلى القمر الصناعي المحدد. يشير العنصر الأخير من كل صف إلى المشتق الجزئي للمدى الزائف بالنسبة لانحياز ساعة جهاز الاستقبال. صِغ المصفوفة Q على أنها مصفوفة التغاير الناتجة عن مصفوفة التوزيع الطبيعي باستخدام طريقة المربعات الصغرى . سؤال=(أتيأ)-1{\displaystyle Q=\left(A^{\mathsf {T}}A\right)^{-1}}

على العموم: سؤال=(جxتي(جدجدجدتي)-1جx)-1{\displaystyle Q=\left(J_{\mathsf {x}}^{\mathsf {T}}\left(J_{d}C_{d}J_{d}^{\mathsf {T}}\right)^{-1}J_{x}\right)^{-1}}

أينجx{\displaystyle J_{x}}يمثل جاكوبيان معادلات البواقي لقياسات المستشعروأنا(x_،د_)=0{\displaystyle f_{i}\left({\underline {x}},{\underline {d}}\right)=0}فيما يتعلق بالمجهولات،x_{\displaystyle {\underline {x}}}؛جد{\displaystyle J_{d}}يمثل جاكوبيان معادلات البواقي لقياسات المستشعر بالنسبة للكميات المقاسةد_{\displaystyle {\underline {d}}}، وجد{\displaystyle C_{d}}هي مصفوفة الارتباط للضوضاء في الكميات المقاسة.

بالنسبة للحالة السابقة المكونة من 4 معادلات متبقية لقياس المدى:x_=(x،y،z،τ)تي{\displaystyle {\underline {x}}=(x,y,z,\tau )^{\mathsf {T}}}،د_=(τ1،τ2،τ3،τ4)تي{\displaystyle {\underline {d}}=\left(\tau _{1},\tau _{2},\tau _{3},\tau _{4}\right)^{\mathsf {T}}}،τ=جت{\displaystyle \tau =ct}،τأنا=جتأنا{\displaystyle \tau _{i}=ct_{i}}،Rأنا=|τأنا-τ|=(τأنا-τ)2{\displaystyle R_{i}=|\tau _{i}-\tau |={\sqrt {(\tau _{i}-\tau )^{2}}}}،وأنا(x_،د_)=(xأنا-x)2+(yأنا-y)2+(zأنا-z)2-(τأنا-τ)2{\displaystyle f_{i}\left({\underline {x}},{\underline {d}}\right)={\sqrt {(x_{i}-x)^{2}+(y_{i}-y)^{2}+(z_{i}-z)^{2}}}-{\sqrt {(\tau _{i}-\tau )^{2}}}}،جx=أ{\displaystyle J_{x}=A}،جد=-أنا{\displaystyle J_{d}=-I}وضوضاء القياس لمختلفτأنا{\displaystyle \tau _{i}}وقد افترض أنها مستقلة مما يجعلجد=أنا{\displaystyle C_{d}=I}.

تنشأ هذه الصيغة لـ Q من تطبيق أفضل تقدير خطي غير متحيز على نسخة خطية من معادلات البواقي لقياسات المستشعر حول الحل الحاليΔx_=-سؤال*(جxتي(جدجدجدتي)-1و){\displaystyle \Delta {\underline {x}}=-Q*\left(J_{x}^{\mathsf {T}}\left(J_{d}C_{d}J_{d}^{\mathsf {T}}\right)^{-1}f\right)}باستثناء حالة اللون الأزرقجد{\displaystyle C_{d}}هي مصفوفة تباين الضوضاء وليست مصفوفة ارتباط الضوضاء المستخدمة في DOP، والسبب في قيام DOP بهذا الاستبدال هو الحصول على خطأ نسبي . عندماجد{\displaystyle C_{d}}هي مصفوفة تباين الضوضاء،سؤال{\displaystyle Q}يمثل هذا تقديرًا لمصفوفة التغاير للضوضاء في المجاهيل الناتجة عن الضوضاء في الكميات المقاسة. وهو التقدير الذي تم الحصول عليه باستخدام تقنية تحديد عدم اليقين من الدرجة الأولى للعزم الثاني (FOSM)، والتي كانت من أحدث التقنيات في ثمانينيات القرن الماضي. ولكي تكون نظرية FOSM قابلة للتطبيق بدقة، يجب أن تكون توزيعات ضوضاء الإدخال غاوسية، أو أن تكون الانحرافات المعيارية لضوضاء القياس صغيرة نسبيًا مقارنةً بمعدل التغير في المخرجات بالقرب من الحل. وفي هذا السياق، عادةً ما يكون المعيار الثاني هو المُتحقق.

يتم حساب هذا (أي بالنسبة لمعادلات البواقي الأربعة لقياس وقت الوصول/المدى) وفقًا لـ [6] حيث مصفوفة الترجيح، P=(جدجدجدتي)-1{\displaystyle P=\left(J_{d}C_{d}J_{d}^{\mathsf {T}}\right)^{-1}}يحدث أن يتبسط إلى مصفوفة الوحدة .

لاحظ أن P لا تتبسط إلا إلى مصفوفة الوحدة لأن جميع معادلات البواقي لقياسات المستشعر هي معادلات زمن الوصول (المدى الزائف). في حالات أخرى، على سبيل المثال عند محاولة تحديد موقع شخص يبث على تردد استغاثة دولي ،P{\displaystyle P}لن يتبسط إلى مصفوفة الوحدة ، وفي هذه الحالة سيكون هناك مكون "تردد DOP" أو FDOP إما بالإضافة إلى مكون TDOP أو بدلاً منه. (بخصوص "بدلاً من مكون TDOP": بما أن ساعات أقمار برنامج Cospas-Sarsat الدولي القديمة في المدار الأرضي المنخفض أقل دقة بكثير من ساعات نظام تحديد المواقع العالمي (GPS)، فإن تجاهل قياساتها الزمنية سيزيد في الواقع من دقة حل تحديد الموقع الجغرافي).

عناصرسؤال{\displaystyle Q}يتم تصنيفها على النحو التالي:

سؤال=[σx2σxyσxzσxتσxyσy2σyzσyتσxzσyzσz2σzتσxتσyتσzتσت2]{\displaystyle Q={\begin{bmatrix}\sigma _{x}^{2}&\sigma _{xy}&\sigma _{xz}&\sigma _{xt}\\\sigma _{xy}&\sigma _{y}^{2}&\sigma _{yz}&\sigma _{yt}\\\sigma _{xz}&\sigma _{yz}&\sigma _{z}^{2}&\sigma _{zt}\\\sigma _{xt}&\sigma _{yt}&\sigma _{zt}&\sigma _{t}^{2}\end{bmatrix}}}

يتم إعطاء PDOP و TDOP و GDOP بواسطة: [ 6 ]PDOP=σx2+σy2+σz2TDOP=σت2GDOP=PDOP2+TDOP2=trسؤال{\displaystyle {\begin{aligned}\operatorname {PDOP} &={\sqrt {\sigma _{x}^{2}+\sigma _{y}^{2}+\sigma _{z}^{2}}}\\\operatorname {TDOP} &={\sqrt {\sigma _{t}^{2}}}\\\operatorname {GDOP} &={\sqrt {\operatorname {PDOP} ^{2}+\operatorname {TDOP} ^{2}}}\\&={\sqrt {\operatorname {tr} Q}}\\\end{aligned}}}

لاحظ أن GDOP هو الجذر التربيعي لأثرسؤال{\displaystyle Q}مصفوفة.

التخفيف الأفقي والرأسي للدقة، HDOP=σن2+σهـ2VDOP=σu2،{\displaystyle {\begin{aligned}\operatorname {HDOP} &={\sqrt {\sigma _{n}^{2}+\sigma _{e}^{2}}}\\\operatorname {VDOP} &={\sqrt {\sigma _{u}^{2}}}\end{aligned}},} يعتمد كلاهما على نظام الإحداثيات المستخدم. لكي يتوافق مع نظام الإحداثيات المحلي شرق-شمال-أعلى ،

EDOP^2 xxx x NDOP^2 xx xx VDOP^2 x xxx TDOP^2

والتخفيفات المشتقة: GDOP=EDOP2+NDOP2+VDOP2+TDOP2HDOP=EDOP2+NDOP2PDOP=EDOP2+NDOP2+VDOP2{\displaystyle {\begin{aligned}\operatorname {GDOP} &={\sqrt {\operatorname {EDOP} ^{2}+\operatorname {NDOP} ^{2}+\operatorname {VDOP} ^{2}+\operatorname {TDOP} ^{2}}}\\\operatorname {HDOP} &={\sqrt {\operatorname {EDOP} ^{2}+\operatorname {NDOP} ^{2}}}\\\operatorname {PDOP} &={\sqrt {\operatorname {EDOP} ^{2}+\operatorname {NDOP} ^{2}+\operatorname {VDOP} ^{2}}}\end{aligned}}}

انظر أيضاً

مراجع

  1. ريتشارد ب. لانغلي (مايو 1999). "تضاؤل ​​الدقة" (ملف PDF) . عالم نظام تحديد المواقع العالمي (GPS ). مؤرشف (ملف PDF) من الأصل بتاريخ 4 أكتوبر 2011. تم الاطلاع عليه بتاريخ 12 أكتوبر 2011 .
  2. دوديك، غريغوري ؛ جينكين، مايكل (2000). المبادئ الحسابية للروبوتات المتنقلة . مطبعة جامعة كامبريدج . ISBN 0-521-56876-5.
  3. بول كينتنر، جامعة كورنيل؛ تود همفريز، جامعة تكساس في أوستن؛ جوانا هينكس، جامعة كورنيل (يوليو-أغسطس 2009). "نظام الملاحة العالمي عبر الأقمار الصناعية والوميض الأيوني: كيفية النجاة من ذروة النشاط الشمسي القادمة" . مجلة Inside GNSS . مؤرشف من الأصل بتاريخ 6 نوفمبر 2011. تاريخ الاسترجاع: 12 أكتوبر 2011 .
  4. "أخطاء نظام تحديد المواقع العالمي (شرح تريمبل)" . مؤرشف من الأصل بتاريخ 7 مارس 2016. تم الاطلاع عليه بتاريخ 8 فبراير 2016 .
  5. 1 2 إيسيك، أوغوز كاغان؛ هونغ، جوهيون؛ بيترونين، إيفان؛ تسوردوس، أنطونيوس (25 أغسطس 2020). "تحليل سلامة الملاحة القائمة على نظام تحديد المواقع العالمي (GPS) للطائرات بدون طيار في البيئة الحضرية" . الروبوتات . 9 (3): 66. doi : 10.3390/robotics9030066 .
  6. القسم 1.4.9 من مبادئ تحديد المواقع عبر الأقمار الصناعية .

للمزيد من القراءة