Местоположение спутника

[xd]

Ну так делайте трассировку и ищите аномалии. Исключения перехватили?
Вы сейчас описываете какую-то непонятную ситуацию и просите участников форума включить режим телепата, читающего мысли Вашего телефона. Хотите конструктива - давайте нормальную обратную связь.

[INA]

Да нет... с режимом телепата я тоже знаком. 
Я просто сейчас хочу точно понять - если модуль SGP4 на выходе выдаёт шесть значений: 3 - это координаты спутника на момент расчёта и ещё 3 - это вектора скорости на тот же момент, то можно ли сделать прогноз, используя эти шесть значений без повторного обращения к модулю SGP4. Поскольку у меня софт начинает сбиваться именно при построении трека, т.к. в этом случае, чтоб отрисовать трек, софт в цикле отрабатывает модуль SGP4 около 90 раз. Вот на этом моменте софт и виснет. Поэтому, не соображая ничего в небесной механике, я и спросил: правильно ли будет, если я при запуске софта один раз в сутки сделаю обращение к модулю SGP4 и по его результатам - этим шести значениям - смогу ли далее, до перезапуска программы, прогнозировать и трек и положение спутника в следующие моменты времени? Я подумал, что если есть координаты и три вектора скорости для какого-то ближайшего момента времени, то возможно ли по ним сделать прогноз движения на последующие моменты времени.
Вот этом момент проясните мне, пожалуйста. Я действительно, не понимаю немного эти тонкости.
Модуль SGP4 очень сильно грузит вычислениями железо. И быть может я неправильно понял сам алгоритм работы подобных программ и сам лишний раз загружаю железо лишними вычислениями. 

[Sleepwalker]

По сишным реализациям помню, что там была отдельно инициализация и отдельно прогноз.

Вот по этому поводу подскажите, пожалуйста подробнее.
В первом своём посте я как раз об этом спросил: нужно ли каждую секунду обрабатывать модуль SGP4 для получения координат спутника или это нужно сделать один раз, получив его координаты на момент времени запуска программы, а потом по этим шести значениям делать прогноз?
По значениям, которые возвращает SGP4.

Конечно, надо запускать прогноз модели sgp для каждого момента времени. Просто первоначальная инициализация делает некоторый предрасчет, который не надо делать при прогнозе. А вообще сам прогноз там аналитический (для низких орбит) и время выполнения прогноза на любом шаге относительно первоначальной эпохи привязки TLE должно быть одинаковым.

[Sleepwalker]

Да нет... с режимом телепата я тоже знаком. 
Я просто сейчас хочу точно понять - если модуль SGP4 на выходе выдаёт шесть значений: 3 - это координаты спутника на момент расчёта и ещё 3 - это вектора скорости на тот же момент, то можно ли сделать прогноз, используя эти шесть значений без повторного обращения к модулю SGP4. Поскольку у меня софт начинает сбиваться именно при построении трека, т.к. в этом случае, чтоб отрисовать трек, софт в цикле отрабатывает модуль SGP4 около 90 раз. Вот на этом моменте софт и виснет. Поэтому, не соображая ничего в небесной механике, я и спросил: правильно ли будет, если я при запуске софта один раз в сутки сделаю обращение к модулю SGP4 и по его результатам - этим шести значениям - смогу ли далее, до перезапуска программы, прогнозировать и трек и положение спутника в следующие моменты времени? Я подумал, что если есть координаты и три вектора скорости для какого-то ближайшего момента времени, то возможно ли по ним сделать прогноз движения на последующие моменты времени.
Вот этом момент проясните мне, пожалуйста. Я действительно, не понимаю немного эти тонкости.
Модуль SGP4 очень сильно грузит вычислениями железо. И быть может я неправильно понял сам алгоритм работы подобных программ и сам лишний раз загружаю железо лишними вычислениями.

Правильно делать как у Вас сейчас, каждый раз вызывать прогноз SGP (вынеся конструктор модели за цикл). Виснет, видимо, из-за ошибки при кодировании.
Можно сделать по-хулигански. Получить xyz vx vy vx и дальше прогнозировать по Кеплеру, только это ничего в принципе не изменит, это будет такой-же прогноз, но орбита будет постоянно уплывать от модельных значения sgp.

[INA]

Правильно делать как у Вас сейчас, каждый раз вызывать прогноз SGP (вынеся конструктор модели за цикл). Виснет, видимо, из-за ошибки при кодировании.
Можно сделать по-хулигански. Получить xyz vx vy vx и дальше прогнозировать по Кеплеру, только это ничего в принципе не изменит, это будет такой-же прогноз, но орбита будет постоянно уплывать от модельных значения sgp.

Спасибо большое Вам за разъяснения!

Разобрался я в чём было дело. Софт заработал как положено.

Да, начальная инициализация у меня сделана при запуске софта, а далее вызывается только модуль SGP4 - он и вычисляет координаты и векторы скорости для данного момента времени.
Я переживал за возможности планшета в вычислительных операциях. Но, как стало видно, напрасно. :-)
Теперь могу просчитать трек на 500(!) минут вперёд и ничего не виснет.

[xd]

Видимо сборщик мусора с вычислительным потоком ссорился.

[INA]

Видимо сборщик мусора с вычислительным потоком ссорился.

Нет... :-)
Подвел JIT компилятор.
Теперь спокойно можно вывести на карту 10 спутников со 100 точками трека для каждого. )

[xd]

Подвел JIT компилятор.

Неожиданно...

[INA]

Да... Я тоже не ожидал.
Поэтому и заблудился в правильности своего алгоритма.

[INA]

Вот тут небольшое видео о работе софта.
За кадром голос Задорного... 

https://youtu.be/SAGbRQAmBtg

[Abdr]

Подвел JIT компилятор.

Неожиданно...

помогите написать прогу на матлабе орбиту навигационного спутника нарисовать надо через формулы движения по орбите плот нерисует круг или эллипс

[xd]

Для начала соизвольте не полениться, и расставьте знаки препинания.

[Abdr]

помогите пжл((( должен получится круг рисует конус. что не так?(((  
function GPS10()
t=405504;
toe=216000;
tk=t-toe;
e=0.00532;
A=26559.45305;
Mo=-176.53961*pi/180;
mu=3.986*10^14;
OMEGA=7.292*10^(-5);
OMEGADOT=-4.65589*10^(-7);
S=23.7702778;
omega=25.84150*pi/180;
io=55.21840*pi/180;
delt=0;
L=-89.76100;
%Расчетное среднее движение
no=sqrt(mu/A^3);
%Время, отсчитываемое от опорной эпохи эфемерид
%n=405504+3600;
%tk=zeros(n,1);
%for i=405504:1:n
%tk(i)=i;
%end
%toe-142214400*delt
%Средняя аномалия
Mk=Mo+no*tk;
%Эксцентрическая аномалия
Ek=Mk;
%Истинная аномалия
vk=atan(((sqrt(1-e^2).*sin(Ek))./(1-e*cos(Ek)))./((cos(Ek)-e)./(1-e*cos(Ek))));
%Аргумент широты
Fk=vk+omega;
OMEGAo=L+S-atan(tan(Fk)*cos(io))-OMEGADOT*t/43200;
%Скорректированный радиус орбиты спутника
rk=A.*(1-e*cos(Ek));
%Координаты навигационного спутника в орбитальной плоскости
xkk=rk.*cos(Fk);
ykk=rk.*sin(Fk);
%Скорректированная долгота восходящего узла
OMEGAk=OMEGAo+(OMEGA-OMEGADOT)*tk-OMEGADOT*toe;
%Координаты спутников в геоцентрической фиксированной системе координат (ECEF)
xk=xkk*cos(OMEGAk)-ykk*cos(io)*sin(OMEGAk)
yk=xkk*sin(OMEGAk)+ykk*cos(io)*cos(OMEGAk)
zk=ykk*sin(io)
hh=sqrt(xk.^2+yk.^2+zk.^2)-6371;

for T=1:1000:85000
    tk=t-toe;
    Mk=Mo+no*tk;
   
    Ek1=Mk;
eps = 10^(-15);
y = e * sin(Mk);
x1 = Mk;
x = y;
for k = 0 : 15 % количество итераций
x2 = x1;
x1 = x;
y1 = y;
y = Mk - (x - e * sin(x));
if (abs(y - y1) < eps)
break
end
x = (x2 * y - x * y1) / (y - y1);
end % kepler
Ek = x;
    vk=atan(((sqrt(1-e^2)*sin(Ek))/(1-e*cos(Ek)))/((cos(Ek)-e)/(1-e*cos(Ek))));
    Fk=vk+omega;
    rk=A*(1-e*cos(Ek));
    xkk=rk*cos(Fk);
    ykk=rk*sin(Fk);
    OMEGAk=OMEGAo+(OMEGA-OMEGADOT)*T-OMEGADOT*toe;
    xk(T)=xkk*cos(OMEGAk)-ykk*cos(io)*sin(OMEGAk);
    yk(T)=xkk*sin(OMEGAk)+ykk*cos(io)*cos(OMEGAk);
    zk(T)=ykk*sin(io);
    %p(T)=sqrt(xk(T)^2+yk(T)^2+zk(T)^2)-6371;
end
figure
%plot(T,p)
%plot(tk,p)
grid on;
%[X,Y,Z]=sphere(30); 
%surfl(X,Y,Z)
%[ X , Y , Z] = meshgrid( xk, yk, zk );
figure
plot3( xk , yk , zk )
grid on;
%hold off
end

[xd]

%Эксцентрическая аномалия
Ek=Mk;

Забыли уравнение Кеплера решить.

%Истинная аномалия
vk=atan(((sqrt(1-e^2).*sin(Ek))./(1-e*cos(Ek)))./((cos(Ek)-e)./(1-e*cos(Ek))));

Здесь неопределённость квадранта угла. Воспользуйтесь формулой
\( \displaystyle tg \frac{\nu}{2} = \sqrt{\frac{1+e}{1-e}} tg \frac{E}{2} \)

OMEGAo=L+S-atan(tan(Fk)*cos(io))-OMEGADOT*t/43200;

Тоже неопределённость квадранта.

    vk=atan(((sqrt(1-e^2)*sin(Ek))/(1-e*cos(Ek)))/((cos(Ek)-e)/(1-e*cos(Ek))));

И тут тоже.

Вообще тяжело читать код без нормальных отступов. Воспользуйтесь BB-тегом [ code ].
И приведите скриншоты графиков, которые Вам не нравятся, будет проще разбираться с происходящим.

[Михаил111]

Добрый день! А кто нибудь знает, где можно скачать альманахи Galileo?

[kodikv]

Добрый день! Подскажите пожалуйста, с чего начать, если необходимо определить координаты по псевдодальности, начально дан ринекс файл.

[Upsilon]

Добрый день! Подскажите пожалуйста, с чего начать, если необходимо определить координаты по псевдодальности, начально дан ринекс файл.

Координаты чего? Если фигурирует слово (псевдо)дальность, то участвуют координаты ИСЗ и координаты наземной точки.
По документине ftp://ftp.unibe.ch/aiub/rinex/rinex301.pdf формат Rinex для разных навигационных ИСЗ разный. О каком типе данных идет речь?

[kodikv]

Необходимо определить координаты точке на земле, дан ринекс файл с расширением nav. Сказали начать с решение системы четырех нелинейных уравнений. Нету ли у Вас алгоритма с наименьшими квадратами или где можно про него посмотреть?

[Upsilon]

http://www.miigaik.ru/library.miigaik.ru/uchebnieposobiya/20150324161133-4015.pdf

[R2_D2]

Скажите, пожалуйста, считал ли кто-нибудь местоположение спутника по Альманаху по алгоритму, предлагаемому в ИКД ГЛОНАСС? Можно ли по нему считать положение на витке, отличном от того, на который Альманах дан? Дело в том, что судя по алгоритму, поправки второй гармоники обнуляются при  каждом прохождении восходящего узла, соответственно траектория будет разрывной.