Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: parse signal recovery problem is the reconstruction of such signals with characteristic of `` sparsity '' from a set of nonadaptive linear measurements . it has great potential of application on various engineering fields such as coding and information theory , signal processing , machine learning and others ( see papers in website @xcite and reference herein ) . sparse signals contains much less information than their ambient dimension suggests . most of entries in its vector representation are zero ( or negligible ) . so it is possible to reconstruct original signal approximately or even accurately using only a small number of linear measurements . the measurements are of the form @xmath0 where @xmath1 is some @xmath2 measurement matrix , @xmath3 and @xmath4 is original signal . clearly the process of sparse signal recovery could be formulated as solving undetermined linear algebraic equation @xmath5 , where @xmath6 is measurement data . it is well - known that undetermined equation has infinite solutions in general . but if we focus on `` sparse '' solution only , situation will be different . although the operation for finding the most sparse solution of undetermined linear equation is np - hard commonly@xcite , theoretical work in compressed sensing has shown that for certain kinds of measurement matrices , it is possible when the number of measurements m is nearly linear in the sparsity of original signal @xcite@xcite . the two major algorithmic approaches to sparse signal recovery are based on @xmath7-minimization and on greedy methods ( matching pursuit ) . finding the most sparse solution of undetermined linear equation is a @xmath8 optimization problem : @xmath9 it could be solved by @xmath7 relaxation for some measurement matrices @xmath1@xcite . that is , solving ( [ label1 ] ) is equivalent to solving the following @xmath7 optimization problem @xmath10 where @xmath11 for @xmath12 . recently , more stronger sufficient condition called restricted isometry property ( rip ) on measurement matrix @xmath1 to guarantee the equivalence of ( [ label1 ] ) and ( [ label2 ] ) was also proposed @xcite . it was widely accepted that @xmath7-minimization ( [ label2 ] ) was normal path to complete sparse signal recovery . ( [ label2 ] ) is essentially a linear programming problem and technique of convex optimization could be utilized to solve it effectively @xcite . the @xmath7-minimization method provides uniform guarantees for sparse recovery . once the measurement matrix @xmath1 satisfies the restricted isometry condition , this method works correctly for all sparse signals @xmath4 . the method is also stable , so it works for non - sparse signals such as those which are compressible , as well as noisy signals . however , the method is based on linear programming , and there is no strongly polynomial time algorithm in linear programming @xcite . but its efficiency was questionable and most popular software package of convex programming , such as cvx@xcite , is hard to be used in practical application for its low rate of convergence , especially when dimension of target signal is large . on the other hand , greedy algorithms are quite fast , both theoretically and experimentally . it runs by iterating in general . typically , on each iteration , the inner products of residue vector @xmath13 with the columns of measurement matrix @xmath1 is computed and a least squares problem is solved to obtain the estimation of original signal on this iteration . it is hoped that the convergence of iterating could be ensured and the estimator could tend to original signal in fewer steps @xcite . the typical case of greedy algorithms for sparse recovery includes orthogonal matching pursuit ( omp ) @xcite , regularized orthogonal matching pursuit ( romp ) @xcite and compressive sampling matching pursuit ( cosamp ) @xcite . it was shown that omp recovered the sparse signal with high probability and had great speed , but it would fail for some sparse signals and matrices @xcite . the development of romp provides a greedy algorithm with uniform guarantees for sparse recovery , just as that provided by @xmath7-minimization method . furthermore , cosamp improves upon these results and provides rigorous runtime guarantees . however , there are one disadvantage for these two algorithms . firstly , sparsity level must be presented as prior parameter for algorithms . but it is unknown in most practical scenario and must be guessed in advance . once the estimated sparsity level has large difference with actual one , the error of algorithm will increase evidently ( although this error could be analyzed theoretically @xcite ) . the problem becomes more severe in noisy environment . romp and cosamp ca nt adapt their running process to noise condition when noise is presented . actually , noise is inevitable in engineering problem and target signal of our recovery algorithms is always buried in it . there exist few `` pure sparse '' signal in real world . because the dimension of target signal is unknown , it is always given with some margin to avoid possible missing . hence it is hard to extract target signal without including certain amount of noise . this not only has influence on accuracy of algorithm , but also reduce the speed of convergence for algorithms . in fact , some calculation is carried through to estimate noise , however , which is useless at all . a new greedy algorithm for sparse recovery is presented in this paper . compared with romp and cosamp , our new algorithm need no any prior information on sparsity level of target signal . furthermore , it is a kind of `` adaptive '' algorithm which can inspect the existence of noise and adjust the halting condition automatically based on detailed state of noise . besides that , it has uniformly guarantee and good efficiency , just as romp or cosamp . hence our new algorithm is a better choice when signal with unknown sparsity level is to be extracted ( such as compressible signal ) under noisy background . this paper is organized as follows : in section 2 we introduce our new algorithm . section 3 describes some consequences of the restricted isometry property that pervade our analysis . the convergence of theorem is also established for sparse signals in sections 3 . practical consideration for algorithm implementation is provided in section 4 . empirical performance and some numerical experiment is described in section 5 . finally , section 7 presents overall conclusion . the most difficult and important part of signal reconstruction is to identify the locations of the components in the target signal . the common approach adopted by most greedy algorithms is `` local approximating '' , that is , computing the inner products between measurement vector @xmath6 and columns of measurement matrix @xmath1 . we will obtain observation vector @xmath14 and use @xmath15 as `` local approximation '' ( or `` proxy '' ) of target signal @xmath4 . note that @xmath1 is a dictionary and @xmath16 is sparse , so @xmath6 has a sparse representation with respect to the dictionary . it is reasonable that only a few entries of @xmath15 are remarkable , which imply the locations of the components of @xmath4 , and most of its entries are comparatively small . of course , the precondition for argument above is that @xmath1 must satisfy some condition such as rip . intuitively , given sparsity level @xmath17 of @xmath4 , every @xmath17 columns form approximately an orthonormal system . therefore , every @xmath17 coordinates of the observation vector u look like correlations of the measurement vector @xmath6 with the orthonormal basis and therefore are close in the euclidean norm to the corresponding coefficients of @xmath4 . popular greedy algorithms , including omp , romp and cosamp , pay much attention to observation vector @xmath15 and build their estimator of location of components in @xmath4 based on @xmath15 . omp uses one biggest coordinate of @xmath15 . it is argued that using only the biggest could nt provides uniformly guarantee . so romp makes use of the @xmath17 biggest coordinates of @xmath15 , rather than just biggest one , and take a further step of regularization to improve the performance of algorithm . it should be noted that sparsity level @xmath17 is always unknown . cosamp employs more coordinates of @xmath15 , the @xmath18 biggest , to avoid the possible leakage of component in @xmath4 . but @xmath17 must be guessed to be input in romp or cosamp as important parameter . if guessed @xmath17 is smaller than its true value , correct result ca nt be found ; on the other hand , if guessed @xmath17 is set to a very large value ( maybe much larger than true value ) to ensure that all of entries of @xmath4 will enter the view of algorithms , certain amount of noise will presented in our calculation inevitably . how can we make a good guess for sparsity level @xmath17 without any knowledge on its true value ? our new approach , named adaptive orthogonal matching pursuit ( amop ) , chooses appropriate number of biggest entries of observation vector @xmath15 by studying the fine feature of @xmath15 . at each step , @xmath15 is computed by @xmath19 where @xmath13 is the residue vector of last step . unlike other algorithms , amop determines the estimation for @xmath17 by analyzing the trend of entries in @xmath15 arranged by descend order of their amplitude . that is , relative amplitude difference of adjacent elements in above queue is calculated and a threshold is set . once the relative amplitude difference between @xmath20th and @xmath21th element in ordered queue of entries in @xmath15 exceed threshold , @xmath20 will be chosen as estimation of @xmath17 . detailed description of amop is proposed as follows :    [ label8 ]    * * input : * measurement matrix @xmath22 , measurement vector @xmath23 , threshold @xmath24 , @xmath25 and @xmath26 . * let @xmath27 , @xmath28 . * calculate @xmath19 and @xmath29 . * arrange @xmath30 by descend order to obtain @xmath31},|u|_{[2]},\cdots,|u|_{[n]}),\label{label3}\ ] ] * determine index @xmath20 as follows , let @xmath32 , @xmath33}-|u|_{[k]}}{|u|_{[k]}}>t*\beta\right.\right\},\ ] ] * if @xmath34 and @xmath35 , set @xmath36 ; else @xmath37 , goto step ( 4 ) ; * update the set of indices by @xmath38,\cdots,[k]\}$ ] . * solve least square problem @xmath39 where @xmath40 is a submatrix of @xmath1 composed of its columns with index in @xmath41 . * calculate the residue vector of @xmath42 . * if @xmath43 , output @xmath44 and @xmath41 , stop ; else go to step ( 2 ) .    as input , the aomp algorithm requires two adjustable parameter @xmath24 and @xmath25 besides matrix @xmath1 and measurement vector @xmath6 . but it does nt need sparsity level of target signal @xmath4 anymore , unlike romp and cosamp . it can be extracted incidentally along with running of algorithm . furthermore , the number of components selected in the step ( 4 ) is determined by algorithm itself automatically . it is easy to understand that this number is critical for performance of algorithm . any manual setting will introduce extra error when mismatch between prescribed value and actual situation of data exists . so it is very necessary to let greedy algorithm of sparse recovery be adaptive , just as aomp . step ( 5 ) should be noted that it give amop algorithm more flexibility and stability . if threshold @xmath24 is set too large so that too much coordinates was selected in one iteration , algorithm is prone to degrade or crash . for this we build a upper bound in step ( 5 ) to prevent the crazy growing of number of chosen components . if this bound is exceed , threshold @xmath24 will be adjusted to smaller value to increase the possibility of components in @xmath30 in step ( 3 ) to satisfy the condition in step ( 4 ) . the importance of step ( 5 ) is also illustrated in following section on analysis of convergence . there are two kinds of iterative invariant of greedy algorithm for sparse recovery deduced in the convergence analysis for romp and cosamp respectively . as to cosamp , assume the sparsity level @xmath45 is preliminary and @xmath46 where @xmath47 is s - sparse approximation for @xmath4 and @xmath48 is additive noise , the following assertion could be proved @xcite . @xmath49 where @xmath50 is result of pruning step in cosamp . hence it is forced to be s - sparse . so this kind of iterative invariant is not suitable for analysis of amop because the sparsity level of intermediate result at each step in amop is nt fixed . however , the iterative invariant in romp simply concerns with the percentage of the coordinates in the newly selected set belong to the support of target signal @xmath4 . it is argued that the ratio above is nt lower than 50% with the help of regularization step @xcite . we will follow the idea in @xcite to derive our result on convergence of amop .      by induction on the iteration of aomp , we study the gain in each iteration . losing no generality , suppose sparsity level of target signal @xmath4 be @xmath51 , and @xmath20 coordinates is selected eventually in this iteration . then its percentage of energy of first @xmath20 components in queue ( [ label3 ] ) is @xmath52}^2}{\sum_{n=1}^k|y|_{[n]}^2+\sum_{n = k+1}^k|y|_{[n]}^2}\ ] ] for the descend order of queue ( [ label3 ] ) , we have @xmath53}^2}{\sum_{n=1}^k|y|_{[n]}^2+(k - k)|y|_{[k]}^2}\nonumber\\ & \geq&\frac{k|y|_{[k]}^2}{\sum_{n=0}^{k-1}(1-t)^{-2n}|y|_{[k]}^2+(k - k)(1-t)^2|y|_{[k]}^2}\nonumber\\ & = & \frac{k}{\sum_{n=0}^{k-1}(1-t)^{-2n}+(k - k)(1-t)^2}\nonumber\\ & = & \frac{k}{\frac{1-(1-t)^{-2k}}{1-(1-t)^{-2}}+(k - k)(1-t)^2}\label{label4}\end{aligned}\ ] ] it is easily seen that [ label4 ] achieves its minimum at @xmath54 , that is @xmath55 here @xmath54 means only the largest coordinates was chosen , which is just the choice of omp . so omp is a special case of amop . according to @xcite lemma 3.6 , a large portion of energy of unidentified part of target signal would be locked by queue ( [ label3 ] ) and certain amount of energy would be reserved by `` regularization step '' in romp by @xcite lemma 3.8 . in amop , the `` regularization step '' is replaced by choosing @xmath20 largest coordinates , so more energy would be identified in amop than in omp because more than one coordinates would be chosen in amop generally . the ability of locking uncovered energy of target signal for amop and romp is compared in fig.[fig1 ]    it is shown that ability of amop is superior to that of romp when sparsity of target signal is small . although the advantage becomes vague when sparsity grows , amop still has relatively good performance considering more coordinates would be chosen and percentage of actual identified energy would be larger than ( [ label5 ] ) .      in @xcite and @xcite , the correctness of support of solution for oga algorithm under different noise scenario were analyzed . mutual coherence of matrix @xmath1 , overcomplete dictionary system , was the key parameter for performance of greedy algorithm in discussion therein . here we will give analogous results for amop under noisy condition based on ric ( restricted isometry constant)@xcite of @xmath1 . in practice the noise is unavoidable and it is always assumed that ideal noiseless signal @xmath56 has sparse representation @xmath57 , where the support of @xmath58 is very small . what we can observe is noisy version @xmath59 where @xmath60 . suppose @xmath58 be solution of @xmath61 @xmath62 be final result of amop , and @xmath63 we argue that @xmath64 under certain conditions on value distribution of target signal . that is , the correctness of support of solution for amop can be guaranteed even in noisy environment . if target signal @xmath65 satisfy @xmath66 here matrix @xmath1 is supposed to has @xmath26-order restricted isometry constant @xmath67 , @xmath26 is twice of sparsity level of target signal @xmath58 . then @xmath64 holds throughout the iteration process of amop unless all the coordinates in @xmath68 were chosen . we proceed by induction . @xmath64 holds at beginning of step 1 in amop initially for @xmath69 . assume it is true at beginning of step 1 in given iteration , we prove it is still true at beginning of step 1 in next iteration . consider case of @xmath70 , we have @xmath71 it is trivial that @xmath72 because @xmath73 is the solution of least square optimization on step 7 in amop , @xmath74 multiply @xmath75 on two side of [ label7 ] , @xmath76 we have @xmath77 where @xmath78 is the identity matrix .    for @xmath79 , the norm of @xmath80 is estimated and bounded . firstly , because @xmath81 is the projection matrix of orthogonal complement of subspace spanned by @xmath82 , its 2-norm is 1 . so @xmath83 secondly , @xmath84 because @xmath85 and @xmath86 , according to @xcite , proposition 3.2 , we have @xmath87 and @xmath88 on the other hand , according to definition of ric , we obtain @xmath89 hence @xmath90 thirdly , by analogous deduction , if @xmath91 , @xmath92 summarize the results above , we have for @xmath91 , @xmath93    lower bound for @xmath94 is considered similarly . for @xmath95 , @xmath96 hence @xmath97 if target signal @xmath98 satisfy ( [ label11 ] ) , we have @xmath99 that is to say , @xmath100 will hold throughout the iteration process of amop unless all the coordinates in @xmath68 were chosen . it is argued that value distribution of target signal , power of noise and ric of matrix @xmath1 are all critical to performance of greedy algorithm . the proposition above gives a general condition for correctness of support of solutions for a large class of greedy algorithms ( not just amop ) which use inner product between residue @xmath13 and dictionary vectors ( columns of @xmath1 ) to obtain information of support of target signal . it seems that condition ( [ label11 ] ) is too restricted . but it is easy to see from ( [ label11 ] ) that `` dynamic scope '' of target signal ( that is , the norm difference between the elements with maximal and minimal norm ) depends on ric @xmath67 of matrix @xmath1 and noise power @xmath25 . consider the requirement on ric in romp , which is @xmath101 with @xmath45 is sparsity level of target signal according to @xcite , theorem 3.1 , we write ( [ label11 ] ) as @xmath102 with maximum is normalized to 1 . it is depicted in fig.[fig2 ] for noise level is 0.1(snr is 20db ) . when sparsity level is small , target signals with considerable dynamic scope are guaranteed to have good performance in greedy algorithms . the restriction on dynamic scope of target signal becomes tighter gently when sparsity level increases . the actual number of chosen coordinates in iteration step of amop in practical scenario is smaller than sparsity level in general . so amop could choose correct coordinates in most cases . this assertion will be illustrated further in numerical experiments . for efficient implementation of amop , it is important to design a appropriate computational scheme with low complexity and good numerical stability to calculate the solution of least square problem in step 7 in amop . it should be noticed that amop is incremental , that is , in each iteration some new coordinates were selected and none of previously chosen coordinates was excluded , unlike cosamp . so it is possible to construct a recursive algorithm to solve the least square problem . assume on the first iteration several coordinates were chosen and denote the set of corresponding columns of matrix @xmath1 as @xmath103 . then observation vector @xmath6 would be linearly approximating with vectors in @xmath103 and coefficients @xmath104are computed by solving the following equation @xmath105 ordinary solver with good numerical performance such as qr decomposition or singular value decomposition could be used to compute @xmath104 . from geometrical point of view , calculation of @xmath104 is equivalent to project @xmath6 onto subspace spanned by @xmath103 . we wrote it intuitively as @xmath106 on the next iteration , a set @xmath107 of columns of matrix @xmath1 with respect to newly chosen coordinates would be added to least square regression . the projection became @xmath108 it is well - known that orthogonalization could simplify the calculation for projection onto subspace spanned by two mutually orthogonal subspace could be regarded as sum of projections on each one . so @xmath107 was written as @xmath109 where @xmath110 was projection of @xmath107 onto @xmath103 and @xmath111 is orthogonal to @xmath103 . hence @xmath112 this could be accomplished with gram - schmidt orthogonalization procedure . without loss of generality , suppose @xmath113 then @xmath114 using following procedure @xmath115 where @xmath116 denotes the inner product of vectors in euclidean space . we can obtain @xmath117 the projection in ( [ label12 ] ) could be written as @xmath118 the numerical stability of gram - schmidt orthogonalization procedure could be improved further @xcite . instead of computing the vector @xmath119 as ( [ label13 ] ) , it is computed as @xmath120 this approach gives the same result as the original formula in exact arithmetic , but it introduces smaller errors in finite - precision arithmetic . there are one points worth mentioning . although other orthogonalization algorithms such as householder transformations or givens rotations are more stable than the stabilized gram - schmidt process , they produce all the orthogonalized vectors only at the end . on the contrary , the gram - schmidt procedure produces the @xmath121th orthogonalized vector after the @xmath121th iteration and this makes it the only choice for iterative algorithm like amop . amop was designed to be a practical method for sparse signal recovery . the main barrier for algorithm efficiency is least square problem in step 7 of amop . using recursive orthogonalization procedure above could mitigate computational burden of amop dramatically . furthermore , the orthogonalization technique has the additional advantage that they only interact with the matrix @xmath1 through its action on vectors . in fact , it only concern with the inner products and additions of columns of matrix @xmath1 . it follows that the algorithm performs better when this sampling matrix has a fast matrix - vector multiply , such as on parallel computational platforms . on the other hand , less memory consumption is another advantage of recursive orthogonalization based least square . in fact , direct method such as svd and qr have storage cost @xmath122 , where @xmath20 is the number of chosen coordinates in each iteration and @xmath123 is row number of matrix @xmath1 . it is too huge for large scale problems . but for amop , only one vector need be put in memory in recursive orthogonalization calculation and storage cost is @xmath124 . it is more suitable for implemented with vlsi circuit . we estimate the time complexity of main steps in amop as follows :    * * step 2 : * in this step , the inner products of residue vector and columns of matrix @xmath1 is computed as proxy for support of target signal . its cost is bounded by matrix - vector multiply @xmath125 , which is @xmath126 with standard multiply operation or @xmath127 for fast matrix - vector multiply . * * step 3 : * according to standard textbook on algorithms @xcite , the expected time for selecting largest @xmath45 entries in vector with dimension @xmath128 is @xmath129 . using efficient schemes such as quicksort or heapsort , a fully sorting of vector could be completed with expected time cost @xmath130 and largest @xmath45 entries could be selected directly , which is faster n some situation . * * step 4 & 5 : * certain amount of support of target signal would be identified in these two steps . although sometimes the threshold needs to be adjusted according to step 5 and several cycles of operations may be necessary , the total cost is still @xmath131 . * * step 7 : * the main advantage of amop is recursive orthogonalization based implementation of least square problem . inner products of vectors are involved in orthogonalization and occupy much of computational resource which can be implemented efficiently by matrix - vector multiply . the cost is @xmath132 with standard multiply operation or @xmath127 for fast matrix - vector multiply . table 1 summarizes the discussion above in standard multiply operation and fast matrix - vector multiply with cost @xmath133 , respectively . .time complexity of amop [ cols="^,^,^ " , ]     storage cost for amop is also considered for showing its practicability . aside from the storage required by the sampling matrix , amop algorithm constructs only one vector of length n as the signal proxy . the sample vectors have length m , so they require @xmath124 storage . the signal approximations require at most @xmath134 storage . similarly , the index sets that appear require only @xmath134 storage . the total storage is at worst @xmath135 . in this section some numerical experiments were conducted to illustrate the performance of signal recovery of aomp . there are three factors to be considered in numerical testing of amop : construction of @xmath1 matrix , value distribution of target signal and snr condition which will be examined in our experiments . the property of measurement matrix @xmath1 is critical to performance of any greedy algorithm for sparse recovery . as indicated in section on theoretical analysis , its ric has direct influence on probability of recovery of algorithms . here , several kinds of matrix @xmath1 were built and utilized to test the performance of amop , including well - known gaussian random matrix , bernoulli random matrix and random fourier matrix . the target signal was set to be flat and no noise was added in , then 500 independent trials were performed . figure [ fig3]-[fig5 ] depicts the percentage ( from the 500 trials ) of sparse flat signals that were reconstructed exactly when @xmath136 gaussian random matrix was chosen as measurement matrix @xmath1 . this plot was generated with @xmath137 for various levels of sparsity @xmath51 . the horizontal axis represents the number of measurements @xmath123 , and the vertical axis represents the exact recovery percentage .    as comparison , recovery percentage of algorithm cosamp is also given under the same setting . standard cosamp needs sparsity of target signal as its important prior knowledge and it is widely regarded as one of main drawbacks of cosamp . in our experiment , sparsity of target signal was given to cosamp as input parameter to guarantee the power of cosamp to be exploited fully . it should be noted that even the sparsity of target signal is known beforehand ( which is impossible in practice ) , recovery percentage of cosamp is lower than that of amop . especially when sparsity was relatively large , performance of cosamp degenerated very rapidly . on the contrary , the behavior of amop was very stable . according to well - known theoretical result of compressed sensing , for gaussian random measurements matrix @xmath1 , if row number @xmath123 , column number @xmath128 and sparsity @xmath51 satisfies @xmath138 where @xmath139 is a constant independent of @xmath51 , then the probability of recovery failure is exponentially small @xcite . our experiment result indicates that for amop , the value of constant @xmath139 is about @xmath140 for gaussian measurement matrix . figure [ fig6]-[fig8 ] depicts corresponding result for bernoulli random measurement matrix , which is analogous to gaussian case . it had been proved that condition ( [ label20 ] is also sufficient for overwhelming probability of successful recovery for binary bernoulli measurement matrix @xcite . it is observed that the constant @xmath139 for amop in bernoulli case is probably the same as that in gaussian case . figure [ fig9]-[fig11 ] depicts corresponding result for fourier random measurement matrix . somewhat surprisingly , it is similar to that of gaussian and bernoulli case . to our knowledge , the best known bounds on size of measurements in fourier case is given by @xcite @xmath141 which is conjectured to be the same as [ label20 ] @xcite but there exists no strict theoretical proof until now . our experiment result verified this conjecture in some extent indirectly . pure flat signal is rarely seen in practical engineering application . so it is necessary to investigate the performance of sparse recovery algorithms on non - flat target signal . there are two cases to be studied . one is piecewise flat signal which is common in various fields of imaging , such as optical , microwave and magnetic resonance . the result is depicted in figure [ fig12 ] to [ fig14 ] . here the measurement matrix is fixed to gaussian random matrix . recovery percentage of cosamp in our experiment is nt ideal . especially when the number of non - zero elements of target signal is large , the successful recovery probability of cosamp becomes more and more ignorable . the capability of cosamp is doubtful in this setting . on the contrary , the behavior of amop is very robust when value distribution of target signal is changed . furthermore , the constant @xmath139 in [ label20 ] is approximately equal to that in flat signal setting . although it is predicted theoretically that constant @xmath139 only depends on the construction of measurement matrix @xmath1 , not rely on nature of target signal , it is common in practice that the detailed value distribution of target signal certainly has influence on performance of recovery algorithms . our experimental result indicates that the actual behavior of amop coincides with theoretical conclusion perfectly . this argument is confirmed by result for the other case of non - flat target signal . here target signal is chosen as compressible signal which is frequently presented in orthogonal representation of signal , such as discrete cosine transform and wavelet transform , and data compression . two kinds of compressible signal are analyzed in our experiment , one is exponential signal , @xmath142 the other is polynomial signal @xmath143 for brevity , only result for exponential signal is depicted in [ fig15 ] and corresponding curve for cosamp is omitted . we also performed the same experiment for polynomial compressible signals and found the results very similar to those in figure [ fig14 ] . random noise was added in target signal to test the performance of sparse recovery algorithms in noisy environment . measurement matrix is fixed to fourier random matrix and target signal is set to flat . the sparsity of target signal is fixed to 20 and size of measurements @xmath123 is fixed to 200 . the relative error of amop and cosamp in gaussian white noise with various level is depicted in figure [ fig16 ]    the capacity of amop in noisy environment is satisfactory . relative recovery error could be controlled within @xmath144 when snr is about 10db . even when snr is as low as 5db , relative error of amop still could be governed within @xmath145 . though the error curve rise very acutely in very low snr region , it is shown that amop works normally in most noisy environment . on the other hand , the relative error of cosamp keeps on a high level when noise is presented . when snr was increased , it did nt exhibit the obvious trend of decreasing . with the language of statistics , cosamp is nt consistent estimator . figure [ fig17 ] illustrates the advantage of amop in noisy environment from other viewpoint . here snr is fixed to 10db , the error with various size of measurements is plotted . it is observed that if number of measurements @xmath17 is small , the performance of amop is heavily abnormal . in fact , relative error of amop is even higher than cosamp when @xmath17 is lower than 100 . but it falls abruptly when @xmath123 increases while that of cosamp remains in narrow range . when @xmath123 is larger than 100 , the relative error of amop tends to be stable . it is well controlled with in @xmath144 , which is similar to figure [ fig16 ] . we would build a measurement matrix @xmath1 with spatial - temporal basis vectors , which is key component in theory and computation of stap ( space - time adaptive processing ) , to investigate the potential of sparse recovery algorithms to be applied in field of modern radar and communication engineering . here @xmath1 is set to @xmath146,\ ] ] and @xmath147^t,\nonumber\end{aligned}\ ] ] where @xmath148 and @xmath149 denotes doppler and spatial frequency , respectively . unlike fourier matrix , the construction of spatial - temporal matrix @xmath1 is more complex . the phase of elements in @xmath1 composed of two parts , one is contributed by doppler frequency and the other by spatial frequency . it lead to following consequence : different @xmath148 and @xmath149 could be combined to form the same ( or approximately equal ) phase . in other words , strong correlation exists in different column vectors in @xmath1 , which corresponding to various points far away with each other on spatial - temporal plane . so restricted isometric constant of @xmath1 is conjectured to be relatively large . it is a challenge for sparse recovery algorithm to be feasible when measurement matrix @xmath1 is chosen as spatial - temporal matrix .    generally speaking , the support ( `` position '' of non - zero elements ) of target signal is much important than its detailed value . it is indeed true in practical engineering discipline . for example , in radar stap processing , sparse recovery is utilized to estimate the energy distribution of clutter and interference on spatial - temporal plane from sample data directly . because echo of clutter and interference is much stronger than radar target , the detailed amplitude and phase of clutter and interference is nt crucial . as long as the accurate support ( `` position '' ) of clutter and interference on spatial - temporal plane is found out , we can design the efficient filter to suppress clutter and interference effectively . hence the most important feature of sparse recovery algorithm applied to stap processing is detecting the support of target signal with great precision . the experiment was performed to test the ability of amop and cosamp to detect signal support . measurement matrix @xmath1 was set to @xmath150 spatial - temporal matrix as [ label21 ] . target signal is noise - free flat signal and its sparsity is fixed to 20 . figure [ fig18 ] depicts the average result from 100 trials . it is clear that amop is much superior to cosamp when applied to stap processing because it could detect a majority of target support with no error while cosamp could only find very few . but even so , the accuracy of algorithms for support , whether amop or cosamp , can hardly satisfy the requirement of practical stap processors . due to ultra - low signal clutter ratio ( generally lower than -50db ) , missing four or five frequency points on spatial - temporal plane would lead to very high false alarm rate and the performance of radar would degrade heavily . so we should detect as much target support as possible to minimize false alarm rate . if tiny error is allowed , the behavior of sparse recovery algorithms becomes better , as depicted in figure [ fig19 ]    it is easily seen that if the `` position '' detected having difference 1 with true `` position '' of target signal is allowed to be counted , the average size of detected target support for amop increases by about 2 and that of cosamp is still negligible . it accounted for that some support of target signal missed by amop was nt really missed . that is to say , their neighborhood , which corresponding to the points adjacent to them on spatial - temporal plane , were discovered instead . this is a good news for high performance filter to suppress clutter could still be designed with amop and carefully chosen notch , without losing much resolution and snr . figure [ fig20 ] depicted the case of error 2 . the behavior of amop continued to be made better and approach the best . it seems that it make little sense to enlarge error tolerance further . in this paper , a novel greedy algorithm for sparse recovery , called amop , was given and examined . its performance was studied by theoretical analysis of simulation experiment . the motivation of this algorithm is two obvious drawbacks in popular methods , such as cosamp : need of sparsity of target signal as prior knowledge , and weak ability of working in noisy environment . with well - designed algorithmic steps , amop can extract the information on sparsity of target signal adaptively and sense the nature of target signal automatically . it can recover the detail value of target signal with very high precision with little prior knowledge . its validity is illustrated by strict deduction . fine stability of performance under random noise is another advantage of amop . furthermore , its robustness for various setting of target signal , flat or compressible , and construction of measurement matrix , such as spatial - temporal matrix , were also shown by thorough numerical experiment . it is argued that amop is a excellent greedy algorithm for sparse recovery and has great potential of widely utilization on signal processing .    1 compressed sensing webpage : http://www.dsp.ece.rice.edu/cs . m.r . garey and d.s . johnson , computers and intractability , a guide to the theory of np - completeness , w. h. freeman and company , new york , 1979 e. candes . compressive sampling . in proc . international congress of mathematics , volume 3 , pages 1433 - 1452 , madrid , spain , 2006 . e. candes and m. wakin , an introduction to compressive sampling . ieee signal processing magazine , 25(2 ) , pp . 21 - 30 , march 2008 . d. l. donoho and p. b. stark . uncertainty principles and signal recovery . siam j. appl . , 49(3):906 - 931 , june 1989 . candes and t. tao . decoding by linear programming , ieee trans . theory 51 , 4203 - 4215 , 2006 . s. boyd , and l. vandenberghe , convex optimization , cambridge university press , 2004 d. needell , r. vershynin , uniform uncertainty principle and signal recovery via regularized orthogonal matching pursuit , foundations of computational mathematics , doi : 10.1007/s10208 - 008 - 9031 - 3 cvx official website : http://www.stanford.edu/  boyd / cvx/ d. needell , j.a . tropp , r. vershynin , greedy signal recovery review , arxiv , 0812:2202 , dec , 2008 j. a. tropp and a. c. gilbert . signal recovery from random measurements via orthogonal matching pursuit . ieee trans . theory , 53(12 ) : 4655 - 4666 , 2007 . d. needell and j. a. tropp . cosamp : iterative signal recovery from incomplete and inaccurate samples . acm technical report 2008 - 01 , california institute of technology , pasadena , july 2008 . h. rauhut . on the impossibility of uniform sparse reconstruction using greedy methods . to appear , sampl . theory signal image process . , 2008 . d. l. donoho , m. elad , and v. n. temlyakov , `` stable recovery of sparse overcomplete representations in the presence of noise , '' ieee trans . theory , vol . 6 - 18 , jan . 2006 . tseng , `` further results on stable recovery of sparse overcomplete representations '' , ieee transactions on information theory , vol . 2 , february 2009 . bau iii , david ; trefethen , lloyd n. , numerical linear algebra , philadelphia : society for industrial and applied mathematics , 1997 . t. cormen , c. lesierson , l. rivest , and c. stein . introduction to algorithms . mit press , cambridge , ma , 2nd edition , 2001 . e. candes and t. tao , near optimal signal recovery from random projections : universal encoding strategies ? ieee trans . on information theory , 52(12 ) , pp . 5406 - 5425 , december 2006 . rudelson , m. , vershynin , r. , sparse reconstruction by convex relaxation : fourier and gaussian measurements . preprint , 2006 . hao zhang , gang li , huadong meng , a class of novel stap algorithms using sparse recovery technique , arxiv:0904.1313 , apr , 2009
Abstract:

Answer: greedy algorithm are in widespread use for sparse recovery because of its efficiency . but some evident flaws exists in most popular greedy algorithms , such as cosamp , which includes unreasonable demands on prior knowledge of target signal and excessive sensitivity to random noise . a new greedy algorithm called amop is proposed in this paper to overcome these obstacles . unlike cosamp , amop can extract necessary information of target signal from sample data adaptively and operate normally with little prior knowledge . the recovery error of amop is well controlled when random noise is presented and fades away along with increase of snr . moreover , amop has good robustness on detailed setting of target signal and less dependence on structure of measurement matrix . the validity of amop is verified by theoretical derivation . extensive simulation experiment is performed to illustrate the advantages of amop over cosamp in many respects . amop is a good candidate of practical greedy algorithm in various applications of compressed sensing . compressed sensing , greedy algorithm , sparse recovery , noisy environment .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the study of the temperature anisotropies in the cosmic microwave background ( cmb ) can provide crucial information on the origin of the universe and help us to discriminate between different cosmological models . recent observations in many cosmological fields tend to support the standard inflationary model , which predicts the anisotropies to represent an isotropic gaussian random field . the wilkinson microwave anisotropy probe ( wmap ) 1-year data ( bennett et al . 2003 ) , measured the cmb fluctuations with high accuracy . in a first approach ( komatsu et al . 2003 ) , these were found to be consistent with the gaussian predictions . further analyses revealed asymmetries or non - gaussian features which have been confirmed in the 3-year wmap data ( hinshaw et al . 2007 ) . some of these anomalies are : low multipole alignments ( de oliveira - costa et al . 2006 , copi et al . 2007 , land & magueijo 2007 , chiang et al . 2007 ) , north - south asymmetries ( eriksen et al . 2007 ) , structure alignment ( vielva et al . 2007 ) , non - gaussian features found with steerable ( wiaux et al . 2008 ) and directional wavelets ( mcewen et al . 2006 ) , and a very cold spot in the southern hemisphere ( cruz et al . this cold spot ( cs ) is found at galactic coordinates ( @xmath1 ) , and has an angular radius in the sky of @xmath2 . the spherical mexican hat wavelet ( smhw ) was used to detect this structure . this spherical wavelet is an optimal tool to enhance non - gaussian features on the sphere . the first smhw analysis ( vielva et al . 2004 ) was performed on 15 different wavelet scales and two estimators , skewness and kurtosis , were used . the kurtosis deviated significantly from the range of values expected from gaussian simulations at scales around @xmath3 . the spot was found to be the main cause of the deviation from gaussianity , using the area estimator ( cruz et al . cruz et al . ( 2006 ) showed the incompatibility of the frequency independent signal of the cs with galactic foregrounds , and cayn , jin & treaster ( 2005 ) introduced two new estimators , namely @xmath4 and higher criticism , redetecting the cs . cruz et al . ( 2007a ) confirmed the detection in the 3-year wmap data , using all the mentioned estimators . in the latter paper , the probability of finding an at least as high deviation in gaussian simulations was calculated without any _ a posteriori _ assumption . a conservative estimate of this probability was 1.85% , and hence this spot is a very special feature whose origin deserves further investigation . rth et al.(2007 ) detect the cs in the 3-year wmap data using scalar indices instead of wavelets . several models were proposed in order to explain this anomaly : voids ( inoue & silk 2006 , 2007 , rudnick , brown & williams 2007 ) , second order gravitational effects ( tomita 2005 , tomita & inoue 2007 ) , anisotropic bianchi vii@xmath5 model ( jaffe et al . 2005 ) , finite cosmology model , ( adler , bjorken & overduin 2006 ) , asymptotically flat lemaitre - tolman - bondi model ( garcia - bellido & haugbolle 2008 ) and brane - world model ( cembranos et al . 2008 ) are some of them . also a large - scale , non - gaussian angular modulation ( naselsky et al . 2007 ) has been suggested although this hypothesis is not based on any physical model . the bianchi vii@xmath5 model was proven by jaffe et al . ( 2006 ) to be ruled out at high significance , whereas there is still no further evidence for the validity of the other explanations . a full bayesian evidence analysis of that bianchi model can be found in bridges et al . ( 2007a ) . bridges et al . ( 2007b ) show that the cold spot was likely to be driving any bianchi vii@xmath5 template detection . recently , cruz et al . ( 2007b ) showed that the cs could be caused by a cosmic texture . the amplitude and scale of the spot were consistent with that interpretation and the kurtosis of the data was compatible with the gaussian cmb plus texture model at all scales . the predicted number of spots at a scale of @xmath3 or larger , is @xmath6 which is consistent with the observed spot .    in this paper we study two other possibilities , namely the rees - sciama ( rees & sciama 1968 ) effect caused by voids as proposed by inoue & silk ( 2006 , 2007 ) , rudnick , brown & williams ( 2007 ) and the sunyaev - zeldovich effect produced by clusters ( hereafter sz effect , sunyaev & zeldovich 1970 ) already considered in cruz et al . both of them could in principle produce temperature decrements as the cs . in addition a dip ( rudnick , brown & williams 2007 ) in the nrao vla sky survey data ( nvss , condon et al . 1998 ) and the eridanus group of galaxies ( brough et al . 2006 ) match approximately the size and position in the sky of the cs . however , the values required by the parameters of either model are not very realistic . a bayesian framework is very convenient to decide which model best describes the data . the prior theoretical and observational knowledge on the parameters is included in the bayesian analysis and can be decisive in some cases . our analysis is performed on the recently released 5-year wmap data ( hinshaw et al . the large - scale structure of the cmb anisotropies remains almost unchanged respect to the 3-year release ( komatsu et al . 2008 ) . we will first review the data and the template fitting method in section 2 , then we recall in section 3 the results obtained for the texture model in order to compare with the sz effect in section 4 and the void model in section 5 . finally we present our conclusions in section 6 . we use the combined , foreground - cleaned q - v - w map ( hereafter wcm ) of the 5-year wmap data release ( hinshaw et al . 2008 ) . the wcm is in the healpix pixelization scheme ( gorski et al . 2005 ) with resolution parameter nside = 64 . since the scale of the spot is around @xmath3 this resolution is sufficient and reduces the number of pixels used in the template fitting . the template fitting is performed in a circular area of @xmath7 radius centered at galactic coordinates ( @xmath1 ) . although the angular radius of the cs is about @xmath3 we have to consider at least a @xmath7 radius patch to take into account the whole neighborhood of the spot for the fit . the smhw convolves all the pixels in this region and they contribute in an important way to the detected structure .    in our template fitting we apply bayesian statistics to the data , @xmath8 , in order to find the best hypothesis , @xmath9 , describing the data through a set of parameters , @xmath10 . bayes theorem states : @xmath11 where @xmath12 is the posterior distribution , @xmath13 the likelihood , @xmath14 the prior and @xmath15 the bayesian evidence . the null hypothesis , @xmath16 , describes the data as an isotropic gaussian random field ( cmb ) plus the noise of the wmap data . the alternative hypothesis , @xmath17 describes the data as cmb plus noise and an additional template , @xmath18 , given by a physical model as for instance a texture spot , the sz effect produced by a cluster or the rees - sciama effect generated by a local void . @xmath19 ( with @xmath20 ) depends on the parameter set @xmath21 .    for @xmath19 , the probability density associated with the observed data is @xmath22 the evidence is hence the average of the likelihood with respect to the prior and has been extensively used for model selection in cosmology ( e.g. jaffe 1996 , hobson , bridle & lahav 2002 , hobson & mclachlan 2003 , mukherjee , parkinson & liddle 2005 , parkinson , mukherjee & liddle 2006 ) . a model having a large evidence shows a high likelihood for its allowed parameter space given the data . we redefine our notation , with the likelihood @xmath23 , and the prior @xmath24 so the previous integral for the evidence , @xmath25 , reads : @xmath26    the likelihood function is @xmath27 , where @xmath28 is defined as : @xmath29 and @xmath30 is the generalised noise matrix including cmb and noise . the template @xmath18 depends on the parameter set @xmath10 . the calculation of the noise contribution to @xmath30 is straightforward since the noise in the wcm is uncorrelated and well known . in order to obtain the cmb part , we calculate the correlation function for the wcm taking into account the pixel and beam effects . comparing the resulting evidence ratios obtained from the analytical expression for the wcm correlation function with those obtained using simulations , the errors are negligible .    the posterior probability ratio : @xmath31 can be used to decide beween both hypotheses . the alternative hypothesis is the more probable one when @xmath32 . if @xmath33 the null hypothesis is the preferred one to describe the data . the a priori probability ratio for the two models , @xmath34 is set to the fraction of the sky , @xmath35 , covered by either texture , clusters or voids . in this manner we compensate the _ a posteriori _ selection of the pixel where the template is centred , namely the center of the cs , ( @xmath1 ) . unified theories of high energy physics predict the production of topological defects ( kibble 1976 , vilenkin & shellard 2000 ) during a symmetry - breaking phase transition in the early universe . depending on the order of the broken symmetry , different types of defects arise , such as domain walls , strings , monopoles or texture . each of these extremely energetic events leave a characteristic pattern in the cmb . for instance cosmic texture ( turok 1989 ) produce hot and cold spots in the cmb sky . cruz et al . ( 2007b ) showed recently that the cs could be due to a collapsing texture . a bayesian template fit was performed on the 3-year wmap data , using an approximated temperature profile based on the analytic result given in turok & spergel ( 1990 ) . the a priori probability ratio for the two models , @xmath34 , was estimated in this case by the fraction of sky covered by texture , @xmath36 . the evidence ratio was @xmath37 , hence the resulting posterior probability ratio was @xmath38 favoring the texture hypothesis versus the null hypothesis . the scale and amplitude of the spot ( @xmath39 and @xmath40 ) were consistent with having being caused by a texture at redshift @xmath41 and energy scale @xmath42 . the number of texture spots of scale @xmath39 or above is around @xmath43 consistent with the single observed spot . the amplitude of the texture spot was shown to be overestimated by a selection bias . texture spots with amplitude @xmath44 were simulated and added to cmb gaussian simulations . performing a template fit on these single extreme events , the inferred amplitude was in average @xmath45 , i.e. almost twice the true amplitude . this means that the amplitude required by the cs is not in conflict with the upper limits on the amplitude calculated from the power spectrum ( urrestilla et al . 2007 ) . moreover , bevis et al . ( 2008 ) and urrestilla et al . ( 2007 ) find a @xmath46 preference for a cosmological model including texture or strings , when considering only cmb data . however this preference is reduced when including hubble key project ( freedman et al . 2001 ) and big bang nuclesosynthesis ( kirkman et al . 2003 ) constraints in the bayesian power spectrum anlysis .    in order to estimate the significance of the @xmath38 result , the same fitting procedure was performed on 10,000 gaussian cmb simulations . the most prominent spot was selected for each simulation and a texture profile was fitted to it . around @xmath47 of the simulations showed @xmath48 and the typical values where @xmath49 . this result was found to be almost independent on the priors or the exact shape of the approximated temperature profile .    using the 5-year data release we obtain @xmath50 and @xmath51 . the lower noise level allows a slightly better detection of the texture template in the new data release . we have not repeated the frequentist estimation of the significance because the result was very similar to the 3-year one . the probability of having an at least as extreme value ( i.e. @xmath52-value ) of @xmath53 assuming that the null hypothesis is true , should be again around @xmath54 . the presence of texture would explain the excess of kurtosis observed at scales around @xmath3 . the @xmath55 texture predicted by the model would be confused with gaussian cmb spots since the latter have their maximum power at these scales . a number of follow - up tests can be performed to check the texture hypothesis . for instance an all - sky template fitting looking for several texture can confirm or refute the alternative hypothesis . also a polarization analysis can help to decide , since gaussian cmb spots have a correlated polarization pattern , whereas late texture do not show associated polarization . in addition the b - mode texture power spectrum is expected to be larger than the typical inflationary one , having also a different shape ( urrestilla et al . the hot gas trapped in clusters modifies the cmb spectrum through the sunyaev - zeldovich effect . the high - energetic electrons present in the gas collide with the cmb photons , changing their original frequency . two different components can be identified . the former is produced by random thermal motion of the electrons ( thermal effect ) , and the second one , by the peculiar motion of the cluster with respect to the cmb ( kinetic effect ) . the thermal sz effect is typically one order of magnitude higher than the kinetic one , and therefore we focus on the first one . the effect on the cmb photons is frequency dependent and at frequencies below 217 ghz it produces a temperature decrement . hence this effect is a possible candidate to explain an anomalously cold spot at the wmap frequencies which range from 23 ghz to 94 ghz . the temperature decrement is given by : @xmath56 where @xmath57 is the frequency in dimensionless units , with @xmath58 and @xmath59 planck and boltzmann constants respectively , and @xmath60 . the cluster comptonization parameter is defined as @xmath61 where @xmath62 and @xmath63 are the intracluster electron density and temperature , @xmath64 is the thomson cross section , and @xmath65 is the electron mass . the integral is performed along the line of sight through the cluster . in hot clusters ( @xmath66 kev ) , the relativistic electrons lead to a slightly modified spectral behaviour of the thermal sz effect , but this is negligible in our analysis . also neglected is the kinetic sz effect of the cluster . cruz et al . ( 2006 ) showed that the flat frequency dependence of the cs is incompatible with being caused by a sz signal alone . however a combination of cmb plus sz effect may explain the spot and could have a sufficiently flat frequency spectrum . cruz et al . ( 2005 ) already reported that a group of local galaxies were present in the direction of the cs . they argued that the amount of gas needed to explain the whole signal was too large for a small group of galaxies . however the group of galaxies in the direction of the spot has been shown recently to be larger than previously thought . brough et al . ( 2006 ) , define it as the eridanus _ supergroup _ with mass @xmath67 . the angular size and position in the sky of eridanus are similar to those of the cs , as can be seen in figure  ( [ fig : patch ] ) . the top panel shows the wcm data and the bottom one , the nearby galaxies catalog ( tully , 1987 ) . both are @xmath68 patches centered at ( @xmath1 ) and were convolved with the smhw at scale @xmath3 . patch of the wmap combined map ( wcm , top ) and the nearby galaxies catalog ( bottom ) , centered at ( @xmath1 ) and convolved with the smhw at scale @xmath3 . the units are wavelet coefficients of temperature in mk and number of counts for the wcm and nearby galaxies catalog respectively.,width=317,height=566 ]    however there are other , more prominent nearby groups which do not have an associated cmb spot , hence this could be a statistical coincidence , but it is worth testing through a template fit whether there is a hidden sz - like template in the wcm data .    the temperature profile as a function of the angle in the sky , @xmath69 , can be calculated assuming a @xmath70-model ( with @xmath71 ) for the gas - density profile ( see e.g. atrio - barandela & mcket , 1999 ) : @xmath72 where @xmath73 is the central electron density and @xmath74 is the core radius of the cluster . expressing the virial radius of the cluster @xmath75 in terms of the core radius as @xmath76 , eq .  [ [ yc ] ] can be rewritten as : @xmath77 where @xmath78r_c t n_c$ ] , and @xmath79 . \label{phi}\ ] ] the parameter @xmath80 is the angle subtended by the core radius of the cluster , whereas @xmath52 is not a parameter but a function of mass @xmath81 and redshift @xmath82 . typical values for clusters yield @xmath83 . summarising , our temperature profile with parameters @xmath80 and @xmath84 reads : @xmath85 . \label{dttfinal}\ ] ] eq .  [ [ dttfinal ] ] gives the effect produced by a single cluster through a particular line of sight . as the wcm is a combination of maps at different frequencies , we choose the central frequency , @xmath86 ghz . we consider two alternative hypotheses in this context . in the first one the spot would be caused by a rich cluster and in the second one by a nearby group of the nearby galaxies catalog . we simulate the entire population of clusters with flux above 10 mjy ( at 353 ghz ) , using a monte carlo algorithm where masses and redshifts were assigned to each cluster according to the press - schechter prescription ( press & schechter 1974 ) . we assume @xmath87 , @xmath88 and @xmath89 ( where @xmath90 and @xmath91 are the matter and vacuum energy density parameters , and @xmath92 the dispersion of the density field at scales of 8 mpc ) . the cut in flux removes nearby small groups of galaxies and distant small clusters . the redshifts of the simulated clusters range from 0.005 to 3.4 . for each cluster with known mass and redshift the core radius is determined from the scaling relation : @xmath93 where @xmath94 is the mass of the cluster in @xmath95 @xmath96 . it is important to note that the cut in flux does not have any negative effect in the results of this paper since the fluxes we are trying to explain are much larger than 10 mjy . since the redshift and core radius are known for each cluster , we calculate the angle subtended by the core , @xmath80 . we find the maximum angle of the 300,000 clusters to be @xmath97 . therefore the cluster hypothesis has to be discarded for generating a @xmath3 spot . returning to the particular case of the eridanus supergroup found in the nearby galaxies catalog , we can consider larger values of @xmath80 , due to its proximity and lower values of @xmath84 since it is not a cluster but a big group or an ensemble of groups . in addition we assume lower values of @xmath52 because the gas could be rather dispersed since the gravitational attraction is lower than in big clusters . now the alternative hypothesis is that a big group from the nearby galaxies catalog is causing the cs . although the center of eridanus does not match exactly the center of the spot , we still center the template at ( @xmath1 ) , since the gas does not need to fit exactly the positions of the galaxies . this assumption is the most favorable for the alternative hypothesis . if the cs was due to the sz effect we would expect high x - ray emission in the considered region . we can use the rosat all - sky survey ( snowden et al . 1997 ) to estimate the prior on @xmath84 . the rosat survey excludes some small regions of the sky , and one of them crosses the considered patch . however the estimated average counts in the rosat r6 band are around @xmath98 cts s@xmath99 arcmin@xmath100 . the @xmath70-models which are compatible with this background and at the same time show an extended sz signal of several degrees , range between models with a core radius of @xmath101 mpc and bolometric luminosities of @xmath102 erg / s and models with @xmath103 mpc and bolometric luminosities of @xmath104 erg / s . these models would predict at most @xmath105 . the nearest groups in the nearby galaxies catalog are at @xmath106 mpc distance and as we have shown , the upper limit for the core radius is @xmath103 mpc . in this case we get @xmath107 rad , which we set as the upper limit for @xmath80 . we set the lower limit on @xmath80 to @xmath108 , since the spot has @xmath3 radius and the pixel size is around @xmath108 for the chosen resolution . hence we consider flat priors for parameter @xmath84 in the interval @xmath109 $ ] and for @xmath80 in @xmath110 $ ] radians . setting @xmath111 we obtain a relative evidence of @xmath112 . the fraction of the sky covered by nearby galaxies in the catalog , is @xmath113 therefore @xmath114 and the sz hypothesis is discarded again . the best fit parameters are the upper limits of the priors , @xmath115 and @xmath107 rad . the amplitude of the temperature decrement does not exceed @xmath116 and is hence negligible compared to the deep decrement . therefore the @xmath112 result is not suggesting that the sz hypothesis is viable for explaining the spot . as the template has a negligible effect on the data , @xmath117 , the likelihood is almost constant in the whole parameter range and as the priors integrate to unity we have @xmath118 . this is a well - known feature of bayesian model selection . in any problem , if @xmath17 has additional parameters compared to @xmath16 , and the likelihood is insensitive to these parameters , then @xmath119 . the results do not significantly depend on the value of @xmath52 , the chosen priors or the resolution . the gravitational effect due to a supercluster can not strongly affect large angular fluctuations because these massive structures are very concentrated in a particular region . on the contrary , voids are better candidates to generate large angular spots . the gravitational effect produced by cosmic voids can produce temperature decrements in the cmb and has been proposed as a possible source of the cs ( inoue & silk 2006 , 2007 ) . in analogy to the cluster case with the eridanus supergroup , there is a dip in the nvss radio source catalogue that approximately matches the position of the spot . rudnick , brown & williams ( 2007 ) , claim that the dip could be associated with a void at redshift @xmath120 which could be the origin of the cs , but affirm that the significance of this coincidence is unclear . as the fraction of the universe filled by voids is around @xmath121 in volume at redshift @xmath122 ( colberg et al . 2005 ) , the fraction of the sky covered in projection , @xmath35 , is close to one and hence such coincidences are not rare . the dip and the cs appear at different scales as can be seen in figure [ fig : wmap_nvss ] where we present a @xmath68 patch of the wmap and nvss data convolved with the smhw at scales @xmath3 and @xmath123 . at scale @xmath3 the cs reaches its maximum power in the wcm , whereas there is no significant dip visible in the nvss data . at scale @xmath123 where the cs is not an outstanding feature , a dip becomes visible but there are other similar dips with no associated cmb spots , therefore the association between dip and spot is unclear and could be a coincidence . moreover , there are several other dips with similar depth outside the selected region hence the dip does not seem to be an outstanding feature of the nvss data . mcewen et al . ( 2007 ) found several regions where both data sets were highly correlated and one of them was located close to the cs . this could be due to the very cold temperature of the cs which gives a high correlation with any underdense region in the nvss data at scale @xmath3 and not to the relation between dip and spot since they have different scales ( see figure [ fig : wmap_nvss ] ) . patch of the wcm ( left column ) and the nvss ( right column ) , centered at ( @xmath1 ) and convolved with the smhw at scale @xmath123 ( top row ) and scale @xmath3 ( bottom row ) . the units are wavelet coefficients of temperature in mk and radio source number counts for the wcm and nvss respectively.,width=317 ]    rudnick , brown & williams ( 2007 ) also performed an order of magnitude calculation to estimate the radius of the void , @xmath75 , needed to generate such a cmb spot through the rees - sciama effect obtaining @xmath124 mpc in concordance with previous calculations ( inoue & silk 2007 ) . the probability of observing such a void in standard theory is very small . for instance , a void with @xmath125 mpc will occur with a probability of @xmath126 in an optimistic case ( rudnick , brown & williams 2007 ) . typical observed values are @xmath127 mpc for a matter density contrast of @xmath128 ( hoyle & vogeley 2004 , patiri et al . similar results are found in n - body simulations ( colberg et al . 2005 ) . however it is worth performing a bayesian template fit for the void hypothesis . martnez  gonzlez et al . ( 1990 ) and martnez  gonzlez & sanz ( 1990 ) , performed a calculation of the temperature anisotropy generated by a compensated void in an einstein - de sitter universe , @xmath129 where @xmath130 is the angle formed between the direction of observation and the direction defined by the centre of the structure , and @xmath131 is the distance to the centre of the structure . they adopt the thin shell approximation , where the void is an empty spherical region surrounded by a thin shell containing the matter which would have been inside the void . outside the void the universe is represented by an einstein - de sitter model and the shell propagates approximately as @xmath132 and @xmath133 . the units are @xmath134 , where @xmath135 is the horizon distance . [ [ eq : dt_t_void ] ] applies to photons crossing the void . the photons that do not cross the void are unaffected since there is no net mass excess . introducing the above expressions for @xmath131 and @xmath141 in eq .  [ [ eq : dt_t_void ] ] we have our temperature template with parameters @xmath82 , @xmath75 in a flat @xmath136-cdm universe . we choose the prior on @xmath75 in concordance with the void radius distribution found in the two - degree field galaxy redshift survey data ( fig . 3 in hoyle & vogeley 2004 ) . we extrapolate these values for large @xmath75 through an exponential decay as shown in figure  ( [ fig : pr_rv1 ] ) . this extrapolation is generous with big voids since the derived decay is faster than the exponential used .        for the redshift , we choose a flat prior between 0 and 1 . this assumption is again rather generous with large radii , since according to n - body simulations ( colberg et al . 2005 ) the void sizes should be smaller at high , unobserved redshifts . we set the a priori probability ratio for the two models to unity since voids cover almost all the sky . the obtained evidence ratio is @xmath118 and therefore the posterior probability ratio is @xmath142 . this result is not suggesting that the void model is viable as the best fit template corresponds to @xmath143 mpc and @xmath144 which produces a temperature anisotropy of @xmath145 . as in the cluster case the template is completely negligible , the likelihood almost constant , and the priors integrate to unity , therefore @xmath118 . the contribution of voids with @xmath146mpc to the evidence integral is negligible and the results remain unchanged varying the upper limit of parameter @xmath75 for values above @xmath147mpc . the big voids ( @xmath148mpc and @xmath149 ) described in rudnick , brown & williams ( 2007 ) which do produce an important temperature decrement are strongly penalized by the prior , despite showing higher likelihood values . we redefine @xmath17 in order to test the big void hypotesis . now only voids with @xmath148mpc and @xmath150 are considered . we renormalise the priors in the new parameter ranges and recalculate the prior probability ratio since the big voids ( in case they exist ) are expected to cover only a small part of the sky . as shown above a void with @xmath125 mpc is expected to occur with a probability of @xmath126 ( rudnick , brown & williams 2007 ) . hence as a valid upper limit we can assume that there is at most one big void with @xmath125mpc and @xmath151 . the fraction of the sky covered by this void is @xmath152 . therefore we assume @xmath153 . recalculating the evidence for the new parameter range , we obtain @xmath154 which yields @xmath155 . this result discards the possibility of explaining the cs through such a void . hence the void hypothesis has to be rejected . inoue & silk ( 2006 , 2007 ) proposed slightly different voids with very high values for the radius of the void , @xmath156 mpc and quasi linear density contrast values @xmath157 , while rudnick , brown & williams ( 2007 ) used @xmath158 . inoue & silk ( 2007 ) recognize that the required values for @xmath75 are inconsistent with observations at low @xmath82 and n - body simulations . such a void is a @xmath159 result in the standard theory with @xmath160 . they consider that a percolation process of the typically observed small voids could form the large quasi linear voids . however these voids have not been observed yet and as already mentioned , according to n - body simulations ( colberg et al . 2005 ) the void sizes should be smaller at high , unobserved redshifts . in this paper we consider different possibilities to explain a non - gaussian cold spot found in the wmap data . this spot is located at galactic coordinates ( @xmath1 ) , covers around @xmath161 in the sky , has been shown to be incompatible with gaussian simulations and is unlikely to be produced by foregrounds ( vielva et al . 2004 , cruz et al . 2005 , cruz et al . 2006 , cruz et al . we compare three alternative hypotheses to explain the cs , namely texture , clusters and voids versus the null hypothesis . the latter describes the data as an isotropic gaussian random field ( cmb ) plus the noise of the wmap data . a bayesian analysis takes into account the prior observational and theoretical knowledge on the parameters and can be used to find the hypothesis best describing the data . the posterior probability ratio , @xmath53 , is the quantity we use to compare the different hypotheses . the alternative hypothesis is preferred if @xmath32 and disfavored otherwise . cruz et al . ( 2007b ) performed such an analysis for the texture model , finding that this model is favored versus the hypothesis of gaussian fluctuations with @xmath162 . the predicted size , amplitude , temperature profile , and number of spots generated by texture is consistent with the observed spot . repeating this analysis for the 5-year wmap data we find @xmath163 due to the lower noise level . the frequentist @xmath52-value associated to this result is around @xmath54 . the finding of a local supergroup ( brough et al . 2006 ) and a dip in the nvss survey ( rudnick , brown & williams 2007 ) , which match the position of the cs could be a hint for the cluster or void hypotheses , but our detailed bayesian analysis reveals that these are statistical coincidences . the fraction of the sky covered by clusters and voids is included in our analysis . the large angular size and temperature decrement of the cs are hard to explain with conventional clusters or voids . we would need either a local , rich cluster or a void with radius @xmath164 mpc located at redshift @xmath149 . in the cluster case only a local group with compton parameter @xmath165 is observed in the required direction . the angular scale is marginally consistent with the spot but the temperature decrement generated is only a few @xmath166 and hence can not account for the cs . our bayesian analysis gives @xmath167 , discarding the cluster hypothesis .    for conventional voids with @xmath127 mpc we obtain @xmath142 because the likelihood is insensitive to the extra parameters . a detailed analysis reveals that these voids produce negligible temperature decrements below @xmath116 . renormalising the priors for voids with radius @xmath164 mpc , we obtain a posterior probability ratio of @xmath168 . therefore the void hypothesis is not viable either . explaining the cold spot by a void with @xmath148 mpc , would introduce a @xmath159 anomaly in the standard model to solve the @xmath169 anomaly of the spot . the solution would be worse than the problem . hence we can conclude that the texture hypothesis is the most plausible explanation for the cs . many scientists are uncomfortable with cosmic defects , because these failed to explain structure formation . however , cosmic defects are a generic consequence of symmetry breaking , which is an essential ingredient in unified theories of particle physics . hence it is very plausible that they exist . they might be formed at any scale below the grand unified scale , and would be consistent with current cosmological observations . they are even compatible with inflation . the only constraint is that they should make a subdominant contribution to structure formation ( see for instance urrestilla et al . 2007 ) . the virtue of the texture hypothesis is that it is fully testable . a texture would not lead to a correlated polarization pattern around the cold spot , but it would lead to a very specific gravitational lensing pattern . in addition there should be more cmb texture spots at scales around @xmath170 . although they are more difficult to detect because the gaussian cmb reaches its maximum power at these scales , we may be able to detect them through an all - sky bayesian analysis , which is already in progress . mc , emg and pv acknowledge financial support from the spanish mcyt project aya2007 - 68058-c03 - 02 . pv acknowledges financial support from the ramn y cajal project . we acknowledge the use of the legacy archive for microwave background data analysis ( lambda ) . support for lambda is provided by the nasa office of space science . this work has used the software package healpix ( hierarchical , equal area and iso - latitude pixelization of the sphere , http://www.eso.org/science/healpix ) , developed by k.m . grski , e. f. hivon , b. d. wandelt , j. banday , f. k. hansen and m. barthelmann ; the visualisation program univiewer , developed by s.m . mingaliev , m. ashdown and v. stolyarov ; and the camb and cmbfast software , developed by a. lewis and a. challinor and by u. seljak and m. zaldarriaga respectively . adler r. j. , bjorken j. d. , overduin j. m. , 2006 , gr - qc/0602102 . atrio - barandela f. , mcket j.  p. , 1999 , apj 515 , 465 . garcia - bellido j. , haugbolle t. , 2008 , jcap , 04 , 003 . bennett c.l . , et al . , 2003 , apjs , 148 , 1 . bevis n. , hindmarsh m. , kunz m. , urrestilla j. , 2008 , phys . lett . , 100 , 021301 . bridges m. , mcewen j.d . , lasenby a.n . , hobson m.p . , 2007 , mnras , 377 , 1473 bridges m. , mcewen j.d . , cruz m. , hobson m.p . , lasenby a.n . , vielva p. , martnez  gonzlez e. , 2008 , submitted to mnras ( arxiv:0712.1789 ) . brough , s. , forbes d.  a. , kilborn v.  a. , couch w. , colless m. , 2006 , mnras , 369 , 1365 . cayn l. , jin l. & treaster a. , 2005 , mnras , 362 , 826 . cembranos j. a. r. , de la cruz - dombriz a. , dobado a. , maroto , a. l. , 2008 , ( arxiv:0803.0694 ) condon j.j . , cotton w.d . , greisen e.w . , yin q.f . , perley r.a . , taylor g.b . , broderick j.j . , 1998 , aj , 115 , 1693 . chiang l .- y . , naselsky p.  d. , coles p. , 2007 , apj , 664 , 8 . cruz m. , martnez  gonzlez e. , vielva p. , cayn l. , 2005 , mnras , 356 , 29 . cruz m. , tucci m. , martnez  gonzlez e. , vielva p. , 2006 , mnras , 369 , 57 . cruz m. , cayn l. , martnez  gonzlez e. , vielva p. , jin j. , 2007a , apj , 655 , 11 . cruz m. , turok n. , vielva p. , martnez  gonzlez e. , hobson m. , 2007b , science , 318 , 1612 . colberg j.m . , sheth r.k . , diaferio a. , liang g. , yoshida n. , 2005 mnras , 360 , 216 . copi c. j. , huterer d. , schwarz d. j. , starkman g. d. , 2007 , phys . rev . d 75 , 023507 . de oliveira - costa a. , tegmark m. , 2006 , phys . rev . d 74 , 023005 . eriksen h. k. , banday a. j. , grski k. m. , hansen f. k. , lilje p. b. , 2007 , apj , 660 , 81 . freedman et al . , 2001 , , apj , 553 , 47 . hinshaw et al . , 2007 , apjs , 170 , 288 . hinshaw et al . , 2008 , submitted to apj , ( arxiv:0803.0732 ) hobson m , bridle s. , lahav o. , 2002 , 335 , 377 . hobson m. , mclachlan , 2003 , mnras , 338 , 765 . m. s. , 2004 , apj , 607,751 . inoue k. t. , silk j. , 2006 , apj , 648 , 23 . inoue k. t. , silk j. , 2007 , apj , 664 , 650 . jaffe a. , 1996 , apj , 471 , 24 . jaffe t. r. , banday a. j. , eriksen h. k. , grski k. m. , hansen f. k. , 2005 , apj , 629 , 1 . jaffe t. r. , hervik s. , banday a. j. , grski k. m. , 2006 , apj 644 , 701 . kibble , j. phys . , 1976 , a9 , 1387 . kirkman d. , tytler d. , suzuki n. , omeara j.m . , lubin d. , 2003 , apjs , 149 , 1 . komatsu e. et al . , 2003 , apjs , 148 , 119 . komatsu e. et al . , 2008 , submitted to apj , ( arxiv:0803.0547 ) . land k. , magueijo j. , 2007 , mnras.378 , 153 . gonzlez e. & sanz j. l. , 1990 , mnras , 247 , 473 . martnez  gonzlez e. , sanz j. l. & silk j. , 1990 , apj , 355 , 5 . mcewen j.d . , hobson m.p . , lasenby a.n . , mortlock d.j . , 2006 , mnras , 371 , l50 mcewen j.d . , vielva p. , hobson m.p . , martnez  gonzlez e. , lasenby a.n . , 2007 , mnras , 376 , 1211 . mukherjee p. , parkinson d. , liddle a. , 2005 , apj , 635 , 51 . naselsky p.  d. , christensen p.  r. , coles p. , verkhodanov o. , novikov d. , kim j. , 2007 , apj submitted , arxiv:0712.1118 . parkinson d. , mukherjee p. , liddle a. , 2006 , phys rev . d 73 , 123523 . patiri s.g . , betancort - rijo j. , prada f. , klypin a. , gottlber s. 2006 , mnras , 369 , 335 . press w.h , schechter p. , 1974 , apj , 187 , 425 . rth c. , schuecker p. , banday a.  j. , 2007 , mnras , 380 , 466 . rees m.j . , . , 1968 , nature , 217 , 511 . rudnick l. , brown s. , williams l. , 2007 , apj accepted , astro - ph/0704.0908 . tomita k. , 2005 , phys . d , 72 , 10 . tomita k. , inoue k.t . , 2007 , arxiv:0712.1291 . tully r. b. , 1988 , nearby galaxies catalog . cambridge university press , 1988 , 221 p. turok n. & spergel d. n. , 1990 , phys . letters , 64 , 2736 . urrestilla j. , bevis n. , hindmarsh m. , kunz m. , liddle a. , 2007 , arxiv:0711.1842 . vielva p. , martnez - gonzlez e. , barreiro r.  b. , sanz j.l . , cayn l. , 2004 , apj , 609 , 22 . vielva p. , wiaux y. , martnez - gonzlez e. , vandergheynst p. , 2007 , mnras , 381 , 932 . vilenkin a. , shellard e.p.s . , 2000 , cosmic defects and other topological defects ( cambridge univ . press ) . wiaux y. , vielva p. , barreiro r.  b. , martnez - gonzlez e. , vandergheynst p. , 2008 , mnras , 385 , 939 .
Abstract:

Answer: the non - gaussian cold spot found in the wmap data has created controversy about its origin . here we calculate the bayesian posterior probability ratios for three different models that could explain the cold spot . a recent work claimed that the spot could be caused by a cosmic texture , while other papers suggest that it could be due to the gravitational effect produced by an anomalously large void . also the sunyaev - zeldovich effect caused by a cluster is taken into account as a possible origin . we perform a template fitting on a @xmath0 radius patch centered at galactic coordinates ( @xmath1 ) and calculate the posterior probability ratios for the void and sunyaev - zeldovich models , comparing the results to those obtained with texture . taking realistic priors for the parameters , the texture interpretation is favored , while the void and sunyaev - zeldovich hypotheses are discarded . the temperature decrement produced by voids or clusters is negligible considering realistic values for the parameters . = = = = = = = =    # 1 # 1 # 1 # 1 @mathgroup@group @mathgroup@normal@groupeurmn @mathgroup@bold@groupeurbn @mathgroup@group @mathgroup@normal@groupmsamn @mathgroup@bold@groupmsamn =  019 =  016 =  040 =  336 =  33e = = = = = = = =    # 1 # 1 # 1 # 1 = = = = = = = =    methods : data analysis - cosmic microwave background