编写计算太阳位置的代码 % 输入参数 time 是时间点 % 输出太阳高度角 sun_altitude 和太阳方位角 sun_azimuth

以下是一个简单的计算太阳位置的 MATLAB 代码：

function [sun_altitude, sun_azimuth] = calculateSunPosition(time, longitude, latitude)
    % 时间转换为儒略日
    J2000 = datenum('2000-01-01 12:00:00');
    J = (time - J2000) / 36525 + 2451545.0;
    
    % 计算太阳平黄经
    L = mod(280.46646 + 36000.76983 * J, 360);
    
    % 计算太阳平近点角
    g = mod(357.52911 + 35999.05029 * J, 360);
    
    % 计算黄道倾角
    e = 0.016708634 - 0.000042037 * J - 0.0000001267 * J^2;
    
    % 计算地球公转轨道离心率
    obliquity = deg2rad(23.439291 - 0.0130042 * J - 0.00000016 * J^2 + 0.000000504 * J^3);
    
    % 计算太阳真黄经
    lambda = deg2rad(L + 1.914602 * sind(g) + 0.019993 * sind(2*g));
    
    % 计算太阳距离
    R = 1.00014 - 0.01671 * cosd(g) - 0.00014 * cosd(2*g);
    
    % 计算太阳黄道经度
    beta = asin(sin(lambda) * sin(obliquity));
    
    % 计算太阳高度角和方位角
    hour_angle = deg2rad((15 * (time - fix(time)) + longitude) - rad2deg(lambda));
    sin_altitude = sind(latitude) * sin(beta) + cosd(latitude) * cos(beta) * cos(hour_angle);
    sun_altitude = asin(sin_altitude);
    cos_azimuth = (sind(latitude) * cos(sun_altitude) - sin(beta)) / (cosd(latitude) * sin(sun_altitude));
    sin_azimuth = -cos(hour_angle) * sin(beta) / cos(sun_altitude);
    sun_azimuth = mod(rad2deg(atan2(sin_azimuth, cos_azimuth)) + 180, 360);
end

其中，`time` 参数是 MATLAB 的日期序列，`longitude` 和 `latitude` 参数是观测点的经度和纬度。函数的输出为太阳高度角 `sun_altitude` 和太阳方位角 `sun_azimuth`，单位均为弧度。