System identification of photovoltaic system based on fractional-order model

System identification of photovoltaic system based on fractional-order model | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article System identification of photovoltaic system based on fractional-order model rabiaa ebead, Belal Abou-Zalam, Essam Nabil This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-1865340/v1 This work is licensed under a CC BY 4.0 License Status: Under Review Version 1 posted 5 You are reading this latest preprint version Abstract A fractional-order model of a photovoltaic (PV) system is proposed in this paper. The system identification method is used to construct an accurate dynamical model for a PV system, a real PV module and boost converter are used to collect experimentally input-output data for the identification process. A black box modeling is considered for system identification to obtain a transfer function, without needing to make mathematical analysis. the input-output data which is measured is fitted by the least-square curve fitting method, and the parameters of the curve fitting are optimized by the Levenberg-Marquardt algorithm. the proposed fractional-order model is employed with the MATLAB and Simulink software to obtain a fractional order PID controller, which is used to extract maximum power from the PV system. fractional-order model System identification PV system least-square curve fitting and PV system. Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Figure 9 I. Introduction The generation of electrical energy has been dependent on fossil fuels for many years, such as coal, natural gas, and oil 1 , which cause environmental issues, and fossil fuel reserves decrease. These problems have made electricity generation go towards renewable energy sources such as gravitational energy, geothermal energy, wind energy, tides, and solar energy 2,3 . renewable energy sources have been taking a significant position in the replacement or complement of conventional power generation, but solar energy source among all renewable energy sources is used widely because of advantages such as having no noise due to the absence of mechanical parts, low maintenance cost, long working life, the great amount of solar irradiation that covers the most of the earth’s surface 4 . The electric energy is generated from solar energy by a small unit called a photovoltaic (PV) cell, which produces a little amount of energy, several PV cells are connected in parallel and series to form a PV module to have more power, and PV panel is a group of PV modules which are joined together in parallel or series form. In a PV system, PV panels are connected in series and parallel to form PV arrays for getting high current and high voltage, to increase the output power of PV arrays 5 . The productivity of PV modules to convert solar energy to electric energy is dependent on the solar irradiance level, the surrounding temperature, and others 6,7 , so the variable environmental conditions cause variation in the voltage and the current of the PV module. therefore, the output power which is generated by the PV module changes also. Hence In a PV module, voltage and current have a nonlinear relation 8 , and in the electrical characteristic of the PV module, there is only one operating point, at which PV output power is maximum, this operating point is called Maximum Power Point (MPP), which makes PV modules are working efficiently. As this point fluctuates with variable environmental conditions, it should design a maximum power point tracking (MPPT)algorithm 9 , MPPT algorithms increase PV system efficiency by moving the PV voltage or current to make PV operate at the maximum power point (MPP). Another problem that has a great effect on the PV efficiency is Partial shading (PS), PS decreases the amount of sunlight that is allocated to the PV array resulting in lower output power from the PV array, under PS multiple local maximum power points (LMPP) emerge on PV curve 10 , but there is one MPP of LMPP is called global maximum power point. Therefore, to overcome the above issues, many MPPT algorithms are proposed to get maximum power from PV array under variable environmental conditions, in 11 the MPPT algorithms are classified into conventional techniques (CTs), soft computing techniques (SCTs), linear and nonlinear control techniques, all these algorithms are used in a simulation environment before real systems, PV model has many important roles in design optimization, fault diagnosis 12 , and prediction performance in the future of any PV system, therefore it should develop an accurate PV model to design MPPT. The mathematical models of real systems can be divided into static and dynamic models 13 , the dynamical model for the PV system is preferred due to the loads always are not constant, and the PWM which is used for switching converters or inverters, the static model for PV module is introduced in many pieces of literature such as in 14–18 , a single diode model is considered as a static PV model and is the most one which is used, the static model can be derived from parameter identification, parameter identification method tray to make approximation and solve analytically the equation of the PV as a polynomial of the first order [ 15 ], this means the static model is simplified and has low accuracy. Alternatively, a second-order dynamic model for PV module is proposed in 19 , the technique which is proposed increases the model accuracy, but it doesn’t consider the power electronics system, another dynamic model is introduced for a PV system based on system identification and considered buck converter, MPPT, and load in 20 , this technique used experimental input-output data and nonlinear Hammerstein-Wiener model to extract dynamic model. This paper proposes a dynamic fractional order model for a PV system, boost converter, MPPT, and load. Ii. Related Work The PV system in this paper is off-grid, and consists of a 60-watt PV module, boost converter, MPPT controller, and load. the circuit diagram of the PV system is shown in Fig. 1 The PV current measurement is used as input data for the system identification and output power of the boost converter as output data. The identification method to minimize the error of model output is based on Nonlinear least-squares estimation, Levenberg-Marquardt algorithm is used to solve nonlinear least-squares problems 21 , and for determining the parameters of the model. The proposed PV model is based on the fractional order model, fractional calculus is a generalization of ordinary differentiation and integration to arbitrary (non-integer) order, it is supposed that fractional order models will generally provide a more accurate description of the system dynamics than those based on conventional differential equations because it is given more degree of freedom for determining parameters of the model, the general definition for an integral differential operator is 22 : Where \(\alpha \in \mathbb{ }\mathbb{R}\) and \(\alpha\) can equal a complex number, there are many multiple definitions of the fractional operator, we consider the definition of Grünwald-Letnikov, which is used throughout this paper to drive numerical solutions for fractional-order differential equations. Fractional order systems modeling and controllers are implemented through the FOMCON toolbox for MATLAB/Simulink, the toolbox Fractional Order Modeling and Control “FOMCON” is developed in 23 , this toolbox will be used for system identification to extract the fractional order model of the PV system. Iii. Static Pv Model Characteristic A single diode model can describe the static behavior of a PV module and is proposed in many kinds of literature such as in 24,25 , the electric circuit for the single diode model is shown in Fig. 2 . From Fig. 2 , the mathematical equation which describes the static PV model can be driven as follow: Where \({I}_{ph}\) is the photocurrent, \({I}_{sat}\) is the reverse saturation current of the diode, q is the electron charge, A is the diode ideality factor, K is the Boltzmann constant, T is the PV module temperature, \({R}_{s}\) is the series resistor, \({R}_{sh}\) is the shunt resistor, and \({N}_{s}\) is the number of connected PV cells in series, the photocurrent ( \({I}_{ph}\) ) which is produced by the PV module is influenced by the solar irradiance falling on the PV module plane and its influenced by the ambient temperature of the PV module as shown below 26,27 : Where \({I}_{phs}\) is the photocurrent at standard conditions (STC), \({I}_{rrs}\) is the irradiation at STC and is 1000 \(\frac{W}{{m}^{2}}\) , \({ K}_{i}\) is the temperature coefficient, \(\varDelta T=T-{T}_{s}\) is the difference between the PV temperature ( \(T\) ) and STD temperature ( \({T}_{s}=25℃\) ), from equations (4) and (5) the output current ( \({I}_{PV}\) ) depends on the amount of sunlight and the PV temperature, the characteristics curves(P-V), and(I-V) for different irradiance and temperature are shown in Fig. 2 , PV System Modeling The modeling process for dynamical systems is used to drive mathematical representation for the system, it depends on two techniques: mathematical representation is derived from physics or measured input-output data, in general, the modeling process can be organized into three techniques: white box, grey box, and black box modeling 28 , in white box modeling the physics of the system, are completely known; and the mathematical model can be constructed entirely from prior knowledge and physical insight 29 , in grey box modeling, some of the physical parameters are known, the other parameters are determined from measured input-output data, but the black box modeling is used when the physical properties of the model are unknown, and the model is constructed from measured input-output data, but the chosen model structure belongs to the candidate models that are known to have good flexibility and have been successful in the past 30 . The back box modeling is chosen for the system identification methodology in this paper, in the back box modeling method, many models are compared with the real data, and the best is determined by an optimization criterion. Four basic stages characterize the system identification technique (I) The measured input-output data, (II) the proposed model, (III) a criterion of fit and parameter estimation, and (IV) Validate the model which is obtained. 4.1. The measured input-output data. The data used in system identification is obtained from the PV system, shown in figure .1, where the input data is the PV current and output data is the PV system output power. 4.2. The proposed model structure . a fractional order model is considered in this paper, A generalized fractional-order dynamic system can be represented by the next Eq. 3 1 : The dynamic system in (6) can be written in the next form: $$\sum _{k=0}^{n}{a}_{k}{\mathcal{D}}^{{\alpha }_{k}}\mathcal{Y}\left(t\right)= \sum _{k=0}^{m}{b}_{k}{\mathcal{D}}^{{\beta }_{k}}\mathfrak{u}\left(t\right) \left(7\right)$$ Where( \({a}_{k}\) , \({ b}_{k})\in \mathcal{ }\mathcal{R}\) , and ( \({\alpha }_{k},{\beta }_{k}\) ) \(\in {\mathcal{R}}_{+}^{2}\) . in (6) if all the derivation orders are integer multiples of a base order q such that \({\alpha }_{k}\) , \({\beta }_{k}=qk\) , \(q\in {\mathcal{R}}_{ }^{+}\) , the system is said to be of commensurate order. If \(q=\frac{1}{\mathcal{g}}\) , \(\mathcal{g}\in {\mathbb{Z}}_{+}\) , the system is said to be of rational order. The Laplace transform of (7) can be written in the form: $$\sum _{k=0}^{n}{a}_{k}{s}^{{\alpha }_{k}}Y\left(s\right)= \sum _{k=0}^{m}{b}_{k}{s}^{{\beta }_{k}}U\left(s\right), \left(8\right)$$ Then the fractional order transfer function can be written in the form: $$G\left(s\right)=\frac{Y\left(s\right)}{U\left(s\right)}=\sum _{k=0}^{m}{b}_{k}{s}^{{\beta }_{k}}/\sum _{k=0}^{n}{a}_{k}{s}^{{\alpha }_{k}}, \left(9\right)$$ By considering the system is commensurate order q , We can write \(\text{ɦ}={s}^{q}\) , then the continuous time rational transfer function is: $$H\left(\lambda \right)= \frac{\sum _{k=0}^{m}{b}_{k}{ɦ}^{k}}{\sum _{k=0}^{n}{a}_{k}{ɦ}^{k}} . \left(10\right)$$ In fractional calculus the differential operator has many definitions, The three most frequently used definitions for the general fractional differ integral are the Grunwald-Letnikov (GL) definition, the Riemann-Liouville (RL), and the Caputo definitions 32 , Grünwald-Letnikov definition is considered in this work for the drive numerical solutions to fractional-order differential equations, which is given by the next equation: $${}_{\text{a}}{}^{\text{ }}{\text{D}}_{\text{t}}^{ \alpha }=\underset{h\to 0}{\text{lim}}\frac{1}{{h}^{\alpha }} \sum _{j}^{\left[\frac{t-a}{h}\right]}{\left(-1\right)}^{k}\left(\begin{array}{c}\alpha \\ k\end{array}\right)f\left(t-kh\right), \left(11\right)$$ Where \(h\) is the step size, and the binomial coefficient can be expressed with Gamma function \(\varGamma \left(.\right)\) in the following expression: $$\left(\begin{array}{c}\alpha \\ k\end{array}\right)=\frac{\alpha !}{k!\left(\propto -k\right)!}=\frac{\varGamma \left(\alpha +1\right)}{\varGamma \left(k+1\right)\varGamma \left(\alpha -k+1\right)} ,$$ 12 Laplace transform of the Grünwald-Letnikov fractional operator is: $$\mathcal{L}\left[{\mathcal{D}}^{\alpha }f\left(t\right)\right]={S}^{\alpha }F\left(S\right)$$ 13 . Oustaloup’s recursive filter is an approximation for fractional differentiator and integrator, in a specified frequency range ( \({\omega }_{b},{\omega }_{h}\) ) and is capable of producing a very good fitting to the fractional-order elements 33,34 , the algorithm for Oustaloup’s recursive filter is shown in the following form: $${S}^{\alpha }\approx \prod _{K=1}^{N}\frac{S+{\omega }_{k}^{\text{'}}}{S+{\omega }_{k}} \left(14\right)$$ Where \({\omega }_{k}^{\text{'}}= {\omega }_{b}* {\omega }_{u}^{\frac{2k-1-\alpha }{N}}\) , \({\omega }_{k}= {\omega }_{b}* {\omega }_{u}^{\frac{2k-1+\alpha }{N}}\) , \(K={\omega }_{h}^{\alpha }\) , \({\omega }_{u}= \sqrt{\frac{{\omega }_{h}}{{\omega }_{b}}}\) And N is the approximation order in the frequency range ( \({\omega }_{b},{\omega }_{h}\) ) 4.3. Criterion of fitting and parameters estimation The system identification method depends on the error minimization between model output and system output, the least-square method is used in this work, and the problem is expressed as follows: $$\mathcal{F}\left(x\right)=\sum _{p=1}^{n}{℮}_{p}$$ 15 Where \({℮}_{p}={y}_{p}-{\widehat{y}}_{p}\) and \({\widehat{y}}_{p}=\widehat{\mathcal{H}}\left({\mu }_{p},\theta \right)\) is the output of the estimated model \(\widehat{\mathcal{H}}\) for the input signal ( \({\mu }_{p}\) ), \(p\) is the iteration, and \(\theta\) is the estimated parameters of the model \(\widehat{\mathcal{H}}\) , from Eq. (6) the model parameter \(\theta\) is formed by: \({a}_{r}=\left[{a}_{0} {a}_{1}\dots { a}_{m}\right]\) , \({\alpha }_{r}=\left[{\alpha }_{0} {\alpha }_{1}\dots { \alpha }_{m}\right]{b}_{\mathcal{F}}=\left[{b}_{0} {b}_{1} \dots {b}_{n}\right]\) , \({\beta }_{\mathcal{F}}=\left[{\beta }_{0} {\beta }_{1} \dots { \beta }_{n}\right]\) . In this work the FOMCON toolbox for MATLAB/Simulink is used, the toolbox is based on fractional calculus for system modeling and control design. The FOMCON toolbox uses the Levenberg-Marquardt algorithm for optimization of the model parameters θ under the nonlinear least squares technique. The Levenberg-Marquardt algorithm is described by the Eq. ( 16 ) According to the references 21,35,36 : $${-J}_{p}^{T}{℮}_{p}=\varDelta {\theta }_{p} \left({J}_{p}^{T}{J}_{p}+\lambda I\right)$$ 16 Where \(\lambda\) is a positive scaler, called a damping factor, \(I\) is the identity matrix, and \(J\) is the Jacobian matrix of the function \(\mathcal{F}\) and is defined as $$J=\left[\begin{array}{ccc}\frac{\partial {℮}_{\text{1,1}}}{\partial {\theta }_{1}}& \frac{\partial {℮}_{\text{1,1}}}{\partial {\theta }_{2}} \dots & \frac{\partial {℮}_{\text{1,1}}}{\partial {\theta }_{N}}\\ \begin{array}{c}\frac{\partial {℮}_{\text{1,2}}}{\partial {\theta }_{1}}\\ \dots \end{array}& \begin{array}{c}\frac{\partial {℮}_{\text{1,2}}}{\partial {\theta }_{1}}\\ \dots \end{array} \dots & \begin{array}{c}\frac{\partial {℮}_{\text{1,2}}}{\partial {\theta }_{N}}\\ \dots \end{array}\\ \frac{\partial {℮}_{1,m}}{\partial {\theta }_{1}}& \frac{\partial {℮}_{1,m}}{\partial {\theta }_{1}} \dots & \frac{\partial {℮}_{1,m}}{\partial {\theta }_{N}}\end{array}\right] \left(17\right)$$ the update parameters are defined in Eq. ( 18 ). $${\theta }_{p+1}={\theta }_{p}+\varDelta {\theta }_{p}$$ 18 From equations ( 16 ) and ( 18 ) the update parameters can be determined by Eq. (19): $${\theta }_{p+1}={\theta }_{p} - \frac{{-J}_{p}^{T}{℮}_{p}}{\left({J}_{p}^{T}{J}_{p}+\lambda I\right)}. \left(19\right)$$ 4.4. The model quality The quality of the identified model will be discussed in this section. let \({y}_{s}\) denote The output data of the system, and \({y}_{d}\) denote The output data of the identified model, for single input single output (SISO)system, both \({y}_{s}\) and \({y}_{d}\) are vectors of dimension N X 1 , and the difference between \({y}_{s}\) and \({y}_{d}\) is the model output error since \({y}_{s}\) and \({y}_{d}\) are vectors, then the error is a vector of size N X 1 which is called residuals $$e= {y}_{s}-{y}_{s}$$ 20 . The percentage fit can be described as $$fit=\left(1-\frac{\parallel e\parallel }{{y}_{s}-\stackrel{-}{{y}_{s}}}\right)*100\% , \left(21\right)$$ where \(\parallel .\parallel\) is the Euclidean norm, and \(\stackrel{-}{{y}_{s}}\) is the mean value of the system output data \({y}_{s}\) 4.5 stability analysis To study the fractional system stability, we consider Matignon’s stability theorem 37 : Theorem 1 ( Matignon’s stability theorem ) The fractional transfer function G(s) = Z(s) = P (s) is stable if and only if the following condition is satisfied in the σ-plane: $$\left|\text{arg}\left(\sigma \right)\right|>q\frac{\pi }{2} , \forall \sigma \in \complement , P\left(\sigma \right)=0$$ 22 , where \(0 < q < 2\) and \(\sigma ={s}^{q}\) . When \(\sigma =0\) is a single root of \(P \left(s\right),\) the system cannot be stable. For q = 1, this is the classical theorem of pole location in the complex plane: no pole is in the closed right plane of the first Riemann sheet. The algorithm for determining the system stability in (10) is shown in figure (4), according to Matignon’s stability theorem, the fractional order system stability in the complex plane is shown in figure.5. Iv. Experiments And Results From Fig. 1 , The PV current \({I}_{PV}\left(n\right)\) and output power \({P}_{PV}\left(n\right)\) ) are used data of the system identification, the measured data is collected under constant irradiation and constant dc load during a 1minute period, for the implementation of the proposed technique, an ARDUINO ATMEGA2560 has been used as a microcontroller, a circuit of the voltage divider and a hall effect current sensor ACS7-12 has been used for reading the analog date for the PV voltage and current respectively, in order to store the signals, a HANTEK 2D42 oscilloscope has been used, to measure the PV current or the output current a shunt resistor of 0.5 ohms is used, the data has been collected with a sample period of \({T}_{s}=1000\mu S\) . The collected data are provided for the FOMCON toolbox in MATLAB to develop the fractional order model. The measured data of the PV current and the output power of the boost converter collected from the previous experiment are shown in figure(6): The PV current and the output power of the boost converter have been used as input data and output data respectively for the FOMCON toolbox in MATLAB, the Levenberg-Marquardt optimization algorithm is used for the optimization problem, and Oustaloup recursive filter approximation is used for model simulation in the time domain, the fractional order model which is estimated by the FOMCON toolbox is shown in the following equation: Conclusion $$G\left(S\right)=\frac{{109.57 S}^{2.5684}+{4374.9.7S}^{1.7303}+3.0654}{{8.0417 S}^{2.8057}+{226.94S}^{1.7615}+{9.347 S}^{1.2891}+1}, \left(23\right)$$ From Eq. (21) The fitting between the system output data and the output of the model (22) is 93.15%, and the error norm is 213.714, as shown in the figure.7, following theorem 1 , the poles of the estimated model in (23) are satisfied with the condition (22), such that the identified model is stable for q = 0.08 as shown in Fig. 8 . Model validation Checking the model validation by comprising the proposed model with another integer model is constructed by MATLAB System Identification TOOLBOX, with the same experimental data. the selection of the integer model was dependent on many iterative processes, the best-fitted integer model shows an accuracy of 88.74%, and the fitting between the measured output data and the output of the integer model is shown in Fig. 9 . Comparing the identified fractional order model (FOM) and the integer order model (IOM), shows the accuracy of FOM is higher than the IOM accuracy. it is clear from figure.7 and figure.9 that the performance of the FOM is better than the IOM, such that the results show the validity of the proposed technique to construct the FOM of the PV system. Conclusion In this paper, a dynamical fractional order model of the PV system is driven based on the system identification technique. The measured input-output data are collected for parameters estimation based on the Levenberg-Marquardt algorithm; the identification process is based on the Least Squares technique for minimization of the errors of the model output. the fractional model stability is checked by Matignon’s stability theorem. an integer order model is constructed and its performance is estimated to demonstrate the fractional order model validation, the fractional order model shows higher accuracy (93.15%) than the other model. Declarations Ethics approval and consent to participate: Not applicable Competing interests: The authors declare that they have no conflicts of interest Funding: This research received no external funding. Consent for publication : The Authors hereby consents to the publication of the Work in the Energy, Sustainability and Society journal. Availability of supporting data: No data were used to support this study. Authors' contributions: R.E. and E.N. conceived and designed the experiments; R.E. and B.A. performed the experiments; R.E. and E.N. prepared figures and tables; R.E and B.A., and E.N. reviewed the manuscript; supervision B.A. E.N.; and R.E. wrote the main manuscript text. Acknowledgments: I wish to thank Prof. Belal Abou-Zalam, Dr. Essam Nabil, and Prof. Mohmed Hamdy (faculty of electronic engineering, Minoufia University) for supporting in research. References Rezk, H., et al.: A novel statistical performance evaluation of most modern optimization-based global MPPT techniques for partially shaded PV system. Renew. Sustain. Energy Rev. 115 , 109372 (2019) Ebrahimi, S.M., Salahshour, E., Malekzadeh, M.: Francisco Gordillo. Parameters identification of PV solar cells and modules using flexible particle swarm optimization algorithm. Energy. 179 , 358–372 (2019) Bollipo, R.B., Mikkili, S., Bonthagorla, P.K.: Critical Review on PV MPPT Techniques: Classical, Intelligent and Optimisation. IET Renewable Power Generation vol. 14 1433–1452 Preprint at (2020). https://doi.org/10.1049/iet-rpg.2019.1163 Pal, R.S., Mukherjee, V.: Metaheuristic based comparative MPPT methods for photovoltaic technology under partial shading condition. Energy. 212 , 118592 (2020) Li, H., Yang, D., Su, W., Lu, J., Yu, X.: An Overall Distribution Particle Swarm Optimization MPPT Algorithm for Photovoltaic System under Partial Shading. IEEE Trans. Industr. Electron. 66 , 265–275 (2019) Abdelaty, A.M., Radwan, A.G., Elwakil, A.S., Psychalinos, C.: Transient and Steady-State Response of a Fractional-Order Dynamic PV Model under Different Loads. J. Circuits Syst. Computers. 27 , 1–25 (2018) da Rocha, M.V., Sampaio, L.P., da Silva, S.A.: O. Comparative analysis of MPPT algorithms based on Bat algorithm for PV systems under partial shading condition. Sustain. Energy Technol. Assess. 40 , 100761 (2020) Ghasemi, M.A., Forushani, M., Parniani, M.: Partial shading detection and smooth maximum power point tracking of PV arrays under PSC. IEEE Trans. Power Electron. 31 , 6281–6292 (2016) Ali, A., et al.: Investigation of MPPT Techniques under Uniform and Non-Uniform Solar Irradiation Condition-A Retrospection. IEEE Access. 8 , 127368–127392 (2020) Yang, B., et al.: Novel bio-inspired memetic salp swarm algorithm and application to MPPT for PV systems considering partial shading condition. J. Clean. Prod. 215 , 1203–1222 (2019) Ali, K., et al.: Robust integral backstepping based nonlinear MPPT control for a PV system. Energies (Basel). 12 , 1–20 (2019) Pillai, D.S., Rajasekar, N.: Metaheuristic algorithms for PV parameter identification: A comprehensive review with an application to threshold setting for fault detection in PV systems. Renew. Sustain. Energy Rev. 82 , 3503–3525 (2018) Rodrigues, M.P., A MISO Nonlinear Model of Photovoltaic Panel Based on System Identification:. IEEE PES Conference on Innovative Smart Grid Technologies, ISGT Latin America 2019 (2019) doi: (2019). 10.1109/ISGT-LA.2019.8895325 Humada, A.M., et al.: Modeling of PV system and parameter extraction based on experimental data: Review and investigation. Sol. Energy. 199 , 742–760 (2020) Dizqah, A.M., Maheri, A., Busawon, K.: An accurate method for the PV model identification based on a genetic algorithm and the interior-point method. Renew. Energy. 72 , 212–222 (2014) Mohamed, M.A.E.H.: The linear model of a PV moduel. Int. J. Power Electron. Drive Syst. 8 , 900–906 (2017) Elazab, O.S., Hasanien, H.M., Elgendy, M.A., Abdeen, A.M.: Parameters estimation of single- and multiplediode photovoltaic model using whale optimisation algorithm. IET Renew. Power Gener. 12 , 1755–1761 (2018) Rasheed, M.S., Shihab, S.: Modelling and Parameter Extraction of PV Cell Using Single-Diode Model. Adv. Energy Convers. Mater. 96–104 (2020). doi: 10.37256/aecm.122020550 di Piazza, M.C., Luna, M., Vitale, G.: Dynamic PV model parameter identification by least-squares regression. IEEE J. Photovolt. 3 , 799–806 (2013) Alqahtani, A., Alsaffar, M., El-Sayed, M., Alajmi, B.: Data-Driven Photovoltaic System Modeling Based on Nonlinear System Identification. Int. J. Photoenergy. 2016 , 9–12 (2016) Gavin, H.P.: The Levenberg-Marquardt Algorithm For Nonlinear Least Squares Curve-Fitting Problems. Duke University 1–19 (2019) Luo, A.C.J., Gheorghe, C.M., Machado, D.J.T.: Introduction to Fractional Dii erential Equations . http://www.springer.com/series/11433 Tepljakov, A.: Fractional-order modeling and control of dynamic systems .Automaticavol. 21 (2016) Teh, C.J.Q., Drieberg, M., Soeung, S., Ahmad, R.: Simple PV Modeling under Variable Operating Conditions. IEEE Access. 9 , 96546–96558 (2021) Blaifi, S., ali, Moulahoum, S., Taghezouit, B., Saim, A.: An enhanced dynamic modeling of PV module using Levenberg-Marquardt algorithm. Renew. Energy. 135 , 745–760 (2019) Nasser, K.W., Yaqoob, S.J., Hassoun, Z.A.: Improved dynamic performance of photovoltaic panel using fuzzy logic-MPPT algorithm. Indonesian J. Electr. Eng. Comput. Sci. 21 , 617–624 (2020) Al-Soeidat, M., Cheng, T., Lu, D.D.C., Agelidis, V.G.: Experimental study of static and dynamic behaviours of cracked PV panels. IET Renew. Power Gener. 13 , 3002–3008 (2019) Celikel, R., Gundogdu, A.: System identification-based MPPT algorithm for PV systems under variable atmosphere conditions using current sensorless approach. Int. Trans. Electr. Energy Syst. 30 , 1–21 (2020) Sjöberg, J., et al.: Nonlinear black-box modeling in system identification: a unified overview. Automatica. 31 , 1691–1724 (1995) Alqahtani, A., Marafi, S., Musallam, B.: Abd El Khalek, N. E. D. Photovoltaic Power Forecasting Model Based on Nonlinear System Identification. Can. J. Electr. Comput. Eng. 39 , 243–250 (2016) Hcheichi, K., Bouani, F.: Comparison between commensurate and non-commensurate fractional systems. Int. J. Adv. Comput. Sci. Appl. 9 , 685–691 (2018) Marouani, H., Hergli, K., Dhahri, H., Fouad, Y.: Implementation and Identification of Preisach Parameters: Comparison Between Genetic Algorithm, Particle Swarm Optimization, and Levenberg–Marquardt Algorithm. Arab. J. Sci. Eng. 44 , 6941–6949 (2019) Li, J., Liu, Y., Shan, C., Dai, C.: Implementation of Simplified Fractional-Order PID Controller Based on Modified Oustaloup’s Recursive Filter. J. Shanghai Jiaotong Univ. (Science). 25 , 44–50 (2020) Colín-Cervantes, J.D., et al.: Rational approximations of arbitrary order: A survey.Fractal and Fractional 5 , (2021) Islam, S.U., et al.: Design of robust fuzzy logic controller based on the levenberg marquardt algorithm and fault ride trough strategies for a grid-connected PV system.Electronics (Switzerland) 8 , (2019) Marouani, H., Hergli, K., Dhahri, H., Fouad, Y.: Implementation and Identification of Preisach Parameters: Comparison Between Genetic Algorithm, Particle Swarm Optimization, and Levenberg–Marquardt Algorithm. Arab. J. Sci. Eng. 44 , 6941–6949 (2019) D, M.: Generalized fractional differential and difference equations: stability properties and modelling issues. Proceedings of the Mathematics Theory of Networks and Systems Symposium 503–506 (1998) Cite Share Download PDF Status: Under Review Version 1 posted Reviewers agreed at journal 23 Jul, 2022 Reviewers invited by journal 22 Jul, 2022 Editor invited by journal 22 Jul, 2022 Editor assigned by journal 18 Jul, 2022 First submitted to journal 16 Jul, 2022 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-1865340","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":123251748,"identity":"ac31c015-bf01-462b-9a51-3287a2b1d4cb","order_by":0,"name":"rabiaa ebead","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABT0lEQVRIie2Qz0sCQRTHHyzspRn2OqDs/gsjC0Uka3/KDAOeVhG67CFoIdBTdw+B/0Ii2HVkIC/bfSEPSmAXDUOIoNBm+4G4pl0j9nMYZh7vw3vfAcjI+JMgUABSX4SE4apso+Qku5UyA7YquzsV41Px6ZoC25SDxm1P1YKBbTX8pyGDYtU53uP3tYDmrdDo3SEoeuG6ko+qTDWjsUuiaYcyKJ8UJO64zYgiIk1xhKAsUgrR+yhcVzyMK13CF4q3Q9zN4fobAon2czqpSC1GrIlWFuqsFfsPLwyWiXL9iusUOdJ61spyQyHJlFAxGvumji95C3DXSBQqkakV6W0oj0yhm3GhneRhIPiVgTs5pLMUlOkeXlLB0opVUXN0OnDsvhjNZuDxVuOiPUcBLdn981E8CbzSDx+9DlXfN+PjCTzc2vuFk+74fUpGRkbGP+cdcSFzH1mo8yoAAAAASUVORK5CYII=","orcid":"https://orcid.org/0000-0002-8723-6803","institution":"Menoufia University","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"rabiaa","middleName":"","lastName":"ebead","suffix":""},{"id":123251749,"identity":"9e5ca0c8-2d56-41fb-9674-9bbd8ac4a72f","order_by":1,"name":"Belal Abou-Zalam","email":"","orcid":"","institution":"Menoufia University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Belal","middleName":"","lastName":"Abou-Zalam","suffix":""},{"id":123251750,"identity":"d00284c6-2bac-4a5f-ab24-92e357930571","order_by":2,"name":"Essam Nabil","email":"","orcid":"","institution":"Menoufia University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Essam","middleName":"","lastName":"Nabil","suffix":""}],"badges":[],"createdAt":"2022-07-16 19:37:13","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-1865340/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-1865340/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":24531439,"identity":"541926d6-840d-441f-bbc9-051912dd5be0","added_by":"auto","created_at":"2022-07-29 20:14:47","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":57863,"visible":true,"origin":"","legend":"\u003cp\u003e\tproposed PV system\u003c/p\u003e\u003cp\u003e\u003cbr\u003e\u003c/p\u003e","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-1865340/v1/2e28b48ddcb4e3227cc4d856.png"},{"id":24531202,"identity":"3787c20f-b869-45aa-b2ee-e00eba8067c6","added_by":"auto","created_at":"2022-07-29 20:09:47","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":47962,"visible":true,"origin":"","legend":"\u003cp\u003e\tPV module equvilant circuit\u0026nbsp;\u0026nbsp;\u003c/p\u003e\u003cp\u003e\u003cbr\u003e\u003c/p\u003e","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-1865340/v1/8213ba9798c24a5727206fc0.png"},{"id":24531598,"identity":"653551b0-f8d3-49f9-9c63-13d3de9a7415","added_by":"auto","created_at":"2022-07-29 20:19:47","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":183466,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003eP-V and I-V characteristics for 60watt PV nodule.\u003c/em\u003e\u003c/p\u003e","description":"","filename":"3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-1865340/v1/6ae9464ab5d9e04d12383e1b.jpg"},{"id":24531199,"identity":"af29d5ae-d114-411b-ba64-8c9f9adb8378","added_by":"auto","created_at":"2022-07-29 20:09:47","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":43625,"visible":true,"origin":"","legend":"\u003cp\u003eThe algorithm for determining the system stability\u003c/p\u003e","description":"","filename":"4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-1865340/v1/9eed1a89ebc80f4165310a5d.jpg"},{"id":24531206,"identity":"2ea25fb4-7eb0-4665-8906-59366cea9d40","added_by":"auto","created_at":"2022-07-29 20:09:47","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":55619,"visible":true,"origin":"","legend":"\u003cp\u003e\tStability of the fractional-order system\u0026nbsp;\u003c/p\u003e\u003cp\u003e\u003cbr\u003e\u003c/p\u003e","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-1865340/v1/96f5705118101d5b74a0fc6d.png"},{"id":24531599,"identity":"3c70cb1c-e103-458a-baac-6e8ebe1b818f","added_by":"auto","created_at":"2022-07-29 20:19:47","extension":"jpeg","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":249406,"visible":true,"origin":"","legend":"\u003cp\u003e\tMeasured input-output data\u003c/p\u003e\u003cp\u003e\u003cbr\u003e\u003c/p\u003e","description":"","filename":"6.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-1865340/v1/617035e1bd90dedb4867b487.jpeg"},{"id":24531441,"identity":"fb272658-9f03-4ab6-81fd-1c39b16e276b","added_by":"auto","created_at":"2022-07-29 20:14:47","extension":"jpg","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":69908,"visible":true,"origin":"","legend":"\u003cp\u003eMeasured and simulated model output\u003c/p\u003e","description":"","filename":"7.jpg","url":"https://assets-eu.researchsquare.com/files/rs-1865340/v1/899c4506c74f3b5f22b9de59.jpg"},{"id":24531634,"identity":"7fdd9e95-6e85-4448-9dd5-018d9015cbe1","added_by":"auto","created_at":"2022-07-29 20:24:47","extension":"jpg","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":19751,"visible":true,"origin":"","legend":"\u003cp\u003e\tStability of the identified model\u003c/p\u003e\u003cp\u003e\u003cbr\u003e\u003c/p\u003e","description":"","filename":"8.jpg","url":"https://assets-eu.researchsquare.com/files/rs-1865340/v1/74e50b73e42a0827284acaa9.jpg"},{"id":24531207,"identity":"372c2da0-07d9-4925-877b-eb4c1f2d82b3","added_by":"auto","created_at":"2022-07-29 20:09:47","extension":"jpg","order_by":9,"title":"Figure 9","display":"","copyAsset":false,"role":"figure","size":32552,"visible":true,"origin":"","legend":"\u003cp\u003e\tMeasured and simulated model output for IOM\u003c/p\u003e\u003cp\u003e\u003cbr\u003e\u003c/p\u003e","description":"","filename":"9.jpg","url":"https://assets-eu.researchsquare.com/files/rs-1865340/v1/634f1804f36c6555d42d5378.jpg"},{"id":24531635,"identity":"f1472e0e-3e87-4917-881a-5661e5ef10ee","added_by":"auto","created_at":"2022-07-29 20:24:50","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":734202,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-1865340/v1/3fc5b908-d2d0-40b4-8e56-2527e016cb19.pdf"}],"financialInterests":"","formattedTitle":"System identification of photovoltaic system based on fractional-order model","fulltext":[{"header":"I. Introduction","content":"\u003cp\u003eThe generation of electrical energy has been dependent on fossil fuels for many years, such as coal, natural gas, and oil \u003csup\u003e1\u003c/sup\u003e, which cause environmental issues, and fossil fuel reserves decrease. These problems have made electricity generation go towards renewable energy sources such as gravitational energy, geothermal energy, wind energy, tides, and solar energy \u003csup\u003e2,3\u003c/sup\u003e. renewable energy sources have been taking a significant position in the replacement or complement of conventional power generation, but solar energy source among all renewable energy sources is used widely because of advantages such as having no noise due to the absence of mechanical parts, low maintenance cost, long working life, the great amount of solar irradiation that covers the most of the earth\u0026rsquo;s surface \u003csup\u003e4\u003c/sup\u003e. The electric energy is generated from solar energy by a small unit called a photovoltaic (PV) cell, which produces a little amount of energy, several PV cells are connected in parallel and series to form a PV module to have more power, and PV panel is a group of PV modules which are joined together in parallel or series form. In a PV system, PV panels are connected in series and parallel to form PV arrays for getting high current and high voltage, to increase the output power of PV arrays \u003csup\u003e5\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eThe productivity of PV modules to convert solar energy to electric energy is dependent on the solar irradiance level, the surrounding temperature, and others \u003csup\u003e6,7\u003c/sup\u003e, so the variable environmental conditions cause variation in the voltage and the current of the PV module. therefore, the output power which is generated by the PV module changes also. Hence In a PV module, voltage and current have a nonlinear relation \u003csup\u003e8\u003c/sup\u003e, and in the electrical characteristic of the PV module, there is only one operating point, at which PV output power is maximum, this operating point is called Maximum Power Point (MPP), which makes PV modules are working efficiently. As this point fluctuates with variable environmental conditions, it should design a maximum power point tracking (MPPT)algorithm \u003csup\u003e9\u003c/sup\u003e, MPPT algorithms increase PV system efficiency by moving the PV voltage or current to make PV operate at the maximum power point (MPP).\u003c/p\u003e \u003cp\u003eAnother problem that has a great effect on the PV efficiency is Partial shading (PS), PS decreases the amount of sunlight that is allocated to the PV array resulting in lower output power from the PV array, under PS multiple local maximum power points (LMPP) emerge on PV curve \u003csup\u003e10\u003c/sup\u003e, but there is one MPP of LMPP is called global maximum power point.\u003c/p\u003e \u003cp\u003eTherefore, to overcome the above issues, many MPPT algorithms are proposed to get maximum power from PV array under variable environmental conditions, in\u003csup\u003e11\u003c/sup\u003e the MPPT algorithms are classified into conventional techniques (CTs), soft computing techniques (SCTs), linear and nonlinear control techniques, all these algorithms are used in a simulation environment before real systems, PV model has many important roles in design optimization, fault diagnosis \u003csup\u003e12\u003c/sup\u003e, and prediction performance in the future of any PV system, therefore it should develop an accurate PV model to design MPPT.\u003c/p\u003e \u003cp\u003eThe mathematical models of real systems can be divided into static and dynamic models \u003csup\u003e13\u003c/sup\u003e, the dynamical model for the PV system is preferred due to the loads always are not constant, and the PWM which is used for switching converters or inverters, the static model for PV module is introduced in many pieces of literature such as in \u003csup\u003e14\u0026ndash;18\u003c/sup\u003e, a single diode model is considered as a static PV model and is the most one which is used, the static model can be derived from parameter identification, parameter identification method tray to make approximation and solve analytically the equation of the PV as a polynomial of the first order [\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e], this means the static model is simplified and has low accuracy.\u003c/p\u003e \u003cp\u003eAlternatively, a second-order dynamic model for PV module is proposed in \u003csup\u003e19\u003c/sup\u003e, the technique which is proposed increases the model accuracy, but it doesn\u0026rsquo;t consider the power electronics system, another dynamic model is introduced for a PV system based on system identification and considered buck converter, MPPT, and load in \u003csup\u003e20\u003c/sup\u003e, this technique used experimental input-output data and nonlinear Hammerstein-Wiener model to extract dynamic model. This paper proposes a dynamic fractional order model for a PV system, boost converter, MPPT, and load.\u003c/p\u003e"},{"header":"Ii. Related Work","content":"\u003cp\u003eThe PV system in this paper is off-grid, and consists of a 60-watt PV module, boost converter, MPPT controller, and load. the circuit diagram of the PV system is shown in Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003eThe PV current measurement is used as input data for the system identification and output power of the boost converter as output data. The identification method to minimize the error of model output is based on Nonlinear least-squares estimation, Levenberg-Marquardt algorithm is used to solve nonlinear least-squares problems\u003csup\u003e21\u003c/sup\u003e, and for determining the parameters of the model.\u003c/p\u003e\n\u003cp\u003eThe proposed PV model is based on the fractional order model, fractional calculus is a generalization of ordinary differentiation and integration to arbitrary (non-integer) order, it is supposed that fractional order models will generally provide a more accurate description of the system dynamics than those based on conventional differential equations because it is given more degree of freedom for determining parameters of the model, the general definition for an integral differential operator is\u003csup\u003e22\u003c/sup\u003e :\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equ1\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equ1\" name=\"EquationSource\"\u003e\u003cimg src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAApIAAABUCAYAAAAxptfmAAAVzElEQVR4nO3d/0sbef4H8Gc+3K+tO8n+VFopGYUTCylnomXXFlowmVUQyu2e6VkWobKp2UVoWbRrLPdDq1tluYXlLmuOCmW5rbNXlx5Hm3MsVLBzZdV0cViPFjRDkV1/csy5/gH5/CAzGI1JjCZR+3xAQGcmk5ejbZ55z/uLLZlMJkFFo+s61tbW4HK5oKoq6uvrS10SERERUV5+U+oC3iSyLOPSpUtwu9344IMP8PTpU4yNjZW6LCIiIqK82NgiWTiSJEFRFDx+/BiNjY2w2+0Ih8M4fvw4mpubEYvF4HQ6S10mERERUV7+r9QFHCSGYWBwcBB+vz+n4zs7OwEAdXV1UFUViUQCfr8fP//8M2praxkiiYiI6EBjkMyRYRiQJAk//vgj+vv7c3rO3NwcRFGEw+GwtkWjUfT29sIwDGiaVqhyiYiIiAqOQTJHra2tEEURsixnbEmMRqOoqKiA3W7HgwcP4PV6AQD19fVoaWnBo0eP8P3338PhcOCXX34pVvlEREREe46DbXIQiUQwPT2N+fn5jMepqoqmpib09PSgr68PNpsN7e3t1n5Zlq2vOciGiIiIDjq2SGZhGAY+++wz3LhxI+UWdTr//ve/IQgC+vr6oOs6AODMmTPFKJOIiIio6BgksxgdHUUikcCVK1eyHvvixQuIogjDMDA6OgoAWFtbQzQaLXSZREREREXH6X+y8Hg8cDgcOd2KjkajaGpqQkdHBz799FO43W74/X6Ew+EiVEqFFIlE8PDhQwBATU0Nrl+/nrWFmoiI6LBjkMxA13WIooiBgQF0dXWVuhzKQyQSwZkzZ+ByufI+hzkf6NDQEADg6tWrEEURU1NTDJNERPRG463tDF69egUAOHXqVIkroXydOHEC7e3t8Hg8kGUZhmHs6PmyLENRFLjdbgQCAQQCAbjdbsTjcav7AhER0ZuKQTKDubk5AMDRo0dLXAkZhoFIJAJJkiBJEqLRqDW3pyRJ1uCmzRobGzEzM4O7d+9icnISlZWVCIVC2x6/2eTkJID1Lg4mc0on81Y3ERHRm4pBMoPXr18DAKqqqkpcyZtN0zRUVlbi6tWruHDhAjo7O9HU1ITh4WEoioKLFy9mXSXI5XIhHA5jfn4e5eXl8Hq9kCQpZUqmdMzAefLkSWtbWVnZ7n8oIiKiQ4DzSGZghohM/eBsNtuevR67q6bX3d2NRCKR0lfV5/Ohu7sbHR0dCAQCOZ/L4XBYt6hVVUVbWxvu3buX17yeO71NTkREdNgwSO4Sw19haZoGRVEAIGUKJkVRIIoibt26ldd5ZVnGvXv3AABtbW1Zj19dXd2yjQNtiIjoTcdb27Sv/fDDDwDWWyDN4Gbejq6oqNhRmNN1HaFQCHa7HZOTk+jt7cXCwgL8fv+2z6mpqQGwPkeoyezyYO4jIiJ6UzFI0r6gaRokSYKmaVmP++mnn6zvzUE4mUSjUfj9fni9XpSVlWF+fh7hcBj19fVZ63rvvfcAANPT0zAMA4ZhYHx8HABymqSeiIjoMGOQpH1hbW0NiqJgbW0tZXtDQwMEQYCiKJAkCf/4xz/Q19cHQRAwPT2N1tbWjMtQ+v1+fPPNN/jwww+xsLCArq6uHbVi1tfXY2RkBIlEAnV1dairq8PKygpGRkayDvAhIiI67DgheQbmRNTxeHxHoUGSJADAt99+mza0GIaB1tZWKIqCx48f4/Lly4jFYm90MFFVFc3NzVhZWdmyzzAMvHz5EseOHbOuka7rWFpaytqqaBjGnvRlNGsAkFNLJhER0ZuAg20yqKmpgaIoWFpa2lHIe+uttxCLxbYNMDdv3gQALC8vo7KyEjdu3HijQyQA3L59GxMTE2n3ORyOLeHN6XTmdM32akBMuhqIiIjedLy1ncFu5gt0u93W16qqQpIka7oYWZbR2dkJAEgkElx+Eeutt7tZxpCIiIiKj0Eyg/LycgDA8+fPsx6rqir8fj+CwSDi8TjOnz8PYD00Njc3Q1EUtLa2Qtd1JBIJHD161GqZJE6lQ0REdBAxSGZw/PjxnI6LRqM4e/YsRFGEy+VCLBZDdXU1AKC2thYA0NLSgt7eXjidTvh8PjQ3N+PcuXMQBMHqU0lERER0kHCwTQa6rkMURfh8vowrn5i3rWdmZhCNRtHU1GRNVG6eY2RkJON8hUREREQHDVskM8h1AIyiKPjggw9gGAYePXoEn8+HaDQKwzDw5MkTAOvT2ESj0UKWS0RERFRUDJJZiKKI6enpjMf4fD7cuXMHw8PDcLlcmJ6exn/+85+Ufn+SJOHXX38tdLlERERERcNb21mYc0lmukzmHIPm9DCqqqZMFaOqasociERERESHAYNkFrkESSIiIqI3EW9tExEREVFeGCSJiIiIKC/7PkgODg7CbrfDbrcjGo3CbrcjEokUtQZRFIv6ekRERHR4qaqKYDAIj8eT07EejwehUAi6rhehup0pSZCUJAk2my3l4fF4MDg4aC0jCKyHyO7ubvz9739HLBbD5cuX4ff7EQgEilrvxuUOiYiIiPIVCoVw9uxZAMDMzEzKPrPBbGPArK+vx3fffYfx8XG43W7IslzUerMpSZBsa2uzvk4mk4jH43A4HOju7k5Z5eXBgwdoaWlBY2MjlpaWkEgkEA6Hi1rr9PS0tdxhrnRdRzAY3He/bNobkUgEkiRBkiSEQqGUDz9ERETbGRwcRH9/P3w+X0qeMQwDfr8fTU1NSCQSW5YNdjqd1sIoly5dgqZpRa07k5IEyZ9++gnA+vyLwPoFMi9oLBazJu6OxWIQRRG6ruO///0vgPUmXlVVi1KnqqpIJBJoaGjI6XjDMBAKhSCKIr7++mssLi4WuEIqNkmScPXqVVy8eBEXL15Ef38/6urqGCaJiCgjwzBw584dAMDAwEDKvtbWVtjtdnR0dGz7fIfDYe3v7u4uXKE7VJIgOT4+DgC4cOGCtW3jHItzc3MAgJ6eHvT392N4eBgNDQ0QBAH3799PmaNxs8HBwS23zdM9cnH79m309PTkNP+jpmmQJAnvvvsu+1QeUrIsQ1EUuN1uBAIBBAIBuN1uxONxjI6Olro8IiLax548eYJEIgFRFOFyuVL2jY2NIRwO4+TJkxnP8Yc//AHA+op6+6W/5G+K/YKGYSAWiwEA3nnnnbTHlJWVAQD6+vrQ19dnbV9ZWcl6/q6uLnR1de26xps3b8IwDFy/fj2n57hcLoyNjcHhcOCrr75CPB7fVQ20PcMwMDo6iocPHwIAOjs7UVdXh9bWVgBAOBxOG/43dpvYTltb27Zrok9OTgJASt8Vr9eLWCyGhw8fFr3vLhERHRzme0hFRUXe59gYQF+9erUvFjopepCcmpqyvt68+oupurq6qDVt9vbbb0MURUxNTW3pp5DJTo6l/GiahvPnzyORSGBgYACnTp1CU1MTBgYGoCgKhoaGtv2H1dvbm/X8x44d23af+elv4ydG80MPERFRJuZ7yMa7sbsxNzeHxsbGPTnXbhQ9SJq3rc3+kab79+8DWJ9qJ9Ot62zMkd7ZZFqpZnl5GTdv3oQkSVYrI+0P3d3dVog0W559Ph+6u7vR0dGRsVVwN39X2bCPJBERvYmKHiSfPn0KAKipqbG2aZqGr7/+GgDw1VdfpRwfiUSgaVrOo7X34ta2w+FAOByGJEn485//nHJ7nUpH0zQoigIAuHLlirVdURSIoohbt25lfP5ub22bVldXt2zjhw0iIspFuveQg6xgQVLXdXzxxRfQdR1vvfUWZFmGYRhWEHj33XdhGAaePHmCYDAIYH0U08Zm2lAohP7+fgiCgJWVlaJPp9Pb24uzZ8/iypUr+6Ifwpvuhx9+ALDeAmkGN/NvoqKiImuY2+2t7ZqaGiiKghcvXljbXr9+be0jIiLajpkjNr6HbGaGzIWFhaznO3Xq1N4UtlvJAojH40lBEJI+ny85NDSUNF+mpaUlCSDlIYpisqWlJfns2bMt53n8+HESQHJgYCA5OztbiFKzEgQhOTQ0tKPniKKYBJDs6OgoUFWH2+zsbNLn8235nZt/Sz6fzzqup6fH2ra8vLzj39VOPHv2LAkgKQhCcnl5Obm8vGz9ruPxeMFel4iIDr6RkREr92zW09OTdLvdW/KR+X5nMt+H9tP7TkGCpPnmbr7Z5hsCzYteyovl8/mSLS0tOR9rBgvz4Xa7kz09PQWu8nAx/6Fs/nBhfkAxg6N5XQVBsD64FPoDx8b/CERRTAqCkBwZGSnoaxIR0cG3vLxsvYelazzLxcbGk/3ClkxmGHWSJ7Mv2tjYGDRNw5EjR/K6Nez3+xGPxzE2NoaXL18WdLDEdiRJwsLCQk7NzNtNlH7kyJEtc0bR9lRVRXNzc9rpngzDwMuXL3Hs2DHrb0rXdSwtLRXt78OsASjsAB4iIjpczAHBPp/PWqkmV7quw+12I5FIYHZ2dt/kioL0kWxra8OlS5cgSRIMw9jxxdooHo9DkiTcvXt3DyvceQ25YKjYG7dv38bExETafQ6HY8t1djqdRe3Dmq4GIiKibLq6urC6uor+/n4Eg8GcBxIbhoGWlhYAwMjIyL4JkQBQkBZJYH2E7dra2q7ecM2Wn6qqqpKNipUkCYqiZJwuiPaWYRgcBU1ERIeWqqq4f/8+ZmZmMDMzk/XYa9euwev17svBvwULkocFgyQRERFReiVZa5uIiIiIDj4GySwWFhYgCEKpyyAiIiLadxgks4jH46itrS11GURERET7DoNkBuYC60RERES0FYNkBktLSwCACxcu7Oo8qqoiEonsRUlERERE+0bB1to+DH7++eddn0OWZQSDQdjtdqyurqKrq2sPKiMiIiIqPbZIZrC4uAgAeOedd/I+x+rqKhKJBNxuN8rLy/eqNCIiIqKSY5DMYHV1dUfHy7KMiooK2Gw2RKNRAMCJEycAAJ999hn8fv+e10hERERUKgySGbx48QIAcOzYsazH6rqO3t5eTE1NoaWlBUePHgUAPHr0CKIowuVypR28o6oq/H4/NE3L6TVyGQA0ODiIYDDIwUJERERUUAySOchlOaLh4WF4vV44HA7IsmwtDWmGOb/fj1evXqU8R5ZlXLt2DX/9618zrpspSRJsNhtEUdxyjnS6urrwxz/+EV6vN6eAmo4sy5AkCcFgEJqmQZKkvM9FREREhxOD5B558eIFTp48uWX7wMAAPvroI3zyySdobGy0tquqimAwiLGxsazrSm9c1P23v/1tTvXU19fj3r17+P3vfw/DMFL2SZKU9mEKhUK4dOkS2traIAgC2tvb4XQ699Ui8fuBruvw+/2w2Wyw2WzweDxQVbXUZRERERUNR23vEUVR0k4T5HK50gaw27dv4/PPP88aIgFgbW0NACCK4o4Wa6+vr4cgCBgeHk4ZLd7b25vxef39/ejp6bH6dI6Pj6eE2cPOMAyMjo5C07Rtf27DMOD1ehGPxzE7OwsAOH36NJqbmxGLxXb0eyIiIjqo2CJZArquQ1EUNDQ05HT8L7/8AgDwer07fq329nb87W9/S9lWX1+f9gHAalF77733YBgG/vnPf8LhcCAajW5p2TxszFbiyspKLC4u4tNPP9322NHRUcTjcfh8PuvDgs/nQyKRwPDwcBGrJiIiKh0GyT2w09uZ09PTALbve2kOwDFvmXZ2dgKA1bJpGAZCoRDsdrt1jPmoqKhIOVd1dTXi8XjOIbCqqgqCIKCtrQ3Dw8M4f/48pqenMTc3l1PrabFtvlahUAi6rsNms8Fut2cdcGQYBmRZhsfjwbVr13Du3DmsrKygr68vY6vixMQEAKCmpsbaZrZIm4O0iIiIDjve2s7ADBKGYexpiFpcXITP50u7LxqNoqmpCYIgYHZ2FsePH8fbb78NADhz5gyA9YE95eXlWFlZgcfjQSwWw+zsbMY+jC9fvrRaHTNxOByIxWJYW1uzzlddXZ3Tc4tt47V69uwZqqqqUFlZiXg8DgD417/+tW0Y1DQNkUgE4+Pj8Hq9uHv37o76gP7vf/8DAJSVlW3Zd9hbbomIiEwMkhmYg2dyDWF7MeG42foYDoe3BBvze7O/o67riMViEARhTwfCbA5f+zFEAqnXyqxRFEV89913GBoaylh3d3c3FhYWMD4+vuf9Gfdjyy0REVEh8NZ2BqdOnQIA/Prrrzkdf/z48ZzPna7VSlVVxONxCIJgDXQxb5una8F88uQJgPz6Th506a6VYRiIxWJwu90IBAIZnz8wMACv1wuv12tNcZSP169f5/U8IiKiw4BBMgNzqp25ubk9PW95eTlisdiW7c+fPwcA1NbWWtv+8pe/AEDaEeFmP73z589nfc2qqqq8at2v0l2r0dFRALm1CLpcLoTDYUxNTeHcuXNob2+Hx+OBLMs5vb75+9jYB/PHH38EAFy8eDG3H4KIiOiAY5DMwOl0wu124+nTp3t6XjP8bB4IYt4aN1srN4eawcHBlIE9Zj+96upqRCKRtIN+nj9/DlEUD+ztVlVVYbPZtvxsm69VJBJJ2a9pGkKhUNbzOxwO+P1+zMzM4Msvv8Tk5CTsdrs1aGc777//PgRBgKIo0DQNmqZhfHwcgiDg/fff3+mPSUREdCAxSGbR3t4ORVFyGkCxcSlFwzCgqmraMOJ0OuHz+awWNFNDQwN8Ph9isRg8Hg+qqqrw4YcfQhAE3LlzB2VlZSn9/sxWsba2NqyurqbtE/jgwQN89NFHOf+8+9WRI0dSvt94rSoqKnDixAkEAgG43W4oioJIJILr16/v6DXq6+sRDocxPz+P8vJyfPHFF9se63Q6MTExAbfbjdOnT+P06dPwer2YmJg4sKGdiIhop2zJZDJZ6iL2M8MwUFlZiRs3bqRM6r2RLMsIBoNYWVkBsD6a+PLly6itrYWiKBgZGbH68ZlUVUVzczPm5+cLFjyi0Sg6OzsxNTV1YMNNNBrFn/70J8zMzJS6FCIiItqELZJZOBwOfP7557hz5862rZKLi4spA146Ozvh9XoxNjaGjo4OBIPBLc8xW78kSSrIdDGapqGzsxPff//9gQ2RAPDNN99gbGys1GUQERFRGgySOQgEAqitrcXHH39sbVNVFR6PBwDw9OlTfPLJJ9a+eDyO3/3udwDWpxBKJBJp+y/6/X58+eWX+Pjjj/MeNZxOJBKx5kg86Otjy7J8oIMwERHRYcZ5JHP07bffQpIk+P1+9Pf348iRI4jFYrDb7bhx40becy1uXJ5wr2Sb+oaIiIhoLzBI5sjhcGBsbAzDw8Po6emBLMvI1L10dXU15fvDNv0OEREREQfbFIAkSVhYWMC9e/fQ1taGiooK9vMjIiKiQ4dBsgAMw8DNmzeh6zqcTidu3brFfn5ERER06DBIEhEREVFeOGqbiIiIiPLCIElEREREeWGQJCIiIqK8MEgSERERUV4YJImIiIgoLwySRERERJQXBkkiIiIiyguDJBERERHlhUGSiIiIiPLy/x++qJEZoo/0AAAAAElFTkSuQmCC\"\u003e\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003eWhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\alpha \\in \\mathbb{ }\\mathbb{R}\\)\u003c/span\u003e\u003c/span\u003e and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\alpha\\)\u003c/span\u003e\u003c/span\u003e can equal a complex number, there are many multiple definitions of the fractional operator, we consider the definition of Gr\u0026uuml;nwald-Letnikov, which is used throughout this paper to drive numerical solutions for fractional-order differential equations.\u003c/p\u003e\n\u003cp\u003eFractional order systems modeling and controllers are implemented through the FOMCON toolbox for MATLAB/Simulink, the toolbox Fractional Order Modeling and Control \u0026ldquo;FOMCON\u0026rdquo; is developed in \u003csup\u003e23\u003c/sup\u003e, this toolbox will be used for system identification to extract the fractional order model of the PV system.\u003c/p\u003e"},{"header":"Iii. Static Pv Model Characteristic","content":"\u003cp\u003eA single diode model can describe the static behavior of a PV module and is proposed in many kinds of literature such as in\u003csup\u003e24,25\u003c/sup\u003e, the electric circuit for the single diode model is shown in Fig. \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003e.\u003c/p\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003eFrom Fig. \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003e, the mathematical equation which describes the static PV model can be driven as follow:\u003c/p\u003e\n\u003cp\u003e\u003cimg src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAf4AAAB+CAYAAADFqHNJAAAe3UlEQVR4nO3d34tTd/4/8Od8+dxO7TG9ckVkzlxUVGapZxTaWKjgJEyhtGzbxHUpgqImuwjuLjM1M1CWttYTaBdKW2fEgpTdJtkKW6SJZgoumLNSx7EkUNGLObkQ8SrH7Dh/wPt7MX2/Pfk5yUx+jXk+IJDknJy88k7OeZ33j/POgBBCgIiIiPrC/+t2AERERNQ5TPxERER9hImfiIiojzDxExER9REmfiIioj7CxE9ERNRHmPiJiIj6CBM/ERFRH2HiJyIi6iNM/ERERH2EiZ+IiKiPMPETERH1ESZ+IiKiPsLET0RE1EeY+ImIiDrAsiyEw2GMjo42tO7o6CimpqaQz+dbGgcTfw8bGBgoufUqy7JK4vT7/S3bdjQaRTQabeo1nSy3ZrcfjUZLYmv2sxHRxjQ1NYX9+/cDAG7fvg1g5dgZDAZLjp25XA4A4PV6kUgkMDc3B8MwEI/HWxeMoJ4Vi8UEAAFAZLPZbodTVygUEgCEruuiUCi0bLumaQrTNJt6TSfLbS27kK7rAoAIBAJtiIiIeo1pmgKA8Pl86rlMJiM0TROZTEYIIUQ2mxWappWsI4QQhUJBaJrW0uMZa/w97MGDBwAAwzAwMjLS5Wjqk01RgUAAHo+nq7H0ernZtg0AeO+997ocCRG1m+M4OHfuHADANE31vNfrxePHj+H1egEAIyMjCIVCSKfTJa/3eDwIhUIAgMnJyZbExMTfw65fvw4AGBsba9t7lDeLV7s10hwtf6yvvPJK22JtVCfKba0sy1L39+3b18VIiKgTfvzxRxSLRei6vmpFpFgswjCMiuffffddACvH2Vb09//furdAbdOJZCqEWPc2ei2Z9dJJSLmbN28CWGmN6HbLCBG1340bNwAAw8PDddfL5XI4f/48stlsxTL3CcP9+/cxNDS0rphY4+9Ra02mfr8ffr+/5aNA62llMvP7/SWtDZOTk5icnGx48GCz5RaPx+H3+xtq1SgfxCgH9pU/546hXC+3RhBR68lj8YEDB2qu4zgOXnvtNcRisVVbBX755Zd1x8Qaf49aazIdGhrC+fPnGz4jbGRUummamJiYqLm8XjILh8Pqh//888/jzTffRDAYrLmta9eulTyWCbne+7s1W2579+7FoUOH6u6UktfrrWghGRgYaKrVpFZrRC6Xq+i/O3XqFMbHxzE1NYU7d+4AWCkfue6ePXvw8ccfN/zeRNR7HMfB4cOHceXKFdXf326s8feotdYMHz9+DJ/P1/D6QohVb6sl3dWa1tPpNKanp6HrOg4dOtTWS9iaLbfl5WUAwMsvv9y2mKR6rRFbt27F4uIiHMfB9PQ0HMfB66+/DgDYtm0b0um0OplbXl5GOp3uya4MIqpuaWmp4jnHceD3+zE9Pa2Sfr2KUasw8fcgx3FUMt29e3fJslQqpZry5X13QllYWMCePXsQjUZLrgltl1Qqpe6/+OKLFcvlNaher1fVTmVybrV65eY4DqampuD3+xEMBlW53Lt3T60TDAYxNTXVltgA4OrVqwAAXdertkbYto1jx47B6/XC4/FA0zQAwM6dOwEA27dvBwCcPn0apmlifHy8bbESUWvIE3bZauf2xz/+EX/5y19U0nccB3Nzc3W3t2vXrnXHxMTfgw4fPqzuX7p0Sd3P5XL473//i3Q6jUgkAmDlh/LRRx8BWOlLsm0bd+7cwa5du5BOpysuDWmlfD6PU6dOqceXL1+uWG7bNt555522xeBWq9xkbIlEAtPT07BtW9X05cCbu3fvYvPmzTh79mxbYkulUjh//jyAlQRfPg7g1q1bAICLFy9ieHgY8/Pz+OqrrwCgpPlvdnYWQONdH0TUXa+++ioAYHFxseT5VCqFRCKBQ4cOqfFBL7zwAorFYsU23MeLahWsZrGPvweV93NLIyMjGBwcxNmzZ/Hee+9hfHwcH3zwgVo+Pz8PYKVv+De/+Q0ANNXs36yhoaGKH7Pb/fv3ATxtRpc/3kb609eiVrlJtm3j7t27SCQS2LRpE4CnLRInTpxAMBhEIBBoS2zj4+N4/PhxzeU//PADNE3D3//+dwCo2te3tLSEc+fO4cqVK22JkYha7+DBg9A0TZ3wy317fHy84fFBsrXQ5/Ote0Q/AM7ct9HIWemEWJn5CYCIxWJCCCEikYhaNjMzo2aFsm27K7G64ykUCsIwDKFpWktn9muGLJNQKKRiAqBmBtQ0TczMzKjy7KRqM3a54deZCCORSAejIqJWqDZzX6Ns2+bMff3uxo0b0DQNfr8fb7zxBkKhkBoMMjc3p2r4S0tLKBaL+OKLL1pzhtgkx3GQSCRUrLIJ6z//+U/Hr1+3LAt+vx///ve/AUD1ncvmddkiUSwWcebMGWzdurWj8YXDYRSLRczPz9e9DFPXdfz5z3/uYGRE1AoTExOIRCJIp9MIh8MNv85xHNUK2cilfo0aEKIFM7hQxwwPD+P48ePw+XwYHBwsSeqWZWHLli0YGhqC4zi4d+9exy4PKZfP5/Ho0SP1eHBwsGvT58qyAKDKxx2jLKNcLldRpp2Qy+XUmINa35f7uyWijcmyLHz77be4ffu2+qOeeuuePn0aY2NjOHr0aEv3fSb+DSSfz0PXdUQiEV6/TUREa8LBfRvIpk2bkMlkuh0GERFtYKzxExER9REO7iMiIuojTPxERER9hImfiIiojzDxExER9REmfiIioj7CxE9ERNRHmPiJiIj6CBM/ERFRH2HiJyIi6iNM/ERERH2EiZ+IiKiPMPETERH1ESZ+IiKiPsLET0RE1EeY+ImIiPoIEz890yzLQjAYRC6X63YobRWNRhEOh5HP57sdChH1OCZ+embF43GcPn0aX375JUZGRrodTltNTEzg97//PcbGxp75kxyijcqyLITDYYyOjja07ujoKKamplp+Qs/E3wOi0SgGBgbULRqNdjukuvx+f0m8lmV1O6QKcge7du0aPB5PyfPu2P1+f8veMxqNNv3duWMZGBhY1/t7vV5cunQJv/vd7+A4zrq2RUStNTU1hf379wMAbt++DeBpi6Tc/0dHR9WJu9frRSKRwNzcHAzDQDweb10wgnqCrusCgAgEAt0OZVW2bQsAAoCYmZnpdjhV+Xy+mrGFQiEBQOi6LgqFQsve0zRNYZpmU6+JxWKqLLPZbEviMAyj6TiIqH1M0xQAhM/nU89lMhmhaZrIZDJCiJXjqq7rQtf1ktcWCgWhaVpLjxGs8fcI27YBAO+9916XI1ndo0eP1P233367i5FUl8/nkU6ncfDgwZrLASAQCJS0BnTDgwcPAACGYbSsO+LYsWO4cOFCS7ZFROvjOA7OnTsHADBNUz3v9Xrx+PFjeL1eAMDQ0BCGh4dhGEbJ6z0eD0KhEABgcnKyJTEx8fcAd1P5vn372vY+5V0KtW6ruXnzJoCVZNXtxFnN/Pw8gJUdqZp0Og0AeOWVVzoWUy3Xr18HAIyNjbVsmzt37oRt22zuJ+oBP/74I4rFInRdr3lyn8vlEI1G8fzzz+PLL7+sWP7uu+8CWDl2taK/n4m/B3QqkU5MTEAIseptNe1IVq304MED+Hy+qss6dZLVqHaehNy7d6/l2ySi5ty4cQMAMDw8XHW53+/Hb3/7W5w7dw4vvfRS1XXcJwz3799fd0xM/D2g1xNpubUkq6mpKfj9fqRSqVXXLR88WO221kF5rTzJKo9zcnISk5OTDcfZ7ElIPB6H3+/v+cGfRPSUrKEfOHCg6vJr165BCIF//OMfOHfuHA4fPlx3e7/88su6Y2Li7wG1Emkul4Pf7y+5ycQpE6lMLHLdqampmu/Tiqb+tdaYd+/ejXQ6jeeee27VdeWOUO927dq1utuo1cxd7yQrHA6rMg0Gg6uOoi2P0zRNmKbZcJzNnoTs3btX/VaI6NkyPj6OUCjUkX2cib/L6iXSrVu3YnFxEY7jYHp6Go7j4PXXXwcAbNu2Del0WvVjLy8vI51O162Ft6Kpf6015qWlJQBQA1naadu2bVhYWKi6bLXWinQ6jenpaei6jkOHDrW1dt1sS8/y8jIA4OWXX25o/R07dqwtMCJqOXkMrGfTpk0diISJv+uuXr0KANB1vWoitW0bx44dg9frhcfjgaZpAFYGcAHA9u3bAQCnT5+GaZoYHx9va7zfffcdAFRMQJHL5RAMBpFKpdT92dnZkuWGYSCVSsHv97f2mtQye/fuBYCKQTDuboYXX3yx4nXyelmv14uPP/4YwNPk3GqO46iTkN27d1csn52drWjlcffZB4PBmq07N2/erPl7IqLOkpWzO3fulDwfjUYxPDysjlP5fB4XLlxQI/hr2bVr17pjYuLvolQqhfPnzwNYSfDlE+HcunULAHDx4kUMDw9jfn4eX331FYDSmrNMsBMTE22NNxqNqpr03Nyc+sE6joPZ2VnYto0PPvgA9+7dw+bNm3Hy5En12rm5ORSLRTx58gSO4+DSpUtti3NoaAg+nw+XL19Wz+XzeZw6dUo9di+Ty23bxjvvvNO2uNzc/XjlZZHL5XDy5EmcOnUKQ0NDqk9PDhK6e/cuNm/ejLNnz1bd9nfffYfjx4+3KXIiasarr74KAFhcXCx5/ujRowgEAhgbG8PAwAAMw0AgEMCHH35YsQ13bqhWaWlaS2YDoLYIhUJqggc5yYMbABGJREomgegmXddFJBIRQggRiUSE/HnJCX/kMk3T2j7xj5wco9EJepLJpACgyjGTyQgATU2Es5YJfKrJZrMCgAiFQsK2bWHbthBipXwNwxBCCBEIBKpO9pRMJls+MRERrZ17Ap61Hqfl8dQ9AdB6MPH3ME3T6n7R+HXGN5lQu0km90wmo37ooVBICFGaVGVSy2QyIplMtjWmWCwmDMNoKAm6T1QKhYIwDKOpE4dWSyaTwjAMNYtXoVAoORGRJ0+xWEy9JpvNCl3XWza7FxG1RrWZ+xpl23bLZ+5j4u9RclpZTdNUja8c2jDt7FrJqWd9Pp+qmcq43ElV1qQ7FXcmkxGBQKDuDlMoFISu6+pES8bXjQRaKBSEz+crKUchKlsk5G9DPp6ZmVEtBETUe+RxUFaIGuGuhLhP8tdrQIgGZmyhjsvlcmoUd62R8JZlYcuWLTVnqOsk+ZewX331FZaXl0smnJB/OiGfsywLO3bs6JnBZ/l8vmQa4sHBwa7+m5/sz3PHIWOUv4VcLofBwcGe+O6JqDGWZeHbb7/F7du31R/11Fv39OnTGBsbw9GjR1u6rzPxU0sMDw9D07RVf8xERNRdTPzUErKW2ks1eSIiqsTET0RE1Ed4HT8REVEfYeInIiLqI0z8REREfYSJn4iIqI8w8RMREfURJn4iIqI+wsRPRETUR5j4iYiI+ggTPxERUR9h4iciIuojTPxERER9hImfiIiojzDxExER9REmfiIioj7CxE9ERNRHmPipL+XzeYTDYeRyuW6HsmbxeByzs7PdDqMEy5Wo9w0IIUS3gyDqJMdxcPjwYfzzn/+Ex+NBPp/Ho0ePAAA7duwAANy7dw8A4PV6S1738OFDjIyMVGwzl8theXlZPd6yZQuGhobqxuF+jft9mmFZFq5evYqPP/54Ta9vpfJyLf98spwHBwdLypDlStRhgqjPmKYpksmkejwzMyMACE3TRDabFUIIoeu60DSt5HWZTEb4fL6q2wwEAgKAMAxDGIYhAJS8R73X6Lq+rs8TCASEbdvr2kYrlJdr+eezbVsAEIFAoOR1LFeizmJTPz3TLMtCPp+H4zjI5/MAgAsXLmDfvn1qnRMnTgAA9u7di5GREViWhcePH2NhYaHh93nzzTcBAIlEArdv3wYA/PLLLw295vPPP2/8A1Xx2muv4fLly+vaRiuUl+uf/vQnAMDx48cBAJcvX4ZhGIjH4w1vk+VK1HpM/PRMyuVyGB4expEjRzA2Nga/36+a823bhsfjKVlf13V1/8iRI/jkk09Uk3Iul4NlWbh79y4cx4FlWbAsC47jqNd8//330HUdjx49QjAYhKZp8Pl8dWO8ceMGAGB8fFw9NzAwoG5+vx+O4yAajZY8Hh4eht/vV6/ZuXPnGkuptcrLdcuWLep+Pp/H5OQkLl68qJ7rdrn6/X4MDAwgGo0il8up+1KvlCtRy3W7yYGoHXRdF4ZhCCFWmqDdP/VqP3ufzyd8Pp8wTbOi2TkSiQifzycMwxCapql1ZbeAEEJomiZ0XVddBI00Eeu6XtHs7fP5RCQSUc3imUxGZLNZdV++zv3emUxGmKbZQKm0V7VyBaDKtDzGbper7OIRYqUbQtd1USgU1Gt6pVyJWo2Jn5458uAei8WEEE+TulR+gJfryOTiTjxutfqik8mkej/3fSmbzQpd1wUAEQqFSmKcmZkp2ZamaSKTyYhYLCYAqEQnk1UymVTbkEzT7IkEVa1c8Wv/fLVl0lrLNRaLCU3TSp5vplwzmYxK/JFIpGTbQvROuRK1Gpv6aU3y+XxJ8+nAwEBTfbftEI/H4ff71YjuBw8eIB6PI51O48CBA2q948eP48cff6x4vW3beP/996uOLq/FcRx888036rHs4/70009Vk3U6nYZhGLBtG9u3bwcAzM/PAyhtTs7lcigWi9i/fz+mp6cRi8UqRrD/8MMP+PDDD0ue+/nnn/H2228DAFKpFILBICzLqht3KpWq+P7WewlerXJdWFjApUuXKrpX6mmkXC9duoRPPvkEsVgMW7duBbD2crVtG8FgsOQ5d7kStYJlWQiHwxgdHW1o3dHRUUxNTanxSa3S0cRffqDpZbJfVd7cfX+9opsxygO8rusQKy1HFQfOTgqHw7hx4wauXbuGkZERRCIR/Pzzz2r5yy+/rO4fPXoUly5dKulL3rNnDwKBACYmJmq+x+DgIPbs2VPy3MOHD/G///0PPp8PS0tL8Hg8CIVC8Hg8ePjwoXrvRCKBy5cvY2JiAo7j4Pvvv4fP58OTJ0/Utn766SdVnouLixXl+fDhQ2zfvr0kgVqWBV3XVSIbHx/H2bNncfr06bonYuPj4xBCYGZmBsDK99jMCU811crV5/PBNM26l9WttVwPHDiAM2fO4LnnnoPX611zucbjcTUQUSovV6L1mpqawv79+wFADVR1C4fDJWN3vF4vEokE5ubmmh4Uu6pONi/IZjYANZtTe4lsni3vL+wl3YpRXloViUQ6+r7VVOuXl9z9uG62bYtIJNL232GhUBChUEj109d7v/IuCTf5Hbuby2OxWEWTtvt963VbSK3+HjtVrslkUn3+8ssuy9UqV9nUX96cX69cidZCjjOqtX/L3FhtudyXW5k3O5r45YeXg656nTxJWe264W7qRoyyH7UXTuBWS3ByAFm34sxms2qnrbXTC/E0CdU6F6+WoFYjP3stvfQ9NstdiahXLvXKVS6TgyaJ2mG1xJ3NZoVhGHX310gksuoxpBkdTfw+n68jtUS5o9e7rXYQdR8wag1K6rZuxShr0eUTpCSTSVWDlMtlzay87N3rue+XPy6v5ZaLxWLrnqjlWSVbGWqNhJff40Y5ESfaiOTxr9pxqlAoCMMwRDabrdtyKfflevtzMzrax59OpwEAr7zySlvfR/za51zvVq8vFwBu3rwJADAMo6lBSZ3UrRjltdhyYhZgpX/q9ddfBwAUCgUEAgHYto3PP/8cwWBQXXttGAaOHj0KYKVfOZPJ4Pbt2zAMAwCgaRr+9re/IZlMAliZuEVOn1vNjRs3MDw83PoP+QyQffb379+vulx+j8eOHetYTET9Rs4rUe049dlnn+HYsWOrjq9xL6+1PzejY4nfPcrYPbtXr7p+/ToAYGxsrOHXyFHljQyysyyrYrBjtVu90dlriXG98vm8mtFOjni2LAvnz5+HpmlIJBLweDxqQpznn38eAGCaJgCgWCzC7/dD13XcunVLzeEut3nlyhWMj4/jxRdfbDge94h9oHIQaT/dqqk20527zA8ePKier/e7rPdb7Pbn7ka5EjVCjsgvP06lUinYtq1mDpUcx6k7UHq1mSsb8X/r3kKD6tVOp6amcOfOHfV4aGgIH374IR4+fIjJyUkAK7OpBYNBta5pmjXPkhrZUU3TrFvrX0vrxN69e3Ho0KGKL7gar9cLsc7/R+pUC4qbnMLUMAw14vmLL74AsJLU3TPgASvTngIrZ6yGYWBhYQGappWMapXb1HVdjf6Wl2Vpmtb0H62st1z7gbwqw/09Amv/XbLMiZrzhz/8AcViEYlEouT5F154YdXZKderYzX+erXTbdu2IZ1O46233sJbb72F8+fP4+uvv8bWrVuxuLiIdDqNvXv3AoA6QZDX7Vaz3qb+tbZOyOvH3ZeOtUurW1Asy1JTmdZz4cIFAKXNw/KHm0wmK8r5xIkTanpUXdeh6zqKxWLJ357KbQYCAfXcp59+CgB4//33V419aWmpwU9JEpv5iTqr/Dj1+PHjkmOlaZrw+XwQQuDatWttjaUjid9xHFU73b17d8XyBw8eQNM0nDhxQk28sWnTJng8HtUvMjQ0hFQqhfn5efW3n+1y9epVACs1UPf7OI6DcDiMaDSq7ofDYbXc3RctWyc6HSOwMmGJvCY0HA6rZG5ZFoLBIILBIFKplEq+8Xgcb7zxBtLpNA4fPlxzsohcLgfbtgGUNg9LT548geM4iMfjGB4eRiqVgmVZ2LdvHxzHwZdffqmS+5kzZ+A4Tsk2gadlvLCwUDIeoJahoaGS1iKqtGvXrpLHuVyuajM/EbWebFFb7Ti1tLSExcXFVbdXvj+vybqHBzZAjuZHjcsR3HN1y3XkSG55CWChUFCjxNspmUyqSy9QdqlPKBQSoVBITREqR0XLdeSymZkZdb/TMQqx8neigUBArSdfg1+vqJBxu6c51TRNBAKBupc2yUtKykeBy2up4RqVn0wm1fpwjfB3rxcKhdQ67ud1XRemaTZ0pQJH9dcmr/ooHwVc63skotarN6q/fB15XCznvoKrFaP6uz5Xv7yWOBKJiEwmU/GhZOKvd6lDJ4VCIXXAlMnU/ecpcplMvt0QCATU3OTyulH5ZyjuuGViLZ/bvhY5WVArJzdZ7zblNbK8FrtSreuC2/E9ElF17uv413qc2tDX8VdTXmsuJxN/o//M1W6yNirE0z92EWLly5UnKEKs/CnIzMzMqsm0HQqFgvqhyPeXscmZ5GStvFAoqO+gUCjUbFFp9XWkQpSexa4ncbfypLA8HtkK1Uh8cgcvn0nO3eIlROlOLH835bf11sZlK075hCHt+B7bwV1msra00SYZIpJWm7mvHrkvY6PO3FfOPatZrUl9ZIH1Qu1E1owNw1DdE/KLKK/9y5OVTtdEZ2Zm1BSlcM3o5/P5hKZpwjRN1TQvy9w9kUutE5VWNw+XN/uvt3VEdsOsVyAQKGlui8ViDccWCASqdvHIZrxkMql24lAoJAqFgshkMuq+3A+q/QNfM+SkINW+y43SzC/LLJvNVpxUE21Ect9rZt+W+7KmaS2tRHY18csDX7Umfkn+fWYvkAcjGbdbeZzZbLYrNSoZRyaTKekjL4+5PP5634EQT5uH6/VDdVssFlt3zd8wDNVVIkT1v2utRo5pkCeqbpFIRJ0kGoZR9fcsf1vr/a1nMpmKsRrVWhV64US6nlAopFrT5IlpvfE98rtnqwD1Mnmi38iJdyaTUVP5tjqXdL2pfyORtblencKX1iebzYpAIFBSQzcMY9WdTg48LRQKVRN/+TTG1bRzMOhG5D7RbKTlTLa4EdHqOjpl70Z34sQJZDIZXjf+jPrpp5/w5ptvIhgMQtM0fPPNNygWiyUT3ESjUXWZpPTZZ5/Btm288MILasIpyXEcLCwswDRNhEIh2LZddQa8ubm5knkM+pm8xHNmZgaxWAzFYrHkb3blJZ9+v79kPgjDMBCNRrv699BEGwETfxNGRkbg9Xr5H93PqIsXL6qJoYLBIBKJREkytiwL586dw/T0tJrrIB6P4+zZs7BtG0IIRCIRAFBzJ9y6dQvAyqROf/3rXwE8nenQvV3btvHSSy+19wNuEP/6178AADt37lQTd33zzTdq+ddff418Po8jR47gwYMHAFamMS0Wi9i2bRsSiUTNuSiICGwbIxKi9LJRIZ6Ofnf3K8vBnTMzM6r5H7823xcKhZIR87FYTA3MkdupdVkPXH3v/c59pYccRyKb/WWZxWIx1fwvvwefz6cuCS6/qoKISg0IwUm2iRoxOzuLnTt34siRI/joo4/YpNwls7Oz2LRpEw4dOgTbtjE0NITNmzfjypUruHnzJq5fv45Tp05hfHy826ES9aSO/UkP0UZ38uRJ+Hw+2LZd978iqH0sy8KZM2eg6zo0TcOmTZuQy+VQLBaxY8cO3Lx5U/3BExFVxxo/UYPy+TwePXqELVu2cJxHF+VyOSwvL2PHjh3weDxwHAf37t2D1+uF4zh4+PDhqv9vTtTPmPiJiIj6CEf1ExER9REmfiIioj7CxE9ERNRHmPiJiIj6CBM/ERFRH2HiJyIi6iNM/ERERH2EiZ+IiKiPMPETERH1ESZ+IiKiPsLET0RE1EeY+ImIiPoIEz8REVEfYeInIiLqI/8f2jSUf135LCoAAAAASUVORK5CYII=\"\u003e\u003c/p\u003e\n\u003cp\u003eWhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({I}_{ph}\\)\u003c/span\u003e\u003c/span\u003e is the photocurrent, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({I}_{sat}\\)\u003c/span\u003e\u003c/span\u003e is the reverse saturation current of the diode, q is the electron charge, A is the diode ideality factor, K is the Boltzmann constant, T is the PV module temperature, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({R}_{s}\\)\u003c/span\u003e\u003c/span\u003e is the series resistor, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({R}_{sh}\\)\u003c/span\u003e\u003c/span\u003e is the shunt resistor, and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({N}_{s}\\)\u003c/span\u003e\u003c/span\u003e is the number of connected PV cells in series, the photocurrent (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({I}_{ph}\\)\u003c/span\u003e\u003c/span\u003e) which is produced by the PV module is influenced by the solar irradiance falling on the PV module plane and its influenced by the ambient temperature of the PV module as shown below\u003csup\u003e26,27\u003c/sup\u003e:\u003c/p\u003e\n\u003cp\u003e\u003cimg src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAksAAAA/CAYAAADwkYz8AAANvElEQVR4nO3dz2sj5R8H8He/eF1qtj2J9pCphyIS6YZd6Q/Bhe6EuhehbLJ7kIVdatODIEjKpi5elCUF8SJpAwp7UJtiRQ8mJl1wYRMFd4MkIOihk0MVT5mN4h8w30N9Hp+0yTSTzOyk7fsFA9NMMn0mTzvzmed55vkMWZZlgYiIiIja+p/fBSAiIiIaZAyWiIiIiGwwWCIiIiKywWCJiIiIyAaDJRetra1haGhILmtra34XiYiIiPo0xKfh3DU+Pg7DMBCNRpHNZv0uDhEREfWJLUsuMwwDAPDGG2/4XBIiIiJyA4MlF5XLZbl+4cIFH0tCREREbmGw5KIff/wRABAOhzEyMuJzaYiIiMgNDJZc9P333wMALl265HNJiIiIyC1P+V2Ak6RYLAIApqenfS4J9cs0Tfz6668AgJmZGZ9LQ0REfmLLkkvsxiutra0hEomgVqshn88jEolgdXX1yG1uW11dbZnaQCz5fN6T33cciekfRkdHMTs7i7ffftvvIrmmVqu1rX+v/t6IiLpRLpexvLyM8fFxmKZ55HuHhoYQiUQAANlsFuPj41hdXUW9XvesjJ4GS+1OzCf1cfrvvvsOAKBpWst4JTHXUrFYxM2bN/H7779jd3cXw8PDttu8sLW1BQAIBAKwLEsu8/Pz8j2nqc7aSSQSiEaj8mc3u1RN08TQ0BDGx8c7viefz7etg06Lk7m8QqGQrPNkMilfv3LlSl/HRUTUC9M0EYvFMDs7i8ePH2NnZ0dePzud82ZnZ1v2EYvFcPfuXWxtbUHTNO/mN7Q8lMvlLABySSaTXv463+RyOSsQCMjjLJVKcluj0ZDfgzj+arVqGYZhu81t1WpVli8ej9sey2moMzvxeFwe/+bmpmv73dzclPutVqtt35NKpSwAVjQatRqNRku9iX/XRqMh/95yuVxPZdE0zQJgaZrW8/EQEfWq0WhY4XC44zVJnPM0TbN0XT+0pFIpR/vrl6fBUqlUajnRNxoNL3/dwEomk4eCqG62eVEGuwu1ZR3vOhP/RP0S/3BHfVe97Ffsu1MQquu6FQgE5Pe+sbEhy6IeWzQatQD0FFirAfFpDIaJyH+6rtvesIlz1MGgyI56c7mxseFWUS3LsixPu+H++OMPuX6aH6ff2dkB0H6gsN02N4kuOE3TEAqFOr6PdQZUKhUA+92Vdt+VE7VaDZVKBR999BGA/+rjoKeffhp37tyR33utVpPbzp07J9cnJycRCAQQDAYdl+Xbb7+V6+yCI6InrVwuywei3n//fdf2GwqFEI/HAQC3bt06cvyTE54GSw8ePJDrp/VxetM0UalUoOu6o21uqtVqcmZxdTxOO6e9ztSB+ufPn5frkUikpd9cfV83MpkMdF3HzMwMwuEwDMNoCYSEbDaLN998U/786NEjuf7iiy/K9UQigcePHzsqg/o7gKMDZyIiL3z88ccA9m9IY7GYq/u+fPkyAKDZbGJ7e9u1/XoaLHU60Z82pVIJ6XTa8Ta3ZDIZuX5US8JprzMxsSgAXLx4Ua6/9dZbcj0cDmNiYqLrfZqmiWw2i+vXrwMAbt68CeC/qSbsPidauYDW4K1X+XwezWYTwNGBMxGR20zTlC3r3ZzTpqamHO1ffWDp66+/dlY4G54FS16c6HslHgd3sojHEvs1MjKCmZmZtt0ldtvc1G1LwiDVmV9+/vlnuS7+SU3TxHvvvQcAiMfjKBQKjron7927BwCYm5sDALz88ssAgC+//NL2c2KeJwA9d7kdxC44IvKTel5Tb0g7uX79urwunz17FpFI5MjpbsLhMICjb0id8CxY8uJE36tEItHyqHw3S6FQ8K28bnLSkjBIdeYXNVicmJhANpvF888/D8MwsLm5iXQ67Xgc14cffohYLCY/FwqFoGkaKpWK7bwgaiuXW4Eru+CIyE/qec1umpxkMglN07C4uIhGo4FSqQRN01AsFvHaa6/ZTmmjnqPdmnvJs2DJzRN9NpvF2toalpeXHY8V8YrTlqp+l145aUk4LnV2cPyQWIrFIorFYs/zEdXrdTm2CwBGR0dx9epVNJtNfPbZZ7Z965lMBmtra4cGFNbrdVQqFVy7dq3ldRG4ilandtRWrm7uwI7CLjgiGiQvvPBCx20ffPABdnd3kUgkZC/MJ598IrcvLy939Tv+/PPPvssJeBgsiTxpQP8n+gcPHmBlZQXr6+s4c+ZMv0VzhdOWqn6XXjlpSTgudVYoFNp+R7quQ9f1ttsSicSR+/3tt9/kejQaRS6Xkz//8ssvtp9dWlrCysrKoVanTz/9FJqmHXraUQSu6j//QeJJScB5v3077IIjouNMtMoD+wO4n2TjiWfB0sOHD+V6vyd6cVfe66Pcfo5Z8pPTloRBqjM/qAHR5ORkS9oau/FF9XodqVQKm5ubh7atr69jcXHx0OuhUAiBQKBjV1y9Xpd1B8DRgPJO2AVHRMedXQYEL3kSLLl9ohdz//TaNXRaxyw5aUkYtDrzg9qyNjU1hZGREXkXYze+KBgMIpFIHOqmy2azaDabWFhYaPs58X41SBXUVq6DKXR6oQbO7YI3IqInTR360a3d3V253k2vxTPPPOP4d7TzlCt7OUA9+Xea2NA0TTkHwsLCAkZGRpDNZrG3t4epqamWbgsx98/FixdRr9exvb2N4eFh+Tlqz0lLgp91ls/nZavO2NgYzp8/78vgcvXJCREsijmRAGB7exsLCwtYXl5GOp1GMBhEuVyW//AHv4O7d+9C1/WOx3L58mWsr6/jm2++ORRo/fDDD3LdjTspNXDuFLwREXmtm16LSCSCYrEITdNagiN1XKnddU0dO+ratcSNacAPEtOYA7DC4fCh7dVqVaamwL/TksfjcSuVSsmcV2oaB5HHKh6PW/F43NP8LyeF05QWftSZyOUTCASsZDIpc7L1mo+tn3Qn6jT56vT7aj43sajHoH7ParkNw7DNa6R+lziQVqZarbbkGmxXH06J/bmxLyKiXjUajbYpnFTquVGcb0ulkryOBAKBjqmoutl/L1wPltQkpGI5mKOlWq1ajUZDJg0Nh8PyYiG+JJErTT1wkSNGTTZK7an1cFR+M7/qTOShUy/gqVSq50TC/QRLag62g39XIq9eIBBom6dIfE4ttzjebpd2312nAM0pNaBzO18SEZFTnW4UhVKp1HKTrd44xuNx22uEV+c7TxPpHkV8YeoduWiREBcPceCBQODQ53ji78yrlgQ368wwDFlOTdOsXC7nalmfBBHwdUoGOQhE0t2DAR0RkR/UFvteexI6Eec7NSG5G3wNlg5GlmpXiHhN3NWLu331znsQTvzJZLKlW8XNLPW98rIlwe06E917YvtxC5jEsXbT1ekH9btnFxwRDQpx3tc0zbWgRty8enHt8zQ3nB11fgQx4PfOnTsA9mfuFK+JuWZeffVVAP/NMq1pmhxg66fp6Wk5GWIwGMSzzz7ra3mA1sG8IsWGG9yss3q9jkwmg1AohEKhgGQyCeDo+YwGjTjW6elp1Go1V7Ncu0Gd9FLkpCMi8tvnn38OTdNgGAZu377d9/5M05T5N+PxeEtCcle4Gno5oI7pEK0z+Lc1QkSZ6l2xaLERkWMgELCi0ajrTXhOqS0rg9DS5WVLgpt1JsYIqYPE7QbtDSp1ILuu6642+7qBXXBENKgajYa8jui63vM5qlQqySEd7caVumHIsvqYHroP4tHAVColX9N1veVRQPHIOYCWGZjF4+oH3++HTCaDpaWlQ484+iWbzeLq1auHXt/Y2Og70nazzur1Oh4+fIi9vT0A+1MGzM3NHbupIGq1GorFIsbGxmzToTwppmlidHT00OvhcBiPHj3yoURERPbK5TK++OIL7Ozs4KeffnJ0Hchms3j33XcRjUZx48YNz6ad8S1YEvnODMM41glbY7EYtra2EI/HkU6n/S6Op05KnRERETnhy5gldeyLXdbh40CMWXnllVd8Lom3TlKdEREROeFLy5I66/GgdF/0olar4aWXXgJw8ltbTkqdEREROeVbN9xJ0Gm8kkgL8txzz2F+fl6O11lYWMDw8HDHbcFgEKZp4t69e9jb28Pw8DDm5uZOdBBGREQ06DzJDXda3L9/H8D+4FnBNE3cvn0bOzs7MAwDyWQShmFga2sLr7/+esdtN27cALA/iLrZbGJxcRG3bt3C33//3TJQmoiIiJ4sBks9yufz2NraAgD89ddfME1TjuB/5513ZMK/K1eu4MyZM5icnMTZs2c7bhOfrVQqiEajSCQSmJqaci1jMhEREfWG3XA9qtVq+Oeff+TPExMTLY87Dg0NdXxc227b8vIy1tfXoWkavvrqK9+nRiAiIjrt2LLUI7sgRjw5dunSJUfbMpkM0uk0rl27htnZWaysrKBQKLhUYiIiIuoFgyUPiKfGpqenHW1bWlrC/fv3oWkaAODcuXMelpKIiIi6wWDJA2NjY0ilUrhw4ULX20zTRC6Xk7nRcrkc5ufnn0h5iYiIqDOOWSIiIiKy4csM3kRERETHBYMlIiIiIhsMloiIiIhsMFgiIiIissFgiYiIiMgGgyUiIiIiGwyWiIiIiGwwWCIiIiKywWCJiIiIyAaDJSIiIiIbDJaIiIiIbDBYIiIiIrLBYImIiIjIBoMlIiIiIhsMloiIiIhsMFgiIiIissFgiYiIiMgGgyUiIiIiGwyWiIiIiGz8H5P7PFQPBBigAAAAAElFTkSuQmCC\"\u003e\u003c/p\u003e\n\u003cp\u003eWhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({I}_{phs}\\)\u003c/span\u003e\u003c/span\u003e is the photocurrent at standard conditions (STC), \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({I}_{rrs}\\)\u003c/span\u003e\u003c/span\u003e is the irradiation at STC and is 1000 \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\frac{W}{{m}^{2}}\\)\u003c/span\u003e\u003c/span\u003e,\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({ K}_{i}\\)\u003c/span\u003e\u003c/span\u003e is the temperature coefficient, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\varDelta T=T-{T}_{s}\\)\u003c/span\u003e\u003c/span\u003e is the difference between the PV temperature (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(T\\)\u003c/span\u003e\u003c/span\u003e) and STD temperature (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({T}_{s}=25℃\\)\u003c/span\u003e\u003c/span\u003e), from equations (4) and (5) the output current (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({I}_{PV}\\)\u003c/span\u003e\u003c/span\u003e ) depends on the amount of sunlight and the PV temperature, the characteristics curves(P-V), and(I-V) for different irradiance and temperature are shown in Fig. \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003e,\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003ePV System Modeling\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe modeling process for dynamical systems is used to drive mathematical representation for the system, it depends on two techniques: mathematical representation is derived from physics or measured input-output data, in general, the modeling process can be organized into three techniques: white box, grey box, and black box modeling\u003csup\u003e28\u003c/sup\u003e, in white box modeling the physics of the system, are completely known; and the mathematical model can be constructed entirely from prior knowledge and physical insight\u003csup\u003e29\u003c/sup\u003e, in grey box modeling, some of the physical parameters are known, the other parameters are determined from measured input-output data, but the black box modeling is used when the physical properties of the model are unknown, and the model is constructed from measured input-output data, but the chosen model structure belongs to the candidate models that are known to have good flexibility and have been successful in the past \u003csup\u003e30\u003c/sup\u003e. The back box modeling is chosen for the system identification methodology in this paper, in the back box modeling method, many models are compared with the real data, and the best is determined by an optimization criterion. Four basic stages characterize the system identification technique (I) The measured input-output data, (II) the proposed model, (III) a criterion of fit and parameter estimation, and (IV) Validate the model which is obtained.\u003c/p\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e4.1. The measured input-output data.\u003c/strong\u003e The data used in system identification is obtained from the PV system, shown in figure .1, where the input data is the PV current and output data is the PV system output power.\u003c/p\u003e\n\u003cp\u003e\u003cspan\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e4.2. The proposed model structure\u003c/strong\u003e. a fractional order model is considered in this paper, A generalized fractional-order dynamic system can be represented by the next Eq.\u0026nbsp;3\u003csup\u003e1\u003c/sup\u003e:\u003c/p\u003e\n\u003cp\u003e\u003c/p\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003e\u003cimg src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAlQAAABCCAYAAACRiEiGAAAes0lEQVR4nO2d32sUWfr/3/Xlc6tJdXs1zAZJ5WKGcciglVHEGVAYqzcuyw7E7VZvBhzG7WbYi8GNpuOwDP4YulcRZEi6FwMirt3ZHWFZtrMmQgS7V9ykDd2oKGyqCa541ZU2zh9wvheZc6zqqurf6XTi84Im6arqqlN13nXqOc95zlMSY4yBIAiCIAiCaJj/t94FIAiCIAiC2OiQQUUQBEEQBNEkZFARBEEQBEE0CRlUBEEQBEEQTUIGFUEQBEEQRJO8dQZVoVBANBpFPB4HAGQyGSSTybaWIRqNwufz1bRdIBBAPp9vQ6k2H8lkEtFoVFy/aDSKQqGwzqWqDuljY5HJZBCNRjE1NQVgVXeZTKZtx69FL6QVolVMTU3B5/PVrKVMJoOBgQFxf2xq2AZE13UWDAaZoigMgPjIssz8fj9LJBKOv0skEkyWZRYOh5nf72exWIwpisLS6XTNx45EIkxVVXFMTdNs2wSDQaaqqm15sVhkqqqycDhc8/HS6TSTZdn1nDY7iUSCaZpmqWcA4jrqum77Db/OmqaxSCTCFEUR+1lrcrkc8/v9TJZlUdZyfaXTaaYoCkulUo7lJn20n3Q6bas3AExRFBYMBl3bCN4O8XYhlUoxWZZZsVisesxisWhrxyKRiG0bN03UqxfSCtEswWCQaZrmqu9EIsH8fj8DYGnfdF1nqqqyYDDYrqKuCxvOoMrlckzTNNHApdNpBoDlcjmWy+VYMBgUho650nVdF9vx39XTuBSLRaZpGkulUmK/qqqyWCxm2w4A8/v9tn1omuba+PEG2e2czWWvB6cH+kYhGAyyWCwmjCZN05jf72fFYpElEgmmKIpjHQaDQcv11zStYiPQKhKJBAsGg6Ke+MO1nHA47FiflfRRLBaZLMs2I4yx5vTBjc5OIBKJtMXoLae83iKRiKi3VColDPryuuEa5Lrixnst9VDejum6zmRZtv02lUoxALZ2hrHKenFb14xWiLebSs+oYrHI/H6/6Lw6tbXFYpEpiuKo5c3ChjOoZFm2PAC4AWUmkUjYGsBwOGxprFOplKs4nHD6PQCbh4QbeOWi4d4xJ6HlcjnhOav1+LWyUQ2qWCxmqVduEJvrnvfQZVkW9cANWvM5c89Bq0in05aHYaXjOvXINE2zGVqV9MHYG527rW9UH2+7QcUNDHO9KYpiKwe//mYdqapquXbhcNhyD/NOmFMHz+m3iqLYyheJRBwNoEp64e2fm+Yb1Qrx9sLbX7dniaZpTFXVqp3Wejy4G5G2xlAZhoHR0VF4PB5IkoRAIIB4PI5oNFrzPhRFQU9Pj/g+MzMDTdMs2wQCAWiahvHxcbHs4cOHOHDggCjHH//4R3i9XvE9mUzC4/GIWANJkkScFQD09PSgu7tbfL9+/ToURUFvb6/l2Pfv3wcA7Nmzx7L80qVLCAQC4pjmsn700UcolUq4cOECJEmCJEm28z5+/Dimp6fbFgMRj8fR19cHSZIwMDCAZDKJQCDQlmMDQFdXFzRNg2EYAIC5uTkAwN69e8U2Xq8Xly9fRqlUwp07dwAAT58+BQDs27cPwGo8SzabxdatW8W+RkdHIUkSkskkBgYGIElSzdd1amoKr1+/Fvsyl0VVVbx+/RoAkM/nkc1m8emnn9r2MT09jYMHD1qWuekjn89DkiSh5W3btkGSJNs90259dBLNtCtbtmyBoijYsmULgNUYS13XRVvBOXv2LGRZxvXr18WybDYr9JjP5zE+Po6uri6hjYmJCXR3d4Mxhu7ubkxMTIjfDgwMWPY/Pj5u0wQAzM7OQpZl9Pf3W5a76cXj8eDIkSMAgEOHDkGSJFt81dusFaIxJiYmoCiKaFfNJJNJTE9P4/LlyzY9ljM4OAiPx4Mff/xxrYq6vrTTelNV1RKzxHt9jY7pc6+AUw+b9+zM37l3iw8l8fiIXC7HUqmUxWUfiUQq9tzhMATAmLP3gVv3bj1G3qN0igcyw2M16gENeKj4kBSvF16+9Rz/5mVywqwBPjTm9/tF7IssyywYDIr64sM43Gtg9jZx3ZR/yq97uYeqHL6f8p6Ykwezmj4YY0KrlWhEH5vBQ9XKdoV7np3qlmvG/J3HMMViMeb3+5mmaSwYDDrqyu3cuCac6h8O4QPV9MLjsirRiFaItxdFUVxHUBRFqWu0h8dhbUba5qGKRqPIZrO4deuWsHKPHj0KAPj4448b2if3RuzYscO2zuzJAFZ7Zd9//z16enowNjaGoaEhfPXVVzhx4gT6+/sxODiIpaUlnDt3Dv39/VhaWnLcLwAxg+fDDz+0rZubm7Odz8uXLwEAW7duddzfo0ePHL1d5fT19WF2drbiNs2SyWRw4cIFxGIx4ZHif528Le3i4cOHUFXVcZ3ZQ+n1evGPf/wDO3fuxMmTJ7Fv3z7cuHED/f39OH/+PIDVXtLi4iLGxsbg9XqxuLiI999/HwAwPDwMtjoUbvkMDw/XVd7Z2VkoimLrsT158gQA8MEHH4hl1fRhGAZ0Xa96/duhj06j1e3K48ePAUDowUy512psbAyHDx/GL3/5S5w4cQIjIyM4cOAAzp49K+qdeywBWDzcZrhX+7333rMs5x6knTt3WpZX08v8/Lyjt8tMLVrJZDLCY17tU8usVGLjous6urq6bMvz+Tx0XcfAwABCoZDwEvf19bnO6tu+fTump6fXusjrQ7ssN0VRbD0tc/BnI3AvgJNnh8cY1FtGvi8+zuvUU+XxPeXrnOJ8zNu7wYOtq1HeQ3ZaDwfvitOn0qyl8utmDvxfL1DBQ1ZvoKOu66IHz4ODc7lcVQ+hmWpxWbIsu84ALddCNX3Uev2r6cPN++b0cfNecA+gU4yROZaNsTfeHidPET+nWj6VzqnV7Yqmaa7enfLJDtXgs4kZY2JWsRN8VlQ53DNcfq9W0wtcgtjNVNMKQXCc4gw5/D42t4c8QL28PeCUjx5tJtrioeJW7P79+y3LFxYWqvakKrGwsODq2Xn06FFdMT88PxHfl8fjgc/nc+yprqysOO6Dx/mUe7bctudMT0/beqGNcPv2bZtnBQDS6bRtudNYOLA6Hl5eJ0+ePHGM42gVhmFY4tXK4T11Jw8N99589tlnNR9vbm7Oco5zc3OIx+NVPYTAm157NpvFoUOHXHvmpVLJcfnMzIzN01ZNH9yD0ez1d/K+aZqGSCTStEeuXvbt22c7ZiQSgaZptuW3b9923Ee97QrPS1YpR1S5NszMz8/jN7/5Tc3n+M0330DXdUiShFevXmFoaMhxu1evXjkuv3fvHgC7t6ySXvi5mT2gBNEMP/30U9VtLl++jMHBQQCrowQjIyMolUrimfi20BaDileI+SYvFAqYnJxsypCYmZlxbfwmJyeF678Went7sbi4KL4vLi5ifn6+apCdmUuXLgGwu+6dXKUc3gCWD1GuF6VSyVYnf/rTnxoelq1GMpnE7t278bvf/c51mwcPHgBwHsK5c+cONE2ryRjiBAIBjI2NAVg1UpaXl8X3apQbAm4PeyempqaEe9xMJX0Aq8OH5RMv1guv14vl5WXbeY+NjWF5edlSD4ODg2CMrdlkhnraFZ/Ph2vXrmFlZQW//vWvHZP55vN5lEolR8OdB6vXY7h7vV7Mz88LndTTlhQKBSSTScdh40p64ca3W4epHmjIjwAgJmzUA+/8PX/+vNXF6WjWLVP6xYsXAdgNiWQyKcZeeQZip9koU1NTKJVKOHnypG1dNBqF3+9vSaPiBH+43bx5E8Bq4xcIBFAqlRw9ZrzBd+oZ85gaJ09YJxCPx7G8vGyLH8lkMsKrxLPP15sJN5lMYuvWrbhy5UrF7a5evYpgMGi7roZh4MyZMzUbQ+3E7/djbm5OaDcej4vzLH9gV9IHsOo12bVr1xqWdvPg1K5MTU1hbm4Ot2/fxvnz53H69GmcOXPG9tu//vWvUBTF0QAMhUK4ceNGXUZRrXz++ecAIIy8TCaDcDiMUqnk2GGspJeFhQXXWMN6cfIiun3q6VgQGwtuHHFj3cw777wDAPjXv/5lWc7bPbc45M1KWwwqftF5hYRCIWzfvl2sj8fjSCaTGB0dxb1793Do0CEkk0ncvHkTs7OzOHXqlG2fV65cQSwWszxkM5kMAoEAVlZWRADyWtDf349wOIzx8XFIkoRwOIyRkREsLy87Nmb8/F+8eOG6z6dPnwrDzO31KHNzczbDptUoiiKCVZPJpKWHkclkMDo6KupmZGQE8Xgc4XAYS0tLOHToUF3HCgQCGBwcdA2u5ccslUo4e/asWMZ77z6fD7du3arLO9UuRkZG4PF48NFHH6Gvrw8rKyviwVluPNeij4cPHwKAuE+caIc+Oola25XHjx9bvJt79+6FruuWfRmGgfHxcdy6dcuyjL9m44svvhBDGq1maGgImqbhyJEj8Hg8uHnzJr7++msAzsO81fSi6zoMwxDtoRNvm1aI5lAUBUtLS7blvb294lnIDfxCoYAvv/wSqqo63jMLCwsd43FvOWsYn2WBB6LxV27wrKmyLIsAymKxKALg+LJIJGIJBE2lUrZXzsD02pn1SmLJy+0WzOuW5JEnp8TPgX1uU70bzXCMOtMm8MRr5gSqPGg2HA6zYrEoUgHwa87Ym8zznEQiIVJPlH/Kz5EHNpbDj1v+4UkRN1pyOLfgY8bc9cHYm+Bknv7BKdCzUX1s9LQJtbQr5fst1xsPYi/XGU9VUc9khVbhNvGF46YXniQYP6dbcPo9ZUsn6sUt8Sxjq8+wcDgsdMfbKbf22XxvbjY6LtQ+FotZclpUmh3TSfAG0K2RqpYJuxqdlruDN8r8fHgeHk4rDKrNhtkALWez6aOTqGZQdSJ8lpQbzeiFtELUS7VM6bXSbDvX6XRcq+L3+y0JxLiRsh69RCfS6bRjin3+Mt5KVHr3ViU6sUdZbvhqmmZ57149bIQHXD0Eg0Gb16dSwkjOZtJHJ8E9fGbjv55EhGtJeUeEscoJi800ohfSCtEold7lVwv0Lr91wPwCWH7zh8PhjnlLNXepq6oqjIdYLOb4YtNy+PBePUMsnfqG+HLDl7t563kJsa7r4ibl9dxp59kIPNcUbzhqfdP6ZtJHJ8GvKx/edHvJ9HpgfgsBH04PBoM1vRetXr2QVohm4bnY6vUw1doGbnQ6yqAqFou22JhYLNZxDQBPBlktrsWNWmNFePxYJ/YmY7GYpVypVIrFYrG6bjRuUFUaDtyomONyKsXGuf12o+uj0ygWiywWi7FIJNJx18scF8pfCVPPfVSLXkgrRKvgr+6qVUt8VKdTOjFricTYz9kfCYIgCIIgiIZYtzxUBEEQBEEQmwUyqAiCIAiCIJqEDCqCIAiCIIgmIYOKIAiCIAiiScigIgiCIAiCaBIyqAiCIAiCIJqEDCqCIAiCIIgmIYOKIAiCIAiiScigIgiCIAiCaBIyqAiCIAiCIJqEDCqCIAiCIIgmIYOKIAiCIAiiScigIgiCIAiCaBIyqAiCIAiCIJqEDCqCIAiCIIgmIYOKIAiCIAiiSf5vvQuwEZiamsLjx4/R1dWFoaEheL3e9S4SsUYYhoE7d+7g+fPn6OnpQSAQWO8iEZsM0hhBbE7IQ1WFUCiE3//+9wCAq1ev4tixY+tcImKtMAwDPp8P165dA7Ba99FodJ1LRWwmSGMEsXkhD1UFkskkkskk/vvf/8Lr9WLv3r345JNP1rtYxBrx7bffwuv14vbt22LZ7OwshoeH17FUxGaCNEYQm5eO8lBFo1FIkgRJkpDP59e7OPj73/+OYDAohvhu3rwJv9/ftuNHo1H4fL6q2wQCgXW5Xj6fD5IkwePxtP3Ya0EymRTeSMMw8Le//Q2ff/75OpeqOp2uE+IN7dRYLbrg25E2CKIFsA4jHA4zWZZbtr9iscgikQhTVZUBsHw0TWOxWIwVi0XH38qyzNLpNEskEkxVVeb3+123bSXFYpGpqsrC4XBN26fTaSbLMkskEmtcMjv8unQKuq6zYDDIFEWx1LUsy8zv97teo1wux/jtEA6HmaIoLBKJtLPodbORdLKZ6HSN1asLxkgbBNEKOspDxTl48GBL9mMYBo4dO4YdO3bg9u3bKBaLkGUZsVgMuq7jwIEDGBkZwe7du229s3w+j1KphH379uH58+cAgJGRkbYEpB87dgwHDx7E+fPnbet8Ph9GR0cty/bt24e7d+/iyJEjDfUyM5kMJElquLz79+9v+LetJJ/PIxQK4ejRo1hcXEQ6nQYA5HI53L17Fx6PB0eOHIHP54NhGJbfPnjwAJqmoVAoAAD6+vpw/Pjxtp9DPbjpxDAMeDweTE1NWZY3q5NaPR7toFnNNspG0Fil9oNft/KyNasNgiDQeR4qVVVZLBZryb78fj/TNE18TyQSDABLp9NiWS6XY7IsM1VVLb+NxWIWz0sqlWqp58yNRCLBZFl29ITx8qdSKcffhsNhy/nWSjqdZo1IoVgsMgAsl8vV/du1QJZlS48/GAzazotfw/Leu9/vt+guHA6vieeNezea9UxU0gk/bzdvaqM6iUQiDf1uLWhUs82yVhorFotM0zThOW/UE15JF9xzVd7WmWlUGwRBdKCHKpvNYs+ePYhGo/B4PI497VrxeDzYtWuX+H7v3j0Aq70xTn9/P06fPo1sNmvpmeXzeezcuVN8f/36tfg/mUzC4/EgHo8jEAhAkiQkk0nxfzweF9sahoFkMom+vj7RMwVWe4oejweZTMZS5kuXLiEQCNg8YbznCwCHDh2CJEk2b8Hx48cxPT3dth7m06dPIcsy3n33XYRCIUiSZDvPdqIoCnp6esT3mZkZaJpm2SYQCEDTNIyPj1uWZ7NZfPDBB+J7qVQS/4+OjkKSJExNTWFgYEBocmBgoK54v0KhgDt37riui0ajGBgYsCyPx+OOnhgnneTzeUiSJM5t27ZtkCTJNous3TrpJAzDwOjoKDweDyRJQiAQQDwer3mm3VppbGJiAt3d3WCMobu7GxMTE2JdPp9HIBCweKYNw0AoFLK1AW7tRzwex7Zt25DNZpHNZkWsann78zZrgyCaZr0tOjN8HD8cDouxfFVVmaIoLdm/qqqOvS/e2zV7rhRFETEO4XCYARBlSqVSLBKJMEVRWC6Xs/1v7sEmEgnxe76cH6/8mLquV/RA8biNSjQSl9Fobz8SiQgvYC6XEx6rYDBY975aDS+L07WIRCKW8+XXXdM0cU6yLAvPWyqVYpqmCc9B+f+8Ds31av6Ua65cI3yZ3++3aIKXs7xuqulEUZSqddCITjaDh4q3J/wacw9TI7FDrdSYWUfl1zkYDIoYUI7Zm1V+DDdd8PXVznUjxA8SRCfSUR6q+/fvA4Al2d3hw4eh63pL9p/NZi0eK47ZYwWs9v50XceVK1cAAF1dXcjlcqJMg4ODWFpawrlz59Df34+lpSX84Q9/EP/v2LFD7CsQCIhYhtnZWXG8WCxmK8fLly8BAFu3bnUs//z8fNX4sr6+PnGctWZ2dhYzMzOIRCLo7++H1+u1xIisJ0+fPgUAS11w9u7da/n+7NkzqKqKL774AsBqTFg2m0V/fz+A1fpeXFzE2NgYvF4vFhcXEYlExP/vv/8+gNV6ZYzZPuYp8m4MDw/j66+/BvDmPhgeHnacVVpJJ1y7n376acXjtVMnnUI0GkU2m8WtW7fEPX/06FEAwMcff1z3/lqpMcDqBe/u7hb/j42N4fDhwwAgPEp/+ctfoCiK5RjV2o9nz54BqH6ub6M2CKIVdFQeqtnZWSiKghMnTrR837wh+vDDD23ruAGwZcsWAMB//vMfKIqCwcFBDA4OOu5vZmYGJ0+eFP/zMieTSZw9exaZTMZmqJkxu/45T548AWA38DjZbBZffvml6z5rxS2Y12k5Y8x1P9PT0wiHw5aHQqfAjZL33nvPtu7FixeQZVl8//e//42BgQHXjNVcH729vSgUClheXkZ/fz/y+Tz6+vrw4sULrKysoLe3t6ayrays1HweO3fuxOTkpGVZJZ3whzw38holk8m45lwr14mmaa5Go8/nw9zcnMjlBqzeI0eOHEEqlRL3V6FQgKqqCAQCGBsbq3rMSsvdNPvnP/8Zfr/fotf79+9DluWa685MKzW2a9cuXL9+HYODg1hYWKg60cPr9aKvr8+yrFr78fjx44bPlSCI6nSUh2pubs5mMHAjq1l4Y+PUO3v27BkURREN7ePHj6Gqquu+3B6wwGqsk8/na+iBVulByw1CJ0OsXso9KHymkpN3xQ0eY/Hb3/7Wsnx6erotDXYmk7HFf5hZWFiAoiiOZXn06JHlwfbw4cOKRuHc3JzwDJr/59/j8XjN5yxJEi5cuIALFy40PEutkk74Q75ZI9fJ2xaJRKBpWkMeuGZpVrP5fB66rtsMlYWFBUevbyaTQTQaRTKZdC1TKzX2zTffQNd1SJKEV69eYWhoyHVbN6oZ6rOzsw154giCqI2OMah4mgKzwWAYBqanp1uSTPPu3buujd8///lPfPXVV+J7T0+PGH5xore3F4uLi+L/5eVlsW5xcRHz8/MNpVfo6upyXccflJW8Xu3kwYMHAKwPbj554Fe/+tWaHbdQKCAUCuGTTz4R18SJmZkZ1+HRyclJMdQDAAcOHMBnn33mui+z1yQQCIiHbH9/P5aXlx09Km7UarBWopJOZmdnbUHS68nt27exvLxsuR8CgQAYYxbvL7+P6rmW9fDTTz8BsHZICoUCJicnLZNPgNWhQT40d+bMGYRCIcd9tlJjXq8X8/PzwkBtdfsBrHZ2Dhw4UPd+CYKojY4xqPhQhdlg4DNdeK4WPhPKMAwYhoF4PC5m1PF1TjMCDcPA5OQkzp07Z1uXyWQwPz9vefVDIBBYE8OF534xDEN4zMzwxt7J87KwsFDRa9Zu8vm87cF95coVMVQKrA7t8PrgPX7u2Uomk6Iu6+HixYs4e/ZsxWsxNTWFUqkkhmTNRKNR+P1+S/0ODw931DAI9zQUCgUsLS3Z1lfSydzcnGOcIGHn4sWLAKzxToZh4NSpU7h16xaGh4dx7do1jI+P2+IC10tjPM4qn8/b7p1KuuD3nVO8F0EQraFjYqh4SoOpqSkMDg4iHo/j1KlTiMVi6O3tRT6fx/fffw9d17GysgJd17Fz506cOnUKKysrWFhYAACcOnXK1vOfmJiA3++3uODz+Tymp6cxOzu75kMWvDHj05UBiGFMcyDqO++8A2A1/sIJXddhGAaePn2KH374wXE4Ym5uDqdPn25p+Z2YmZnB8vIyCoUCurq68O2332J6ehq5XA7AaqqBUqmE8fFxJBIJ3Lt3D4VCAbOzs+jt7YUsy2JaeT3vMeMejEo9+CtXrgjdcDKZDH744QcoiuKY8LAT4EY2HxIE3ujEMAxxztV08vDhQwCrRuujR48cz7ddOukU+DW7f/8+9u3bh1AohO3bt4v18XgcXV1dePfddwG88bxyo+jly5cWPbVbY7x9O3ToEABY4rM41XQBrIYz7N69GxMTE5bJP2aqaSMUCmFmZkZ46QmC+Jm1m0BYH6qqile8AGCKolim//JEdX6/nymKIr7DNE2fp13g6Lru+MoZ/DzduF2vWeBTnHnSUl3Xma7rTJZl27RrVVUdp7zzBKQAmN/vt6RbMG+DBhJt1jsFnV/XWCwmXr/BUydwisWiKA9PZhiJRGxTxfk6XddFOgGnj67rljLw6edmUqmU7XUgML0SxOmadRK8vPz1JcVikaVSKccUG2464UklZVlmwWDQdt0Ya1wnGz1tAk9lwNuWYrHIFEVhsiwLHTrt13zt10NjvEy8Tnm7yNOxmMvrpgvG3qSI4GkRnJJ/1qINnn6BIAgrEmMNBnKsEx6PBzdu3MDg4CAKhQIURUGxWITX60U8Hkc+n1+zOIx2kEwmEQqFLLOiaiUUCqFQKLQlSLgW4vE4rl69ivn5eQCr5ZNlWfTePR4PstmsCO7/8ccfXfc1NDRk8Qb4fD4cOHCgLu/WZmIz6aST4LMbzc2iJElIp9MdE79YiWZ0AdSmDZ4E9W299wjClXU26OqC9554zyqRSFh6zNzL5dQr30homlbXi00Za9zrsJb4/X7LeZgTKvKePmOsofpy8lC9bWwWnXQS5ckxy9ucjUAjumCsNm0kEomOSNxLEJ3IhjKoYrGY5T1U5Q9s/Dz818y7sDoB/s6tWg2GTn1TvCzLrg8mnl2eD2/VijnzPL9GG92AbpTNopNOg7+RIBKJMFVVGzJO1pN6dcEYaYMgWkHHBKXXwi9+8Qt899134vv+/fuxZ88e8T2VSuF///sfzp4925C7u1PgU6ij0Sh8Pl9V9/vCwgLu3r3bUQk2DcPA6dOnsXv3brEsFouJeuEzN/fu3Vv3UEpXVxcikUjrCrtB2Qw66UTGxsaQTCbx/PlzfPfdd67JfTuVenQBkDYIolVsuBgqgiAIgiCITqNj8lARBEEQBEFsVMigIgiCIAiCaBIyqAiCIAiCIJqEDCqCIAiCIIgmIYOKIAiCIAiiSf4/jFoRf0EQjg8AAAAASUVORK5CYII=\"\u003e\u003c/p\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003eThe dynamic system in (6) can be written in the next form:\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equb\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equb\" name=\"EquationSource\"\u003e$$\\sum _{k=0}^{n}{a}_{k}{\\mathcal{D}}^{{\\alpha }_{k}}\\mathcal{Y}\\left(t\\right)= \\sum _{k=0}^{m}{b}_{k}{\\mathcal{D}}^{{\\beta }_{k}}\\mathfrak{u}\\left(t\\right) \\left(7\\right)$$\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003eWhere( \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({a}_{k}\\)\u003c/span\u003e\u003c/span\u003e,\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({ b}_{k})\\in \\mathcal{ }\\mathcal{R}\\)\u003c/span\u003e\u003c/span\u003e, and (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\alpha }_{k},{\\beta }_{k}\\)\u003c/span\u003e\u003c/span\u003e) \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\in {\\mathcal{R}}_{+}^{2}\\)\u003c/span\u003e\u003c/span\u003e. in (6) if all the derivation orders are integer multiples of a base order q such that \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\alpha }_{k}\\)\u003c/span\u003e\u003c/span\u003e, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\beta }_{k}=qk\\)\u003c/span\u003e\u003c/span\u003e, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(q\\in {\\mathcal{R}}_{ }^{+}\\)\u003c/span\u003e\u003c/span\u003e, the system is said to be of commensurate order. If \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(q=\\frac{1}{\\mathcal{g}}\\)\u003c/span\u003e\u003c/span\u003e, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\mathcal{g}\\in {\\mathbb{Z}}_{+}\\)\u003c/span\u003e\u003c/span\u003e, the system is said to be of rational order.\u003c/p\u003e\n\u003cp\u003eThe Laplace transform of (7) can be written in the form:\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equc\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equc\" name=\"EquationSource\"\u003e$$\\sum _{k=0}^{n}{a}_{k}{s}^{{\\alpha }_{k}}Y\\left(s\\right)= \\sum _{k=0}^{m}{b}_{k}{s}^{{\\beta }_{k}}U\\left(s\\right), \\left(8\\right)$$\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003eThen the fractional order transfer function can be written in the form:\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equd\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equd\" name=\"EquationSource\"\u003e$$G\\left(s\\right)=\\frac{Y\\left(s\\right)}{U\\left(s\\right)}=\\sum _{k=0}^{m}{b}_{k}{s}^{{\\beta }_{k}}/\\sum _{k=0}^{n}{a}_{k}{s}^{{\\alpha }_{k}}, \\left(9\\right)$$\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003eBy considering the system is commensurate order \u003cem\u003eq\u003c/em\u003e, We can write \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\text{ɦ}={s}^{q}\\)\u003c/span\u003e\u003c/span\u003e, then the continuous time rational transfer function is:\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Eque\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Eque\" name=\"EquationSource\"\u003e$$H\\left(\\lambda \\right)= \\frac{\\sum _{k=0}^{m}{b}_{k}{ɦ}^{k}}{\\sum _{k=0}^{n}{a}_{k}{ɦ}^{k}} . \\left(10\\right)$$\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003eIn fractional calculus the differential operator has many definitions, The three most frequently used definitions for the general fractional differ integral are the Grunwald-Letnikov (GL) definition, the Riemann-Liouville (RL), and the Caputo definitions\u003csup\u003e32\u003c/sup\u003e, Gr\u0026uuml;nwald-Letnikov definition is considered in this work for the drive numerical solutions to fractional-order differential equations, which is given by the next equation:\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equf\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equf\" name=\"EquationSource\"\u003e$${}_{\\text{a}}{}^{\\text{ }}{\\text{D}}_{\\text{t}}^{ \\alpha }=\\underset{h\\to 0}{\\text{lim}}\\frac{1}{{h}^{\\alpha }} \\sum _{j}^{\\left[\\frac{t-a}{h}\\right]}{\\left(-1\\right)}^{k}\\left(\\begin{array}{c}\\alpha \\\\ k\\end{array}\\right)f\\left(t-kh\\right), \\left(11\\right)$$\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003eWhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(h\\)\u003c/span\u003e\u003c/span\u003e is the step size, and the binomial coefficient can be expressed with Gamma function \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\varGamma \\left(.\\right)\\)\u003c/span\u003e\u003c/span\u003ein the following expression:\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equ2\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equ2\" name=\"EquationSource\"\u003e$$\\left(\\begin{array}{c}\\alpha \\\\ k\\end{array}\\right)=\\frac{\\alpha !}{k!\\left(\\propto -k\\right)!}=\\frac{\\varGamma \\left(\\alpha +1\\right)}{\\varGamma \\left(k+1\\right)\\varGamma \\left(\\alpha -k+1\\right)} ,$$\u003c/div\u003e\n \u003cdiv class=\"EquationNumber\"\u003e12\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003eLaplace transform of the Gr\u0026uuml;nwald-Letnikov fractional operator is:\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equ3\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equ3\" name=\"EquationSource\"\u003e$$\\mathcal{L}\\left[{\\mathcal{D}}^{\\alpha }f\\left(t\\right)\\right]={S}^{\\alpha }F\\left(S\\right)$$\u003c/div\u003e\n \u003cdiv class=\"EquationNumber\"\u003e13\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003e.\u003c/p\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003eOustaloup\u0026rsquo;s recursive filter is an approximation for fractional differentiator and integrator, in a specified frequency range (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\omega }_{b},{\\omega }_{h}\\)\u003c/span\u003e\u003c/span\u003e) and is capable of producing a very good fitting to the fractional-order elements \u003csup\u003e33,34\u003c/sup\u003e, the algorithm for Oustaloup\u0026rsquo;s recursive filter is shown in the following form:\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equg\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equg\" name=\"EquationSource\"\u003e$${S}^{\\alpha }\\approx \\prod _{K=1}^{N}\\frac{S+{\\omega }_{k}^{\\text{\u0026apos;}}}{S+{\\omega }_{k}} \\left(14\\right)$$\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003eWhere\u003c/p\u003e\n\u003cp\u003e\u003cspan class=\"InlineEquation\"\u003e\u0026nbsp;\u003cspan class=\"mathinline\"\u003e\\({\\omega }_{k}^{\\text{\u0026apos;}}= {\\omega }_{b}* {\\omega }_{u}^{\\frac{2k-1-\\alpha }{N}}\\)\u003c/span\u003e\u0026nbsp;\u003c/span\u003e, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\omega }_{k}= {\\omega }_{b}* {\\omega }_{u}^{\\frac{2k-1+\\alpha }{N}}\\)\u003c/span\u003e\u003c/span\u003e, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(K={\\omega }_{h}^{\\alpha }\\)\u003c/span\u003e\u003c/span\u003e, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\omega }_{u}= \\sqrt{\\frac{{\\omega }_{h}}{{\\omega }_{b}}}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003eAnd \u003cem\u003eN\u003c/em\u003e is the approximation order in the frequency range (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\omega }_{b},{\\omega }_{h}\\)\u003c/span\u003e\u003c/span\u003e)\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e4.3. Criterion of fitting and parameters estimation\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe system identification method depends on the error minimization between model output and system output, the least-square method is used in this work, and the problem is expressed as follows:\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equ4\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equ4\" name=\"EquationSource\"\u003e$$\\mathcal{F}\\left(x\\right)=\\sum _{p=1}^{n}{℮}_{p}$$\u003c/div\u003e\n \u003cdiv class=\"EquationNumber\"\u003e15\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003eWhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({℮}_{p}={y}_{p}-{\\widehat{y}}_{p}\\)\u003c/span\u003e\u003c/span\u003e and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\widehat{y}}_{p}=\\widehat{\\mathcal{H}}\\left({\\mu }_{p},\\theta \\right)\\)\u003c/span\u003e\u003c/span\u003e is the output of the estimated model\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\widehat{\\mathcal{H}}\\)\u003c/span\u003e\u003c/span\u003e for the input signal (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\mu }_{p}\\)\u003c/span\u003e\u003c/span\u003e), \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(p\\)\u003c/span\u003e\u003c/span\u003e is the iteration, and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\theta\\)\u003c/span\u003e\u003c/span\u003e is the estimated parameters of the model \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\widehat{\\mathcal{H}}\\)\u003c/span\u003e\u003c/span\u003e, from Eq.\u0026nbsp;(6) the model parameter\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\theta\\)\u003c/span\u003e\u003c/span\u003e is formed by:\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({a}_{r}=\\left[{a}_{0} {a}_{1}\\dots { a}_{m}\\right]\\)\u003c/span\u003e\u003c/span\u003e, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\alpha }_{r}=\\left[{\\alpha }_{0} {\\alpha }_{1}\\dots { \\alpha }_{m}\\right]{b}_{\\mathcal{F}}=\\left[{b}_{0} {b}_{1} \\dots {b}_{n}\\right]\\)\u003c/span\u003e\u003c/span\u003e, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({\\beta }_{\\mathcal{F}}=\\left[{\\beta }_{0} {\\beta }_{1} \\dots { \\beta }_{n}\\right]\\)\u003c/span\u003e\u003c/span\u003e.\u003c/p\u003e\n\u003cp\u003eIn this work the FOMCON toolbox for MATLAB/Simulink is used, the toolbox is based on fractional calculus for system modeling and control design. The FOMCON toolbox uses the Levenberg-Marquardt algorithm for optimization of the model parameters \u0026theta; under the nonlinear least squares technique. The Levenberg-Marquardt algorithm is described by the Eq.\u0026nbsp;(\u003cspan class=\"InternalRef\"\u003e16\u003c/span\u003e) According to the references\u003csup\u003e21,35,36\u003c/sup\u003e :\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equ5\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equ5\" name=\"EquationSource\"\u003e$${-J}_{p}^{T}{℮}_{p}=\\varDelta {\\theta }_{p} \\left({J}_{p}^{T}{J}_{p}+\\lambda I\\right)$$\u003c/div\u003e\n \u003cdiv class=\"EquationNumber\"\u003e16\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003eWhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\lambda\\)\u003c/span\u003e\u003c/span\u003e is a positive scaler, called a damping factor, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(I\\)\u003c/span\u003e\u003c/span\u003e is the identity matrix, and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(J\\)\u003c/span\u003e\u003c/span\u003e is the Jacobian matrix of the function \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\mathcal{F}\\)\u003c/span\u003e\u003c/span\u003e and is defined as\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equh\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equh\" name=\"EquationSource\"\u003e$$J=\\left[\\begin{array}{ccc}\\frac{\\partial {℮}_{\\text{1,1}}}{\\partial {\\theta }_{1}}\u0026amp; \\frac{\\partial {℮}_{\\text{1,1}}}{\\partial {\\theta }_{2}} \\dots \u0026amp; \\frac{\\partial {℮}_{\\text{1,1}}}{\\partial {\\theta }_{N}}\\\\ \\begin{array}{c}\\frac{\\partial {℮}_{\\text{1,2}}}{\\partial {\\theta }_{1}}\\\\ \\dots \\end{array}\u0026amp; \\begin{array}{c}\\frac{\\partial {℮}_{\\text{1,2}}}{\\partial {\\theta }_{1}}\\\\ \\dots \\end{array} \\dots \u0026amp; \\begin{array}{c}\\frac{\\partial {℮}_{\\text{1,2}}}{\\partial {\\theta }_{N}}\\\\ \\dots \\end{array}\\\\ \\frac{\\partial {℮}_{1,m}}{\\partial {\\theta }_{1}}\u0026amp; \\frac{\\partial {℮}_{1,m}}{\\partial {\\theta }_{1}} \\dots \u0026amp; \\frac{\\partial {℮}_{1,m}}{\\partial {\\theta }_{N}}\\end{array}\\right] \\left(17\\right)$$\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003ethe update parameters are defined in Eq.\u0026nbsp;(\u003cspan class=\"InternalRef\"\u003e18\u003c/span\u003e).\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equ6\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equ6\" name=\"EquationSource\"\u003e$${\\theta }_{p+1}={\\theta }_{p}+\\varDelta {\\theta }_{p}$$\u003c/div\u003e\n \u003cdiv class=\"EquationNumber\"\u003e18\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003eFrom equations (\u003cspan class=\"InternalRef\"\u003e16\u003c/span\u003e) and (\u003cspan class=\"InternalRef\"\u003e18\u003c/span\u003e) the update parameters can be determined by Eq.\u0026nbsp;(19):\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equi\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equi\" name=\"EquationSource\"\u003e$${\\theta }_{p+1}={\\theta }_{p} - \\frac{{-J}_{p}^{T}{℮}_{p}}{\\left({J}_{p}^{T}{J}_{p}+\\lambda I\\right)}. \\left(19\\right)$$\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003e\u003cstrong\u003e4.4. The model quality\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe quality of the identified model will be discussed in this section. let \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({y}_{s}\\)\u003c/span\u003e\u003c/span\u003e denote The output data of the system, and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({y}_{d}\\)\u003c/span\u003e\u003c/span\u003e denote The output data of the identified model, for single input single output (SISO)system, both \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({y}_{s}\\)\u003c/span\u003e\u003c/span\u003e and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({y}_{d}\\)\u003c/span\u003e\u003c/span\u003e are vectors of dimension \u003cem\u003eN X 1\u003c/em\u003e, and the difference between \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({y}_{s}\\)\u003c/span\u003e\u003c/span\u003e and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({y}_{d}\\)\u003c/span\u003e\u003c/span\u003e is the model output error since \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({y}_{s}\\)\u003c/span\u003e\u003c/span\u003e and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({y}_{d}\\)\u003c/span\u003e\u003c/span\u003e are vectors, then the error is a vector of size \u003cem\u003eN X 1\u003c/em\u003e which is called residuals\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equ7\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equ7\" name=\"EquationSource\"\u003e$$e= {y}_{s}-{y}_{s}$$\u003c/div\u003e\n \u003cdiv class=\"EquationNumber\"\u003e20\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003e.\u003c/p\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003eThe percentage fit can be described as\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equj\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equj\" name=\"EquationSource\"\u003e$$fit=\\left(1-\\frac{\\parallel e\\parallel }{{y}_{s}-\\stackrel{-}{{y}_{s}}}\\right)*100\\% , \\left(21\\right)$$\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003cp\u003ewhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\parallel .\\parallel\\)\u003c/span\u003e\u003c/span\u003e is the Euclidean norm, and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\stackrel{-}{{y}_{s}}\\)\u003c/span\u003e\u003c/span\u003e is the mean value of the system output data\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({y}_{s}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e4.5 stability analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo study the fractional system stability, we consider Matignon\u0026rsquo;s stability theorem\u003csup\u003e37\u003c/sup\u003e :\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eTheorem 1\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e(\u003cem\u003eMatignon\u0026rsquo;s stability theorem\u003c/em\u003e) The fractional transfer function G(s)\u0026thinsp;=\u0026thinsp;Z(s)\u0026thinsp;=\u0026thinsp;P (s) is stable if and only if the following condition is satisfied in the \u0026sigma;-plane:\u003c/p\u003e\n\u003cdiv class=\"Equation\" id=\"Equ8\"\u003e\n \u003cdiv class=\"mathdisplay\" id=\"FileID_Equ8\" name=\"EquationSource\"\u003e$$\\left|\\text{arg}\\left(\\sigma \\right)\\right|\u0026gt;q\\frac{\\pi }{2} , \\forall \\sigma \\in \\complement , P\\left(\\sigma \\right)=0$$\u003c/div\u003e\n \u003cdiv class=\"EquationNumber\"\u003e22\u003c/div\u003e\n\u003c/div\u003e\n\u003cp\u003e,\u003c/p\u003e\n\u003cp\u003e\u003cem\u003ewhere\u003c/em\u003e \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(0 \u0026lt; q \u0026lt; 2\\)\u003c/span\u003e\u003c/span\u003e\u003cem\u003eand\u003c/em\u003e \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\sigma ={s}^{q}\\)\u003c/span\u003e\u003c/span\u003e. \u003cem\u003eWhen\u003c/em\u003e \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\sigma =0\\)\u003c/span\u003e\u003c/span\u003e \u003cem\u003eis a single root of\u003c/em\u003e \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(P \\left(s\\right),\\)\u003c/span\u003e\u003c/span\u003e \u003cem\u003ethe system cannot be stable. For q\u003c/em\u003e = 1, \u003cem\u003ethis is the classical theorem of pole location in the complex plane: no pole is in the closed right plane of the first Riemann sheet.\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003eThe algorithm for determining the system stability in (10) is shown in figure (4), according to Matignon\u0026rsquo;s stability theorem, the fractional order system stability in the complex plane is shown in figure.5.\u003c/p\u003e"},{"header":"Iv. Experiments And Results","content":"\u003cp\u003eFrom Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, The PV current \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({I}_{PV}\\left(n\\right)\\)\u003c/span\u003e\u003c/span\u003e and output power\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({P}_{PV}\\left(n\\right)\\)\u003c/span\u003e\u003c/span\u003e) are used data of the system identification, the measured data is collected under constant irradiation and constant dc load during a 1minute period, for the implementation of the proposed technique, an ARDUINO ATMEGA2560 has been used as a microcontroller, a circuit of the voltage divider and a hall effect current sensor ACS7-12 has been used for reading the analog date for the PV voltage and current respectively, in order to store the signals, a HANTEK 2D42 oscilloscope has been used, to measure the PV current or the output current a shunt resistor of 0.5 ohms is used, the data has been collected with a sample period of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\({T}_{s}=1000\\mu S\\)\u003c/span\u003e\u003c/span\u003e. The collected data are provided for the FOMCON toolbox in MATLAB to develop the fractional order model.\u003c/p\u003e \u003cp\u003eThe measured data of the PV current and the output power of the boost converter collected from the previous experiment are shown in figure(6):\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe PV current and the output power of the boost converter have been used as input data and output data respectively for the FOMCON toolbox in MATLAB, the Levenberg-Marquardt optimization algorithm is used for the optimization problem, and Oustaloup recursive filter approximation is used for model simulation in the time domain, the fractional order model which is estimated by the FOMCON toolbox is shown in the following equation:\u003cdiv id=\"Equk\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equk\" name=\"EquationSource\"\u003eConclusion\n$$G\\left(S\\right)=\\frac{{109.57 S}^{2.5684}+{4374.9.7S}^{1.7303}+3.0654}{{8.0417 S}^{2.8057}+{226.94S}^{1.7615}+{9.347 S}^{1.2891}+1}, \\left(23\\right)$$\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003eFrom Eq.\u0026nbsp;(21) The fitting between the system output data and the output of the model (22) is 93.15%, and the error norm is 213.714, as shown in the figure.7, following theorem \u003cspan refid=\"FPar1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, the poles of the estimated model in (23) are satisfied with the condition (22), such that the identified model is stable for q\u0026thinsp;=\u0026thinsp;0.08 as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eModel validation\u003c/b\u003e \u003c/p\u003e \u003cp\u003eChecking the model validation by comprising the proposed model with another integer model is constructed by MATLAB System Identification TOOLBOX, with the same experimental data. the selection of the integer model was dependent on many iterative processes, the best-fitted integer model shows an accuracy of 88.74%, and the fitting between the measured output data and the output of the integer model is shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e9\u003c/span\u003e. Comparing the identified fractional order model (FOM) and the integer order model (IOM), shows the accuracy of FOM is higher than the IOM accuracy. it is clear from figure.7 and figure.9 that the performance of the FOM is better than the IOM, such that the results show the validity of the proposed technique to construct the FOM of the PV system.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e "},{"header":"Conclusion","content":"\u003cp\u003eIn this paper, a dynamical fractional order model of the PV system is driven based on the system identification technique. The measured input-output data are collected for parameters estimation based on the Levenberg-Marquardt algorithm; the identification process is based on the Least Squares technique for minimization of the errors of the model output. the fractional model stability is checked by Matignon\u0026rsquo;s stability theorem. an integer order model is constructed and its performance is estimated to demonstrate the fractional order model validation, the fractional order model shows higher accuracy (93.15%) than the other model.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate:\u003c/strong\u003e Not applicable\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests:\u003c/strong\u003e The authors declare that they have no conflicts of interest\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding:\u0026nbsp;\u003c/strong\u003eThis research received no external funding.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e: The Authors hereby consents to the publication of the Work in the Energy, Sustainability and Society journal.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of supporting data:\u0026nbsp;\u003c/strong\u003eNo data were used to support this study.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors\u0026apos; contributions:\u003c/strong\u003e R.E. and E.N. conceived and designed the experiments; R.E. and B.A. performed the experiments; R.E. and E.N. \u0026nbsp;prepared figures and tables; R.E and B.A., and E.N. reviewed the manuscript; supervision B.A. E.N.; and R.E. \u0026nbsp; wrote the main manuscript text.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgments:\u0026nbsp;\u003c/strong\u003eI wish to thank Prof. Belal Abou-Zalam, Dr. Essam Nabil, and Prof. Mohmed Hamdy (faculty of electronic engineering, Minoufia University) for supporting in research.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eRezk, H., et al.: A novel statistical performance evaluation of most modern optimization-based global MPPT techniques for partially shaded PV system. Renew. Sustain. Energy Rev. \u003cb\u003e115\u003c/b\u003e, 109372 (2019)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eEbrahimi, S.M., Salahshour, E., Malekzadeh, M.: Francisco Gordillo. Parameters identification of PV solar cells and modules using flexible particle swarm optimization algorithm. Energy. \u003cb\u003e179\u003c/b\u003e, 358\u0026ndash;372 (2019)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBollipo, R.B., Mikkili, S., Bonthagorla, P.K.: Critical Review on PV MPPT Techniques: Classical, Intelligent and Optimisation. \u003cem\u003eIET Renewable Power Generation\u003c/em\u003e vol.\u0026nbsp;14 1433\u0026ndash;1452 Preprint at (2020). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1049/iet-rpg.2019.1163\u003c/span\u003e\u003cspan address=\"10.1049/iet-rpg.2019.1163\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePal, R.S., Mukherjee, V.: Metaheuristic based comparative MPPT methods for photovoltaic technology under partial shading condition. Energy. \u003cb\u003e212\u003c/b\u003e, 118592 (2020)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi, H., Yang, D., Su, W., Lu, J., Yu, X.: An Overall Distribution Particle Swarm Optimization MPPT Algorithm for Photovoltaic System under Partial Shading. IEEE Trans. Industr. Electron. \u003cb\u003e66\u003c/b\u003e, 265\u0026ndash;275 (2019)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAbdelaty, A.M., Radwan, A.G., Elwakil, A.S., Psychalinos, C.: Transient and Steady-State Response of a Fractional-Order Dynamic PV Model under Different Loads. J. Circuits Syst. Computers. \u003cb\u003e27\u003c/b\u003e, 1\u0026ndash;25 (2018)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eda Rocha, M.V., Sampaio, L.P., da Silva, S.A.: O. Comparative analysis of MPPT algorithms based on Bat algorithm for PV systems under partial shading condition. Sustain. Energy Technol. Assess. \u003cb\u003e40\u003c/b\u003e, 100761 (2020)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGhasemi, M.A., Forushani, M., Parniani, M.: Partial shading detection and smooth maximum power point tracking of PV arrays under PSC. IEEE Trans. Power Electron. \u003cb\u003e31\u003c/b\u003e, 6281\u0026ndash;6292 (2016)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAli, A., et al.: Investigation of MPPT Techniques under Uniform and Non-Uniform Solar Irradiation Condition-A Retrospection. IEEE Access. \u003cb\u003e8\u003c/b\u003e, 127368\u0026ndash;127392 (2020)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang, B., et al.: Novel bio-inspired memetic salp swarm algorithm and application to MPPT for PV systems considering partial shading condition. J. Clean. Prod. \u003cb\u003e215\u003c/b\u003e, 1203\u0026ndash;1222 (2019)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAli, K., et al.: Robust integral backstepping based nonlinear MPPT control for a PV system. Energies (Basel). \u003cb\u003e12\u003c/b\u003e, 1\u0026ndash;20 (2019)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePillai, D.S., Rajasekar, N.: Metaheuristic algorithms for PV parameter identification: A comprehensive review with an application to threshold setting for fault detection in PV systems. Renew. Sustain. Energy Rev. \u003cb\u003e82\u003c/b\u003e, 3503\u0026ndash;3525 (2018)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRodrigues, M.P., A MISO Nonlinear Model of Photovoltaic Panel Based on System Identification:. \u003cem\u003eIEEE PES Conference on Innovative Smart Grid Technologies, ISGT Latin America 2019\u003c/em\u003e (2019) doi: (2019). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1109/ISGT-LA.2019.8895325\u003c/span\u003e\u003cspan address=\"10.1109/ISGT-LA.2019.8895325\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHumada, A.M., et al.: Modeling of PV system and parameter extraction based on experimental data: Review and investigation. Sol. Energy. \u003cb\u003e199\u003c/b\u003e, 742\u0026ndash;760 (2020)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDizqah, A.M., Maheri, A., Busawon, K.: An accurate method for the PV model identification based on a genetic algorithm and the interior-point method. Renew. Energy. \u003cb\u003e72\u003c/b\u003e, 212\u0026ndash;222 (2014)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMohamed, M.A.E.H.: The linear model of a PV moduel. Int. J. Power Electron. Drive Syst. \u003cb\u003e8\u003c/b\u003e, 900\u0026ndash;906 (2017)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eElazab, O.S., Hasanien, H.M., Elgendy, M.A., Abdeen, A.M.: Parameters estimation of single- and multiplediode photovoltaic model using whale optimisation algorithm. IET Renew. Power Gener. \u003cb\u003e12\u003c/b\u003e, 1755\u0026ndash;1761 (2018)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRasheed, M.S., Shihab, S.: Modelling and Parameter Extraction of PV Cell Using Single-Diode Model. Adv. Energy Convers. Mater. 96\u0026ndash;104 (2020). doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.37256/aecm.122020550\u003c/span\u003e\u003cspan address=\"10.37256/aecm.122020550\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003edi Piazza, M.C., Luna, M., Vitale, G.: Dynamic PV model parameter identification by least-squares regression. IEEE J. Photovolt. \u003cb\u003e3\u003c/b\u003e, 799\u0026ndash;806 (2013)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAlqahtani, A., Alsaffar, M., El-Sayed, M., Alajmi, B.: Data-Driven Photovoltaic System Modeling Based on Nonlinear System Identification. Int. J. Photoenergy. \u003cb\u003e2016\u003c/b\u003e, 9\u0026ndash;12 (2016)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGavin, H.P.: The Levenberg-Marquardt Algorithm For Nonlinear Least Squares Curve-Fitting Problems. \u003cem\u003eDuke University\u003c/em\u003e 1\u0026ndash;19 (2019)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLuo, A.C.J., Gheorghe, C.M., Machado, D.J.T.: \u003cem\u003eIntroduction to Fractional Dii erential Equations\u003c/em\u003e. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.springer.com/series/11433\u003c/span\u003e\u003cspan address=\"http://www.springer.com/series/11433\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTepljakov, A.: \u003cem\u003eFractional-order modeling and control of dynamic systems\u003c/em\u003e.Automaticavol. 21 (2016)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTeh, C.J.Q., Drieberg, M., Soeung, S., Ahmad, R.: Simple PV Modeling under Variable Operating Conditions. IEEE Access. \u003cb\u003e9\u003c/b\u003e, 96546\u0026ndash;96558 (2021)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBlaifi, S., ali, Moulahoum, S., Taghezouit, B., Saim, A.: An enhanced dynamic modeling of PV module using Levenberg-Marquardt algorithm. Renew. Energy. \u003cb\u003e135\u003c/b\u003e, 745\u0026ndash;760 (2019)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNasser, K.W., Yaqoob, S.J., Hassoun, Z.A.: Improved dynamic performance of photovoltaic panel using fuzzy logic-MPPT algorithm. Indonesian J. Electr. Eng. Comput. Sci. \u003cb\u003e21\u003c/b\u003e, 617\u0026ndash;624 (2020)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAl-Soeidat, M., Cheng, T., Lu, D.D.C., Agelidis, V.G.: Experimental study of static and dynamic behaviours of cracked PV panels. IET Renew. Power Gener. \u003cb\u003e13\u003c/b\u003e, 3002\u0026ndash;3008 (2019)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCelikel, R., Gundogdu, A.: System identification-based MPPT algorithm for PV systems under variable atmosphere conditions using current sensorless approach. Int. Trans. Electr. Energy Syst. \u003cb\u003e30\u003c/b\u003e, 1\u0026ndash;21 (2020)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSj\u0026ouml;berg, J., et al.: Nonlinear black-box modeling in system identification: a unified overview. Automatica. \u003cb\u003e31\u003c/b\u003e, 1691\u0026ndash;1724 (1995)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAlqahtani, A., Marafi, S., Musallam, B.: Abd El Khalek, N. E. D. Photovoltaic Power Forecasting Model Based on Nonlinear System Identification. Can. J. Electr. Comput. Eng. \u003cb\u003e39\u003c/b\u003e, 243\u0026ndash;250 (2016)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHcheichi, K., Bouani, F.: Comparison between commensurate and non-commensurate fractional systems. Int. J. Adv. Comput. Sci. Appl. \u003cb\u003e9\u003c/b\u003e, 685\u0026ndash;691 (2018)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMarouani, H., Hergli, K., Dhahri, H., Fouad, Y.: Implementation and Identification of Preisach Parameters: Comparison Between Genetic Algorithm, Particle Swarm Optimization, and Levenberg\u0026ndash;Marquardt Algorithm. Arab. J. Sci. Eng. \u003cb\u003e44\u003c/b\u003e, 6941\u0026ndash;6949 (2019)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi, J., Liu, Y., Shan, C., Dai, C.: Implementation of Simplified Fractional-Order PID Controller Based on Modified Oustaloup\u0026rsquo;s Recursive Filter. J. Shanghai Jiaotong Univ. (Science). \u003cb\u003e25\u003c/b\u003e, 44\u0026ndash;50 (2020)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCol\u0026iacute;n-Cervantes, J.D., et al.: Rational approximations of arbitrary order: A survey.Fractal and Fractional\u003cb\u003e5\u003c/b\u003e, (2021)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eIslam, S.U., et al.: Design of robust fuzzy logic controller based on the levenberg marquardt algorithm and fault ride trough strategies for a grid-connected PV system.Electronics (Switzerland)\u003cb\u003e8\u003c/b\u003e, (2019)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMarouani, H., Hergli, K., Dhahri, H., Fouad, Y.: Implementation and Identification of Preisach Parameters: Comparison Between Genetic Algorithm, Particle Swarm Optimization, and Levenberg\u0026ndash;Marquardt Algorithm. Arab. J. Sci. Eng. \u003cb\u003e44\u003c/b\u003e, 6941\u0026ndash;6949 (2019)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eD, M.: Generalized fractional differential and difference equations: stability properties and modelling issues. \u003cem\u003eProceedings of the Mathematics Theory of Networks and Systems Symposium\u003c/em\u003e 503\u0026ndash;506 (1998)\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":true,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"journal-of-computational-electronics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"jcel","sideBox":"Learn more about [Journal of Computational Electronics](https://www.springer.com/journal/10825)","snPcode":"10825","submissionUrl":"https://submission.nature.com/new-submission/10825/3","title":"Journal of Computational Electronics","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"fractional-order model, System identification, PV system, least-square curve fitting, and PV system.","lastPublishedDoi":"10.21203/rs.3.rs-1865340/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-1865340/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eA fractional-order model of a photovoltaic (PV) system is proposed in this paper. The system identification method is used to construct an accurate dynamical model for a PV system, a real PV module and boost converter are used to collect experimentally input-output data for the identification process. A black box modeling is considered for system identification to obtain a transfer function, without needing to make mathematical analysis. the input-output data which is measured is fitted by the least-square curve fitting method, and the parameters of the curve fitting are optimized by the Levenberg-Marquardt algorithm. the proposed fractional-order model is employed with the MATLAB and Simulink software to obtain a fractional order PID controller, which is used to extract maximum power from the PV system.\u003c/p\u003e","manuscriptTitle":"System identification of photovoltaic system based on fractional-order model","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2022-07-29 20:09:45","doi":"10.21203/rs.3.rs-1865340/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"reviewerAgreed","content":"","date":"2022-07-23T10:45:36+00:00","index":0,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2022-07-22T16:56:36+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"Journal of Computational Electronics","date":"2022-07-22T14:44:41+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2022-07-18T09:27:48+00:00","index":"","fulltext":""},{"type":"submitted","content":"Journal of Computational Electronics","date":"2022-07-16T15:35:57+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"journal-of-computational-electronics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"jcel","sideBox":"Learn more about [Journal of Computational Electronics](https://www.springer.com/journal/10825)","snPcode":"10825","submissionUrl":"https://submission.nature.com/new-submission/10825/3","title":"Journal of Computational Electronics","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"914de916-7547-4cf8-a725-6d566ce87b04","owner":[],"postedDate":"July 29th, 2022","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"under-review","subjectAreas":[],"tags":[],"updatedAt":"2022-12-03T15:09:36+00:00","versionOfRecord":[],"versionCreatedAt":"2022-07-29 20:09:45","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-1865340","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-1865340","identity":"rs-1865340","version":["v1"]},"buildId":"cBFmMYwuxLRRLfASyISRj","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}