import numpy as np
import [Link] as plt
import control as ctrl
# Thermal Reheat Unit Parameters
Kt = 1 # Turbine gain
Tt = 0.3 # Main turbine time constant
Tr = 10 # Reheat time constant
Kp = 1 # Proportional gain
Ki = 0.5 # Integral gain
Kd = 0.1 # Derivative gain
# Turbine Transfer Function
Gt_num = [Kt * (Tr + Tt), Kt] # Numerator coefficients
Gt_den = [Tt * Tr, (Tt + Tr), 1] # Denominator coefficients
Gt = [Link](Gt_num, Gt_den)
# PID Controller
Gc_num = [Kd, Kp, Ki]
# PID numerator (Derivative, Proportional, Integral)
Gc_den = [1, 0]
Gc = [Link](Gc_num, Gc_den)
# Closed-Loop System
Gol = [Link](Gc, Gt)
Gcl = [Link](Gol, 1)
# Step Response
time = [Link](0, 50, 500) # Time range
time, response = ctrl.step_response(Gcl, time)
# Plotting the Response
[Link](figsize=(8, 5))
[Link](time, response, label="Step Response")
[Link]('Python Step Response for LFC', fontsize=14)
[Link]('Time (s)', fontsize=12)
[Link]('Frequency Deviation (Hz)', fontsize=12)
[Link]()
[Link]()
[Link]()