嘿，作为系统建模新手，从R转到Simulink确实会遇到时间逻辑上的小障碍，咱们一步步把你的太阳辐射模型搬过去，顺便搞懂Simulink的时间怎么玩。

先搞懂Simulink的时间逻辑 #

和R里手动生成时间序列（你代码里的hours向量）不一样，Simulink是基于仿真步长的自动时间推进环境：

• 你只需要在仿真配置里设定起始时间、结束时间、步长，Simulink会自动在每个步长生成当前时间，不需要手动创建时间数组。
• 注意：Simulink默认时间单位是秒，而你的R代码用的是小时，这是最容易踩坑的地方，一定要做好单位转换！

第一步：把R代码转成Simulink兼容的MATLAB逻辑 #

先把你的R计算逻辑转换成MATLAB函数（Simulink原生支持MATLAB代码），处理好边界条件：

function irradiance = solarRadiance(t_hour)
    % t_hour: 当前仿真时间（已转换为小时）
    % 计算时角（弧度）
    h_radians = (pi/12) * (t_hour - (12 + 26/60));
    % 年积日（可改成输入参数实现动态调整）
    doy = 268;
    % 计算赤纬角（弧度）
    decl_radians = 23.45 * sind(2*(284 + doy)/365) * pi/180;
    % 纬度（弧度）
    lat_radians = 23.45 * pi/180;
    % 计算sin(gamma)
    sin_gamma = sin(lat_radians)*sin(decl_radians) + cos(lat_radians)*cos(decl_radians)*cos(h_radians);
    
    % 处理夜晚/太阳在地平线以下的情况
    if sin_gamma <= 0
        irradiance = 0;
    else
        m = 1/sin_gamma;
        irradiance = 1353 * sin_gamma * (0.687)^(m^0.678);
    end
end

⚠️ 说明：

• MATLAB里用sind()处理角度输入的正弦，和R的sin()（弧度输入）逻辑对应，避免混淆。
• 主动处理sin_gamma <= 0的情况，防止出现NaN或无穷大的计算错误。

第二步：在Simulink里搭建模型 #

跟着这些步骤拖模块就行：

1. 拖一个MATLAB Function模块到画布，双击打开，把上面的代码粘贴进去，确认输入端口是t_hour，输出是irradiance。
2. 拖一个Clock模块：它会输出当前仿真时间（单位是秒），所以需要加一个Gain模块，把增益设为1/3600，把秒转换成小时，再连接到MATLAB Function的输入端口。
3. 拖一个Scope模块，连接到MATLAB Function的输出端口，用来查看辐照度的仿真结果。
4. 设置仿真参数：
   • 点击菜单栏的「Simulation」→「Model Configuration Parameters」
   • 把Stop Time设为(18+31/60 - 6-20/60)*3600（把你的时间范围转换成秒）
   • 选择Fixed-step solver，把Fixed-step size设为0.01*3600（对应你R里的0.01小时步长，转成秒）

第三步：验证结果 #

运行仿真后，Scope里的曲线应该和你在R里得到的结果一致。如果需要调整年积日、纬度等参数，只需要修改MATLAB函数里的变量，或者把这些参数做成Constant模块作为输入即可。

内容的提问来源于stack exchange，提问作者sametsokel