1

Есть точка А(ш,д) и точка Б(ш,д). Мы знаем расстояние между ними (с единичным радиусом), вычисленное по формуле

arccos( sin(д1)*sin(д2) + cos(д1)*cos(д2)*cos(ш2-ш1)).

Далее я сдвигаю точку А. Теперь, чтобы вычислить новое расстояние, нужно снова просчитывать 6 тригонометрических функций. И я вот думаю, а не проще ли будут расчеты, если новое расстояние считать, зная старые координаты и новые координаты? Помогите составить формулу

2
  • Я думаю, вы экономите на спичках. Проще формула не будет (она не проще даже на плоскости, а уж на сфере и подавно), а тригонометрические функции типа косинуса — одна инструкция процессора.
    – VladD
    Commented 12 окт. 2015 в 8:58
  • 1
    @VladD, я тоже так думаю. Но вдруг бы это сократило формулу с 6 функций до 3х. Далее умножаем на количество вычисленных расстояний... В общем попытаться стоило )
    – iRumba
    Commented 12 окт. 2015 в 9:17

1 ответ 1

2

A(x1,y1,z1), B(x2,y2,z2), и dA(dx, dy, dz)

!!A и B вычисляются переводом сферических координат в декартову!!

старое расстояние R1=Корень(квадрат(x1-x2)+квадрат(y1-y2)+квадрат(z1-z2))

новое расстояние R2=Корень(квадрат(x1-x2+dx)+квадрат(y1-y2+dy)+квадрат(z1-z2+dz))

разница расстояний dR = R2-R1;

прочитайте про перевод из сферических координат в декартову тут, и будет Вам счастье

По другому никак

UPD:

Есть вектор A(ϴa,φa), Вектор B(ϴb,φb), и вектор A2(ϴa+dϴ,φa+dφ)

они же в декартовой системе A(Xa,Ya,Za) B(Xb,Yb,Zb),

Где Xi=sin(ϴi)*cos(φi), Yi=sin(ϴi)*sin(φi), Zi=cos(ϴi)

Найдем угол между ними αab=ArcCos(ab'/(a'*b'))

где ab'=XaXb+YaYb+Za*Zb, a'=Корень(Xa^2+Ya^2+Za^2), b'=Корень(Xb^2+Yb^2+Zb^2)

и длинну дуги Lab=αab*R

т.е. примерно так:

public double GetSphereSectorLen(double ϴ1,double φ1,double ϴ2,double φ2, double R)
{
  double X1=Math.Sin(ϴ1)*Math.Cos(φ1);
  double Y1=Math.Sin(ϴ1)*Math.Sin(φ1);
  double Z1=Math.Cos(ϴ1);
  double X2=Math.Sin(ϴ2)*Math.Cos(φ2);
  double Y2=Math.Sin(ϴ2)*Math.Sin(φ2);
  double Z2=Math.Cos(ϴ2);
  double ab_=X1*X2+Y1*Y2+Z1*Z2;
  double a_=Math.Sqrt(Math.Pow(X1,2)+Math.Pow(Y1,2)+Math.Pow(Z1,2));
  double b_=Math.Sqrt(Math.Pow(X2,2)+Math.Pow(Y2,2)+Math.Pow(Z2,2));
  double αab=Math.Acos(ab_/(a_*b_));
  return αab*R;
}

double LAB=GetSphereSectorLen(ϴa,φa,ϴb,φb);
double LA2B=GetSphereSectorLen(ϴa+dϴ,φa+dφ,ϴb,φb);
double dL = LAB-LA2B;
8
  • В данной задаче уже нет x и y и расстояние вычисляется иначе
    – iRumba
    Commented 9 окт. 2015 в 10:53
  • ответил в посте Commented 9 окт. 2015 в 10:59
  • А вообще вы зря взялись за такую задачу с Вашими знаниями линейной алгебры (судя по предыдущему посту) Commented 9 окт. 2015 в 11:08
  • Меня никто не спрашивал, хочу я за нее браться или нет. Заодно заполню дыры в знаниях. По крайней мере теперь знаю для чего все это преподавалось, а значит будет веселее ))
    – iRumba
    Commented 9 окт. 2015 в 11:10
  • 1
    дык надож описывать задачу полностью, и перемещение у вас, звените будет не dX, dY, dZ, а dθ, dφ, и искать вам надо дугу окружности, а не расстояние между точками, что является корнем сумм квадратов разностей координат. Формулируйте задачу правильно Commented 9 окт. 2015 в 11:28

Ваш ответ

Нажимая «Отправить ответ», вы соглашаетесь с условиями пользования и подтверждаете, что прочитали политику конфиденциальности.

Всё ещё ищете ответ? Посмотрите другие вопросы с метками или задайте свой вопрос.