Development of Frequency Weighted Model Reduction Algorithm with Error Bound: Application to Doubly Fed Induction Generator Based Wind Turbines for Power System

: The state-space representations grant a convenient, compact, and elegant way to examine the induction and synchronous generator-based wind turbines, with facts readily available for stability, controllability, and observability analysis. The state-space models are used to look into the functionality of different wind turbine technologies to fulﬁll grid code requirements. This paper deals with the model order reduction of the Variable-Speed Wind Turbines model with the aid of improved stability preserving a balanced realization algorithm based on frequency weighting. The algorithm, which is in view of balanced realization based on frequency weighting, can be utilized for reducing the order of the system. Balanced realization based model design uses a full frequency spectrum to perform the model reduction. However, it is not possible practically to use the full frequency spectrum. The Variable-Speed Wind Turbines model utilized in this paper is stable and includes various input-output states. This brings a complicated state of affairs for analysis, control, and design of the full-scale system. The proposed work produces steady and precise outcomes such as in contrast to conventional reduction methods which shows the efﬁcacy of the proposed algorithm.


Introduction
Among the available sources of energy such as oil, gas coal, and nuclear, the wind is recognized as an infinite source of clean energy. Wind energy has reduced the sole reliance on fossil fuel means. During the last two decades, the rapid growth of wind power plants is observed all over the world [1][2][3][4]. According to the global wind energy council (GWEC), by the completion of the year 2019, the total cumulative installed capacity of wind energy was 60.4 GW [5] around the globe, which is a 19 % globally increase from installation in 2018 [6]. Moreover, the second-best year for wind energy historically. The total cumulative installed capacity for wind energy is 651 GW around the globe, a 10 % increase compared to 2018 (see GWEC for more detail). Continual change in speed, density, and temperature round the clock is the main issue with wind energy. Therefore the integration of the wind masts with the grid needs to be governed by certain policies, known as grid codes, to avoid the undesirable impact on the grid power [1,[7][8][9]. The developed grid codes are different for the different countries depending upon their environmental conditions and operational 1.
Frequency operating range.
Communications and external control.
• Power quality; • Wind farm modeling and verification.
Traditionally three types of power generation machines are used in wind turbines to convert wind energy into electrical energy. These are Squirrel (single/double) Cage Induction Generator (SCIG), Doubly Fed (wound rotor) Induction Generator (DFIG), and Permanent Magnet Synchronous Generator (PMSG). Among the said three generators, the DFIG has shown a good performance during LVRT and remained connected to the power grid. DFIGs are frequently being used in the wind energy conversion systems due to the distinguished features of handling of variable speed, easy to control, higher energy efficiency, and improved power quality [13]. Nevertheless, Squirrel Cage Induction Machines (SCIM) equipped with reactive power compensators are also used in wind energy conversion systems because of their simplicity [14], low cost, and robustness. The SCIM with an additional cage, also called double squirrel cage, has higher-order complexity as compared to the single cage machines.
Hence, the computer simulation-based study, investigations and research demands a realization of a perfect depiction of Double Cage Induction that can address their significant issues, especially, concerning the integration of wind energy conversion systems to the grid. Thus, the use of wind energy and its integration with the grid is now an essential subject to learn about [15].
Comprehensive studies to discover the interplay between wind farms and the power system are necessary. Different researches are carried out whenever a wind farm is designed, in a comparable way to other new technology amenities [13,15]. In general, the influences of wind technology are assessed in the course of model planning concerning voltage profile, electricity flow, short-circuit currents, reactive power capability, Low Voltage Ride Through (LVRT), and transient stability [16][17][18][19]. A detailed representation of each unique unit and connections between the units and the system is every so often considered. Alternatively, the wind farm can also be modeled as a lumped equivalent model as viewed from the system [17]. The other related work on transient stability improvement and spinning reserve dispatch in wind-thermal power system is also given in [20,21].
Nevertheless, the response of the system will rely on the type of tools being used. Several different dynamic models have been derived and are entirely documented in the literature [22][23][24][25][26][27]. However, some of them use significant simplification to make use of control system design techniques [28].
To provide a high-quality tool for energy system studies, the reduced-order modeling of state-space representations of synchronous and induction machines for wind turbine applications is given [29]. The reduced models are presented as state-space representations, allowing a convenient, compact and elegant way to check the induction and synchronous generator-based wind turbines [30][31][32], with facts effectively accessible for stability, controllability, and observability analyses [33,34]. Natural and steady-state responses of reduced-order model (ROM) effectively compared with the original system which shows the effectiveness of ROM by using balanced based realization [35].
Model order reduction (MOR) is a process of reducing the higher-order system to lower-order systems. A Higher-order system involves lots of complexity for the analysis, design, and simulations. In control system theory it is a very useful way to reduce the higher-order model and perform tasks [36][37][38][39][40][41][42][43][44][45][46].
One of the most extensively used technique of MOR is balance truncation (BT) [47]. It is used to acquire the ROM. It grants the stable ROM and also grants error bound formula. It keeps the properties like input/output behavior and passivity etc while the use of the full frequency spectrum. The BT approach [47] yields ROM which not only preserves stability but additionally gives easily computable error bounds [47]. Ideally, the BT [47] approach tends to approximate the original model and offers ROM with less approximation error for the entire frequency range. However, there is a certain scenario where approximation error is required to be small for certain frequency weights alternatively than the entire frequency range. This motivates to use of frequency weights in a balanced realization algorithm [48]. Enns [48] has extended the BT [47] approach by way of introducing frequency weights in the given model realization. This technique may additionally include input, output, and two-sided weighting. However, this method sometime yields unstable reduced order systems when both sided weighting is present [49]. To overcome this main drawback, Wang et al's [50] and Imran and Ghafoor (IG) [51] presented stability preserving methods. Gawronski and Juang [52] and Wang and Zilouchian (WZ) [53] proposed a frequency confined balanced MOR technique for the continuous-time and the discrete-time systems respectively, where weights are not explicitly predefined, but the approximation is considered in certain frequency intervals. In this technique, Gramians were described for preferred frequency intervals. However, it can also yield unstable ROM for a stable original system. Moreover, there is no error bound exists. The instability issues appeared in [52,53] were carried out by [54][55][56][57][58][59].
To overcome the instability issue of Enns [48], the proposed scheme effectively performs a model reduction of the DFIG based Variable-Speed Wind Turbines model [24] and produces stable ROM by ensuring the positive/semi-positive definiteness of some input and output related matrices respectively, simulation results that are comparable with [48] and existing stability preserving techniques ( [50,51]). Moreover, the proposed technique provides computable a priori error bounds formula for frequency weighting and limited frequency-intervals respectively. Simulation results are given to show the usefulness of the proposed technique when compared with other methods.
The main contributions of this paper are as follows

Grid Connection Configuration of Induction Machines
This section discusses the grid connection configuration of Double Fed Induction Generator (DFIG) and Squirrel Cage Induction Generator (SCIG).

Double Fed Induction Generator (DFIG)
DFIGs are integrated with the grid in a pattern as shown in Figure 1. A low-speed shaft of the wind turbine is connected to the rotor of the DFIG through a gearbox. The gearbox enhances the speed to a value required by the generator for power generation. DFIG uses a wound type rotor which is connected to the grid through, AC-DC-AC, two back-to-back voltage source converters of partial rating (usually 30 percent of generator power). The converter connected to the rotor end is called Rotor Side Converter (RSC) while the converter connected to the grid end is known as Grid Side Converter (GSC). Both converters are separated by a DC-link capacitor as energy storage. These converters address the problem of the variable wind speed, during both sub-synchronous and super-synchronous modes, and produce fixed output frequency as required at the grid [16]. The stator is connected to the grid through a step-up transformer [13]. A control system is incorporated to control input power at the shaft of the wind turbine, reactive power, and the voltage at the grid terminal. This system generates separate voltage commands V r and V gc for RSC and GSC respectively. The RSC provides the control of both active and reactive power while GSC manages the voltage at DC-link capacitor in between RSC and GSC for its operation at unity power factor.

Squirrel Cage Induction Generator (SCIG)
Direct grid-connected wind energy turbines use SCIGs (single or double cage) and operate at a fixed speed [14]. Similar to the other wind energy conversion systems, the turbine is linked to SCIG through a gearbox to achieve a nominal speed for power generation while the generator is directly coupled to the grid as displayed in Figure 2. As speed fluctuations are due to the variations in the rotor slip, the rotor speed changes are minor and the wind turbine is normally used to drive at a fixed speed. At the instant of voltage variations at the grid, SCIG acts as an induction motor and absorbs reactive power. Therefore, the power system is augmented with a reactive power compensator, usually a capacitor bank, to improve the power factor. The pitch angle control is incorporated for the optimum value of the wind power to address the rotor speed instability of the generator as the wind speed varies. The Squirrel cage systems are also used in variable-speed wind energy systems.

Mathematical Model for DFIG and SCIG
The variables indicated by prime are referred to the stator. The stator and rotor quantities are in arbitrary 'dq' (d-axis and q-axis) reference frame. Figure 1 depicts circuit diagram of DFIG and SCIG.

Mathematical Model for Double Cage Induction Machines (DCIM)
Mathematical relations for the conversion of wind energy into electrical energy is specified in [17][18][19]. Electrical variables and parameters designated by primes, in the mathematical model, are referred to as the stator. Parameters for the Double Cage Induction Machine (DCIM) are shown in Table 1.

Electrical System
The electrical system comprises a stator and a double squirrel cage rotor. The machine is represented by a sixth-order model while the mechanical part is a second-order system [22]. The electrical circuit is represented in Figure 3. Stator voltages on quadrature (q) axis and direct (d) axis can be obtained as following Rotor equations for the first cage quadrature (q) axis and direct (d) axis can be obtained as following Rotor equations for the second cage quadrature (q) axis and direct (d) axis can be obtained as following Electrical torque equations double squirrel cage induction machines be T e = 1.5p(ωϕ d s i q s − ωϕ q s i d s ) Figure 3. Electrical circuit of a double cage induction machines in the arbitrary two-axis, dqreference frame.

Remark 1.
During the conversion process of electromechanical energy, some of the energy constitutes as heat due to the resistive loss in the conductors of rotor and stator as well as core losses (hysteresis and eddy current) and dielectric loss. Technically, core losses are taken into account by placing a small resistor in parallel with the magnetizing inductance Lm in the per phase equivalent circuit of the model. However, these losses are ignored in the circuit diagram, as shown in Figure 3, due to the following reasons [22]: 1.
Selected core is made up of ferromagnetic material organized in lamination to minimize the core losses 2.
The core losses are negligible when compared to the stator and rotor copper winding losses

Mechanical System
A second-order mechanical system of double cage induction machine is represented as [25]

Mathematical Model for DFIG and SCIG Systems
For a computer simulations based study it is necessary to understand mathematical models of the machines for their electrical and mechanical systems [19]. Table 1 displays all the mathematical symbols and their meaning in the context of the mathematical models discussed.

Wind Energy Conversion
The kinetic energy of the wind is converted through the blades of turbine into mechanical power P m and computed as follows [19] where T m is the torque exerted on the mechanical shaft and ω m (rad/sec) is rotational speed of the turbine. In steady state operation at a fixed speed for a lossless generator, mechanical torque is equal to electrical torque i.e., T m = T e = P m /ω m [17,18]. The power produced by the wind is given as follows where ζ p = Performance coefficient of the turbine λ = Tip speed ratio of the rotor blade tip speed to wind speed = ω m /ν β = Blade pitch angle (deg) ρ = Density of Air (kg/m) r = Radius of the turbine blades (m) ν = Wind speed (m/s)

Electrical Systems for Wound Rotor (DFIG) and Squirrel Cage (SCIG) Machine
This section discusses mathematical model of electrical system. The variables indicated by prime are referred to the stator. The stator and rotor quantities are in arbitrary 'dq' (d-axis and q-axis) reference frame. Figure 4 depicts circuit diagram of DFIG and SCIG. The stator voltages are given as follow The stator voltages are given as folloẃ The active and reactive stator and rotor powers are related as follow [19]

Balancing Related Model Order Reduction Schemes
Consider a linear time invariant continuous time system where A ∈ n×n , B ∈ n×p , C ∈ q×n , D ∈ q×p and {A, B, C, D} is nth order minimal realization with p inputs and q outputs. The problem of model reduction is to find which approximates the actual system (in the full frequency band), where A tr ∈ r×r , B tr ∈ r×p , C tr ∈ q×r , D tr ∈ q×p with r << n.

Balance Truncation Technique [47]
Let P cg and Q og be the full frequency controllability and observability Gramians respectively are the solution of following Lyapunov equations: Let T B be a contragredient matrix obtained as where σ j ≥ σ j+1 , j = 1, 2, 3, . . . , n − 1 and σ r > σ r+1 , and r is the order of ROM. By applying the transformation and then partitioning the original system, the ROM G tr (s) = C tr (sI − A tr ) −1 B tr + D tr is obtained as: Remark 2. For the minimal and stable realization the transformed balanced realization is minimal and the stability of ROM (A tr , B tr , C tr ) is also guaranteed and yields frequency response error bounds. However, the balance truncation technique [47] uses a full frequency spectrum to perform MOR which not practically feasible all the time.

Enns's Technique [48]
Consider a stable input weighting stable model of continuous time systems where Let a stable output weighting model of continuous time systems where The augmented systems are given by satisfy the following Lyapunov equations: A T so Q so + Q so A so + C T so C so = 0 (10) Remark 3. In controller reduction, scenario pole-zero cancellation may occur that leads to loss of controllability and observability of realizations input-augmented {A si , B si , C si , D si } and outputaugmented {A so , B so , C so , D so } respectively.
By considering the first and fourth block of Equations (9) and (10) respectively, we have: where By eigenvalues decomposition of X en and Y en we have following X en = U en S en 1 0 0 S en 2 U T en = U en S en U T en (15) B en = U en S en 1 where k and q are the number of positive eigenvalues of X en and Y en respectively. Let T en be a contragredient matrix obtained as where σ j ≥ σ j+1 , j = 1, 2, 3, . . . , n − 1 and σ r > σ r+1 and r is the order of ROM. By applying the transformation and then partitioning the original system, the ROM G tr (s) = C tr (sI − A tr ) −1 B tr + D tr is obtained as as similar as way as in (5)-(6).
Remark 4. P en and Q og are used to obtain the balancing (contragredient) transformation T en in (19) subject to only input weights are present. Likewise, P cg and Q en are used to obtain balancing (contragredient) transformation T en in (19) subject to only output weights are present, where P cg and Q og are un-weighted Gramians computed as in (3)-(4).

Remark 5.
Enns's technique [48] do not ensure the stability of the ROM since X en ≤ 0 and Y en ≤ 0, it may cause sometimes unstable ROMs for both sided weighting case. Furthermore, the transformed realization may not be minimal due to pole-zero cancellation, therefore, this method is not useful especially in the controller reduction scenario [49].

Wang and Sreeram's Technique [50]
Wang and Sreeram [50] provides the solution to the instability issue that appeared in [48]. This method [50] ensured the positive/semi-positive definiteness of input and output related matrices X en (13) and Y en (14) respectively by doing some variation. Let the improvised controllability and observability Gramians P WS and Q WS (calculated by solving the following Lyapunov equations) respectively.

Imran and Ghafoor's Technique [51]
Imran and Ghafoor [51] technique also provided the solution to instability issue appeared in [48], by ensuring positive/semi-positive definiteness of input and output related matrices X en and Y en respectively. However, this technique produces large approximation error due to large variation in X en and Y en . The controllability and observability Gramians P IG and Q IG (calculated by solving the following Lyapunov equations) respectively can be obtained as: the fictitious input and output related matrices B IG and C IG respectively are defined as: The terms U IG , S IG , V IG , and R IG are calculated by the orthogonal eigenvaluesdecomposition of symmetric matrices B IG B T IG = U IG S IG U T IG and C T IG C IG = V IG R IG V T IG , where S IG = S en − s n I = diag(s 1 ,s 2 , s 3 , · · · ,s n−1 , 0), R IG = R en − r n I = diag(r 1 ,r 2 , r 3 , · · · , r n−1 , 0),s 1 ≥s 2 ≥ · · · ≥s n−1 ≥ 0, andr 1 ≥r 2 ≥ · · · ≥r n−1 ≥ 0. Let T IG be contragredient matrix obtained as whereσ j ≥σ j+1 , j = 1, 2, 3, . . . , n − 1 andσ r >σ r+1 and r is the order of the ROM. By applying the transformation and then partitioning the original system, the ROM G tr (s) =C tr (sI −Ā tr ) −1B tr +D tr is obtained as similar way as in (5-6).

Remark 9.
Since X en ≤ B IG B T IG ≥ 0, Y en ≤ C T IG C IG ≥ 0, P IG > 0 and Q IG > 0. Therefore, the realization (A, B IG , C IG , D) is minimal. Moreover, ROMs are stable.

Remark 10.
The stability of ROMs in the presence of both input and output weighting is guaranteed and the following error bound holds [51]

Main Results
In [48], input related and output related matrices X en and Y en respectively are not ensured to be positive/semi-positive definite due to negative eigenvalues obtained by using eigenvalue decomposition of X en and Y en . Whereas, the proposed technique ensured the positive/semi-positive definite of the matrices X en and Y en by introducing small variation in input related and output related matrices X en and Y en respectively. This is achieved by subtracting all the eigenvalues with the negative sum of eigenvalues which results in stable ROM. The proposed technique also provides computable a priori error bounds formula.
Let controllability P SB and observability Q SB Gramians satisfying the following Lyapunov equations: where X SB = B SB B T SB and Y SB = C T SB C SB By eigenvalues decomposition of X SB = U SB S SB U T SB and Y SB = VR SB V T . The new imaginary input and output related matrices B SB and C SB respectively are given as following where S SB 2 = (S en 2 −ŝI) 1/2 for s n < 0 S 1/2 en 2 for s n ≥ 0 The terms U SB , S SB , V SB , and R SB are calculated by the orthogonal eigenvalues-decomposition of input and output related matrices B SB B T SB = U SB S SB U T SB and C T SB C SB = V SB R SB V T SB respectively, where S SB = S en −ŝ n I = diag(ŝ 1 ,ŝ 2 ,ŝ 3 , · · · ,ŝ n−1 ,ŝ n ), R SB = R en −r n I = diag(r 1 ,r 2 ,r 3 , · · · ,r n−1 ,r n ),ŝ 1 ≥ŝ 2 ≥ · · · ≥ŝ n ≥ 0, and r 1 ≥r 2 ≥ · · · ≥r n ≥ 0. Note that, the matrices B SB and C SB are constructed by ensuring similar effect on all eigenvalues of symmetric matrices X en and Y en . This will guarantee the preserving of eigenvalue structure of matrices X en and Y en and therefore, better results are obtained in proposed technique as compared to Wang and Sreeram [50] and Imran and Ghafoor [51] technique. A contragradient transformation matrix T SB is obtained as whereσ j ≥σ j+1 andσ r ≥σ r+1 and r is the order of the ROM. By applying the transformation and then partitioning the original system, the ROM G tr (s) =Ĉ tr (sI −Â tr ) −1B tr +D tr is obtained as: Lemma 1. X en ≤ B SB B T SB ≥ 0, Y en ≤ C T SB C SB ≥ 0, likewise, P en < P SB > 0 and Q en < Q SB > 0. Therefore, minimal and stable realization (A, B SB , C SB ) is obtained which guaranteed the stability of the ROM.
Proof of Lemma 1. We will first show that the realization {A, B SB , C SB } is minimal. Since P SB and Q SB are solution of Lyapunov Equations (26) and (27) respectively, so Since P en is positive definite, P SB must be positive definite. Similarly, we can say Q SB is positive definite. Since P SB and Q SB are positive definite and A is stable, it follows immediately that the pair (A, B SB ) is controllable and (A, C SB ) is observable or in other words the realization {A, B SB , C SB } is minimal.
If {Â tr , B SB 1 , C SB 1 ,D tr } is approximated model obtained by partitioning a balanced realization {Â, B SB , C SB , D}, we have from [48] For the case G(s)V i (s), the error bound expression becomes  [50,54]) are satisfied.
Remark 11. When input and output related matrices X en ≥ 0 and Y en ≥ 0 respectively, then P en = P SB and Q en = Q SB . Otherwise P en < P SB and Q en < Q SB . Moreover, frequency weighted Hankel singular values satisfies: (λ j [P en Q en ]) 1/2 ≤ (λ j [P SB Q SB ]) 1/2 Remark 12. When input V i (s) and output W o (s) weights are co-inner and inner respectively [51], then P cg = P en = P SB and Q og = Q en = Q SB .

Remark 13.
When input and output related matrices X en ≥ 0 and Y en ≥ 0 respectively, then approximated models obtained using [48] and proposed technique are the equivalent.

Remark 15.
Since the input matrix B SB and output matrix C SB ensure positive (semi-positive) definiteness of input and output related matrices, consequently positive definiteness of P SB and Q SB in a way leads to transformation matrix T SB which subsequently results in stability preserving MOR technique. Moreover, L SB and K SB form bases for the derivation of the error bound for the proposed technique.

Computational Aspects
Let S si and R so be the Cholesky factors of the augmented system Gramians matrices P si and Q so of Equations (9) and (10) respectively, By making use of the Cholesky factors S so and R so calculated above, the Cholesky factors corresponding to Gramians in frequency weighted model reduction techniques like [48] and proposed technique can be received as follows: 1.

Enns's Technique [48]:
The Cholesky factors S en = [S 11 S 12 ] and R en = R 12 R 22 satisfy [54] P en = S en S T en = S 11 S T 11 + S 12 S T 12 = S 11 S 12 S T

22
and

Wang and Sreeram's Technique [50]:
The Cholesky factors S WS and R WS satisfy P WS = S WS S T WS and Q WS = R T WS R WS , where P WS (20) and Q WS (21) [51]. Compute from the realization {A, B WS , C WS } using method of Hammarling [60].

3.
Imran and Ghafoor's Technique [51]: The Cholesky factors S IG and R IG satisfy P IG = S IG S T IG and Q IG = R T IG R IG , where P IG (24) and Q IG (25).

4.
Proposed Technique: The Cholesky factors S SB and R SB satisfy P SB = S SB S T SB and Q SB = R T SB R SB , where P SB (26) and Q SB (27). Next we establish a relationship between Cholesky factors Gramian matrices of Enns and proposed technique. Equations (26) and (27) can be expressed as: A(P en + P (ext) ) + (P en + P (ext) )A T + (X en + X (ext) ) = 0, for s n < 0 AP en + P en A T + X en = 0, for s n ≥ 0 A T (Q en + Q (ext) ) + (Q en + Q (ext) )A + (Y en + Y (ext) ) = 0, for r n < 0 By using the Hammarling's technique [60] to calculate the Cholesky factors of the Gramians P (ext) and Q (ext) from the realization {A, B (ext) , C (ext) , D}, we can write P (ext) = S (ext) S T (ext) and Q (ext) = R T (ext) R (ext) . Therefore, P SB (26) and Q SB (27) can be expressed as: Note that, the Cholesky factors for the Enns and the proposed technique respectively, are computed directly from the augmented system realization using Hammarling technique without calculating the augmented system realization Gramian matrices P si and Q so . However, Cholesky factorization for Wang et al's technique is computed directly from corresponding frequency weighted realization using the Hammarling technique without calculating associated frequency weighted Gramians [61].

Analysis and Discussion
It can be seen that in the given input/output frequency weights in the proposed technique compare well with Enns's technique [48]. From Table 3 it can be seen that 3rd, 4th and 5th order ROMs of example 1 and 3rd, 4th and 5th order ROMs of Example 2 obtained by using Enns's technique [48] yields unstable in the given both-sided frequency weights. However, the proposed technique yield stable ROMs in given frequency weights. Figures 5 and 6) of Example 1 and Figures 7 and 8 of Example 2 respectively provide a comparison for the frequency response Bode plot (magnitude, phase) in the given frequency weights for (1st, 2nd) and (1st, 2nd) order ROMs respectively, which indicates that the proposed technique provide better approximation results when compared with existing frequency weighted stability preserving MOR approaches ( [50,51]). Moreover, the proposed technique also carry error bounds formula with low frequency-response approximation error along-with error bound value. From Table 2 it can be observed that frequency response approximation error attained by the proposed technique is mostly low as compared with other MOR techniques.

Conclusions
In this work, the frequency weighted MOR technique is presented which approximates the Variable-Speed Wind Turbines model in the presence of frequency weights. The proposed method is used to derived a stable ROMs and frequency response error comparison at different order is provided which is well comparable with the existing stability preserving MOR techniques (frequency weighted MOR). The frequency response error and error bound comparison of the proposed technique with existing methods is provided which shows that the proposed method provides better results. Numerical analysis of DCIM has proven that ROMs obtained by using the proposed method are stable. However, analysis, control, and design of the original large scale model is not an easy task to do, whereas, it can be performed easily for approximated systems. This paper produces steady and precise outcomes by using the proposed method which shows the effectiveness of the proposed method.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: