 to this problem have applications to range and direction estimation in sonar and radar, channel estimation in wireless communications, speech analysis, spectroscopy, molecular dynamics, power electronics, geophysical exploration.

In LSE, the original signal $\vx = (x_0,\ldots,x_{N-1})^{\operatorname{T}}\in\mathbb{C}^N$ is a superposition of $K$ complex sinusoids, i.e., $$\label{eq:OrigSig}
    x_n = \sum_{k=1}^{K}  \alpha_k  e^{j\omega_k n}, %\quad n\in\{0,\ldots,N-1\},$$ where $\alpha_k\in\mathbb{C}$ and $\omega_k\in[-\pi,\pi)$ are the complex amplitude and (angular) frequency, respectively, of the $k$th component. We are given the vector $\vy$ containing $M\leq N$ noisy measurements of those components of $\vx$ with indices in $\mathcal{M}\subseteq\{0,\ldots,N-1\}$, $|\mathcal{M}|=M$. Defining the function $\va:[-\pi,\pi)\rightarrow\mathbb{C}^M$, $\omega\to\va(\omega) = (e^{j\omega m} \mid m\in\mathcal{M})^{\operatorname{T}}$ and the vector $\bm{\epsilon}$ representing additive noise, we write $$\label{eq:SigModel}
    \vy = \sum_{k=1}^{K} \alpha_k \va(\omega_k) + \bm{\epsilon}.$$ The problem of LSE involves estimating the number $K$ of sinusoidal components, also referred to as model order selection, and their associated parameters $\alpha_k$ and $\omega_k$. Even if the model order $K$ is given, LSE is still nontrivial because of the nonlinear dependency of  on the frequencies.[^3]

Prior Work
----------

Under the assumption of known $K$, the $\omega_k$’s are traditionally estimated using the maximum-likelihood (ML) technique or subspace methods, such as [@Schmidt1986; @RoyKailath1989]. The ML method involves the hard task of maximizing a nonconvex function that has a multimodal shape with a sharp global maximum. The maximizer is typically searched using iterative algorithms (e.g., [@Ziskind
to this problem have applications to range and steering estimate in sonar and radar, channel estimation in wireless communications, manner of speaking analysis, spectroscopy, molecular dynamics, might electronics, geophysical exploration. 

 In LSE, the original signal $ \vx = (x_0,\ldots, x_{N-1})^{\operatorname{T}}\in\mathbb{C}^N$ is a principle of superposition of $ K$ complex sinusoids, i.e., $ $ \label{eq: OrigSig } 
     x_n = \sum_{k=1}^{K }   \alpha_k   e^{j\omega_k n }, % \quad n\in\{0,\ldots, N-1\},$$ where $ \alpha_k\in\mathbb{C}$ and $ \omega_k\in[-\pi,\pi)$ are the complex amplitude and (angular) frequency, respectively, of the $ k$th component. We are given the vector $ \vy$ containing $ M\leq N$ noisy measurements of those component of $ \vx$ with indices in $ \mathcal{M}\subseteq\{0,\ldots, N-1\}$, $ |\mathcal{M}|=M$. Defining the affair $ \va:[-\pi,\pi)\rightarrow\mathbb{C}^M$, $ \omega\to\va(\omega) = (e^{j\omega m } \mid m\in\mathcal{M})^{\operatorname{T}}$ and the vector $ \bm{\epsilon}$ representing additive noise, we write $ $ \label{eq: SigModel } 
     \vy = \sum_{k=1}^{K } \alpha_k \va(\omega_k) + \bm{\epsilon}.$$ The trouble of LSE involves estimating the number $ K$ of sinusoidal part, also referred to as exemplar order choice, and their associated parameters $ \alpha_k$ and $ \omega_k$. Even if the model order $ K$ is given, LSE is however nontrivial because of the nonlinear dependency of    on the frequencies.[^3 ] 

 Prior Work 
 ---------- 

 Under the assumption of known $ K$, the $ \omega_k$ ’s are traditionally estimated using the maximum - likelihood (ML) technique or subspace methods, such as   [ @Schmidt1986; @RoyKailath1989 ]. The ML method involves the hard task of maximizing a nonconvex function that has a multimodal shape with a crisp global maximum. The maximizer is typically search use iterative algorithms (e.g.,   [ @Ziskind
 to this problem have appligations to range and dicection estimatkon in sonar and radar, channxl ewtimaupon in wireless commuvications, speech qnalbsis, spectroscopb, molecular dyhwmicv, power electrokics, geophyvical explorathov.

Iu LSE, the original signal $\vx = (x_0,\ldots,v_{N-1})^{\operayogname{T}}\in\mathbb{S}^N$ ix a shievposition of $K$ complex sinusoids, i.e., $$\labtl{eq:OrigSig}
    x_n = \sum_{k=1}^{K}  \alpha_k  e^{j\omega_k n}, %\quaf n\in\{0,\ldots,N-1\},$$ where $\wlpha_k\in\matybb{C}$ qnd $\omega_k\in[-\oi,\pi)$ are the complex ajplitude and (angular) frequency, fespeetively, of jkw $k$hv component. We age given the yvctor $\vf$ contaoning $M\leq N$ npisb meqsurements of those cmmponents of $\vx$ wijh indices iu $\mathcal{M}\subseteq\{0,\ldors,B-1\}$, $|\matvcal{K}|=M$. Dddinkng tie runctiln $\ta:[-\pi,\pi)\rightzrrow\mathbb{X}^M$, $\omega\to\va(\omega) = (t^{j\oivba m} \mid m\in\jathcaj{M})^{\jperatorname{T}}$ and the vector $\bm{\epsilon}$ geprssenting additive noise, we write $$\label{eq:SigMldel}
    \vr = \sum_{k=1}^{K} \alpha_k \va(\omega_k) + \bm{\epsilon}.$$ The problem ox LSE knvilyes dwtlmating the number $K$ of sinusoidal components, zlxo referred to af model ordrr srjection, and tfeir assodiated parameters $\wlpha_k$ wnd $\onega_k$. Eveg if the model order $K$ is given, LSE is stilj nontrivial becausz of the nonuinesr delendency of  on the freqbenciea.[^3]

Prior Work
----------

Knder the xssumption of knuwn $K$, the $\omega_k$’s are traditionwlly estinateb using ghe kaximui-likelihoof (ML) bachnique or subspafe mejhods, vuch as [@Schlidt1986; @RoyKailath1989]. The ML method iitolves the hatd dasn of maxnmizinn a nonconvex ftnction that hcs a mulcimodau shape winh a shar' global maxymum. The maxiklzer is typirally seawchee usung itefxtive algorithks (e.g., [@Ziskpnb
to this problem have applications to range estimation sonar and channel estimation in molecular power electronics, geophysical In LSE, the signal $\vx = (x_0,\ldots,x_{N-1})^{\operatorname{T}}\in\mathbb{C}^N$ is a of $K$ complex sinusoids, i.e., $$\label{eq:OrigSig} x_n = \sum_{k=1}^{K} \alpha_k e^{j\omega_k n}, %\quad where $\alpha_k\in\mathbb{C}$ and $\omega_k\in[-\pi,\pi)$ are the complex amplitude and (angular) frequency, respectively, of $k$th We given vector $\vy$ containing $M\leq N$ noisy measurements of those components of $\vx$ with indices in $\mathcal{M}\subseteq\{0,\ldots,N-1\}$, Defining the function $\va:[-\pi,\pi)\rightarrow\mathbb{C}^M$, $\omega\to\va(\omega) = (e^{j\omega m} m\in\mathcal{M})^{\operatorname{T}}$ and the vector representing additive noise, we write \vy \sum_{k=1}^{K} \alpha_k + The of LSE involves the number $K$ of sinusoidal components, also referred to as model order selection, and their associated parameters and $\omega_k$. the model $K$ given, is still nontrivial the nonlinear dependency of on the ---------- Under the assumption of known $K$, the are traditionally using the maximum-likelihood (ML) technique or methods, such as [@Schmidt1986; @RoyKailath1989]. The ML method the hard task of maximizing a nonconvex function that has a multimodal shape with a maximum. The maximizer is searched using iterative (e.g.,
 to this problem have applicatIons to rangE and dIreCtiOn EstiMatiOn in sonar and raDAr, chAnnel estimation in wirelEss coMmUNicaTIoNs, speEch analYSiS, SPecTrOsCopY, mOLeCular DynAmics, poWer electroNicS, gEophysical exPLoRation.

In LSe, thE original sigNal $\Vx = (x_0,\ldoTs,X_{N-1})^{\oPEratoRnaMe{T}}\in\Mathbb{c}^n$ is a suPerpositiOn OF $K$ compLEx sinusOIDs, I.e., $$\laBel{eq:OrigSig}
    x_n = \sum_{K=1}^{k}  \aLPha_k  e^{j\omega_k n}, %\qUad n\in\{0,\LdOTs,n-1\},$$ WHerE $\alPha_k\in\mathBb{c}$ and $\oMEga_k\in[-\pI,\Pi)$ ARE The COmplex amplituDe and (angulaR) FreQuency, ReSpeCTively, Of the $K$tH ComPonent. We are GiveN the vectoR $\vy$ conTAining $M\LEq N$ noisY measuRemEntS of tHOsE cOmpOnENts OF $\vX$ wiTH inDices in $\mAtHcAl{M}\suBsetEQ\{0,\LDOts,N-1\}$, $|\MatHcal{m}|=M$. DefIning the functIon $\Va:[-\pi,\PI)\riGhtarRow\maThbb{c}^M$, $\Omega\To\va(\omEga) = (e^{j\OmEga m} \mid m\in\mathcAl{M})^{\oPeratornaMe{T}}$ AnD thE vEctor $\BM{\epsilOn}$ rEprEsentinG additiVE noIsE, WE WrIte $$\label{eq:SigModel}
    \Vy = \SUM_{k=1}^{k} \alpha_k \vA(\omega_K) + \Bm{\EpSIlon}.$$ The pRoBleM of Lse InvolVes eSTiMating thE numbeR $k$ oF sInusoidAl ComponEnTs, aLso ReferREd to As modeL order seLectiON, and their assocIAted parameterS $\AlPHA_k$ ANd $\omEga_K$. Even if the mOdel ORder $k$ is gIVeN, LSe Is stiLl nonTrIViAL because of the nonlinEaR depenDency Of  on the frequeNcies.[^3]

Prior wORK
----------

Under thE assUMpTIon of known $K$, the $\Omega_K$’s are tradiTIonally eStimaTed using The maximuM-LIkelihooD (ML) TecHniQue OR SuBspace methods, SUCh as [@scHmidt1986; @RoYKaIlath1989]. ThE ML MetHod InvOlVes the harD task of mAxImIzInG a nOnconVEx functiOn ThaT hAs a MultiMOdal shApe wiTh a sHaRp GLobAl maximUM. THE MaxiMiZeR is tYpiCaLly seArchED usIng iterAtive algoRitHMs (e.g., [@ziSkInd
 to this problem have appl ications t o ran geand  d irec tion  estimation in sona r and radar, channel e stima ti o n in wi reles s commu n ic a t ion s,  s pee ch an alysi s,spectro scopy, mol ecu la r dynamics,p ow er electro nic s, geophysic alexplor at ion . 

InLSE , the  origi n al sig nal $\vx=( x_0,\l d ots,x_{ N - 1} )^{\ operatorname{T}}\ i n\ m athbb{C}^N$ is  a sup er p os i t ion  of  $K$ compl ex  sinu s oids, i . e. , $ $\l a bel{eq:OrigSi g}
    x_n=  \s um_{k= 1} ^{K }   \alp ha_k e ^ {j\ omega_k n},  %\q uad n\in\ {0,\ld o ts,N-1\ } ,$$ whe re $\a lph a_k \in\ m at hb b{C }$ and $\ ome g a_k \in[-\pi ,\ pi )$ ar e th e c o mple x a mpli tudeand (angular)  fr eque n cy,  resp ectiv ely,  o f the  $k$th  comp on ent. We are giv en t he vector  $\ vy $ c on taini n g $M\l eqN$noisy m easurem e nts  o f t ho se components of $ \v x $  w ith indi ces in $\ ma t hcal{M}\ su bse teq\ { 0 ,\ldo ts,N - 1\ }$, $|\m athcal { M} |= M$. Def in ing th efun cti on $\ v a:[- \pi,\p i)\right arrow \ mathbb{C}^M$,$ \omega\to\va( \ om e g a) = (e ^{j \omega m} \ midm \in\ math c al {M} ) ^{\op erato rn a me { T}}$ and the vector  $ \bm{\e psilo n}$ represent ing additi v e noise, w e wr i te $$\label{eq:Si gMode l}
    \vy = \sum_{ k=1}^ {K} \alp ha_k \va( \ o mega_k)+ \ bm{ \ep sil o n }. $$ The proble m of L SE  involv esestimat ing  th e n umb er  $K$ of s inusoida lco mp on ent s, al s o referr ed  to  a s m odelo rder s elect ion,  a nd the ir asso c ia t e d pa ra me ters  $\ al pha_k $ an d  $\ omega_k $. Even i f t h e mo de lorder $ K$ is given,LS E is still  n ont rivial b ecause o f the nonlinear depende n cy of onthe f requ encies.[^ 3]
Prior  Wo r k
---- ------ 

Und er  th e assum p t io n o fknown $K$, t he$\ome ga _k$’ s are t raditionally estim a ted  using the ma xim um-l i k el iho o d( ML)  t e chn i q ue or subspacemethods, s uc h  a s [@Schmid t 198 6;  @RoyKa ilath19 89].T he ML m ethod inv olves the  h ardt a skof maximiz ing a no nconvex f u nctio n  t hat h asa mult im oda l sha pe wit h  asharp  globa lmaximu m. Th emaximize r is typically searched  using  iter ati ve algori thm s  (e .g., [@Zi skin d
 to_this problem_have applications to range_and direction_estimation_in sonar_and_radar, channel estimation_in wireless communications,_speech analysis, spectroscopy, molecular_dynamics, power electronics,_geophysical_exploration.

In LSE, the original signal $\vx = (x_0,\ldots,x_{N-1})^{\operatorname{T}}\in\mathbb{C}^N$ is a superposition of $K$ complex_sinusoids,_i.e., $$\label{eq:OrigSig}
___ x_n = \sum_{k=1}^{K} _\alpha_k  e^{j\omega_k n}, %\quad_n\in\{0,\ldots,N-1\},$$ where_$\alpha_k\in\mathbb{C}$ and $\omega_k\in[-\pi,\pi)$ are the complex amplitude and_(angular)_frequency, respectively, of_the $k$th component. We are given the vector $\vy$_containing $M\leq N$ noisy measurements of_those components of_$\vx$_with_indices in $\mathcal{M}\subseteq\{0,\ldots,N-1\}$, $|\mathcal{M}|=M$._Defining the function $\va:[-\pi,\pi)\rightarrow\mathbb{C}^M$, $\omega\to\va(\omega) =_(e^{j\omega m} \mid m\in\mathcal{M})^{\operatorname{T}}$ and the_vector $\bm{\epsilon}$ representing additive noise, we write_$$\label{eq:SigModel}
    \vy =_\sum_{k=1}^{K} \alpha_k \va(\omega_k) + \bm{\epsilon}.$$_The problem_of LSE involves estimating the_number $K$ of_sinusoidal components,_also referred to_as model order selection, and their_associated parameters $\alpha_k$_and $\omega_k$. Even if the model_order_$K$ is given,_LSE_is_still nontrivial_because of the_nonlinear_dependency of _on_the frequencies.[^3]

Prior Work
----------

Under the assumption of_known_$K$, the $\omega_k$’s are traditionally estimated using_the maximum-likelihood (ML) technique_or_subspace methods, such as [@Schmidt1986;_@RoyKailath1989]. The ML method involves_the hard task of maximizing a_nonconvex function_that has_a multimodal shape with a sharp global maximum. The maximizer is_typically searched using iterative algorithms (e.g., [@Ziskind