Novel Heuristic Optimization Technique to Solve Economic Load Dispatch and Economic Emission Load Dispatch Problems

: The fundamental objective of economic load dispatch is to operate the available generating units such that the needed load demand satisﬁes the lowest generation cost and also complies with the various constraints. With proper power system operation planning using optimized generation limits, it is possible to reduce the cost of power generation. To fulﬁll the needs of such objectives, proper planning and economic load dispatch can help to plan the operation of the electrical power system. To optimize the economic load dispatch problems, various classical and new evolutionary optimization approaches have been used in research articles. Classical optimization techniques are outdated due to many limitations and are also unable to provide a global solution to the ELD problem. This work uses a new variant of particle swarm optimization techniques called modiﬁed particle swarm optimization, which is effective and efﬁcient at ﬁnding optimum solutions for single as well as multi-objective economic load dispatch problems. The proposed MPSO is used to solve single and multi-objective problems. This work considers constraints like power balance and power generation limits. The proposed techniques are tested for three different case studies of ELD and EELD problems. (1) The ﬁrst case is tested using the data of 13 generating unit systems along with the valve point loading effect; (2) the second case is tested using 15 generating unit systems along with the ramp rate limits; and (3) the third case is tested using the economic emission dispatch (EELD) as a multi-objective problem for 6 generating unit systems. The outcomes of the suggested procedures are contrasted with those of alternative optimization methods. The results show that the suggested strategy is efﬁcient and produces superior optimization outcomes than existing optimization techniques.


Introduction
Since fuel resources (coal) are limited but the demand for electrical power is increasing, part of this research is to find the best solution to reduce the cost of power generation and reduce emissions so that the global warming problem also reduces.The above-discussed issues are the main motivation for considering this topic in our work [1].
The primary objective of the proposed study is to identify the ideal producing point where the power production system will function and fuel costs may be reduced.If the generating stations operate at the optimum generation point, the generation cost should be minimal for any particular load demand.The second motive is to reduce emissions; if the burning of fuel definitely reduces, the emission of toxic gases also decreases, reducing the pollution level [2,3].
Various classical approaches are listed for finding the optimum solution to the ELD problem, such as linear programming (LP) [3], pattern search PS [4], quadratic programming (QP) [5], etc. Classical methods have many drawbacks, like the LP and QP methods, which need linear mathematical objectives.Pattern search methods are sensitive to changes in the parameters, and the DS (direct search) method's performance for large data is not good.
The objective of the ELD problem is to maintain the generating output power according to the load demand and reduce the cost of power generation [1].ELD problems are also associated with constraints like power demand, generation limits, line losses, etc. Due to the presence of such constraints, economic load dispatch is a nonlinear computation problem.So, the optimization of such a problem requires efficient evolutionary techniques [2].Some classical and heuristic methods are not effective for solving nonlinear optimization problems because they have many limitations and take a long time to compute.Most of the modern heuristic techniques are more effective and provide the best solution to the problem [6].Heuristic optimizations techniques are the newest optimization techniques, as explained in various articles, and are very efficient, simple, and fast at computation.Such techniques effectively optimize the ELD data within seconds and cannot deviate in the direction of a local solution [7].
The latest optimization techniques have many features; they can optimize big data and provide a global solution without deviation from the constraints.These methods are unaffected by initial conditions or changes in the variables [8].Optimization is the task of finding the most suitable solution to a problem.The optimization process starts randomly to find the solution of the objective function (constraints are also defined) in the search space within the defined limits and with a number of agents, or swarms, and finally finds the solution to the problem, which is the global optimum solution, by exploring as small a set of solutions as possible [9].
The proposed work considered a modified PSO (MPSO), which is the latest variant of particle swarm optimization and is very effective for single objective (ELD) and multiobjective (EELD) problems.This is a new variant of the PSO that has an attraction factor that helps the particle move in the direction of a solution; it controls the movement of particles in such a way that they cannot deviate from the search area.
Other variants of the PSO need to update the velocity and position of the particle continuously until they can find a solution to the problem.This work proposes a new modified PSO, where the velocity of the particles is controlled by the attraction factor in such a way that they cannot stop during the iteration process; hence, there is no need to update the velocity of the particles; only the position of the particles is required to be updated.That is the big advantage of the proposed MPSO.Hence, the speed of computation increases and the particles are also not moved to the local solution.The MPSO always provides a global solution to the proposed ELD and EELD problems.This paper's remaining sections are structured as follows: the second part demonstrates the works published by the researchers.Here, we considered the last 20 years' papers for the study.The third part of the article presents the mathematical formulation of the economic load dispatch and economic emission dispatch.The fourth section of the article presents the mathematical model and algorithm of the new PSO.The fifth section shows the different case studies and the results of the new PSO compared with other optimization techniques.In the last section, the conclusion of the work is provided.

Literature Review
The economic load dispatch is formulated and optimized with various constraints using different optimization techniques to find out the minimum generation cost and fulfill the load demand.This section discusses the various optimization techniques used for finding the optimum solution to the ELD problem.
To address the ELD issues, a biogeography-based optimization technique was proposed.The recommended approach primarily uses the two processes of migration and mutation to obtain the global optimum [1].A new evolutionary technique is suggested for the optimization of multi-objective ELD problems [2].
The linear programming method is used to obtain the real and reactive power of the electrical generation system; however, such methods take a long computation time and are sometimes unable to provide a global solution for large data sets [3].
The pattern search method was proposed for finding the optimum solution to the ELD problem and the valve loading effects were also considered.The proposed algorithm was tested on various test data for validation of the results and compared with other optimization techniques [4].Quadratic programming was used for the solution of the ELD problem; this consists of a DC load flow and network security constraints [5].The PSO techniques used the ELD problem along with transmission losses, dynamic operation limits, and restricted operating zones [6].
Discrete and continuous ELD problem solutions are obtained by using a new variant of the PSO (EPSO).In this article, the ELD problem is formulated with penalties and constraints [7].
Article [8] solved the ELD problem, including the ramp rate limits and the transient stability, by using a novel particle swarm optimization called the chaotic quantum-behaved PSO algorithm.
The authors of [9] considered class-topping optimization techniques that are based on human intelligence.They solve the ELD problem as well as the CEED using four different case studies.
Article [10] proposes QPSO and HQPSO to solve the ELD problem.They solve a case ELD problem consisting of data from 13 generating unit systems along with the valve loading effect, and the results are compared to classical PSO.
To resolve the pollution and economic dispatch issues, a new combined genetic-tabu search technique is provided.The suggested technique is designed such that a basic GA serves as the first search to guide the search towards the ideal region [11].
Article [12] optimized the ELD problem along with the valve loading using the EPSO.The research also includes constraints like the prohibited operating zone.
The QPSO technique was used to find the optimum solution to the ELD problem.The suggested PSO particles move in the search area with good speed and explore the solution quickly [13].Two test cases of 10 and 15 generating units were taken and optimized for the ELD problem using the PPSO [14].Article [15] used a combination of differential evolution and a PSO called quantum particle swarm optimization (QPSO) as a solution for economic load dispatch optimization.
An artificial bee colony approach is discussed in the article to resolve the economic dispatch problem [16].The suggested technique is used to solve the economic dispatch problem along with the power balance and network loss constraints.
The BAT method is demonstrated for the optimization of multi-objective problems.This work considered the economic emission dispatch as a multi-objective problem and found the overall generation cost of the plant [17].The Chameleon Swarm Algorithm is proposed for the minimization of emissions as well as the generation cost of the power plant [18].
The ELD problem with a valve-point effect is solved using a hybrid method [19].A novel evolutionary optimization approach is proposed to resolve the economic load dispatch issue.The standard harmony search random selection process is replaced with a mutation process based on wavelet theory to increase the performance of the suggested strategy [20].
An optimization technique that is based on the teaching-learning approach is used to solve ELD problems without taking transmission losses into account.The suggested solution can handle an ELD while taking a nonlinearity, such as valve point loading, into account [21].A single objective ELD problem along with a prohibited zone and ramp rate and a multi-objective ELD problem are optimized using C andGPSO [22].
An effective evolutionary technique is using the novel adaptive particle swarm optimization for a limited ED issue [23].The differential evolution algorithm was investigated for the purpose of resolving power system economic load dispatch issues.Five ELD issues with various features were utilized to test the proposed methods [24].A contemporary search and rescue optimization approach that is motivated by the behavior of humans during search and rescue operations was used to solve the combined emission and economic dispatch [25].
To overcome the ELD difficulties, an enhanced arithmetic optimization is suggested.Two crucial variables in the original AOA are math optimizer acceleration and math optimizer probability [26].To optimize multi-objective ELDs, data mining technology is used.The proposed techniques provide the best decision for the sample test data [27].
The "Grey Wolf" approach is taken to solve the challenging nonconvex ELD optimization problems effectively and reliably.The grey wolf optimizer algorithm's search agents are helped by the sine and cosine functions in such a way that the solution cannot move to the local optima [28].
Harris Hawks' techniques are used to address the issues of economic load dispatch.The Harris Hawks Optimizer is used to find the number of possible solutions in the search space; each of these regions can be thoroughly searched for the best local solution via adaptive hill climbing [29].
An optimization algorithm based on the theory of demand and supply was used to optimize the ELD model.The proposed algorithm efficiently moves the particles in the search space between the local and global searches [30].
To address the complex economic load dispatch issue, a new hybrid grey wolf optimization method with a strong learning mechanism has been developed [31].The foraging optimization method is based on the Spiral Foraging approach given for the solution of the ELD problem.Spiral foraging strategies improve global search capability and convergence velocity [32].A multi-objective ELD problem was formulated using a wind energy system and a distributed generation system [33].
A hybrid generation system can help reduce the cost of power generation as well as carbon emissions [34].An adaptive multi-population-based differential evolutionary technique is proposed that improves ELD global solutions [35].
A diffusion model with heterogeneous coefficients and a generic nonlinear incidence rate for brucellosis was used for the solution of economic load dispatch [36].For the solution of economic load dispatch with emissions subject to power balance and generation limit constraints, a new particle swarm optimization technique was proposed [37].The model proposed is a nonlinear distributed delayed periodic AG-ecosystem with competition on the time scales.Our approach integrates and generalizes the discrete and continuous situations in terms of the time scale [38].
The Bumble Bee Mating Optimization algorithm is used for economic load dispatch optimization.The proposed BBMO works in three different modes: the queen, the workers, and the drones (males) [39].A multi-objective optimization method is used in this article for the optimal planning of the distributed generators in electric distribution networks.Reducing the annual cost and network loss, along with improving the reliability of the network, are the main objectives of a multi-objective algorithm [40].

Formulation of ELD Problem
The mathematical model for minimizing the cost of generation is given as

Formulation of Economic Dispatch along with Valve Point
When using a multi-valve steam turbine system, the fuel cost function of the generating units may vary [4].The cost of generation increases with the use of multi-valve steam turbines.The generator heat rate curve is the result of the multi-valve steam turbine valve opening operation.The impact of valve-point loading is shown in Figure 1 [22].In Figure 1,  A, B, C, D, and E are the points where the valve can be opened to control the speed of the generator if the load suddenly increases.Due to the opening of the valve, the normal fuel ignition characteristic of the generating plant changes.

Formulation of Economic Dispatch along with Valve Point
When using a multi-valve steam turbine system, the fuel cost function of th ating units may vary [4].The cost of generation increases with the use of mu steam turbines.The generator heat rate curve is the result of the multi-valve st bine valve opening operation.The impact of valve-point loading is shown in [22].In Figure 1, A, B, C, D, and E are the points where the valve can be opened t the speed of the generator if the load suddenly increases.Due to the opening of t the normal fuel ignition characteristic of the generating plant changes.
The importance of this impact lies in the fact that the huge steam-produ curve's function is, in reality, nonlinear rather than continuous.In reality put-output curves for generating units are not a smooth cost function when con a multi-valve steam operating system.The representation of the cost functio generating units is changed by the inclusion of the valve-point loading.Hence nomic load dispatch cost function is modified when considering the valve poi [22]; it is mathematically formulated as given in Equation (3).

Formulation of the ELD with a Ramp Rate Limit
When generating units are operated online, the actual operating ranges de the ramp rate limits  The importance of this impact lies in the fact that the huge steam-producing cost curve's function is, in reality, nonlinear rather than continuous.In reality, the input-output curves for generating units are not a smooth cost function when considering a multi-valve steam operating system.The representation of the cost functions of the generating units is changed by the inclusion of the valve-point loading.Hence, the economic load dispatch cost function is modified when considering the valve point effects [22]; it is mathematically formulated as given in Equation (3).

Formulation of the ELD with a Ramp Rate Limit
When generating units are operated online, the actual operating ranges depend on the ramp rate limits [14].The operation of the generating units at the different intervals is shown in Figure 2. The steady-state condition of generation is shown in Figure 2a; when generation increases, it is represented in Figure 2b; Figure 2c shown when generation decreases.When power generation increases, as shown in Figure 2b, it is represented as and when power generation decreases, as shown in Figure 2c, it is given as Now the economic load dispatch, with an up rate and down rate of power, is defined as Where URi and DRi are the upper and down ramp rate limits, respectively.

Formulation of the EELD Multi-Objective Problem
When generating power in a thermal power plant, coal is used as fuel; due to the burning of coal, toxic gases like COx, NOx, and SOx are emitted.Such toxic gases not only pollute the environment but also affect human health.So, reducing the emission of gases is a great achievement of the ELD solution.So, when considering the minimization of generation costs and environmental emissions, it is formulated as a multi-objective EELD (economic emission load dispatch) problem [2].
The first objective is the minimization of the cost function already discussed, given in Equations ( 1) and (2); the minimization of the environmental emission considered as the second objective is given as To obtain the economic emission dispatch, both objectives are combined using the penalty factor given in Equation (11).The ratio of the power produced by the plant to the actual power requirement of the load being satisfied after line loss is known as the penalty factor (hi) where FC (P .) = a P .+ b P .+ c and ( 12) When power generation increases, as shown in Figure 2b, it is represented as and when power generation decreases, as shown in Figure 2c, it is given as Now the economic load dispatch, with an up rate and down rate of power, is defined as P i,min where UR i and DR i are the upper and down ramp rate limits, respectively.

Formulation of the EELD Multi-Objective Problem
When generating power in a thermal power plant, coal is used as fuel; due to the burning of coal, toxic gases like COx, NOx, and SOx are emitted.Such toxic gases not only pollute the environment but also affect human health.So, reducing the emission of gases is a great achievement of the ELD solution.So, when considering the minimization of generation costs and environmental emissions, it is formulated as a multi-objective EELD (economic emission load dispatch) problem [2].
The first objective is the minimization of the cost function already discussed, given in Equations ( 1) and (2); the minimization of the environmental emission considered as the second objective is given as To obtain the economic emission dispatch, both objectives are combined using the penalty factor given in Equation (11).The ratio of the power produced by the plant to the actual power requirement of the load being satisfied after line loss is known as the penalty factor (h i ) Fc i (P i.max where FC i (P i.max ) = a i P where E T is the total economic emission cost.

Constraints
When optimized, the ELD problem is required to follow some restrictions or limitations, called constraints.During the formulation of the ELD mathematical model, keep in mind that the constraints should be inserted and then optimized.This work included the following constraints.

Power Balance
The power balancing restrictions make sure that the total amount of power generated by all the different kinds of generating unit is sufficient to satisfy the power load in each time period.

Power Limits
The generator's output power must not exceed its nominal value and must not fall below the level required for the boiler to operate steadily.As a result, the production is capped at a certain minimum and maximum range.The following equation may be used to represent each productivity unit with an expected production in the circuit.P min ≤ P i ≤ P max (16)

Heuristic Optimization Technique
The ELD problem is nonlinear in nature because of the presence of different constraints.For the optimization of such nonlinear ELD problems, classical methods were found to be inefficient, and hence, new heuristic optimization techniques were required.The PSO was first introduced in 1995 [1].After that, many new variants of PSO were investigated and implemented in optimization problems.The PSO is used for the optimization of various nonlinear problems; however, it has some disadvantages, such as difficulties in premature convergence, lagging performance for large test data, diversity in the local solutions, and the appropriate tuning of its parameters.
This work is considered a modified particle swarm optimization technique.A modified PSO is an effective variant of a PSO and can easily optimize single and multi-objective ELD problems.In the MPSO, we are using an attraction factor as shown in Equation ( 25) that has the following features:

•
It controls the movements of the particles in such a way that they cannot move away from the search area; • It continuously changes the particles' positions so the particles cannot stop during the iteration process until the final results are obtained; therefore, there is no need to update the velocity of the particles after every iteration; • It also attracts the particles so they cannot move to the local solution.

Initialization of the Swarm
To optimize the proposed technique, the swarm is required to be initialized randomly within the effective real power minimum and maximum limits of the power generated (P min and P max ) using the generating units, as given in Equation ( 17), and the velocity of the initialized particles, as given in Equation (18).S k i = P min + rand(P max − P min ) where P initial and V initial are the initial values of the swarm's position and its velocities are randomly generated using a random positive number (rand) between zero and one using MATLAB software when the algorithm is executed.Now, initialize the limits of the velocity of the particles using Equations ( 19) and ( 20).

Updating the Velocity of the Particle
The velocity of the particle is updated using the personal best and the global best value of the swarm is given as Inertia weight helps to accelerate the particles in the search space and is given as During the considered optimization, there were 100 iterations: W max and W min were considered 0.9 and 0.4, respectively, and c 1 and c 2 were 2.05.

Updating the Particles' Positions
This PSO used one attraction factor (S d ) that can attract the particles in the search area and find the global solution to the problem.N is the number of generating units.It is represented as follows.
where rand is a random variable generated using MATLAB between 0 and 1, P best is the best value of the particle, and g i,best is the global value.The positions of the particles are updated by using the best global position of the particles and an attraction factor, as follows.

Initialization of the Swarm
To optimize the proposed technique, the swarm is required to be initialized randomly within the effective real power minimum and maximum limits of the power generated (Pmin and Pmax) using the generating units, as given in Equation ( 17), and the velocity of the initialized particles, as given in Equation ( 18).

S = P + rand(P − P )
where Pinitial and Vinitial are the initial values of the swarm's position and its velocities are randomly generated using a random positive number (rand) between zero and one using MATLAB software when the algorithm is executed.Now, initialize the limits of the velocity of the particles using Equations ( 19) and (20).

Updating the Velocity of the Particle
The velocity of the particle is updated using the personal best and the global best value of the swarm is given as Inertia weight helps to accelerate the particles in the search space and is given as During the considered optimization, there were 100 iterations: Wmax and Wmin were considered 0.9 and 0.4, respectively, and c1 and c2 were 2.05.

Updating the Particles' Positions
This PSO used one attraction factor (Sd) that can attract the particles in the search area and find the global solution to the problem.N is the number of generating units.It is represented as follows.
where rand is a random variable generated using MATLAB between 0 and 1, Pbest is the best value of the particle, and gi,best is the global value.The positions of the particles are updated by using the best global position of the particles and an attraction factor, as follows.

Initialization of the Swarm
To optimize the proposed technique, the swarm is required to be initialized randomly within the effective real power minimum and maximum limits of the power generated (Pmin and Pmax) using the generating units, as given in Equation ( 17), and the velocity of the initialized particles, as given in Equation (18).S = P + rand(P − P ) (17) where Pinitial and Vinitial are the initial values of the swarm's position and its velocities are randomly generated using a random positive number (rand) between zero and one using MATLAB software when the algorithm is executed.Now, initialize the limits of the velocity of the particles using Equations ( 19) and (20).V = (P − P )/10 (19)

Updating the Velocity of the Particle
The velocity of the particle is updated using the personal best and the global best value of the swarm is given as Inertia weight helps to accelerate the particles in the search space and is given as During the considered optimization, there were 100 iterations: Wmax and Wmin were considered 0.9 and 0.4, respectively, and c1 and c2 were 2.05.

Updating the Particles' Positions
This PSO used one attraction factor (Sd) that can attract the particles in the search area and find the global solution to the problem.N is the number of generating units.It is represented as follows.
where rand is a random variable generated using MATLAB between 0 and 1, Pbest is the best value of the particle, and gi,best is the global value.The positions of the particles are updated by using the best global position of the particles and an attraction factor, as follows.
where R1 and R2 are random variables within [0, 1] and R3 within

Algorithm of the Proposed MPSO
 Consider the number of particles;  Define the number of iterations; Define the number of iterations; Electronics 2023, 12, x FOR PEER REVIEW 8 of 16

Initialization of the Swarm
To optimize the proposed technique, the swarm is required to be initialized randomly within the effective real power minimum and maximum limits of the power generated (Pmin and Pmax) using the generating units, as given in Equation ( 17), and the velocity of the initialized particles, as given in Equation (18).S = P + rand(P − P ) where Pinitial and Vinitial are the initial values of the swarm's position and its velocities are randomly generated using a random positive number (rand) between zero and one using MATLAB software when the algorithm is executed.Now, initialize the limits of the velocity of the particles using Equations ( 19) and (20).

Updating the Velocity of the Particle
The velocity of the particle is updated using the personal best and the global best value of the swarm is given as Inertia weight helps to accelerate the particles in the search space and is given as During the considered optimization, there were 100 iterations: Wmax and Wmin were considered 0.9 and 0.4, respectively, and c1 and c2 were 2.05.

Updating the Particles' Positions
This PSO used one attraction factor (Sd) that can attract the particles in the search area and find the global solution to the problem.N is the number of generating units.It is represented as follows.
where rand is a random variable generated using MATLAB between 0 and 1, Pbest is the best value of the particle, and gi,best is the global value.The positions of the particles are updated by using the best global position of the particles and an attraction factor, as follows.
where R1 and R2 are random variables within [0, 1] and R3 within

Algorithm of the Proposed MPSO
 Consider the number of particles;  Define the number of iterations; Initialize the swarm using the minimum and maximum values of the generated power, as given in Equations ( 17) and (18); Electronics 2023, 12, x FOR PEER REVIEW 8 of 16

Initialization of the Swarm
To optimize the proposed technique, the swarm is required to be initialized randomly within the effective real power minimum and maximum limits of the power generated (Pmin and Pmax) using the generating units, as given in Equation ( 17), and the velocity of the initialized particles, as given in Equation (18).S = P + rand(P − P ) where Pinitial and Vinitial are the initial values of the swarm's position and its velocities are randomly generated using a random positive number (rand) between zero and one using MATLAB software when the algorithm is executed.Now, initialize the limits of the velocity of the particles using Equations ( 19) and (20).

Updating the Velocity of the Particle
The velocity of the particle is updated using the personal best and the global best value of the swarm is given as Inertia weight helps to accelerate the particles in the search space and is given as During the considered optimization, there were 100 iterations: Wmax and Wmin were considered 0.9 and 0.4, respectively, and c1 and c2 were 2.05.

Updating the Particles' Positions
This PSO used one attraction factor (Sd) that can attract the particles in the search area and find the global solution to the problem.N is the number of generating units.It is represented as follows.
S = P * rand + g , * (1 − rand) where rand is a random variable generated using MATLAB between 0 and 1, Pbest is the best value of the particle, and gi,best is the global value.The positions of the particles are updated by using the best global position of the particles and an attraction factor, as follows.

Initialization of the Swarm
To optimize the proposed technique, the swarm is required to be initialized randomly within the effective real power minimum and maximum limits of the power generated (Pmin and Pmax) using the generating units, as given in Equation ( 17), and the velocity of the initialized particles, as given in Equation (18).S = P + rand(P − P ) where Pinitial and Vinitial are the initial values of the swarm's position and its velocities are randomly generated using a random positive number (rand) between zero and one using MATLAB software when the algorithm is executed.Now, initialize the limits of the velocity of the particles using Equations ( 19) and (20).

Updating the Velocity of the Particle
The velocity of the particle is updated using the personal best and the global best value of the swarm is given as Inertia weight helps to accelerate the particles in the search space and is given as During the considered optimization, there were 100 iterations: Wmax and Wmin were considered 0.9 and 0.4, respectively, and c1 and c2 were 2.05.

Updating the Particles' Positions
This PSO used one attraction factor (Sd) that can attract the particles in the search area and find the global solution to the problem.N is the number of generating units.It is represented as follows.
S = P * rand + g , * (1 − rand) where rand is a random variable generated using MATLAB between 0 and 1, Pbest is the best value of the particle, and gi,best is the global value.The positions of the particles are updated by using the best global position of the particles and an attraction factor, as follows.The optimum results obtained with the proposed MPSO are listed in Table 4.The results obtained with the proposed algorithm are compared for validation with recently developed optimization algorithms like the PPSO [14], Fuzzy and APSO [23], GPSO-w [22], GA [24], G-SCNHGWO [28], and ESCSDO [30].Results shown in the table clearly show the effectiveness of the modified PSO; it performs better and has a minimum generation cost.The minimum fuel cost given by the MPSO is USD 32,465.69/hand the computation time is 0.824 s.

Case 3
This case study is considered multi-objective test data.In this case, we evaluated the economic load dispatch and environmental emissions.Then, both objectives are combined to find the overall economic emissions dispatch.The test data for this case study are shown in Table 5.The case data are tested for a load demand of 900 MW for six generating units along with the emission and fuel cost coefficients.Now, the six generating units' data are analyzed for the load demand of 900 MW.The case data are tested with environmental emissions.The results obtained with the proposed modified PSO technique are shown in Table 6.The MSPO results were compared, for validation, with the BAT [17], ABC [16], RGA [17], and ACB [18] and the proposed algorithm was found to gives the best results compared to the other techniques, as shown in Table 6.The MPSO gives a minimum generation cost of USD 47,889.45/h, a minimum emission discharge of 669.3217T/h, and a total economic emission dispatch cost of 81,489.12for a load demand of 900 MW.Before combining the fuel cost and emission required values of the penalty factor shown in Equation (11) using the Equations ( 12) and ( 13), the penalty factor (hi) is calculated for six generating units, as given in Table 7.

Conclusions
Policies are required for the generation of power in such a way that it fulfils load demands, generates the power at a low cost, and emits fewer toxic gases into the environment.ELD is one of the best optimization approaches that can help generate power with a minimum fuel cost.This work proposes a new PSO for the optimization of economic load dispatch as a single objective for thirteen generating units with valve loading effects and for fifteen generating units with ramp rate limits.The performance of the proposed PSO is important in terms of the minimum generation cost and the computation time taken.According to the results shown for both single objective cases, the performance of the proposed new PSO is better than the other optimization techniques; it also has a

Figure 1 .
Figure 1.Effect of the valves' loading on generation cost function.
[14].The operation of the generating units at the different in shown in Figure 2. The steady-state condition of generation is shown in Figure generation increases, it is represented in Figure 2b; Figure 2c shown when ge decreases.

Figure 1 .
Figure 1.Effect of the valves' loading on generation cost function.

Figure 2 .
Figure 2. Power generation limits due to ramp rate limits.

Figure 2 .
Figure 2. Power generation limits due to ramp rate limits.

Figure 3 .
Figure 3. Pseudo code for the NPSO.Figure 3. Pseudo code for the NPSO.

Figure 3 .
Figure 3. Pseudo code for the NPSO.Figure 3. Pseudo code for the NPSO.
Initialize the velocity of the particles, as shown in Equations (19) and (20);

Table 3 .
Fuel cost coefficients, ramp rate limits and the generation limits of 15 generating unit systems for a load demand of 2650 MW.

Table 4 .
Test results of case 2 for 15 generating unit systems.

Table 5 .
Test data for case 3 for a load demand of 900 MW.

Table 6 .
Results of the six generating units for the economic emission dispatch for a demand of 900 MW.

Table 7 .
Penalty factor h values.