MATLAB
JORGELALDocumentos de Investigación31 de Mayo de 2021
13.199 Palabras (53 Páginas)86 Visitas
RESPALDO DIAGRAMAS
[pic 1]
[pic 2][pic 3][pic 4][pic 5][pic 6][pic 7][pic 8][pic 9][pic 10][pic 11][pic 12][pic 13][pic 14][pic 15][pic 16][pic 17][pic 18][pic 19][pic 20][pic 21][pic 22][pic 23]
RESPALDOS MATLAB – GRAFICAS
Fig 4
close all
hold on
x1=[2015 2020 2025 2030 2035 2040 2045 2050];
y1=[0 16 38 66 88 93 98 99];
x2=[2015 2020 2025 2030 2035 2040 2045 2050];
y2=[ 0 16 37 51 71 84 96 97];
x3=[2015 2020 2025 2030 2035 2040 2045 2050];
y3=[0 9 13 17 21 26 30 35];
x4=[2015 2020 2025 2030 2035 2040 2045 2050];
y4=[0 16 37 65 88 96 98 99];
x5=[2015 2020 2025 2030 2035 2040 2045 2050];
y5=[0 16 31 40 60 83 96 97];
grid on
xlabel('Year')
ylabel('VRE penetracion [%de demanda anual]')
plot(x1, y1, 'k*-', x2,y2, 'gx--', x3,y3, 'bo--', x4,y4, 'kd--', x5,y5, 'g*--')
lgd = legend('BC','NCC', 'CP', 'BCWC', 'NCCWC');
axis([2015 2050 0 100])
[pic 24]
close all
hold on
x1=[2015 2020 2025 2030 2035 2040 2045 2050];
y1=[0 2.5 2.5 9.2 11 12.5 11.4 11];
x2=[2015 2020 2025 2030 2035 2040 2045 2050];
y2=[ 0 2.5 2 5.5 8 11.8 12 12.4];
x3=[2015 2020 2025 2030 2035 2040 2045 2050];
y3=[0 4.8 3 2.5 2.2 2 1.8 1.7];
x4=[2015 2020 2025 2030 2035 2040 2045 2050];
y4=[0 0.5 0.5 0.3 0.2 0.3 0.2 0.2];
x5=[2015 2020 2025 2030 2035 2040 2045 2050];
y5=[0 0.6 0.6 0.4 0.2 0.3 0.2 0.2];
x6=[2015 2020 2025 2030 2035 2040 2045 2050];
y6=[0.1 3.4 3.1 12 15.2 18.2 20 20.5];
x7=[2015 2020 2025 2030 2035 2040 2045 2050];
y7=[0 3.4 3.2 7.5 11.5 13.4 18.7 19 ];
x8=[2015 2020 2025 2030 2035 2040 2045 2050];
y8=[0.1 6.4 4.6 4.9 4 3.8 3 2.8 ];
x9=[2015 2020 2025 2030 2035 2040 2045 2050];
y9=[0 1.8 2 5.8 8.8 11.8 12.4 12.1];
x10=[2015 2020 2025 2030 2035 2040 2045 2050];
y10=[0 1.8 2 2.7 5.4 7.3 9.3 10.4];
x11=[2015 2020 2025 2030 2035 2040 2045 2050];
y11=[0 0.8 0.8 2.1 4.2 5 8.9 9.3 ];
x12=[2015 2020 2025 2030 2035 2040 2045 2050];
y12=[0 2 2 2.2 2.5 2.4 6.5 6.4];
x13=[2015 2020 2025 2030 2035 2040 2045 2050];
y13=[0 2.5 2.2 2.6 2.7 2.7 2 1.8];
x14=[2015 2020 2025 2030 2035 2040 2045 2050];
y14=[0 1.9 2.3 6 9 11.8 12.2 12];
x15=[2015 2020 2025 2030 2035 2040 2045 2050];
y15=[0 1.8 2.4 2.5 5.2 7.5 9 10.2];
grid on
xlabel('Year')
ylabel('Pérdida [%de la generación total de RE]')
plot(x1, y1, 'k*-', x2,y2, 'Kx-', x3,y3, 'Ko--', x4,y4, 'kd--', x5,y5, 'Kp--', x6, y6, 'b*-', x7,y7, 'bx-', x8,y8, 'bo--', x9,y9, 'bd-', x10,y10, 'bp--',x11, y11, 'g*-', x12,y12, 'gx-', x13,y13, 'go--', x14,y14, 'gd-', x15,y15, 'gp--')
lgd = legend('Curtailment, BC','Curtailment, NCC', 'Curtailment, CP', 'Curtailment, BCWC', 'Curtailment, NCCWC','Total loss, BC','Total loss, NCC', 'Total loss, CP', 'Total loss, BCWC', 'Total loss, NCCWC','Store loss, BC','Store loss, NCC', 'Store loss, CP', 'Store loss, BCWC', 'Store loss, NCCWC');
axis([2015 2050 0 25])
[pic 25]
Fig 9
close all
hold on
x1=[110 800 2000 9000 11900];
y1=[65 88 93 97 99];
x2=[0 30 50 60 90 110 130];
y2=[9 11 15 21 25 30 35 ];
x3=[150 5000 8900 13000 15000];
y3=[64 88 96 98 99.5];
x4=[0 120 100 4500 5500];
y4=[1 70 85 94.5 95 ];
x5=[0 30 30 125 3900 8500 11000];
y5=[15 32 40 60 83 95 97];
x6=[0 100 2000 11900 ];
y6=[0 85 93 99];
x7=[0 150 3900 5000 8900 11000 13000 15000];
y7=[1 64 83 88 96 97 98 99.5];
grid on
xlabel('Capacidad de bateria [GWh]')
ylabel('VRE penetracion [% Demanda anual]')
plot(x1, y1, 'b+', x2,y2, 'kd', x3,y3, 'yx', x4,y4, 'mo', x5,y5, 'g*', x6,y6, 'c--',x7,y7, 'r-')
lgd = legend('BC','CP', 'BCWC', 'NCC', 'NCCWC','Curtailment allowed combined fit', 'No-curtailent combined fit');
axis([0 16000 0 100])
[pic 26]
close all
hold on
x1=[0 0.32 0.6 4.2 6 24.9 27];
y1=[1 38 65 89.5 93 98 99];
x2=[0.12 0.13 0.15 0.17 0.18 0.19 0.2];
y2=[9 11 15 21 25 30 35 ];
x3=[0 1 18 27 36 37];
y3=[16 64 88 96 97.5 98];
x4=[0 0.8 0.9 12 14.5];
y4=[15 71 85 96 97];
x5=[0 0.5 12 24 27];
y5=[15 60 84 95 96 ];
x6=[0 0.9 4.2 6 14.5 24.9 27];
y6=[0 85 89.5 93 97 98 99];
x7=[0 0.5 1 12 18 24 27 28 37];
y7=[1 60 64 84 88 95 96 95 98];
grid on
xlabel('Capacidad de bateria [demanda media diaria]')
ylabel('VRE penetracion [% Demanda anual]')
plot(x1, y1, 'b+', x2,y2, 'kd', x3,y3, 'yx', x4,y4, 'mo', x5,y5, 'g*', x6,y6, 'c--',x7,y7, 'r-')
lgd = legend('BC','CP', 'BCWC', 'NCC', 'NCCWC','Curtailment allowed combined fit', 'No-curtailent combined fit');
...