معادلات مشاهداتی شبه‌فاصلهتعیین موقعیت تک‌نقطه‌ای (SPP)

در هر مسئله اندازه‌گیری، یک معادله مشاهداتی (Observation Equation) رابطه‌ای ریاضی بین اندازه‌گیری‌ها (مقادیر معلوم) و مجهولات (مقادیری که می‌خواهیم تعیین کنیم) برقرار می‌کند. بدون این معادله، هیچ الگوریتم تعیین موقعیتی قابل طراحی نیست.

در تعیین موقعیت تک‌نقطه‌ای (SPP — Single Point Positioning):

• اندازه‌گیری: شبه‌فاصله (Pseudorange) — فاصله اندازه‌گیری‌شده بین گیرنده و ماهواره
• مجهولات: مختصات سه‌بعدی گیرنده (X, Y, Z) + خطای ساعت گیرنده (δt_r)

چرا درک معادله مهم است؟ وقتی نرم‌افزار تعیین موقعیت جوابی با دقت ۵ متر می‌دهد، بدون درک معادله مشاهداتی نمی‌توانید بفهمید چرا دقت بهتر نشده، کدام منبع خطا غالب است، و چگونه می‌توان نتیجه را بهبود داد. معادله مشاهداتی نقشه راه شماست.

معادله اصلی شبه‌فاصله که پایه تمام محاسبات SPP است:

هر جمله در این معادله معنای فیزیکی مشخصی دارد. بیایید یک‌به‌یک بررسی کنیم:

اندازه معمول هر جمله خطا:

معادله شبه‌فاصله غیرخطی است — زیرا فاصله هندسی ρ شامل جذر مجموع مربعات مختصات مجهول است. برای حل آن باید معادله را خطی‌سازی کنیم.

بسط تیلور مرتبه اول فاصله هندسی:

مشتقات جزئی (که به آن‌ها کسینوس‌های هادی یا Direction Cosines می‌گویند):

پس از تصحیح خطاهای مدل‌شده (ساعت ماهواره، یونوسفر، تروپوسفر)، معادله خطی‌شده به این شکل در می‌آید:

(اختلاف مشاهده و محاسبه)

با n ماهواره، n معادله خواهیم داشت. این سیستم معادلات را به شکل ماتریسی می‌نویسیم:

      | a_x1   a_y1   a_z1   1 |       | δX    |       | ΔP1 |
      | a_x2   a_y2   a_z2   1 |       | δY    |       | ΔP2 |
A =   | a_x3   a_y3   a_z3   1 |   x = | δZ    |   l = | ΔP3 |
      | ...    ...    ...    .. |       | c·δt_r|       | ... |
      | a_xn   a_yn   a_zn   1 |                       | ΔPn |

۴ مجهول داریم: δX, δY, δZ, c·δt_r — بنابراین حداقل ۴ ماهواره نیاز است. با بیش از ۴ ماهواره، سیستم فوق‌تعیین (Over-determined) می‌شود و از روش کمترین مربعات استفاده می‌کنیم.

با n ماهواره (که n ≥ 4) و ۴ مجهول، سیستم معادلات فوق‌تعیین است. روش کمترین مربعات (Least Squares) جوابی پیدا می‌کند که مجموع مربعات باقی‌مانده‌ها را مینیمم کند.

ماتریس وزن W یک ماتریس قطری است که اهمیت نسبی هر مشاهده را مشخص می‌کند:

فرآیند حل تکراری (Iterative) است:

بیایید با یک مثال ساده‌شده، تمام مراحل را عملی کنیم. فرض کنید ۴ ماهواره GPS در دید هستند و موقعیت تقریبی گیرنده را داریم.

Sat   X (m)             Y (m)             Z (m)
S1    15600000.0        7540000.0         20140000.0
S2    18760000.0        2750000.0         18610000.0
S3    17610000.0        14630000.0        13480000.0
S4    19170000.0        610000.0          18390000.0

X₀ = 3859000.0 m
Y₀ = 3170000.0 m
Z₀ = 3467000.0 m
c·δt₀ = 0.0 m (فرض اولیه)

مرحله ۱: محاسبه فاصله هندسی تقریبی ρ₀ برای هر ماهواره:

ρ₀_1 = √((15600000-3859000)² + (7540000-3170000)² + (20140000-3467000)²)
     = 20462187.9 m

ρ₀_2 = √((18760000-3859000)² + (2750000-3170000)² + (18610000-3467000)²)
     = 21114825.8 m

ρ₀_3 = √((17610000-3859000)² + (14630000-3170000)² + (13480000-3467000)²)
     = 18622262.5 m

ρ₀_4 = √((19170000-3859000)² + (610000-3170000)² + (18390000-3467000)²)
     = 21255338.1 m

مرحله ۲: محاسبه کسینوس‌های هادی و ساختن ماتریس A:

      |  a_x      a_y      a_z     1 |
A =   | -0.5737  -0.2136  -0.8152  1 |    (S1)
      | -0.7058  +0.0199  -0.7172  1 |    (S2)
      | -0.7382  -0.6155  -0.5378  1 |    (S3)
      | -0.7203  +0.1204  -0.7025  1 |    (S4)

مرحله ۳: فرض کنید شبه‌فاصله‌های اندازه‌گیری‌شده (پس از تصحیح ساعت ماهواره، یونوسفر و تروپوسفر) اینها باشند:

P1_corrected = 20462315.0 m     ΔP1 = 20462315.0 - 20462187.9 = +127.1 m
P2_corrected = 21114948.0 m     ΔP2 = 21114948.0 - 21114825.8 = +122.2 m
P3_corrected = 18622387.0 m     ΔP3 = 18622387.0 - 18622262.5 = +124.5 m
P4_corrected = 21255464.0 m     ΔP4 = 21255464.0 - 21255338.1 = +125.9 m

مرحله ۴: حل سیستم (بدون وزن‌دهی برای سادگی):

x = | δX     |   =   | -30.2 m  |
    | δY     |       | +18.7 m  |
    | δZ     |       | -42.5 m  |
    | c·δt_r |       | +95100.0 m | (≈ 317 μs clock error)

Updated Position:
X = 3859000.0 + (-30.2) = 3858969.8 m
Y = 3170000.0 + (+18.7) = 3170018.7 m
Z = 3467000.0 + (-42.5) = 3466957.5 m

تعیین موقعیت SPP با معادله شبه‌فاصله، ساده‌ترین و در دسترس‌ترین روش است — اما محدودیت‌های جدی دارد. دقت معمول آن ۲ تا ۵ متر است.

روش‌های دقیق‌تر چگونه این محدودیت‌ها را رفع می‌کنند؟