Scribd
Upload
36 views3 pages

CEP MATLAB Codes

This document discusses solving an economic load dispatch problem using a lambda iteration method. It provides MATLAB code to implement the lambda iteration method to calculate optimal generator outputs both without and with transmission line losses considered. The codes calculate the incremental costs, power outputs of generators, and check for convergence of the solution within a tolerance level over multiple iterations.

Uploaded by

WAQAR BHATTI
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
36 views3 pages

CEP MATLAB Codes

This document discusses solving an economic load dispatch problem using a lambda iteration method. It provides MATLAB code to implement the lambda iteration method to calculate optimal generator outputs both without and with transmission line losses considered. The codes calculate the incremental costs, power outputs of generators, and check for convergence of the solution within a tolerance level over multiple iterations.

Uploaded by

WAQAR BHATTI
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

Name: Hafiz Waqar Ahmad Bhatti

Reg. No: L1F20BSEE0053

Section: EA (Power)

Complex Engineeing problem (PSOC)

(b) Part 2: (WP2,WP3, WP4, EP2) - 10 marks • Design a lambda iteration program that can
solve the Lagrange function you have formulated in Q3 to find the incremental cost and the
power generated by each plant.

MATLAB Code (EDP Without Losses):

% Given data
Pd = 200; % Total load demand

% Cost characteristics [a, b, c] for each unit


A = [105 2.450 0.005;
44.1 3.510 0.005;
40.6 3.890 0.005];

% Maximum and minimum generation values for each unit


Pmin = [10; 20; 20];
Pmax = [160; 80; 50];

% Lambda Iteration Method


epsilon = 1e-5; % Convergence tolerance
maxIter = 1000; % Maximum number of iterations

% Initialize Lambda
lambda_min = min(A(:, 2)); % Minimum b value
lambda_max = max(A(:, 2)); % Maximum b value
lambda = (lambda_min + lambda_max) / 2;

iter = 0;
while iter < maxIter
% Calculate power generation for each unit
P = max(Pmin, min(Pmax, (lambda - A(:, 2)) ./ (2 * A(:, 3))));

% Calculate total generation and mismatch


Ptotal = sum(P);
deltaP = Pd - Ptotal;

% Check for convergence


if abs(deltaP) < epsilon
break;
end

% Update Lambda
if deltaP > 0
lambda_min = lambda;
else
lambda_max = lambda;
end
lambda = (lambda_min + lambda_max) / 2;

iter = iter + 1;
end

% Display results
if iter == maxIter
disp('Lambda Iteration did not converge.');
else
fprintf('Lambda Iteration converged in %d iterations.\n', iter);
fprintf('Lambda = %.6f\n', lambda);
fprintf('Power generation for each unit (MW):\n');
for i = 1:length(P)
fprintf('P%d = %.4f\n', i, P(i));
end
end

The Output:

(c) Part 3: (WP2, WP3, WP4, WP5) - 20 marks • Compute the Transmission Line Loss
Coefficients, reformulate the Economic Load Dispatch problem along with the Lagrange
function and use lambda iterative method to find the new set of generated powers.

MATLAB Code (EDP With Losses):

% Given cost characteristics (A matrix)


A = [105 2.450 0.005; % Cost coefficients for generator 1
44.1 3.510 0.005; % Cost coefficients for generator 2
40.6 3.890 0.005]; % Cost coefficients for generator 3

% Power demand
Pd = 200;

% Maximum and minimum generation values


Pmax = [160; 80; 50];
Pmin = [10; 20; 20];
% Calculated B-coefficients matrix for transmission losses
B = [0.0002 0.0001 0.00015;
0.0001 0.0003 0.0002;
0.00015 0.0002 0.00025];

% Lambda iteration parameters


epsilon = 1e-5;
max_iterations = 1000;

% Initial guess for lambda


lambda = 10;
d_lambda = 0.01;

% Function to calculate power output based on lambda


calculate_P = @(lambda) max(min((lambda - A(:,2)) ./ (2 * A(:,3)), Pmax), Pmin);

% Lambda iteration method


for iter = 1:max_iterations
P = calculate_P(lambda);
P_total = sum(P);
P_loss = P' * B * P;

% Check if total generation equals demand plus losses


if abs(P_total - Pd - P_loss) < epsilon
break;
end

% Adjust lambda
if P_total > Pd + P_loss
lambda = lambda - d_lambda;
else
lambda = lambda + d_lambda;
end
end

% Display results
fprintf('Lambda: %.5f\n', lambda);
fprintf('Iterations: %d\n', iter);
fprintf('Generator Outputs (MW):\n');
fprintf('P1: %.2f\n', P(1));
fprintf('P2: %.2f\n', P(2));
fprintf('P3: %.2f\n', P(3));
fprintf('Total Generation: %.2f MW\n', P_total);
fprintf('Transmission Losses: %.2f MW\n', P_loss);
fprintf('Total Demand + Losses: %.2f MW\n', Pd + P_loss);

The Output:

You might also like