<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JPEE</journal-id><journal-title-group><journal-title>Journal of Power and Energy Engineering</journal-title></journal-title-group><issn pub-type="epub">2327-588X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jpee.2014.23001</article-id><article-id pub-id-type="publisher-id">JPEE-44281</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Dynamic State Forecasting in Electric Power Networks
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>ideig</surname><given-names>A. Dowi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Amar</surname><given-names>Ibrahim Hamza</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>School of Electrical and Electronics Engineering, North China Electric Power University, Beijing, China</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>sideigdowi@yahoo.com(IAD)</email>;<email>amarhamaza2010@hotmail.com(AIH)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>28</day><month>03</month><year>2014</year></pub-date><volume>02</volume><issue>03</issue><fpage>1</fpage><lpage>11</lpage><history><date date-type="received"><day>3</day>	<month>February</month>	<year>2014</year></date><date date-type="rev-recd"><day>3</day>	<month>March</month>	<year>2014</year>	</date><date date-type="accepted"><day>10</day>	<month>March</month>	<year>2014</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   The real time monitoring and control have become very important in electric power system in order to achieve a high reliability in the system. So, improvement in Energy Management System (EMS) leads to improvement in the monitoring and control functions in the control center. In this paper, DSE is proposed based on Weighted Least Squares (WLS) estimator and Holt’s exponential smoothing to state predicting and Extended Kalman Filter to state filtering. The results viewing the dynamic state the estimator performance under normal and abnormal operating conditions. 
 
</p></abstract><kwd-group><kwd>Dynamic State Estimation; Extended Kalman Filter; State Estimation; Electric Power Flow</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Recently, the power system has begun to grow very largely and more complex, so real time monitoring and control become very important in order to fulfill a reliable operation. Energy Management System (EMS) is responsible for this mission, and it forms the basis for efficient operating and control. State Estimation (SE) forms the spine of the EMS by providing the information of the real time state of the system which can be used in other EMS functions. Hence, an accurate and efficient state estimation is necessary for a reliable and efficient operation of the power system [<xref ref-type="bibr" rid="scirp.44281-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.44281-ref2">2</xref>] .</p><p>The state estimator computes the voltage magnitudes and voltage angles at the buses of the power system. We know that, power system is not a static system, but it changes very slowly with time and continuously. That means, when the load on the buses changes, the generations also have to change to overcome these changes in load. This in turn causes the change in power flows and injections at the buses, also leads to change in voltage angular at the buses and perhaps change in voltage magnitude at some buses depending on the size of this change; therefore, change the nature of the power system from static state to dynamic state nature. These dynamic behaviors of the power system are difficult to overcome by the conventional Static State Estimation (SSE). This led to the development of a new algorithm called Dynamic State Estimation (DSE) [<xref ref-type="bibr" rid="scirp.44281-ref1">1</xref>] -[<xref ref-type="bibr" rid="scirp.44281-ref5">5</xref>]. The DSE technique possesses a mathematical model for the time variation of the power system. The DSE uses this mathematical model depending on the previous state of the system at time (t), to predict the state vectors at the next state of the system (one step ahead) at time (t + 1). This capability of predicting the state one step ahead is a very important advantage in the control center, because the state forecasting gives a longer decision time to the system operator, and because the security assessment, economic dispatching and the other functions can be performed in advance. So, Dynamic State Estimation has the ability to represent an important role in the modern-day control center [<xref ref-type="bibr" rid="scirp.44281-ref1">1</xref>] .</p><p>In this paper, we will not only describe the dynamic model for the time behavior of the system state, but will show more details about the DSE, mainly the state predicting and state filtering. When the state variables are estimated at time k by state estimation technique, we will use these state variables to forecast the state vectors at time k + 1 using linear exponential smoothing. The state vectors are filtered based on Extended Kalman Filter and weighted least squares method. The proposal is tested using IEEE 14 bus test system. The test includes normal and abnormal operations.</p></sec><sec id="s2"><title>2. Mathematical Models</title><p>The measurement vector consists of active power and reactive power flows and injection’s power as well as some voltage magnitudes, is denoted by an m-dimensional vector z. The power equation s is expressed by [<xref ref-type="bibr" rid="scirp.44281-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.44281-ref7">7</xref>] .</p><disp-formula id="scirp.44281-formula2305"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x5.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.44281-formula2306"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x6.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.44281-formula2307"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x7.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.44281-formula2308"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x8.png"  xlink:type="simple"/></disp-formula><p>The measurement and state variable <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x9.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x10.png" xlink:type="simple"/></inline-formula> at time instant k are given by equation</p><disp-formula id="scirp.44281-formula2309"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x11.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x12.png" xlink:type="simple"/></inline-formula> is measurement vector m-dimensional<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x13.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x14.png" xlink:type="simple"/></inline-formula>is state variables dimensioned<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x15.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x16.png" xlink:type="simple"/></inline-formula>is nonlinear function relating measurement to the state vector .dimensioned <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x17.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x18.png" xlink:type="simple"/></inline-formula> is the measurements error with zero mean and standard deviation of R.</p><p>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x19.png" xlink:type="simple"/></inline-formula>, σ is the standard deviation of error.</p><disp-formula id="scirp.44281-formula2310"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x20.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x21.png" xlink:type="simple"/></inline-formula> is Jacobian matrix dimensional (m&#180;n), m is number of measurements and<sup> </sup>n is number of state vectors.</p><p>The general model for DSE is given by.</p><disp-formula id="scirp.44281-formula2311"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x22.png"  xlink:type="simple"/></disp-formula><p>where x<sub>k</sub> and x<sub>k</sub><sub>+1 </sub>are the state vector at instants k and k+1 respectively, F<sub>k</sub> is nonzero diagonal matrix dimensioned (n &#180; n), a function represent the state transition between two instant of time, G<sub>k</sub><sub> </sub>is nonzero vector associated with trend behavior of the state trajectory dimensional (n &#180; 1) and w<sub>k</sub> is white Gaussian noise with zero mean and covariance matrix Q [<xref ref-type="bibr" rid="scirp.44281-ref1">1</xref>] -[<xref ref-type="bibr" rid="scirp.44281-ref4">4</xref>].</p><sec id="s2_1"><title>2.1. Parameters Identification</title><p>The parameters F<sub>k</sub> and G<sub>k</sub> are identified using Holt’s two-parameter linear exponential smoothing method [<xref ref-type="bibr" rid="scirp.44281-ref1">1</xref>] , [<xref ref-type="bibr" rid="scirp.44281-ref8">8</xref>] -[<xref ref-type="bibr" rid="scirp.44281-ref10">10</xref>]. This method is very simple, used when the data shows a trend. In this method, two components must be updated each period of time, level and trend.</p><p>・ The level is a smoothed estimate of the value of the data at the end of each period represented by a<sub>k</sub>. as shown in Equation (8).</p><p>・ The trend is a smoothed estimate of average growth at the end of each period. represented by<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x23.png" xlink:type="simple"/></inline-formula>.as shown in Equation (8).</p><p>The specific model for simple exponential smoothing is written as:</p><disp-formula id="scirp.44281-formula2312"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x24.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x25.png" xlink:type="simple"/></inline-formula> &amp; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x26.png" xlink:type="simple"/></inline-formula> 0 &lt; β &lt; 1, 0 &lt; α &lt; 1</p><disp-formula id="scirp.44281-formula2313"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x27.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.44281-formula2314"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x28.png"  xlink:type="simple"/></disp-formula><p>α and β represent the smoothing parameters [<xref ref-type="bibr" rid="scirp.44281-ref2">2</xref>] , [<xref ref-type="bibr" rid="scirp.44281-ref8">8</xref>] -[<xref ref-type="bibr" rid="scirp.44281-ref10">10</xref>].</p></sec><sec id="s2_2"><title>2.2. State Predicting or State Forecasting</title><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x29.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x30.png" xlink:type="simple"/></inline-formula> be the estimated state and estimated state covariance respectively at time k the forecasted state vector and its covariance matrix is obtained by</p><disp-formula id="scirp.44281-formula2315"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x31.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.44281-formula2316"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x32.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x33.png" xlink:type="simple"/></inline-formula> is covariance matrix of the forecasted state.</p></sec><sec id="s2_3"><title>2.3. State Filtering</title><p>The forecasted state would use to forecast new measurements Z<sub>k</sub><sub>+1</sub> at the time k + 1 based on the data at instant k; the predicted state vector at k + 1 will be filtered to obtain new estimates (filtered states) <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x34.png" xlink:type="simple"/></inline-formula>with its error’s covariance matrix<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x35.png" xlink:type="simple"/></inline-formula>. Then, the objective function for the filtering process at the instant of time (k + 1) is</p><disp-formula id="scirp.44281-formula2317"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x36.png"  xlink:type="simple"/></disp-formula><p>Note that, the time index (k + 1) has been omitted.</p><p>Extended Kalman Filter (EKF) used for minimizing the objective function and getting the final filtering state.</p><disp-formula id="scirp.44281-formula2318"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x37.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.44281-formula2319"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x38.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.44281-formula2320"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x39.png"  xlink:type="simple"/></disp-formula><p>where K is called the gain matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x40.png" xlink:type="simple"/></inline-formula> dimensioned. Perform only one iteration in Equation (15).</p></sec></sec><sec id="s3"><title>3. The Implementation</title><p>The steps of the dynamic state estimation algorithm are described above. The covariance matrix R of the measurement error is assumed to be calculated online.</p><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x41.png" xlink:type="simple"/></inline-formula></p><p>where a<sub>m</sub> &amp; b<sub>m</sub> are the manufacture factors. s<sub>t</sub> represents the real value of the measurement. s<sub>f</sub> is the maximum value of the measurement. As mentioned before, for predicting state we used Holt’s 2-exponential smoothing, the values of the smoothing parameters α and β, are fixed at 0.7 and 0.45, and for the filtering state we used Extended Kalman Filter. The elements of the covariance matrix Q of the system, is set at 10<sup>‒6</sup>. The load curve at each bus was composed of a linear trend and random fluctuation (jitter).</p><p>In this paper, the Dynamic State simulation is studied over a period of 20 time sample intervals. with increasing of constant value 5% of the load at all the buses at each period. Once the load is changed the load flow is ready to update all the real and reactive power and injection power on the lines, voltage magnitudes and angles on the system. In this paper, the actual values of the state vectors are obtained from the 14-bus IEEE standard data for the base case and load flow for rest of the time samples [<xref ref-type="bibr" rid="scirp.44281-ref11">11</xref>] . The test cases and their results are shown in tables and figures in the next sections.</p><sec id="s3_1"><title>3.1. Case Study</title><p>In this paper, we used standard IEEE 14-bus test system [<xref ref-type="bibr" rid="scirp.44281-ref11">11</xref>] . The measurement vector consists of 34 observations distributed as in <xref ref-type="table" rid="table1">Table 1</xref>.</p><p>The measurement value was simulated by adding random errors to the true values represented by normally distribution with zero mean and standard deviation.</p></sec><sec id="s3_2"><title>3.2. Performance Indices</title><p>The performance of the algorithm in the simulation studies was obtained by comparing the forecasted and estimated values at time k + 1 with the actual values. The average performance indices for voltage magnitude and voltage angular forecasted and estimated are given as.</p><disp-formula id="scirp.44281-formula2321"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x42.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.44281-formula2322"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x43.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.44281-formula2323"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1770021x44.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x45.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x46.png" xlink:type="simple"/></inline-formula> are the true value of voltage magnitude and angle, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x47.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x48.png" xlink:type="simple"/></inline-formula> are the forecasted and filtered voltage magnitude, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x49.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1770021x50.png" xlink:type="simple"/></inline-formula> are the forecasted and filtered voltage angle, J<sub>k</sub> is the performance index of the overall achievement, the ratio of estimated at time k and actual error of the measurement.</p></sec><sec id="s3_3"><title>3.3. Results Discussion</title><p>In this paper, three test cases are performed, the normal operation, bad data and sudden load change.</p><p>Test 1: Normal operation case:</p><p>The normal operation case was illustrated by <xref ref-type="table" rid="table2">Table 2</xref>, <xref ref-type="table" rid="table3">Table 3</xref>, Figures 1 and 2. <xref ref-type="table" rid="table2">Table 2</xref> shows the perfor-</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> The Measurements prepared by the load flow</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Type</th><th align="center" valign="middle" >Measurement vector</th></tr></thead><tr><td align="center" valign="middle" >Active power flow</td><td align="center" valign="middle" >p(1-5), p(4-5), p(4-9), p(6-11), p(6-12), p(7-8), p(7-9), p(9-10) , p(13-14)</td></tr><tr><td align="center" valign="middle" >Reactive power flow</td><td align="center" valign="middle" >q(1-5), q(4-5), q(4-9), q(6-11), q(6-12), q(7-8), q(7-9), q(9-10), q(13-14)</td></tr><tr><td align="center" valign="middle" >Active and Reactive power injection</td><td align="center" valign="middle" >P<sub>1</sub>, P<sub>3</sub>, P<sub>6</sub>, P<sub>10</sub>, P<sub>12.</sub> Q<sub>1</sub>, Q<sub>3</sub><sub> </sub>, Q<sub>6</sub>, Q<sub>10</sub>, Q<sub>12</sub></td></tr><tr><td align="center" valign="middle" >Voltages magnitude</td><td align="center" valign="middle" >V<sub>1</sub>, V<sub>3</sub>, V<sub>8</sub>, V<sub>11</sub>, V<sub>12</sub>, V<sub>14.</sub></td></tr></tbody></table></table-wrap><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Performance indices under normal operation</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  colspan="2"   rowspan="2"  >cases</th><th align="center" valign="middle"  colspan="2"  >Predicted</th><th align="center" valign="middle"  colspan="2"  >Filtered</th><th align="center" valign="middle"  rowspan="2"  >J_k</th></tr></thead><tr><td align="center" valign="middle" >voltage</td><td align="center" valign="middle" >angle</td><td align="center" valign="middle" >voltage</td><td align="center" valign="middle" >angle</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Normal operation</td><td align="center" valign="middle" >Max</td><td align="center" valign="middle" >0.510</td><td align="center" valign="middle" >1.7546</td><td align="center" valign="middle" >0.499</td><td align="center" valign="middle" >1.0455</td><td align="center" valign="middle" >0.970</td></tr><tr><td align="center" valign="middle" >Ave</td><td align="center" valign="middle" >0.1845</td><td align="center" valign="middle" >1.3450</td><td align="center" valign="middle" >0.178</td><td align="center" valign="middle" >0.6385</td><td align="center" valign="middle" >0.957</td></tr></tbody></table></table-wrap><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> The errors per bus under normal operation</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >Bus number</th><th align="center" valign="middle"  colspan="2"  >predicted</th><th align="center" valign="middle"  colspan="2"  >Filtered</th></tr></thead><tr><td align="center" valign="middle" >voltage</td><td align="center" valign="middle" >Angle</td><td align="center" valign="middle" >voltage</td><td align="center" valign="middle" >angle</td></tr><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.2286</td><td align="center" valign="middle" >0.2198</td><td align="center" valign="middle" >0.000</td><td align="center" valign="middle" >0.0000</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >0.1977</td><td align="center" valign="middle" >0.1986</td><td align="center" valign="middle" >1.440</td><td align="center" valign="middle" >0.6792</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >0.2969</td><td align="center" valign="middle" >0.2954</td><td align="center" valign="middle" >1.666</td><td align="center" valign="middle" >0.9037</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >0.2480</td><td align="center" valign="middle" >0.2389</td><td align="center" valign="middle" >1.530</td><td align="center" valign="middle" >0.7662</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >0.2257</td><td align="center" valign="middle" >0.2168</td><td align="center" valign="middle" >1.458</td><td align="center" valign="middle" >0.6974</td></tr><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle" >0.1560</td><td align="center" valign="middle" >0.1530</td><td align="center" valign="middle" >1.409</td><td align="center" valign="middle" >0.6483</td></tr><tr><td align="center" valign="middle" >7</td><td align="center" valign="middle" >0.2090</td><td align="center" valign="middle" >0.2080</td><td align="center" valign="middle" >1.457</td><td align="center" valign="middle" >0.6967</td></tr><tr><td align="center" valign="middle" >8</td><td align="center" valign="middle" >0.1621</td><td align="center" valign="middle" >0.1568</td><td align="center" valign="middle" >1.436</td><td align="center" valign="middle" >0.6754</td></tr><tr><td align="center" valign="middle" >9</td><td align="center" valign="middle" >0.2036</td><td align="center" valign="middle" >0.1995</td><td align="center" valign="middle" >1.444</td><td align="center" valign="middle" >0.6838</td></tr><tr><td align="center" valign="middle" >10</td><td align="center" valign="middle" >0.1953</td><td align="center" valign="middle" >0.1907</td><td align="center" valign="middle" >1.438</td><td align="center" valign="middle" >0.6780</td></tr><tr><td align="center" valign="middle" >11</td><td align="center" valign="middle" >0.1228</td><td align="center" valign="middle" >0.1207</td><td align="center" valign="middle" >1.395</td><td align="center" valign="middle" >0.6350</td></tr><tr><td align="center" valign="middle" >12</td><td align="center" valign="middle" >0.1420</td><td align="center" valign="middle" >0.1426</td><td align="center" valign="middle" >1.398</td><td align="center" valign="middle" >0.6375</td></tr><tr><td align="center" valign="middle" >13</td><td align="center" valign="middle" >0.0901</td><td align="center" valign="middle" >0.0870</td><td align="center" valign="middle" >1.389</td><td align="center" valign="middle" >0.6294</td></tr><tr><td align="center" valign="middle" >14</td><td align="center" valign="middle" >0.1055</td><td align="center" valign="middle" >0.0650</td><td align="center" valign="middle" >1.369</td><td align="center" valign="middle" >0.6093</td></tr></tbody></table></table-wrap><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Performance index under normal operation</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1770021x51.png"/></fig><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> V9 under normal operation</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1770021x52.png"/></fig><p>mance indices of the case. The maximum and average percentage error of voltage magnitudes and angles are made over 20 time sample. From this table, the algorithm has achieved very high performance.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref> describes the Performance index of the system dynamic nature under normal operation case corresponding to filtered indices calculated at each time sample.</p><p><xref ref-type="table" rid="table3">Table 3</xref> shows the percentage for estimated errors per bus for predicted and filtered state with regard to the voltage magnitudes and voltage angles.</p><p>In both states, the maximum error occurred in busbar 3 for voltage magnitudes and voltage angles. Furthermore the average errors are equal to the average error in <xref ref-type="table" rid="table2">Table 2</xref>. Bus 1 is the reference bus, so it is angle is outside of state variables. <xref ref-type="fig" rid="fig3">Figure 3</xref> represents the time behavior for true, forecasted and filtered values of voltage magnitude at bus 9. Figures 3 and 4 show the graphs of the <xref ref-type="table" rid="table3">Table 3</xref>.</p><p>Test 2: Bad data case:</p><p>In this test, the simulation was carried out under bad data conditions with three different cases.</p><p>・ Single bad data was considered. Active power flow pf(4 - 9) was suspected in error of 10% at the 6th time sample.</p><p>・ Two measurements were suspected in error pf(9 - 10) of 20% at the 11th time sample. And also pf(13 - 14) of 50% decrement at the 11th time sample.</p><p>・ Single bad data was considered. Reactive power flow qf(4 - 5) was suspected in error of 20% at the 18<sup>th</sup> time sample.</p><p>Suppose that no work is taken to eliminate these errors.</p><p>The result of these tests is shown in Tables 4 and 5, Figures 5-7.</p><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Average voltage error per bus *100</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1770021x53.png"/></fig><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Average voltage angle error per bur*100</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1770021x54.png"/></fig><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Performance indices under bad data condition</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  colspan="2"   rowspan="2"  >cases</th><th align="center" valign="middle"  colspan="2"  >Predicted</th><th align="center" valign="middle"  colspan="2"  >Filtered</th><th align="center" valign="middle"  rowspan="2"  >J_k</th></tr></thead><tr><td align="center" valign="middle" >voltage</td><td align="center" valign="middle" >angle</td><td align="center" valign="middle" >voltage</td><td align="center" valign="middle" >angle</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Bad date</td><td align="center" valign="middle" >Max</td><td align="center" valign="middle" >0.501</td><td align="center" valign="middle" >3.2594</td><td align="center" valign="middle" >0.484</td><td align="center" valign="middle" >2.539</td><td align="center" valign="middle" >4.3631</td></tr><tr><td align="center" valign="middle" >Ave</td><td align="center" valign="middle" >0.210</td><td align="center" valign="middle" >2.501</td><td align="center" valign="middle" >0.206</td><td align="center" valign="middle" >1.786</td><td align="center" valign="middle" >1.6755</td></tr></tbody></table></table-wrap><table-wrap id="table5" ><label><xref ref-type="table" rid="table5">Table 5</xref></label><caption><title> The percentage error per bus under bad data conditions</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >Bus number</th><th align="center" valign="middle"  colspan="2"  >Predicted</th><th align="center" valign="middle"  colspan="2"  >Filtered</th></tr></thead><tr><td align="center" valign="middle" >Voltage</td><td align="center" valign="middle" >angle</td><td align="center" valign="middle" >voltage</td><td align="center" valign="middle" >angle</td></tr><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.2264</td><td align="center" valign="middle" >0.0000</td><td align="center" valign="middle" >0.2283</td><td align="center" valign="middle" >0.0000</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >0.2172</td><td align="center" valign="middle" >1.5183</td><td align="center" valign="middle" >0.2179</td><td align="center" valign="middle" >0.7570</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >0.3045</td><td align="center" valign="middle" >1.6960</td><td align="center" valign="middle" >0.3039</td><td align="center" valign="middle" >0.9337</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >0.2500</td><td align="center" valign="middle" >1.5674</td><td align="center" valign="middle" >0.2563</td><td align="center" valign="middle" >0.8061</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >0.2730</td><td align="center" valign="middle" >1.4960</td><td align="center" valign="middle" >0.2578</td><td align="center" valign="middle" >0.7351</td></tr><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle" >0.1571</td><td align="center" valign="middle" >3.1423</td><td align="center" valign="middle" >0.1533</td><td align="center" valign="middle" >2.3690</td></tr><tr><td align="center" valign="middle" >7</td><td align="center" valign="middle" >0.2581</td><td align="center" valign="middle" >3.1045</td><td align="center" valign="middle" >0.2637</td><td align="center" valign="middle" >2.3318</td></tr><tr><td align="center" valign="middle" >8</td><td align="center" valign="middle" >0.2271</td><td align="center" valign="middle" >3.0803</td><td align="center" valign="middle" >0.2144</td><td align="center" valign="middle" >2.3103</td></tr><tr><td align="center" valign="middle" >9</td><td align="center" valign="middle" >0.2418</td><td align="center" valign="middle" >2.8826</td><td align="center" valign="middle" >0.2536</td><td align="center" valign="middle" >2.1113</td></tr><tr><td align="center" valign="middle" >10</td><td align="center" valign="middle" >0.2171</td><td align="center" valign="middle" >2.9726</td><td align="center" valign="middle" >0.2266</td><td align="center" valign="middle" >2.2005</td></tr><tr><td align="center" valign="middle" >11</td><td align="center" valign="middle" >0.1176</td><td align="center" valign="middle" >3.1732</td><td align="center" valign="middle" >0.1227</td><td align="center" valign="middle" >2.3997</td></tr><tr><td align="center" valign="middle" >12</td><td align="center" valign="middle" >0.1352</td><td align="center" valign="middle" >3.3294</td><td align="center" valign="middle" >0.1352</td><td align="center" valign="middle" >2.5548</td></tr><tr><td align="center" valign="middle" >13</td><td align="center" valign="middle" >0.1606</td><td align="center" valign="middle" >3.6942</td><td align="center" valign="middle" >0.1470</td><td align="center" valign="middle" >2.9168</td></tr><tr><td align="center" valign="middle" >14</td><td align="center" valign="middle" >0.1462</td><td align="center" valign="middle" >3.5880</td><td align="center" valign="middle" >0.1025</td><td align="center" valign="middle" >2.5840</td></tr></tbody></table></table-wrap><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Performance index under bad data conditio</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1770021x55.png"/></fig><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Average voltage error per bus *100</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1770021x56.png"/></fig><p><xref ref-type="table" rid="table4">Table 4</xref> shows the performance indices of these tests. From this results the average and maximum error becomes larger comparing with normal operation test due to these bad data in measurements.</p><p><xref ref-type="fig" rid="fig5">Figure 5</xref> describes the Performance index of the system dynamic nature under bad data condition’s corresponding to filtered indices calculated at each time sample.</p><p>Figures 6 and 7 show the percentage error per bus of the voltage magnitudes and voltage angle’s respectively based on <xref ref-type="table" rid="table5">Table 5</xref>.</p><p>Test 3: Sudden load change:</p><p>In the case, the injected load at busbar 3, busbar 9 and busbar 13 are assumed to be sudden changed.</p><p>Busbar 3 and busbar 13, 50% of their values are cut at the 6th time sample and 15th time sample respectively. For busbar 9, we assumed that the load is increased to 40% at the 10th time sample.</p><p>The result of this case is shown in Tables 6-8, and also shown in Figures 8-10. <xref ref-type="table" rid="table6">Table 6</xref> shows the performance indices of the case of sudden load change. From this results the average and maximum error has become difference of normal operation and bad data test due to these changes in load.</p><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Average voltage angle error per bus *100</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1770021x57.png"/></fig><table-wrap id="table6" ><label><xref ref-type="table" rid="table6">Table 6</xref></label><caption><title> Performance indices under sudden load change</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  colspan="2"   rowspan="2"  >cases</th><th align="center" valign="middle"  colspan="2"  >Predicted</th><th align="center" valign="middle"  colspan="2"  >Filtered</th><th align="center" valign="middle"  rowspan="2"  >J_k</th></tr></thead><tr><td align="center" valign="middle" >voltage</td><td align="center" valign="middle" >angle</td><td align="center" valign="middle" >voltage</td><td align="center" valign="middle" >angle</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Sudden Load change</td><td align="center" valign="middle" >Max</td><td align="center" valign="middle" >0.5145</td><td align="center" valign="middle" >1.6942</td><td align="center" valign="middle" >0.4951</td><td align="center" valign="middle" >0.9857</td><td align="center" valign="middle" >0.9866</td></tr><tr><td align="center" valign="middle" >Ave</td><td align="center" valign="middle" >0.195</td><td align="center" valign="middle" >1.35</td><td align="center" valign="middle" >0.1922</td><td align="center" valign="middle" >0.642</td><td align="center" valign="middle" >0.9748</td></tr></tbody></table></table-wrap><table-wrap id="table7" ><label><xref ref-type="table" rid="table7">Table 7</xref></label><caption><title> The percentage error per bus under sudden load change</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >Bus number</th><th align="center" valign="middle"  colspan="2"  >Predicted</th><th align="center" valign="middle"  colspan="2"  >Filtered</th></tr></thead><tr><td align="center" valign="middle" >voltage</td><td align="center" valign="middle" >angle</td><td align="center" valign="middle" >voltage</td><td align="center" valign="middle" >angle</td></tr><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.2546</td><td align="center" valign="middle" >0.0000</td><td align="center" valign="middle" >0.2465</td><td align="center" valign="middle" >0.0000</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >0.2258</td><td align="center" valign="middle" >1.4344</td><td align="center" valign="middle" >0.2268</td><td align="center" valign="middle" >0.6735</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >0.2860</td><td align="center" valign="middle" >1.6345</td><td align="center" valign="middle" >0.2910</td><td align="center" valign="middle" >0.8728</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >0.2525</td><td align="center" valign="middle" >1.5383</td><td align="center" valign="middle" >0.2541</td><td align="center" valign="middle" >0.7771</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >0.2318</td><td align="center" valign="middle" >1.4830</td><td align="center" valign="middle" >0.2341</td><td align="center" valign="middle" >0.7226</td></tr><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle" >0.1710</td><td align="center" valign="middle" >1.4125</td><td align="center" valign="middle" >0.1683</td><td align="center" valign="middle" >0.6522</td></tr><tr><td align="center" valign="middle" >7</td><td align="center" valign="middle" >0.2270</td><td align="center" valign="middle" >1.4660</td><td align="center" valign="middle" >0.2310</td><td align="center" valign="middle" >0.7057</td></tr><tr><td align="center" valign="middle" >8</td><td align="center" valign="middle" >0.1733</td><td align="center" valign="middle" >1.4453</td><td align="center" valign="middle" >0.1680</td><td align="center" valign="middle" >0.6847</td></tr><tr><td align="center" valign="middle" >9</td><td align="center" valign="middle" >0.2155</td><td align="center" valign="middle" >1.4520</td><td align="center" valign="middle" >0.2180</td><td align="center" valign="middle" >0.6912</td></tr><tr><td align="center" valign="middle" >10</td><td align="center" valign="middle" >0.2058</td><td align="center" valign="middle" >1.4460</td><td align="center" valign="middle" >0.2072</td><td align="center" valign="middle" >0.6852</td></tr><tr><td align="center" valign="middle" >11</td><td align="center" valign="middle" >0.1300</td><td align="center" valign="middle" >1.4002</td><td align="center" valign="middle" >0.1306</td><td align="center" valign="middle" >0.6400</td></tr><tr><td align="center" valign="middle" >12</td><td align="center" valign="middle" >0.1546</td><td align="center" valign="middle" >1.4025</td><td align="center" valign="middle" >0.1563</td><td align="center" valign="middle" >0.6422</td></tr><tr><td align="center" valign="middle" >13</td><td align="center" valign="middle" >0.0934</td><td align="center" valign="middle" >1.3900</td><td align="center" valign="middle" >0.0938</td><td align="center" valign="middle" >0.6296</td></tr><tr><td align="center" valign="middle" >14</td><td align="center" valign="middle" >0.1092</td><td align="center" valign="middle" >1.3712</td><td align="center" valign="middle" >0.0651</td><td align="center" valign="middle" >0.6113</td></tr></tbody></table></table-wrap><table-wrap id="table8" ><label><xref ref-type="table" rid="table8">Table 8</xref></label><caption><title> V13 under sudden load change</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Time sample</th><th align="center" valign="middle" >True</th><th align="center" valign="middle" >predicted</th><th align="center" valign="middle" >filtered</th></tr></thead><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >1.0501</td><td align="center" valign="middle" >1.0482</td><td align="center" valign="middle" >1.0478</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >1.0489</td><td align="center" valign="middle" >1.0485</td><td align="center" valign="middle" >1.0481</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >1.0474</td><td align="center" valign="middle" >1.0482</td><td align="center" valign="middle" >1.0479</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >1.0459</td><td align="center" valign="middle" >1.0467</td><td align="center" valign="middle" >1.0464</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >1.0246</td><td align="center" valign="middle" >1.0237</td><td align="center" valign="middle" >1.0235</td></tr><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle" >1.0336</td><td align="center" valign="middle" >1.034</td><td align="center" valign="middle" >1.0338</td></tr><tr><td align="center" valign="middle" >7</td><td align="center" valign="middle" >1.0173</td><td align="center" valign="middle" >1.0165</td><td align="center" valign="middle" >1.0164</td></tr><tr><td align="center" valign="middle" >8</td><td align="center" valign="middle" >1.0056</td><td align="center" valign="middle" >1.0058</td><td align="center" valign="middle" >1.0058</td></tr><tr><td align="center" valign="middle" >9</td><td align="center" valign="middle" >0.9937</td><td align="center" valign="middle" >0.991</td><td align="center" valign="middle" >0.991</td></tr><tr><td align="center" valign="middle" >10</td><td align="center" valign="middle" >0.9827</td><td align="center" valign="middle" >0.9823</td><td align="center" valign="middle" >0.9824</td></tr><tr><td align="center" valign="middle" >11</td><td align="center" valign="middle" >0.9800</td><td align="center" valign="middle" >0.9791</td><td align="center" valign="middle" >0.9793</td></tr><tr><td align="center" valign="middle" >12</td><td align="center" valign="middle" >0.9769</td><td align="center" valign="middle" >0.9769</td><td align="center" valign="middle" >0.977</td></tr><tr><td align="center" valign="middle" >13</td><td align="center" valign="middle" >0.9738</td><td align="center" valign="middle" >0.9737</td><td align="center" valign="middle" >0.9739</td></tr><tr><td align="center" valign="middle" >14</td><td align="center" valign="middle" >0.9706</td><td align="center" valign="middle" >0.9723</td><td align="center" valign="middle" >0.9725</td></tr><tr><td align="center" valign="middle" >15</td><td align="center" valign="middle" >0.9798</td><td align="center" valign="middle" >0.9804</td><td align="center" valign="middle" >0.9805</td></tr><tr><td align="center" valign="middle" >16</td><td align="center" valign="middle" >0.9771</td><td align="center" valign="middle" >0.9745</td><td align="center" valign="middle" >0.9747</td></tr><tr><td align="center" valign="middle" >17</td><td align="center" valign="middle" >0.9744</td><td align="center" valign="middle" >0.9728</td><td align="center" valign="middle" >0.973</td></tr><tr><td align="center" valign="middle" >18</td><td align="center" valign="middle" >0.971</td><td align="center" valign="middle" >0.9708</td><td align="center" valign="middle" >0.971</td></tr><tr><td align="center" valign="middle" >19</td><td align="center" valign="middle" >0.9675</td><td align="center" valign="middle" >0.9681</td><td align="center" valign="middle" >0.9683</td></tr><tr><td align="center" valign="middle" >20</td><td align="center" valign="middle" >0.9641</td><td align="center" valign="middle" >0.9649</td><td align="center" valign="middle" >0.9652</td></tr></tbody></table></table-wrap><fig id="fig8"  position="float"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> Average voltage error per bus *100%</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1770021x58.png"/></fig><fig id="fig9"  position="float"><label><xref ref-type="fig" rid="fig9">Figure 9</xref></label><caption><title> Average voltage angle error per bus *100</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1770021x59.png"/></fig><fig id="fig10"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>0</label><caption><title> V13 under sudden load change condition</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1770021x60.png"/></fig><p><xref ref-type="fig" rid="fig1">Figure 1</xref>0 represents the time behavior for true, forecasted and filtered values of voltage magnitude at busbar 13 based on <xref ref-type="table" rid="table8">Table 8</xref>. As they are shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>0, at time sample 6, 10 and 15 when sudden load change, the oscillation in the voltage is appeared very clear.</p></sec></sec><sec id="s4"><title>4. Conclusion</title><p>The dynamic state estimator technique has been made based on Holt’s exponential smoothing and Extended Kalman Filter for forecasting and filtering state respectively. The system dynamic was simulated over 20 time samples, with increasing the load at all busbars 5% during any time sample. The algorithm of DSE has been simulated through these 20 time intervals used standard IEEE 14_bus test system under normal and abnormal operation. In this paper, the algorithm gave very good performance results through the normal operation and abnormal operation (bad data and sudden load change) conditions. The error per bus for voltage magnitude and voltage angle has been calculated over these 20 time intervals in both predicted and filtered states.</p></sec></body><back><ref-list><title>References</title><ref id="scirp.44281-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Shivakumar, N.R. and Jain, A.M. (2008) A Review of Power System Dynamic State Estimation Techniques. Joint International Conference on Power System Technology and IEEE Power India Conference, New Delhi, 15 October 2008, 1-6.</mixed-citation></ref><ref id="scirp.44281-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Jegatheesa, R. and Leean, K. (1991) Dynamic State Estimation. TENCON’91.1991 IEEE Region 10 International Conference on EC3-Energy, Computer, Communication and Control Systems, New Delhi, 28-30 Aug 1991, 492-496.</mixed-citation></ref><ref id="scirp.44281-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Debs, A.S. and Larson, R.E. (1970) A Dynamic Estimator for Tracking the State of a Power System. IEEE Transactions on Power Apparatus and Systems, PAS-89, 1670-1678.</mixed-citation></ref><ref id="scirp.44281-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Jain, A.M. and Shivakumar, N.R. (2008) Phasor Measurement in Dynamic State Estimation of a Power System. TENCON 2008, 2008 IEEE Region 10 Conference, Hyderabad, 19-21 November 2008, 1-6.</mixed-citation></ref><ref id="scirp.44281-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Huang, S.J. and Shih, K.S. (2002) Dynamic-State-Estimation Scheme including Nonlinear Measurement-Function Considerations. Generation, Transmission and Distribution, IEE Proceedings, 149, 673-678.</mixed-citation></ref><ref id="scirp.44281-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Ali, A. and Antonio, G.E. (2004) Power System State Estimation Theory and Implementation. CRC Press, New York.</mixed-citation></ref><ref id="scirp.44281-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Mukhtar, A. (2013) Power System State Estimation. Artech House, London.</mixed-citation></ref><ref id="scirp.44281-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Leite da silva, A.M., Coutto Filho, M.B. and de Queiroz, J.F. (1976) State Forecasting in Electric Power Systems. Generation, Transmission and Distribution, IEE Proceedings C, 130, 79-87.</mixed-citation></ref><ref id="scirp.44281-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Li, H. and Li, W.G. (2009) Estimation and Forecasting of Dynamic State Estimation in Power Systems. International Conference on Sustainable Power Generation and Supply, Nanjing, 6-7 April 2009, 1-6.</mixed-citation></ref><ref id="scirp.44281-ref10"><label>10</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Prajakta</surname><given-names> S.K. </given-names></name>,<etal>et al</etal>. (<year>2004</year>)<article-title>Time Series Forecasting Using Holt-Winters Exponential Smoothing</article-title><source> Kanwal Rekhi School of Information Technology</source><volume> 4329008</volume>,<fpage> 1</fpage>-<lpage>13</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.44281-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">University of Washington Electrical Engineering (1993) Power Systems Test Case Archive. http://www.ee.washington.edu/research/pstca/</mixed-citation></ref></ref-list></back></article>