風(fēng)力發(fā)電機(jī)功率曲線統(tǒng)計(jì)MATLAB代碼實(shí)現(xiàn)

風(fēng)力發(fā)電機(jī)功率曲線統(tǒng)計(jì)MATLAB代碼實(shí)現(xiàn)function[windspeed_avractive_power_devwindspeed_fin_avractive_power_fin_avrCP]=Powercurve_cal_v03(filename1,filename2,filename3,filename4)%讀取excel的xlsx文件%B列為功率C列為風(fēng)速data1=xlsread(filename1,1,'B2:C52000');data2=xlsread(filename2,1,'B2:C52000');data3=xlsread(filename3,1,'B2:C52000');data4=xlsread(filename4,1,'B2:C52000');%將四個(gè)Excel表合成一個(gè)表格data=[data1;data2;data3;data4];%對(duì)數(shù)據(jù)進(jìn)行篩選風(fēng)速<3m/s功率<10kw的數(shù)據(jù)不要n=length(data(:,2));forj=1:1:nifdata(j,2)>=3&&data(j,1)>=10data_opt(j,1)=data(j,2);data_opt(j,2)=data(j,1);endend%對(duì)風(fēng)速和功率進(jìn)行排序并返回索引[windspeed,ind]=sort(data_opt(:,1));%求0.5m/s的時(shí)間間隔的平均風(fēng)速和相對(duì)應(yīng)的平均功率,標(biāo)準(zhǔn)差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)^2;forwindspeed_num=3:0.5:19.5fori=1:1:m%初步計(jì)算平均風(fēng)速和平均功率active_power(i)=data_opt(ind(i),2);ifwindspeed(i)>=windspeed_num&&windspeed(i)<windspeed_num+0.5k=k+1;windspeed_sum=windspeed_sum+windspeed(i);active_power_sum=active_power_sum+active_power(i);endq=(windspeed_num-3)/0.5+1;windspeed_avr(q)=windspeed_sum/k;active_power_avr(q)=active_power_sum/k;%計(jì)算標(biāo)準(zhǔn)差ifwindspeed(i)>=windspeed_num&&windspeed(i)<windspeed_num+0.5h=h+1;active_power_error(i)=(active_power(i)-active_power_avr(q))^2;active_power_error_sum=active_power_error_sum+active_power_error(i);endactive_power_dev(q)=sqrt(active_power_error_sum/h);%剔除散點(diǎn)和平均功率差的絕對(duì)值與標(biāo)準(zhǔn)差比值大于4的點(diǎn)，并計(jì)算最終的平均風(fēng)速/功率曲線/發(fā)電機(jī)功率系數(shù)ifwindspeed(i)>=windspeed_num&&windspeed(i)<windspeed_num+0.5e=abs(active_power(i)-active_power_avr(q))/active_power_dev(q);ife<=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;%計(jì)算發(fā)電機(jī)功率系數(shù)CP(q)=1000*active_power_fin_avr(q)/(0.5*P*A*(windspeed_fin_avr(q)^3));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%畫(huà)出初步平均風(fēng)速和平均功率的曲線figureplot(windspeed_avr,active_power_avr,'r.-');set(gca,'xtick',0:1:20)set(gca,'ytick',0:100:1600)gridonxlabel('windspeed_avrm/s')ylabel('active_power_avr(kW)')title('powercurve')%畫(huà)出標(biāo)準(zhǔn)差的曲線figureplot(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)gridonxlabel('windspeed_avrm/s')ylabel('active_power_dev')title('powerstandarddeviation')%畫(huà)出最終的平均風(fēng)速和平均功率曲線figureplot(windspeed_fin_avr,active_power_fin_avr,'k.-');set(gca,'xtick',0:1:20)set(gca,'ytick',0:100:1600)gridonxlabel('windspeed_fin_avrm/s')ylabel('active_power_fin_avr(kW)')title('finalpowercurve')%畫(huà)出最終的發(fā)電機(jī)功率系數(shù)figureplot(windspeed_fin_avr,CP,'m.-');set(gca,'xtick',0:1:20)set(gca,'ytick',0:0.1:0.6)gridonxlabel('windspeed_fin_avrm/s')ylabel('active_power_dev')title('CP系數(shù)')%畫(huà)出采集的散點(diǎn)圖figurepl