<?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">JSIP</journal-id><journal-title-group><journal-title>Journal of Signal and Information Processing</journal-title></journal-title-group><issn pub-type="epub">2159-4465</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jsip.2017.83008</article-id><article-id pub-id-type="publisher-id">JSIP-77747</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject></subj-group></article-categories><title-group><article-title>
 
 
  Fading Channels Parametric Data Simulation Supported by Real Data from Outdoor Experiments
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Azra</surname><given-names>Kapetanovic</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>Mohamed</surname><given-names>A. Zohdy</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Redhwan</surname><given-names>Mawari</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Electrical and Computer Engineering, Oakland University, Rochester, MI, USA</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>azra.kapetanovic@gm.com(AK)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>18</day><month>07</month><year>2017</year></pub-date><volume>08</volume><issue>03</issue><fpage>113</fpage><lpage>131</lpage><history><date date-type="received"><day>April</day>	<month>2,</month>	<year>2017</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>July</month>	<year>16,</year>	</date><date date-type="accepted"><day>July</day>	<month>19,</month>	<year>2017</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>
 
 
  Optimizing the estimates of received power signals is important as it can improve the process of transferring an active call from one base station in a cellular network to another base station without any interruptions to the call. The lack of effective techniques for estimation of shadow power in fading mobile wireless communication channels motivated the use of Kalman Filtering (KF) as an effective alternative. In our research, linear second-order state space Kalman Filtering was further investigated and tested for applicability. We first created simulation models for two KF-based estimators designed to estimate local mean (shadow) power in mobile communications corrupted by multipath noise. Simulations were used extensively in the initial stage of this research to validate the proposed method. The next challenge was to determine if the models would work with real data. Therefore, in [1] we presented a new technique to experimentally characterize the wireless small-scale fading channel taking into consideration real environmental conditions. The two-dimensional measurement technique enabled us to perform indoor experiments and collect real data. Measurements from these experiments were then used to validate simulation models for both estimators. Based on the indoor experiments, we presented new results in [2], where we concluded that the second-order KF-based estimator is more accurate in predicting local shadow power profiles than the first-order KF-based estimator, even in channels with imposed non-Gaussian measurement noise. In the present paper, we extend experiments to the outdoor environment to include higher speeds, larger distances, and distant large objects, such as tall buildings. Comparison was performed to see if the system is able to operate without a failure under a variety of conditions, which demonstrates model robustness and further investigates the effectiveness of this method in optimization of the received signals. Outdoor experimental results are provided. Findings demonstrate that the second-order Kalman filter outperforms the first-order Kalman filter.
 
</p></abstract><kwd-group><kwd>Kalman Filtering</kwd><kwd> Small-Scale Fading</kwd><kwd> Large-Scale Fading</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Because wireless technology and smart cell phones are experiencing dramatic growth, the accurate estimation of local mean (shadow) power in a cell phone is becoming a popular area of challenge for engineers in both industry and academia. Researchers are encouraged to find ways to enhance device performance in power control and handoff, particularly to address mobility-induced fading in metropolitan areas.</p><p>Wireless cell phones operate by transferring information over a distance between two or more stations that are not connected by cables. Instead, cell phones use radio waves to carry information, such as sound, by systematically modulating some property of electromagnetic energy waves transmitted through space, such as their amplitude, frequency, phase, or pulse width [<xref ref-type="bibr" rid="scirp.77747-ref3">3</xref>] . A transmitted signal from a cell tower will undergo changes while traveling through the propagation path to the cell phone. These changes may fluctuate with time, geographical position, and radio frequency. The term fading is used to describe the effect of these changes. As a result, the quality of communications decreases.</p><p>Two significant forms of fading in cellular communications are multipath and shadow fading. Since cell phone users tend to move a lot, received signal strength fluctuates with these two multiplicative forms of fading [<xref ref-type="bibr" rid="scirp.77747-ref4">4</xref>] . In the outdoor environment, shadow fading in cell phones causes long-term variation primarily caused by nearby mountains or tall buildings. Tall building structures shadow the radio signal, which results in a power drop at a receiver input. Multipath, as illustrated in <xref ref-type="fig" rid="fig1">Figure 1</xref>, results in fluctuations of the signal amplitude because of the addition of signals that arrive with different phases [<xref ref-type="bibr" rid="scirp.77747-ref2">2</xref>] .</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Versions of the same signal can take many different paths between cell tower and cell phone which causes fading effect in the resultant received signal</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-3400502x2.png"/></fig><p>Having an accurate estimate of the shadowing component of a received power signal will allow the mobile communication system to efficiently compensate for the signal degradation that will occur. As a result, it can help the system perform handoff at the most effective times (predict when and where to handoff user). In [<xref ref-type="bibr" rid="scirp.77747-ref5">5</xref>] , the authors have presented and discussed the limitation of different types of windows-based estimators that are utilized to filter multipath noise from the instantaneous received power signal to estimate the local mean shadow power. Unfortunately, window-based estimators work well only under the assumption that the shadowing component is relatively constant during the window period. However, shadowing components are not constant; the fluctuations can vary and at times can significantly decrease the performance of the windows-based estimate. In [<xref ref-type="bibr" rid="scirp.77747-ref2">2</xref>] , we proposed a second-order KF-based estimator as an alternative method to windows-based estimation for estimating local mean (shadow) power.</p><p>Kalman filtering is a very effective algorithm that uses a series of measured observations and produces optimal estimates of states as explained in [<xref ref-type="bibr" rid="scirp.77747-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.77747-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.77747-ref8">8</xref>] . We applied this method to derive equations for an estimator that can estimate local mean (shadow) power profiles. In the initial stage of research, we used simulation models to validate the proposed method. The next challenge was to see if the model would work with real data. Therefore, in [<xref ref-type="bibr" rid="scirp.77747-ref1">1</xref>] we presented a new technique to experimentally characterize the wireless small-scale fading channel, taking into consideration real environmental conditions. This new two dimensional measurement technique provided essential information regarding the constructive and destructive interference patterns caused by the interaction between the mobile station (while in motion) and surrounding obstacles. The two dimensional measurement technique enabled us to perform indoor experi- ments and collect real data, which we then used to confirm validation of the simulation model for the second-order KF-based estimator.Based on results from the indoor experiments, we concluded that the second-order KF-based estimator is more accurate in predicting local mean (shadow) power profiles than the first-order KF-based estimator, even in channels with imposed non-Gaussian measurement noise.</p><p>The next challenge was to find a way to measure cell phone signal strength outside of the lab environment and to test our second-order KF-based estimator. In this paper, we looked at mobility-induced fading and present experimental results from the outdoor environment that further confirmed validation of the proposed method. We will explain how Kalman filter method can be applied in optimization of received signals in mobile communications. The system was able to operate without a failure under a variety of conditions, which demonstrates model robustness. In subsequent sections, it will be demonstrated that the second-order KF based estimator we designed exceeds the performance of the first-order KF-based estimator, even in the outdoor environment where parameters for mobile velocity varied.</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Model for the received power at cell phone</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-3400502x3.png"/></fig></sec><sec id="s2"><title>2. Model for Multipath Signal</title><p>The description of Shadow Power Signal and its models as they pertain to our problem are presented in this section. In a wireless cellular radio environment, a model for an instantaneous received power signal l(t) at a cell phone is given in <xref ref-type="fig" rid="fig2">Figure 2</xref>, where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x4.png" xlink:type="simple"/></inline-formula> represents fast power fluctuation due to multipath and x(t) represents slow power fluctuation due to shadowing. Many common shadow power estimation methods in industry rely on an accurate model for multipath. Multipath is often modeled as Rayleigh noise for modeling purposes. It is customary to express power measurements in decibels. Handoff algorithms rely on estimates of shadow power in decibels [<xref ref-type="bibr" rid="scirp.77747-ref9">9</xref>] .</p><p>To solve the problem, we start with the multipath model shown in Equation (2).</p><disp-formula id="scirp.77747-formula1"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x5.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula2"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x6.png"  xlink:type="simple"/></disp-formula><p>where:</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x7.png" xlink:type="simple"/></inline-formula>is the magnitude of the mobile velocity [10 m/s - 30 m/s],</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x8.png" xlink:type="simple"/></inline-formula>is the wavelength corresponding to the carrier frequency, which is typically<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x9.png" xlink:type="simple"/></inline-formula>,</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x10.png" xlink:type="simple"/></inline-formula>is the Doppler spread equal to<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x11.png" xlink:type="simple"/></inline-formula>. Range used here is [<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x12.png" xlink:type="simple"/></inline-formula>-<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x13.png" xlink:type="simple"/></inline-formula>],</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x14.png" xlink:type="simple"/></inline-formula>is the number of paths in the multipath power [20 - 30],</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x15.png" xlink:type="simple"/></inline-formula>is gain [0 - 20],</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x16.png" xlink:type="simple"/></inline-formula>represents angles between incoming waves and mobile antenna. The value range is uniformly distributed<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x17.png" xlink:type="simple"/></inline-formula>. Other distributions like Normal can be used,</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x18.png" xlink:type="simple"/></inline-formula>represents phase random variables whose values are also uniformly distributed<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x19.png" xlink:type="simple"/></inline-formula>.</p></sec><sec id="s3"><title>3. Model for Shadow Power Signal</title><p>Equation (3) shows a widely accepted first-order state space model for the shadow process given by [<xref ref-type="bibr" rid="scirp.77747-ref10">10</xref>] . We derived a second-order state space model for the shadow process as shown in Equation (4).</p><disp-formula id="scirp.77747-formula3"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x20.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula4"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x21.png"  xlink:type="simple"/></disp-formula><p>First shadow power coefficient <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x22.png" xlink:type="simple"/></inline-formula> is given by Equation (5) and second shadow power coefficient <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x23.png" xlink:type="simple"/></inline-formula> is given by Equation (6):</p><disp-formula id="scirp.77747-formula5"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x24.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula6"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x25.png"  xlink:type="simple"/></disp-formula><p>where, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x26.png" xlink:type="simple"/></inline-formula>is the magnitude of the mobile velocity, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x27.png" xlink:type="simple"/></inline-formula>is time sample, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x28.png" xlink:type="simple"/></inline-formula> is effective correlation distance. The effective correlation distance is key attribute of the wireless environment. In urban area it can be as low as 10 m while in suburban areas it can be as high as 500 m. The effective correlation distance is given in Equation (7), where variable <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x29.png" xlink:type="simple"/></inline-formula> is the distance between cell tower and cell phone measured in meters. Term <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x30.png" xlink:type="simple"/></inline-formula> is the correlation coefficient of shadow process between two points separated by distance D.</p><disp-formula id="scirp.77747-formula7"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x31.png"  xlink:type="simple"/></disp-formula><p>Finally, system noise covariance is given in Equation (8), where term <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x32.png" xlink:type="simple"/></inline-formula> denotes the shadow variance which depends on environment. In urban areas, typical value for shadow variance is 4 dB while in suburban areas typical value is 8 dB. Term <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x33.png" xlink:type="simple"/></inline-formula> in Equations (3) and (4) denotes zero mean white Gaussian noise with variance<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x34.png" xlink:type="simple"/></inline-formula>.</p><disp-formula id="scirp.77747-formula8"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x35.png"  xlink:type="simple"/></disp-formula></sec><sec id="s4"><title>4. Kalman Filter Algorithm</title><p>Kalman Filter theory was developed and introduced in 1960 by Dr. Rudolf Kalman. This led to the use of the Kalman filter during the Apollo program, carried out by NASA, which accomplished landing the first humans on the Moon. Since then, his contributions and thoughts educated and inspired inventors across many disciplines. As a result, the Kalman filter has been the subject of extensive research and application, especially in digital computing.</p><p>In this research, we applied the Kalman filter algorithm to estimate power signal in a mobile communication corrupted by multipath noise. The Kalman filter is a form of a linear algorithm for optimal recursive estimation of a system state with a specific set of output equations. The estimates are calculated every time a new measurement is received. Data received is processed sequentially, so it is not necessary to store the complete data set or to reprocess existing data when new measurement data is received.</p><sec id="s4_1"><title>4.1. Derivation of the Discrete-Time Linear Kalman Filter</title><p>This section derives the equations of the discrete-time Kalman filter. This filter is applied as a recursive solution to the estimation problem studied in this research. To use the Kalman filter to estimate signal of interest, one must first create matrices to fix the system model into a Kalman filter format. The following sets of equations describe the format of the linear discrete-time system:</p><disp-formula id="scirp.77747-formula9"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x36.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula10"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x37.png"  xlink:type="simple"/></disp-formula><p>The Kalman filter is a great tool, but its computation is complex and requires some explanation. An optimal value for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x38.png" xlink:type="simple"/></inline-formula> in Equation (9) is calculated based on the available knowledge of the system dynamics and the noise measurement<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x39.png" xlink:type="simple"/></inline-formula>. In Equation (10), <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x40.png" xlink:type="simple"/></inline-formula>represents the measured output of the system (measurement of system state) with measurement noise.</p><p>The k’s on the subscripts are states and can be treated as discrete time intervals. In general, when applying the Kalman Filter, the goal is to estimate state<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x41.png" xlink:type="simple"/></inline-formula>. For example, in signal processing, it is basically the estimate of some signal x that we want to find for each subsequent k. During this process, the Kalman filter forms ana priori estimate and an a posterior estimate denoted as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x42.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x43.png" xlink:type="simple"/></inline-formula>. Equation (11) computes the expected value of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x44.png" xlink:type="simple"/></inline-formula> conditioned on all of the measurement up to time k. Similarly, Equation (12) computes the expected value of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x45.png" xlink:type="simple"/></inline-formula> conditioned on all of the measurements after time k.</p><disp-formula id="scirp.77747-formula11"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x46.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula12"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x47.png"  xlink:type="simple"/></disp-formula><p>Each system has to have initial values. Notation <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x48.png" xlink:type="simple"/></inline-formula> denotes an initial estimate of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x49.png" xlink:type="simple"/></inline-formula> before any measurements are taken. The first measurement in this algorithm is taken at time k = 1. During this time period, the system does not have any measurements available to estimate<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x50.png" xlink:type="simple"/></inline-formula>, and therefore <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x51.png" xlink:type="simple"/></inline-formula> is formed as the best expected value of the initial state <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x52.png" xlink:type="simple"/></inline-formula></p><disp-formula id="scirp.77747-formula13"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x53.png"  xlink:type="simple"/></disp-formula><p>This algorithm takes into account the measurement noise, process noise, and the previous estimated output values so that it can minimize the prediction error upon a continuous cycle of prediction and filtering. It looks at the error between the true state and the estimated state. Therefore, the next step is to derive an error equation. Equations (14) and (15) define an priori estimate error and an posteriori estimate error, respectively.</p><disp-formula id="scirp.77747-formula14"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x54.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula15"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x55.png"  xlink:type="simple"/></disp-formula><p>Then, Equations (14) and (15) are used to compute covariance of the estimation error, which is denoted as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x56.png" xlink:type="simple"/></inline-formula>. The term <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x57.png" xlink:type="simple"/></inline-formula> denotes the covariance of the estimation error of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x58.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x59.png" xlink:type="simple"/></inline-formula> denotes the covariance of the estimation error of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x60.png" xlink:type="simple"/></inline-formula>:</p><disp-formula id="scirp.77747-formula16"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x61.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula17"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x62.png"  xlink:type="simple"/></disp-formula><p>After the measurement at time k − 1 is processed, the estimate of the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x63.png" xlink:type="simple"/></inline-formula> is computed, which is denoted as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x64.png" xlink:type="simple"/></inline-formula>. Also, the covariance of that estimate is computed at the same time, and it is denoted as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x65.png" xlink:type="simple"/></inline-formula>. Then, at time k, before the measurement is processed, estimate of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x66.png" xlink:type="simple"/></inline-formula> is computed and denoted as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x67.png" xlink:type="simple"/></inline-formula> Then, the covariance of these estimates are computed and denoted as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x68.png" xlink:type="simple"/></inline-formula> Then, at time k = 1, the measurement is processed to refine the estimate of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x69.png" xlink:type="simple"/></inline-formula>. The resulting improved estimates of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x70.png" xlink:type="simple"/></inline-formula> and its covariance are denoted as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x71.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x72.png" xlink:type="simple"/></inline-formula>.</p><p>To begin the estimation process, initial values of the system determined in Equation (13) must be initialized. Then, with<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x73.png" xlink:type="simple"/></inline-formula>, variable <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x74.png" xlink:type="simple"/></inline-formula> is computed using Equation (18). Based on this equation, the time update equation for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x74.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x75.png" xlink:type="simple"/></inline-formula> is computed as indicated in Equation (19).</p><disp-formula id="scirp.77747-formula18"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x76.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula19"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x77.png"  xlink:type="simple"/></disp-formula><p>From time <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x78.png" xlink:type="simple"/></inline-formula> to time<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x78.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x79.png" xlink:type="simple"/></inline-formula>, the state estimate and its mean propagates the same way. As there are no additional measurements available between these two time-steps, the state estimate has to be updated based on a knowledge of system dynamics.</p><p>Next step is to derive time-update equation for the covariance of the state estimation error. The term <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x80.png" xlink:type="simple"/></inline-formula> represents the covariance of the initial estimate of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x81.png" xlink:type="simple"/></inline-formula> (the uncertainty in initial estimates of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x81.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x82.png" xlink:type="simple"/></inline-formula>). If the exact value of the initial state is known, then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x81.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x82.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x83.png" xlink:type="simple"/></inline-formula> can be set to zero. However, if the initial value is not known, then then<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x81.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x82.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x84.png" xlink:type="simple"/></inline-formula>. The mathematical from for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x81.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x82.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x85.png" xlink:type="simple"/></inline-formula> is shown in Equation (20). The general description of how the covariance of the state of a linear discrete-time system propagates with time is given by Equation (21). The term <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x81.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x82.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x86.png" xlink:type="simple"/></inline-formula> can be computed by substituting <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x81.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x82.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x86.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x87.png" xlink:type="simple"/></inline-formula> value from Equation (20). The time equation for term P in general form is given by Equation (23).</p><disp-formula id="scirp.77747-formula20"><label>(20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x88.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula21"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x89.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula22"><label>(22)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x90.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula23"><label>(23)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x91.png"  xlink:type="simple"/></disp-formula><p>Final step requires derivation of measurement-update equations for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x92.png" xlink:type="simple"/></inline-formula> and P. Given <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x93.png" xlink:type="simple"/></inline-formula> from Equation (19), we need to find a method to compute<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x94.png" xlink:type="simple"/></inline-formula>, the estimate of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x95.png" xlink:type="simple"/></inline-formula>, which takes the measurement <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x95.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x96.png" xlink:type="simple"/></inline-formula> into account. Recall <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x95.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x96.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x97.png" xlink:type="simple"/></inline-formula> in Equation (10) represents the measurement of the system state with measurement noise<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x95.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x96.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x98.png" xlink:type="simple"/></inline-formula>. Measurement noise is usually caused by the measurement instrument. Based on the recursive least square estimation theory, we know that the availability of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x95.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x96.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x98.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x99.png" xlink:type="simple"/></inline-formula> changes the estimate of a constants <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x95.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x96.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x98.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x99.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x100.png" xlink:type="simple"/></inline-formula> as follows:</p><disp-formula id="scirp.77747-formula24"><label>(24)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x101.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula25"><label>(25)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x102.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula26"><label>(26)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x103.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x105.png" xlink:type="simple"/></inline-formula> are estimates before the measurement is processed, while <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x106.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x107.png" xlink:type="simple"/></inline-formula> take the measurements <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x108.png" xlink:type="simple"/></inline-formula> into account. The next step is to rewrite Equations (24) through (26) in a format that Kalman used when he derived his estimation theory. To formulate the measurement-update equations for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x109.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x110.png" xlink:type="simple"/></inline-formula>, simply perform the following substitutions: substitute <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x111.png" xlink:type="simple"/></inline-formula> with<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x111.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x112.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x111.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x112.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x113.png" xlink:type="simple"/></inline-formula>with<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x111.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x112.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x113.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x114.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x111.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x112.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x113.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x114.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x115.png" xlink:type="simple"/></inline-formula>with<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x111.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x112.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x113.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x114.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x115.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x116.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x111.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x112.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x113.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x114.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x115.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x116.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x117.png" xlink:type="simple"/></inline-formula> with<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x111.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x112.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x113.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x114.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x115.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x116.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x117.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x118.png" xlink:type="simple"/></inline-formula>. These substitutions lead to the following equations:</p><disp-formula id="scirp.77747-formula27"><label>(27)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x119.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula28"><label>(28)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x120.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula29"><label>(29)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x121.png"  xlink:type="simple"/></disp-formula><p>The matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x122.png" xlink:type="simple"/></inline-formula> given in Equation (27) is called the Kalman filter gain. This gain is a blending factor that minimizes the a posteriori error covariance. If the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x123.png" xlink:type="simple"/></inline-formula> is a constant, then<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x124.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x125.png" xlink:type="simple"/></inline-formula>, and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x126.png" xlink:type="simple"/></inline-formula>.</p><p>The random variable <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x127.png" xlink:type="simple"/></inline-formula> represents process noise and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x127.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x128.png" xlink:type="simple"/></inline-formula> represents measurement noise. Terms <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x127.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x128.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x129.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x127.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x128.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x129.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x130.png" xlink:type="simple"/></inline-formula> are independent of each other. They are assumed to be white and with normal probability distributions</p><disp-formula id="scirp.77747-formula30"><label>(30)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x131.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula31"><label>(31)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x132.png"  xlink:type="simple"/></disp-formula><p>In outdoor experiments, the process noise covariance matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x133.png" xlink:type="simple"/></inline-formula> and measurement noise covariance matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x133.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x134.png" xlink:type="simple"/></inline-formula> can change for every measurement.</p></sec><sec id="s4_2"><title>4.2. First-Order Kalman Filter Application to Fading Channels</title><p>The Kalman Filter is a form of a linear algorithm for optimal recursive estimation of a system state with a specific set of output equations. To build a simulator, understanding of system model and its dynamic behaviors is necessary. Then, the system must be represented in the state space format to be able to apply Kalman filtering. In other words, we need to mathematically model its states and parameters. This section presents set of equations used to create first-order KF-based estimator. References [<xref ref-type="bibr" rid="scirp.77747-ref10">10</xref>] and [<xref ref-type="bibr" rid="scirp.77747-ref11">11</xref>] were useful for programming in MATLAB during the initial stages of research.</p><p>To build an estimation model in MATLAB, we started with equations introduced in Section 4.1 and substituted suitable entries from this problem to reflect the linear channel model [<xref ref-type="bibr" rid="scirp.77747-ref12">12</xref>] . Equations (9) and (10) can be rewritten as Equations (32) and (33). This assumes a linear time invariant system with a mean of zero and white noise on both the state and output.</p><disp-formula id="scirp.77747-formula32"><label>(32)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x135.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula33"><label>(33)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x136.png"  xlink:type="simple"/></disp-formula><p>where:</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x137.png" xlink:type="simple"/></inline-formula>is a symbol value of shadow signal state that needs to be estimated.</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x138.png" xlink:type="simple"/></inline-formula>is the control signal for handoff.</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x139.png" xlink:type="simple"/></inline-formula>is process white noise.</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x140.png" xlink:type="simple"/></inline-formula>is measurement noise.</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x141.png" xlink:type="simple"/></inline-formula>is the measurement value for both shadow and multipath. In this document, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x141.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x142.png" xlink:type="simple"/></inline-formula>or L(k) is the measurement value used to update shadow power estimate.</p><p>After initializing Kalman filter using initial values <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x143.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x143.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x144.png" xlink:type="simple"/></inline-formula> as explained in Section 4.1, time-update and measurement-update equations had to be identified. These equations are computed for each time step k = 1, 2, 3. Equations (34) and (35) represent the “Time Update” state of the Kalman Filter, also known as the “Prediction States.” Equations (34) and (35) are derived by substituting suitable entries from this problem into Equations (19) and (22). For first-order KF, matrix A was set to<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x143.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x144.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x145.png" xlink:type="simple"/></inline-formula>, which represents the first shadow power coefficient given by Equation (5.5). Term for system noise covariance <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x143.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x144.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x145.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x146.png" xlink:type="simple"/></inline-formula> in Equation (22), was substituted with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x143.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x144.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x145.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x146.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x147.png" xlink:type="simple"/></inline-formula> in Equation (35) and it denotes variance defined in Equation (34). Here we project the current state estimate forward in time with Equation (34) projecting the state ahead and Equation (35) projecting the error covariance ahead as represented below:</p><disp-formula id="scirp.77747-formula34"><label>(34)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x148.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula35"><label>(35)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x149.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x150.png" xlink:type="simple"/></inline-formula> is the rough estimate before the measurement <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x150.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x151.png" xlink:type="simple"/></inline-formula> is processed at time k.</p><p>Similarly, Equations (37), (38) and (39) belong to the “Measurement Update” state of the linear Kalman Filter, also known as the “Correction State.” Equations (37), (38) and (39) are derived by substituting suitable entries from this problem into Equations (27), (28) and (29). Here we adjust the projected estimate by an actual measurement at time k. Equation (37) was derived by substituting environment noise covariance <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x152.png" xlink:type="simple"/></inline-formula> with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x152.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x153.png" xlink:type="simple"/></inline-formula> defined in Equation (36).</p><disp-formula id="scirp.77747-formula36"><label>(36)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x154.png"  xlink:type="simple"/></disp-formula><p>Equation (37) computes the Kalman Gain, Equation (38) adjusts the projected estimate by an actual measurement<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x155.png" xlink:type="simple"/></inline-formula>, and Equation (39) updates the error covariance, as follows:</p><disp-formula id="scirp.77747-formula37"><label>(37)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x156.png"  xlink:type="simple"/></disp-formula><p>If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x157.png" xlink:type="simple"/></inline-formula> is small and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x157.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x158.png" xlink:type="simple"/></inline-formula> is close to Identity, then<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x157.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x158.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x159.png" xlink:type="simple"/></inline-formula>, which is a well-known Pseudo inverse.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x160.png" xlink:type="simple"/></inline-formula>: Estimate of x after the actual measurement <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x160.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x161.png" xlink:type="simple"/></inline-formula> at time k.</p><disp-formula id="scirp.77747-formula38"><label>(38)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x162.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula39"><label>(39)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x163.png"  xlink:type="simple"/></disp-formula><p>The next step is to represent these estimates over a period of sufficient time. The output estimate in the previous step will be the input estimate in the next step. The main goal is to find an optimal value for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x164.png" xlink:type="simple"/></inline-formula></p></sec><sec id="s4_3"><title>4.3. Second-Order Kalman Filter Application to Fading Channels</title><p>In this research, we assumed that the first-order state space model can be used to model Shadow Power. To extend the first-order state space model equations [<xref ref-type="bibr" rid="scirp.77747-ref1">1</xref>] presented in Section 4.2, a second-order state space model for the linear Kalman Filter was formulated and applied as suggested in the equations below. With this notation, we can describe an algorithm for the second-order KF as follows:</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x165.png" xlink:type="simple"/></inline-formula>is the first Shadow power coefficient as defined in Equation (5).</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x166.png" xlink:type="simple"/></inline-formula>is the second Shadow power coefficient as defined in Equation (6).</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x167.png" xlink:type="simple"/></inline-formula>is the effective correlation distance as defined in Equation (7).</p><p>・ <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x168.png" xlink:type="simple"/></inline-formula>is the correlation coefficient of the shadow process between two points separated by a distance D as measured in meters,</p><disp-formula id="scirp.77747-formula40"><label>(40)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x169.png"  xlink:type="simple"/></disp-formula><p>Equation (40) shows <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x170.png" xlink:type="simple"/></inline-formula> expressed in matrix format for second-order state space.</p><p>Next set of equations present prediction states for second-order linear Kalman filter. Equation (34) can be rewritten as Equation (41). Therefore, Equation (41) in this section projects the state ahead, and Equation (43) projects the error covariance ahead. For the second-order KF, the matrices are defined as follows:</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x171.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x171.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x172.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x171.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x172.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x173.png" xlink:type="simple"/></inline-formula></p><disp-formula id="scirp.77747-formula41"><label>(41)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x174.png"  xlink:type="simple"/></disp-formula><p>B is input matrix that relates the control input to the state <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x175.png" xlink:type="simple"/></inline-formula> Matrix H is output equation whose function is to relate state to the measured output <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x175.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x176.png" xlink:type="simple"/></inline-formula></p><p>Parameter Q in Equation (42) represents the predicted process noise. Term <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x177.png" xlink:type="simple"/></inline-formula> denotes the shadow variance with range from 4 dB to 8 dB. The notation <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x177.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x178.png" xlink:type="simple"/></inline-formula> in Equation (42) denotes the variance of the zero mean white Gaussian noise.</p><disp-formula id="scirp.77747-formula42"><label>(42)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x179.png"  xlink:type="simple"/></disp-formula><p>Equation (43) can then be expressed in the following state space format:</p><disp-formula id="scirp.77747-formula43"><label>(43)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x180.png"  xlink:type="simple"/></disp-formula><p>Next set of equations present correction states. Equation (44) was used to compute Kalman gain, which takes into consideration measurement noise due to multipath.</p><disp-formula id="scirp.77747-formula44"><graphic  xlink:href="http://html.scirp.org/file/1-3400502x181.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula45"><graphic  xlink:href="http://html.scirp.org/file/1-3400502x182.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula46"><label>(44)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x183.png"  xlink:type="simple"/></disp-formula><p>Equation (45) updates the estimate via<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x184.png" xlink:type="simple"/></inline-formula>, a measured value, and Equation (46) updates the error covariance.</p><disp-formula id="scirp.77747-formula47"><label>(45)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x185.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.77747-formula48"><label>(46)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x186.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-3400502x187.png" xlink:type="simple"/></inline-formula>in Equation (47) is the optimal estimate of the second order shadow process.</p><disp-formula id="scirp.77747-formula49"><label>(47)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-3400502x188.png"  xlink:type="simple"/></disp-formula><p>It is assumed that channel variation is mainly due to the changing mobile velocity and the correlation distance. Therefore, only the variation of the shadow process coefficient is considered. The smaller the sample period, the closer the shadow process coefficient is to one.</p><p>When the channel is nonlinear, the Unscented Kalman Filter also can be applied to the state space model optimize the shadow power presented in this section. In our problem, distribution of multipath is non-Gaussian. However, even when the white Gaussian noise assumption is not valid, the linear Kalman filter is still the optimal LMMSE estimator if the driving and measurement noises are white [<xref ref-type="bibr" rid="scirp.77747-ref5">5</xref>] .</p><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Block diagram of wireless communication system used in the outdoor experi- ment</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-3400502x189.png"/></fig></sec></sec><sec id="s5"><title>5. Outdoor Experiments</title><p>In [<xref ref-type="bibr" rid="scirp.77747-ref1">1</xref>] , we presented a new technique to experimentally characterize the wireless small-scale fading channel, taking into consideration real environmental conditions. Then, using the technique that we described in [<xref ref-type="bibr" rid="scirp.77747-ref1">1</xref>] , laboratory experiments were performed to collect real data and to validate the simulation model for the second-order KF-based estimator (designed to estimate power signal in cell phones). In [<xref ref-type="bibr" rid="scirp.77747-ref2">2</xref>] , we showed that our simulation results in the MATLAB environment and laboratory experiment results validate the proposed algorithm and the theoretical analysis. We concluded that the second-order KF-based estimator is more accurate in predicting local shadow power profiles than the first-order KF-based estimator, even in channels with imposed non-Gaussian measurement noise. Next, we wanted to see how robust a second-order KF-based estimator is in the outdoor environment where we have different large-scale fading configurations. Therefore, outdoor experiments were set up as shown in <xref ref-type="fig" rid="fig3">Figure 3</xref> to test a second-order KF-based model designed to estimate shadow power where wireless communication takes place while the user is in motion.</p><sec id="s5_1"><title>5.1. Measurements</title><p>A cell phone or portable phone uses radio waves to establish connection with its base station. Radio waves can travel long distances, but they easily get interrupted. As the transmitted signals travel from tower station to cell phone, they penetrate the atmosphere, and some signals are scattered, reflected, or observed. Objects obstructing the propagation path between the transmitter and receiver can cause variations in the received signal. All this can have a significant impact on signal strength in the cell phone device.</p><p>In this experiment, a mobile phone signal refers to signal strength received by a mobile antenna from a cellular network. There are several ways to measure mo- bile signal strength. The two most common units of measurement used in radio signals are dBm (decibels) and RSSI (Received Signal Strength Indicator). RSSI is</p><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Cell phone signal strength range</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-3400502x190.png"/></fig><p>a measurement of the power present in a received radio signal. The higher the RSSI number, the stronger the signal. These values allow users to know when they are receiving a stronger signal or a weaker signal. <xref ref-type="fig" rid="fig4">Figure 4</xref> shows that decibel values are expressed as negative numbers, which implies that the closer the value is to zero, the stronger the received signal. In Michigan, the standard frequency range for cellular phone operation is between −50 dB to −120 dB. User will get the best signal at −50 dB as it is considered full strength. On the other hand, −120 dB is considered a dead zone and the user will have no phone service. It also indicates that an ideal signal strength for optimum performance of a cell phone is about −65 dBm.</p></sec><sec id="s5_2"><title>5.2. Field Test Scenario</title><p>Area or a region can impact signal strength or path loss. Therefore, as part of the experiment, we collected data in suburban and urban areas. Measurements have been conducted in two different environments while the user was driving a vehicle at different speeds:</p><p>1) Suburban environment, Oakland University campus in Auburn Hills, Michigan.</p><p>2) Urban environment, downtown Detroit, Michigan.</p></sec><sec id="s5_3"><title>5.3. Outdoor Experiment Prerequisites and Setup</title><p>Valid research experiment must meet the certain criteria. To satisfy terrain requirement, an experiment must be conducted in an area that has good wireless cell phone coverage. Presence of large obstructing objects such as tall building structures are essential for forming a fading channel. Finally, equipment required for power signal data acquisition and processing include: base station, mobile station, vehicle, telecom tool, and laptop with built-in Bluetooth model and MATLAB software.</p><p>Cell phones work by communicating via radio waves using a system of cell towers that send and receive calls. A base station, also known as a cell tower, is placed on a big metal pole about 300 ft. high. Cell towers have triangular platforms on the top of the pole for cellular providers to keep their equipment. The process of a cell phone tower transmission requires the following equipment: radios, antennas for receiving and transmitting radio frequency signals, computerized switching control equipment, GPS receivers, power sources, and some kind of protective cover. In this experiment, Verizon was the cellular provider and the location of the base station is shown in <xref ref-type="fig" rid="fig5">Figure 5</xref>.</p><p>A mobile station consists of the physical equipment (radio transceiver, display and digital signal processors) and software package needed for communication with a mobile network. In this experiment, the Samsung Galaxy S5 smartphone was used as a mobile station. Any type of legal vehicle is acceptable to perform a driving test on public roads. As cell phone user moves around while using a cell phone, tall buildings will shadow the radio signal, which can result in a power drop at the receiver input. In this research, initial experiments were performed next to large buildings on the Oakland University campus to create a shadow fading phenomena in the outdoor environment. Supplementary experiments were conducted in downtown Detroit.</p><p>Telecom tool that was released to the market by Wylisis in March 2017 is recommended for recording captured data (<xref ref-type="fig" rid="fig6">Figure 6</xref>) and vehicle movement. Telecom software has the capability to save logging data, which can be imported into MATLAB. Measured data includes cell tower location markers, signal strength, position, velocity, and time. Alternatively, the Data Acquisition Toolbox provides functions for connecting MATLAB to data acquisition hardware. Data can be analyzed as it is acquired or it can be saved for post-processing. Block diagram of wireless communication system used in the experiment setup outside of the lab environment is illustrated in <xref ref-type="fig" rid="fig3">Figure 3</xref>. The algorithm flow chart is presented in <xref ref-type="fig" rid="fig7">Figure 7</xref>.</p><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Showing distance between mobile station and base station during outdoor experiment at Oakland University</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-3400502x191.png"/></fig><fig-group id="fig6"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Signal Strength measurement during outdoor experiment.</title></caption><fig id ="fig6_1"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-3400502x192.png"/></fig><fig id ="fig6_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-3400502x193.png"/></fig></fig-group><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> The algorithm flow chart</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-3400502x194.png"/></fig></sec><sec id="s5_4"><title>5.4. Outdoor Experiment Results</title><p>This section presents outdoor experiment results for shadow process estimation and pertinent performance analysis. The purpose of these experiments was to study and analyze output results of the first-order state space model and to compare them to the second-order state space model while applying a Kalman Filter technique to determine shadow power signal in mobile communications from measurements that have impinged Rayleigh fast fading noise. As stated before, we were able to validate this concept through laboratory experiments with data from real scenarios, but those experiments performed in the indoor environment were limited by lower speed and obstacle contribution. The outdoor experiment allowed us to conduct tests that include higher mobile velocity, exact shadow variance values, and large-scale fading configurations.</p><p>Measurements have been conducted outside while the cell phone user was driving a vehicle at different speeds, which caused variation in default parameters, such as mobile receiver velocity, shadow variance, and effective correlation distance. Multiple experiment trials were performed to collect sufficient amount of data, but in this paper we include results from driving the vehicle at 36 mph in urban area as shown in <xref ref-type="fig" rid="fig7">Figure 7</xref>. The second-order KF-based estimator performed equally well when we varied the parameters.</p><p>The plots of outdoor experiment results supported by the field data are shown in Figures 8-10. These plots show results of the actual shadow power signal and estimations with Kalman Filtering. In <xref ref-type="fig" rid="fig8">Figure 8</xref> and <xref ref-type="fig" rid="fig9">Figure 9</xref>, predicted power signal with the second-order KF-based estimator (marked in blue color) is very close to the measured signal (marked in black color). However, <xref ref-type="fig" rid="fig1">Figure 1</xref>0 shows noticeable disparity between measured signal (marked in black color) and first- order KF estimate (marked in red color). These results clearly show that the second-order KF-based estimator tracks the actual shadow power more accu- rately than the first-order KF-based estimator. Average Error for the second- order KF-based estimator is lower than the first-order KF-based estimator. Also,</p><fig id="fig8"  position="float"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> Power estimation in amobile station with asecond-order KF-based estimator using real data from the outdoor experiment</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-3400502x195.png"/></fig><fig id="fig9"  position="float"><label><xref ref-type="fig" rid="fig9">Figure 9</xref></label><caption><title> Zoomed out version of <xref ref-type="fig" rid="fig8">Figure 8</xref> for better visualization</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-3400502x196.png"/></fig><fig id="fig10"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>0</label><caption><title> Comparison of system performance with the integrated second-order KF versus the same system with the first-order KF, which shows that the implementation of the second-order KF results in better estimation</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-3400502x197.png"/></fig><p>the second-order Kalman filter output has less lag from the actual shadow power.</p><p>Authors in [<xref ref-type="bibr" rid="scirp.77747-ref5">5</xref>] presented results demonstrating that the first-order Kalman Filter method is superior to conventional window-based estimators like the sample average estimator, the uniformly minimum variance unbiased estimator, and the maximum likelihood estimator. Our results show that the second-order KF-based estimator improves the signal estimate significantly over the first-or- der KF estimate.</p></sec></sec><sec id="s6"><title>6. Conclusions</title><p>In this work, a second-order KF-based estimator has been further investigated in the outdoor environment, which is able to estimate local mean shadow power in mobile communications corrupted by multipath noise. In our experiments, we mainly explored how the second-order KF-based estimator compares to the first-order KF-based estimator. Based on our results from the indoor experiments of small-scale fading presented in [<xref ref-type="bibr" rid="scirp.77747-ref2">2</xref>] , we concluded that the second-order KF-based estimator is more accurate in predicting local mean (shadow) power profiles than the first-order KF-based estimator, even in channels with imposed non-Gaussian measurement noise. To fully complete the proposed theory, we recently extended the research to the outdoor environment and compared how these two estimators handled variability due to higher vehicle speed, larger distances, and distant large objects in the outdoor environment, such as mountains or large buildings. A Telecom tool/software released to market in 2017 was used to measure cell phone signal strength and other key parameters outside of the lab environment. Signal measurements have been conducted in typical environments like urban and suburban areas.</p><p>In this paper, we presented results from outdoor experiments which further confirmed validation of the proposed method and the theoretical analysis. The results supported by field data are provided in Figures 8-10. These plots clearly show that the second-order Kalman filter tracks the actual shadow power more accurately than the first-order Kalman filter. The system was able to operate without a failure under variety of conditions, which demonstrates model robustness. With MATLAB software, we were able to efficiently explore, analyze, and visualize measured data from the outdoor experiment. Comparison analysis was performed as explained in [<xref ref-type="bibr" rid="scirp.77747-ref11">11</xref>] [<xref ref-type="bibr" rid="scirp.77747-ref12">12</xref>] . Simulations in the MATLAB environment, laboratory, and outdoor experiment results have been consistent in showing that our implementation of the second-order KF results in better estimation.</p></sec><sec id="s7"><title>7. Future Work</title><p>Math Works currently offers some basic examples of Kalman Filter theory. Therefore, we will most likely share our code for a first-order KF-based estimator and second-order KF-based estimator by deploying an Application with MATLAB, so others can use it too. According to MathWorks’ web site, there is a wide range of options for deploying and sharing an application that was developed in MATLAB. As future work, we will look into these options.</p><p>When the channel is nonlinear, the Unscented Kalman filter also can be applied to the state space model to further improve and optimize the shadow power presented in this paper. The Unscented Kalman filter is popular due to its superiority in approximating and estimating nonlinear systems and its ability to handle non-Gaussian noise environments [<xref ref-type="bibr" rid="scirp.77747-ref13">13</xref>] . We may consider this optimization in the future.</p><p>As future work, we also are considering designing a third-order KF-based estimator. When the order of the filter is higher, we predict that there will be better noise repair. However, there is a tradeoff between three things: order of filter, computational difficulty of filter, and accuracy. Therefore, we need to look at these to determine if higher order estimators are practical.</p></sec><sec id="s8"><title>Acknowledgements</title><p>The authors would like to express appreciation to Eric Yaharmatter from Autoliv Inc. for his initial thoughts on this subject.</p></sec><sec id="s9"><title>Conflict of Interest</title><p>The authors declare that they have no conflict of interest.</p></sec><sec id="s10"><title>Cite this paper</title><p>Kapetanovic, A., Zohdy, M.A. and Mawari, R. (2017) Fading Channels Parametric Data Simulation Su- pported by Real Data from Outdoor Experiments. Journal of Signal and Information Processing, 8, 113-131. https://doi.org/10.4236/jsip.2017.83008</p></sec></body><back><ref-list><title>References</title><ref id="scirp.77747-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Mawari, R., Henderson, A., Akbar, M., Dargin, G. and Zohdy, M. (2016) ) An Improved Characterization of Small Scale Fading Based on 2D Measurements and Modeling of a Moving Receiver in an Indoor Environment. Journal of Signal and Information Processing, 7, 160-174. https://doi.org/10.4236/jsip.2016.73016</mixed-citation></ref><ref id="scirp.77747-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Kapetanovic, A., Mawari, R. and Zohdy, M. (2016) Second-Order Kalman Filtering Application to Fading Channels Supported by Real Data. Journal of Signal and Information Processing, 7, 61-74. https://doi.org/10.4236/jsip.2016.72008</mixed-citation></ref><ref id="scirp.77747-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Graf, Z. (1974) Dictionary of Electronics.</mixed-citation></ref><ref id="scirp.77747-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Yarhmatter, E. and Kapetanovic, A. (2012) Power Estimation in Mobile Communications: Comparison of the First Order AR Model to Second Order AR Model. Oakland University, Rochester, MI. (Unpublished)</mixed-citation></ref><ref id="scirp.77747-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Jiang, T., Sidiropoulos, N.D. and Giannakis, G.B. (2003) Kalman Filtering for Power Estimation in Mobile Communications, IEEE Transactions on Wireless Communications, 2,151-161. https://doi.org/10.1109/TWC.2002.806386</mixed-citation></ref><ref id="scirp.77747-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Kalman, R.E. (1960) A New Approach to Linear Filtering and Prediction Problems. Research Institute for Advanced Study, Baltimore.</mixed-citation></ref><ref id="scirp.77747-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Simon, D. (2006) Optimal State Estimation: Kalman, H∞ and Nonlinear Approaches. 1st Edition, John Wiley &amp; Sons Inc., New Jersey.  
https://doi.org/10.1002/0470045345</mixed-citation></ref><ref id="scirp.77747-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Nsour, A., Abdallah, A.-S. and Zohdy, M. (2013) An Investigation into Using Kalman Filtering for Phase Estimation in Bluetooth Receivers for Gaussian and Non-Gaussian Noise. 2013 IEEE International Conference on Electro/Information Technology, Rapid City, 9-11 May 2013, 1-5.  
https://doi.org/10.1109/EIT.2013.6632644</mixed-citation></ref><ref id="scirp.77747-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Rappaport, T.S. (2010) Wireless Communications Principles and Practice. 2nd Edition, Persons Education, London.</mixed-citation></ref><ref id="scirp.77747-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Gudmundson, M. (1991) Correlation Model for Shadowing Fading in Mobile Radio Systems. Electronics Letters, 27, 2145-2146. https://doi.org/10.1049/el:19911328</mixed-citation></ref><ref id="scirp.77747-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Grewal, M.S. and Andrews, A.P. (2014) Kalman Filtering Theory and Practice Using MATLAB. 4th Edition, John Wiley &amp; Sons Inc., New York.</mixed-citation></ref><ref id="scirp.77747-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Brown, R.G. and Hwang, P.Y.C. (2012) Introduction to Random Signals and Applied Kalman Filtering with Matlab Exercises. 4th Edition, John Wiley &amp; Sons Inc., Hoboken.</mixed-citation></ref><ref id="scirp.77747-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Li, L. and Xia, Y. (2013) Unscented Kalman Filter over Unreliable Communication Networks with Markovian Packet Dropouts. IEEE Transactions of Automatic Control, 58, 3224-3230. https://doi.org/10.1109/TAC.2013.2263650</mixed-citation></ref></ref-list></back></article>