I need help with my Senior Project I am trying to implement a PID Controller that is able to detect when a building power demand is exceeding beyond its power rating it is using a Battery energy storage system to help supply the power towards the building. I am trying to Tune a PID controller using the PID tuner app to get discrete values that are above 5% of the given setpoint 93.75 (KVA). I am trying to get the kp, ki, and kd values and I don't have much experience using the PID tuner application. Attached is the excel file that I was working on to be able to get the values:
Code I was trying to use:
% === STEP 1: Load Power Data from Excel ===
data = readtable('battery_capacity_calculation_with_totals.xlsx');
% Display the column names to verify structure
disp("Column Names:")
disp(data.Properties.VariableNames)
% Extract the power column safely (fix column name if needed)
power_kw = data{:, "Power_kW"}; % Replace "Power_kW" if name differs
% Clean data: ensure numeric, finite, and valid
power_kw = fillmissing(power_kw, 'linear'); % Fill gaps
power_kw = power_kw(~isnan(power_kw)); % Remove remaining NaNs
power_kw = double(power_kw); % Ensure numeric
% Time vector in minutes (15-minute intervals)
time_minutes = (0:length(power_kw)-1) * 15;
% === STEP 2: Define Setpoint ===
setpoint = 93.75; % kW
dt = 15 * 60; % 15 minutes in seconds
% === STEP 3: Define First-Order Plant (or later ID from data) ===
G = tf(1, [300 1]); % First-order plant (adjust as needed)
% === STEP 4: Launch PID Tuner ===
pidTuner(G, 'PID'); % GUI will open
% === STEP 5: Export Tuned PID Controller (replace with your tuned values) ===
% After tuning, you’ll export something like:
% >> C = pid(Kp, Ki, Kd);
C = pid(0.6, 0.05, 0.01); % Example gains; replace after tuning
% === STEP 6: Closed-Loop Simulation ===
T = feedback(C * G, 1); % Closed-loop transfer function
% Simulate step response
figure;
step(T)
title('Step Response of Tuned PID Control System')
xlabel('Time (seconds)')
ylabel('Building Power (kW)')
grid on;
% === STEP 7: Optional – Apply PID to real power data ===
error = setpoint - power_kw;
integral = 0;
prev_error = 0;
pid_output = zeros(size(error));
simulated_power = zeros(size(error));
for i = 1:length(error)
e = error(i);
integral = integral + e * dt;
derivative = (i > 1) * (e - prev_error) / dt;
u = C.Kp * e + C.Ki * integral + C.Kd * derivative;
pid_output(i) = u;
simulated_power(i) = power_kw(i) + u;
prev_error = e;
end
% Plot actual vs. PID-controlled power
figure;
plot(time_minutes, power_kw, 'b-', 'LineWidth', 2); hold on;
plot(time_minutes, simulated_power, 'r--', 'LineWidth', 2);
yline(setpoint, 'k:', 'Setpoint');
xlabel('Time (minutes)');
ylabel('Power (kW)');
title('PID-Controlled Power vs Original Power');
legend('Original Power', 'PID-Controlled Power');
grid on;
