1. 项目背景与核心价值风力发电作为清洁能源的重要组成部分其开发前期的资源评估至关重要。气象塔测量的历史风力数据是评估风电场选址和发电潜力的第一手资料但原始数据往往存在噪声、缺失和格式不统一等问题。这个项目正是针对这一痛点通过Matlab实现了一套完整的风能数据预处理和分析流程。我在新能源行业工作八年处理过上百个风电项目的评估数据。实际工作中发现很多工程师把60%的时间花在数据清洗和格式转换上真正用于分析的时间反而不足。这套代码的价值就在于自动化处理常见数据问题单位转换、时间对齐、异常值检测内置行业标准计算方法Weibull分布拟合、风玫瑰图生成可视化模块直接生成符合IEC标准的报告图表2. 数据准备与导入2.1 原始数据格式解析典型的气象塔数据包含以下关键字段时间戳 | 风速(m/s) | 风向(°) | 温度(℃) | 气压(hPa) | 湿度(%)实际接收的数据往往存在三种问题不同设备导出的时间格式不一致UTC时间/本地时间/自定义格式风速单位可能是节(knots)、公里/小时(km/h)等非标准单位数据记录可能存在跳变或设备故障导致的异常值2.2 Matlab数据导入实战推荐使用readtable函数处理CSV/TXT文件相比xlsread有更好的兼容性% 处理含中文表头的CSV文件 opts detectImportOptions(wind_data.csv); opts.VariableNames {Time,WS,WD,Temp,Press,Humidity}; rawData readtable(wind_data.csv, opts); % 时间格式转换处理Excel奇怪的日期格式 rawData.Time datetime(rawData.Time, ConvertFrom, excel);重要提示遇到混合格式数据时先用detectImportOptions预览数据结构避免直接导入导致类型识别错误。3. 数据预处理关键技术3.1 异常值检测与处理采用三级过滤机制物理范围过滤风速40m/s或0m/s视为无效变化率过滤相邻记录风速差5m/s需检查统计过滤3σ原则剔除离群值% 风速有效性检查 validIdx (rawData.WS 0) (rawData.WS 40) ... (abs(diff([0; rawData.WS])) 5); % 使用filloutliers处理剩余异常值 cleanData.WS filloutliers(rawData.WS(validIdx), linear);3.2 数据重采样与对齐不同传感器可能以不同频率记录数据需要统一到相同时间戳% 创建标准时间轴10分钟间隔 timeVector dateshift(rawData.Time(1), start, hour):... minutes(10):... dateshift(rawData.Time(end), end, hour); % 线性插值重采样 resampledWS interp1(rawData.Time, cleanData.WS, timeVector);4. 核心分析模块实现4.1 Weibull分布拟合风资源评估的核心是计算Weibull分布参数形状参数k和尺度参数c% 最大似然估计法计算参数 pd fitdist(cleanData.WS, Weibull); k pd.B; % 形状参数 c pd.A; % 尺度参数 % 可视化验证 histogram(cleanData.WS, Normalization, pdf); hold on x linspace(0, max(cleanData.WS), 100); plot(x, wblpdf(x, c, k), LineWidth, 2);4.2 风玫瑰图绘制风向频率分布可视化是选址的关键参考% 风向数据分箱16方位 windDirEdges linspace(0, 360, 17); [counts, ~] histcounts(cleanData.WD, windDirEdges); % 极坐标绘图 polarhistogram(BinEdges, windDirEdges*pi/180,... BinCounts, counts,... FaceColor, blue);5. 工程应用扩展5.1 发电量预估模型结合风机功率曲线进行发电量模拟% 典型风机功率曲线示例值 windSpeed 0:25; powerCurve [0 0 50 150 300 500 800 1200 1700 2000 2100... 2150 2200 2250 2300 2300 2300 2300 0 0 0 0 0 0 0]; % 计算理论发电量 energyOutput sum(interp1(windSpeed, powerCurve, cleanData.WS, linear, 0));5.2 数据报告自动生成使用mlreportgen包创建PDF报告import mlreportgen.dom.*; doc Document(WindReport, pdf); % 添加标题和摘要 title Paragraph(风能资源评估报告); title.Style {FontSize(18pt), Bold(true)}; append(doc, title); % 插入分析图表 img Image(which(weibull_fit.png)); img.Style {Width(100%)}; append(doc, img); close(doc);6. 实战经验与避坑指南时区处理陷阱气象数据通常用UTC时间而国内项目常用北京时间解决方法统一转换为本地时间后标注时区信息rawData.Time.TimeZone UTC; rawData.Time.TimeZone 08:00;缺失数据处理原则连续缺失2小时可用线性插值连续缺失2小时建议整段标记为无效总缺失率15%需重新评估数据质量性能优化技巧对于多年数据先用tall数组处理避免在循环中反复操作时间数据预分配数组内存特别是处理秒级数据时常见报错解决Invalid date/time format检查原始文件是否含隐藏字符Array indices must be positive通常是时间序列未正确排序导致Out of memory改用datastore增量读取大数据文件这套代码在实际风电项目中已验证过以下场景复杂地形下的多测风塔数据融合极端气候条件下的数据修正结冰、沙尘等不同型号测风设备的交叉校准对于想进一步深入的研究者建议扩展以下方向添加湍流强度分析模块集成测风激光雷达(LiDAR)数据对比开发基于机器学习的长期修正模型