<?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">JCC</journal-id><journal-title-group><journal-title>Journal of Computer and Communications</journal-title></journal-title-group><issn pub-type="epub">2327-5219</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jcc.2016.43007</article-id><article-id pub-id-type="publisher-id">JCC-64114</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>
 
 
  Joint Noise Reduction and lp-Norm Minimization for Enhancing Time Delay Estimation in Colored Noise
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Jingxian</surname><given-names>Tu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Youshen</surname><given-names>Xia</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>College of Mathematics and Computer Science, Fuzhou University, Fuzhou, China</addr-line></aff><pub-date pub-type="epub"><day>02</day><month>03</month><year>2016</year></pub-date><volume>04</volume><issue>03</issue><fpage>46</fpage><lpage>53</lpage><history><date date-type="received"><day>24</day>	<month>December</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>26</month>	<year>February</year>	</date><date date-type="accepted"><day>2</day>	<month>March</month>	<year>2016</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><html>
 <head></head>
 
   Time delay estimation (TDE) is an important issue in signal processing. Conventional TDE algorithms are usually efficient under white noise environments. In this paper, a joint noise reduction and <img src="Edit_0963a419-05bb-4c51-85d5-1d51f3ad1fed.bmp" alt="" />
   
    <!--[if gte mso 9]><xml>
 <o:oleobject type="Embed" progid="Equation.DSMT4" shapeid="_x0000_i1025" drawaspect="Content" objectid="_1518423965">
 </o:oleobject>
</xml><![endif]-->-norm minimization method is presented to enhance TDE in colored noise. An improved subspace method for colored noise reduction is first performed. Then the time delay is estimated by using an <img src="Edit_0e0ae937-e412-4096-afde-1cdffa089129.bmp" alt="" />
 
    <!--[if gte mso 9]><xml>
 <o:oleobject type="Embed" progid="Equation.DSMT4" shapeid="_x0000_i1026" drawaspect="Content" objectid="_1518423966">
 </o:oleobject>
</xml><![endif]-->-norm minimization method. Because the clean speech signal form the noisy signal is well extracted by noise reduction and the <img src="Edit_46eea1a8-3730-4f8e-bcae-60c1a58153ac.bmp" alt="" />
  
    <!--[if gte mso 9]><xml>
 <o:oleobject type="Embed" progid="Equation.DSMT4" shapeid="_x0000_i1027" drawaspect="Content" objectid="_1518423967">
 </o:oleobject>
</xml><![endif]-->-norm minimization method is robust, the TDE accuracy can be enhanced. Experiment results confirm that the proposed joint estimation method obtains more accurate TDE than several conventional algorithms in colored noise, especially in the case of low signal-to-noise ratio.  
     
 
</html></p></abstract><kwd-group><kwd>Time Delay Estimation</kwd><kwd> Speech Enhancement</kwd><kwd> Noise Reduction</kwd><kwd> Subspace</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Estimating the time delay from two received signals at spatially separated sensors is of important significance in signal processing [<xref ref-type="bibr" rid="scirp.64114-ref1">1</xref>]. It has many practical applications, such as multichannel speech enhancement, echo cancellation and wireless communications. The basic problem of the time delay estimation (TDE) is estimating accurately the time delay of interfering signals, aiming to exclude the influence of noise and interference. Many TDE approaches have been proposed. They mainly include the generalized correlation method [<xref ref-type="bibr" rid="scirp.64114-ref2">2</xref>], the statistical method [<xref ref-type="bibr" rid="scirp.64114-ref3">3</xref>], the parametric estimation method [<xref ref-type="bibr" rid="scirp.64114-ref4">4</xref>], the adaptive estimation method [<xref ref-type="bibr" rid="scirp.64114-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.64114-ref6">6</xref>], the combinational estimation method [<xref ref-type="bibr" rid="scirp.64114-ref7">7</xref>], and the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x7.png" xlink:type="simple"/></inline-formula>-norm minimization-based estimation method [<xref ref-type="bibr" rid="scirp.64114-ref8">8</xref>]. Among them, the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x8.png" xlink:type="simple"/></inline-formula>-norm minimization-based estimation method can find the time delay by minimizing an <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x9.png" xlink:type="simple"/></inline-formula>-norm objective function. It was reported in [<xref ref-type="bibr" rid="scirp.64114-ref8">8</xref>] that this method can obtain more robust results than other several conventional approaches against impulse noise. However, these conventional TDE approaches do not consider the influence of noise, specially under low SNR conditions. Moreover, in general they are only efficient in white noise.</p><p>Speech enhancement techniques have been applied in speech recognition and voice communication. They can recover the clean speech signal from the noisy signal by noise reduction. Speech enhancement algorithms can be classified as single channel and multichannel speech enhancement algorithms. The multichannel speech enhance- ment algorithms usually mix multiple noisy signals for noise reduction. By contrast, the single speech enhance- ment algorithms utilize only one noisy signal and thus do not change the time delay of the noisy signal at each channel. So, the single speech enhancement algorithms have the potential to improve the performance of the TDE algorithm. At present, the single-channel speech enhancement algorithms mainly include the spectral subtraction-based methods [<xref ref-type="bibr" rid="scirp.64114-ref9">9</xref>], the Kalman filtering-based parametric method [<xref ref-type="bibr" rid="scirp.64114-ref10">10</xref>], the statistic-based approach [<xref ref-type="bibr" rid="scirp.64114-ref11">11</xref>] [<xref ref-type="bibr" rid="scirp.64114-ref12">12</xref>], and the subspace-based method [<xref ref-type="bibr" rid="scirp.64114-ref13">13</xref>]. These conventional algorithms are suitable for white noise reduction. To deal with colored noise, one conventional approach is that the noisy speech signal is multiplied by the square root of the noise covariance matrix’s inverse [<xref ref-type="bibr" rid="scirp.64114-ref14">14</xref>]. Another conventional approach is the prewhitening covariance matrix of the colored noise. These prewhitening approaches all require to estimate the covariance matrix of the colored noise in advance. Recently, to avoid disadvantage of estimating the covariance matrix of the colored noise, an improved subspace method was proposed in [<xref ref-type="bibr" rid="scirp.64114-ref15">15</xref>]. It was reported that the improved subspace method outperforms conventional speech enhancement methods in colored noise reduction.</p><p>In this paper, a new method for enhancing time delay estimation (TDE) in colored noise is presented by joining noise reduction and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x10.png" xlink:type="simple"/></inline-formula>-norm minimization. We first perform the improved subspace method for enhanced signals corrupted by colored noise and then we use the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x11.png" xlink:type="simple"/></inline-formula>-norm minimization based TDE method to estimate the time delay from the enhanced signals. Experiment results show that the proposed joint algorithm can obtain more accurate TDE than several conventional algorithms in colored noise, especially in the case of low signal-to-noise ratio.</p></sec><sec id="s2"><title>2. TDE Signal Model and Estimation</title><sec id="s2_1"><title>2.1. Signal Model</title><p>Consider the following TDE signal model:</p><disp-formula id="scirp.64114-formula56"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/64114x12.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x13.png" xlink:type="simple"/></inline-formula> is the unknown random source signal, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x14.png" xlink:type="simple"/></inline-formula>is the attenuation factor, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x15.png" xlink:type="simple"/></inline-formula>is the time delay to be estimated, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x16.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x17.png" xlink:type="simple"/></inline-formula> are uncorrelated noise observations which are independent of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x18.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x19.png" xlink:type="simple"/></inline-formula> is the sampling length of the noisy signals. The goal of TDE is to estimate the delay <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x20.png" xlink:type="simple"/></inline-formula> from <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x21.png" xlink:type="simple"/></inline-formula> noise observations.</p></sec><sec id="s2_2"><title>2.2. l<sub>p</sub>-Norm Minimization-Based Estimation</title><p>For robust TDE, Ma and Nikias introduced [<xref ref-type="bibr" rid="scirp.64114-ref16">16</xref>] the following <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x22.png" xlink:type="simple"/></inline-formula>-norm cost function about the delay D and the attenuation factor<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x23.png" xlink:type="simple"/></inline-formula>:</p><disp-formula id="scirp.64114-formula57"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/64114x24.png"  xlink:type="simple"/></disp-formula><p>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x25.png" xlink:type="simple"/></inline-formula>, M is the order parameter and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x26.png" xlink:type="simple"/></inline-formula> is the sinc function. (2) can be written in a</p><p>matrix form as:</p><disp-formula id="scirp.64114-formula58"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/64114x27.png"  xlink:type="simple"/></disp-formula><p>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x28.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x29.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x30.png" xlink:type="simple"/></inline-formula>, and</p><disp-formula id="scirp.64114-formula59"><graphic  xlink:href="http://html.scirp.org/file/64114x31.png"  xlink:type="simple"/></disp-formula><p>To minimize<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x32.png" xlink:type="simple"/></inline-formula>, Zeng et al. presented an efficient two-steps procedure [<xref ref-type="bibr" rid="scirp.64114-ref8">8</xref>]. In the first step, the global optimum <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x33.png" xlink:type="simple"/></inline-formula> is estimated for each given D. The estimation for the global optimum <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x34.png" xlink:type="simple"/></inline-formula> has the following three cases.</p><p>Case 1:<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x35.png" xlink:type="simple"/></inline-formula>. The <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x36.png" xlink:type="simple"/></inline-formula> cost function is a one-dimensional quadratic function of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x37.png" xlink:type="simple"/></inline-formula>. Its optimal solution is given by</p><disp-formula id="scirp.64114-formula60"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/64114x38.png"  xlink:type="simple"/></disp-formula><p>Case 2:<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x39.png" xlink:type="simple"/></inline-formula>. The <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x40.png" xlink:type="simple"/></inline-formula> cost function is the least absolute deviation function:</p><disp-formula id="scirp.64114-formula61"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/64114x41.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x42.png" xlink:type="simple"/></inline-formula> is the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x43.png" xlink:type="simple"/></inline-formula>th element of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x44.png" xlink:type="simple"/></inline-formula>. Let<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x45.png" xlink:type="simple"/></inline-formula>, then the optimal</p><p>solution of the cost function is the weighted median of the sequence <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x46.png" xlink:type="simple"/></inline-formula> with the weights<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x47.png" xlink:type="simple"/></inline-formula>. The procedure of computation of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x48.png" xlink:type="simple"/></inline-formula> is listed in Algorithm 1.</p><p>Case 3:<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x58.png" xlink:type="simple"/></inline-formula>. The <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x59.png" xlink:type="simple"/></inline-formula> cost function has derivative. Thus the following fixed-point iteration is used to find the optimal solution<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x60.png" xlink:type="simple"/></inline-formula>:</p><disp-formula id="scirp.64114-formula62"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/64114x61.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x62.png" xlink:type="simple"/></inline-formula> is the estimate of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x63.png" xlink:type="simple"/></inline-formula> in the kth iteration, and</p><disp-formula id="scirp.64114-formula63"><graphic  xlink:href="http://html.scirp.org/file/64114x64.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x65.png" xlink:type="simple"/></inline-formula> is the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x66.png" xlink:type="simple"/></inline-formula>th element of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x67.png" xlink:type="simple"/></inline-formula>. In the second step, a search range <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x68.png" xlink:type="simple"/></inline-formula> and a step size <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x69.png" xlink:type="simple"/></inline-formula> are first determined, then the value of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x70.png" xlink:type="simple"/></inline-formula> increases from <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x71.png" xlink:type="simple"/></inline-formula> to <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x72.png" xlink:type="simple"/></inline-formula> with the step size being<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x73.png" xlink:type="simple"/></inline-formula>, after that the delay profile is computed by substituting each given <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x74.png" xlink:type="simple"/></inline-formula> and the corresponding <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x74.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x75.png" xlink:type="simple"/></inline-formula> into</p><p>the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x76.png" xlink:type="simple"/></inline-formula> cost function, finally the minimum <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x76.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x77.png" xlink:type="simple"/></inline-formula> is used to estimate the<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x76.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x77.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x78.png" xlink:type="simple"/></inline-formula>.</p></sec></sec><sec id="s3"><title>3. Proposed TDE Algorithm</title><sec id="s3_1"><title>3.1. Improved Subspace Method for Colored Noise Reduction</title><p>Without the loss of generality, we consider the following noise signal model:</p><disp-formula id="scirp.64114-formula64"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/64114x79.png"  xlink:type="simple"/></disp-formula><p>Our goal is to restore <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x80.png" xlink:type="simple"/></inline-formula> from <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x81.png" xlink:type="simple"/></inline-formula> by colored noise reduction.</p><p>Recently, for colored noise reduction, an improved subspace method was presented in [<xref ref-type="bibr" rid="scirp.64114-ref15">15</xref>]. Let the colored noise be modeled as the pth order autoregressive signal process</p><disp-formula id="scirp.64114-formula65"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/64114x82.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x83.png" xlink:type="simple"/></inline-formula> are the AR noise model parameters, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x84.png" xlink:type="simple"/></inline-formula>is the drive noise which is assumed to be white with variance<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x85.png" xlink:type="simple"/></inline-formula>. Let K denote the length of one frame signal, and let<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x86.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x86.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x87.png" xlink:type="simple"/></inline-formula>, and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x86.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x87.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x88.png" xlink:type="simple"/></inline-formula>. Then (7) can be written in a vector form:</p><disp-formula id="scirp.64114-formula66"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/64114x89.png"  xlink:type="simple"/></disp-formula><p>and (8) can be written in a vector form:</p><disp-formula id="scirp.64114-formula67"><graphic  xlink:href="http://html.scirp.org/file/64114x90.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x91.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x91.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x92.png" xlink:type="simple"/></inline-formula> is the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x91.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x93.png" xlink:type="simple"/></inline-formula> whitening matrix:</p><disp-formula id="scirp.64114-formula68"><graphic  xlink:href="http://html.scirp.org/file/64114x94.png"  xlink:type="simple"/></disp-formula><p>Multiplying (9) by<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x95.png" xlink:type="simple"/></inline-formula>, we have</p><disp-formula id="scirp.64114-formula69"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/64114x96.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x97.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x98.png" xlink:type="simple"/></inline-formula>. Since <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x98.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x99.png" xlink:type="simple"/></inline-formula> is the white noise, the conventional subspace method for white noise reduction can be directly used to estimate<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x98.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x99.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x100.png" xlink:type="simple"/></inline-formula>, given as</p><disp-formula id="scirp.64114-formula70"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/64114x101.png"  xlink:type="simple"/></disp-formula><p>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x102.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x102.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x103.png" xlink:type="simple"/></inline-formula>is the Lagrangian multiplier, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x102.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x103.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x104.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x102.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x103.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x105.png" xlink:type="simple"/></inline-formula>is the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x102.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x103.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x106.png" xlink:type="simple"/></inline-formula></p><p>covariance matrix of the whitening signal vector<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x107.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x108.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x109.png" xlink:type="simple"/></inline-formula> consist of the eigenvector and eigenvalue of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x110.png" xlink:type="simple"/></inline-formula>, respectively. Then the clean speech signal <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x111.png" xlink:type="simple"/></inline-formula> can be estimated by</p><disp-formula id="scirp.64114-formula71"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/64114x112.png"  xlink:type="simple"/></disp-formula><p>For each signal frame, the improved subspace algorithm (denoted as ISS) is summarized as:</p></sec><sec id="s3_2"><title>3.2. Proposed TDE Algorithm</title><p>In this section, we introduce a new method for enhancing time delay estimation (TDE) in colored noise, based on joint noise reduction and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x119.png" xlink:type="simple"/></inline-formula>-norm minimization. An improved subspace method for colored noise reduction is first performed. The time delay is then estimated by using the enhanced signal, based on the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x119.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x120.png" xlink:type="simple"/></inline-formula>-norm minimization. The proposed TDE algorithm is listed in Algorithm 3. Compared with conventional TDE algorithms, the proposed TDE algorithm can greatly reduce the interference of colored noise such that the TDE accuracy is enhanced.</p><disp-formula id="scirp.64114-formula72"><graphic  xlink:href="http://html.scirp.org/file/64114x121.png"  xlink:type="simple"/></disp-formula></sec></sec><sec id="s4"><title>4. Experimental Results</title><p>In this section, we conduct numerical simulations to demonstrate the effectiveness of the proposed algorithm. We compare the proposed TDE algorithm with the TDE algorithm based <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x122.png" xlink:type="simple"/></inline-formula>-norm minimization without noise reduction. We also compare with other two TDE algorithms based <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x123.png" xlink:type="simple"/></inline-formula>-norm minimization with noise reduction, where the minimum mean square error(MMSE) and maximum a posterior(MAP) estimators of the magnitude- squared spectrum(denoted as MMSE-MSS and MAP-MSS) are used for noise reduction, respectively. The source signal and noisy signal are taken from the NOIZEUS [<xref ref-type="bibr" rid="scirp.64114-ref17">17</xref>] and NOISEX [<xref ref-type="bibr" rid="scirp.64114-ref18">18</xref>] corpora, respectively. We randomly select twenty different speech sentences from the NOIZEUS corpora. Babble and factory noises are selected from the NOISEX corpora. Each speech sentence is corrupted by these two noises with different input SNRs. The true delay is set to<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x124.png" xlink:type="simple"/></inline-formula>, the attenuation factor is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x125.png" xlink:type="simple"/></inline-formula>, the approximation order parameter is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x126.png" xlink:type="simple"/></inline-formula>, the delay search range is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x126.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x127.png" xlink:type="simple"/></inline-formula>, and the search step size is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x126.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x127.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x128.png" xlink:type="simple"/></inline-formula>. To evaluate the performance of our proposed methods, we use the root mean square error(RMSE), which is defined as:</p><disp-formula id="scirp.64114-formula73"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/64114x129.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x130.png" xlink:type="simple"/></inline-formula> is the number of speech sentences and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x130.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x131.png" xlink:type="simple"/></inline-formula> is the delay estimate of the mth speech sentence. By the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x130.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x131.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x132.png" xlink:type="simple"/></inline-formula>-norm minimization method we see that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x130.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x131.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x133.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x130.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x131.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x133.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x134.png" xlink:type="simple"/></inline-formula> are the best choice in Babble noise and factory noise, respectively.</p><p>In the first test, we perform the four algorithms for different values of the input SNRs. <xref ref-type="fig" rid="fig1">Figure 1</xref> and <xref ref-type="fig" rid="fig2">Figure 2</xref> display the RMSE results of the four algorithms with different values of the input SNRs (From 0 dB to 10 dB) in factory noise and babble noise, respectively. From the two figures, we first see that the four algorithms obtain higher value of RMSE when increasing the input SNRs. This indicates that the noise decrease the performance of TDE. Second, we see that the proposed algorithm can outperform the TDE algorithm based <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x135.png" xlink:type="simple"/></inline-formula>-norm minimization without noise reduction for all input SNRs in terms of RMSE. Third, the proposed TDE algorithm can get a lower value of RMSE than the other two TDE algorithms with noise reduction, based on the MMSE- MSS and MAP-MSS estimators, respectively. This also indicates that the proposed algorithm can obtain the best accurate TDE in terms of RMSE.</p><p>In the second test, we perform the four algorithms with the input SNRs being 5 dB via different values of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x136.png" xlink:type="simple"/></inline-formula>. <xref ref-type="fig" rid="fig3">Figure 3</xref> and <xref ref-type="fig" rid="fig4">Figure 4</xref> display their RMSE results of the four algorithms in factory noise and babble noise,</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> RMSE of TDE based on four algorithms with different input SNRs and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x138.png" xlink:type="simple"/></inline-formula> in factory noise</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/64114x137.png"/></fig><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> RMSE of TDE based on four algorithms with different input SNRs and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x140.png" xlink:type="simple"/></inline-formula> in babble noise</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/64114x139.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> RMSE of TDE based on four algorithms via different values of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x142.png" xlink:type="simple"/></inline-formula> with the input SNRs being 5 dB and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x142.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x143.png" xlink:type="simple"/></inline-formula> in factory noise</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/64114x141.png"/></fig><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> RMSE of TDE based on four algorithms via different values of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x145.png" xlink:type="simple"/></inline-formula> with the input SNRs being 5 dB and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x145.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x146.png" xlink:type="simple"/></inline-formula> in babble noise</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/64114x144.png"/></fig><p>respectively. From the two figures, we first see that the proposed algorithms outperform TDE without speech enhancement in any value of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x147.png" xlink:type="simple"/></inline-formula>. Second, the proposed TDE algorithm can get a lower value of RMSE than the other two TDE algorithms with noise reduction, based on the MMSE-MSS and MAP-MSS estimators, respec- tively. This indicates that the proposed algorithm can obtain the best accurate TDE in any value of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x147.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/64114x148.png" xlink:type="simple"/></inline-formula>.</p></sec><sec id="s5"><title>Cite this paper</title><p>Jingxian Tu,Youshen Xia, (2016) Joint Noise Reduction and lp-Norm Minimization for Enhancing Time Delay Estimation in Colored Noise. Journal of Computer and Communications,04,46-53. doi: 10.4236/jcc.2016.43007</p></sec></body><back><ref-list><title>References</title><ref id="scirp.64114-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Xing, H.Y. and Tang, J. (2008) Analysis and Survey of Algorithms for Time-Delay Estimation. Technical Acoustics, 27, 110-114.</mixed-citation></ref><ref id="scirp.64114-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Knapp, C.H. and Carter, G.C (1976) The Generalized Correlation Method for Estimation of Time Delay. IEEE Trans. Acoust. Speech Signal Process, 24, 320-327. http://dx.doi.org/10.1109/TASSP.1976.1162830</mixed-citation></ref><ref id="scirp.64114-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Torrieri, D.J. (1984) Statistical Theory of Passive Location Systems. IEEE Transactions on Aerospace and Electronic Systems, 20, 183-197. http://dx.doi.org/10.1109/TAES.1984.310439</mixed-citation></ref><ref id="scirp.64114-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Kenneth, W.K., Lui, B., Frankie, K.W. and So Chan, H.C. (2009) Accurate Time Delay Estimation Based Passive Localization. Signal Processing, 89, 1835-1838. http://dx.doi.org/10.1016/j.sigpro.2009.03.009</mixed-citation></ref><ref id="scirp.64114-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Palanisamy, P. and Kalyanasundaram, N. (2010) Sonar Target Detection by Modified Adaptive Noise Cancellation Using Correlating Filter. International Journal of Electronics, 98, 41-60.  
http://dx.doi.org/10.1080/00207217.2010.497670</mixed-citation></ref><ref id="scirp.64114-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Widrow, B. and Stearn, S.D. (1993) Adaptive Signal Processing. Prentice-Hall, Englewood Cliffs, NJ.</mixed-citation></ref><ref id="scirp.64114-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Lotf, S., Mourad, T. and Sabeur, A.A.C. (2013) Performance of Wavelet Analysis and Neural Networks for Pathological Voices Identification. International Journal of Electronics, 98, 1129-1140.</mixed-citation></ref><ref id="scirp.64114-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Zeng, W.J., So, H.C. and Abdelhak, M.Z. (2013) An   Minimization Approach to Time Delay Estimation in Impulse Noise. Digital Signal Processing, 23, 1247-1254. http://dx.doi.org/10.1016/j.dsp.2013.03.013</mixed-citation></ref><ref id="scirp.64114-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Boll, S. (1979) Suppression of Acoustic Noise in Speech Using Spectral Subtraction. IEEE Trans. Acoust., Speech, Signal Process, 27, 113-120. http://dx.doi.org/10.1109/TASSP.1979.1163209</mixed-citation></ref><ref id="scirp.64114-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">So, S. and Paliwal, K.K. (2011) Suppressing the Influence of Additive Noise on the Kalman Gain for Low Residual Noise Speech Enhancement. Speech Communication, 53, 355-378. http://dx.doi.org/10.1016/j.specom.2010.10.006</mixed-citation></ref><ref id="scirp.64114-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Lu, Y. and Loizou, P.C. (2011) Estimators of the Mag-nitude-Squared Spectrum and Mehthods for Incorporating SNR Uncertainty. IEEE Transactions on Audio, Speech and Language Processing, 19, 1123-1137.  
http://dx.doi.org/10.1109/TASL.2010.2082531</mixed-citation></ref><ref id="scirp.64114-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Martin, R. (2005) Speech Enhancement Based on Minimum Mean-Square Error Estimation and Supergaussian Priors. IEEE Transactions on Speech Audio Processing, 13, 845-856. http://dx.doi.org/10.1109/TSA.2005.851927</mixed-citation></ref><ref id="scirp.64114-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Ephraim, Y. and Van Trees, H.L. (1995) A Signal Subspace Ap-proach for Speech Enhancement. IEEE Transactions on Speech Audio Processing, 3, 251-266. http://dx.doi.org/10.1109/89.397090</mixed-citation></ref><ref id="scirp.64114-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Lev-ari, H. and Ephraim, Y. (2003) Extension of The Signal Subspace Speech Enhancement Approach to Colored Noise. IEEE Signal Processing Letters, 10, 104-106. http://dx.doi.org/10.1109/LSP.2003.808544</mixed-citation></ref><ref id="scirp.64114-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Wei, Q. and Xia, Y.S. (2013) A Novel Prewhitening Subspace Method for Enhancing Speech Corrupted by Colored Noise. Proc. IEEE CISP, 3, 1282-1286. http://dx.doi.org/10.1109/cisp.2013.6743870</mixed-citation></ref><ref id="scirp.64114-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Ma, X. and Nikias, C.L. (1996) Joint Estimation of Time Delay and Frequent Delay in Impulsive Noise Using Fractional Lower Order Statistics. IEEE Trans. Signal Process., 44, 2669-2687. http://dx.doi.org/10.1109/78.542175</mixed-citation></ref><ref id="scirp.64114-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Loizou, P.C. (2007) Speech Enhancement: Theory and Practice. CRC, Boca Raton, FL, USA.</mixed-citation></ref><ref id="scirp.64114-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Varga, A. and Steeneken, H.J.M. (1993) Assessment for Automatic Speech Recog-nition: II, NOISEX-92: A Database and An Experiment to Study The Effect of Additive Noise on Speech Recognition Systems. Speech Communication, 12, 247-251. http://dx.doi.org/10.1016/0167-6393(93)90095-3</mixed-citation></ref></ref-list></back></article>