Weighted average algorithm adjusted a novel (1 + FOPI)-FOPI-TID controller structure for AGC with integration of non-linearities and cyber-attack

Moses Awal; Michael Robson Atim; Jimmy Nabende Wanzala; Johnes Obungoloch; Mohamed Barakat

PMC · DOI:10.1038/s41598-026-37004-0·February 2, 2026

Weighted average algorithm adjusted a novel (1 + FOPI)-FOPI-TID controller structure for AGC with integration of non-linearities and cyber-attack

Moses Awal, Michael Robson Atim, Jimmy Nabende Wanzala, Johnes Obungoloch, Mohamed Barakat

PDF

Open Access

TL;DR

This paper introduces a new controller for power systems that improves frequency stability and handles cyber-attacks more effectively than existing methods.

Contribution

A novel (1 + FOPI)-FOPI-TID controller optimized with a weighted average algorithm for robust automatic generation control in smart grids.

Findings

01

The proposed controller reduces ITAE by 45.3% and improves settling times by 47.7% and 32.8%.

02

It shows robustness under ±25% parameter variations and random load changes.

03

The controller maintains the lowest RoCoF during cyber-attacks, demonstrating superior resilience.

Abstract

The integration of diverse energy sources and the advent of smart grids have intensified the challenges in load frequency management (LFM). Modern power systems are increasingly vulnerable to inherent nonlinearities, such as generation rate constraints, governor dead bands, boiler dynamics, and communication delays, as well as sophisticated cyber-attacks, which collectively threaten frequency stability and tie-line power balance. To address these challenges, this study proposes a novel cascade controller, designated as (1 + FOPI)-FOPI-TID, for robust automatic generation control in hybrid two-area power systems. The controller uniquely combines fractional-order (FO) dynamics with a tilt-integral-derivative stage and is optimized using a green metaheuristic, the weighted average algorithm (WAA). The WAA effectively balances exploration and exploitation to achieve superior parameter…

Linked entities

Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.

Species1

Homo sapiens(human · species)

Chemicals3

BD water PI

Diseases8

ITAE AGC CTD WAA PID SS BD LFC

Figures18

Click any figure to enlarge with its caption.

Transient reactions for case 2 with 4% SLC at area 1: (a) Δ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${F}_{1}$$\end{document}$ , (b) Δ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${F}_{2}$$\end{document}$ , and (c) ΔP_tie_.

Random step loads $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta P}_{L1}$$\end{document}$ in area 1.Fig. 12Transient response with various algorithm techniques.

Nuclear power plant transfer function model.

The literature controller structure: (a) PIDF-(1 + FOD), (b) PIFOD-(1 + PI), (c) PDF-(1 + PI), (d) PD-PI, and (e) (1 + PI)-PI-PID.

Proposed structure of (1 + FOPI)-FOPI-TID CC.

Diagram of search agent position updating based on exploitation strategies.

Cyber-attack signal model applied to the PS.

Controller performances under the cyber-attack: ΔF_1_, ΔF_2_, and ΔP_tie._Table 6RoCoF (Hz/s) for different controllers.ControllerMax. Amplitude (Area 1)Max. Amplitude (Area 2)WAA: (1 + FOPI)-FOPI-TID0.06541.3276WAA: PD-PI0.09671.3835WAA: (1 + PI)-PI-PID0.11581.7025WAA: PIFOD-(1 + PI)0.09301.5391WAA: PDF-(1 + PI)0.08801.7652PIDF(1 + FOD)0.09591.5442Significant values are in [bold].

RoCoF results of PS areas for cyber-attack for Area-1 and Area-2.

Transfer function model of the tested PSI.

The convergence characteristic curves of PSO, GWO, and the proposed FHO.Table 2ITAE values for WAA, MA, and GWO at case1.AlgorithmMinimumMaximumAverageSt. deviationWAA0.0138100.0154430.0142870.00085502MA0.0161090.0228110.0188630.00336873GWO0.0169670.0230630.0223990.00334436

Transient response with various controllers but same optimizer at case 2.

Keywords

Automatic generation control(1 + FOPI)-FOPI-TID controllerWeighted average algorithmCyber-attackPower systems nonlinearitiesRandom step loadEnergy science and technologyEngineeringMathematics and computing

Peer Reviews

No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.

Videos

No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.

Taxonomy

TopicsFrequency Control in Power Systems · Power System Optimization and Stability · HVDC Systems and Fault Protection

Full text

Introduction

The stability and efficiency of modern electrical power systems are fundamentally governed by the precise balance between generation and demand, a critical task managed by load frequency control (LFC) or automatic generation control (AGC). As global energy consumption rises and power grids evolve into increasingly interconnected networks, incorporating a diverse mix of conventional thermal plants, hydroelectric facilities, gas turbines, and renewable sources, the challenge of maintaining system frequency within narrow operational limits has become markedly more complex^1,2^. These interconnected systems, while enhancing reliability and enabling economic power exchange, exhibit heightened dynamic sensitivity. Even minor load perturbations can trigger frequency deviations and unscheduled tie-line power flows, potentially leading to instability, equipment stress, and, in extreme cases, cascading failures^3^. Consequently, the development of advanced, robust control strategies for AGC is paramount for ensuring the security and reliability of future power networks.

This imperative is further intensified by two distinct yet converging sets of challenges: inherent physical nonlinearities and emerging cyber-physical vulnerabilities. Traditional linear models often fail to capture the true dynamic behavior of practical systems. Critical physical constraints, including generation rate constraints (GRC) that limit how fast a plant can ramp output, the nonlinear governor dead band (GDB) that causes insensitivity to small speed changes, slow-acting boiler dynamics (BD) in thermal units, and inevitable communication time delays (CTD) in signal transmission, collectively degrade control performance^4–6^. These nonlinearities manifest as increased overshoot, prolonged settling times, and persistent oscillations following disturbances. Simultaneously, the digital transformation of the grid through smart technologies introduces a new frontier of risk. The very communication channels and digital controllers that enable advanced AGC also expose the system to sophisticated cyber-attacks. Threats such as stealthy resonance-based attacks (ResA), which subtly manipulate load signals to drive the system into instability at its natural oscillatory modes, pose a serious and plausible danger^7–9^. This dual challenge, managing complex physical dynamics while fortifying the system against intelligent cyber threats, defines a critical research gap in contemporary AGC design.

In response to the need for improved regulation, control strategies have progressed significantly beyond the ubiquitous Proportional-Integral-Derivative (PID) controller. While PID variants, often tuned via classical or heuristic methods, remain widely used for their simplicity^10^, their performance is frequently inadequate in handling the multifaceted nonlinearities and high variability of modern, multi-source grids^11,12^. This has spurred interest in cascade control (CC) architectures, where two or more controllers are arranged in series, offering superior disturbance rejection and set-point tracking by correcting errors at multiple stages^13^. This evolution highlights the need for more structured and multi-stage control architectures in modern LFM design. The search for optimal parameters for these advanced structures has naturally aligned with the rise of metaheuristic optimization. Algorithms such as TLBO^2^, PSO^14^, Whale optimization algorithm^15^, Jaya optimization^16^, differential evolution^17^, hybrid PSO-BFOA algorithm^18,19^, grey wolf optimization^20^, and others, have been extensively applied to tune AGC controllers, demonstrating the ability to find high-performance solutions in complex, non-convex search spaces.

A parallel and transformative advancement has been the incorporation of Fractional-Order (FO) calculus into control theory. Unlike integer-order controllers, FO controllers possess non-integer differentiation and integration orders, providing additional tuning parameters that grant greater flexibility to shape the dynamic response. This has led to the development of effective strategies like the FOPID controller^21,22^, and more sophisticated cascade designs such as the PIFOD–(1 + PI) and FOPD–FOPI^23,24^ for multi-area systems. The quest for further performance gains continues to yield novel hybrid structures, including tri-stage controllers like (1 + PI)–PI–PID^25^ optimized via enhanced artificial bee colony–particle swarm optimization (IABC–PSO)^25^, hybrid 2DOF–PID–TD controllers^26^, fuzzy-FOPI + PIDN schemes for deregulated two-area systems^27^, and fuzzy-PIDF–(1 + PI) controllers^28^. These developments collectively illustrate the growing shift toward multi-stage and FO-based control architectures for advanced AGC applications.

Despite these considerable advances, a synthesis of the literature reveals several persistent limitations that this work seeks to address. First, many studies validate controller performance under a simplified set of assumptions, often incorporating only one or two nonlinearities. There is a lack of comprehensive evaluation under the simultaneous action of GRC, GDB, BD, and CTD within a hybrid (thermal-hydro-gas-nuclear) generation environment. Second, while cybersecurity is recognized as a paramount concern, very few AGC studies integrate realistic, dynamic cyber-attack models into their simulation-based validation, leaving a gap in understanding controller resilience. Third, the proliferation of complex cascade controller structures is sometimes driven primarily by numerical optimization outcomes, with insufficient explanation of the control-theoretic rationale behind the specific signal flow and component selection. Finally, the claimed superiority of newer metaheuristic optimizers often lacks rigorous benchmarking against established algorithms in the specific context of the AGC tuning problem.

To bridge these identified gaps, this paper proposes a novel, carefully justified fractional-order cascade controller, designated as the (1 + FOPI)-FOPI-TID, specifically tailored for automatic generation control of two-area power systems. The proposed controller simultaneously utilizes the area control error (ACE) and frequency deviation ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta F$$\end{document}$ ) signals, incorporating fractional-order dynamics and a tilting integral block $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$s^{{ - 1_{{/_{n} }} }}$$\end{document}$ to enhance transient response, disturbance rejection, and robustness^29^. Optimal parameter tuning is achieved using the weighted average algorithm (WAA)^30^. The straightforward yet efficient search process of the Weighted Average Algorithm (WAA), which is based on deterministic weighted average of potential solutions based on their fitness scores, sets it apart from other metaheuristic techniques. WAA avoids complicated population dynamics and requires fewer control parameters than evolutionary and swarm-based algorithms, which rely on several stochastic operators and algorithm-specific parameters. WAA is ideally suited for controller tuning in large-scale interconnected power systems because the adaptive weighting technique allows for a natural balance between exploration and exploitation, leading to quick convergence, numerical robustness, and less computational effort.

In accordance with the “no-free-lunch” principle^31^, the dual exploration–exploitation capability of WAA ensures robust regulation in the presence of GRC, GDB, BD, CTD, and stealthy cyber threats. Its novelty lies in the dynamic adjustment of search space via weighted population positions, which enhances convergence stability.

The Contributions of this study are summarized as follows:

i.A novel tri-stage FO CC, designed as (1 + FOPI)-FOPI-TID, is proposed for AGC of PSI to enhance frequency regulation and tie-line power stability.
ii.The proposed controller employs dual input signals, $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ACE$$\end{document}$ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta F$$\end{document}$ , enabling faster disturbance rejection and improved transient performance compared to conventional single-input CC structures.
iii.The WAA, a recent metaheuristic algorithm, is utilized for the optimal tuning of the proposed controller parameters, providing an effective balance between exploration and exploitation and ensuring robust convergence characteristics.
iv.A realistic multi-source two-area power system model is developed by incorporating GRC, GDB, BD, and CTD nonlinearities.
v.The cyber-resilience of the proposed LFC scheme is explicitly investigated by subjecting the system to resonance-based cyber-attack scenarios, assessing frequency deviation and rate-of-change-of-frequency (RoCoF) performance.

Materials and methods

Modelling of hybrid power system

The investigated system in Fig. 1. is a two-area interconnected power system with multi-source generation in each control area. Area-1 consists of thermal reheat, hydro, gas, and nuclear power generation units, while Area-2 comprises an identical combination of generation sources.Fig. 1. Transfer function model of the tested PSI.

The power demand is represented in a step load form, with an increment of $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{l1}$$\end{document}$ for area1 and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{l2}$$\end{document}$ for area-2. The system includes two frequency variations $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta F}_{i}$$\end{document}$ , $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta F}_{1}$$\end{document}$ for area-1 and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta F}_{2}$$\end{document}$ for area-2, and tie-line power change $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{tie}$$\end{document}$ . Frequency regulation is achieved through feedback mechanisms utilizing the speed regulation (R) and frequency bias factors (B), which stabilize system dynamics by adjusting tie-line power and frequency deviations^32^.

The effect of combining tie-line power error and frequency fluctuation is attributed to the ACE^17,33^. The area-specific ACEs are:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} ACE_{1} = & \Delta P_{tie} + B_{1} \Delta F_{1} \\ ACE_{2} = & \Delta P_{tie} + B_{2} \Delta F_{2} \\ \end{aligned}$$\end{document}

In the event of a load disturbance on the PS, ACEs are employed as a regulating signal to reduce $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{tie}$$\end{document}$ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta F}_{\mathrm{i}}$$\end{document}$ to zero at the steady state.

The integration of the HVDC link in the PSI allows for additional investigation. The generation rates from thermal, hydro, gas, and nuclear units are managed by independent secondary controllers in each sector, with AC-DC parallel lines connecting these units, as illustrated in Fig. 1. The DC link is modeled using a first-order transfer function, characterized by a gain of $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{K}}_{\mathrm{DC}}$$\end{document}$ and a time constant of $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{T}}_{\mathrm{DC}}$$\end{document}$ . This first-order transfer function effectively describes the behavior of the DC link.

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{\mathrm{tie,DC}}=\frac{{\mathrm{K}}_{\mathrm{DC}}}{1+{\mathrm{T}}_{\mathrm{DC}}\mathrm{s}}$$\end{document}

The HVDC links are increasingly being utilized for power exchange across control areas alongside existing high-voltage AC links. Research indicates that HVDC links, when operating in tandem with conventional AC systems, provide cost-effective and environmentally friendly benefits while also ensuring higher-quality electric power delivery^34^. Moreover, it reinforces the system’s resistance to minor disruptions^35^.

Nonlinear test models

Several physical nonlinearities are incorporated into the test model of the PSI to enhance realism. As shown in Figs. 1 and 2, the test model includes nonlinearities such as GDB, GRC, BD, and CTD. To limit the rate of generation change, a significant nonlinear component known as GRC is added. The typical benchmark for thermal plants is a GRC value of 3% per minute. To prevent excessive generation, the GRC for the thermal system is specified as $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left|\Delta {P}_{g}\right|$$\end{document}$ ≤ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.0005 p.u.MW/s$$\end{document}$ . To efficiently control generation levels, the saturation block of the thermal unit is restricted to ± 0.0005. A GRC of 360% per minute for decreasing generation and 270% per minute for increasing generation was taken into consideration for hydropower plants. This means that for decreasing generation, |ΔP_g_(s)|≤ 0.06 p.u.MW/s, and for increasing generation, |ΔP_g_(s)|≤ 0.045 p.u.MW/s. According to^36^ and^37^, the system model includes a saturation block in the hydro turbine that is limited between -0.06 and 0.045. According to Gheisarnejad M.^38,39^, the GDB nonlinearity is spelled out as:Fig. 2. Boiler dynamic nonlinearity.

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${G}_{\mathrm{d}}\left(s\right)=\frac{0.8-\frac{0.2}{\pi }s}{1+{T}_{\mathrm{g}}s}$$\end{document}

Figure 3. presents a simplified dynamic model of a nuclear power plant for load frequency control analysis. It includes the governor, high- and low-pressure turbine stages with reheater dynamics, and a nuclear participation factor. Frequency deviation is processed through droop and governor dynamics, while turbine-reheater blocks capture the slow response of nuclear units. The aggregated turbine output determines the incremental nuclear power contribution, enabling assessment of nuclear participation in frequency regulation under load disturbances.Fig. 3. Nuclear power plant transfer function model.

Similarly, a boiler dynamic nonlinearity is added to the thermal reheat system, as illustrated in Fig. 2. ^11,40^. The figure presents a dynamic boiler-steam turbine model used in load frequency control studies of thermal power plants, capturing the interactions among the boiler, fuel system, pressure control loop, and steam turbine. The control input signal regulates steam flow to the turbine through appropriate gain factors, while feedback from drum and throttle pressure deviations ensures stable operation. Boiler dynamics are represented using a storage model, and the pressure control loop incorporates proportional-integral with lead-lag characteristics. The fuel system is modeled as a first-order process with time delay to account for fuel transport and combustion effects. Overall, the model effectively represents the dominant delayed and nonlinear dynamics of thermal units for frequency regulation analysis under load disturbances.

Additionally, a CTD is incorporated after the controller design to enhance the realism of the system. According to^41^ and ^42^, this study used a time delay parameter ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\uptau }_{d}$$\end{document}$ ) value of 0.25 s, which is normally in the [0:1] range:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${e}^{-s{\uptau }_{d}} = \frac{1-\frac{{\uptau }_{d }}{2}\mathrm{s}}{1+\frac{{\uptau }_{d }}{2}\mathrm{s}}$$\end{document}

Cyber attack

Cyber-attacks in LFM are increasingly important to consider in order to ensure the reliable and secure operation of modern PSs, as such attacks may disrupt control actions and lead to severe operational consequences. In this work, a resonance-based cyber-attack (ResA) is adopted as a physics-informed and literature-supported representation of stealthy load manipulation scenarios reported in recent studies on cyber–physical power systems. In a ResA, an attacker introduces load variations that follow a resonance source while remaining within permissible operating limits, thereby avoiding immediate detection by conventional protection mechanisms^9,43^. Although these variations are small in magnitude, their persistent nature can induce pronounced frequency deviations and elevated RoCoF, which may threaten system stability.

The attack parameters, including amplitude and frequency, are selected based on established literature to reflect realistic adversarial behavior rather than idealized worst-case conditions. This choice ensures that the injected disturbances remain credible and allows a fair comparative assessment, as all controllers are evaluated under identical attack conditions. In practical power systems, RoCoF-based protection relays are widely deployed to mitigate such disturbances. Typical RoCoF relay thresholds lie within the range of 0.1–1.0 Hz/s, depending on system inertia characteristics^9^, and exceeding these limits may trigger protective actions that result in major system disruptions. Under the considered ResA, the manipulated load signal is verified for plausibility by the power generation system; once accepted, the attacker can influence system dynamics in a sustained manner^9,43^.

In this study, RoCoF is employed as the primary cyber-resilience indicator due to its direct relevance to frequency protection schemes and operational security in interconnected power systems. While additional cyber-resilience metrics, such as observability degradation and control signal distortion, are recognized as important complementary indicators, their detailed investigation is beyond the scope of the present work and is identified as a direction for future research. The mathematical representation of the resonance-based cyber-attack is given by:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{L2}= -0.1.sign\left[\Delta {f}_{2}{\prime}\right]$$\end{document}

which relates the injected load disturbance to the frequency variation in the second control area, consistent with previously reported models^43^.

Design of proposed (1 + FOPI)-FOPI-TID controller structure

As shown in Table 1, the literature controller structures include two-stage configurations such as PIDF-(1 + FOD), PIFOD-(1 + PI), and PDF-(1 + PI), where the $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ACE$$\end{document}$ is the sole input. In contrast, the PD-PI controller utilizes both $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ACE$$\end{document}$ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta F$$\end{document}$ as inputs; however, $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta F$$\end{document}$ is directly fed into the secondary PI controller, which limits disturbance rejection. Additionally, the tri-stage (1 + PI)-PI-PID controller incorporates $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ACE$$\end{document}$ , $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta F$$\end{document}$ , and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{\mathrm{tie}}$$\end{document}$ , but retains a conventional PID structure, with $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{\mathrm{tie}}$$\end{document}$ passing through only one controller, the PID. These controller structures are shown in Fig. 4. A comparative analysis of controller attributes provided in terms of controller structure, number of stages, scale of PSI, and limitations and advantages are summarized in Table 1.Table 1. Controller structure, scale, input, limitations, and advantages.RefScaleCC structureInput SignalsLimitationsAdvantages^44^Maritime microgridPI-(1 + PD) $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta F$$\end{document}$ - PID controller and conventional structure- It depends on ACE or $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \mathrm{F}$$\end{document}$ only- Reduces high-frequency noise better than single PID- Better stability and faster disturbance rejection^45^Two-area hydrothermal IPS in deregulated environmentPIDF-(1 + FOD)ACE^46^Three-area system with GRC &FACTSPDF + (1 + PI)ACE^47^Two-area of different test modelsPD-PIACE, and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta F$$\end{document}$ - $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \mathrm{F}$$\end{document}$ only passes to PI controller, losing disturbance rejection since it does not pass through two controllers- Using PID controller- Fast stabilization with minimal peaks- Enhances robustness and flexibility in high-variability conditions^25^Two area: reheat, thermal, and hydropower PSI(1 + PI)-PI-PIDACE, $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta F$$\end{document}$ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{\mathrm{tie}}$$\end{document}$ - Using PID controller- $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {\mathrm{P}}_{\mathrm{tie}}$$\end{document}$ only passes to PID controller, losing disturbance rejection as it does not pass through two controllers^23^Two-area microgrid with tidal PSPIFOD-(1 + PI)ACE-Tested on small PSs and one input ACE- Cascade controller attain good transient responseFig. 4The literature controller structure: (a) PIDF-(1 + FOD), (b) PIFOD-(1 + PI), (c) PDF-(1 + PI), (d) PD-PI, and (e) (1 + PI)-PI-PID.

This paper introduces a novel controller structure of (1 + FOPI)-FOPI-TID, as illustrated in Fig. 5. This configuration allows for the swift rejection of disturbance sources before they can propagate through the system. The unique controller incorporates two input signals: the main $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ACE$$\end{document}$ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta F$$\end{document}$ . This feature distinguishes it from other CC configurations in the field, as the $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ACE$$\end{document}$ serves as the primary input to the (1 + FOPI) controller, in parallel, the $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta F$$\end{document}$ is the input to the FOPI controller, and the combined output of both feeds into the TID controller.Fig. 5. Proposed structure of (1 + FOPI)-FOPI-TID CC.

Although the $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ACE$$\end{document}$ already contains frequency deviation information ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ACE={\Delta P}_{tie}+ {\Delta F}_{i}$$\end{document}$ ), using both $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta F}_{i}$$\end{document}$ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ACE$$\end{document}$ as inputs in the proposed controller is intentional and meaningful. The two signals represent different aspects of system behavior. The frequency deviation $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta F}_{i}$$\end{document}$ reflects fast local changes caused by sudden load disturbances and is useful for taking quick corrective action. On the other hand, the $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ACE$$\end{document}$ signal represents the overall control objective, as it includes both frequency deviation and tie-line power deviation, and is responsible for coordinating power exchange between areas and restoring steady-state operation. By assigning $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta F}_{i}$$\end{document}$ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ACE$$\end{document}$ two different stages of the CC, the proposed structure separates fast local response from slower inter-area regulation. Even though these signals are mathematically related, they act on different time scales, and their combined use does not lead to conflicting control actions. Instead, this structured design improves disturbance rejection, damping, and robustness, especially in the presence of nonlinearities and uncertainties.

Therefore, the proposed controller is not unnecessarily complex but rather designed to handle the practical dynamic characteristics of interconnected power systems more effectively. The S-domain representations of the (1 + FOPI), FOPI, and TID controllers are as follows:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${G}_{1+FOPI}\left(\mathrm{s}\right)=1+{K}_{p1}+ {K}_{i1}{s}^{-\lambda 1}$$\end{document}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${G}_{FOPI}\left(\mathrm{s}\right)={K}_{p2}+ {K}_{i2}{s}^{-\lambda 2}$$\end{document}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G_{{TID}} \left( s \right) = K_{t} s^{{ - 1_{{/_{n} }} }} + \frac{{K_{{i3}} }}{S} + K_{d} s$$\end{document}

where $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{pi}$$\end{document}$ , $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{i}$$\end{document}$ , $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{t}$$\end{document}$ , $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{d}$$\end{document}$ are proportional, integral, tilt, derivative, while $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document}$ , and n are the fractional parameters of integral and tilt factors. Achieving optimal performance for this controller requires the simultaneous optimization of these ten parameters. To fulfill this objective and attain fair comparison, a WAA is selected for this study to tune the proposed and the literature controller structures.

In this study, The FO operators were realized using the Oustaloup recursive approximation method, which provides accurate finite‑dimensional representations of fractional differ‑integrators over the selected frequency band [0.01–100] rad/s. This choice ensures numerical stability and efficient implementation in MATLAB/Simulink. The controller parameters are optimized offline, and the resulting implementation requires only low-order filters and simple algebraic operations, making the proposed control scheme computationally efficient and suitable for real-time deployment in large-scale interconnected power systems, even in the presence of communication delays.

The proposed (1 + FOPI)–FOPI–TID CC architecture is motivated by the physical structure and different response time requirements of automatic generation control rather than mathematical complexity alone. Each stage is assigned a distinct control role: the (1 + FOPI) loop, driven by the area control error, ensures accurate steady-state regulation and coordinated inter-area power exchange, with fractional-order integration enhancing robustness to uncertainties and low-frequency oscillations. The FOPI loop directly processes local frequency deviation, enabling rapid disturbance rejection and improved damping under sudden load changes and nonlinear operating conditions. The final TID stage functions as a dynamic compensator, where the tilt component provides wide-band damping, the integral term removes residual steady-state offsets, and the derivative term improves transient shaping without excessive noise amplification. This cascade configuration deliberately separates fast local frequency control from slower inter-area regulation, reducing control interaction and avoiding aggressive gain tuning required in simpler single- or dual-loop fractional controllers. As a result, the proposed structure closely reflects the physical control objectives of interconnected power systems, thereby enhancing robustness, interpretability, and practical applicability beyond purely performance-based validation.

Performance index of LFM

Customer demand changes affect the PS’s frequency, which causes $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{i}$$\end{document}$ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{tie}$$\end{document}$ to deviate from their nominal values. With an appropriate Performance function, it is advised to first successfully use contemporary heuristic optimization-based controller techniques in order to improve LFM performance^39,48,49^. The ISE, ITAE, IAE, and ITSE are some of the common performance functions that have been literature to address AGC issues and restore nominal performance^43,50^.

The ITAE criterion was chosen as the primary objective because of its proven effectiveness in minimizing sustained deviations and improving damping in AGC studies, thereby ensuring faster convergence and reduced steady‑state error^51^. ITAE provides a balanced trade‑off between transient suppression and long‑term stability, which aligns with the focus of this work. The application of the ITAE criterion will be essential in this investigation in order to modify the (1 + FOPI)-FOPI-TID CC’s scaling parameters. This objective is represented by the symbol ξ as follows:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi =\mathrm{ITAE}={\int }_{0}^{{\mathrm{t}}_{\mathrm{sim}}}\mathrm{t}\cdot \left[{\Delta F}_{1} +{\Delta F}_{2} +\Delta {P}_{\mathrm{tie}}\right]\mathrm{dt}$$\end{document}

In order to find out optimistic controller gains on a large scale, the gains range are chosen in the range $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left[0, 2\right]$$\end{document}$ as in^24,36^ for non-linear MSTA model. The n and λ in the range $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left[2, 3\right]$$\end{document}$ ^52,53^ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left[0, 2\right]$$\end{document}$ ^11^ respectively. In conclusion, by exploiting the WAA metaheuristic algorithm, optimal or near optimal set of (1 + FOPI)-FOPI-TID CC parameters is derived within the recommended edges for the smallest value of $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi .$$\end{document}$

Weighted average algorithm (WAA)

The two main steps of population-based meta-heuristic algorithms are exploitation and exploration. Exploration is concerned with avoiding local optima in the search space (SS), whereas exploitation is concerned with bettering solutions found during exploration. This approach is best shown by the WAA, which strikes a balance between exploration and exploitation. To keep this equilibrium, it uses movement strategies based on random constants and iteration numbers and calculates the population’s weighted average position (WAP) at each iteration.

The weighted average position’s inspiration

Each element is given a certain weight in order to calculate the WAP; the final average is more affected by the higher weights. The SS directs the updating of each person’s location for the subsequent iteration during each iteration. Individual positions, their personal best position (PBP), and the global best position (GBP) form the basis of the new SS. The global position is the best fitness value for the entire population, whereas a person’s PBP is the maximum fitness value attained during its movement history. Together with the PBP and GBP, the WAP is used to update future positions and represents the current population distribution.

Algorithm and mathematical model

Matrix S, which represents a collection of potential solutions, is used to randomly begin the optimization process:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S=\left[\begin{array}{c}\begin{array}{c}{S}_{1}\\ {S}_{2}\end{array}\\ \begin{array}{c}\vdots \\ {S}_{i}\\ \vdots \end{array}\\ {S}_{n}\end{array}\right]=\left[\begin{array}{cccccc}{s}_{1}^{1}& {s}_{1}^{2}& \cdots & {s}_{1}^{j}& \cdots & {s}_{1}^{d}\\ {s}_{2}^{1}& {s}_{2}^{2}& \cdots & {s}_{2}^{j}& \cdots & {s}_{2}^{d}\\ \vdots & \vdots & \cdots & \vdots & \cdots & \vdots \\ {s}_{i}^{1}& {s}_{i}^{2}& \cdots & {s}_{i}^{j}& \cdots & {s}_{i}^{d}\\ \vdots & \vdots & \cdots & \vdots & \cdots & \vdots \\ {s}_{n}^{1}& {s}_{n}^{2}& \cdots & {s}_{n}^{j}& \cdots & {s}_{n}^{d}\end{array}\right], \mathrm{where} \left\{\begin{array}{c}i=1, 2, \dots ,n.\\ j=1, 2, \dots ,d.\end{array}\right.$$\end{document}

where $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d$$\end{document}$ is the dimension of the deemed problem, $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n \text{is the number}$$\end{document}$ of solutions candidate along the SS, and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{i}$$\end{document}$ is the $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${i}^{th}$$\end{document}$ candidate solution, , ^30^. The starting positions $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${s}_{i}^{j}$$\end{document}$ are disclosed randomly as:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${s}_{i}^{j}={LB}_{j}+rand.\left({UB}_{j}-{LB}_{j}\right)$$\end{document}

where $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${LB}_{j}$$\end{document}$ refer to the $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${j}^{th}$$\end{document}$ lower bound value, $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$rand$$\end{document}$ is a random number within [0: 1] range, and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${UB}_{j}$$\end{document}$ expresses the $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${j}^{th}$$\end{document}$ upper bound value.

Weighted average position (WAP)

First, each participant’s fitness is determined, then the population is rearranged according to the cost function to determine the WAP. The first N candidates in the population are then chosen in order to compute the WAP in the manner described below:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${N}_{candidate}=\left(nP-4\right)\frac{it-1}{1-{Max}_{It}}+nP$$\end{document}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${SUM}_{Fitness}=\sum_{i=1}^{{N}_{candidate}}Fitness ({S}_{i})$$\end{document}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${X}_{miu}=\frac{\sum_{i=1}^{{N}_{candidate}}{S}_{i}({SUM}_{Fitness}-Fitness \left({S}_{i}\right))}{{SUM}_{Fitness}({N}_{candidate}-1)}$$\end{document}

where $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{i}$$\end{document}$ is the i^th^ candidate solution, $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${SUM}_{Fitness}$$\end{document}$ is the total fitness values of the chosen candidates, $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$nP$$\end{document}$ is the populations number, $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Fitness$$\end{document}$ is the function used to calculate the fitness value, $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${N}_{candidate}$$\end{document}$ is the number of chosen candidate solutions, $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${N}_{candidate}$$\end{document}$ is the applicant solutions number, $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${X}_{miu}$$\end{document}$ is WAP. Besides, $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$it$$\end{document}$ is the current iteration number, and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Max}_{It}$$\end{document}$ is the limit number of iterations.

Exploitation or exploration process

The search phase for the potential solution, i.e., exploration or exploitation, should be identified during the iteration process:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f\left(it\right)=\left(\alpha .rand-1\right)\mathrm{sin}\left(\pi \frac{it}{{Max}_{It}}\right)$$\end{document}

The balance between exploration and exploitation phases is controlled by the constant $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}$ . An equilibrium is reached if $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f\left(it\right)= 0.5$$\end{document}$ . Candidate solutions with $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f\left(it\right)\ge 0.5$$\end{document}$ move in accordance to their exploitation capabilities, whereas those with $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f\left(it\right)<0.5$$\end{document}$ move in accordance to their exploration abilities^30^

(b)Exploitation phase

The exploitation strategy in the WAA directs search agents towards areas likely to yield new global best values. To efficiently exploit the SS, the WAA employs three movement strategies that incrementally move search agents closer to the optimal solution as shown in Fig. 6.Fig. 6. Diagram of search agent position updating based on exploitation strategies.

(i)The first strategy exploits the SS using the WAP of the entire population, PBP, and GBPs.

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{i}\left(it+1\right)={r}_{11}.\left({S}_{miu}\left(it\right)-{S}_{G{lobal}_{best}}\left(it\right)\right)+{r}_{12}.\left({S}_{miu}\left(it\right)-{S}_{{Personal}_{best}}\left(it\right)\right)+{r}_{13}.{S}_{miu}\left(it\right)$$\end{document}

where $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${r}_{ij}$$\end{document}$ are random values within [0–1], to modify the expansion of the search space around the personal best $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{{Personal}_{best}}$$\end{document}$ and GBP s $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{G{lobal}_{best}}$$\end{document}$ to express the PBP and GBP in the iteration numbers ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$it)$$\end{document}$ .

(ii)The second technique ignores the GBP and concentrates on the region between the PBP and the population’s WAP in order to narrow the SS. This method, in contrast to the first technique, prioritizes a narrower area for exploitation. Using the following, the new SS is produced^34^:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{i}\left(it+1\right)={r}_{21}.\left({S}_{miu}\left(it\right)-{S}_{{Personal}_{best}}\left(it\right)\right)+{r}_{22}.{S}_{{Personal}_{best}}\left(it\right)$$\end{document}

This approach shows a greater bias towards improving accuracy and convergence rate by guiding movement from the agent’s PBP towards the population’s WAP.

(iii)The third trick deals with the SS between the GBP and the population’s WAP. Each candidate’s position update is determined for this purpose using the formula below:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{i}\left(it+1\right)={r}_{31}.\left({S}_{miu}\left(it\right)-{S}_{{Global}_{best}}\left(it\right)\right)+{r}_{32}.{S}_{{Global}_{best}}\left(it\right)$$\end{document}

This approach has a narrow SS that concentrates on the area between the GBP and the WAP. Both the accuracy and the pace of convergence are improved by this focused movement^34^.

(c)Exploration phase

The exploration of the solution space utilizes two strategies: the Levy Flight method and dynamic movement distances. The Levy flight method focuses on the GBP, creating a broad search area around the WAP. This involves short movements for local optimization and longer movements to explore distant regions. Candidates dynamically move in various directions from the GBP, enhancing exploration and balancing the exploitation of current solutions. The Levy flight model follows a Levy distribution, characterized by mostly small strides with occasional larger ones, reflecting natural patterns observed in birds and insects, expressed as^30^:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X=\frac{U}{|V{|}^{1/\beta }}, U=normal\left(0, {\sigma }_{u}^{2}\right)\& V=normal\left(0, {\sigma }_{v}^{2}\right)$$\end{document}

where $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X$$\end{document}$ is the step length of the Levy flight influenced by the parameter β, and Γ stands for the Gamma distribution function. Besides, V and U satisfy normal distributions with standard deviations equal to $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\sigma }_{v}$$\end{document}$ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\sigma }_{u}$$\end{document}$ respectively, and zero mean values. The first exploration technique, which is based on a controlled adjustment of step length X, is as follows:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{i,j}\left(it+1\right)={S}_{{Global}_{best}}\left(it\right)+X$$\end{document}

where $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{i,j}\left(it+1\right)$$\end{document}$ denotes the $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$j$$\end{document}$ ^th^ position of the $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i$$\end{document}$ ^th^ solution at the next iteration. In certain instances, the algorithm’s global optimum might be located close to the GBP yet too far from the ideal GBP value. When using the first method in such circumstances, the WAA runs the danger of convergent to a local optimum. The second movement method is used to modify the SS in order to reduce this danger^34^:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{i}(it+1)={LB}_{min}+rand.\left({UB}_{min}-{LB}_{min}\right)$$\end{document}

where the lower and upper boundaries’ minimum values are denoted by $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${LB}_{min}$$\end{document}$ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${UB}_{min}$$\end{document}$ , respectively. By using this tactic, the search agents can shift to different locations within the SS in an effort to locate possibly superior areas.

The WAA system optimizes gains in (1 + FOPI)-FOPI-TID controllers using the ITAE criteria. It iteratively generates new gains by minimizing ITAE values, aiming for the lowest ITAE to reduce system deviations. Research indicates that larger initial populations do not necessarily outperform smaller ones in finding optimal solutions, making the approach suitable for real-time LFM scenarios^54^. The operational framework of the proposed WAA is depicted in the flowchart presented in Fig. 7, which outlines the sequential logic and decision-making steps involved in the algorithm.Fig. 7WAA flowchart.

Results and discussion

This study investigates the effectiveness of a novel (1 + FOPI)-FOPI-TID controller using a combination of hybrid power sources and non-linearities, GRC, GDB, BD, and CTD, alongside cyber-attacks. To achieve the LFM goal, the controller coefficients are optimized with the WAA optimizer in MATLAB, integrated with Simulink on an Intel Core i-5 processor. The WAA-based controller settings for the cases study after 100 iterations of optimization are conducted. The robustness of the proposed controller is compared to other regulators, including PIDF(1 + FOD), PIFOD-(1 + PI), (1 + PI)-PI-PID CC, PDF-(1 + PI), and PD-PI CCs, all under fixed load changes of 1% SLP. Results from the analyzed multiarea PSI are thoroughly examined in subsequent case studies.

Case 1: Effectiveness of the WAA algorithm

To demonstrate the adaptability of the WAA, the convergence characteristics of WAA, mayfly algorithm (MA)^43^, and grey wolf optimization (GWO)^55^ are analyzed in the context of the (1 + FOPI)-FOPI-TID CC using a test model as in Fig. 1. Each optimizer is configured with 100 iterations. The convergence graph of WAA, extracted from MSTA with GDB, GRC, BD, and CTD nonlinearities is depicted in Fig. 8, showcasing WAA’s swift convergence and effectiveness in enhancing AGC control strategies to optimize system performance. Furthermore, a comparison of the minimum, maximum, average, and standard deviation of the ITAE values, $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\xi }_{ITAE}$$\end{document}$ , in Table 2 illustrates that WAA outperforms MA and GWO in addressing AGC challenges.Fig. 8. The convergence characteristic curves of PSO, GWO, and the proposed FHO.Table 2ITAE values for WAA, MA, and GWO at case1.AlgorithmMinimumMaximumAverageSt. deviationWAA0.0138100.0154430.0142870.00085502MA0.0161090.0228110.0188630.00336873GWO0.0169670.0230630.0223990.00334436

Case 2: Analyses of controller’s performance

This study conducted a comprehensive performance comparison between the proposed (1 + FOPI)-FOPI-TID controller and various two-stage and three-stage control approaches, including PIDF(1 + FOD), PIFOD-(1 + PI), (1 + PI)-PI-PID CC, PDF-(1 + PI), and PD-PI CC under WAA scheme. The gains of these controllers are listed in Table 3. The evaluation metrics focused on the ITAE value, minimum undershoot (Ush), and settling time (T_s_), considering changes in $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{1}$$\end{document}$ , * $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{tie}$$\end{document}$ , and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{2}$$\end{document}$ . Overall, Table 4 compares the performance of each optimization procedure using transient metrics. Figure 9 illustrates the dynamic responses of a two-area interconnected power system subjected to a load disturbance, presenting the frequency deviations in Area-1 ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{1}$$\end{document}$ ), Area-2 ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{2}$$\end{document}$ ), and the tie-line power deviation ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{tie}$$\end{document}$ ) for various WAA-based cascade controller structures. A comparative evaluation is carried out between the proposed WAA-(1 + FOPI)–FOPI–TID controller and other existing controllers, including PIDF-(1 + FOD), (1 + PI)–PI–PID, PIFOD-(1 + PI), PDF-(1 + PI), and PD–PI.Table 3WAA optimized controller gains under 1% SLP in area 1.WAA based CC structureArea 1Area 2 $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{K}}}_{{\boldsymbol{P}}}/{{\boldsymbol{K}}}_{{\boldsymbol{t}}}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{K}}}_{{\boldsymbol{I}}}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{K}}}_{{\boldsymbol{D}}}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\mu}}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\uplambda}}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{n}}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{N}}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{K}}}_{{\boldsymbol{P}}}/{{\boldsymbol{K}}}_{{\boldsymbol{t}}}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{K}}}_{{\boldsymbol{I}}}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{K}}}_{{\boldsymbol{D}}}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\mu}}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\uplambda}}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{n}}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{N}}$$\end{document}$ PIDF (1 + FOD)PIDF0.63332.00001.3e-4––––1.50190.14952.0000––––(1 + FOD)–––0.2611––23.53–––1.0e-6––0.0126WAA: PIFOD-(1 + PI)PIFOD0.65832.00007.02e-050.9757–––2.00000.27690.38160.3802–––(1 + PI)2.00001.129e-7–––––0.05270.00373–––––PD-PI CCPD2.0000–0.4935––––2.0000–2.0000––––PI0.27091.9308–––––0.05690.0873–––––(1 + FOPI)-FOPI-TID CC(1 + FOPI)1.27772.0000––1e-5––2.00000.0786––1.9977––FOPI0.67310.0769––0.9628––1.0e-40.0031––0.4319––TID1.99951.63990.2435––2.5071–1.0e-60.00000.0000––2.999–(1 + PI)-PI-PID CC*1 + PI1.69191.8777–––––1.67521.2e-05–––––PI0.32191.9545–––––0.35821.3421–––––PID0.88230.00000.1154––––1.24191.08610.1967––––*PDF-(1 + PI)*PDF2.0000–0.2135–––0.21350.0311–0.5343–––100(1 + PI)0.00002.0000–––––0.00002.0000–––––Table 4. Comparison performance of several controllers for case 2.WAA based: $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi$$\end{document}$ ST (0.0002 band) (s)Undershoot(-ve)ITAE $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{1}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{2}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta P$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{1}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times {10}^{-3}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{2}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times {10}^{-3}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta P$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times {10}^{-4}$$\end{document}$ (1 + FOPI)-FOPI-TID0.013811.7071.7351.35513.01.62114.15PIFOD-(1 + PI)0.025263.2643.3931.72313.02.31616.22PD-PI0.026443.4942.5822.61713.01.87115.62PIDF(1 + FOD)0.031704.4603.8902.86513.02.41616.33(1 + PI)-PI-PID0.033904.1574.9154.17713.01.67515.34PDF-(1 + PI)0.034934.6612.6102.85813.02.38816.17% Reduction**45.3%****47.7%****32.8%21.4%—****3.2%****7.8%**Significant values are in [bold].Fig. 9. Transient response with various controllers but same optimizer at case 2.

As demonstrated in Fig. 9 and summarized in Table 4, the proposed (1 + FOPI)-FOPI-TID control scheme outperformed all comparative techniques. Table 4 presents a comparative performance evaluation of several WAA-based cascade controllers for Case 2 using key dynamic performance indices, including ITAE, settling time (within a ± 0.0002 band), and maximum undershoot of frequency and tie-line power deviations.

Notably, it achieved an ITAE value of 0.01381—representing a 45.3% improvement over the PIFOD-(1 + PI) controller. The scheme also showed significant enhancements in settling times: 47.7% improvement for $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{1}$$\end{document}$ , 21.4% for $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{tie}$$\end{document}$ , both compared to PIFOD-(1 + PI) CC, and 32.8% improvement for $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta$$\end{document}$ F₂ (compared to PD-PI CC). While no improvement was observed in the percentage undershoot for $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{1}$$\end{document}$ , the proposed controller achieved reductions of 3.2% for $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{2}$$\end{document}$ and 7.8% for $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{tie}$$\end{document}$ when compared to the (1 + PI)-PI-PID CC.

As a result, the proposed (1 + FOPI)-FOPI-TID CC appears to be a potential solution for improving LFM performance in various PSIs. The transient responses of the proposed WAA: (1 + FOPI)-FOPI-TID CC and other controllers were further examined under 4% SLC in area 1, as shown in Fig. 10. The WAA: (1 + FOPI)-FOPI-TID CC outperformed all competing controllers under 4% SLC in area 1, confirming its superiority and demonstrating MSTA stability with nonlinearities.Fig. 10. Transient reactions for case 2 with 4% SLC at area 1: (a) Δ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${F}_{1}$$\end{document}$ , (b) Δ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${F}_{2}$$\end{document}$ , and (c) ΔP_tie_.

Case 3: Performance analyses under random load perturbation

The dynamic analysis of the proposed WAA: (1 + FOPI)-FOPI-TID CC was conducted using a random step load (RSL) in area 1. The RSL perturbation involves a sequence of sudden generating unit outages or abrupt load switches, representing a mix of random load variations as depicted in Fig. 11. Figure 12. depicts the transient responses of a two-area interconnected power system under successive load disturbances, illustrating the frequency deviations in Area-1 ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{1}$$\end{document}$ ), Area-2 ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{2}$$\end{document}$ ), and the tie-line power deviation ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{tie}$$\end{document}$ ) for different WAA-based cascade controller schemes. The highlighted zoomed-in regions emphasize the initial disturbance intervals, allowing a clear comparison of dynamic performance.Fig. 11. Random step loads $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta P}_{L1}$$\end{document}$ in area 1.Fig. 12. Transient response with various algorithm techniques.

It is apparent from Fig. 12 that, when compared to other controllers in the existing literature, the WAA: (1 + FOPI)-FOPI-TID CC demonstrates consistent dynamic performance under randomized loading patterns in area 1. It shows reduced undershoot and overshoot, along with faster convergence, highlighting its robustness and effectiveness in handling dynamic variations Fig. 13.Fig. 13(a): Variation of $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{g}$$\end{document}$ , $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B$$\end{document}$ , $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R$$\end{document}$ , and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{12}$$\end{document}$ PS parameters for Δ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${F}_{1}$$\end{document}$ . (b): Variation of $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{g}$$\end{document}$ , $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B$$\end{document}$ , $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R$$\end{document}$ , and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{12}$$\end{document}$ PS parameters for Δ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${F}_{2}$$\end{document}$ . (c): Variation of $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{g}$$\end{document}$ , $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B$$\end{document}$ , $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R$$\end{document}$ , and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{12}$$\end{document}$ PS parameters for ΔP.

Case 4: sensitivity analysis/robustness

A sensitivity analysis was conducted to evaluate the robustness of the WAA optimized (1 + FOPI)-FOPI-TID CC. The stability of the proposed controller was assessed by comparing the responses of the original parameters to various metrics when parameters such as $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{g}$$\end{document}$ *, * $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B$$\end{document}$ *, * $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R$$\end{document}$ , and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{12}$$\end{document}$ are varied by approximately ± 25%, as outlined in Table 5 under $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {P}_{L1}=1\%$$\end{document}$ SLP at area 1. Table 5 and Fig. 14 demonstrate the suggested controller’s reliability under altered system settings. The impact of parameter fluctuations on the ITAE criterion, settling time, and undershoot achieved by the proposed approach under modifications of system parameters was found to be negligible. The analysis reveals that even with significant variations of up to 25% in parameters, the proposed scheme exhibits robustness and maintains excellent performance, with minimal effects on the key evaluation metrics. This underscores the robust nature of the controller and its ability to perform effectively under varying system conditions.Table 5. Percentage variations in hybrid PS parameters.Parameters%of ChangeITAESettling time (0.0002 band) T_s_ (s)Undershoot(-ve) $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{1} \Delta {F}_{2} \Delta P$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{1}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times {10}^{-3}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{2}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times {10}^{-3}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta P$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times {10}^{-4}$$\end{document}$ $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{g}$$\end{document}$ + 250.014811.9242.3121.35513.01.80414.09Nominal0.013811.6931.7521.35513.01.62514.13− 250.016912.1692.3332.37613.01.68714.87 $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B$$\end{document}$ + 250.022493.8832.4322.33313.01.71313.49Nominal0.013811.6931.7521.35513.01.62514.13− 250.020273.5951.4792.16313.01.74415.60 $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R$$\end{document}$ + 250.016612.8311.9081.31213.01.86014.31Nominal0.013811.6931.7521.35513.01.62514.13− 250.017362.7882.3332.46113.01.61214.50 $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{12}$$\end{document}$ + 250.016523.7822.3541.22713.02.17816.98Nominal0.013811.6931.7521.35513.01.62514.13− 250.019623.9702.4543.01413.01.32011.70Fig. 14Cyber-attack signal model applied to the PS.

Case 5: Resonance Cyber Attack (ResA)

In this case, the ResA is simulated by introducing load disturbances at arbitrary frequency values, without any prior knowledge of the system’s operational state. Figure 14. illustrates a cyber-attack injection signal applied to the power system model, represented as a bounded high-frequency square-wave disturbance varying between ± 0.1 p.u. over the simulation duration. This signal is used to emulate false data injection or malicious control signal manipulation, periodically perturbing system measurements or control inputs to evaluate the robustness of the load frequency control scheme under cyber-physical threats.

The cyber-attack signal is synthetically generated within the simulation environment using a periodic switching (square-wave) function with fixed amplitude and duty cycle.

The dynamic response of the system to this attack is analyzed using various control strategies, with corresponding results illustrated in Fig. 15. Among the compared controllers, the WAA: (1 + FOPI)-FOPI-TID configuration demonstrates the most robust performance, exhibiting the lowest RoCoF amplitude in area 1 (0.0654) and maintaining a relatively moderate amplitude in area 2 (1.3276) summarized in Table 6. This indicates superior disturbance rejection and frequency stability under cyber-attack conditions.Fig. 15. Controller performances under the cyber-attack: ΔF_1_, ΔF_2_, and ΔP_tie._Table 6RoCoF (Hz/s) for different controllers.ControllerMax. Amplitude (Area 1)Max. Amplitude (Area 2)WAA: (1 + FOPI)-FOPI-TID0.0654****1.3276WAA: PD-PI0.09671.3835WAA: (1 + PI)-PI-PID0.11581.7025WAA: PIFOD-(1 + PI)0.09301.5391WAA: PDF-(1 + PI)0.08801.7652PIDF(1 + FOD)0.09591.5442Significant values are in [bold].

The WAA: PDF-(1 + PI) controller ranks second, with a slightly higher RoCoF amplitude in area 1 (0.0880) and a comparable value in area 2 (1.7652), suggesting effective damping characteristics. In contrast, the (1 + PI)-PI-PID controller performs the poorest, recording the highest RoCoF amplitude in area 2 (1.7025), alongside a relatively elevated value in area 1 (0.1158). Similarly, the WAA: PIFOD-(1 + PI) controller shows inadequate performance in area 2 with an amplitude of 1.5391, further emphasizing its reduced resilience to cyber-attacks.

The RoCoF profiles, presented in Fig. 1, align with these observations, while most controllers effectively limit RoCoF within acceptable thresholds in area 1 and area 2 exhibits higher susceptibility to instability.

Conclusion

This study presents a novel (1 + FOPI)-FOPI-TID cascade control strategy, optimized using a WAA, for the frequency regulation of heterogeneous two-area PSIs incorporating reheat thermal, gas, hydro, and nuclear generation units. The system model integrates critical nonlinearities, including GRC, BD, GDB, and CTD, as well as cyber-attack scenarios to evaluate real-world robustness. Extensive dynamic simulations, including random step load perturbations, are conducted to assess the effectiveness of the proposed control scheme. Comparative evaluations with contemporary control configurations, such as WAA: PIFOD-(1 + PI) and WAA: PD-PI CCs, demonstrate the superior dynamic performance of the proposed method. Specifically, the (1 + FOPI)-FOPI-TID controller achieved a 45.3% improvement in ITAE and reduced frequency settling times by 47.7% ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{1}$$\end{document}$ ) and 32.8% ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {F}_{2}$$\end{document}$ ) in the first and second areas, respectively, under test case 2. Moreover, sensitivity analysis confirms the robustness of the proposed controller, maintaining consistent performance despite substantial parameter variations. Under cyber-attack scenarios, the controller exhibited minimal RoCoF amplitudes in both areas, further underscoring its resilience Fig. 16.Fig. 16. RoCoF results of PS areas for cyber-attack for Area-1 and Area-2.

In summary, while the fractional-order cascade structure is more complex than classical controllers, its fixed-parameter implementation post-tuning is feasible for real-time application, and the design-time optimization cost is justified by the achieved performance and robustness benefits. This research contributes significantly to the advancement of load frequency control strategies through the development of robust, WAA-based CC architecture. The WAA: (1 + FOPI)-FOPI-TID approach not only enhances system stability and dynamic responsiveness but also ensures operational resilience under nonlinearities and cybersecurity threats.

Future work will focus on implementing the controller on real-time hardware (like OPAL-RT) to empirically validate its computational performance, alongside the integration of fault-detection mechanisms capable of rapidly identifying and diagnosing malicious intrusions before they can destabilize the power system.

Bibliography19

The reference list from the paper itself. Each links out to its DOI / PubMed record.

1Daraz, A., Malik, S. A., Basit, A., Aslam, S. & Zhang, G. Modified FOPID controller for frequency regulation of a hybrid interconnected system of conventional and renewable energy sources. Fractal Fract.7 (1). 10.3390/fractalfract 7010089 (2023).
2Ali, T. et al. Terminal voltage and load frequency control in a real four-area multi-source interconnected power system with nonlinearities via OOBO algorithm. IEEE Access, (2024).
3Liu, S., Liu, X. P. & Saddik, A. E. Denial-of-Service (dos) attacks on load frequency control in smart grids, in 2013 ieee pes innovative smart grid technologies conference (isgt), 1–6. (2013).
4Ye, Y., Daraz, A., Basit, A., Khan, I. A. & Alqahtani, S. A. Cascaded Fractional-Order Controller-Based Load Frequency Regulation for Diverse Multigeneration Sources Incorporated with Nuclear Power Plant, Int. J. Energy Res., 10.1155/2024/7939416 (2024).
5Ogar, V. N., Hussain, S. & Gamage, K. A. A. Load frequency control using the particle swarm optimisation algorithm and PID controller for effective monitoring of transmission line. Energies 16 (15). 10.3390/en 16155748 (2023).
6Singh, K. & Arya, Y. Jaya-ITDF control strategy-based frequency regulation of multi-microgrid utilizing energy stored in high-voltage direct current-link capacitors. Soft Comput. Fusion Found. Methodol. & Appl, 27, 9, (2023).
7Shalaby, R., El-Hossainy, M., Abo-Zalam, B. & Mahmoud, T. A. Optimal fractional-order PID controller based on fractional-order actor-critic algorithm. 3 Springer Lond.3510.1007/s 00521-022-07710-7 (2023).
8Khan, I. A. et al. November., Load frequency control in power systems with high renewable energy penetration: A strategy employing PIλ(1 + PDF) controller, hybrid energy storage, and IPFC-FACTS, Alexandria Eng. J., 106, 337–366. 10.1016/j.aej.2024.06.087 (2024).