% STATCOM (Static Synchronous Compensator) Implementation in MATLAB
% Clear the workspace and command window
clear;
clc;
% System parameters
Vload_nominal = 1.0; % Nominal load voltage (pu)
Vload = 0.9; % Actual load voltage (pu)
Qload = 0.2; % Load reactive power (pu)
% STATCOM parameters
Vstatcom = 1.0; % STATCOM voltage magnitude (pu)
Qstatcom = 0.2; % STATCOM reactive power injection (pu)
% Calculate STATCOM current injection
Istatcom = (Qload + Qstatcom) / Vstatcom;
% Update load voltage with STATCOM compensation
Vload_compensated = Vload_nominal + 1j * Istatcom;
% Display the results
disp('Load Voltage without STATCOM:');
disp(Vload);
disp('Load Voltage with STATCOM:');
disp(abs(Vload_compensated));
% Calculate Voltage Stability Index (VSI) without STATCOM
VSI_without_STATCOM = abs(Vload_nominal - Vload) / Vload_nominal;
% Calculate Voltage Stability Index (VSI) with STATCOM
VSI_with_STATCOM = abs(Vload_nominal - abs(Vload_compensated)) / Vload_nominal;
disp('Voltage Stability Index (VSI) without STATCOM:');
disp(VSI_without_STATCOM);
disp('Voltage Stability Index (VSI) with STATCOM:');
disp(VSI_with_STATCOM);
% Plot the voltage phasors
figure;
hold on;
plot([0, real(Vload)], [0, imag(Vload)], 'b', 'LineWidth', 2);
plot([0, real(Vload_compensated)], [0, imag(Vload_compensated)], 'r', 'LineWidth', 2);
hold off;
axis equal;
xlim([-1.5, 1.5]);
ylim([-1.5, 1.5]);
xlabel('Real (pu)');
ylabel('Imaginary (pu)');
legend('Load Voltage without STATCOM', 'Load Voltage with STATCOM', 'Location', 'northeast');
title('Voltage Phasors with STATCOM Compensation');
