資源簡介
GPS歷書解算及天空視圖繪制MATLAB代碼,經測試通過。
該程序可以讀取YUMA格式的歷書,然后根據(jù)設定的時間和地點,解算出GPS衛(wèi)星的位置,繪制出當前位置處的GPS衛(wèi)星分布圖。
代碼片段和文件信息
function?[ELAZ]=calAzEl(Sat?Pos?PosLat?PosLon)
R2D?=?180?/?pi;
for(i=[1:32])
dX?=?Sat(1i)?-?Pos(1);
dY?=?Sat(2i)?-?Pos(2);
dZ?=?Sat(3i)?-?Pos(3);
????
????
????M00?=?-sin(PosLon);
M01?=?cos(PosLon);
M02?=?0.0;
M10?=?-sin(PosLat)?*?cos(PosLon);
M11?=?-sin(PosLat)?*?sin(PosLon);
M12?=?cos(PosLat);
M20?=?cos(PosLat)?*?cos(PosLon);
M21?=?cos(PosLat)?*?sin(PosLon);
M22?=?sin(PosLat);
dE?=?M00?*?dX?+?M01?*?dY?+?M02?*?dZ;
dN?=?M10?*?dX?+?M11?*?dY?+?M12?*?dZ;
dU?=?M20?*?dX?+?M21?*?dY?+?M22?*?dZ;
r?=?sqrt(dE*dE?+?dN*dN);
if?(r?1E-20)
AZ(i)?=?0;
EL(i)?=?90;
else
AZ(i)?=?atan2(dE?dN)?*?R2D;
EL(i)?=?atan2(dU?r)?*?R2D;
if?(AZ(i)?0)
AZ(i)?=?AZ(i)?+?360;
????????end
end???
end
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????163??2016-07-27?20:21??歷書下載網址.txt
?????文件??????18410??2016-07-26?21:24??20160617.ALM
?????文件????????753??2016-07-27?17:06??calAzEl.m
?????文件???????5293??2016-07-28?11:26??calsatpos.m
?????文件????????343??2016-07-28?13:20??getSOW.m
?????文件???????1217??2016-07-27?16:07??llh2xyz.m
?????文件????????631??2016-07-28?11:56??main.m
?????文件???????2408??2016-07-27?17:46??ReadAlm.m
?????文件???????6699??2016-07-27?17:31??skyPlot.m
?????文件???????1565??2013-08-02?20:31??xyz2enu.m
?????文件???????1615??2013-08-02?20:31??xyz2llh.m
-----------?---------??----------?-----??----
????????????????39097????????????????????11
- 上一篇:圖論著色問題matlab程序
- 下一篇:matlab課后題及實驗答案
評論
共有 條評論