<?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.2017.53008</article-id><article-id pub-id-type="publisher-id">JCC-74656</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>
 
 
  Parameter Estimation of Maneuvering Target Using Maximum Likelihood Estimation for MIMO Radar with Colocated Antennas
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Hai</surname><given-names>Li</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>Ling</surname><given-names>Xu</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>Zhe</surname><given-names>Zhang</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Tianjin Key Lab. for Advanced Signal Processing, Civil Aviation University of China, Tianjin, China</addr-line></aff><pub-date pub-type="epub"><day>08</day><month>03</month><year>2017</year></pub-date><volume>05</volume><issue>03</issue><fpage>69</fpage><lpage>74</lpage><history><date date-type="received"><day>March</day>	<month>2,</month>	<year>2017</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>March</month>	<year>10,</year>	</date><date date-type="accepted"><day>March</day>	<month>13,</month>	<year>2017</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   
   The estimation of target parameters in MIMO radar signal processing is one of the most important research topics. An efficient implementation of the Maximum Likelihood estimator is presented in this paper to estimate the DOA (Direction of Arrival), initial velocity and acceleration of a maneuvering target in colocated MIMO radar. The target’s DOA is estimated in the first place, then a Maximum-Likelihood (ML) estimation based on peak search is applied to a two-dimensional grids providing estimation of initial velocity and acceleration. Simulations show that the MIMO radar has a better performance in DOA estimation than the phased array radar. By means of Monte Carlo simulations, the estimation error of initial velocity and acceleration on different SNRs are calculated. The results also suggest the effectiveness of this method. 
  
 
</p></abstract><kwd-group><kwd>MIMO Radar</kwd><kwd> Parameter Estimation</kwd><kwd> Maximum Likelihood Estimation</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The Multiple-Input Multiple-Output (MIMO) Radio Transmission Technology, has the potential to enhance system capacity and improve spectral efficiency without requiring extra emission signal bandwidth or transmission power [<xref ref-type="bibr" rid="scirp.74656-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.74656-ref2">2</xref>]. The overwhelming popularity of MIMO systems draws radar researchers’ attention to engage MIMO techniques in radar systems implementation.</p><p>MIMO radar employs multiple elements at both transmitting and receiving antennas as the phased array radar does, on flip side, MIMO radar transmits orthogonal or partially correlated waveforms. According to the antenna conﬁgura- tion, MIMO radar systems can be classified into two categories: collocated (coherent) MIMO radars are those with closely spaced antennas [<xref ref-type="bibr" rid="scirp.74656-ref3">3</xref>]. While distributed MIMO radars, also known as non-coherent MIMO radars, whose antennas are placed far from each other [<xref ref-type="bibr" rid="scirp.74656-ref4">4</xref>]. The latter conﬁguration exploits the random ﬂuctuation of the targets’ Radar Cross Section (RCS). This ﬂuctuation causes spatial decorrelation between target echoes, which prevents coherent processing. It is known that spatial resolution can be improved by taking advantage of coherent processing introduced by coherent MIMO radar. Several methods have been proposed for DOA estimation offering significant estimate performance without considering Doppler effects [<xref ref-type="bibr" rid="scirp.74656-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.74656-ref6">6</xref>]. As shown in [<xref ref-type="bibr" rid="scirp.74656-ref7">7</xref>] and [<xref ref-type="bibr" rid="scirp.74656-ref8">8</xref>], parameters of moving target can be measured by the maximum likelihood estimator for distributed MIMO radar. There are many efforts have been made on joint DOA and Doppler frequency estimation [<xref ref-type="bibr" rid="scirp.74656-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.74656-ref10">10</xref>], while few studies have yet been conducted on maneuvering targets.</p><p>This paper considers parameters estimation of maneuvering target in colocated MIMO radar by developing a maximum-likelihood estimator. The MUltiple SIgnal Classification (MUSIC) is applied in the first place to estimate DOA, followed by a two-dimensional search applied to estimate the velocity and acceleration. Performance of estimation for phased array radar and MIMO radar are compared under the same condition. The root mean square errors (RMSE) of the ML estimation of target parameters are obtained by Monte Carlo simulations. Numerical results demonstrate that MIMO radar has better angle estimation performance compared with conventional phased array radar and verify the effectiveness of this method as well.</p></sec><sec id="s2"><title>2. MIMO Radar Signal Model</title><p>Let’s consider a MIMO radar system in which there are <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x2.png" xlink:type="simple"/></inline-formula> transmitting antennas and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x3.png" xlink:type="simple"/></inline-formula> receiving antennas, with inter-element spacing of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x4.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x5.png" xlink:type="simple"/></inline-formula>, respectively. Radar working wavelength is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x6.png" xlink:type="simple"/></inline-formula>.The signal emitts from the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x7.png" xlink:type="simple"/></inline-formula>-th antenna, is defined as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x8.png" xlink:type="simple"/></inline-formula>, where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x9.png" xlink:type="simple"/></inline-formula> is the length of coded pulse signal, with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x10.png" xlink:type="simple"/></inline-formula> denoting the transpose of matrices. Assuming the different transmitting signals are orthogonal, then transmit steering vector of the array towards direction <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x11.png" xlink:type="simple"/></inline-formula> can be denoted as</p><disp-formula id="scirp.74656-formula78"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/74656x12.png"  xlink:type="simple"/></disp-formula><p>and the receive steering vector is</p><disp-formula id="scirp.74656-formula79"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/74656x13.png"  xlink:type="simple"/></disp-formula><p>Assuming that the target has a reflection coefficient<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x14.png" xlink:type="simple"/></inline-formula>, moves with an initial radial velocity <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x15.png" xlink:type="simple"/></inline-formula> and an acceleration<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x16.png" xlink:type="simple"/></inline-formula>.The received signal at antenna <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x17.png" xlink:type="simple"/></inline-formula> and time index <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x18.png" xlink:type="simple"/></inline-formula> can be described by</p><disp-formula id="scirp.74656-formula80"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/74656x19.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x20.png" xlink:type="simple"/></inline-formula> represents Doppler shift caused by the maneuvering target, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x21.png" xlink:type="simple"/></inline-formula> is the noise components approximated by complex white Gaussian samples at the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x22.png" xlink:type="simple"/></inline-formula>-th receiver. From (3) the received space-time matrix can be expressed as</p><disp-formula id="scirp.74656-formula81"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/74656x23.png"  xlink:type="simple"/></disp-formula><p>The total <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x24.png" xlink:type="simple"/></inline-formula> samples can be collected into a single matrix</p><disp-formula id="scirp.74656-formula82"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/74656x25.png"  xlink:type="simple"/></disp-formula><p>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x26.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x27.png" xlink:type="simple"/></inline-formula>is an <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x28.png" xlink:type="simple"/></inline-formula> diagonal matrix,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x29.png" xlink:type="simple"/></inline-formula> ,where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x30.png" xlink:type="simple"/></inline-formula> is the pulse width, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x31.png" xlink:type="simple"/></inline-formula> is an <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x32.png" xlink:type="simple"/></inline-formula> noise matrix with a variance of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x33.png" xlink:type="simple"/></inline-formula>.</p></sec><sec id="s3"><title>3. Maximum Likelihood Estimation</title><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x34.png" xlink:type="simple"/></inline-formula> be the vector denoting the unknown parameters that need to be estimated. We use data <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x35.png" xlink:type="simple"/></inline-formula> as stated in (5) and be vectorized, then the received data vector <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x36.png" xlink:type="simple"/></inline-formula> can be obtained, which follows a Gaussian probability density function whose mean is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x37.png" xlink:type="simple"/></inline-formula>. Since a maneuvering target is discussed, the Doppler frequency <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x38.png" xlink:type="simple"/></inline-formula> involves two variables: initial velocity <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x39.png" xlink:type="simple"/></inline-formula> and acceleration<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x40.png" xlink:type="simple"/></inline-formula>.</p><p>The probability density function of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x41.png" xlink:type="simple"/></inline-formula> is</p><disp-formula id="scirp.74656-formula83"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/74656x42.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x43.png" xlink:type="simple"/></inline-formula> denotes the conjugate transpose. The corresponding log-Likelihood function for radar data can then be described as</p><disp-formula id="scirp.74656-formula84"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/74656x44.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x45.png" xlink:type="simple"/></inline-formula> is a constant which does not depend on<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x46.png" xlink:type="simple"/></inline-formula>. To estimate the value of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x47.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x48.png" xlink:type="simple"/></inline-formula>, the cost-function to be maximized is given by</p><disp-formula id="scirp.74656-formula85"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/74656x49.png"  xlink:type="simple"/></disp-formula><p>The first term <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x50.png" xlink:type="simple"/></inline-formula> in (8) does not depend on<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x51.png" xlink:type="simple"/></inline-formula>. The final term <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x52.png" xlink:type="simple"/></inline-formula> is constant with respect to<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x53.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x54.png" xlink:type="simple"/></inline-formula>. By exploiting the relation <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x55.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x56.png" xlink:type="simple"/></inline-formula>, where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x57.png" xlink:type="simple"/></inline-formula> denotes Kronecker product, we have</p><disp-formula id="scirp.74656-formula86"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/74656x58.png"  xlink:type="simple"/></disp-formula><p>It can be easily shown that the product in (9) yields a diagonal matrix with diagonal elements maintain a constant which equals to<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x59.png" xlink:type="simple"/></inline-formula>. Therefore the cost-function in (8) can be simplified as</p><disp-formula id="scirp.74656-formula87"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/74656x60.png"  xlink:type="simple"/></disp-formula><p>In order to find the maximum likelihood estimations of the unknown parameters, it is necessary to perform a three-dimensional (3-D) search in the parameter space. From a practical point of view, such a 3-D search, for all possible values of DOA, velocity and acceleration of a target, can cause tremendous computational burden. Therefore, we propose a two-stage sequential estimation process.</p><p>As we simplify the target moving scenario by assuming the target is moving along a certain direction, a classical sub-space method, multiple signal classification, is used to estimate the DOA of the target which can be treated as static or relatively static. Then the estimation of initial velocity and acceleration can be obtained when the cost function in (10) achieves its maximum.</p></sec><sec id="s4"><title>4. Simulations</title><p>We consider a MIMO radar with 6 antennas both on transmitting and receiving ends. Antennas are uniformly spaced with an distance of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x61.png" xlink:type="simple"/></inline-formula>, where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x62.png" xlink:type="simple"/></inline-formula>. Each transmitter transmits different Hadamard codes with the same carrier frequency. The code length is set to <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x63.png" xlink:type="simple"/></inline-formula> and code duration is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x64.png" xlink:type="simple"/></inline-formula>. The azimuth of the target is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x65.png" xlink:type="simple"/></inline-formula>, the initial velocity of the maneuvering target is <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x66.png" xlink:type="simple"/></inline-formula> and acceleration is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x67.png" xlink:type="simple"/></inline-formula>. Initial velocity searching range is from <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x68.png" xlink:type="simple"/></inline-formula> to<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x69.png" xlink:type="simple"/></inline-formula>, acceleration searching range is from <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x70.png" xlink:type="simple"/></inline-formula> to<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x71.png" xlink:type="simple"/></inline-formula>, the step for searching of initial velocity and acceleration are set to <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x72.png" xlink:type="simple"/></inline-formula>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x68.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/74656x73.png" xlink:type="simple"/></inline-formula>, respectively.</p><p>Monte Carlo simulations with 200 trials per signal-to-noise ratio (SNR) were conducted to acquire the RMSE of the ML estimation. The RMSE of DOA is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>, which indicates that MIMO radar has better angular estimation performance than phased array radar. The RMSE of initial velocity is shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>, and that of acceleration is shown in <xref ref-type="fig" rid="fig3">Figure 3</xref>. Under various SNRs, MIMO radar and phased array radar both can obtain excellent estimation performance in estimating initial velocity and acceleration of maneuvering targets.</p></sec><sec id="s5"><title>5. Conclusions</title><p>This paper studies the performance of parameters estimation of maneuvering target in coherent MIMO radars with the application of maximum likelihood estimation. Numerical simulation results show that MIMO radar has a better angular resolution than conventional phased array radar. The effectiveness of maximum likelihood estimator is also verified through the simulations.</p></sec><sec id="s6"><title>Acknowledgements</title><p>This work is supported by National Nature Science Foundation of China (NSFC)</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> RMSE of DOA</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/74656x74.png"/></fig><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> RMSE of initial velocity</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/74656x75.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> RMSE of acceleration</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/74656x76.png"/></fig><p>under Grants 61471365, U1633107 and 61231017, National University’s Basic Research Foundation of China under Grant No. 2000300446. The work is also supported by the Foundation for Sky Young Scholars of Civil Aviation University of China.</p></sec><sec id="s7"><title>Cite this paper</title><p>Li, H., Xu, L. and Zhang, Z. (2017) Parameter Estimation of Maneuvering Target Using Maximum Like- lihood Estimation for MIMO Radar with Colocated Antennas. Journal of Computer and Communications, 5, 69-74. https://doi.org/10.4236/jcc.2017.53008</p></sec></body><back><ref-list><title>References</title><ref id="scirp.74656-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Goldsmith, A., Jafar, S.A., Jindal, N. and Vishwanath, S. (2003) Capacity Limits of MIMO Channels. IEEE Journal on Selected Areas in Communications, 21, 684-702.  
https://doi.org/10.1109/JSAC.2003.810294</mixed-citation></ref><ref id="scirp.74656-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Cui, S.G., Goldsmith, A.J. and Bahai, A. (2004) Energy-Efficiency of MIMO and Cooperative MIMO Techniques in Sensor Networks. IEEE Journal on Selected Areas in Communications, 22, 1089-1098.  
https://doi.org/10.1109/JSAC.2004.830916</mixed-citation></ref><ref id="scirp.74656-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Li, J. and Stoica, P. (2007) MIMO Radar with Colocated Antennas. IEEE Signal Processing Magazine, 24, 106-114. https://doi.org/10.1109/MSP.2007.904812</mixed-citation></ref><ref id="scirp.74656-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Haimovich, A.M., Blum, R.S. and Cimini, L. (2008) MIMO Radar with Widely Separated Antennas. IEEE Signal Processing Magazine, 25, 116-129.   
https://doi.org/10.1109/MSP.2008.4408448</mixed-citation></ref><ref id="scirp.74656-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Dang, X.F., Ma, K., Chen, B.X. and Yang, M.L. (2016) Direction of Arrival Esti-mation for Monostatic MIMO Radar Using Reduced-dimension RISR Algorithm. 11th International Symposium on Antennas, Prop-agation and EM Theory (ISAPE), Guilin, China, 18-21 October 2016, 607-611.</mixed-citation></ref><ref id="scirp.74656-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Hokiguchi, M., Kikuma, N. and Sakakibara, K. (2016) Target Direction Estimation by MIMO Radar Using Root-MUSIC with Minimum Redundancy Array, International Symposium on Antennas and Propagation (ISAP), Okinawa, Japan, 24-28 October 2016, 1036-1037.</mixed-citation></ref><ref id="scirp.74656-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">He, Q., Blum, R.S., Godrich, H. and Haimovich, A.M. (2010) Target Velocity Estimation and Antenna Placement for MIMO Radar with Widely Separated Antennas. IEEE Journal of Selected Topics in Signal Processing, 4, 79-100.   
https://doi.org/10.1109/JSTSP.2009.2038974</mixed-citation></ref><ref id="scirp.74656-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Jiang, M., Huang, J.G. and Hou, Y.S. (2010) MIMO Radar Joint Estimation of Target Location and Velocity with Multiple Subcarrier Signals. IEEE International Symposium on Phased Array Systems and Tech-nology. Boston, 12-15 October 2010, 398-402. https://doi.org/10.1109/ARRAY.2010.5613338</mixed-citation></ref><ref id="scirp.74656-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Sajjadieh, M.H.S. and Asif, A. (2015) Compressive Sensing Time Reversal MIMO Radar: Joint Direction and Doppler Frequency Estimation. IEEE Signal Processing Letters, 22, 1283-1287. https://doi.org/10.1109/LSP.2015.2396650</mixed-citation></ref><ref id="scirp.74656-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Wang, N., Wang, W.G., Zhang, F. and Yuan, Y.N. (2014) The PARAFAC-MUSIC Algorithm for DOA Estimation with Doppler Frequency in a MIMO Radar System. International Journal of Antennas and Propagation, 2014, 1-5.   
https://doi.org/10.1155/2014/684591</mixed-citation></ref></ref-list></back></article>