简介面向光伏系统设计与太阳能利用研究这份 MATLAB 程序包以双轴太阳跟踪和太阳辐射计算为核心适合需要模拟太阳位置、估算倾斜面辐射量或开发跟踪控制算法的工程师和科研人员。资源围绕辐射模型与跟踪策略展开涉及太阳赤纬角、时角计算并考虑地面反射、散射对接收辐射的影响可用于分析双面电池与双轴跟踪结合的发电增益。压缩包内含 17 个文件全部为 .m 源文件仅 33KB主程序与适应度计算、方向解算、参数更新等功能模块分开便于按需调用和二次开发。目前已有 125 人学习下载。通过这套代码可快速搭建辐射计算与跟踪角度求解的仿真框架理解双轴系统与双面电池协同工作的原理并方便修改地理、季节参数拓展实验节省从零编写算法的时间。1. 从 try_1.rar 说起双轴太阳跟踪要解决的三个问题双轴太阳跟踪器听起来像机械问题但拿到 try_1.rar 这类工程包后最先让人挠头的往往是数据太阳辐射记录仪的曲线和跟踪角度显示完全不在同一时间轴上两轴电机明明在转辐射量却时高时低。把压缩包名字拆开看其实就是三个要素双轴太阳跟踪负责“对准”太阳辐射负责“评价”后缀的 rar 说明这是一份要解包、要解析、要复现的原始资料。这里想讲清楚的内容是为什么双轴跟踪需要天文算法和辐射反馈配合如何从一份压缩包里的数据表格还原出跟踪误差以及现场调试时那些出厂文档不会写的坑。适合正在做光伏支架、聚光跟踪器或气象辐射观测的工程师。2. 双轴太阳跟踪的轴系选择与太阳位置算法2.1 高度角-方位角轴系的建模基础双轴太阳跟踪常见的机械结构有两种地平坐标系的高度角-方位角轴系和赤道坐标系的极轴-赤纬轴系。项目包里出现“方位角”“高度角”两个字段通常就是地平轴系。第一根轴垂直于地面负责把反射镜或光伏板转到太阳所在的方位角第二根轴平行于地面负责调整出射方向的高度角。优点是控制量直观安装时不需要对准北极轴缺点是方位角在正午前后变化太快尤其在低纬度地区一秒可能转十几角分普通步进电机容易丢步。建模时先定义本地太阳位置。给定观测点纬度 φ、太阳赤纬角 δ 和时角 ω高度角 α 的表达式为α arcsin(sin φ sin δ cos φ cos δ cos ω)方位角 γ 的符号和时角相关常见实现避免用 arctan 单独判断象限而是用 atan2。工程里我一般直接用坐标变换写成向量形式避免公式在过顶时出现奇点。这一步看似单纯后续跟踪误差、辐射修正全都基于这个坐标系。2.2 太阳位置算法与“跟踪太阳”的双重含义标题里的“跟踪太阳_辐射”值得拆一下。双轴太阳跟踪在控制上经常混用两种“跟踪”一是开环的天文跟踪根据时间计算太阳位置二是闭环的太阳跟踪用太阳辐射计或四象限传感器判断当前平面是否正对太阳。真正可靠的做法是前者做前馈、后者做修正。天文坐标告诉你正午 12 点 10 分太阳高度角约 65 度辐射传感器告诉你当前背光面还有 30 W/m² 的杂散这时要把修正量叠加到轴角指令上。Python 实现太阳位置计算时可以用 pvlib 的 solar_position 函数也可以用手动公式。我一般会在原型阶段用手动公式便于把每一步中间量写进日志import math import datetime def solar_position(lat, lon, ts, tz_hour): # 简化版天文算法只用于调试正式项目建议用 pvlib 或天文年历 day_of_year ts.timetuple().tm_yday # 赤纬角近似公式单位度 decl 23.45 * math.sin(math.radians(360 * (284 day_of_year) / 365)) # 当地时角以15度/小时换算午前为负 hour_angle 15 * (ts.hour ts.minute / 60 tz_hour - 12) lat_rad math.radians(lat) decl_rad math.radians(decl) h_rad math.radians(hour_angle) # 高度角 alt math.degrees(math.asin( math.sin(lat_rad) * math.sin(decl_rad) math.cos(lat_rad) * math.cos(decl_rad) * math.cos(h_rad))) # 方位角以正南为0向西为正 az math.degrees(math.atan2( -math.sin(h_rad), math.tan(decl_rad) * math.cos(lat_rad) - math.sin(lat_rad) * math.cos(h_rad))) return alt, az这段代码输入的 ts 是本地时间 datetimetz_hour 是时区相对 UTC 的小时偏移输出高度角单位度方位角以正南为 0向西为正。注意这里所有三角函数都先转成弧度最后输出再转回度否则正午附近的高度角误差会非常大。用这个函数跑一天再和辐射数据里的实际功率曲线对比就能看出问题在算法还是执行机构。2.3 双轴太阳跟踪项目的 6 个关键参数解包 try_1.rar 后第一步不是看代码而是找描述文件里这六个参数参数符号取值示例影响纬度φ39.9°N高度角和方位角整体偏移经度λ116.3°E地方时与标准时差时区偏移UTCi8太阳时计算起点磁偏角Δ-7°指南针校准方位零位海拔高度h50m大气折射修正量级辐射计采样周期Δt1s跟踪响应带宽上限磁偏角这个参数最容易漏。很多原型机用电子罗盘找零位罗盘给的是磁北极而太阳位置算法算的是真北方向高纬度地区磁偏角能到十几度不修正相当于跟踪系统一开机就带着十几度系统误差。海拔高度只在太阳高度角低于 5 度时有明显影响但如果项目包含聚光系统这个影响会被聚光倍数放大不能省略。3. 解包 try_1.rar 后辐射数据与跟踪控制闭环的复现路径3.1 先看数据太阳辐射记录与跟踪误差的对应关系拿到 try_1.rar常见做法是先解开按后缀名和体积分分类。通常里面会有几类内容控制器的配置文本、按日期命名的 CSV 或 DAT 辐射数据、可能还有几张轴角日志。把数据表头打出来head -5 radiation_hour.csv unzip -l try_1.rar先确认每列的时间戳是 UTC 还是本地时间辐射单位是 W/m² 还是 MJ/m²。辐射记录和轴角日志如果时间戳差很多后面所有分析都会失真。压缩包里常见的文件类型和处理动作可以靠下面这张表快速对齐文件类型常见后缀需要做的事配置文件.cfg/.ini提取纬度、经度、时区偏移辐射日志.csv/.dat统一时间戳、单位换算角度日志.log/.txt检查指令角度和实际角度说明文档.pdf/.md查找量程和采样周期3.2 控制闭环从太阳辐射传感器到电机驱动的信号链双轴太阳跟踪系统里辐射数据不只是“结果”还参与闭环。常见链路是太阳辐射计输出微伏级电压经过变送器转成 4-20mA 电流采样板按 1Hz 采样控制器根据当前太阳位置给出角度指令同时读取辐射传感器两翼的偏差信号生成校正量校正量送进步进驱动器。辐射数据在这个闭环里的作用类似伺服系统里的编码器反馈只不过被控量不是位置而是辐射通量。这个设计的关键参数有三个辐射计响应时间、变送器量程和采样同步时钟。一般热电堆辐射计的响应时间约 5s如果跟踪器两分钟才做一次位置修正反馈太慢如果压缩包里的辐射曲线看起来有相位滞后多半是采样没有和轴角指令打同一时间戳。3.3 用 Python 快速解析 try_1 的辐射数据并绘制日轨迹下面的代码演示怎么把辐射日志和计算出的太阳高度角放在一张图里判断跟踪是否到位。假设压缩包解出的文件是 solar_log.csv包含 timestamp、direct_irradiance 两列。import pandas as pd import matplotlib.pyplot as plt # 读取原始日志时间列不是UTC就转成带时区的本地时间 df pd.read_csv(solar_log.csv, parse_dates[timestamp]) df[timestamp] df[timestamp].dt.tz_localize(Asia/Shanghai) # 对每个时间点计算理论太阳高度角 df[alt_theory] df[timestamp].apply( lambda t: solar_position(39.9, 116.3, t.to_pydatetime(), 8)[0]) # 辐射曲线和理论高度角曲线共享时间轴 fig, ax1 plt.subplots(figsize(10, 5)) ax1.plot(df[timestamp], df[direct_irradiance], colororange) ax1.set_ylabel(Direct irradiance (W/m²)) ax2 ax1.twinx() ax2.plot(df[timestamp], df[alt_theory], b--, alpha0.6) ax2.set_ylabel(Solar altitude (deg)) plt.xticks(rotation45) plt.tight_layout()解析逻辑里最要紧的是时间戳规范化如果 CSV 里时间格式是“20230801 12:00:05”pandas 还能猜但如果是“2023-08-01T12:00:050800”就需要在 read_csv 时指定 strict 时间格式。双轴太阳跟踪的日志通常每秒一条一条错误的时间记录会让整段曲线平移几度比机械误差更难排查。4. 双轴太阳跟踪系统调试中的 5 个必查点与常见坑必查点典型故障关键参数零位早晚偏差大home offset辐射表读数偏低安装倾角、遮阴天气切换电机频繁启停辐射阈值丢步阶梯状曲线驱动器电流时间戳曲线平移采集延迟4.1 机械零位与限位开关的校准双轴太阳跟踪的零位不在地面而在“高度角 0 度、方位角对准正南”的位置。用水平尺校准完高度轴后要把方位轴开到软件认为的 0 度再用影子法核验插一根铅垂杆影子最短的时刻是太阳中天方位角应为 0。很多 try_1.rar 里的寄存器地址表会写“home offset”字段这个值就是用来把机械限位和软件零位对上的。限位开关不能只装末端还要在正南 0 度和天顶 90 度附近预留软限位防止雨天大风时电机越过机械限位。4.2 太阳辐射表的安装角度与遮阴影响辐射表必须放在跟踪平台中央感光面与光伏板或聚光镜面平行有些支架为了排水会让板面倾斜但辐射表的底座如果没垫平读数就会系统性偏低。调试时检查两个现象一是早晚低角度时辐射表周围有栏杆或线缆阴影二是辐射表线缆被风吹动造成电压波动。用两套辐射表交叉验证是最省力的办法把第二块表装成固定倾角比较它的日累计量和移动平台的累计量。4.3 天气突变时的跟踪策略切换阴天时太阳辐射被云层折射直接辐射分量很小跟踪器如果还按照晴天算法追高精度反而会让电机频繁启停。我一般会在控制流程里加一条切换规则当辐射传感器测到总辐射低于 120 W/m² 且变化率超过 10 W/m²/s就切换到“平均太阳位置”模式只做慢速追踪。try_1.rar 里如果包含天气数据可以按这个阈值回放一遍看看切换时电机指令有没有突变。4.4 电机步进丢步与传感器漂移双轴太阳跟踪的两个轴向受力不同方位轴要克服线缆扭力和风载高度轴要克服重力矩。步进电机驱动器如果电流设置偏小在大风天会丢步表现为高度曲线上出现阶梯状缺口。排查时可以对比编码器回读角度和指令角度如果没有编码器就在每天太阳中天时记录一次辐射峰值峰值慢慢变小就说明累计误差在增加。4.5 数据记录时间戳与太阳时的对齐这是最隐蔽的坑。辐射计调理电路通常有模拟滤波采样值相对真实辐射会延迟 1 到 2 秒控制器把角度写入日志时又可能带进一个延时。两边时间戳不齐画出来就是一条曲线比另一条晚几秒看似毫不起眼但用于标定跟踪误差时最多能引入 0.15 度的角度偏移。建议把所有日志统一转成 Unix 时间戳打印出每个数据点的采集和处理延迟。下面这段代码用来检查两个日志的计时偏差import pandas as pd rad pd.read_csv(radiation.csv, parse_dates[ts]) pos pd.read_csv(position.csv, parse_dates[ts]) # 匹配最近时间戳输出最大间隔 merged pd.merge_asof(rad.sort_values(ts), pos.sort_values(ts), onts, tolerancepd.Timedelta(1s)) gap (merged[ts_pos] - merged[ts]).dt.total_seconds() print(gap.describe())merge_asof会按照最近时间戳连接gap 如果出现几百毫秒以上的系统性偏差就要回去检查采集线程的调度优先级。5. 用辐射数据验证跟踪效果直接辐射与散射辐射分离技巧在缺少直接辐射表的情况下可以用遮挡法把总辐射表中的散射分量分离出来从而验证双轴太阳跟踪效果。标准做法是拿一根遮光环挡住太阳圆面测量散射辐射 D不遮挡时测总辐射 I。直接辐射 B I - D。如果跟踪器的平台相对太阳有偏差遮光环并不能完全挡住太阳分离出的直接辐射会比真实值减少。具体到验证流程我会在控制器的日志里加两个通道一个通道是普通总辐射表另一个通道固定在同一平台但装有遮光环采样间隔 1s。在中午连续记录 1 小时计算每个五分钟窗口的平均直接辐射和理论直接辐射把大气质量、太阳高度角代入然后观察两者的偏差。如果偏差稳定的单向变大说明跟踪漂移如果偏差是随机的则是云或气溶胶造成。下面是一个简单的散占比统计脚本用来量化跟踪误差对辐射分离的干扰import pandas as pd import numpy as np df pd.read_csv(two_radiation_log.csv, parse_dates[ts]) df[direct] df[total] - df[shaded] # 计算散射占比 df[scatter_ratio] df[shaded] / df[total] # 太阳高度角大于20度时统计 valid df[df[alt] 20] print(valid[scatter_ratio].describe()) # 如果散射占比稳定在30%以上高度角又很高需要怀疑遮光环未对准太阳 high_ratio valid[valid[scatter_ratio] 0.3] if len(high_ratio) 10: print(跟踪偏差疑似大于2度检查遮光环安装和太阳位置计算)这里 scatter_ratio 的数值可以作为现场告警阈值晴天散射占比通常在 10% 到 20% 之间如果超过 30% 且持续 5 分钟说明遮光环相对太阳已经偏移。利用这个技巧不需要高精度气象站就能把“跟踪准不准”变成一个数字指标也算把 try_1.rar 里辐射数据用到了实处。本文还有配套的精品资源点击获取