风力发电机功率曲线统计MATLAB代码实现function windspeed_avr active_power_dev windspeed_fin_avr active_power_fin_avr CP= Powercurve_cal_v03(filename1,filename2,filename3,filename4)%读取 excel 的 xlsx 文件%B列为功率C列为风速data1=xlsread(filename1,1,B2:C52000);data2=xlsread(filename2,1,B2:C52000);data3=xlsread(filename3,1,B2:C52000);data4=xlsread(filename4,1,B2:C52000);%将四个 Excel 表合成一个表格data=data1;data2;data3;data4;%对数据进行筛选风速3m/s功率=3 & data(j,1)=10data_opt(j,1)=data(j,2);data_opt(j,2)=data(j,1);endend%对风速和功率进行排序并返回索引windspeed,ind=sort(data_opt(:,1);% 求0.5m/s的时间间隔的平均风速和相对应的平均功率,标准差m=length(ind);k=0;u=0;h=0;windspeed_sum=0;active_power_sum=0;windspeed_avr=0;active_power_avr=0;active_power_error_sum=0;windspeed_fin_sum=0;active_power_fin_sum=0;P=1.062;A=pi*(88/2)T;for windspeed_num=3:0.5:19.5for i=1:1:m%初步计算平均风速和平均功率 active_power(i)=data_opt(ind(i),2);if windspeed(i)=windspeed_num & windspeed(i)=windspeed_num & windspeed(i)=windspeed_num & windspeed(i)windspeed_num+0.5 e=abs(active_power(i)-active_power_avr(q)/active_power_dev(q);if e=4u=u+1;windspeed_fin_sum=windspeed_fin_sum+windspeed(i); active_power_fin_sum=active_power_fin_sum+active_power(i);endendwindspeed_fin_avr(q)=windspeed_fin_sum/u; active_power_fin_avr(q)=active_power_fin_sum/u;%计算发电机功率系数CP(q)=1000*active_power_fin_avr(q)/(0.5*P*A*(windspeed_fin_avr(q)A3); endwindspeed_sum=0;active_power_sum=0;active_power_error_sum=0;windspeed_fin_sum=0;active_power_fin_sum=0;k=0;h=0;u=0;end%画出初步平均风速和平均功率的曲线 figureplot(windspeed_avr, active_power_avr,r .-); set(gca,xtick,0:1:20)set(gca,ytick,0:100:1600)grid on xlabel(windspeed_avr m/s) ylabel(active_power_avr (kW) title(power curve)%画出标准差的曲线figure plot(windspeed_avr,active_power_dev,g .-); %bar(2*windspeed_num+0.5)/2,active_power_dev); set(gca,xtick,0:1:20)set(gca,ytick,0:20:200) grid onxlabel(windspeed_avr m/s) ylabel(active_power_dev) title(power standard deviation)%画出最终的平均风速和平均功率曲线 figureplot(windspeed_fin_avr, active_power_fin_avr,k .-); set(gca,xtick,0:1:20)set(gca,ytick,0:100:1600)grid on xlabel(windspeed_fin_avr m/s) ylabel(active_power_fin_avr (kW) title(final power curve)%画出最终的发电机功率系数 figureplot(windspeed_fin_avr,CP,m .-); set(gca,xtick,0:1:20)set(gca,ytick,0:0.1:0.6)grid on xlabel(windspeed_fin_avr m/s) ylabel(active_power_dev) title(CP 系数)%画出采集的散点图 figure plot(data(:,2),data(:,1),.); set(gca,xtick,0:1:22) set(gca,ytick,0:100:1600) grid on xlabel(windspeed m/s) ylabel(active_power) title(power point )