Здравствуйте. Высылаю вам свое ДЗ.
В файле DZ_Sharov. m сделано:
1. Построение доступности спутников для выбранных станций.
2. Построены графики для C1, P1, P2, L1, L2, S1, S1, D1, D2, зависимость S1 от P1 и S2 от P2 для выбранных станций и различных спутников.
3. Построены линейные комбинации: Ionosphere Linear Combination (LI и PI), Widelane linear combination (PW и LW) и Melbourne-W Aubbena линейная комбинация.
В файле PathPlot. m должны строиться траектории спутников GPS, но получается не так, как должно быть, хотя что-то похожее на правду наблюдается. Скачал с сайта Mathworks GPS Navigation Toolbox GNT08.1.2 http://www. /matlabcentral/fileexchange/20578 (прикрепил к письму - GNT08.1.2.zip). Там тоже есть определение earth-fixed координат (в файле SV_Ephemeris_Model. m), разница только в том, что, видимо, исходные данные в файле project_data. mat имеют немного другой формат и угловая скорость везде домножается на пи. Определение Ек (Eccentric Anomaly) из уравнения Кеплера (заданного в Kepler_Eq. m) осуществляется с помощью функции fzero. Преобразование декартовых координат в сферические я делал также, как было и у вас, с помощью функции cart2sph. Полученные значения широты и долготы переводятся из радиан в градусы домножением на 180/pi. Расчет производится за сутки (86400 секунды), разделенные на 100 частей.
Но могу понять, что там не так.
P. S. GPS Navigation Toolbox для MatLab'а строит орбиты спутников в аксонометрии, Земля выполняется в виде эллипсоида, не понял как там подключить свой Navigation Message File.
С уважением,
Шаров Алексей
