 reconstruction capability, and boosts the computational complexity. To alleviate the mismatch issues, several approaches are conceived, e.g.: in [@Duarte2013], the concept of structured sparsity is utilized to inhibit closely-spaced frequency estimates; the method in [@Malioutov2005] starts with a coarse grid and heuristically iterates between estimating the weights and placing a finer grid around the location of the non-zero weight estimates; in [@Ekanadham2011; @Yang2013; @Hu2013; @Fyhn2015], a less fine grid is used as a baseline and the dictionary matrix is modified to include auxiliary interpolation functions.

In the quest for gridless methods which work directly with continuously parameterized dictionaries, i.e., dictionaries whose parameter ranges in $[-\pi,\pi)$, several works depart from using a static dictionary given by a fixed grid. By including the parameters that dictate the dictionary in the estimation problem, they obtain dynamic dictionary algorithms in which the candidate frequencies and hence the dictionary columns are gradually refined. In [@Austin2013], two such algorithms are designed based on the $\ell_p$ regularized least squares objective by adding a penalty term to prohibit closely spaced frequencies and respectively imposing a hard constraint on the minimum distance between frequencies. The algorithms approximately solve the involved nonlinear estimation and still require an initial grid [@Austin2013]. A different line of works adopts the Bayesian framework and augments the probabilistic model of sparse Bayesian learning (SBL) [@Tipping2001; @Wipf2004] to incorporate the candidate frequencies. In SBL, a sparse weight vector is promoted by selecting a parameterized/hierarchical prior model for its entries [@Tipping2001; @Wipf2004]. Estimation in the augmented model is performed using variational inference methods [@ShutinFleury2011; @Hu2012; @Shutin2013] or maximization of the marginalized posterior pdf [@Hansen2014]. Common to all existing SBL-based approaches is that they restrict to compute point estimates of the frequencies (i.e., MAP/ML estimates), which implies nontrivial maximization of highly multimodal functions (similar to classical ML frequency estimation) in each iteration. The maximization is accomplished approximately by using a grid followed by refinement with Newton’s method or interpolation. Another limitation is that, while providing good reconstruction performance, the SBL-based methods reportedly overestimate the model order, i.e., they consistently output additional spurious
reconstruction capability, and boosts the computational complexity. To alleviate the mismatch issues, several access are conceive, e.g.: in   [ @Duarte2013 ], the concept of structured sparsity is utilize to inhibit closely - space frequency estimates; the method in   [ @Malioutov2005 ] starts with a coarse power system and heuristically repeat between estimating the weight and placing a finer grid around the localization of the non - zero system of weights estimates; in   [ @Ekanadham2011; @Yang2013; @Hu2013; @Fyhn2015 ], a less fine grid is used as a baseline and the dictionary matrix is modify to include auxiliary interpolation function. 

 In the quest for gridless methods which work directly with endlessly parameterized dictionaries, i.e., dictionary whose parameter ranges in $ [ -\pi,\pi)$, several works deviate from using a static dictionary given by a fixed grid. By include the parameters that dictate the dictionary in the estimation problem, they obtain dynamic dictionary algorithms in which the candidate frequencies and hence the dictionary column are gradually refined. In   [ @Austin2013 ], two such algorithms are design based on the $ \ell_p$ regularized least squares objective by adding a punishment term to forbid closely spaced frequencies and respectively levy a hard constraint on the minimum distance between frequencies. The algorithms approximately solve the involved nonlinear estimation and still require an initial grid   [ @Austin2013 ]. A unlike line of work espouse the Bayesian framework and augments the probabilistic model of sparse Bayesian eruditeness (SBL)   [ @Tipping2001; @Wipf2004 ] to incorporate the candidate frequencies. In SBL, a sparse weight vector is promoted by selecting a parameterized / hierarchical prior model for its entries   [ @Tipping2001; @Wipf2004 ]. estimate in the augmented mannequin is performed using variational inference method   [ @ShutinFleury2011; @Hu2012; @Shutin2013 ] or maximization of the marginalized later pdf   [ @Hansen2014 ]. Common to all existing SBL - based approaches is that they qualify to compute point estimates of the frequency (i.e., MAP / ML calculate), which implies nontrivial maximization of highly multimodal functions (similar to classical ML frequency estimation) in each iteration. The maximization is accomplished approximately by using a grid followed by polish with Newton ’s method or interpolation. Another limit is that, while providing good reconstruction operation, the SBL - based methods reportedly overestimate the model order, i.e., they consistently output additional spurious
 refonstruction capability, xnd boosts the eimputavional domplexigy. To alleviate the mismatch iwsues, several approaches ard conceivvd, e.g.: in [@Dyartt2013], the concept of structurcb spaddity ms utilized to lnhibit clovely-spaced frexudney estimates; the method in [@Malioutov2005] ftarts eihh a coarse gryd amq hehgiwtically iterates between estjmating the weights snd placing a finer grid agounf the location of hhe non-zero weidyt estimates; in [@Ekanadham2011; @Yang2013; @Hu2013; @Fihn2015], a less fine grid is used as x baszline and tye difjionary matrmx is iodified to lmclude auxilisry interpolatlon fnnctuons.

In the quest for jridless methods whish work dhrzctly with continuousoy paraketesizea diztikneriss, i.e., ficvionaries wgose paramerer ranges in $[-\pi,\pi)$, xedvtal works delart fwoi using a static dictionary given by a xixsd grid. By including thw parameters that dicjate the dyctionary in the estimation problem, they obtain dfnamir aicuiinary qlhorithms in which the candidate frequencies ags nekce the dictionavy columns are gracuwlki refined. In [@Aurtin2013], tck auch algorithms arf desigged bqsed on tre $\ekl_p$ regularized least squarws objective vy adding a penaltv term to pruhibot clpsely spaced frequencier ans respectivfly imposjvg a hard constrxinn on the minimum distance betwqen frequxnciex. The augorothms wpproximatfly solve the involved nlnlingar esdimation ajd still require an initial griv [@Eustin2013]. A diffgrett nine of corks sdopts the Baresian framewotk and auyments the probagilistir model of s[arse Bayesiat learning (SBN) [@Tipping2001; @Wipd2004] to incorpufate the candicate freqleucies. In WBL, a sparse weighb vecjod is promoted bv welecting a parsmegeryzvd/hmerarsvical prior kodeu fut its entries [@Tipiine2001; @Wilf2004]. Estimation in the augjented model is petfjrmed usung variwtional inferrnce methods [@ShutinVleurb2011; @Hu2012; @Siutin2013] pr iaximization of the marginalizsd posterlor pdf [@Hansen2014]. Cjmmok to all existnng SBL-based approaches is that they resvrict to compute point wstimates of the frgqutncies (i.e., MAP/KL eseimates), wvich implies nontrivual maximization pf highly multimodal runctimns (slmilar to classical ML frequency estimation) in each iteration. The maximizqtion ms accomplishsd alproxpmacelv by usyng e jrid followed by gefinement with Newton’s method oc interpoldtnon. Another limitation is thay, dhile providivg good reconstruction psrformanve, the SBL-based methods reportecly overestimate the modcl order, u.e., thsy consistxntly output additiomal spbrious
reconstruction capability, and boosts the computational complexity. the issues, several are conceived, e.g.: structured is utilized to closely-spaced frequency estimates; method in [@Malioutov2005] starts with a grid and heuristically iterates between estimating the weights and placing a finer grid the location of the non-zero weight estimates; in [@Ekanadham2011; @Yang2013; @Hu2013; @Fyhn2015], a fine is as baseline and the dictionary matrix is modified to include auxiliary interpolation functions. In the quest for methods which work directly with continuously parameterized dictionaries, dictionaries whose parameter ranges $[-\pi,\pi)$, several works depart from a dictionary given a grid. including the parameters dictate the dictionary in the estimation problem, they obtain dynamic dictionary algorithms in which the candidate frequencies hence the are gradually In two algorithms are designed the $\ell_p$ regularized least squares objective penalty term to prohibit closely spaced frequencies and imposing a constraint on the minimum distance between The algorithms approximately solve the involved nonlinear estimation still require an initial grid [@Austin2013]. A different line of works adopts the Bayesian framework the probabilistic model of Bayesian learning (SBL) @Wipf2004] incorporate candidate In SBL, sparse weight vector is promoted by selecting a parameterized/hierarchical prior model its entries [@Tipping2001; @Wipf2004]. Estimation in the augmented model is variational methods [@ShutinFleury2011; @Hu2012; or maximization of the posterior [@Hansen2014]. Common to all approaches that compute estimates the frequencies (i.e., MAP/ML which implies nontrivial maximization of multimodal functions (similar to each iteration. The maximization is accomplished approximately by a grid followed by refinement with Newton’s or interpolation. Another limitation is that, while providing good reconstruction performance, the methods reportedly model order, i.e., they consistently output additional spurious
 reconstruction capability, aNd boosts thE compUtaTioNaL comPlexIty. To alleviate THe miSmatch issues, several appRoachEs ARe coNCeIved, e.G.: in [@DuarTE2013], tHE ConCePt Of sTrUCtUred sParSity is uTilized to iNhiBiT closely-spacED fRequency esTimAtes; the methoD in [@maliouToV2005] stARts wiTh a CoarsE grid aND heuriStically iTeRAtes beTWeen estIMAtIng tHe weights and placiNG a FIner grid around The locAtIOn OF The Non-Zero weight EsTimatES; in [@EkanADhAM2011; @yAng2013; @hU2013; @Fyhn2015], a less finE grid is used AS a bAselinE aNd tHE dictiOnary MaTRix Is modified tO incLude auxilIary inTErpolatIOn functIons.

In The QueSt foR GrIdLesS mEThoDS wHicH WorK directlY wItH contInuoUSLY ParaMetErizEd dicTionaries, i.e., diCtiOnarIEs wHose pArameTer rAnGes in $[-\Pi,\pi)$, seVeral WoRks depart from usIng a Static dicTioNaRy gIvEn by a FIxed grId. BY inCluding The paraMEteRs THAT dIctate the dictionarY iN THe EstimatiOn probLEm, ThEY obtain dYnAmiC dicTIOnary AlgoRItHms in whiCh the cANdIdAte freqUeNcies aNd HenCe tHe dicTIonaRy coluMns are grAdualLY refined. In [@AustIN2013], two such algorIThMS ArE DesiGneD based on the $\Ell_p$ REgulArizED lEasT SquarEs objEcTIvE By adding a penalty terM tO prohiBit clOsely spaced frEquencies aND REspectivEly iMPoSIng a hard constrAint oN the minimuM Distance BetweEn frequeNcies. The aLGOrithms aPprOxiMatEly SOLvE the involved nONLineAr EstimatIon And stilL reQuiRe aN inItIal grid [@AuStin2013]. A difFeReNt LiNe oF workS Adopts thE BAyeSiAn fRamewORk and aUgmenTs thE pRoBAbiListic mODeL OF spaRsE BAyesIan LeArninG (SBL) [@tIppIng2001; @Wipf2004] To incorpoRatE The cAnDiDate freQuencies. In SBL, A sParse weighT vEctOr is prOMOted by seLecting a parameterized/hiERarchicAl pRior mOdel For its entRieS [@TippiNg2001; @WIPf2004]. EstiMation In the AuGmeNTEd modEL Is PerFoRmed using vARIatIonal InFereNce methOds [@ShutinFleury2011; @Hu2012; @SHUtiN2013] or maximizatiOn oF the MARgInaLIzED poStERioR PDf [@Hansen2014]. Common tO all existiNg sbL-Based approACheS iS that thEy restrIct to COmpute pOint estimAtes of the FrEqueNCIes (I.e., MAP/ML estImates), whIch implieS NontrIViAl maxImiZation Of HigHly muLtimodAL fuNctioNs (simiLaR to claSsicaL Ml frequenCy estimation) in each iteraTion. ThE maxiMizAtion is acComPLisHed approxImatEly by using A grId fOllowEd bY RefinEmenT WiTh NEWton’s MethOD or interpOLaTioN. aNoTher limitatION Is tHat, whIle PRovidiNg goOd reconstruction pERformance, the SBl-basED MetHodS RepoRtEdly overestimaTe tHe MODel order, I.e., They consistEntly outPuT AdditIonal sPuriouS
 reconstruction capability , and boos ts th e c omp ut atio nalcomplexity. To alle viate the mismatch iss ues,se v eral ap proac hes are co n c eiv ed ,e.g .: in  [@Du art e2013],  the conce ptof  structureds pa rsity is u til ized to inhi bit  close ly -sp a ced f req uency  estim a tes; t he method  i n  [@Mal i outov20 0 5 ]star ts with a coarseg ri d  and heuristic ally i te r at e s  be twe en estimat in g the weights an d p lac i ng a finer gr id around t h e l ocatio noft he non -zero  w e igh t estimates ; in  [@Ekanad ham201 1 ; @Yang 2 013; @H u2013;  @F yhn 2015 ] ,ales sf ine gr idi s u sed as a  b as eline  and t h e  dic tio nary  matr ix is modifie d t o in c lud e aux iliar y in te rpola tion f uncti on s.

In the ques t fo r gridles s m et hod swhich work d ire ctl y withcontinu o usl yp a r am eterized dictionar ie s ,  i .e., dic tionar i es  w h ose para me ter  ran g e s in$[-\ p i, \pi)$, s everal wo rk s depar tfrom u si nga s tatic dict ionary  given b y a f i xed grid. By i n cluding the p a ra m e te r s th atdictate the  dic t iona ry i n  t hee stima tionpr o bl e m, they obtain dyna mi c dict ionar y algorithmsin which t h e candidat e fr e qu e ncies and henc e the  dictionar y  columns  aregraduall y refined . In [@Aus tin 201 3],  tw o su ch algorithms a re d es igned b ase d on th e $ \el l_p $ r eg ularizedleast sq ua re sob jec tiveb y adding  a  pe na lty  term to pro hibit  clo se ly spa ced fre q ue n c iesan dresp ect iv ely i mpos i nga hardconstrain t o n  the  m in imum di stance betwee nfrequencie s.  Th e algo r i thms app roximately solve the in v olved n onl inear  est imation a ndstillreq u ire an  initi al gr id  [@ A u stin2 0 1 3] . A  d ifferent l i n e o f wor ks  ado pts the  Bayesian framewor k  an d augments th e p roba b i li sti c  m o del  o f  sp a r se Bayesian lea rning (SBL )[ @T ipping2001 ;  @W ip f2004]to inco rpora t e the c andidatefrequenci es . In S BL,  a sparseweight v ector isp romot e dby se lec ting a  p ara meter ized/h i era rchic al pri or  model  forit s entrie s [@Tipping2001; @Wipf2 004].Estim ati on in the  au g men ted model  isperformedusi ngvaria tio n al in fere n ce  me t hods[@Sh u tinFleury 2 01 1;@ H u2 012; @Shuti n 2 0 13]  or m axi m izatio n of  the marginalized posterior pdf[@Ha n s en2 014 ] . Co mm on to all exis tin gS B L-basedap proaches is  that th ey restr ict to  compu te poin t es t imates  ofthe  frequenc ies  ( i .e., MA P/ ML estima tes) ,whichimplie s  non t r ivial maximizati on of h ighly mul timod al  functi o ns ( similar to  classicalML fre quen cy es timatio n)  in ea chit eration. T h e maximiz ation  is acc om plis hed  appro xima t e ly by  usi ng  agrid foll o w ed by  r e fin emen t wit hNewt on’s meth o d or int erp o lation.  A not h e r limi t at i o n is that,  wh ile p r o viding goo d  rec o ns t ructi on per forman ce, the SBL -b ased me tho d s  reported ly overes t ima te  the  model o rd er,i. e.,  t hey consistentl y output  additio n a l spur iou s 
 reconstruction_capability, and_boosts the computational complexity._To alleviate_the_mismatch issues,_several_approaches are conceived,_e.g.: in [@Duarte2013], the_concept of structured sparsity_is utilized to_inhibit_closely-spaced frequency estimates; the method in [@Malioutov2005] starts with a coarse grid and heuristically iterates_between_estimating the_weights_and_placing a finer grid around_the location of the non-zero_weight estimates;_in [@Ekanadham2011; @Yang2013; @Hu2013; @Fyhn2015], a less fine grid_is_used as a_baseline and the dictionary matrix is modified to include_auxiliary interpolation functions.

In the quest for_gridless methods which_work_directly_with continuously parameterized dictionaries,_i.e., dictionaries whose parameter ranges in_$[-\pi,\pi)$, several works depart from using_a static dictionary given by a fixed_grid. By including the parameters that_dictate the dictionary in the_estimation problem,_they obtain dynamic dictionary algorithms_in which the_candidate frequencies_and hence the_dictionary columns are gradually refined. In [@Austin2013],_two such algorithms_are designed based on the $\ell_p$_regularized_least squares objective_by_adding_a penalty_term to prohibit_closely_spaced frequencies_and_respectively imposing a hard constraint on_the_minimum distance between frequencies. The algorithms approximately_solve the involved nonlinear_estimation_and still require an_initial grid [@Austin2013]. A different line_of works adopts the Bayesian framework_and augments_the probabilistic_model of sparse Bayesian learning (SBL) [@Tipping2001; @Wipf2004] to incorporate the candidate_frequencies. In SBL, a sparse weight_vector is promoted by_selecting a_parameterized/hierarchical_prior model for_its_entries [@Tipping2001; @Wipf2004]._Estimation in the augmented model is performed_using variational_inference methods [@ShutinFleury2011; @Hu2012; @Shutin2013] or maximization_of the marginalized posterior_pdf [@Hansen2014]._Common to all existing SBL-based approaches_is that they restrict to compute_point estimates of the frequencies_(i.e.,_MAP/ML_estimates), which implies nontrivial maximization_of highly multimodal functions (similar to_classical ML frequency_estimation) in each iteration. The maximization is_accomplished_approximately by using a grid followed_by_refinement with Newton’s method or interpolation._Another_limitation_is that, while providing good_reconstruction performance, the SBL-based methods reportedly_overestimate the model order, i.e., they consistently output additional_spurious