Определить координаты точки (x, y, z) в сферической системе координат.( С введением: Ветвления) На языке Си.

[Stanislav09]

Сообщение от type_Oleg

Stanislav09, и разберитесь со сферическими координатами.
На примере географических ( они правда не совсем сферические, но приблизительно можно считать так).
Например на полюсе долгота не определена , определена только широта фи=90 и расстояние до центра. Потому что там x=0 и y=0.
Ну и в центре Земли ни широта ни долгота не определены, только расстояние =0.

Вот с фи я разобрался...

#include<stdio.h>        // Программа по нахождению координат (x,y,z) в сферической системе координат //
#include<math.h>
#include<locale.h>

void main()
{
    setlocale(LC_ALL, "Rus");
    float x,y,z,p,fi,alpha;
    printf("Введите координату x:");    // Ввод координаты x //
    scanf("%f", &x);
    printf("Введите координату y:");    // Ввод координаты y //
    scanf("%f", &y);
    printf("Введите координату z:");    // Ввод координаты z //
    scanf("%f", &z);
    if (x > 0, y > 0, z > 0)
    {
        p = sqrt( pow(x,2) + pow(y,2) + pow(z,2) );      // Нахождение координаты p //
        fi = atan ( x / y );            // Нахождение координаты fi //
        alpha = acos( z / p );                          // Нахождение координаты alpha //
        printf("Координаты x,y,z в сферической системе координат p,fi,alpha равны: %f, %f, %f", p, fi, alpha);  // Вывод координат p,fi,alpha //
    }
    else
    {
        printf("Ошибка!!!");  // Вывод ошибки //
    }
}

[Stanislav09]

Сообщение от type_Oleg

if (x > 0, y > 0, z > 0)

а почему только в одном из 8 квадрантов ? По моему в пространстве x,y,z могут быть любые

А вот перед этим -

fi = x / sqrt( pow(x,2) + pow(y,2) );
alpha = acos( z / p );

как раз можно проверку знаменателя на 0. Вот тогда как раз сферические координаты и не определить.

Кстати, зачем эти pow(x,2), не проще x*x.

Да, и фи у вас неправильно. Арктангенс забыли. Или что-то напутали..

А вот это я не понял, что именно нужно сделать?

[digitalis]

Это звучит так:

if ( (x != 0) || ( y != 0) || ( z != 0) )

(достаточно, чтобы хоть одна координата была ненулевой, тогда неопределённости не будет)
Действительно, точка может быть в одном из 8 квадрантов октантов, так что x,y,z не обязательно положительные.

А вот это я не понял, что именно нужно сделать?

Дык это - к вышеприведённой строке - если все x,y,z = 0, то мы в центре сферы, и определение сфер. координат лишено смысла.

[type_Oleg]

Посмотрел ради интереса, как на сайте IAU SOFA сделано.
Вот -

void iauC2s(double p[3], double *theta, double *phi)
{
   double x, y, z, d2;
   x  = p[0];
   y  = p[1];
   z  = p[2];
   d2 = x*x + y*y;
   *theta = (d2 == 0.0) ? 0.0 : atan2(y, x);
   *phi = (z == 0.0) ? 0.0 : atan2(z, sqrt(d2));
   return;
}

У них там в векторах. Вам нужно оттуда всего 3 строчки

  d2 = x*x + y*y;
  theta = (d2 == 0.0) ? 0.0 : atan2(y, x);
  phi = (z == 0.0) ? 0.0 : atan2(z, sqrt(d2));

Плюс расстояние найти sqrt(d2+z*z)