MERCADOS ELÉCTRICOS TAREA 1 JOHAN R. RIVERO V.
Enviado por Albert • 15 de Enero de 2019 • 529 Palabras (3 Páginas) • 424 Visitas
...
[pic 10]
[pic 11]
[pic 12]
[pic 13]
[pic 14]
[pic 15]
[pic 16]
[pic 17]
El código realizado en Matlab para resolver la minimización utilizando linprog es el siguiente:
%Código para calcular el mÍnimo de la funcion de costos
%con su restricciones correspondientes usando LINPROG
%Cp1 = 500+15*P1
%Cp2 = 600+25*P2
%Cp3 = 700+35*P3
% Minimizar mincost= Cp1+Cp2+Cp3
% sujeta a las restricciones siguientes:
% P1+P2+P3=310
% 10
% 20
% 30
clc; clear all; format compact;
CP=[15 25 35]+500+600+700; % coeficientes de la función objetivo
P0=[1 1 1]; %valores de arranque de la iteración
A=[];
b=[];
Aeq=[1 1 1]; % coeficientes de las ecuaciones de restricciones de igualdad
beq=[310]; % término independiente de restricciones de igualdad
Lb=[10 20 30]; % límite inferior para las variables P1, P2 y P3
Ub=[100 200 300]; % límite superior para las variables P1, P2 y P3
options=optimset('LargeScale','off','Display','iter'); %Para minimizar con las opciones de optimizacion de matlab
[P,CPmin,lambda] = linprog(CP,A,b,Aeq,beq,Lb,Ub,P0,options);
P1=P(1); P2=P(2); P3=P(3); % reasignacion de valores
% despliegue de resultados
fprintf(' -----------------------------------------------------------\n')
fprintf(' ----------- despliegue de resultados ----------------------\n')
fprintf(' -----------------------------------------------------------\n')
fprintf('P1 = %10.4f\n',P1)
fprintf('P2 = %10.4f\n',P2)
fprintf('P3 = %10.4f\n',P3)
fprintf('CPmin = %10.4f\n',CPmin)
LA SOLUCIÓN OBTENIDA DE LA MINIMIZACIÓN DE LOS COSTOS EMPLEANDO EL CÓDIGO ANTERIOR ES LA SIGUIENTE:
-------------------------------------------------------------
------------ Despliegue de resultados ----------------
-------------------------------------------------------------
P1 = 100.0000
P2 = 180.0000
P3 = 30.0000
CPmin = 565050.0000
La solución de los problemas anteriores manualmente, empleando el método de los multiplicadores de Lagrange, es el siguiente:
λ[pic 18][pic 19]
λ[pic 20][pic 21]
λ([pic 22][pic 23]
λ[pic 24]
λ[pic 25]
λ[pic 26]
[pic 27]
Resolviendo el sistema:
P1=179,21 MW
P2=264,04 MW
P3=-133,25MW
Ninguna de las potencias cumple con las restricciones de generación, por lo cual, basado en el hecho de que la ecuación de costo de producción de P1 parecer ser la más económica, se asumirá que P1 operará en su valor máximo de potencia permitido:
P1=100 MW
Por lo cual:
P3=30 MW
P2=180 MW
...