![Bit畉pm担nV畉nhnhvi畛ukhi畛nh畛th畛ngi畛n
SV: V探Ch鱈Li棚m MSSV 40901369
Bi 1. T鱈nhto叩ni畛u畛kinht畉choh畛th畛ngg畛m 3 m叩ynhi畛tc坦c叩cth担ngs畛nh動sau:
H1 = 225 + 7P1 + 0.0049P12 , 45MW P1 350 MW, FuelCost = 0.8 R/MBtu
H2 = 792 + 6.9P1 + 0.005P12 , 45MW P2 350 MW, FuelCost = 1.02 R/MBtu
H3 = 400 + 7.3P1 + 0.0059P12 , 47.5MW P3 450 MW, FuelCost = 0.9 R/MBtu
C担ngsu畉ty棚uc畉u l 433MW. S畛d畛ngph動董ngph叩pl畉p Lambda v畛i 了0 = 7.0 R/MWh.
Gi畉i:
Chi ph鱈ri棚ngc畛at畛ngm叩y :
F1 = H1.0,8 = 180 + 5.6P1 + 0.00392P12
F2 = H2.1,02 = 807.84 + 7.038P2 + 0.0051P22
F3 = H3.0,9 = 360 + 6.57P3 + 0.00531P32
畉ohm :
dF1/dP1 = 5.6 + 0.00784P1 = 了
dF2/dP2 = 7.038 + 0.0102P2 = 了
dF3/dP3 = 6.57 + 0.01062P3 = 了
P1 + P2 + P3 = Pyc = 433MW
Suyra : + + = 433
了 = 7.680728
P1 = 265.39898, P2 = 63.01255, P3 = 104.5883
Sai s畛 : 竜0 = Pyc - Pi = 433 ( 265.39898 + 63.01255 + 104.5883 ) = 0.00018 MW 0
V畉yk畉tqu畉 l 了 = 7.680728
Bi 2 :
Code :
gendata = [ 180 6.72 0.00248
743.58 6.834 0.007038
360 6.12 0.00369];
power = [100, 300, 100 ];
Pload = 500;
%form H
n = length ( gendata ) ;
H = zeros (n+1,n+1);
for i = 1 : n
H(i,i) = gendata (i,3) * 2;
H(i,n+1) = -1, H(n+1,i) = -1; end
x0 = zeros(n+1,1);
x0(1:n,1) = transpose( power );
%calculate the gradient and Hessian matrices
for kk = 1 : 10
disp(kk)
gradient = zeros(n+1,1);
gradient(n+1,1) = Pload;
for i =1 : n x0(n+1,1);
gradient(i,1)=gendata(i,2) + 2 * gendata(i,3) * x0(i,1) -
x0(n+1,1);](https://image.slidesharecdn.com/bt-130312015133-phpapp02/85/Bt-1-320.jpg)
![gradient(n+1,1)= gradient(n+1,1) - x0(i,1) ; end
dx = H gradient;
cost = 0;
for i = 1 : n
cost = cost +gendata(i,1) +gendata(i,2) * x0(i) +gendata(i,3) *
x0(i) *x0(i); end
disp( [x0', cost/1000])
x0 = x0 - dx;
end
K畉t qu畉 ta 動畛c 叩p 叩n B](https://image.slidesharecdn.com/bt-130312015133-phpapp02/85/Bt-2-320.jpg)
Bt
- 1. Bit畉pm担nV畉nhnhvi畛ukhi畛nh畛th畛ngi畛n SV: V探Ch鱈Li棚m MSSV 40901369 Bi 1. T鱈nhto叩ni畛u畛kinht畉choh畛th畛ngg畛m 3 m叩ynhi畛tc坦c叩cth担ngs畛nh動sau: H1 = 225 + 7P1 + 0.0049P12 , 45MW P1 350 MW, FuelCost = 0.8 R/MBtu H2 = 792 + 6.9P1 + 0.005P12 , 45MW P2 350 MW, FuelCost = 1.02 R/MBtu H3 = 400 + 7.3P1 + 0.0059P12 , 47.5MW P3 450 MW, FuelCost = 0.9 R/MBtu C担ngsu畉ty棚uc畉u l 433MW. S畛d畛ngph動董ngph叩pl畉p Lambda v畛i 了0 = 7.0 R/MWh. Gi畉i: Chi ph鱈ri棚ngc畛at畛ngm叩y : F1 = H1.0,8 = 180 + 5.6P1 + 0.00392P12 F2 = H2.1,02 = 807.84 + 7.038P2 + 0.0051P22 F3 = H3.0,9 = 360 + 6.57P3 + 0.00531P32 畉ohm : dF1/dP1 = 5.6 + 0.00784P1 = 了 dF2/dP2 = 7.038 + 0.0102P2 = 了 dF3/dP3 = 6.57 + 0.01062P3 = 了 P1 + P2 + P3 = Pyc = 433MW Suyra : + + = 433 了 = 7.680728 P1 = 265.39898, P2 = 63.01255, P3 = 104.5883 Sai s畛 : 竜0 = Pyc - Pi = 433 ( 265.39898 + 63.01255 + 104.5883 ) = 0.00018 MW 0 V畉yk畉tqu畉 l 了 = 7.680728 Bi 2 : Code : gendata = [ 180 6.72 0.00248 743.58 6.834 0.007038 360 6.12 0.00369]; power = [100, 300, 100 ]; Pload = 500; %form H n = length ( gendata ) ; H = zeros (n+1,n+1); for i = 1 : n H(i,i) = gendata (i,3) * 2; H(i,n+1) = -1, H(n+1,i) = -1; end x0 = zeros(n+1,1); x0(1:n,1) = transpose( power ); %calculate the gradient and Hessian matrices for kk = 1 : 10 disp(kk) gradient = zeros(n+1,1); gradient(n+1,1) = Pload; for i =1 : n x0(n+1,1); gradient(i,1)=gendata(i,2) + 2 * gendata(i,3) * x0(i,1) - x0(n+1,1);
- 2. gradient(n+1,1)= gradient(n+1,1) - x0(i,1) ; end dx = H gradient; cost = 0; for i = 1 : n cost = cost +gendata(i,1) +gendata(i,2) * x0(i) +gendata(i,3) * x0(i) *x0(i); end disp( [x0', cost/1000]) x0 = x0 - dx; end K畉t qu畉 ta 動畛c 叩p 叩n B