gold-music: a variation on music to accurately determine peaks of the spectrum
TRANSCRIPT
![Page 1: gold-MUSIC: A Variation on MUSIC to Accurately Determine Peaks of the Spectrum](https://reader037.vdocument.in/reader037/viewer/2022100512/575092bf1a28abbf6baa0332/html5/thumbnails/1.jpg)
Copyright (c) 2013 IEEE. Personal use is permitted. For any other purposes, permission must be obtained from the IEEE by emailing [email protected].
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication.
100 IEEE TRANS ANTENNA & PROPAGATION
gold-MUSIC: A Variation on MUSIC to AccuratelyDetermine Peaks of the Spectrum
Kaluri V. Rangarao E-mail: [email protected] Venkatanarasimhan E-mail: [email protected] of Electronics and Communication Engineering
Jawaharlal Nehru Technological UniversityHYDERABAD
INDIA
Abstract— This paper presents a new way of accuratelydetermining peaks of the MUSIC [1] spectrum, here consid-ered from the point of view of estimating the directions ofarrival (DOAs) of narrow-band signals. It can be used, withany smart antenna geometry and for any purpose whereMUSIC is applicable.
The MUSIC algorithm for DOA estimation evaluates theMUSIC spectrum for various angles and chooses the max-ima or peaks as the angles of arrival. The values obtaineddepend on the interval at which the spectrum is evaluated.The coarser the interval, the less accurate are the resultsin case of MUSIC. To improve accuracy and not depend onthe interval, Root-MUSIC [2], which involves finding theroots of a polynomial, is available. However, Root-MUSICis applicable, in its original form, only to uniform lineararrays (ULA).
The gold-MUSIC algorithm proposed in this paper is atwo-stage process. The first stage evaluates the objec-tive function at coarse intervals and determines peaks fol-lowed by an iterative approach based on gold-section uni-variate(GSU) minimization [3] to find accurate values ofthe peaks. If the number of peaks found so far is equal tothe number of estimated peaks, the algorithm stops withthis first stage. The second stage is an iterative step for fineresolution using finer intervals around the peaks found so farfor finding peaks that were missing in previous iterations.
This paper, also presents a method, based on a partition-ing algorithm for estimating the number of emitters.
The performance of gold-MUSIC is described, includingits advantages and comparison of time complexities for MU-SIC, Root-MUSIC and gold-MUSIC. The proposed algo-rithm gives good results even when the number of snap-shots is small. This gives it an additional computationaladvantage. It does not compromise on the resolving powerof MUSIC.
Index Terms: gold-MUSIC, Root-MUSIC, Directions ofArrival (DOAs), Gold Section Univariate Minimization,Smart Antenna, Sensor Arrays, Partitioning Algorithm,Resolution Technique.
I. Introduction
Directions of arrival for narrow-band signals have beencomputed using high and super-resolution algorithms suchas those discussed in [1], [2], [4], [5] and [6] over the last fewdecades. In this paper, we present a variation on MUSIC.In this method, which we call gold-MUSIC, the peaks ofthe MUSIC spectrum are evaluated quite accurately usinggold-section univariate(GSU) minimization.
GSU minimization [3], which is an iterative method has
been similarly and successfully used by the authors in de-lay estimation [7] and DOA estimation using spatial tunedfilter [8]. The spatial tuned filter approach [8] is based ona filter described in the book [9].
Stage I of gold-MUSIC finds the peaks in two passes.In the first pass, the spectrum is evaluated as in MUSICfor −90◦ to 90◦ using coarse intervals. This provides theintervals with left and right points within which accuratepeaks are located. The second pass, then estimates veryaccurate measures of the peak DOAs using gold-sectionunivariate minimization. If all the emitters are not foundin Stage I, Stage II finds them by repeating the two passesusing fine sweep angles around the peaks found so far.
II. MUSIC Spectrum
This section revises the concepts of MUSIC spectrum asused in [1]. Let M be the number of sensors in a smartantenna array. Let there be K snapshots at each sensor.For any k ∈ {1, . . . , K}, the kth snapshot at each sen-sor is for the same instant in time. Let Y represent theK × M matrix of snapshots. This matrix has complexvalues in general, representing in-phase and quadraturecomponents. The covariance matrix is given by:
S = YHY (1)
where H denotes Hermitian. The eigenanalysis of S yieldsthe eigenvalues L and the corresponding eigenvectors E.It has been shown that L will contain only positive realvalues as S is positive definite.
If there are D < M independent signals, M − D of theeigenvalues will ideally be 0 under no noise condition, but,will be close to 0 depending on the signal-to-noise ratio(SNR). After sorting the eigenvalues in L in ascendingorder, the M × (M − D) matrix EN , the matrix of theM − D eigenvectors corresponding to the M − D lowesteigenvalues is found. EN contains the noise eigenvectors.E can be written as:
E = {EN |ES} where ES are the signal eigenvectors.
The array manifold for a direction of arrival θ, is theM × 1 column vector a(θ). It depends on the geometry.
![Page 2: gold-MUSIC: A Variation on MUSIC to Accurately Determine Peaks of the Spectrum](https://reader037.vdocument.in/reader037/viewer/2022100512/575092bf1a28abbf6baa0332/html5/thumbnails/2.jpg)
Copyright (c) 2013 IEEE. Personal use is permitted. For any other purposes, permission must be obtained from the IEEE by emailing [email protected].
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication.
KALURI: GOLD-MUSIC: A VARIATION ON MUSIC TO ACCURATELY DETERMINE PEAKS OF THE SPECTRUM 101
For example, for a ULA, it is given by:
at(θ) = [1, e−jω, e−j2ω, . . . , e−j(M−1)ω ] (2)= [1, z−1, z−2, . . . , z−(M−1)] (3)
z = ejω and ω is given by:
ω = 2πd
λsin θ (4)
where d is the spacing between the ULA elements and λis the wavelength corresponding to the center frequencyof the narrow-band signal. The filtered vector VMUSIC
which is analogous to the output of the spatial tuned filteras described in [8] is here given by:
VMUSIC(θ) = at(θ)EN (5)
which is a 1 × (M − D) row vector. Then, the energy inthe filtered output is:
EMUSIC(θ) = VMUSIC(θ)VHMUSIC (θ) (6)
= at(θ)ENENHa∗(θ) (7)
For MUSIC, this energy is minimized, i.e., almost close tozero, when the angle θ corresponds to a DOA. Thus, theobjective function for maximization is:
JN (θ) =1
EMUSIC(θ)(8)
It is interesting to note that the evaluation of JN (θ) is lin-ear in θ while non-linear in ω unlike conventional frequencyresponse. The spectrum can be thought of as the inverseof the magnitude squared of the frequency response of thefilter given by Eq 5.
III. MUSIC, Root-MUSIC and gold-MUSIC
We begin with the signal snapshots Y, which is a K ×M matrix. Covariance matrix S is computed via Eq 1.The next step is eigenanalysis, getting the eigenvalues Land the corresponding eigenvectors E. In this paper, wedescribe a partitioning method in section IV on the setof the normalized eigenvalues, which gives us an estimateof D. We separate the eigenvectors into M − D noiseeigenvectors EN and the remaining D signal eigenvectorsES.
MUSIC, Root-MUSIC and gold-MUSIC differ after thisstep as depicted in Fig 1. In MUSIC, the spectrum JN (θ)(Eq 8) is evaluated for different angles and the peaks arechosen. Numerically, the estimate is limited by the evalu-ation interval. In Root-MUSIC, the complex polynomialsassociated with EN (Eq 5) are solved and the D most com-mon roots are treated as the DOA angles. gold-MUSIC isexplained in section V.
IV. Estimating the Number of Signals
We estimate the number of signals, D, as the numberof significant eigenvalues in L, i.e., the cluster of eigenval-ues which carries maximum mean energy. The size of thecluster provides an estimate of D.
Root-MUSIC
Solve Polynomialin EN. Find DCommon Roots
gold-MUSIC
Coarse JN (θ)SpectrumGSU Minimization
MUSIC
Evaluate JN (θ)SpectrumFind Peaks.
Partition Eigenvectors{EN |ES}
Estimate Emitters D
Eigenvalues LEigenvectors EEigenanalysis of S
Covariance S = YHY
Signal Snapshots Y
�
�
��
�
� ��
Fig. 1. Steps involved in MUSIC, Root-MUSIC and gold-MUSIC
A. Partitioning Algorithm
The partitioning method we propose has a number ofiterations starting from k = 1. After iteration k, the resultis a collection of clusters {p1, p2, ..., pk} which have beenformed and a remainder set Pk, which has to be brokendown into further clusters if possible. Each cluster pk isa set consisting of normalized eigenvalues of covariancematrix S. The algorithm is started with P0 in Eq 9 as theset of normalized eigenvalues for k = 1.
∀x ∈ Pk−1 pk = {x ≥ σk−1} and Pk = {x < σk−1} (9)
where σk−1 =
√∑(x − µ)2
n − 1∀x ∈ Pk−1
in which n = size[Pk−1] and µ = E(Pk−1)
The choice of the dynamic threshold, σk(second moment)is heuristic [10]. The recursive set of operations as given inEq 9 is performed till either set pk or Pk becomes a NULLset, in which k is an integer. In case pk is NULL, the setpk is replaced by Pk.
B. Numerical Example of the Method
A numerical example of how this algorithm works isdemonstrated below. Only for this purpose, narrow-bandsignals from three emitters arriving from 20◦, −35◦ and25◦ were chosen. The SNR was about −2 dB. The numberof snapshots was K = 50 and the number of sensors in thereceiving ULA was M = 20.
{P0} Initial Set of Normalized Eigenvalues0.018377 0.020537 0.024351 0.030716 0.0408480.045091 0.059093 0.060620 0.068971 0.0702510.088147 0.096373 0.104326 0.118996 0.1323630.151811 0.162140 0.233112 0.499889 1.000000Threshold σ0 = 0.226902
![Page 3: gold-MUSIC: A Variation on MUSIC to Accurately Determine Peaks of the Spectrum](https://reader037.vdocument.in/reader037/viewer/2022100512/575092bf1a28abbf6baa0332/html5/thumbnails/3.jpg)
Copyright (c) 2013 IEEE. Personal use is permitted. For any other purposes, permission must be obtained from the IEEE by emailing [email protected].
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication.
102 IEEE TRANS ANTENNA & PROPAGATION
k = 1 size[P0] = 20 size[p1] = 3 and size[P1] = 17{p1} First Partition E[p1] = 0.577667 size = 30.233112 0.499889 1.000000{P1} First Remaining Part0.018377 0.020537 0.024351 0.030716 0.0408480.045091 0.059093 0.060620 0.068971 0.0702510.088147 0.096373 0.104326 0.118996 0.1323630.151811 0.162140Threshold σ1 = 0.045736
k = 2 size[P1] = 17 size[p2] = 11 and size[P2] = 6{p2} Second Partition E[p2] = 0.10118999 size = 110.059093 0.060620 0.068971 0.070251 0.0881470.096373 0.104326 0.118996 0.132363 0.1518110.162140{P2} Second Remaining Part0.018377 0.020537 0.024351 0.030716 0.0408480.045091Threshold σ2= 0.009065
k = 3 size[P2] = 6 size[p3] = 6 and size[P3] = NULL{p3} Third Partition E[p3] = 0.02998 size = 60.018377 0.020537 0.024351 0.030716 0.0408480.045091
Now, having formed the clusters, we choose D as the sizeof the cluster with maximum E[p]. In this case D = 3corresponding to size of p1.
2 4 6 8 10 12 14 16 180
0.1
0.2
0.3
0.4
0.5
0.6
0.7CLUSTER DIAGRAM Snap Shots = 15 at SNR = −2.1042
Estimated Emitters = 3
Str
engt
h
Clusters using Algorithm Eq 9
Fig. 2. Partitioning Algorithm
For a better visual understanding clusters are depictedin Fig 2. Ideally( SNR ≥ +10 dB ) there must be only twoclusters one with full strength (1.0) and second with (0.0).The reason for the pop-up of two other non-zero clustersis due to the presence of noise ( SNR about −2 dB ).
The estimation of number of emitters has been testedexhaustively over different angles in the [−90◦, 90◦] rangeand for various SNR values as found in Fig 6 in sectionVI.
V. DOA Estimation Using gold-MUSIC
The objective function JN (θ) as given by Eq (8) peaksat each DOA. The sharpness of JN (θ) is dependent onthree factors:
1. K the number of snap-shots2. M the number of ULA elements3. SNR Signal to Noise Ratio
The resolving power of the MUSIC algorithm is also re-lated to the above parameters. To combat poor SNR con-ditions parameters K and M are chosen to fit the reality.
A. gold-MUSIC Algorithm
The gold-MUISC algorithm has two stages:
1. Finding a number of emitters using a coarse scan in[−90◦, 90◦]. Let these be the initial D1 ≤ D DOAs.
2. Determining the remaining D−D1 angles using finerscanning.
The second part is not required if all the initial emittersfound are equal to the total number of estimated DOAs.
A.1 Stage I : Algorithm for Initial Angles
By evaluating at suitable points, intervals are foundwithin which the peaks, i.e., the DOAs are present. Then,by using Gold Section Univariate (GSU) minimizationtechnique peaks are located very accurately.
For demonstrating the gold-MUSIC algorithm, wechoose three emitters with relative strengths 1.0, 1.3 and1.2, arriving from 30◦, −10◦ and 26◦ respectively. In thefirst pass, peaks or energy clusters are found by evaluat-ing JN (θ) at suitable intervals in [−90◦, 90◦]. Fig 3(a)(toppart) shows how 2 peaks (of the red curve) have been foundin a MUSIC spectrum evaluated between −90◦ and 90◦ atevery 10◦ interval. The corresponding gold-MUSIC peaksfound are also marked in the same figure by blue stemswith a bubble.
GSU minimization method [3] finds the extreme val-ues of a univariate function accurately. The algorithmworks by successively narrowing the range of values in-side which the extremum is known to exist. To locate thepeaks of JN (θ), we need to consider −JN (θ) as our uni-modal function. We require 4 points which include thelimits of the interval in which the maximum is to be ob-tained. Let [θa, θb] represent the coarse interval such thatθa < θb for a possible peak. This interval is obtained forall possible peaks by performing a coarse sweep in the in-terval [−90◦, 90◦]. We choose intermediate points θc andθd such that θc = θa + g(θb − θa) and θd = θc + g(θb − θc)where g = 3−√
52 , called the Golden Ratio . The objec-
tive function, −JN (θ), is evaluated at θa ,θb, θc and θd. If−JN (θc) ≤ −JN (θd), θb takes the value of θa and θa takesthe value of θd in the next iteration. Otherwise, θa takesthe value of θc and θc takes the value of θd in the next it-eration. Typically a fixed number of iterations is used andthe minimum value is quickly obtained quite accuratelyresulting in DOA value. Using a simple logic this processis repeated till all peaks,i.e., DOA values are found.
In Fig 3(a), for each of the 2 peaks, θa and θb are chosen2 × 10 = 20 degrees apart, keeping the local maximumin between them. In the same figure, θa and θb pointto these positions for one of the peaks. For this figure,noise was present in the signal with SNR about −2 dB.The accurate angles are found corresponding about −10.43(corresponding to −10◦) and 28.67◦ (corresponding to 26◦
![Page 4: gold-MUSIC: A Variation on MUSIC to Accurately Determine Peaks of the Spectrum](https://reader037.vdocument.in/reader037/viewer/2022100512/575092bf1a28abbf6baa0332/html5/thumbnails/4.jpg)
Copyright (c) 2013 IEEE. Personal use is permitted. For any other purposes, permission must be obtained from the IEEE by emailing [email protected].
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication.
KALURI: GOLD-MUSIC: A VARIATION ON MUSIC TO ACCURATELY DETERMINE PEAKS OF THE SPECTRUM 103
−100 −80 −60 −40 −20 0 20 40 60 80 1000
0.5
1
1.5
Nor
mal
ized
JN
Coarse: ULA 32 SNR −2.0593dB Samples/Cycle 5 Step 10 Snap Shots 15
gold−MUSIC peaks θa θ
b
20 22 24 26 28 30 32 34 36 380
0.5
1
1.5
DOA in degrees
Nor
mal
ized
JN
Fine: ULA 32 Snap Shots 400 Sweep 20 to 38 in Step 0.9
gold−MUSIC peaks
20 22 24 26 28 30 32 34 36 380
0.5
1
1.5
Nor
mal
ized
JN
Fine: ULA 32 Snap Shots 100 Sweep 20 to 38 in Step 0.9
gold−MUSIC peaks
−10.4315328.672252
26.23872329.723267
26.98689329.669132
UNABLE toRESOLVEblue color
Fig. 3. gold-MUSIC Spectrum
or 30◦). Even though the actual and estimated numberof emitters is 3, we find only two. This is due to thehigher sampling interval and lower snap-shots, which mayselect only one peak whereas several may be present. Closeresolution is obtained by the technique mentioned in sub-section A.2 below. If the number of peaks found is equalto the number of emitters, there is no need for furtherresolution.
A.2 Stage II: Resolution Technique
If two emitters are too close for the chosen sweep, i.e.,sampling angle, the above procedure may give less numberof peaks as it has been shown for the two angles around30◦. This is a single value of 28.67◦ (corresponding toeither one of 26◦ and 30◦). However, as in this case, therecould be two (or more) peaks such as θ1 and θ2 and onlyone of them, say θ1 will show up by gold-MUSIC usingcoarse scanning. However, here, we know apriori by usingthe partitioning algorithm of section IV that there can bemore (in this case 3) peaks.
To resolve these peaks, we repeat the gold-MUSIC al-gorithm near each of the 2 peaks, but, this time using afine sweep angle around on either side of the respectivepeak and an increased number of snapshots. We obtainthe close resolution as shown in Fig 3(b) for K = 100, giv-ing values about 26.24◦ (corresponding to original angle of26◦) and 29.72◦ (30◦).
For K = 400, where 27◦ was used instead of 26◦ asan original DOA, Fig 3(c) shows how we get a superiorresolution of about 26.99◦ (27◦)and 29.67◦ (30◦) due to thehigher number of snapshots. However, using K = 100 and27◦ instead of 26◦, the fine scan fails to resolve between 27◦
and 30◦. This is shown by the blue line in Fig 3(b). Thisvariation of the results with K shows how increasing thenumber of snapshots (from 100 to 400 in this case) givesmore resolution power with accurate results. However, italso shows that increasing K from 100 to 400 has improvedresolving power only from 4◦ to 3◦. Also, as Fig 4 shows,
increasing the number of snapshots improves resolution.The finer sweep angle used in both of these is about 0.9
degrees. This way gold-MUSIC does not compromise onthe resolution but is limited by the choice of the numberof sensors M , the number of snapshots K and the SNRconditions as in Root-MUSIC.
If the number of angles found so far is still less than D,the algorithm of section A.2 is repeated using still finersweep angles till all the D emitters are found.
B. Snapshots and Coarseness
As the number of snapshots, K is increased, peaks be-come sharper. This effect is depicted in Fig 4. Highercoarse scanning will tend to miss sharp and close peaks.To counter this it is preferable to use less number of snap-shots which impacts the resolving power. This is overcomeby using the fine sweep of stage-II as explained in sub-section A.2 above. Improved resolving ability can be seenin the Fig 4(b), with a higher number of snapshots (in thiscase 500).
−100 −80 −60 −40 −20 0 20 40 60 80 1000
0.2
0.4
0.6
0.8
1
Nor
mal
ized
JN
Sharpness Variations ULA 32 SNR −2.1232dB
20 25 30 35 40 45 500
0.2
0.4
0.6
0.8
1 ZOOM from 21 to 47 degrees
DOA in degrees
Nor
mal
ized
JN
Snapshots 15100500
Fig. 4. Sharpness Variations with Snapshots
A unique dual computational advantage is thus achievedby the gold-MUSIC method by lowering K and the numberof angles to evaluate the spectrum at the same time.
VI. Performance
gold-MUSIC has some advantages which are mentionedin this section. As mentioned in the sub-sections below,gold-MUSIC has quick convergence, so that objective func-tion needs to be evaluated overall less number of times.There is trade-off between sharpness of peaks and the num-ber of snapshots which affects the way the scanning inter-val of the first pass and the number of snapshots are tobe chosen. The time-complexity of gold-MUSIC is derivedand compared with MUSIC and Root-MUSIC algorithms.The accuracy of gold-MUSIC is shown for various SNRconditions in Fig 6.
![Page 5: gold-MUSIC: A Variation on MUSIC to Accurately Determine Peaks of the Spectrum](https://reader037.vdocument.in/reader037/viewer/2022100512/575092bf1a28abbf6baa0332/html5/thumbnails/5.jpg)
Copyright (c) 2013 IEEE. Personal use is permitted. For any other purposes, permission must be obtained from the IEEE by emailing [email protected].
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication.
104 IEEE TRANS ANTENNA & PROPAGATION
0 5 10 15 20 25 30 35−15
−10
−5
0ULA 32 SNR −2.1042dB Samples/Cycle 5
1st D
OA
val
ue
0 5 10 15 20 25 30 3529
29.5
30
30.5
2nd D
OA
val
ue
0 5 10 15 20 25 30 3526.5
27
27.5
28
3rd D
OA
val
ue
Number of iterations
Fig. 5. Convergence of GSU Minimization
A. Convergence of gold-MUSIC
As can be seen from Fig 5, it does not take many iter-ations for gold-section to converge. The first part of thefigure is for −10 degrees and the remaining two below itcorrespond to finer angles around 26 and 30 degrees. Thefixed number of iterations used was 32. However, Fig 5shows the quick convergence for all the three angles, sothat about 10 iterations would have sufficed.
B. Comparison of Complexity
Till the isolation of the noise eigenvectors EN, Root-MUSIC and gold-MUSIC follow the same steps. Thetime-complexity is that of computing covariance and doingeigendecomposition. The detailed analysis of these stepsgive
∆ = O(KM2 + M3)
.After that, conventional Root-MUSIC solves M − D
polynomials of degree M −1 and finds D common roots ineach of the polynomials. The time complexity for polyno-mial rooting is O(Mn) where n >> 2 is a relatively largeinteger for known algorithms. So, the overall complexityof these steps is ∆ + O((M − D)Mn).
In gold-MUSIC, the objective function is evaluated onceper coarse spectrum point and a fixed number of timesper peak. The time complexity for objective functionevaluation is that of multiplying a 1 × M matrix withan M × (M − D) followed by multiplying 1 × (M − D)matrix with its Hermitian. That is O(M2) ignoring D.Considering evaluating the objective funtion as a unit ofcomputation (of complexity O(M2)), 180
CAunits are spent
when the spectrum is evaluated with coarse angle of CA
degrees. If all the D1 initial emitters are scanned finerwith a finer coarse angle FA degrees, the total numberof units added for this is D1 × CA
FA+ D2NF , where D2 is
the number of new peaks for which the objective func-tion is evaluated NF times in the gold-MUSIC for de-
termining the peaks. Hence the overall complexity is:O(KM2 + M3 + M2 × [ 180CA
+ D1NF + D2CA
FA+ D2NF ]).
Assuming NF , CA and FA are fixed values and lettingD = D1 + D2 as the total number of emitters, the overallasymptotic time-complexity is O(KM2 +M3 +DM2).LetCA degrees be the coarse scanning angle of the first passand let NF be the number of iterations used. We havedemonstrated gold-MUSIC with less number of snapshots.K = 15 was used in the 3-emitter example. This improvesthe computational speed of gold-MUSIC.
The following table summarizes the time-complexitiesof MUSIC, Root-MUSIC and gold-MUSIC:
Algorithm Complexity
MUSIC ∆ + O(DM2 × 180CA
)
Root-MUSIC ∆ + O((M − D)Mn)
gold-MUSIC ∆ + O(DM2)
If the number of emitters that show up initially are small,gold-MUSIC also has an advantage in less number of eval-uations of the objective function during fine resolution cor-responding to the less number of initial angles.
C. Accuracy
For evaluating the algorithm, consider θ as true DOAand the corresponding estimate by gold-MUSIC as θg. Letus define error Eg for a given ULA size, number of snap-shots K and SNR conditions as:
Eg(θ) = θg − θ ∀θ ∈ {−90◦ < θ < +90◦} (10)
The standard deviation of the error Eg(θ) is:
σg =
√∑θ Eg(θ)2
Nθ − 1(11)
where Nθ is the number of angles for which θg is estimated.
−80 −60 −40 −20 0 20 40 60 80−100
−50
0
50
100ULA 10 SNR −1.0056dB at 5 Samples/Cycle Step 15 Snap Shots 15
DOA in Degrees
Est
imat
ed D
OA
True DOA Estimated DOA σ
g 0.96357
−3 −2 −1 0 1 2 3 4 50
0.5
1
1.5
2
SNR in dB
σg
mean of σg (100) = 0.3850 .. mean of σ
g (15) = 0.9871
Snap Shots 100 Snap Shots 15
Fig. 6. Performance of Gold-MUSIC
![Page 6: gold-MUSIC: A Variation on MUSIC to Accurately Determine Peaks of the Spectrum](https://reader037.vdocument.in/reader037/viewer/2022100512/575092bf1a28abbf6baa0332/html5/thumbnails/6.jpg)
Copyright (c) 2013 IEEE. Personal use is permitted. For any other purposes, permission must be obtained from the IEEE by emailing [email protected].
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication.
KALURI: GOLD-MUSIC: A VARIATION ON MUSIC TO ACCURATELY DETERMINE PEAKS OF THE SPECTRUM 105
The Fig 6a shows θg against the actual DOA angles fora single emitter case with DOA values ranging from −90◦
to 90◦ for a fixed value of SNR.Fig 6(b) depicts σg for K = 15 and K = 100 for various
SNR values in dB. Each value of σg for a given SNR isgenerated using Eq 10 and Eq 11 in which Nθ was chosenas 180. SNR values are varied from −3 dB to 5 dB insteps of 0.1 dB. As can be seen from the Fig 6(b), themean values are µ(σg)|K=100 = 0.3850o and µ(σg)|K=15 =0.9871o
The bias in DOA estimation is present and dependsupon the values of the number of snapshots K, numberof sensors M , the number of sources D and the SNR asdiscussed in [11].
Discussion of the asymptotic accuracy of Root-MUSICcan be found in [12]. Asymptotic efficiency of the MUSICalgorithm is given in [13].
[14] gives a detailed derivation of the Cramer-Raobounds for MUSIC and explains its relationship to theMaximum Likelihood method. The performance analysisis based on the minimization of the objective function. Itis clear from the paper that the asymptotic bound reachedby different forms of the MUSIC algorithm are the sameand depend on the objective function. Thus the Cramer-Rao bound is reached by gold-MUSIC also, under the con-ditions mentioned in that paper.
VII. Conclusions
We have presented a new variation on MUSIC calledgold-MUSIC. This method is useful for getting accurateresults for different kinds of array geometries where MU-SIC can be applied. gold-MUSIC makes use of gold-sectionunivariate minimization, which is an iterative method. Afixed number of iterations gives good results. Results pre-sented in this paper show that it is a very effective methodto get quite accurate results under different SNR condi-tions. We get resolution measured by the standard devi-ation σg better than 1◦ on the average and always within2◦ where SNR is varied from −3 dB to 5 dB for K = 15and better than 0.5◦ for K = 100. We have analyzedthe time complexity of gold-MUSIC and compared it withRoot-MUSIC.
We have presented a novel method for estimating thenumber of emitters, based on a partitioning approachwhich forms clusters of normalized eigenvalues. Thismethod works well under a wide range of SNR conditions.
The fact that gold-MUSIC works well with less numberof snapshots gives it a computational advantage.
VIII. Acknowledgements
The authors wish to thank Defence Electronics ResearchLaboratories, (DLRL), Defence Research and Develop-ment Organization, (DRDO), Ministry of Defence, Gov-ernment of India, Hyderabad, India for sponsoring the ini-tial parts of this research work. Special thanks are toDr. M. Laxminarayana and Dr. Bhoopathi, of DLRL.We thank Dr. L. Pratap Reddy and Dr. D. Srinivasa
Rao of the Electronics and Communication EngineeringDepartment, Jawaharlal Nehru Technological University,Hyderabad (JNTU H) for helping us by providing supportat the Department.
References
[1] Ralph O. Schmidt; Multiple Emitter Location and Signal Pa-rameter Estimation, IEEE Transactions on Antennas and Prop-agation, Vol. AP-34, No.3, March 1986.
[2] Anthony J. Weiss and Benjamin Friedlander; Direction Findingfor Diversely Polarized Signals using Polynomial Rooting IEEETrans. SP, 1993, 41,(5).
[3] Richard Brent, Algorithms for Minimization without Deriva-tives. Prentice - Hall, Inc. (1973).
[4] J. Capon; High Resolution Frequency-Wavenumber SpectrumAnalysis, Proc. IEEE, Vol. 57, No. 8, August 1969.
[5] R. Roy, A. Paulraj and T. Kailath; Estimation of Signal Param-eters via Rotational Invariance Techniques-ESPRIT, c©1986IEEE.
[6] Richard Roy and Thomas Kailath; ESPRIT-Estimation of Sig-nal Parameters via Rotational Invariance Techniques, IEEETransactions on Acoustics, Speech and Signal Processing, Vol.37, No. 7, July 1989.
[7] Kaluri V. Rangarao and Shridhar V.; An Accurate Numer-ical Method for Estimating the Delay Between Two Omni-Directional Receiving Elements IEEE 30th Digital Avionics Sys-tems Conference 2011, Seattle, Washington, U.S.A., October2011, pp 8C3-4-8C3-7.
[8] Kaluri V. Rangarao, Shridhar V., V. Prathyusha, T. Neela Pa-vani; Estimating the Directions of Arrival of Narrow Band Sig-nals using Spatial Tuned Filter IEEE International Conferenceon Radar 2011, Chengdu, China, 24-27 October 2011, pp 619-622.
[9] Kaluri V. Rangarao and Ranjan K. Mallik; Digital Signal Pro-cessing: A Practitioner’s Approach Wiley.
[10] Harry L. Van Trees; Detection, Estimation, and ModulationTheory, Part 1 John Wiley & Sons.
[11] Fu Li and Yang Lu; Unified Bias Analysis for DOA Estima-tion Algorithms IEEE International Conference on Acoustics,Speech, and Signal Processing, 1993, Volume: 4, pp: 376 - 379
[12] Rao, B.D.; Hari, K.V.S.; Performance Analysis of Root-MusicIEEE Transactions on Acoustics, Speech and Signal ProcessingVolume: 37, Issue: 12 1989, pp: 1939 - 1949
[13] Porat, B. and B. Friedlander, Analysis of the Asymptotic Rel-ative Efficiency of MUSIC Algorithm, IEEE Transactions onAcoustics, Speech and Signal Processing, Vol. 36, No. 4, 532-544, 1988.
[14] Petre Stoica and Arye Nehorai; MUSIC, Maximum Likely-hood,and Cramer-Rao Bound, IEEE Transactions on Acoustics,Speech and Signal Processing. Vol 37. No. 5. May 1989.
Kaluri V. Ranga Rao:
Earned his BSEE from Andhra University in 1974 and MSCS from
the Indian Institute of Science, Bangalore in 1977. He received an
MSSE in 1991 from the US Naval Post Graduate School and a Ph.D.
from the Indian Institute of Technology (Madras) in 1994. He is a
senior member of the IEEE. He is the author of a text on DSP appli-
cations, published by John Wiley(UK). Currently he is a Professor
Emeritus at the Dept of ECE, JNTU-H, Hyderabad, and Director
at Voxta Communications, India.
![Page 7: gold-MUSIC: A Variation on MUSIC to Accurately Determine Peaks of the Spectrum](https://reader037.vdocument.in/reader037/viewer/2022100512/575092bf1a28abbf6baa0332/html5/thumbnails/7.jpg)
Copyright (c) 2013 IEEE. Personal use is permitted. For any other purposes, permission must be obtained from the IEEE by emailing [email protected].
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication.
106 IEEE TRANS ANTENNA & PROPAGATION
Before that, he was the Chief Scientist of Tanla Solutions, VP
Engineering at Sarayu Softech, Chief Engineer at General Elec-
tric (Hyderabad), Vice President Satyam Computer Services, CEO
of GMR Vasavi Infotech and Project Director,Senior Scientist in
DRDO. Principal Investigator for DRDO project. He has been a
Visiting Professor at IITD, Visiting Scholar at Oklahoma State Uni-
versity, Visiting Professor at IIIT-H, and Visiting Scholar at CMU.
Prof Rangarao published 25 IEEE transaction and conference pa-
pers, and holds several US and Indian patents.
Shridhar Venkatanarasimhan: Shridhar received BSCS from
the Indian Institute of Technology (Mumbai) in 1991 and MSCS
from Rutgers University in 1993. Currently he is a Research Scholar
at JNTU-H.