GPS历书解算及天空视图绘制MATLAB

GPS历书解算及天空视图绘制MATLAB代码，经测试通过。 该程序可以读取YUMA格式的历书，然后根据设定的时间和地点，解算出GPS卫星的位置，绘制出当前位置处的GPS卫星分布图。

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