how to modify this code to estimate LCOE using the parameters for offshore wind?
CAPEX ESTIMATION (KUSD)
% CAPEX OF WIND TURBINE
USD = 1000; Pwt = 400; dp = 30; Nwt = 50; Lcc = 10620; Acc = 6691; Loffs = 10; Lonsh = 2;
CAPwt = 1497.9 * Pwt^0.87; T = 20; WACC = 0.10;
% CAPEX OF BUILDING
CAbui = (0.9181 * dp^2 – 31.43 * dp + 747.4) * 1.09 * Pwt;
% CAPEX OF PROJECT MANAGEMENT
CApm = 130.82 * Pwt;
% CAPEX OF SUPERVISORY CONTROL AND DATA ACQUISITION
CAscda = 555.1 * Nwt;
% CAPEX OF ACQUISITION OF WIND POWER CONNECTION CABLES
CAawpc = ((4.26 * 10^-4 * Acc) + (2.324 * 10^-1)) * 1.09 * Lcc;
% INSTALLATION OF WIND POWER CONNECTION CABLES
CAiwpc = ((-2.2684 * 10^-3 * Nwt) + (3.8018 * 10^-1)) * 1.09 * Lcc;
% ACQUISITION OF OFFSHORE TRANSMISSION CABLES
CAaotc = ((3.3565 * 10^-3 * Pwt) + (8.3872 * 10^-2)) * Loffs * 1.09;
% CAPEX OF INSTALLATION OF OFFSHORE TRANSMISSION CABLES
CAiotc = ((-6.318 * 10^-4 * Pwt) + (3.8125 * 10^-1)) * 1.09 * Loffs;
% CAPEX OF ACQUISITION OF ONSHORE TRANSMISSION CABLES
CAaontc = ((-3.3565 * 10^-3 * Pwt) + (8.3872 * 10^-2)) * 0.27 * Lonsh;
% CAPEX OF INSTALLATION OF ONSHORE TRANSMISSION CABLES
CAiontc = 0.578 * Lonsh;
% Total CAPEX
CApexT = (CAPwt + CAbui + CApm + CAscda + CAawpc + CAiwpc + CAaotc + CAiotc + CAaontc + CAiontc);
CAXNEW = CApexT – 0.2 * CApexT;
% OPEX
Dexex = 0.003 * CApexT;
% TOTAL CAPEX IN NAIRA
CApexTN = CApexT * USD * 10^3;
% TOTAL OPEX AND DEXEX IN NAIRA
DexexN = Dexex * USD * 10^3;
Opex = 0.035 * CApexT;
OpexN = Opex * USD * 10^3;
n = 20; Rb = 10; AEPv = 28298304; Civ = 668860000; Cdex = 2006600; Cpex = 23410000;
% Calculate LCOE
AEP = 28298304;
lcoe = calculate_lcoe(CAXNEW, Opex, WACC, AEP, T);
% Display the LCOE
disp([‘The LCOE is ‘, num2str(lcoe), ‘ currency units per kWh’]);
% Function definition
function lcoe = calculate_lcoe(CAXNEW, Opex, WACC, AEP, T)
% Calculate the numerator of the LCOE formula
numerator = CAXNEW;
for t = 1:T
numerator = numerator + Opex / (1 + WACC)^t;
end
% Calculate the denominator of the LCOE formula
denominator = 0;
for t = 1:T
denominator = denominator + AEP / (1 + WACC)^t;
end
% Calculate the LCOE
lcoe = numerator / denominator;
endCAPEX ESTIMATION (KUSD)
% CAPEX OF WIND TURBINE
USD = 1000; Pwt = 400; dp = 30; Nwt = 50; Lcc = 10620; Acc = 6691; Loffs = 10; Lonsh = 2;
CAPwt = 1497.9 * Pwt^0.87; T = 20; WACC = 0.10;
% CAPEX OF BUILDING
CAbui = (0.9181 * dp^2 – 31.43 * dp + 747.4) * 1.09 * Pwt;
% CAPEX OF PROJECT MANAGEMENT
CApm = 130.82 * Pwt;
% CAPEX OF SUPERVISORY CONTROL AND DATA ACQUISITION
CAscda = 555.1 * Nwt;
% CAPEX OF ACQUISITION OF WIND POWER CONNECTION CABLES
CAawpc = ((4.26 * 10^-4 * Acc) + (2.324 * 10^-1)) * 1.09 * Lcc;
% INSTALLATION OF WIND POWER CONNECTION CABLES
CAiwpc = ((-2.2684 * 10^-3 * Nwt) + (3.8018 * 10^-1)) * 1.09 * Lcc;
% ACQUISITION OF OFFSHORE TRANSMISSION CABLES
CAaotc = ((3.3565 * 10^-3 * Pwt) + (8.3872 * 10^-2)) * Loffs * 1.09;
% CAPEX OF INSTALLATION OF OFFSHORE TRANSMISSION CABLES
CAiotc = ((-6.318 * 10^-4 * Pwt) + (3.8125 * 10^-1)) * 1.09 * Loffs;
% CAPEX OF ACQUISITION OF ONSHORE TRANSMISSION CABLES
CAaontc = ((-3.3565 * 10^-3 * Pwt) + (8.3872 * 10^-2)) * 0.27 * Lonsh;
% CAPEX OF INSTALLATION OF ONSHORE TRANSMISSION CABLES
CAiontc = 0.578 * Lonsh;
% Total CAPEX
CApexT = (CAPwt + CAbui + CApm + CAscda + CAawpc + CAiwpc + CAaotc + CAiotc + CAaontc + CAiontc);
CAXNEW = CApexT – 0.2 * CApexT;
% OPEX
Dexex = 0.003 * CApexT;
% TOTAL CAPEX IN NAIRA
CApexTN = CApexT * USD * 10^3;
% TOTAL OPEX AND DEXEX IN NAIRA
DexexN = Dexex * USD * 10^3;
Opex = 0.035 * CApexT;
OpexN = Opex * USD * 10^3;
n = 20; Rb = 10; AEPv = 28298304; Civ = 668860000; Cdex = 2006600; Cpex = 23410000;
% Calculate LCOE
AEP = 28298304;
lcoe = calculate_lcoe(CAXNEW, Opex, WACC, AEP, T);
% Display the LCOE
disp([‘The LCOE is ‘, num2str(lcoe), ‘ currency units per kWh’]);
% Function definition
function lcoe = calculate_lcoe(CAXNEW, Opex, WACC, AEP, T)
% Calculate the numerator of the LCOE formula
numerator = CAXNEW;
for t = 1:T
numerator = numerator + Opex / (1 + WACC)^t;
end
% Calculate the denominator of the LCOE formula
denominator = 0;
for t = 1:T
denominator = denominator + AEP / (1 + WACC)^t;
end
% Calculate the LCOE
lcoe = numerator / denominator;
end CAPEX ESTIMATION (KUSD)
% CAPEX OF WIND TURBINE
USD = 1000; Pwt = 400; dp = 30; Nwt = 50; Lcc = 10620; Acc = 6691; Loffs = 10; Lonsh = 2;
CAPwt = 1497.9 * Pwt^0.87; T = 20; WACC = 0.10;
% CAPEX OF BUILDING
CAbui = (0.9181 * dp^2 – 31.43 * dp + 747.4) * 1.09 * Pwt;
% CAPEX OF PROJECT MANAGEMENT
CApm = 130.82 * Pwt;
% CAPEX OF SUPERVISORY CONTROL AND DATA ACQUISITION
CAscda = 555.1 * Nwt;
% CAPEX OF ACQUISITION OF WIND POWER CONNECTION CABLES
CAawpc = ((4.26 * 10^-4 * Acc) + (2.324 * 10^-1)) * 1.09 * Lcc;
% INSTALLATION OF WIND POWER CONNECTION CABLES
CAiwpc = ((-2.2684 * 10^-3 * Nwt) + (3.8018 * 10^-1)) * 1.09 * Lcc;
% ACQUISITION OF OFFSHORE TRANSMISSION CABLES
CAaotc = ((3.3565 * 10^-3 * Pwt) + (8.3872 * 10^-2)) * Loffs * 1.09;
% CAPEX OF INSTALLATION OF OFFSHORE TRANSMISSION CABLES
CAiotc = ((-6.318 * 10^-4 * Pwt) + (3.8125 * 10^-1)) * 1.09 * Loffs;
% CAPEX OF ACQUISITION OF ONSHORE TRANSMISSION CABLES
CAaontc = ((-3.3565 * 10^-3 * Pwt) + (8.3872 * 10^-2)) * 0.27 * Lonsh;
% CAPEX OF INSTALLATION OF ONSHORE TRANSMISSION CABLES
CAiontc = 0.578 * Lonsh;
% Total CAPEX
CApexT = (CAPwt + CAbui + CApm + CAscda + CAawpc + CAiwpc + CAaotc + CAiotc + CAaontc + CAiontc);
CAXNEW = CApexT – 0.2 * CApexT;
% OPEX
Dexex = 0.003 * CApexT;
% TOTAL CAPEX IN NAIRA
CApexTN = CApexT * USD * 10^3;
% TOTAL OPEX AND DEXEX IN NAIRA
DexexN = Dexex * USD * 10^3;
Opex = 0.035 * CApexT;
OpexN = Opex * USD * 10^3;
n = 20; Rb = 10; AEPv = 28298304; Civ = 668860000; Cdex = 2006600; Cpex = 23410000;
% Calculate LCOE
AEP = 28298304;
lcoe = calculate_lcoe(CAXNEW, Opex, WACC, AEP, T);
% Display the LCOE
disp([‘The LCOE is ‘, num2str(lcoe), ‘ currency units per kWh’]);
% Function definition
function lcoe = calculate_lcoe(CAXNEW, Opex, WACC, AEP, T)
% Calculate the numerator of the LCOE formula
numerator = CAXNEW;
for t = 1:T
numerator = numerator + Opex / (1 + WACC)^t;
end
% Calculate the denominator of the LCOE formula
denominator = 0;
for t = 1:T
denominator = denominator + AEP / (1 + WACC)^t;
end
% Calculate the LCOE
lcoe = numerator / denominator;
end my research MATLAB Answers — New Questions