Concentration of Measure Inequalities for Compressive Toeplitz Matrices with Applications to Detection and System Identification
|
|
- Pierce Shields
- 6 years ago
- Views:
Transcription
1 Concentration of Measure Inequalities for Compressive Toeplitz Matrices with Applications to Detection an System Ientification Borhan M Sananaji, Tyrone L Vincent, an Michael B Wakin Abstract In this paper, we erive concentration of measure inequalities for compressive Toeplitz matrices having fewer rows than columns with entries rawn from an inepenent an ientically istribute ii Gaussian ranom sequence These inequalities show that the norm of a vector mappe by a Toeplitz matrix to a lower imensional space concentrates aroun its mean with a tail probability boun that ecays exponentially in the imension of the range space ivie by a factor that is a function of the sample covariance of the vector Motivate by the emerging fiel of Compressive Sensing CS, we apply these inequalities to problems involving the analysis of high-imensional systems from convolution-base compressive measurements We iscuss applications such as system ientification, namely the estimation of the impulse response of a system, in cases where one can assume that the impulse response is high-imensional, but sparse We also consier the problem of etecting a change in the ynamic behavior of a system, where the change itself can be moele by a system with a sparse impulse response I INTRODUCTION We live in an era that one might call the century of ata explosion inee, we have now reache a point where all of the ata generate by intelligent sensors, igital cameras, an so on excees the worl s total ata storage capacity [] Motivate to reuce the burens of acquiring, transmitting, storing, an analyzing such vast quantities of ata, signal processing researchers have over the last few ecaes evelope a variety of techniques for ata compression an imensionality reuction Unfortunately, many of these techniques require a raw, high-imensional ata set to be acquire before its essential low-imensional structure can be ientifie, extracte, an exploite In contrast, what woul be truly esirable are sensors that require fewer raw measurements yet still capture the essential information in a ata set One positive outcome of this past work, however, has been the evelopment of sparse an compressible representations as concise moels for high-imensional ata sets Builing on these principles, Compressive Sensing CS has very recently emerge as a powerful paraigm for combining measurement with compression First introuce by Canès, Romberg an Tao [] [4], an Donoho [5], the CS problem can be viewe as recovery of an s-sparse signal a R n from just < n or even n observations y Xa R, where X R n is a matrix representing All authors are with Division of Engineering, Colorao School of Mines, Golen, CO 8040, USA {bmolazem,tvincent,mwakin}@mineseu This work was partially supporte by AFOSR Grant FA , NSF Grant CCF , DARPA Grant HR , NSF Grant CNS an Department of Energy, Office of Energy Efficiency an Renewable Energy Grant DE-FG36-08GO8800 the linear measurement process An s-sparse signal a R n is a signal of length n with s < n nonzero significant entries The notation s : a 0 enotes the sparsity level of a Since the null space of X is non-trivial, exact recovery of the signal epens on whether the solution set for y Xa contains only one solution of the require sparseness It has been shown that if X is a fully populate ranom matrix with entries rawn from an inepenent an ientically istribute ii Gaussian ranom sequence, then exact recovery of a can be guarantee with high probability from O s log n s measurements [6] In aition, recovery is possible by solving a convex optimization problem an is robust to measurement noise [7] As an extension of CS, there is also recent work on the etection of sparse signals from compressive measurements [8] Concentration of measure inequalities are one of the leaing techniques use in the theoretical analysis of ranomize compressive operators [9] A typical concentration inequality takes the following form [0] For any fixe signal a R n, an a suitable ranom matrix X, the norm of the projecte signal by X will be highly concentrate aroun the norm of the original signal with high probability In other wors, there exist constants c an c, such that for any fixe a R n, P { Xa E [ Xa ] E [ Xa ]} c e cc0, where c 0 is a function of 0, In this paper, we erive concentration of measure inequalities for ranom, compressive Toeplitz matrices an iscuss applications to etection an system ientification A Compressive Toeplitz Matrices While the recovery results iscusse in the previous section are for a measurement matrix X with inepenent elements, many applications will require X to have a particular structure In particular, when ynamical systems are involve, the measurement process may involve the convolution of an unknown system impulse response with a known input Consier ientifying the impulse response of a linear time-invariant LTI system from its input-output observations Let {x k } n+ k be the applie input sequence to an LTI system characterize by its finite impulse response {a k } n k Then the corresponing output y is calculate from the time-omain convolution y x a Consiering the x k an a k sequences to be zero-pae from both sies, each
2 output sample y i can be written as n y i a j x i j j If we only keep observations of the system, {y k } n+ kn+, then can be written in a matrix-vector multiplication format as y Xa, where x n x n x x n+ x n x X x n+ x n+ x is a n Toeplitz matrix Supposing the system has an s-sparse impulse response, we are intereste in efficiently acquiring an retrieving a via its multiplication by the n compressive Toeplitz matrix X We call this a compressive Toeplitz matrix because < n In aition to the application to convolution measurement problems, it is worth noting that while fully populate ranom matrices require n elements to be generate, Toeplitz matrices require only n + istinct ranom entries, which may provie an avantage in other CS applications B Relate Results Compressive Toeplitz matrices an convolution have been previously stuie in the context of compressive sensing in [] [7], with applications incluing channel estimation an synthetic aperture raar The work in this paper is most closely relate to [7], which also erives a concentration of measure boun for compressive Toeplitz matrices an utilizes this boun to stuy recovery A comparison between our work an this work is given in the next section, with further iscussion of relate work postpone to later sections C Main Result In this paper, we erive a concentration of measure boun for compressive Toeplitz matrices as given in 3 with entries rawn from an ii Gaussian ranom sequence Our main result, etaile in Theorem, states that the upper an lower probability tail bouns epen on the number of measurements an on the eigenvalues of the measurement covariance matrix P P a efine as: R a 0 R a R a R a R a 0 R a P, 4 R a R a R a 0 where n τ R a τ : a i a i+τ enotes the un-normalize sample autocorrelation function of a R n 3 Theorem : Let a R n be fixe Define two quantities ρa an µa associate with the eigenvalues of the measurement covariance matrix P P a as: an ρa max i λ i a µa λ i a 4 where {λ i } are the eigenvalues of P Suppose X is a ranom compressive Toeplitz matrix with ii Gaussian entries having zero mean an unit variance Then for any 0,, the upper tail probability boun is P { y > a + } + e an the lower tail probability boun is P { y < a ] } e ρa 7 µa 8 In [7], an upper an lower boun essentially ientical to 7 is given, but only for a range of boune away from 0 Our result 8 an proof technique appear to be new D Organization The rest of the paper is organize as follows In Section II, we state the proof of the main result In Section III, we provie bouns on the spectral norm of the covariance matrix P A system ientification application is consiere in Section IV In Section V, we explore the problem of etection with compressive sensing an present experimental results We conclue in Section VI II PROOF OF MAIN RESULT Our analysis begins with the observation that, for fixe a an ranom Gaussian X, y Xa will be a Gaussian ranom vector Lemma : If {x k } n+ k is a zero mean, unit variance ii Gaussian ranom sequence, then y Xa is a Gaussian ranom vector with zero mean an covariance matrix P P a given in 4 Proof: See Appenix A From Lemma, it quickly follows that E [ y ] a It is also easy to see that P is a symmetric Toeplitz matrix with tr P a where {λ i P } enote the eigenvalues of P Note that λ i P 0 for all i,,, The proof of the main theorem utilizes Markov s inequality along with a suitable boun on the moment generating function of y This boun epens on the following lemma Lemma : If y R is a Gaussian ranom vector with covariance matrix P, then E [ e ±ty y ] et I tp, 9 where we require t 0, max i λ i for the case involving et I tp, an y enotes the transpose of y
3 Proof: [ ] E e ±ty y e ±ty y e π et P y P y y e π et P y P tiy y P et ti et P et P ti et P et I tp Remark : As a special case of Lemma, if y R is a scalar Gaussian ranom variable of unit variance, then we obtain the well known result an E [ e ±ty] t 0 The main objective of this paper is to fin bouns on P { y > a + } a P { y < a } b We use Markov s inequality an Chernoff s bouning metho for computing the upper an lower probability tail bouns a an b Specifically, for a ranom variable Z, an all t > 0, P {Z > } P { e tz > e t} E [ e tz] e t see eg [8] Applying to a yiels [ ] P { y > a + } E e ty y 3 e a +t Using Lemma, we rewrite the right han sie of 3 as [ ] E e ty y e et I tp a +t e a +t 4 In 4, t > 0 is a free variable which can be varie to fin the tightest possible boun, resulting in the Chernoff boun approximation of the inequality Claim : The value of t that minimizes 4 satisfies λ i a tλ + 5 i Proof: See Appenix B Unfortunately, 5 is ifficult to use for further analysis However, since 4 is vali for all t, we are free to choose a suboptimal, but still useful value We propose to use t + fa a, where f is a function of a that will be chosen later in a way that leas to tighter bouns Now we are reay to prove the main theorem We state the upper tail probability boun in Lemma 3 an the lower tail probability boun in Lemma 4 Lemma 3: Let y Xa R be a zero mean Gaussian ranom vector with covariance matrix P, where X is a n compressive Toeplitz matrix populate with ii zero mean, unit variance Gaussian ranom variables, an a is an s-sparse vector of length n Then, for any 0,, P { y > a + } + e ρa 6 Proof: From 3 an 4, we get P { y > a + } et I tp e a +t 7 Choosing t as t + ρ a a, 8 the right han sie of 7 can be written as et I P + ρ a a e ρa 9 This expression can be simplifie by applying the previous lemmas Note that P λ i et I + ρ a a + ρ a a e log λ i + ρa a Using the fact that log c c c log c for any c, c [0, ] an tr P a, we have e log + λ i ρa a e e trp ρa a λ i ρa a log + e log ρa + + Combining 7, 9, an 0 gives us P { y > a + } + log + ρa 0 ρa e + e ρa ρa Lemma 4: Using the same assumptions as in Lemma 3, for any 0,, P { y < a } e 3 µa +
4 Proof: obtain Applying Markov s inequality to b, we P { y < a } P { y > a } [ ] E e ty y e a t Using Lemma, this implies P { y < a } et I + tp e a t In this case, we choose t + µa a Plugging t into an following similar steps as for the upper tail boun, we get P eti + tp et I + + µa a λ i + + µa a e + log λ i + µa a 3 Since log + c c c for c > 0, λ i log + + µa a λ i + µa a λ i + µa a Thus, λ i log + + µa a λ i + µa a λ i + µa a λ i + µa a + µa a + µa + µa + µa + Using this in 3 gives the boun eti + tp e + µa + e + µa + 4 λ i f upper ε f lower ε ε Fig Comparison between upper an lower tail boun functions By substituting 4 in, we obtain P { y < a } e + e + e 3 µa + µa III CONCENTRATION BOUNDS FOR SPARSE SIGNALS As can be seen the statement of Theorem, the upper an lower tail probability bouns are functions of the quantities ρ a an µ a For any 0,, we have the upper tail probability boun P { y > a + } f upper an the lower tail probability boun where P { y < a } f lower f upper + e ρa µa, an f lower e 3 + Figure shows a comparison between f upper an f lower for 0, These functions are similar, an behave as e 4 for small Typically, we are seeking to fin values of for which the boun is small for a particular class of signals, which in our case will be s-sparse signals Consequently, it is necessary to obtain an upper boun for the quantities ρ a an µ a for the class of signals of interest It is easy to show that for all a, µ a ρ a Therefore, we limit our problem to fining the sharpest boun for ρ a Note that P is a symmetric Toeplitz matrix which can be ecompose as: P A A,
5 where A is an n + Toeplitz matrix forme by a, a a n a A 0 0, a n a a n an A is the transpose of A By the Cauchy Interlacing Theorem [9], max λ i P max λ i A A max λ i à Ã, where à is an n + n + circulant matrix with ã [a, a,, a n, 0,, 0], as its first row Since λ i à à λ i Ã, an upper boun for ρ a is provie by the maximum eigenvalue of à Since à is circulant, λ i à simply equals the un-normalize n + -length Discrete Fourier Transform DFT of the first row of à As a result, λ i à n k a k e jπi k/n+ 5 When a is an s-sparse signal, only s < n + terms in this summation are nonzero, therefore λ i à s a This results in the boun ρ a s Remark : Although µa ρa s for all s-sparse signals, this appears to be a highly pessimistic boun for most signals Numerical experiments have shown that for Gaussian ranom s-sparse signals, E [ρa] logs However, it remains to be seen if this non-uniformity can be exploite IV SYSTEM IDENTIFICATION / TOEPLITZ RECOVERY In this section, we aress problem of recovering a sparse signal after convolution We seek to fin a number of measurements of an s-sparse signal a using a ranom, compressive Toeplitz measurement matrix X that are sufficient to ensure that a can be recovere with high probability A Relate Work Consiering structure in the measurement matrix X was aresse by Tropp et al, as one of the early papers which consiere reconstruction of a signal from its convolution with a fixe Finite Impulse Response FIR filter having ranom taps [] Bajwa et al [] stuie Toeplitzstructure compresse sensing matrices with applications to sparse channel estimation They allowe the entries of the measurement matrix be rawn from a symmetric Bernoulli istribution Later they extene this to ranom matrices whose entries are boune or Gaussian-istribute It is shown that O s log n measurements are sufficient in orer to have exact signal recovery [4], [5], which appears to be the best result to ate This compares to O s log n/s measurements when X is unstructure, see eg [0] Rauhut consiere circulant an Toeplitz matrices whose entries are inepenent Bernoulli ± ranom variables Furthermore, they constraine the signs of the non-zero entries of the s-sparse signal to be rawn from a Bernoulli istribution as well Consiering this configuration, they propose that the require number of measurements for exact recovery scales linearly with sparsity up to a log-factor of the signal length [6] Xiang et al propose convolution with a white noise waveform followe by eterministic subsampling as a framework for CS Imposing a uniformly ranom sign pattern on the non-zero entries of the signal, they require O s log 5 n/δ measurements to have exact recovery with probability exceeing δ [7] B Recovery Conition: The Restricte Isometry Property Introuce by Canès an Tao [], the Restricte Isometry Property RIP is an isometry conition on measurement matrices use in CS Definition : X satisfies the RIP of orer s if there exists a δ s 0, such that δ s a Xa + δ s a 6 hols for all s-sparse signals a Given a matrix X R n an any set T of column inices, we enote by X T the T submatrix of X compose of these columns Equation 6 can be interprete as a requirement that all of the eigenvalues of the Grammian matrix X T X T lie in the range [ δ s, + δ s ] for all sets T with T s When the RIP is satisfie of orer s with δ s <, then correct recovery of an s-sparse signal can be obtaine by fining the minimizer of the following optimization problem [7]: min a a R n subject to y Aa In [7], it is also shown that robust recovery is possible in the presence of measurement noise Verifying the RIP is generally a ifficult task This property requires boune conition number for all submatrices built by selecting s T arbitrary columns Therefore, there are n s such submatrices to check, an computing the spectral norm of a matrix is not generally an easy task However, for certain ranom constructions of X, the RIP follows in a simple way from concentration of measure inequalities An important result in this vein concerns the n matrix X whose entries x i,j are inepenent realizations of Gaussian ranom variables, x i,j N 0, All Gaussian ranom variables mentione in this section will have this same mean an variance Note that in this realization of X, there are n fully inepenent Gaussian ranom entries
6 {x k } n+ k {a k } n k {z k } n+ kn+ + Fig FIR filter with impulse response a k y Xa + z We call such realizations simply Unstructure X comparing to Toeplitz X which have the form of 3 In previous work, a concentration of measure inequality has been use to show that Unstructure X satisfy the RIP with high probability The relevant concentration of measure inequality for Unstructure X is as follows see eg [0] For any 0, P { Xa a } a e c 0 7 with c Base on this this concentration inequality, the following result was proven in [0] Theorem : An Unstructure X with ii Gaussian entries satisfies the RIP of orer s with high probability for O s log n s An approach ientical to the one taken in [0] can be use to establish the RIP for Toeplitz X base on the concentration of measure inequalities given in Theorem In particular, some basic bouning gives us the following result for Toeplitz X For any 0,, P { Xa a } a e ρp c0 8 Comparing inequality 8 with 7 inicates that these concentration inequalities iffer only by the factor ρa Since for sparse signals, ρa is boune by s, we have the following result, whose proof is omitte for space Theorem 3: A Toeplitz X with ii Gaussian entries satisfies the RIP of orer s with high probability if O s log n s Note that this result is essentially ientical to the bouns given previously in [4], [5] However, given the extremely non-uniform istribution of ρa over the set of all s-sparse signals a, it may be that a moifie approach to the RIP conition will allow the recovery results to be tightene In particular, we note that if the sign pattern for the non-zero entries of the signal a is chosen to be ranom, ρ a logs This remains an area of current research V DETECTION WITH COMPRESSIVE SENSING A Problem Setup In another application, we consier solving a CS etection problem Consier an FIR filter with impulse response {a k } n k The response of this filter to a test signal x k is as escribe in Moreover, suppose the observations are corrupte by a ranom aitive measurement noise z Figure shows the schematic of this filter Now a etection problem can be consiere as follows: the impulse response a k is known, but we wish to etect whether the ynamics of the system change to b k, where the ifference c k b k a k is known The response {y k } n+ kn+ to a known input x k is monitore an the goal is to etect when the impulse response of the system changes using < n measurements Since the the nominal impulse response is known, the expecte response Xa can be subtracte off, an thus without loss of generality, we can consier a 0 The etection problem can be formulate as follows [8] Define two events E 0 an E as: E 0 y z E y Xc + z where z is a vector of ii Gaussian noise However, ue to the presence of ranom noise, the algorithm can result in false etection This leas to efining the probabilities P F A an P D as: P F A P {E chosen when E 0 } P D P {E chosen when E } where P F A enotes the false-alarm probability an P D enotes the etection probability A Receiver Operating Curve ROC is a plot of P D as a function of P F A A Neyman- Pearson NP etector maximizes P D for a given limit on failure probability, P F A α The NP test is base on the likelihoo ratio which in our case can be written as y Xc E E 0 γ 9 where the threshol γ is chosen to meet the constraint P F A α Consequently, we consier the compressive etector t as t : y Xc 30 By evaluating t an consiering the threshol γ, we are now able to istinguish between the occurrence of two events, E 0 an E To fix the failure limit, we set P F A α which leas to: P D α Q Q α Xc 3 σ where Qq π q e u u 3 Since P D α is a function of Xc, we postulate that the ultimate performance of the etector will epen on µc, ρc an consequently on P which is the sample covariance matrix for c Furthermore, this epenence woul not occur if the measurement utilize an Unstructure measurement matrix which, of course, woul not apply to the FIR filter measurement consiere here but is a useful comparison to the structure measurement impose by convolution B Experiments an ROCs In these experiments, we test ifferent signals c with ifferent ρ c values Several signal classes are esigne base on the location, sign an the value of their non-zero elements For a given fixe signal c, 000 Unstructure an Toeplitz X matrices are generate For each X, a curve of P D over P F A is compute using 3 Figure 3 shows the
7 Unstructure X Toeplitz X PD PD α Toeplitz X α PD ρ 5, µ ρ 50, µ 3334 ρ 3 75, µ ρ 4 00, µ ρ 5 50, µ ρ 6 00, µ α Fig 3 ROCs for 000 X Unstructure an Toeplitz matrices for a fixe signal c ρc 50 The soli black curve is the average of 000 curves Fig 5 Average ROCs with 000 Toeplitz X for 6 ifferent signals c The curves escen the same orer they appear in the legen box result for a special class of signals whose non-zero elements are put in a block position in front with same value an sign In this experiment signal length is n 50 an number of measurements 5 As can be seen, the ROCs associate with Toeplitz X are more scattere comparing to the ROCs when we apply Unstructure X This result is ue to the weaker concentration of Xc for Toeplitz X In orer to compare the ROCs of ifferent signals with ifferent ρ c values when projecte by ranom n Toeplitz matrices, we esign an experiment with 6 ifferent signals Figures 4 an 5 show how these signals get treate by Unstructure an Toeplitz X matrices PD Unstructure X ρ 5, µ 668 ρ 50, µ 3334 ρ 3 75, µ ρ 4 00, µ ρ 5 50, µ ρ 6 00, µ α Fig 4 Average ROCs with 000 Unstructure X for 6 ifferent signals c Note that all curves are overlapping As can be seen, the average ROCs o not change with respect to ifferent signals when we apply Unstructure X However, they change when applying Toeplitz X More importantly, it can be conclue that ρc has a irect influence on these curves Generally signals with higher spectral norm ρc have lower ROCs, although this relation is not firm In all of these experiments, we fix c an let the noise stanar eviation σ 03 VI CONCLUSION Concentration of measure inequalities for a compressive Toeplitz matrix are erive in this work We have consiere two important quantities ρa an µa associate with the eigenvalues of the covariance matrix P Two applications were aresse for these types of matrices A Toeplitz recovery application was consiere where we want to recover an s-sparse signal a from its compressive measurements via y Xa, where X is a ranom, compressive Toeplitz matrix We showe that for Toeplitz X, the require number of measurements shoul scale with sρ a for X to satisfy RIP of orer s for all s-sparse signals a In another application, a CS etection problem is consiere Experimental results show that signals with ifferent ρ a values have ifferent ROCs when measure by a Toeplitz matrix ACKNOWLEDGMENT The authors gratefully acknowlege Chris Rozell, Han Lun Yap, Alejanro Weinstein an Luis Tenorio for helpful conversations uring the evelopment of this work A Proof of Lemma APPENDIX Proof: It is easy to show that y is a zero mean vector, since E [y] E [Xa] E [X] a 0 The covariance matrix P associate with y can be calculate as follows: E [y y ] E [y y ] E [y y ] P E [yy E [y y ] E [y y ] E [y y ] ] E [y y ] E [y y ] E [y y ]
8 Now consier one of the elements eg E [y l y m ] Without loss of generality, assume m < l Then E [y l y m ] a E x n+l x n+l x l [ ] xn+m x n+m x m a a a a n l m 0 0 a a a n l m a i a i+l m R a l m, where the ones start at row l m + Clearly, with E [y l y m ] R a l m, the covariance matrix P has the form as in 4 which completes the proof Note that the covariance matrix is a symmetric Toeplitz matrix B Proof of Claim Proof: We start by taking the erivative of 4 with respect to t an then putting it equal to zero Using Jacobi s formula for an invertible matrix A, eta A etatr A, α α we have et I tp e +t a t et I tp e a +t et I tp tr I tp P a + e a +t et I tp et I tp e a +t tr I tp P a + Now setting the above equation equal to zero, we get tr I tp P a + Knowing that P is a symmetric matrix, we can ecompose it using the Singular Value Decomposition SVD P UΣU T, where U is an orthogonal matrix an Σ is a iagonal matrix Now for 0 < t < 05, we get tr I tp P tr U I tσ U T UΣU T tr tr UU T tuσu T UΣU T I tσ Σ From here, we get the following polynomial equation: λ i tλ i a + where λ i s are the eigenvalues of the P matrix which lie along the iagonal of Σ REFERENCES [] Data, ata everywhere, The Economist, 5 Feb 00 [] E Canès, Compressive sampling, in Proceeings of the International Congress of Mathematicians, vol, no 3 Citeseer, 006 [3] E Canès an T Tao, Near-optimal signal recovery from ranom projections: Universal encoing strategies? IEEE Transactions on Information Theory, vol 5, no, pp , 006 [4] E Canès, J Romberg, an T Tao, Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information, IEEE Transactions on information theory, vol 5, no, pp , 006 [5] D Donoho, Compresse sensing, IEEE Transactions on Information Theory, vol 5, no 4, pp , 006 [6] E Canès an M Wakin, An introuction to compressive sampling, IEEE Signal Processing Magazine, vol 5, no, pp 30, 008 [7] E Canès, The restricte isometry property an its implications for compresse sensing, Comptes renus-mathématique, vol 346, no 9-0, pp , 008 [8] M A Davenport, P T Boufounos, M B Wakin, an R G Baraniuk, Signal processing with compressive measurements, IEEE Journal of Selecte Topics in Signal Processing, vol 4, no, pp , 00 [9] M Leoux, The concentration of measure phenomenon Amer Mathematical Society, 00 [0] D Achlioptas, Database-frienly ranom projections: Johnsonlinenstrauss with binary coins, Journal of Computer an System Sciences, vol 66, no 4, pp , 003 [] J Tropp, M Wakin, M Duarte, D Baron, an R Baraniuk, Ranom filters for compressive sampling an reconstruction, in Proc Int Conf Acoustics, Speech, Signal Processing ICASSP, 006 [] W Bajwa, J Haupt, G Raz, S Wright, an R Nowak, Toeplitzstructure compresse sensing matrices, in IEEE/SP 4th Workshop on Statistical Signal Processing, 007 SSP 07, 007, pp [3] J Romberg, Compressive sensing by ranom convolution, SIAM Journal on Imaging Sciences, vol, no 4, pp 098 8, 009 [4] W Bajwa, J Haupt, G Raz, an R Nowak, Compresse channel sensing, in Proc 4n Annu Conf Information Sciences an Systems CISS08, 008, pp 5 0 [5] J Haupt, W Bajwa, G Raz, an R Nowak, Toeplitz compresse sensing matrices with applications to sparse channel estimation, IEEE Trans Inform Theory, 008 [6] H Rauhut, Circulant an Toeplitz matrices in compresse sensing, Proc SPARS, vol 9, 009 [7] Y Xiang, L Li, an F Li, Compressive sensing by white ranom convolution, 009 [Online] Available: [8] G Lugosi, Concentration-of-measure inequalities, Lecture notes, 004 [9] S Hwang, Cauchy s Interlace Theorem for Eigenvalues of Hermitian Matrices, The American Mathematical Monthly, vol, no, pp 57 59, 004 [0] R Baraniuk, M Davenport, R DeVore, an M Wakin, A simple proof of the restricte isometry property for ranom matrices, Constructive Approximation, vol 8, no 3, pp 53 63, 008 [] E Canès an T Tao, Decoing via linear programming, IEEE Trans Inform Theory, vol 5, no, pp , 005
Least-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationConcentration of Measure Inequalities for. Toeplitz Matrices with Applications
Concentration of Measure Inequalities for 1 Toeplitz Matrices with Applications Borhan M. Sanandaji, Tyrone L. Vincent, and Michael B. Wakin Abstract We derive Concentration of Measure (CoM) inequalities
More informationu!i = a T u = 0. Then S satisfies
Deterministic Conitions for Subspace Ientifiability from Incomplete Sampling Daniel L Pimentel-Alarcón, Nigel Boston, Robert D Nowak University of Wisconsin-Maison Abstract Consier an r-imensional subspace
More informationLecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012
CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration
More informationOn combinatorial approaches to compressed sensing
On combinatorial approaches to compresse sensing Abolreza Abolhosseini Moghaam an Hayer Raha Department of Electrical an Computer Engineering, Michigan State University, East Lansing, MI, U.S. Emails:{abolhos,raha}@msu.eu
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationTopic 7: Convergence of Random Variables
Topic 7: Convergence of Ranom Variables Course 003, 2016 Page 0 The Inference Problem So far, our starting point has been a given probability space (S, F, P). We now look at how to generate information
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More informationExact Topology Identification of Large-Scale Interconnected Dynamical Systems from Compressive Observations
Exact Topology Identification of arge-scale Interconnected Dynamical Systems from Compressive Observations Borhan M Sanandaji, Tyrone Vincent, and Michael B Wakin Abstract In this paper, we consider the
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationRobust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k
A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine
More informationCapacity Analysis of MIMO Systems with Unknown Channel State Information
Capacity Analysis of MIMO Systems with Unknown Channel State Information Jun Zheng an Bhaskar D. Rao Dept. of Electrical an Computer Engineering University of California at San Diego e-mail: juzheng@ucs.eu,
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationSurvey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013
Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing
More informationALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS
ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationLower Bounds for the Smoothed Number of Pareto optimal Solutions
Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.
More informationarxiv: v1 [cs.lg] 22 Mar 2014
CUR lgorithm with Incomplete Matrix Observation Rong Jin an Shenghuo Zhu Dept. of Computer Science an Engineering, Michigan State University, rongjin@msu.eu NEC Laboratories merica, Inc., zsh@nec-labs.com
More informationA New Minimum Description Length
A New Minimum Description Length Soosan Beheshti, Munther A. Dahleh Laboratory for Information an Decision Systems Massachusetts Institute of Technology soosan@mit.eu,ahleh@lis.mit.eu Abstract The minimum
More informationAPPLICATION of compressed sensing (CS) in radar signal
A Novel Joint Compressive Single Target Detection an Parameter Estimation in Raar without Signal Reconstruction Alireza Hariri, Massou Babaie-Zaeh Department of Electrical Engineering, Sharif University
More informationImage Denoising Using Spatial Adaptive Thresholding
International Journal of Engineering Technology, Management an Applie Sciences Image Denoising Using Spatial Aaptive Thresholing Raneesh Mishra M. Tech Stuent, Department of Electronics & Communication,
More informationAn Optimal Algorithm for Bandit and Zero-Order Convex Optimization with Two-Point Feedback
Journal of Machine Learning Research 8 07) - Submitte /6; Publishe 5/7 An Optimal Algorithm for Banit an Zero-Orer Convex Optimization with wo-point Feeback Oha Shamir Department of Computer Science an
More informationLower bounds on Locality Sensitive Hashing
Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,
More informationA PAC-Bayesian Approach to Spectrally-Normalized Margin Bounds for Neural Networks
A PAC-Bayesian Approach to Spectrally-Normalize Margin Bouns for Neural Networks Behnam Neyshabur, Srinah Bhojanapalli, Davi McAllester, Nathan Srebro Toyota Technological Institute at Chicago {bneyshabur,
More informationEntanglement is not very useful for estimating multiple phases
PHYSICAL REVIEW A 70, 032310 (2004) Entanglement is not very useful for estimating multiple phases Manuel A. Ballester* Department of Mathematics, University of Utrecht, Box 80010, 3508 TA Utrecht, The
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationTEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE
TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS Yannick DEVILLE Université Paul Sabatier Laboratoire Acoustique, Métrologie, Instrumentation Bât. 3RB2, 8 Route e Narbonne,
More informationSparse Reconstruction of Systems of Ordinary Differential Equations
Sparse Reconstruction of Systems of Orinary Differential Equations Manuel Mai a, Mark D. Shattuck b,c, Corey S. O Hern c,a,,e, a Department of Physics, Yale University, New Haven, Connecticut 06520, USA
More informationTOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH
English NUMERICAL MATHEMATICS Vol14, No1 Series A Journal of Chinese Universities Feb 2005 TOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH He Ming( Λ) Michael K Ng(Ξ ) Abstract We
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationMulti-View Clustering via Canonical Correlation Analysis
Technical Report TTI-TR-2008-5 Multi-View Clustering via Canonical Correlation Analysis Kamalika Chauhuri UC San Diego Sham M. Kakae Toyota Technological Institute at Chicago ABSTRACT Clustering ata in
More informationThe effect of dissipation on solutions of the complex KdV equation
Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,
More informationRobustness and Perturbations of Minimal Bases
Robustness an Perturbations of Minimal Bases Paul Van Dooren an Froilán M Dopico December 9, 2016 Abstract Polynomial minimal bases of rational vector subspaces are a classical concept that plays an important
More informationSpace-time Linear Dispersion Using Coordinate Interleaving
Space-time Linear Dispersion Using Coorinate Interleaving Jinsong Wu an Steven D Blostein Department of Electrical an Computer Engineering Queen s University, Kingston, Ontario, Canaa, K7L3N6 Email: wujs@ieeeorg
More informationNew Coherence and RIP Analysis for Weak. Orthogonal Matching Pursuit
New Coherence and RIP Analysis for Wea 1 Orthogonal Matching Pursuit Mingrui Yang, Member, IEEE, and Fran de Hoog arxiv:1405.3354v1 [cs.it] 14 May 2014 Abstract In this paper we define a new coherence
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationEfficiently Decodable Non-Adaptive Threshold Group Testing
Efficiently Decoable Non-Aaptive Threshol Group Testing arxiv:72759v3 [csit] 3 Jan 28 Thach V Bui, Minoru Kuribayashi, Mahi Cheraghchi, an Isao Echizen SOKENDAI The Grauate University for Avance Stuies),
More informationOn the Observability of Linear Systems from Random, Compressive Measurements
On the Observability of Linear Systems from Random, Compressive Measurements Michael B Wakin, Borhan M Sanandaji, and Tyrone L Vincent Abstract Recovering or estimating the initial state of a highdimensional
More informationANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS
ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics
More informationNew Bounds for Distributed Storage Systems with Secure Repair
New Bouns for Distribute Storage Systems with Secure Repair Ravi Tanon 1 an Soheil Mohajer 1 Discovery Analytics Center & Department of Computer Science, Virginia Tech, Blacksburg, VA Department of Electrical
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More informationNecessary and Sufficient Conditions for Sketched Subspace Clustering
Necessary an Sufficient Conitions for Sketche Subspace Clustering Daniel Pimentel-Alarcón, Laura Balzano 2, Robert Nowak University of Wisconsin-Maison, 2 University of Michigan-Ann Arbor Abstract This
More informationMultipath Matching Pursuit
Multipath Matching Pursuit Submitted to IEEE trans. on Information theory Authors: S. Kwon, J. Wang, and B. Shim Presenter: Hwanchol Jang Multipath is investigated rather than a single path for a greedy
More informationLeaving Randomness to Nature: d-dimensional Product Codes through the lens of Generalized-LDPC codes
Leaving Ranomness to Nature: -Dimensional Prouct Coes through the lens of Generalize-LDPC coes Tavor Baharav, Kannan Ramchanran Dept. of Electrical Engineering an Computer Sciences, U.C. Berkeley {tavorb,
More informationDelocalization of boundary states in disordered topological insulators
Journal of Physics A: Mathematical an Theoretical J. Phys. A: Math. Theor. 48 (05) FT0 (pp) oi:0.088/75-83/48//ft0 Fast Track Communication Delocalization of bounary states in isorere topological insulators
More informationAnalyzing Tensor Power Method Dynamics in Overcomplete Regime
Journal of Machine Learning Research 18 (2017) 1-40 Submitte 9/15; Revise 11/16; Publishe 4/17 Analyzing Tensor Power Metho Dynamics in Overcomplete Regime Animashree Ananumar Department of Electrical
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationFLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction
FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS ALINA BUCUR, CHANTAL DAVID, BROOKE FEIGON, MATILDE LALÍN 1 Introuction In this note, we stuy the fluctuations in the number
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationSpurious Significance of Treatment Effects in Overfitted Fixed Effect Models Albrecht Ritschl 1 LSE and CEPR. March 2009
Spurious Significance of reatment Effects in Overfitte Fixe Effect Moels Albrecht Ritschl LSE an CEPR March 2009 Introuction Evaluating subsample means across groups an time perios is common in panel stuies
More informationOn the number of isolated eigenvalues of a pair of particles in a quantum wire
On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationOptimal CDMA Signatures: A Finite-Step Approach
Optimal CDMA Signatures: A Finite-Step Approach Joel A. Tropp Inst. for Comp. Engr. an Sci. (ICES) 1 University Station C000 Austin, TX 7871 jtropp@ices.utexas.eu Inerjit. S. Dhillon Dept. of Comp. Sci.
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationWe G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies
We G15 5 Moel Reuction Approaches for Solution of Wave Equations for Multiple Frequencies M.Y. Zaslavsky (Schlumberger-Doll Research Center), R.F. Remis* (Delft University) & V.L. Druskin (Schlumberger-Doll
More informationApplications of the Wronskian to ordinary linear differential equations
Physics 116C Fall 2011 Applications of the Wronskian to orinary linear ifferential equations Consier a of n continuous functions y i (x) [i = 1,2,3,...,n], each of which is ifferentiable at least n times.
More informationTIME-DELAY ESTIMATION USING FARROW-BASED FRACTIONAL-DELAY FIR FILTERS: FILTER APPROXIMATION VS. ESTIMATION ERRORS
TIME-DEAY ESTIMATION USING FARROW-BASED FRACTIONA-DEAY FIR FITERS: FITER APPROXIMATION VS. ESTIMATION ERRORS Mattias Olsson, Håkan Johansson, an Per öwenborg Div. of Electronic Systems, Dept. of Electrical
More informationAcute sets in Euclidean spaces
Acute sets in Eucliean spaces Viktor Harangi April, 011 Abstract A finite set H in R is calle an acute set if any angle etermine by three points of H is acute. We examine the maximal carinality α() of
More informationQubit channels that achieve capacity with two states
Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March
More informationPermanent vs. Determinant
Permanent vs. Determinant Frank Ban Introuction A major problem in theoretical computer science is the Permanent vs. Determinant problem. It asks: given an n by n matrix of ineterminates A = (a i,j ) an
More informationRelative Entropy and Score Function: New Information Estimation Relationships through Arbitrary Additive Perturbation
Relative Entropy an Score Function: New Information Estimation Relationships through Arbitrary Aitive Perturbation Dongning Guo Department of Electrical Engineering & Computer Science Northwestern University
More informationSystems & Control Letters
Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen
More informationDesign of an Industrial Distributed Controller Near Spatial Domain Boundaries
Design of an Inustrial Distribute Controller Near Spatial Domain Bounaries Stevo Mijanovic, Greg E. Stewart, Guy A. Dumont, an Michael S. Davies ABSTRACT This paper proposes moifications to an inustrial
More informationMEASURES WITH ZEROS IN THE INVERSE OF THEIR MOMENT MATRIX
MEASURES WITH ZEROS IN THE INVERSE OF THEIR MOMENT MATRIX J. WILLIAM HELTON, JEAN B. LASSERRE, AND MIHAI PUTINAR Abstract. We investigate an iscuss when the inverse of a multivariate truncate moment matrix
More informationarxiv:hep-th/ v1 3 Feb 1993
NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø,
More informationarxiv: v4 [math.pr] 27 Jul 2016
The Asymptotic Distribution of the Determinant of a Ranom Correlation Matrix arxiv:309768v4 mathpr] 7 Jul 06 AM Hanea a, & GF Nane b a Centre of xcellence for Biosecurity Risk Analysis, University of Melbourne,
More informationTechnion - Computer Science Department - M.Sc. Thesis MSC Constrained Codes for Two-Dimensional Channels.
Technion - Computer Science Department - M.Sc. Thesis MSC-2006- - 2006 Constraine Coes for Two-Dimensional Channels Keren Censor Technion - Computer Science Department - M.Sc. Thesis MSC-2006- - 2006 Technion
More informationSignal Recovery from Permuted Observations
EE381V Course Project Signal Recovery from Permuted Observations 1 Problem Shanshan Wu (sw33323) May 8th, 2015 We start with the following problem: let s R n be an unknown n-dimensional real-valued signal,
More informationOn the Surprising Behavior of Distance Metrics in High Dimensional Space
On the Surprising Behavior of Distance Metrics in High Dimensional Space Charu C. Aggarwal, Alexaner Hinneburg 2, an Daniel A. Keim 2 IBM T. J. Watson Research Center Yortown Heights, NY 0598, USA. charu@watson.ibm.com
More informationSlide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13)
Slie10 Haykin Chapter 14: Neuroynamics (3r E. Chapter 13) CPSC 636-600 Instructor: Yoonsuck Choe Spring 2012 Neural Networks with Temporal Behavior Inclusion of feeback gives temporal characteristics to
More informationExperimental Robustness Study of a Second-Order Sliding Mode Controller
Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans
More informationOn Characterizing the Delay-Performance of Wireless Scheduling Algorithms
On Characterizing the Delay-Performance of Wireless Scheuling Algorithms Xiaojun Lin Center for Wireless Systems an Applications School of Electrical an Computer Engineering, Purue University West Lafayette,
More informationA Modification of the Jarque-Bera Test. for Normality
Int. J. Contemp. Math. Sciences, Vol. 8, 01, no. 17, 84-85 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1988/ijcms.01.9106 A Moification of the Jarque-Bera Test for Normality Moawa El-Fallah Ab El-Salam
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationA Review of Multiple Try MCMC algorithms for Signal Processing
A Review of Multiple Try MCMC algorithms for Signal Processing Luca Martino Image Processing Lab., Universitat e València (Spain) Universia Carlos III e Mari, Leganes (Spain) Abstract Many applications
More informationSYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS
SYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS GEORGE A HAGEDORN AND CAROLINE LASSER Abstract We investigate the iterate Kronecker prouct of a square matrix with itself an prove an invariance
More informationA simple model for the small-strain behaviour of soils
A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:
More informationParameter estimation: A new approach to weighting a priori information
Parameter estimation: A new approach to weighting a priori information J.L. Mea Department of Mathematics, Boise State University, Boise, ID 83725-555 E-mail: jmea@boisestate.eu Abstract. We propose a
More informationLecture 2: Correlated Topic Model
Probabilistic Moels for Unsupervise Learning Spring 203 Lecture 2: Correlate Topic Moel Inference for Correlate Topic Moel Yuan Yuan First of all, let us make some claims about the parameters an variables
More informationAll s Well That Ends Well: Supplementary Proofs
All s Well That Ens Well: Guarantee Resolution of Simultaneous Rigi Boy Impact 1:1 All s Well That Ens Well: Supplementary Proofs This ocument complements the paper All s Well That Ens Well: Guarantee
More informationSINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
Communications on Stochastic Analysis Vol. 2, No. 2 (28) 289-36 Serials Publications www.serialspublications.com SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationOptimal Signal Detection for False Track Discrimination
Optimal Signal Detection for False Track Discrimination Thomas Hanselmann Darko Mušicki Dept. of Electrical an Electronic Eng. Dept. of Electrical an Electronic Eng. The University of Melbourne The University
More informationMulti-View Clustering via Canonical Correlation Analysis
Keywors: multi-view learning, clustering, canonical correlation analysis Abstract Clustering ata in high-imensions is believe to be a har problem in general. A number of efficient clustering algorithms
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationExact Signal Recovery from Sparsely Corrupted Measurements through the Pursuit of Justice
Exact Signal Recovery from Sparsely Corrupted Measurements through the Pursuit of Justice Jason N. Laska, Mark A. Davenport, Richard G. Baraniuk Department of Electrical and Computer Engineering Rice University
More informationGeneralizing Kronecker Graphs in order to Model Searchable Networks
Generalizing Kronecker Graphs in orer to Moel Searchable Networks Elizabeth Boine, Babak Hassibi, Aam Wierman California Institute of Technology Pasaena, CA 925 Email: {eaboine, hassibi, aamw}@caltecheu
More informationPerfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs
Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs Ashish Goel Michael Kapralov Sanjeev Khanna Abstract We consier the well-stuie problem of fining a perfect matching in -regular bipartite
More informationMulti-View Clustering via Canonical Correlation Analysis
Kamalika Chauhuri ITA, UC San Diego, 9500 Gilman Drive, La Jolla, CA Sham M. Kakae Karen Livescu Karthik Sriharan Toyota Technological Institute at Chicago, 6045 S. Kenwoo Ave., Chicago, IL kamalika@soe.ucs.eu
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationTractability results for weighted Banach spaces of smooth functions
Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March
More information