风电出力曲线生成代码

clc
clear
close all 
load('turbin5');
P_Wind1=sum(turbin5,2);    
%%  将功率等级折算成2348*0.65MW

P_Wind1=P_Wind1'/max(P_Wind1)*2348*0.65;  %MW   
figure(1)
plot(P_Wind1);
xlabel('时刻/min');       
ylabel('功率/MW');        
title('5月份某月风电1分钟分度值的出力');
%%  将5月份某月风电1分钟分度值的出力折算成15min的
P_Wind15 = zeros(1,30*24*4);
for i = 1:30*24*4
    P_Wind15(1,i) = sum(P_Wind1(1,15*i-14:15*i))/15;
end
figure(2)
plot(P_Wind15);
xlabel('时刻/min');       
ylabel('功率/MW');        
title('5月份某月风电15分钟分度值的出力');

%%   将15min与1min相比的波动值绘画出来
P_Wind15_1 = zeros(1,30*24*60);
for i = 1:30*24*4
    P_Wind15_1(1,15*i-14:15*i) = P_Wind15(1,i);
end
P_Wind1_bodong = P_Wind1(1,1:43200)-P_Wind15_1;
figure(3)
plot(P_Wind1_bodong);
xlabel('时刻/min')