Если у вас есть векторы положения A=(xa,ya) и B=(xb,yb) и соответствующие скорости Va и Vb, вы пытаетесь решить A+s*Va=B+t*Vb. Это решение для любой системы координат, но сначала вы должны выбрать систему координат.
Давайте сначала решим это для евклидова пространства ...
Запишите это для каждого компонента x и y. Теперь у вас есть два уравнения одновременно для двух переменных, и вы можете решить как s, так и t. Они столкнутся, если s==t (или «близко»). Не упустите деление на ноль, когда Va и Vb параллельны, и числовую нестабильность.
Для сферического решения рассмотрим две точки, движущиеся по окружностям окружностей, уравнения имеют одинаковую форму, хотя A и B будут углами, а Va и Vb - угловыми скоростями. Чтобы получить круги, в любой момент можно мгновенно вычислить евклидовы Va и Vb, взятые с помощью евклидовых A и B (при условии, что Земля центрирована в 0), чтобы сказать вам, в какой плоскости вы работаете, спроецируйте в эту плоскость, чтобы получить двумерную задачу для каждой плоскости отдельно.