clc;
clear all;
load('FS_1_1.mat')
load('FS_1_2.mat')
load('FS_1_3.mat')
load('FS_1_4.mat')
load('FS_1_6.mat')
load('FS_1_8.mat')
load('FS_1_9.mat')
load('FS_1_10.mat')
load('FS_1_11.mat')
load('FS_1_12.mat')
load('FS_1_13.mat')
load('FS_1_14.mat')
load('FS_1_15.mat')
load('FS_1_16.mat')
load('FS_1_17.mat')
load('FS_1_18.mat')
load('FS_1_20.mat')
load('FS_1_21.mat')
load('FS_1_22.mat')
load('FS_1_23.mat')
load('FS_1_24.mat')
load('FS_1_25.mat')
load('FS_1_26.mat')
load('FS_1_27.mat')
load('FS_1_28.mat')
load('FS_1_29.mat')
load('FS_1_30.mat')
load('FS_1_31.mat')
load('FS_1_32.mat')
load('FS_1_33.mat')
FS_Z=FS_1_1(1:4464)+FS_1_2(1:4464)+FS_1_3(1:4464)+FS_1_4(1:4464)+FS_1_15(1:4464)+FS_1_6(1:4464)+FS_1_7(1:4464)+FS_1_8(1:4464)+FS_1_9(1:4464)+FS_1_10(1:4464);
FS_Z=FS_Z+FS_1_11(1:4464)+FS_1_12(1:4464)+FS_1_13(1:4464)+FS_1_14(1:4464)+FS_1_15(1:4464)+FS_1_16(1:4464)+FS_1_17(1:4464)+FS_1_18(1:4464)+FS_1_19(1:4464)+FS_1_20(1:4464);
FS_Z=FS_Z+FS_1_21(1:4464)+FS_1_22(1:4464)+FS_1_23(1:4464)+FS_1_24(1:4464)+FS_1_25(1:4464)+FS_1_26(1:4464)+FS_1_27(1:4464)+FS_1_28(1:4464)+FS_1_29(1:4464)+FS_1_30(1:4464);
FS_Z=FS_Z+FS_1_31(1:4464)+FS_1_32(1:4464)+FS_1_33(1:4464);
FS_Z_av=FS_Z/33;
figure(1)
plot(FS_Z_av(1+144*7:144*8),'b-');
set(gca,'xtick',0:12:144)
set(gca,'xticklabel',{'0','2','4','6','8','10','12','14','16','18','20','22','24'})
xlabel('\bf小时');
ylabel('\bf风电场平均风速曲线/(m/s)');
FS_Z_d_av=sum(FS_Z_av(1+144*7:144*8))/144;
FS_Z_av_LeiJi=FS_Z_av(1:144);
for i=2:31
FS_Z_av_LeiJi=FS_Z_av(144*i+1-144:144*i)+FS_Z_av_LeiJi;
end
FS_Z_av_LeiJi_d=FS_Z_av_LeiJi/31;
figure(2)
plot(FS_Z_av_LeiJi_d,'b-');
set(gca,'xtick',0:12:144)
set(gca,'xticklabel',{'0','2','4','6','8','10','12','14','16','18','20','22','24'})
xlabel('\bf小时');
ylabel('\bf风电场累计平均风速曲线/(m/s)');
FS_Z_av_LeiJi_d(1,1)=4.700;
FS_Z_dl_av=sum(FS_Z_av_LeiJi_d)/144;
%统计整个风电场的风速概率分布
figure(3)
h = histogram(FS_Z_av,'Normalization','probability')
xlabel('\bf风速/(m/s)');
ylabel('\bf风电场平均风速概率分布');
%单台风电机组风速概率统计
figure(4)
h = histogram(FS_1_1,'Normalization','probability')
xlabel('\bf风速/(m/s)');
ylabel('\bf风电机组风速概率分布');
%统计持续风速概率
FS_Z_av_1=FS_Z_av-4.5;
FSshiduan_1=zeros();
y_1=find(FS_Z_av_1>0);
n_1=length(y_1);
j_1=1;m_1=1;
FSshiduan_1(1,1)=FS_Z_av_1(y_1(1));
shichang_1=zeros();
for i=1:1:m_1
s_1=find(FSshiduan_1(i,:)>0);
shichang_1(i)=length(s_1);
end
%求解时段大于6的总时长
shichang_Z=0;
for i=1:1:m_1
if shichang_1(i)>6
shichang_Z=shichang_1(i)+shichang_Z;
end
end
figure(5)
t=3:10;
h=[70.2 59.8 47.3 30.7 17.9 8.4 4.5 2];
bar(t,h);
xlabel('\bf风速');
ylabel('\bf持续风速概率/%');
评论0