UTM และ Remote ID

ตัวอย่างที่ 1 เส้นทางส่งเวชภัณฑ์สองเส้นที่ตัดกัน

เที่ยวบิน A บินตะวันออก 6 km ที่ความสูง 30–90 m เที่ยวบิน B บินขึ้นเหนือตัดผ่านกลางเส้นทาง A ที่ 50–90 m ทั้งคู่บิน 15 m/s แบ่งช่วงละ 1 km กว้างข้างละ 100 m เผื่อเวลา ±30 s

import math

def intent(start, end, alt, speed, depart, seg=1000, half_w=100, buf=30):
    (x0, y0), (x1, y1) = start, end
    length = math.dist(start, end)
    vols = []
    for i in range(math.ceil(length / seg)):
        a, b = i * seg / length, min((i + 1) * seg / length, 1)
        xa, ya = x0 + a * (x1 - x0), y0 + a * (y1 - y0)
        xb, yb = x0 + b * (x1 - x0), y0 + b * (y1 - y0)
        vols.append((min(xa, xb) - half_w, max(xa, xb) + half_w, min(ya, yb) - half_w, max(ya, yb) + half_w,
                     alt[0], alt[1], depart + a * length / speed - buf, depart + b * length / speed + buf))
    return vols

def overlap(p, q):
    return all(p[2 * k] < q[2 * k + 1] and q[2 * k] < p[2 * k + 1] for k in range(4))

A = intent((0, 0), (6000, 0), (30, 90), 15, depart=0)
for delay in range(0, 601, 30):
    B = intent((3000, -3000), (3000, 3000), (50, 90), 15, depart=delay)
    clashes = [(i + 1, j + 1) for i, p in enumerate(A) for j, q in enumerate(B) if overlap(p, q)]
    print(f"B departs +{delay:>3} s: {len(clashes)} conflicting segment pairs {clashes}")
    if not clashes:
        break

B departs +  0 s: 4 conflicting segment pairs [(3, 3), (3, 4), (4, 3), (4, 4)]
B departs + 30 s: 4 conflicting segment pairs [(3, 3), (3, 4), (4, 3), (4, 4)]
B departs + 60 s: 3 conflicting segment pairs [(3, 3), (4, 3), (4, 4)]
B departs + 90 s: 3 conflicting segment pairs [(3, 3), (4, 3), (4, 4)]
B departs +120 s: 3 conflicting segment pairs [(3, 3), (4, 3), (4, 4)]
B departs +150 s: 1 conflicting segment pairs [(4, 3)]
B departs +180 s: 1 conflicting segment pairs [(4, 3)]
B departs +210 s: 0 conflicting segment pairs []

ถ้า B ออกพร้อม A จะขัดแย้งสี่คู่บริเวณจุดตัด การเลื่อน B ออกไป 210 วินาทีทำให้ไม่มีคู่ใดทับกัน อีกทางหนึ่งคือแยกความสูงไม่ให้ทับกันเลย ซึ่งต้องอยู่ในเพดานที่ได้รับอนุญาต

ตัวอย่างที่ 2 ตรวจ log Remote ID ตามเกณฑ์ Part 89

ใช้เกณฑ์ของสหรัฐฯ เป็นตัวอย่างวิธีตรวจ ข้อมูลเป็น log สมมติ 20 ข้อความ

FT100 = 100 * 0.3048
meas = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12.4, 13.4, 14.4, 15.4, 16.4, 17.4, 18.4, 19.4]
lat = [0.3] * 20
lat[7] = 1.2
errs = [4, 6, 5, 8, 7, 12, 9, 6, 5, 11, 14, 10, 8, 7, 34, 9, 6, 5, 7, 8]

sent = [m + l for m, l in zip(meas, lat)]
gaps = [b - a for a, b in zip(sent, sent[1:])]
print(f"broadcasts {len(sent)}, largest gap {max(gaps):.1f} s, gaps over 1 s: {sum(g > 1.0 for g in gaps)}")
print(f"latency over 1 s: {sum(x > 1.0 for x in lat)}")
within = sum(e <= FT100 for e in errs) / len(errs)
print(f"position within {FT100:.2f} m: {within:.0%} (needs at least 95%)")

broadcasts 20, largest gap 1.9 s, gaps over 1 s: 2
latency over 1 s: 1
position within 30.48 m: 95% (needs at least 95%)

ความแม่นตำแหน่งผ่านพอดี แต่มีช่วงที่ไม่ได้ส่งเกิน 1 วินาทีสองช่วง และข้อความหนึ่งส่งช้าเกิน 1 วินาที ต้องหาสาเหตุ เช่น งานอื่นแย่งเวลาประมวลผลของโมดูลส่งสัญญาณ