太阳方位角的matlab代码

以下是一段用 MATLAB 编写的计算太阳方位角的简单代码：

% 设置日期及位置信息（纬度、经度）
date = datetime(2019, 6, 21, 10, 0, 0); % 示例为2019年6月21日上午10点
lat = 39.9; % 北京市的纬度
lon = 116.4; % 北京市的经度

% 计算太阳位置相关参数
D = day(date, 'dayofyear');
decl = 23.45 * sind(360/365 * (284 + D));
t = 72.83 - 0.229 * B - 2.83 * sind(B) - 1.56 * sind(2*B);
B = 360/365 * (D - 81);
h = hour(date) + minute(date)/60 + second(date)/3600 - t/60;
omega = -120 / 180 * pi * h;

% 计算太阳方位角
cos_a = sind(decl) * sind(lat) + cosd(decl) * cosd(lat) * cos(omega - lon / 180 * pi);
a = acosd(cos_a);
if omega > lon / 180 * pi
    az = 360 - a;
else
    az = a;
end

以上代码仅供参考，实际使用中需要根据具体需求进行修改。