<?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.2015.61002</article-id><article-id pub-id-type="publisher-id">JSIP-53580</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>
 
 
  Robust Non-Coherent Demodulation Scheme for Bluetooth Voice Transmission Using Linear, Extended, and Unscented Kalman Filtering
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>li</surname><given-names>S. Alghamdi</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>Mahdi</surname><given-names>N. Ali</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</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><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Hyundai Motor Company, Superior Township, MI, USA</addr-line></aff><aff id="aff1"><addr-line>Electrical and Computer Engineering Department, Oakland University, Rochester, MI, USA</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>asalgham@oakland.edu(LSA)</email>;<email>zohdyma@oakland.edu(MNA)</email>;<email>mali@hatci.com(MAZ)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>21</day><month>01</month><year>2015</year></pub-date><volume>06</volume><issue>01</issue><fpage>9</fpage><lpage>27</lpage><history><date date-type="received"><day>10</day>	<month>January</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>26</month>	<year>January</year>	</date><date date-type="accepted"><day>28</day>	<month>January</month>	<year>2015</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>
 
 
  This paper presents a novel and cost effective method to be used in the optimization of the Gaussian Frequency Shift Keying (GFSK) at the receiver of the Bluetooth communication system. The proposed method enhances the performance of the noncoherent demodulation schemes by improving the Bit Error Rate (BER) and Frame Error Rate (FER) outcomes. Linear, Extended, and Unscented Kalman Filters are utilized in this technique. A simulation model, using Simulink, has been created to simulate the Bluetooth voice transmission system with the integrated filters. Results have shown improvements in the BER and FER, and that the Unscented Kalman Filters (UKF) have shown superior performance in comparison to the linear Kalman Filter (KF) and the Extended Kalman Filter (EKF). To the best of our knowledge, this research is the first to propose the usage of the UKF in the optimization of the Bluetooth System receivers in the presence of additive white Gaussian noise (AWGN), as well as interferences.
 
</p></abstract><kwd-group><kwd>Bluetooth System Unscented Kalman Filter (UKF)</kwd><kwd> Interferences</kwd><kwd> Gaussian Frequency Shift Keying (GFSK)</kwd><kwd> Noncoherent Demodulation</kwd><kwd> Additive White Gaussian Noise (AWGN)</kwd><kwd> Bit Error Rate (BER)</kwd><kwd> Matlab Simulink</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>In recent years, we have seen that the Bluetooth technology [<xref ref-type="bibr" rid="scirp.53580-ref1">1</xref>] becomes usable in diverse applications, in addition to being a standard feature in cell-phone devices. Due to its low costs, it has been used in most of the short distance communication systems, and there is a significant reason for that: the Bluetooth technology has free wireless bandwidth for short distance connection within 2.4 GHz frequencies.</p><p>The proposed method in this research is employed at the Gaussian Frequency Shift Keying (GFSK) modulation level in the physical layer of the Bluetooth system [<xref ref-type="bibr" rid="scirp.53580-ref2">2</xref>] . Bluetooth receivers implement on-coherent demodulation, which results in low cost products, but unreliable and poor performance in most applications. On the other hand, in order to optimize a non-coherent demodulation system, some complex Maximum Likelihood Sequence Estimators (MLSE) are used. <xref ref-type="fig" rid="fig1">Figure 1</xref> shows a top level diagram of the Bluetooth system with our contribution of the robust Kalman Filter scheme. In the system block diagram, each block will be explained in detail in subsequent sections of this paper. We can see that the Kalman Filter (KF) is integrated in the receiver before the demodulation stage, so that it optimizes the detection processes. Furthermore, since we are dealing with a nonlinear system, the Extended Kalman Filter (EKF) and the Unscented Kalman Filter (UKF) are used in the optimization of the detected signals.</p><p>In 1960, Rudolf E. Kalman presented his idea of Kalman Filtering in [<xref ref-type="bibr" rid="scirp.53580-ref3">3</xref>] . KF is an estimator that estimates the states in linear dynamic systems by minimizing the Mean Square Error (MMSE). Kalman invented his filter based on modifying the Wiener-Kolmogorov filter, which was derived in 1940. While Wiener and Kolmogorov developed their filters based on frequency domain analysis, Kalman proposed his filter based on time domain analysis. The KF originated from Baye’s probabilistic theory, in which, if a good number of past samples are known, then the future samples can be predicted and updated based on continuously collected results. Since most practical systems in real life are nonlinear, many research methods have been proposed to modify KF to be applied to nonlinear systems. The EKF has been proposed to be used in nonlinear systems by linearizing the states and measurement variables under study. In [<xref ref-type="bibr" rid="scirp.53580-ref4">4</xref>] , the EKF has been used to optimize the received signal in the Bluetooth system in the presence of interferences, and produced better MSE in different Signal to Noise Ratio (SNR).</p><p>Literature [<xref ref-type="bibr" rid="scirp.53580-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.53580-ref5">5</xref>] has presented solid improvements, in which the KF proved to be useful in optimizing the Bluetooth system; however, the research in [<xref ref-type="bibr" rid="scirp.53580-ref5">5</xref>] investigated the performance of the system with the UKF. It concluded that the UKF produced a better outcome than the EKF. However, that paper performed the analysis without considering the interferences in the 2.4 GHz band in the AWGN. The work of this paper will continue on their work, which was listed in [<xref ref-type="bibr" rid="scirp.53580-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.53580-ref5">5</xref>] . In this research, the KF, EKF, and UKF are integrated into the Bluetooth receivers. The mathematical theory and derivation will be presented. The simulation model and results analysis will be analyzed and compared.</p></sec><sec id="s2"><title>2. Literature Study</title><p>There has been a great deal of literature, which has focused on the performance of the Bluetooth receivers. A</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Bluetooth system with the integrated robust Kalman Filtering technique</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x5.png"/></fig><p>significant amount of research has been performed to improve the non-coherent detection in the area of signal processing. Most of the research has developed techniques to improve the synchronization performance in low complexity demodulation schemes. The following sections describe the contributions that have been done in this related work.</p><p>Since 2003, up to the present, a GFSK has been optimized by many different techniques, which resulted in BER improvements. In [<xref ref-type="bibr" rid="scirp.53580-ref6">6</xref>] , the authors investigated the performance of frequency hopped on GFSK for the Rayleigh fading channel; they showed an improvement in BER over the fading channel, which was corrupted by AWGN by the pilot assisted differential detection. Meanwhile, the authors in [<xref ref-type="bibr" rid="scirp.53580-ref7">7</xref>] proposed a new digital receiver for a BT, where they combined a zero crossing demodulation with the de-correlating matched filter, and it also was tested on the BER performance versus<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x6.png" xlink:type="simple"/></inline-formula>, the modulation index for GFSK.</p><p>The paper in [<xref ref-type="bibr" rid="scirp.53580-ref8">8</xref>] presents an improvement on the GFSK modulation method by using maximum likelihood sequence estimation (MLSE), based on a linearized model of the modulation. Furthermore, the authors were dedicated to evolving a scheme that had MLSE-oriented detection using the Viterbi equalizer considering a GFSK inherent inter symbol interference (ISI), with the channel gain sensitivity of 0 - 7 dB. In the article [<xref ref-type="bibr" rid="scirp.53580-ref9">9</xref>] , the author proposed a new receiver design, studying the sequence detection and discovering an algorithm, which they called non-coherent sequence detection (NSD). The new algorithm is unique in low complexity, as it operates distinctions caused by low-cost local oscillators. The NSD is combined with a forward error correction decoder, which yields optimal results in real Bluetooth devices. Another proceeding [<xref ref-type="bibr" rid="scirp.53580-ref10">10</xref>] developed a simulator for the Bluetooth scatternet, which had been tested in BER performance with an ad-hoc network using the ARQ scheme.</p><p>Some of these researches were proposed for the general purpose of optimization of signal detection for the GFSK modulation; hence, it can be considered for the Bluetooth systems improvement method. Digital carrier synchronization in [<xref ref-type="bibr" rid="scirp.53580-ref11">11</xref>] used a decision-aided data cancellation algorithm to compute the carrier frequency offset (CFO) for GFSK, where they measured the error difference due to the decision output between the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x7.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x8.png" xlink:type="simple"/></inline-formula> receiver channels. The CFO will be taken into account, if the orthogonal frequency division multiplexing (OFDM) multi-user’s technique is desirable, yet it is useful only for the multi-Piconet connection in a fading environment. Later on, a new optimization was published in a differential GFSK demodulator [<xref ref-type="bibr" rid="scirp.53580-ref12">12</xref>] . The author proposed the splitting of the transmitting data into a sequence of segments assigned different coefficients by setting up the first half of the sequence as linear and the other half as non-linear.</p><p>Some literature, which is the focus of our research, was published during the last four years. The article in [<xref ref-type="bibr" rid="scirp.53580-ref13">13</xref>] investigated the Bluetooth interference by enabling the Enhanced Data Rates (EDR) to run on three types of modulation techniques: Differential Quadratic Phase Shift Keying (π/4-DQPSK, 8-DPSK, and GFSK). Their result and analysis were impressive, because they applied to a real Bluetooth experimental application. The article produced significant results for the three types of modulations in the packet delay, which was stable and unaffected when the BER was 10<sup>−3</sup> or higher. On the other hand, a study had been done on a similar topic in [<xref ref-type="bibr" rid="scirp.53580-ref14">14</xref>] , in which the authors studied the interference of 802.11b WLAN, in which they developed a new Bluetooth voice packet synchronous connection oriented with a repeated transmission (SCORT) scheme. Their analysis focused on the interferences in the 2.4 GHz frequency unlicensed bandwidth, and their application was applied in a medium access control (MAC) layer. This paper focused on the future work as it is presented in [<xref ref-type="bibr" rid="scirp.53580-ref5">5</xref>] , in which the authors proposed an optimization method of the GFSK and DPSK demodulation schemes using UKF, without considering the 802.11b interference factors in their work. The authors came up with an interactive Kalman Filtering technique to optimize the non-coherent GFSK and DPSK to produce better BER results. Finally, in [<xref ref-type="bibr" rid="scirp.53580-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.53580-ref15">15</xref>] -[<xref ref-type="bibr" rid="scirp.53580-ref17">17</xref>] , the authors made decent contributions using an EKF for GFSK, in order to estimate the phase of the signal at the receiver of the Bluetooth system, while considering the 802.11b interference factors in Gaussian noise and non-Gaussian noise environments.</p></sec><sec id="s3"><title>3. GFSK Signal Description</title><p>The GFSK modulation method is widely used in digital communications, including Bluetooth systems. <xref ref-type="fig" rid="fig2">Figure 2</xref> illustrates the GFSK modulation system, and its modulated waveform can be represented as follows [<xref ref-type="bibr" rid="scirp.53580-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.53580-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.53580-ref18">18</xref>] :</p><disp-formula id="scirp.53580-formula160"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x9.png"  xlink:type="simple"/></disp-formula><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> GFSK modulation block diagram</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x10.png"/></fig><disp-formula id="scirp.53580-formula161"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x11.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.53580-formula162"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x12.png"  xlink:type="simple"/></disp-formula><p>Equation (2) is the formal GFSK signal generated by shifting the carrier within the amount of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x13.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x14.png" xlink:type="simple"/></inline-formula>, to reflect the digital information that will be transmitted. This GFSK modulation is different from the Frequency Shift Keying (FSK), because we need to avoid using signals having large spectral side lobes. The information-bearing signal frequency modulates a signal carrier whose frequency is changed continuously. This is the result of the modulated signal, which is phase-continuous, and is called the Continuous-Phase Frequency Shift Keying (CPFSK); hence, this type of signal has memory, and the phase of the carrier is constrained to be continuous [<xref ref-type="bibr" rid="scirp.53580-ref18">18</xref>] .</p><p>A complete presentation of CPFSK will start to drive the signal representation from the Pulse Amplitude Modulation (PAM):</p><disp-formula id="scirp.53580-formula163"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x15.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x16.png" xlink:type="simple"/></inline-formula> denoted the sequence of the amplitudes obtained by mapping <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x17.png" xlink:type="simple"/></inline-formula> block in<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x18.png" xlink:type="simple"/></inline-formula>, information symbols <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x19.png" xlink:type="simple"/></inline-formula> can be represented in the format of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x20.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x21.png" xlink:type="simple"/></inline-formula> is a rectangular</p><p>pulse of the amplitude<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x22.png" xlink:type="simple"/></inline-formula>, and duration <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x23.png" xlink:type="simple"/></inline-formula> seconds. The signal <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x24.png" xlink:type="simple"/></inline-formula> is used to modulate the frequency of</p><p>the carrier. Consequently, the complex baseband waveform <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x25.png" xlink:type="simple"/></inline-formula> can be expressed as:</p><disp-formula id="scirp.53580-formula164"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x26.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x27.png" xlink:type="simple"/></inline-formula> is the peak frequency deviation and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x28.png" xlink:type="simple"/></inline-formula> is the initial phase of the carrier. Therefore, by using Euler’s formula in order to extend an exponential complex function Equation (5) can be expressed as:</p><disp-formula id="scirp.53580-formula165"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x29.png"  xlink:type="simple"/></disp-formula><p>The <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x30.png" xlink:type="simple"/></inline-formula> is the time-varying phase of the carrier, which is distinct as:</p><disp-formula id="scirp.53580-formula166"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x31.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula167"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x32.png"  xlink:type="simple"/></disp-formula><p>Carrying out the integral of the phase carrier in the interval <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x33.png" xlink:type="simple"/></inline-formula> can be represented in form:</p><disp-formula id="scirp.53580-formula168"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x34.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula169"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x35.png"  xlink:type="simple"/></disp-formula><p>where:</p><disp-formula id="scirp.53580-formula170"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x36.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula171"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x37.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula172"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x38.png"  xlink:type="simple"/></disp-formula><p>The parameter <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x39.png" xlink:type="simple"/></inline-formula> is called the modulation index. However, we show the mathematical GFSK modulation proof, which has already been discussed in [<xref ref-type="bibr" rid="scirp.53580-ref18">18</xref>] . The GFSK has the property of a Gaussian filter applied to a square signal waveform. Since the square signal has either 0 s or 1 s for output, which represent bits, and when we multiply the Gaussian filter, the frequency increases when it has symbol 1, and decreases when the symbol is 0, which all describes the frequency modulation (FM) criteria. The Gaussian filter can be described as follows:</p><disp-formula id="scirp.53580-formula173"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x40.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula174"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x41.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x42.png" xlink:type="simple"/></inline-formula> is related to the filters 3 dB bandwidth<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x43.png" xlink:type="simple"/></inline-formula>.</p></sec><sec id="s4"><title>4. Noisy Channel</title><p>The employed channel is an additive white Gaussian noise (AWGN), which is widely and extensively used in wireless communication, especially with short distance connections [<xref ref-type="bibr" rid="scirp.53580-ref9">9</xref>] . The receiving signal can be represented as:</p><disp-formula id="scirp.53580-formula175"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x44.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x45.png" xlink:type="simple"/></inline-formula> is described in Equation (6), <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x46.png" xlink:type="simple"/></inline-formula>represents the 802.11b interference with 40 dB free space loss, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x47.png" xlink:type="simple"/></inline-formula> is the white Gaussian noise AWGN carrying <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x48.png" xlink:type="simple"/></inline-formula> noise power. <xref ref-type="fig" rid="fig3">Figure 3</xref> describes the channel structure.</p><sec id="s4_1"><title>802.11 Interference Model</title><p>IEEE 802.11b is an exponential model developed for a 2.4 GHz indoor channel. It is used for small coverage areas, such as a Wi-Fi network inside a building. However, the channel conditions may vary with time and location, due to the environment being enclosed by walls, which still requires a power delay profile (PDP) to characterize their channel delay and power average. The channel impulse response can be calculated through the</p><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Noisy channel structure</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x49.png"/></fig><p>output of a finite impulse response (FIR) filter; of course, each channel is modeled by a complex Gaussian random variable with its average power. In order to determine the maximum number of paths, we need to calculate the RMS delay<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x50.png" xlink:type="simple"/></inline-formula>, and sampling period<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x51.png" xlink:type="simple"/></inline-formula>, as [<xref ref-type="bibr" rid="scirp.53580-ref23">23</xref>] :</p><disp-formula id="scirp.53580-formula176"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x52.png"  xlink:type="simple"/></disp-formula><p>Let us consider that we have a Gaussian distribution<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x53.png" xlink:type="simple"/></inline-formula>, with mean and variance respectively. The</p><p>maximum excess delay, according to the specification of the 802.11 channels, is 10 times the RMS delay. So, the power of each channel can be represented as:</p><disp-formula id="scirp.53580-formula177"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x54.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x55.png" xlink:type="simple"/></inline-formula> is:</p><disp-formula id="scirp.53580-formula178"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x56.png"  xlink:type="simple"/></disp-formula><p><xref ref-type="fig" rid="fig4">Figure 4</xref> is an illustration of an example of the 802.11 indoor channel impulse response, and its power average, when<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x57.png" xlink:type="simple"/></inline-formula>, and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x58.png" xlink:type="simple"/></inline-formula>.</p></sec></sec><sec id="s5"><title>5. Non-Coherent Detection</title><p>The detection is part of the demodulation stage at the receiver side of the BT in <xref ref-type="fig" rid="fig1">Figure 1</xref>. In BT, we do not have a reliable receiver, or, in other words, the transmitter does not have a distinguished receiver, which has none of the transmitting signal property. Therefore, they are called a non-coherent detection, due to the invalid assumption of the signals<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x59.png" xlink:type="simple"/></inline-formula>, or the availability of an orthonormal basis <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x60.png" xlink:type="simple"/></inline-formula> [<xref ref-type="bibr" rid="scirp.53580-ref18">18</xref>] , caused by either a random attenuation, or a random phase shift impact, of the signals. Hence, the non-coherent GFSK is modulated in the identity, when the signal phase matter for detection, <xref ref-type="fig" rid="fig5">Figure 5</xref>, showing how the phase trajectory for each bit can be acknowledged [<xref ref-type="bibr" rid="scirp.53580-ref18">18</xref>] .</p><p>Obviously, <xref ref-type="fig" rid="fig5">Figure 5</xref> provides a proof of our statement, when we address the GFSK modulation as a product of the periodic square signal and the Gaussian filter, so the phase increases <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x61.png" xlink:type="simple"/></inline-formula> (amplitude) when (binary 1) is transmitted, and decreases <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x62.png" xlink:type="simple"/></inline-formula> when (binary 0) is transmitted. The optimal non-coherent detection of GFSK is</p><fig-group id="fig4"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> 802.11b channel impulse example with their power.</title></caption><fig id ="fig4_1"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x63.png"/></fig><fig id ="fig4_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x64.png"/></fig></fig-group><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Phase trajectory for binary GFSK</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x65.png"/></fig><p>the real convolution of the complex baseband waveform Equation (3):</p><disp-formula id="scirp.53580-formula179"><label>(20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x66.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x67.png" xlink:type="simple"/></inline-formula> is denoted as the energy per transmitted symbol, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x68.png" xlink:type="simple"/></inline-formula>is the frequency separation between signals, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x69.png" xlink:type="simple"/></inline-formula> is the pulse period. Therefore, the general BER for GFSK can be evaluated using the general form of BER for any digital communication [<xref ref-type="bibr" rid="scirp.53580-ref2">2</xref>] .</p><disp-formula id="scirp.53580-formula180"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x70.png"  xlink:type="simple"/></disp-formula><p>In this case <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x71.png" xlink:type="simple"/></inline-formula> is the threshold voltage, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x72.png" xlink:type="simple"/></inline-formula> is the test statistic, which is a basic random variable. Then,</p><p>with the probability of transmitting<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x73.png" xlink:type="simple"/></inline-formula>, and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x74.png" xlink:type="simple"/></inline-formula>, and an optimum threshold<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x74.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x75.png" xlink:type="simple"/></inline-formula>, Equation (21) yields:</p><disp-formula id="scirp.53580-formula181"><label>(22)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x76.png"  xlink:type="simple"/></disp-formula></sec><sec id="s6"><title>6. GFSK State Space Model</title><p>A GFSK module signal needs to be modified according to the KF procedure, because KF operates in two steps: a prediction and an update, so GFSK signals have a magnitude and a phase angle. For non-coherent detection, we proposed the phase estimation, in order to detect the information using the phase trajectory illustrated in <xref ref-type="fig" rid="fig5">Figure 5</xref>. Prediction and update processes have the relation between the previous phase and current phase estimations, finding the error among them, and then predicting the next phase.</p>GFSK Description<p>In the KF development procedure, the implementation took place to evaluate or estimate in two steps: predict and update. The prediction has been computed by KF to evaluate the current phase and the previous state phase, finding the error, to update the future state, using the residual calculated error. GFSKs have a magnitude and a phase angle, which are unpredictable, due to the noise effect. Our KF are involved in two states to estimate the phase<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x77.png" xlink:type="simple"/></inline-formula>, and the difference equation<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x77.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x78.png" xlink:type="simple"/></inline-formula>, where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x77.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x78.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x79.png" xlink:type="simple"/></inline-formula> input data is a modulated phase [<xref ref-type="bibr" rid="scirp.53580-ref5">5</xref>] . The dynamic system model is represented as follows:</p><disp-formula id="scirp.53580-formula182"><label>(23)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x80.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula183"><label>(24)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x81.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x82.png" xlink:type="simple"/></inline-formula>in Equation (16) (25)</p><p>Now, our desirable result is to estimate the phase<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x83.png" xlink:type="simple"/></inline-formula>, and that can be computed using Euler’s identity for complex exponential functions. Therefore, the difference in the phase <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x84.png" xlink:type="simple"/></inline-formula> is evaluated by summing the current phase <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x85.png" xlink:type="simple"/></inline-formula> with the previous change in the phase, where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x86.png" xlink:type="simple"/></inline-formula> is the GFSK modulation index, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x86.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x87.png" xlink:type="simple"/></inline-formula> is defined in Equation (24) summation of AWGN and 802.11b interference noise.</p><disp-formula id="scirp.53580-formula184"><label>(26)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x88.png"  xlink:type="simple"/></disp-formula><p>Equation (23) shows that the system begins the module in nonlinear format, and the receiving signal can be written as:</p><disp-formula id="scirp.53580-formula185"><label>(27)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x89.png"  xlink:type="simple"/></disp-formula></sec><sec id="s7"><title>7. State Model</title><p>By substituting Equation (26) into Equations (23) and (24), the state difference equations can be formalized as:</p><disp-formula id="scirp.53580-formula186"><label>(28)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x90.png"  xlink:type="simple"/></disp-formula>Noise Measurement Model<p>According to the non-coherent detection, the typical receiver of AWGN has a correlated process when the receiving signal is divided in two phases: in-phase and quadratic phase components, which are the more familiar phase <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x91.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x91.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x92.png" xlink:type="simple"/></inline-formula> receiver channel when<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x91.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x93.png" xlink:type="simple"/></inline-formula>.</p><disp-formula id="scirp.53580-formula187"><label>(29)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x94.png"  xlink:type="simple"/></disp-formula></sec><sec id="s8"><title>8. Proposed Scheme Analysis</title><p>Referring to <xref ref-type="fig" rid="fig1">Figure 1</xref>, the first stage at the receiver, which we have called the “Robust Kalman Filter Scheme”, yields our contributed scheme, yet we are using a KF, which is one of the adaptive filters for signal processing manner. Although in this paper we adjust our results in three types of filters: linear KF, EKF, and UKF, all of them have been extensively used in stochastic random signal processing. The result could be a rich optimization for a wireless communication connection. Moreover, the author in [<xref ref-type="bibr" rid="scirp.53580-ref5">5</xref>] implemented an interactive KF, using a two-modulation technique, and simultaneously concluded his article by proving the UKF had a better result, but his effort was in the system with no interference, assuming the BT system was operated in an area such as a vehicle, where there is not much 802.11b interference. Nonetheless, the authors in [<xref ref-type="bibr" rid="scirp.53580-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.53580-ref15">15</xref>] -[<xref ref-type="bibr" rid="scirp.53580-ref17">17</xref>] have also contributed meaningfully to the same topic, to improve the non-coherent detection for GFSK proposed in their work for the BT system, by involving 802.11b interference in their results. Hence, their efforts were published, concluding with using the EKF for the purpose of a nonlinear optimization minimizing the mean square error (MES).</p><p>Beginning with our scheme description, we are designing three types of stochastic filters, as we mentioned earlier: linear KF, EKF, and UKF, while we have fully considered the latest literature review, as applied to the BT voice transmission system [<xref ref-type="bibr" rid="scirp.53580-ref19">19</xref>] . Consequently, we will be presenting our schemes analysis in three sections, each one discussing the details of our filters, and <xref ref-type="fig" rid="fig6">Figure 6</xref> gives an image of our robust scheme that cascaded with the demodulation stage.</p><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Slave non-coherent demodulation black diagram</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x95.png"/></fig><p>In the early 1960s, R. E. Kalman considered the same problem that Wiener had dealt with. In his 1960 paper, he proposed the noise measurement to be a discrete sequence in time, applied to a continuous-time signal [<xref ref-type="bibr" rid="scirp.53580-ref20">20</xref>] . In addition, the KF recursive equation was treated as a computational loop between an input and an output of the dynamic system, so the loop had two major tasks: predict and update. The following is the general form of the linear KF equations [<xref ref-type="bibr" rid="scirp.53580-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.53580-ref21">21</xref>] :</p><disp-formula id="scirp.53580-formula188"><label>(30)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x96.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula189"><label>(31)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x97.png"  xlink:type="simple"/></disp-formula><p>where:</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x98.png" xlink:type="simple"/></inline-formula>.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x99.png" xlink:type="simple"/></inline-formula>.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x100.png" xlink:type="simple"/></inline-formula>.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x101.png" xlink:type="simple"/></inline-formula>.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x102.png" xlink:type="simple"/></inline-formula>.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x103.png" xlink:type="simple"/></inline-formula>.</p><p>Now, for our system, the general Kalman state equations need to modify the base on the signal representation of GFSK, so we let<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x104.png" xlink:type="simple"/></inline-formula>. We assume the initial estimation of some point in time, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x105.png" xlink:type="simple"/></inline-formula>, and that is based on our knowledge of the process prior to<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x106.png" xlink:type="simple"/></inline-formula>. In the process of prediction, the phase <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x107.png" xlink:type="simple"/></inline-formula> takes the difference from the previous phase angle, and now we can drive from the Equation (24):</p><disp-formula id="scirp.53580-formula190"><label>(32)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x108.png"  xlink:type="simple"/></disp-formula><p>Nonetheless, we need to have<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x109.png" xlink:type="simple"/></inline-formula>, so then we can compute<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x110.png" xlink:type="simple"/></inline-formula>, using the joint estimation property<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x111.png" xlink:type="simple"/></inline-formula>; the Equation (32) can be written as:</p><disp-formula id="scirp.53580-formula191"><label>(33)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x112.png"  xlink:type="simple"/></disp-formula><p>The prior (or a priori) estimate will be denoted as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x113.png" xlink:type="simple"/></inline-formula> where the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x113.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x114.png" xlink:type="simple"/></inline-formula> denotes the estimate, and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x113.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x114.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x115.png" xlink:type="simple"/></inline-formula>, with “super minus”, is denoted as our best estimate prior to assimilating the measurement at <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x113.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x114.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x115.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x117.png" xlink:type="simple"/></inline-formula> . We can define the estimation error to be:</p><disp-formula id="scirp.53580-formula192"><label>(34)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x118.png"  xlink:type="simple"/></disp-formula><sec id="s8_1"><title>8.1. Linear Kalman Filter (KF)</title><p>First of all, we have to rewrite a KF state space model, corresponding to our GFSK modulated signal, carried in the general KF format in Equations (30) and (31).</p><disp-formula id="scirp.53580-formula193"><label>(35)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x119.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula194"><label>(36)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x120.png"  xlink:type="simple"/></disp-formula><p>For linear estimation, we have two types of covariance matrices that need to be considered for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x121.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x121.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x122.png" xlink:type="simple"/></inline-formula> vectors, which are given as [<xref ref-type="bibr" rid="scirp.53580-ref20">20</xref>] :</p><disp-formula id="scirp.53580-formula195"><label>(37)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x123.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula196"><label>(38)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x124.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula197"><label>(39)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x125.png"  xlink:type="simple"/></disp-formula><p>Corresponding to error Equation (34), the error covariance matrix is:</p><disp-formula id="scirp.53580-formula198"><label>(40)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x126.png"  xlink:type="simple"/></disp-formula><p>Now, we are looking to use the measurement <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x127.png" xlink:type="simple"/></inline-formula> to improve the prior estimate, so we assume a linear blending of the noisy measurement and the prior estimation in the equation:</p><disp-formula id="scirp.53580-formula199"><label>(41)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x128.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x129.png" xlink:type="simple"/></inline-formula> is the update estimation, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x129.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x130.png" xlink:type="simple"/></inline-formula> is the blending coefficient. Thus, by using the minimum mean square error (MMSE), we can find an optimal<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x129.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x130.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x131.png" xlink:type="simple"/></inline-formula>. In order to determine the MMSE, we need to handle the expression of the error covariance in matrix form in relation to (a posteriori) the estimate:</p><disp-formula id="scirp.53580-formula200"><label>(42)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x132.png"  xlink:type="simple"/></disp-formula><p>By substituting the Equation (36) into Equation (41), we will yield a result that has been substituted in Equa-</p><p>tion (42), and illustrate that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x133.png" xlink:type="simple"/></inline-formula> is the a priori estimation error, which is uncorrelated with measurement</p><p>error . Although, from the algebraic substitutions, we will be end up with a general expression for the update er- ror covariance matrix.</p><disp-formula id="scirp.53580-formula201"><label>(43)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x135.png"  xlink:type="simple"/></disp-formula><p>Expanding Equation (43) becomes:</p><disp-formula id="scirp.53580-formula202"><label>(44)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x136.png"  xlink:type="simple"/></disp-formula><p>Optimizing <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x137.png" xlink:type="simple"/></inline-formula> by process, two differentiations produce the trace of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x137.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x138.png" xlink:type="simple"/></inline-formula> with respect to<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x137.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x138.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x139.png" xlink:type="simple"/></inline-formula>, then solve for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x137.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x138.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x139.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x140.png" xlink:type="simple"/></inline-formula>:</p><disp-formula id="scirp.53580-formula203"><label>(45)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x141.png"  xlink:type="simple"/></disp-formula><p>Thus, the Equation (45) gives the suitable value of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x142.png" xlink:type="simple"/></inline-formula>, which, in a sense, will control minimizing the MSE estimation, and is called Kalman gain. Finalizing the computation of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x142.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x143.png" xlink:type="simple"/></inline-formula> by having a routine substitution Equation (45) in Equation (44) leads to:</p><disp-formula id="scirp.53580-formula204"><label>(46)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x144.png"  xlink:type="simple"/></disp-formula><p>Hence, we still need to find <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x145.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x145.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x146.png" xlink:type="simple"/></inline-formula>, which are our desirable best estimations; this has to be done in a similar process, due to the use of the measurement <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x145.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x146.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x147.png" xlink:type="simple"/></inline-formula> in Equation (36). From Equation (35), neglecting the term<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x145.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x146.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x147.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x148.png" xlink:type="simple"/></inline-formula>, since we have the AWGN with zero mean, and, of course, not correlated, so we have:</p><disp-formula id="scirp.53580-formula205"><label>(47)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x149.png"  xlink:type="simple"/></disp-formula><p>Back to using the error estimation Equation (34), according to Equation (44), we have:</p><disp-formula id="scirp.53580-formula206"><label>(48)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x150.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula207"><label>(49)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x151.png"  xlink:type="simple"/></disp-formula><p>So far, we derived the entire KF algorithm processes associated with the modulated BT system, and <xref ref-type="fig" rid="fig7">Figure 7</xref> shows the block diagram of the KF algorithm to prove the statement of recursive equation.</p></sec><sec id="s8_2"><title>8.2. The Extended Kalman Filter (EKF)</title><p>It was mentioned earlier, in Section 2, that [<xref ref-type="bibr" rid="scirp.53580-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.53580-ref15">15</xref>] -[<xref ref-type="bibr" rid="scirp.53580-ref17">17</xref>] have made meaningful contributions in the same subject that with which we are working, especially in using EKF for a non-coherent BT receiver. Furthermore, proposing an EKF for this article becomes more a matter of comparison, and some of the better achievements will be illustrated in Section 10. The EKF has been developed for nonlinear dynamic systems, and nonlinear measurement relationships, certainly, will require a linearization process for prediction and update measurements. Thus, the linearization will accommodate the general linear KF system module, so, starting with our system, we have an associated measurement relationship that may be written as [<xref ref-type="bibr" rid="scirp.53580-ref20">20</xref>] :</p><disp-formula id="scirp.53580-formula208"><label>(50)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x152.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula209"><label>(51)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x153.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x154.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x154.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x155.png" xlink:type="simple"/></inline-formula> are known functions, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x154.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x155.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x156.png" xlink:type="simple"/></inline-formula>is a forcing function, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x154.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x155.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x156.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x157.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x154.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x155.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x156.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x157.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x158.png" xlink:type="simple"/></inline-formula> are uncorrelated white noise processes. For linearization, we may choose <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x154.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x155.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x156.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x157.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x158.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x159.png" xlink:type="simple"/></inline-formula> to be an approximate trajectory determined by the mean.</p><disp-formula id="scirp.53580-formula210"><label>(52)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x160.png"  xlink:type="simple"/></disp-formula><p>So, Equations (50) and (51) can be written as:</p><disp-formula id="scirp.53580-formula211"><label>(53)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x161.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula212"><label>(54)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x162.png"  xlink:type="simple"/></disp-formula><p>Now, we approximate the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x163.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x163.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x164.png" xlink:type="simple"/></inline-formula> functions by using Tylor’s series expansions:</p><disp-formula id="scirp.53580-formula213"><label>(55)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x165.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula214"><label>(56)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x166.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula215"><label>(57)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x167.png"  xlink:type="simple"/></disp-formula><p>Equation (57) is the Jacobian matrix, and it allows having the nominal trajectory of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x168.png" xlink:type="simple"/></inline-formula> to satisfy the deterministic differential equation. From that, we can give the linearized equation as:</p><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Kalman recursive algorithm</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x169.png"/></fig><disp-formula id="scirp.53580-formula216"><label>(58)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x170.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula217"><label>(59)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x171.png"  xlink:type="simple"/></disp-formula><p>Finally, the EKF can apply to the linear KF algorithm shown in <xref ref-type="fig" rid="fig7">Figure 7</xref>; once the linearization process is taken care of, the sequences and the coefficients of the measurements <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x172.png" xlink:type="simple"/></inline-formula> can be evaluated as the following:</p><disp-formula id="scirp.53580-formula218"><label>(60)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x173.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula219"><label>(61)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x174.png"  xlink:type="simple"/></disp-formula></sec><sec id="s8_3"><title>8.3. The Unscented Kalman Filter (UKF)</title><p>For random variables (RV), or Gaussian random variables (GRV), the UKF proposed an extension development of EKF, in particular, for nonlinear statistical estimation regarding the main principles of linear KF, which are prediction and update. Obviously, the name of UKF came after the unscented transformation (UT), which developed the solution; in another words, the UT evaluates the mean and covariance (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x175.png" xlink:type="simple"/></inline-formula>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x175.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x176.png" xlink:type="simple"/></inline-formula>) respectively, of the RV, based on discrete samples projected through the associated nonlinear transformation [<xref ref-type="bibr" rid="scirp.53580-ref20">20</xref>] . Furthermore, the UT carries the mean and covariance of the data with the last two information samples, where these samples related to each other on the standard deviation and approximate dimension; this relation creates the sigma points from the mean, in <xref ref-type="fig" rid="fig8">Figure 8</xref>, showing the UT [<xref ref-type="bibr" rid="scirp.53580-ref21">21</xref>] [<xref ref-type="bibr" rid="scirp.53580-ref22">22</xref>] .</p><p>Assuming a nonlinear function <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x177.png" xlink:type="simple"/></inline-formula> that propagates an RV<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x177.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x178.png" xlink:type="simple"/></inline-formula>, to yield a result in another RV<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x177.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x178.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x179.png" xlink:type="simple"/></inline-formula>, so let us derive the UT from the sample point to estimate the mean and covariance in <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x177.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x178.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x179.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x180.png" xlink:type="simple"/></inline-formula> (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x177.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x178.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x179.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x180.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x181.png" xlink:type="simple"/></inline-formula>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x177.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x178.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x179.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x180.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x181.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x182.png" xlink:type="simple"/></inline-formula>) respectively.</p><p>(i) Scaling the unscented samples:</p><disp-formula id="scirp.53580-formula220"><label>(62)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x183.png"  xlink:type="simple"/></disp-formula><p>where:</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x184.png" xlink:type="simple"/></inline-formula>.</p><fig id="fig8"  position="float"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> The principle of UT</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x185.png"/></fig><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x186.png" xlink:type="simple"/></inline-formula>.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x187.png" xlink:type="simple"/></inline-formula>.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x188.png" xlink:type="simple"/></inline-formula>.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x189.png" xlink:type="simple"/></inline-formula>.</p><p>(ii) Transforming the unscented samples:<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x190.png" xlink:type="simple"/></inline-formula>.</p><p>(iii)</p><disp-formula id="scirp.53580-formula221"><label>(63)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x191.png"  xlink:type="simple"/></disp-formula><p>(iiii)</p><disp-formula id="scirp.53580-formula222"><label>(64)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x192.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula223"><label>(65)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x193.png"  xlink:type="simple"/></disp-formula><p>Equations (64) and (65) describe how UT computes the mean and covariance; also, <xref ref-type="fig" rid="fig9">Figure 9</xref> visualizes the UT data process. Now, it is necessary to run the linear KF through the output of UT, starting with the basic form of the linear KF dynamic system, involving the process update and measurement update.</p><p>The UKF algorithm starts from the original linear KF measurement model, and then takes inside an initial state where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x194.png" xlink:type="simple"/></inline-formula> for the UT processes as follows:</p><disp-formula id="scirp.53580-formula224"><label>(66)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x195.png"  xlink:type="simple"/></disp-formula><fig id="fig9"  position="float"><label><xref ref-type="fig" rid="fig9">Figure 9</xref></label><caption><title> UT mean and covariance processes</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x196.png"/></fig><disp-formula id="scirp.53580-formula225"><label>(67)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x197.png"  xlink:type="simple"/></disp-formula><p>Equations (66) and (67) are the nonlinear UKF equations with noise inputs <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x198.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x198.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x199.png" xlink:type="simple"/></inline-formula>, respectively for the BT system, which will be considered white Gaussian noise [<xref ref-type="bibr" rid="scirp.53580-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.53580-ref21">21</xref>] .</p><p>Initialize:</p><disp-formula id="scirp.53580-formula226"><label>(68)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x200.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula227"><label>(69)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x201.png"  xlink:type="simple"/></disp-formula><p>Calculating the sigma points:</p><disp-formula id="scirp.53580-formula228"><label>(70)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x202.png"  xlink:type="simple"/></disp-formula><p>Time update processes:</p><p>The prediction of UKF projects the set of samples through the nonlinear Equation (66) as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x203.png" xlink:type="simple"/></inline-formula>, and then reconstitutes the a priori state estimation. For the record, as mentioned earlier, a “super minus” script means the optimal, or the best, estimation.</p><disp-formula id="scirp.53580-formula229"><label>(71)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x205.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula230"><label>(72)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x206.png"  xlink:type="simple"/></disp-formula><p>Thus, after computing the optimal mean Equation (71), and its error covariance Equation (72), we can re- sample a new sigma point for the remainder from the following relationship:</p><disp-formula id="scirp.53580-formula231"><label>(73)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x207.png"  xlink:type="simple"/></disp-formula><p>The next steps are the measurement update processes, which required the sigma points operated through the nonlinear Equation (67).</p><disp-formula id="scirp.53580-formula232"><label>(74)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x208.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula233"><label>(75)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x209.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula234"><label>(76)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x210.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula235"><label>(77)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x211.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula236"><label>(78)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x212.png"  xlink:type="simple"/></disp-formula><p>Finally, the Kalman gain <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-3400384x213.png" xlink:type="simple"/></inline-formula> is used to update the state estimation and error covariance:</p><disp-formula id="scirp.53580-formula237"><label>(79)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x214.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.53580-formula238"><label>(80)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-3400384x215.png"  xlink:type="simple"/></disp-formula></sec></sec><sec id="s9"><title>9. Simulations and Numerical Results</title><p>The BT voice transmission system has been created on Matlab/Simulink environments [<xref ref-type="bibr" rid="scirp.53580-ref19">19</xref>] ; the demo version has been modified for the purpose of our optimization. We set the transmitting signal power to 10 mW, which is sufficient for an indoor radius of 20 meters. The high level Simulink model is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>0 for the BT system [<xref ref-type="bibr" rid="scirp.53580-ref19">19</xref>] .</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>1 shows the linear KF system model that has been implemented as a cascaded stage at the slave receiver side, and, as it is illustrated, it has two embedded functions: one for imaginary, and the other for real estimation. Furthermore, <xref ref-type="fig" rid="fig1">Figure 1</xref>2 and <xref ref-type="fig" rid="fig1">Figure 1</xref>3 will illustrate the modulation scheme for the EKF and UKF respectively, while all of the three Kalman filters have been identified in <xref ref-type="fig" rid="fig6">Figure 6</xref> as a Kalman filter subsystem block.</p><p>The simulation results have been analyzed in two different aspects: Bit Error Rate (BER) and Frame Error Rate (FER). The significance of the BER shows the performance of these types of filters in BT system application, and, as a result, it is concluded that the UKF contributes the better result between the linear KF and the EKF. Moreover, the linear KF gives the worst result, which is proof of our expectation about the nonlinear propagation of the transmitting signals, because it is only suitable for linear purposes. According to the theoretical AWGN, shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>4, the UKF gave the best result at a high noise measurement, and at a low noise measurement, versus the AWGN theoretical result.</p><p>Second, in an analysis approach based on the FER, or, in other words, on a package transmitting estimation, we divide every 366 bits into a frame, or a package. The package transmitting concept was more useful on broadcasting data transmission, and, from <xref ref-type="fig" rid="fig1">Figure 1</xref>5, where the UKF wins the computation of the three types of Kalman filters, according to expectation; the UKF will produce the best results in the matter of the FER. Nonetheless, the FER does not show a theoretical result, but the UKF in <xref ref-type="fig" rid="fig1">Figure 1</xref>4 illustrates better than do the linear KF and the EKF.</p><p>In summary, the experiment simulation is designed and run in different SNR values to come up with the previous results shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>3, and <xref ref-type="fig" rid="fig1">Figure 1</xref>4, and each result was computed out of 100 bits of error. Even</p><fig id="fig10"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>0</label><caption><title> Top level Simulink Bluetooth system model</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x216.png"/></fig><fig id="fig11"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>1</label><caption><title> Linear Kalman Filter Scheme implemented at the robust stage at the receiver of the Bluetooth</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x217.png"/></fig><fig id="fig12"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>2</label><caption><title> Extended Kalman Filter Scheme implemented at the robust stage at the receiver of the Bluetooth</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x218.png"/></fig><fig id="fig13"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>3</label><caption><title> Unscented Kalman Filter Scheme implemented at the robust stage at the receiver of the Bluetooth</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x219.png"/></fig><fig id="fig14"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>4</label><caption><title> Bit Error Rate (BER)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x220.png"/></fig><p>so, the results will be summarized in <xref ref-type="table" rid="table1">Table 1</xref>, denoting the meaningful outcomes of this experiment.</p></sec><sec id="s10"><title>10. Conclusion and Future Work</title><p>In this paper, a novel achievement has been delivered through the design of the UKF, which shows a high qual- ity for the BT system. The challenge was not only filtering a signal, but also utilizing the optimal contributions from among the related works that have been done previously. Therefore, the non-coherent GFSK has been</p><fig id="fig15"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>5</label><caption><title> Frame Error Rate (FER)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-3400384x221.png"/></fig><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Summarizing the Kalman filters performance verses signal to noise ratio (SNR)</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >SNR</th><th align="center" valign="middle"  colspan="2"  >Linear KF</th><th align="center" valign="middle"  colspan="2"  >EKF</th><th align="center" valign="middle"  colspan="2"  >UKF</th></tr></thead><tr><td align="center" valign="middle" >BER</td><td align="center" valign="middle" >FER</td><td align="center" valign="middle" >BER</td><td align="center" valign="middle" >FER</td><td align="center" valign="middle" >BER</td><td align="center" valign="middle" >FER</td></tr><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.25</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.27</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.19</td><td align="center" valign="middle" >0.97</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >0.23</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.19</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.11</td><td align="center" valign="middle" >0.70</td></tr><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle" >0.21</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.10</td><td align="center" valign="middle" >0.69</td><td align="center" valign="middle" >0.032</td><td align="center" valign="middle" >0.075</td></tr><tr><td align="center" valign="middle" >9</td><td align="center" valign="middle" >0.20</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.034</td><td align="center" valign="middle" >0.15</td><td align="center" valign="middle" >0.0065</td><td align="center" valign="middle" >0.025</td></tr><tr><td align="center" valign="middle" >12</td><td align="center" valign="middle" >0.19</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.010</td><td align="center" valign="middle" >0.15</td><td align="center" valign="middle" >0.0029</td><td align="center" valign="middle" >0.025</td></tr></tbody></table></table-wrap><p>validated using the Bluetooth voice transmission model, which had already been created in Matlab/Simulink software. As a result, an improvement has been executed on the physical layer, which is where we dealt with non-coherent GFSK modulation and demodulation; definitely, our robust scheme has been designed in order to use a faster sampling rate for the execution time of the filter.</p><p>There is always room for future works in this area, and some of them are going to be our next tasks, to continue our work. One future task will be applying the same experiment from this paper to an interactive Kalman filter, where the signals have to be modulated into two different techniques simultaneously for tunable purposes. A second task will be to develop the UKF for a multi-path fading channel, which is extensively valid for the 802.11 indoor channel models. The third task will be to consider the multi-user OFDM system, which is useful for master transmitting to Piconet and scatternet; in other words, using a multi-slave receivers connection. Finally, the experiment for the in-vehicle channel mode should be created, which would be a benefit for a tremendous number of the in-vehicle multimedia applications.</p></sec><sec id="s11"><title>Acknowledgements</title><p>I really would like to thank to Saudi Red Crescent Authority for its help in providing me the time, and the financial support, to establish my Ph.D. research.</p><p>I would also like to express my appreciation to the Saudi Arabian Cultural Mission for its infinite support, and supervising us in the United States, in order to achieve our dream.</p><p>Additionally, this research paper could not have been done without Debbie Barry, for her editing and linguistics help.</p></sec></body><back><ref-list><title>References</title><ref id="scirp.53580-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">B. S. I. Group (2014) Bluetooth SIG. Bluetooth SIG. https://www.bluetooth.org/en-us</mixed-citation></ref><ref id="scirp.53580-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Leon, I. and Couch, W. (2007) Digital and Analog Communication Systems. 7th Edition, Pearson Prentice Hall, Upper Saddle River.</mixed-citation></ref><ref id="scirp.53580-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Brown, R.G. and Hwang, P.Y. (2012) Introduction to Random Signals and Applied Kalman Filtering. 4th Edition, John Wiley &amp; Sons, Inc., Hoboken.</mixed-citation></ref><ref id="scirp.53580-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Nsour, A., Abdallah, A.-S. and Zohdy, M. (2013) Phase Tracking in Bluetooth Receivers Using Extended Kalman Filtering. 3rd International Conference on Wireless Communications, Vehicular Technology, Information Theory and Aerospace &amp; Electronic Systems (VITAE), Atlantic City, 24-27 June 2013, 1-5.</mixed-citation></ref><ref id="scirp.53580-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Ali, M. and Zohdy, M. (2012) Interactive Kalman Filtering for Differential and Gaussian Frequency Shift Keying Modulation with Application in Bluetooth. Journal of Signal and Information Processing, 3, 63-76. http://dx.doi.org/10.4236/jsip.2012.31009</mixed-citation></ref><ref id="scirp.53580-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Liu, H., Venkatesan, V., Nilsen, C., Kyker, R. and Maga&amp;ntilde;a, M. (2003) Performance of Frequency Hopped Noncoherent GFSK in Correlated Rayleigh Fading Channels. ICC’03. IEEE International Conference, 11-15 May 2003, 2779-2783.</mixed-citation></ref><ref id="scirp.53580-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Scholand, T. and Jung, P. (2003) Novel Receiver Structure for Bluetooth Based on Modified Zero-Crossing Demodulation. Global Telecommunications Conference, GLOBECOM’03. IEEE, 1-5 December 2003, 729-733.</mixed-citation></ref><ref id="scirp.53580-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Lampe, L., Schober, R. and Jain, M. (2005) Noncoherent Sequence Detection Receiver for Bluetooth Systems. Selected Areas in Communications, IEEE Journal, 23, 1718-1727. http://dx.doi.org/10.1109/JSAC.2005.853791</mixed-citation></ref><ref id="scirp.53580-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Roh, J.S. (2006) Performance Analysis and Evaluation of Bluetooth Networks in Wireless Channel Environment. Systems and Networks Communications, ICSNC’06, October 2006, Tahiti, 61 p.</mixed-citation></ref><ref id="scirp.53580-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Chang, D.-C. and Shiu, T.-H. (2006) Digital GFSK Carrier Synchronization. Circuits and Systems, APCCAS 2006, IEEE Asia Pacific, Singapore, 4-7 December 2006, 1523-1526.</mixed-citation></ref><ref id="scirp.53580-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Yu, B., Yang, L. and Chong, C.-C. (2011) Optimized Differential GFSK Demodulator. IEEE Transactions on Communications, 59, 1497-1501. http://dx.doi.org/10.1109/TCOMM.2011.041111.100010A</mixed-citation></ref><ref id="scirp.53580-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Luque, J.R., Morón, M.J. and Casilari, E. (2012) Analytical and Empirical Evaluation of the Impact of Gaussian Noise on the Modulations Employed by Bluetooth Enhanced Data Rates. EURASIP Journal on Wireless Communications and Networking, 2012, 94. http://dx.doi.org/10.1186/1687-1499-2012-94</mixed-citation></ref><ref id="scirp.53580-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Mathew, A., Chandrababu, N., Elleithy, K. and Rizvi, S. (2010) Interference of 802.11B WLAN and Bluetooth: Analysis and Performance Evaluation. International Journal of Computer Networks &amp; Communications (IJCNC), 2, 140-150.</mixed-citation></ref><ref id="scirp.53580-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">John, P.G. and Masoud, S. (2008) Digital Communications. McGraw-Hill, New York.</mixed-citation></ref><ref id="scirp.53580-ref15"><label>15</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. Proceedings of the 2013 IEEE International Conference on Electro/Information Technology (EIT), Rapid City, 9-11 May 2013, 1-5.</mixed-citation></ref><ref id="scirp.53580-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Nsour, A., Abdallah, A.-S. and Zohdy, M. (2013) GFSK Phase Estimation Using Extended Kalman Filtering for Non-Gaussian Noise. Proceedings of the Wireless Telecommunications Symposium (WTS), Phoenix, 17-19 April 2013, 1-6.</mixed-citation></ref><ref id="scirp.53580-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Nsour, A., Abdallah, A.-S. and Zohdy, M. (2014) Non-Coherent Detection of GFSK Using Extended Kalman Filtering for Non-Gaussian Noise. Proceedings of the Wireless Telecommunications Sym-posium (WTS), Washington DC, 9-11 April 2014, 1-6.</mixed-citation></ref><ref id="scirp.53580-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Mathworks, Inc. (2014) Bluetooth&amp;reg; Voice Transmission.http://www.mathworks.com/help/comm/examples/bluetooth-voice-transmission.html</mixed-citation></ref><ref id="scirp.53580-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Wan, E. and Van der Merwe, R. (2000) The Unscented Kalman Filter for Nonlinear Estimation. Proceedings of the IEEE 2000 Adaptive Systems for Signal Processing, Communications, and Control Symposium, Lake Louise, 1-4 October 2000, 153-158.</mixed-citation></ref><ref id="scirp.53580-ref20"><label>20</label><mixed-citation publication-type="book" xlink:type="simple">Brown, R.G. and Hwang, P.Y.C. (2012) Introduction to Random Signals and Applied Kalman Filtering. In: Brown, R.G. and Hwang, P.Y.C., Eds., Introduction to Random Signals and Applied Kalman Filtering with Matlab Exercises, 4th Edition, John Wiley &amp; Sons, Inc., Hoboken, 141-148.</mixed-citation></ref><ref id="scirp.53580-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Yang, W.Y., Cho, Y.S., Choo, C.Y., Shim, J., Jeon, W.G., Kim, J., Yu, S., Park, K.W. and Wee, J. (2012) MATLAB/Simulink for Digital Signal Processing. Hongrung Publishing Company, Seoul.</mixed-citation></ref><ref id="scirp.53580-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Julier, S. and Uhlmann, J. (2004) Unscented Filtering and Nonlinear Estimation. Proceedings of the IEEE, 92, 401- 422.</mixed-citation></ref><ref id="scirp.53580-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Cho, Y.S., Kim, J., Yang, W.Y. and Kang, C.G. (2010) MIMO-OFDM Wireless Communications with MATLAB. John Wiley&amp; Sons (Asia), Singapore.</mixed-citation></ref></ref-list></back></article>