SPECTRAL ANALYSIS OF NON-UNIFORMLY SAMPLED DATA: A NEW APPROACH VERSUS THE PERIODOGRAM
|
|
- Claud Nash
- 5 years ago
- Views:
Transcription
1 SPECTRAL ANALYSIS OF NON-UNIFORMLY SAMPLED DATA: A NEW APPROACH VERSUS THE PERIODOGRAM Hao He, Jian Li and Petre Stoica Dept. of Electrical and Computer Engineering, University of Florida, Gainesville, FL, USA Dept. of Information Technology, Uppsala University, Uppsala, Sweden ABSTRACT We begin by revisiting the plain least-squares periodogram () for real-valued data. Then we introduce a new method for spectral analysis of non-uniformly sampled data by iteratively weighting, and we name the new method realvalued iterative adaptive approach (). and are most suitable for data sequences with discrete spectra. For such type of data, we present a procedure to obtain a parametric spectral estimate, from the or nonparametric estimate, by means of the Bayesian information criterion (BIC). We also discuss a possible strategy for designing the sampling pattern of future measurements. Several numerical examples are provided to illustrate the performance of our proposed approaches. Index Terms Spectral analysis, non-uniformly sampled data, the periodogram, the least-squares method, the iterative adaptive approach, BIC. 1. INTRODUCTION Spectral analysis of non-uniformly sampled (NUS) data is widely used in many areas, like economics, biomedicine, and especially astronomy. can be readily applied to the NUS data case (see, e.g., [1]), but it suffers from both local and global leakage problems. Many new approaches have been proposed (see, e.g., [][3][]) to outperform the periodogram. A recent detailed discussion can be found in [5], where the well-known Capon, MUSIC and ESPRIT methods are extended from uniform to NUS data cases. The amplitude and phase estimation (APES) method, proposed in [] for uniformly sampled data, has significantly less leakage than the periodogram. We follow here the ideas in [7] to extend APES to the NUS data case, and the resulting new method is referred to as. Among all methods mentioned above, Capon (when used with the same covariance matrix as ) is essentially the same as the first iteration of, and other parametric methods all have special requirements (e.g., a guess of the number of sinusoidal components This work was supported in part by the National Science Foundation under Grant No. CCF-37 and ECCS-7977, and by the Swedish Research Council (VR). or further assumptions on noise) which are not as favorable as the non-parametric. We also combine with BIC (see, e.g., []) to provide parametric spectral estimates in the form of a number of estimated sinusoidal components that are deemed to fit the data well. The use of BIC bypasses the need for testing the significance of the periodogram peaks, for which, to the best of our knowledge, no satisfactory procedure exists for the NUS data case (see [9] and the references therein). Finally, we present a method for designing an optimal sampling pattern (see, e.g., [1][]) that imizes an objective function based on the spectral window. In doing so, we assume that a sufficient number of observations are already available, from which we can get a reasonably accurate spectral estimate. We make use of this spectral estimate to design the sampling times when future measurements should be performed.. AND THE SPECTRAL WINDOW Let {y(t n )} N denote the zero mean, real-valued sequence whose spectral analysis is our main goal. The classical Fourier transform-based periodogram (FP) P F (ω) = 1 N y(t n )e jωtn can be written as P F (ω) = ˆβ(ω), where ˆβ(ω) comes from the following least-squares (LS) data fitting problem: ˆβ(ω) =arg β(ω) y(t n ) β(ω)e jωtn. (1) Because {y(t n )} N R, it is better to use α cos(ωt n + φ) instead of βe jωtn to fit {y(t n )} N and the so obtained LS criterion is: a,b [y(t n ) a cos(ωt n ) b sin(ωt n )], () where a = α cos φ and b = α sin φ (we omit the dependence of α and φ on ω, for notational simplicity). With the
2 following additional notations denote the number of the grid points ( x is the largest integer x). Thus the frequency grid is ω k = kδω, k =1,...,K. y(t 1 ) [ ] a(ω) y =., θ(ω) =, A(ω) = [ c(ω) s(ω) ] The notations in Eq. (3) are accordingly modified to A k =, b(ω) A(ω k ), θ k = θ(ωk ). Now consider again the LS problem: y(t N ) c(ω) = cos(ωt 1 ). cos(ωt N ), s(ω) = sin(ωt 1 ). sin(ωt N ), (3) the LS fitting criterion in Eq. () can be written in the following vector form: θ y Aθ, () where denotes the Euclidean norm (and the dependence of A and θ on ω is omitted). The solution to the imization problem in Eq. () is well-known to be: ˆθ =(A T A) 1 A T y, (5) and the LS periodogram is accordingly given by P LS (ω) = 1 N (Aˆθ) T (Aˆθ) = 1 N ˆθ T (A T A)ˆθ. () To analyze the performance of, consider a possibly existing sinusoidal componet ȳ =[ȳ(t 1 ) ȳ(t N )] T, where ȳ(t n ) = ā cos( ωt n )+ b sin( ωt n ),n = 1,...,N. Its effect on Eq. (5) is to introduce an error term (A T A) 1 A T ȳ. For non-pathological sampling patterns, R N I and N cos((ω + ω)t n) N sin((ω + ω)t n). Then the squared Euclidean norm of the error term can be well approximated by 1 N P ( ω) e j(ω ω)tn, (7) where P ( ω) =ā (ω) + b (ω). The part of Eq. (7) that depends on the sampling pattern viz. W (ω) = e jωtn, () is called the spectral window. Then we follow [11] to define the maximum frequency interval [,ω max ] that can be dealt with unambiguously. The spectral window attains its maximum value of N at ω =: W () = N W (ω), ω. We evaluate W (ω) for ω> to find the smallest frequency value ω, at which the spectral window has a peak whose height is close to N. Then ω max = ω / can be claimed as the largest frequency of the data sequence in question. 3. AND BIC Let Δω denote the frequency grid size (small enough as to not affect the potential resolution), and let K = ω max /Δω y A k θ k, k =1,...,K. (9) With respect to the sinusoidal component with frequency ω k, all other sinusoidal components can be considered as noise and their contributions to the spectrum can be described by the covariance matrix defined as Q k = p=1,p k A p D p A T p, D p = a (ω p )+b (ω p ) [ 1 1 () Assug that Q k is available, and that it is invertible (a necessary condition is that (K 1) N, which is easily satisfied in general), it would make sense to consider the following weighted LS (WLS) criterion (which is more accurate than the corresponding LS problem, under quite general conditions), instead of Eq. (9), [y A k θ k ] T Q 1 k [y A kθ k ]. (11) The vector θ k that imizes Eq. (11) is given by: ˆθ k =(A T k Q 1 k A k) 1 (A T k Q 1 y). (1) Using the matrix inversion lemma, we obtain the following alternative expression for the WLS estimate: where Γ = ˆθ k =(A T k Γ 1 A k ) 1 (A T k Γ 1 y), (13) A p D p A T p = Q k + A k D k A T k. () p=1 Eq. (13) is computationally more appealing than Eq. (1), as Γ 1 in Eq. (13) needs to be computed only once for all values of k =1,...,K. Because Γ in Eq. (13) depends on the very quantities that we want to estimate, viz. the {θ k } K k=1, the only apparent solution is some form of iterative process, for example the algorithm outlined in Table 1. In most applications, the algorithm is expected to require no more than iterations (see, e.g., the numerical examples in Section 5). Next we present BIC as an alternative to test the significance of the doant peaks. Let {P (ω k )} K k=1 denote the values taken by either the or periodogram, and let { ω m, ă m, b m } M m=1 denote the frequency, amplitude and phase related parameters corresponding to the M largest peaks of {P (ω k )} K k=1, arranged in a decreasing order: P ( ω 1 ) P ( ω ) P ( ω M ), (17) k ].
3 Table 1. The Algorithm Initialization: Use the LS method (see Eq. (5)) to obtain initial estimates of {θ k }, denoted by {ˆθ k } Iteration: Let {ˆθ i k} denote the estimates of {θ k } at the i th iteration, and let ˆΓ i denote the estimate of Γ obtained from {ˆθ i k}. Fori =, 1,,..., compute: ˆθ i+1 k =[A T k (ˆΓ i ) 1 A k ] 1 [A T k (ˆΓ i ) 1 y], k =1,...,K, (15) until a given number of iterations is performed. Periodogram Calculation: Let {ˆθ I k} denote the estimates of {θ k } obtained by the above iterative process. The periodogram is computed as (see Eq. ()): P (ω k )= 1 N (ˆθ I k) T (A T k A k )(ˆθ I k), k=1,...,k. (1) for M 1. Under the idealizing assumptions that the data sequence consists of a finite number of sinusoidal components and of normal white noise, and that { ω m, ă m, b m } M m=1 are the maximum likelihood (ML) estimates, the BIC rule estimates M as follows: where ˆM =arg M=,1,,... BIC(M), (1) BIC(M) =ρm ln N [ M +N ln y(t n ) (ă m cos( ω m t n )+ b m sin( ω m t n )) m=1 ] (19) and where ρ = ρ ML =5(see, e.g., [1][]). It follows from the derivation of BIC (see, e.g., the cited works) that the large complexity penalty ρ ML is mainly due to the high accuracy of the ML estimates. By comparison, { ω m, ă m, b m } M m=1 estimated with and are expected to be less accurate than the ML estimates. Consequently, we suggest the use of BIC with smaller values of ρ in the case of and than the ρ ML =5, namely: ρ =1and ρ =. () Admittedly, these choices are somewhat ad-hoc, but the corresponding BIC rules are simple to use and they appear to provide accurate estimates of M (see Section 5 for details).. SAMPLING PATTERN DESIGN Suppose that the data samples {y(t n )} N have already been collected and a spectral estimate {P (ω k )} K k=1 has been obtained based on them; {τ g } G g=1 = S is the set of times when future measurements of y(t) could in principle be performed. Our goal is to choose the sampling times {x f } F f=1 from {τ g } G g=1 for future measurements in an optimal manner. The previous sections have evidenced that the performance of both and depends on the spectral window. With the {P (ω k )} K k=1 imposed as weights, the integrated spectral window is given at ω = ω p (see Eq. (7)) by G P (ω k ) e j(ωp ω k)t n + e j(ωp ω k)τ g. (1) k=1 g=1 A well-known fact is that the integrated spectral window is a constant in the uniform sampling case. This is also approximately true in the non-uniformly sampling case, which is the reason why we have to assume the prior information {P (ω k )} K k=1 available (otherwise the integrated spectral window would be insensitive to the choice of {x f } F f=1 ). For a given ω p, the larger the P (ω k ), the more emphasis will be put on imizing the sidelobe level of the spectral window corresponding to ω k. The function in Eq. (1) still depends on ω p. We eliate the frequency dependence by integrating Eq. (1) weighted by P (ω p ), with respect to ω p as well, and obtain the following imization problem: {μ g} G g=1 p=1 k=1 P (ω p ) G P (ω k ) e j(ωp ω k)t n + μ g e j(ωp ω k)τ g, s.t. μ g {, 1} and g=1 G μ g = F, () g=1 where {μ g } G g=1 {, 1} are Boolean variables. Because the Boolean constraint makes the problem hard, we relax it to: μ g 1, g =1,...,G, (3) with which Eq. () becomes a linearly constrained quadratic program that can be solved efficiently, e.g., by the MATLAB 1 function quadprog. Once the solution {μ g } G g=1 is found, we choose its F largest elements and set them to 1 (and the rest to ). Letting g 1 <g < <g F denote the indices of these elements, we use the following sampling times x f = τ gf, f =1,...,F () 1 MATLAB is a registered trademark of The MathWorks, Inc.
4 as an approximate solution to the imization problem of Eq. (). (See, e.g., [13] for a general discussion on this type of approximation to a Boolean optimization problem.) 5.1. Simulated Data 5. NUMERICAL EXAMPLES 3 Consider a data sequence consisting of M = 3 sinusoidal components with frequencies.1,. and.1 Hz, and amplitudes, and 5, respectively. The phases of the 3 sinusoids are independently and uniformly distributed over [, π] and the additive noise is white and normally distributed with mean and variance σ =.1. The sampling pattern follows a Poisson process with parameter λ =.1s 1, that is, the sampling times are exponentially distributed with mean μ = 1/λ = s. We generate N =samples and the sampling times are rounded off to ten decimals (see [11]). Figure 1 presents the spectral estimates averaged over independent realizations while Figure shows the overlapped estimates from the first 15 Monte-Carlo trials. nearly misses the smallest sinusoid while successfully resolves all three sinusoids. Note that suffers from much less variability than from one trial to another. Next, we use BIC to select the significant peaks of the spectra obtained by and. Let M BIC denote the number of sinusoids picked up by BIC. Then the probabilities of correct detection, false alarm and miss are defined as P D = Prob.(M BIC = 3), P FA = Prob.(M BIC > 3) and P M = Prob.(M BIC < 3). Figure 3 shows the scatter plots of the estimates for the doant sinusoids obtained via using +BIC and +BIC. The amplitude estimates of the two closely-spaced sinusoids are biased and the smallest sinusoid is frequently missed. For, P D is only.3, P FA is., and P M is as high as.97. Compared to, shows much better stability and accuracy, as illustrated in Figure 3. For, P D =.9, P FA =., and P M =.. We also note that usually after 15 iterations s performance does not improve visibly. So in the above and all subsequent examples, we terate after 15 iterations. 5.. Sampling Pattern Design Suppose we have already collected data at {t n } N=3, which are the sampling times of the measured astronomical data HD195 (whose original analysis led to the discovery of an extrasolar planet named ET-1, see []). The unit of the sampling times is changed from day to second for simplicity. Figure shows the sampling pattern and Figure shows the corresponding spectral window. The smallest frequency f > at which the spectral window has a peak close to N is approximately 1 Hz. According to the discussion in Section, f max = f /.5 Hz. The set of sampling times of possible future measurements {τ g } G= g=1 is randomly selected Fig. 1. Average spectral estimates from Monte-Carlo trials. The solid line is the estimated spectrum and the circles represent the true frequencies and amplitudes of the 3 sinusoids. and Fig.. Spectral estimates from 15 Monte-Carlo trials. The solid lines are the estimated spectra and the circles represent the true frequencies and amplitudes of the 3 sinusoids. and.
5 by and the amplitude estimates are also quite accurate Fig. 3. Scatter plot of the doant sinusoidal components selected by BIC in Monte-Carlo trials. Dots symbolize the estimates while circles represent the true frequencies and amplitudes of the 3 sinusoids. +BIC, and +BIC. from the interval [, ] s, which covers approximately the same period of time as {t n } N=3. Consider the same data sequence as in the previous subsection, except that the amplitudes of the two closely-spaced sinusoids (at. and.1 Hz) are now set to and 5, respectively. Figure 5 shows the scatter plots of the estimated frequencies and amplitudes using this data sequence sampled at {t n } N=3. In most trials +BIC catches only the strongest sinusoid. +BIC, on the other hand, has a much better detection rate, although the amplitude estimates are somewhat biased. The corresponding probabilities P D, P FA and P M are.1,. and.99 for and.73,.1 and. for, respectively. Next we use the optimal sampling design technique outlined in Section to choose {x f } F f=1 =3 from {τ g } G= g=1. In Eq. () the {P (ω k )} K k=1 are equal to the spectrum estimates obtained from the given 3 samples and N is chosen as because we want to achieve better estimation from only the designed samples. Figure shows the designed sampling pattern and the corresponding spectral window (although the inferred f max is much larger than.5 Hz, we still choose f max =.5Hz for a fair comparison). With the data sequence sampled at {x f } F f=1 =3, we obtain the scatter plots of the new estimates in Figure 7. The new probabilities P D, P FA and P M are.5,. and.95 for ; and.9,. and. for. All 3 sinusoidal components are now clearly resolved. CONCLUDING REMARKS We have presented a new spectral estimation method named for real-valued non-uniformly sampled data. is non-parametric, has better resolution than, and is most suitable for sinusoidal data, in which case we combine with BIC to obtain a parametric spectral estimate. We have also shown a sampling pattern design strategy, which makes use of the spectral window and the already obtained spectral estimates for designing the future sampling times in an optimal manner. Several numerical examples have been provided to show the properties of, especially compared with the benchmark periodogram. 7. REFERENCES [1] J. D. Scargle, Studies in astronomical time series analysis II: Statistical aspects of spectral analysis of unevenly spaced data, The Astrophysical Journal, vol. 3, pp , December 19. [] R. H. Jones, Fitting a continuous time autoregression to discrete data, Appl. Time Series Anal. II, pp. 51, 191. [3] R. J. Martin, Autoregression and irregular sampling: Spectral estimation, Signal Processing, vol. 77, pp , [] A. Rivoira and G. A. Fleury, A consistent nonparametric spectral estimator for randomly sampled signals, IEEE Transactions on Signal Processing, vol. 5, pp ,. [5] P. Stoica and N. Sandgren, Spectral analysis of irregularlysampled data: Paralleling the regularly-sampled data approaches, Digital Signal Processng, vol. 1, no., pp , November. [] J. Li and P. Stoica, An adaptive filtering approach to spectral estimation and SAR imaging, IEEE Transactions on Signal Processing, vol., no., pp. 9, June 199. [7] T. Yardibi, J. Li, P. Stoica, M. Xue, and A. B. Baggeroer, Source localization and sensing: A nonparametric iterative adaptive approach based on weighted least squares, submitted to IEEE Transactions on Aerospace and Electronic Systems, 7. [] P. Stoica and Y. Selén, Model-order selection: a review of information criterion rules, IEEE Signal Processing Magazine, vol. 1, no., pp. 3 7, July. [9] F. A. M. Frescura, C. A. Engelbrecht, and B. S. Frank, Significance tests for periodogram peaks, E-print arxiv:7.5, June 7. [] E. S. Saunders, T. Naylor, and A. Allan, Optimal placement of a limited number of observations for period searches, Astronomy & Astrophysics, vol. 55, pp , May. [11] L. Eyer and P. Bartholdi, Variable stars: Which Nyquist frequency?, Astronomy and Astrophysics, Supplement Series, vol. 135, pp. 1 3, February [1] P. Stoica and R. L. Moses, Spectral Analysis of Signals, Prentice-Hall, Upper Saddle River, NJ, 5. [13] S. Boyd and L. Vandenberghe, Convex Optimization, Cambridge University Press, Cambridge, United Kingdom,. [] J. Ge, J. van Eyken, S. Mahadevan, C. DeWitt, S. R. Kane, R. Cohen, A. V. Heuvel, S. W. Fleg, and P. Guo, The first extrasolar planet discovered with a new-generation highthroughput Doppler instrument, The Astrophysical Journal, vol., pp. 3 95, September.
6 Sampling Pattern Sampling Pattern Sampling Times (sec) Sampling Times (sec) Spectral Window (db) Spectral Window (db) Fig.. Sampling pattern and spectral window for the given sampling times in the sampling design example. The sampling pattern used for the Monte-Carlo trials in Figures 5. The distance between two consecutive bars represents the sampling interval. The corresponding spectral window. Fig.. Sampling pattern and spectral window for the designed sampling times in the sampling design example. The sampling pattern used for the Monte-Carlo trials in Figures 7. The distance between two consecutive bars represents the sampling interval. The corresponding spectral window Fig. 5. Scatter plot of the doant sinusoidal components selected by BIC in Monte-Carlo trials, using the given 3 sampling times. Dots symbolize the estimates while circles represent the true frequencies and amplitudes of the 3 sinusoids. +BIC, and +BIC. Fig. 7. Scatter plot of the doant sinusoidal components selected by BIC in Monte-Carlo trials, using the designed 3 sampling times. Dots symbolize the estimates while circles represent the true frequencies and amplitudes of the 3 sinusoids. +BIC, and +BIC.
Spectral Analysis of Non-Uniformly Sampled Data: A New Approach Versus the Periodogram
IEEE TRANSACTIONS ON SIGNAL PROCESSING 1 Spectral Analysis of Non-Uniformly Sampled Data: A New Approach Versus the Periodogram Petre Stoica, Fellow, IEEE, Jian Li, Fellow, IEEE, Hao He Abstract We begin
More informationIEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 57, NO. 3, MARCH X/$ IEEE
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL 57, NO 3, MARCH 2009 843 Spectral Analysis of Nonuniformly Sampled Data: A New Approach Versus the Periodogram Petre Stoica, Fellow, IEEE, Jian Li, Fellow, IEEE,
More informationOn Moving Average Parameter Estimation
On Moving Average Parameter Estimation Niclas Sandgren and Petre Stoica Contact information: niclas.sandgren@it.uu.se, tel: +46 8 473392 Abstract Estimation of the autoregressive moving average (ARMA)
More informationsine wave fit algorithm
TECHNICAL REPORT IR-S3-SB-9 1 Properties of the IEEE-STD-57 four parameter sine wave fit algorithm Peter Händel, Senior Member, IEEE Abstract The IEEE Standard 57 (IEEE-STD-57) provides algorithms for
More informationUnstable Oscillations!
Unstable Oscillations X( t ) = [ A 0 + A( t ) ] sin( ω t + Φ 0 + Φ( t ) ) Amplitude modulation: A( t ) Phase modulation: Φ( t ) S(ω) S(ω) Special case: C(ω) Unstable oscillation has a broader periodogram
More informationSPICE : a sparse covariance-based estimation method for array processing
SPICE : a sparse covariance-based estimation method for array processing Petre Stoica, Fellow, IEEE, Prabhu Babu, Jian Li, Fellow, IEEE 1 Abstract This paper presents a novel sparse iterative covariance-based
More informationPeriodogram of a sinusoid + spike Single high value is sum of cosine curves all in phase at time t 0 :
Periodogram of a sinusoid + spike Single high value is sum of cosine curves all in phase at time t 0 : X(t) = µ + Asin(ω 0 t)+ Δ δ ( t t 0 ) ±σ N =100 Δ =100 χ ( ω ) Raises the amplitude uniformly at all
More informationPARAMETER ESTIMATION AND ORDER SELECTION FOR LINEAR REGRESSION PROBLEMS. Yngve Selén and Erik G. Larsson
PARAMETER ESTIMATION AND ORDER SELECTION FOR LINEAR REGRESSION PROBLEMS Yngve Selén and Eri G Larsson Dept of Information Technology Uppsala University, PO Box 337 SE-71 Uppsala, Sweden email: yngveselen@ituuse
More informationPeriodogram of a sinusoid + spike Single high value is sum of cosine curves all in phase at time t 0 :
Periodogram of a sinusoid + spike Single high value is sum of cosine curves all in phase at time t 0 : ( ) ±σ X(t) = µ + Asin(ω 0 t)+ Δ δ t t 0 N =100 Δ =100 χ ( ω ) Raises the amplitude uniformly at all
More informationSource Localization and Sensing: A Nonparametric Iterative Adaptive Approach Based on Weighted Least Squares
IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS Source Localization and Sensing: A i Nonparametric Iterative Adaptive Approach Based on Weighted Least Squares Tarik Yardibi, Student Member, IEEE,
More informationGRIDLESS COMPRESSIVE-SENSING METHODS FOR FREQUENCY ESTIMATION: POINTS OF TANGENCY AND LINKS TO BASICS
GRIDLESS COMPRESSIVE-SENSING METHODS FOR FREQUENCY ESTIMATION: POINTS OF TANGENCY AND LINKS TO BASICS Petre Stoica, Gongguo Tang, Zai Yang, Dave Zachariah Dept. Information Technology Uppsala University
More informationREGULARIZED ESTIMATION OF LINE SPECTRA FROM IRREGULARLY SAMPLED ASTROPHYSICAL DATA. Sébastien Bourguignon, Hervé Carfantan, Loïc Jahan
REGULARIZED ESTIMATION OF LINE SPECTRA FROM IRREGULARLY SAMPLED ASTROPHYSICAL DATA Sébastien Bourguignon, Hervé Carfantan, Loïc Jahan Laboratoire d Astrophysique de Toulouse - Tarbes, UMR 557 CNRS/UPS
More informationSPICE and LIKES : Two hyperparameter-free methods for sparse-parameter estimation
SPICE and LIKES : Two hyperparameter-free methods for sparse-parameter estimation Petre Stoica, Prabhu Babu Department of Information Technology, Uppsala University, SE-751 5 Uppsala, Sweden Abstract SPICE
More informationarxiv:astro-ph/ v1 14 Sep 2005
For publication in Bayesian Inference and Maximum Entropy Methods, San Jose 25, K. H. Knuth, A. E. Abbas, R. D. Morris, J. P. Castle (eds.), AIP Conference Proceeding A Bayesian Analysis of Extrasolar
More informationProblem Sheet 1 Examples of Random Processes
RANDOM'PROCESSES'AND'TIME'SERIES'ANALYSIS.'PART'II:'RANDOM'PROCESSES' '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''Problem'Sheets' Problem Sheet 1 Examples of Random Processes 1. Give
More informationHigh-resolution Parametric Subspace Methods
High-resolution Parametric Subspace Methods The first parametric subspace-based method was the Pisarenko method,, which was further modified, leading to the MUltiple SIgnal Classification (MUSIC) method.
More informationRobust Implementation of the MUSIC algorithm Zhang, Johan Xi; Christensen, Mads Græsbøll; Dahl, Joachim; Jensen, Søren Holdt; Moonen, Marc
Aalborg Universitet Robust Implementation of the MUSIC algorithm Zhang, Johan Xi; Christensen, Mads Græsbøll; Dahl, Joachim; Jensen, Søren Holdt; Moonen, Marc Published in: I E E E International Conference
More informationAN ITERATIVE ADAPTIVE APPROACH FOR BLOOD VELOCITY ESTIMATION USING ULTRASOUND
18th European Signal Processing Conference (EUSIPCO-1) Aalborg Denmark August 3-7 1 AN ITERATIVE ADAPTIVE APPROACH FOR BLOOD VELOCITY ESTIMATION USING ULTRASOUND Erik Gudmundson Andreas Jakobsson Jørgen
More informationFrequency estimation by DFT interpolation: A comparison of methods
Frequency estimation by DFT interpolation: A comparison of methods Bernd Bischl, Uwe Ligges, Claus Weihs March 5, 009 Abstract This article comments on a frequency estimator which was proposed by [6] and
More informationJ. Liang School of Automation & Information Engineering Xi an University of Technology, China
Progress In Electromagnetics Research C, Vol. 18, 245 255, 211 A NOVEL DIAGONAL LOADING METHOD FOR ROBUST ADAPTIVE BEAMFORMING W. Wang and R. Wu Tianjin Key Lab for Advanced Signal Processing Civil Aviation
More informationEstimating Periodic Signals
Department of Mathematics & Statistics Indian Institute of Technology Kanpur Most of this talk has been taken from the book Statistical Signal Processing, by D. Kundu and S. Nandi. August 26, 2012 Outline
More informationEEE508 GÜÇ SİSTEMLERİNDE SİNYAL İŞLEME
EEE508 GÜÇ SİSTEMLERİNDE SİNYAL İŞLEME Signal Processing for Power System Applications Frequency Domain Analysis Techniques Parametric Methods for Line Spectra (Week-5-6) Gazi Üniversitesi, Elektrik ve
More informationLIKELIHOOD-BASED ESTIMATION OF PERIODICITIES IN SYMBOLIC SEQUENCES. Dept. of Mathematical Statistics, Lund University, Sweden
LIKELIHOOD-BASED ESTIMATION OF PERIODICITIES IN SYMBOLIC SEQUENCES Stefan Ingi Adalbjörnsson, Johan Swärd, and Andreas Jaobsson Dept. of Mathematical Statistics, Lund University, Sweden ABSTRACT In this
More informationTransient Disturbance Detection for Power Systems with a General Likelihood Ratio Test
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Transient Disturbance Detection for Power Systems with a General Likelihood Ratio Test Song, X.; Sahinoglu, Z.; Guo, J. TR23-24 May 23 Abstract
More informationDesign Criteria for the Quadratically Interpolated FFT Method (I): Bias due to Interpolation
CENTER FOR COMPUTER RESEARCH IN MUSIC AND ACOUSTICS DEPARTMENT OF MUSIC, STANFORD UNIVERSITY REPORT NO. STAN-M-4 Design Criteria for the Quadratically Interpolated FFT Method (I): Bias due to Interpolation
More informationA ROBUST BEAMFORMER BASED ON WEIGHTED SPARSE CONSTRAINT
Progress In Electromagnetics Research Letters, Vol. 16, 53 60, 2010 A ROBUST BEAMFORMER BASED ON WEIGHTED SPARSE CONSTRAINT Y. P. Liu and Q. Wan School of Electronic Engineering University of Electronic
More informationDigital Signal Processing
Digital Signal Processing Introduction Moslem Amiri, Václav Přenosil Embedded Systems Laboratory Faculty of Informatics, Masaryk University Brno, Czech Republic amiri@mail.muni.cz prenosil@fi.muni.cz February
More informationAn Efficient and Fast Algorithm for Estimating the Parameters of Two-Dimensional Sinusoidal Signals
isid/ms/8/ November 6, 8 http://www.isid.ac.in/ statmath/eprints An Efficient and Fast Algorithm for Estimating the Parameters of Two-Dimensional Sinusoidal Signals Swagata Nandi Anurag Prasad Debasis
More informationStatistical and Adaptive Signal Processing
r Statistical and Adaptive Signal Processing Spectral Estimation, Signal Modeling, Adaptive Filtering and Array Processing Dimitris G. Manolakis Massachusetts Institute of Technology Lincoln Laboratory
More informationParameter estimation in epoch folding analysis
ASTRONOMY & ASTROPHYSICS MAY II 1996, PAGE 197 SUPPLEMENT SERIES Astron. Astrophys. Suppl. Ser. 117, 197-21 (1996) Parameter estimation in epoch folding analysis S. Larsson Stockholm Observatory, S-13336
More informationDirection of Arrival Estimation: Subspace Methods. Bhaskar D Rao University of California, San Diego
Direction of Arrival Estimation: Subspace Methods Bhaskar D Rao University of California, San Diego Email: brao@ucsdedu Reference Books and Papers 1 Optimum Array Processing, H L Van Trees 2 Stoica, P,
More information2.161 Signal Processing: Continuous and Discrete
MIT OpenCourseWare http://ocw.mit.edu.6 Signal Processing: Continuous and Discrete Fall 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. M MASSACHUSETTS
More information8.2 Harmonic Regression and the Periodogram
Chapter 8 Spectral Methods 8.1 Introduction Spectral methods are based on thining of a time series as a superposition of sinusoidal fluctuations of various frequencies the analogue for a random process
More information! # % & () +,.&/ 01),. &, / &
! # % & () +,.&/ ),. &, / & 2 A NEW METHOD FOR THE DESIGN OF ENERGY TRANSFER FILTERS Xiaofeng Wu, Z Q Lang and S. A. Billings Department of Automatic Control and Systems Engineering The University of Sheffield
More informationBispectral resolution and leakage effect of the indirect bispectrum estimate for different types of 2D window functions
Bispectral resolution and leakage effect of the indirect bispectrum estimate for different types of D window functions Teofil-Cristian OROIAN, Constantin-Iulian VIZITIU, Florin ŞERBAN Communications and
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 25 Feb 2003
arxiv:cond-mat/0302507v1 [cond-mat.stat-mech] 25 Feb 2003 SIGNIFICANCE OF LOG-PERIODIC SIGNATURES IN CUMULATIVE NOISE HANS-CHRISTIAN GRAF V. BOTHMER Abstract. Using methods introduced by Scargle we derive
More informationLasso Maximum Likelihood Estimation of Parametric Models with Singular Information Matrices
Article Lasso Maximum Likelihood Estimation of Parametric Models with Singular Information Matrices Fei Jin 1,2 and Lung-fei Lee 3, * 1 School of Economics, Shanghai University of Finance and Economics,
More informationDetection theory 101 ELEC-E5410 Signal Processing for Communications
Detection theory 101 ELEC-E5410 Signal Processing for Communications Binary hypothesis testing Null hypothesis H 0 : e.g. noise only Alternative hypothesis H 1 : signal + noise p(x;h 0 ) γ p(x;h 1 ) Trade-off
More informationOn the Frequency-Domain Properties of Savitzky-Golay Filters
On the Frequency-Domain Properties of Savitzky-Golay Filters Ronald W Schafer HP Laboratories HPL-2-9 Keyword(s): Savitzky-Golay filter, least-squares polynomial approximation, smoothing Abstract: This
More informationPart III Spectrum Estimation
ECE79-4 Part III Part III Spectrum Estimation 3. Parametric Methods for Spectral Estimation Electrical & Computer Engineering North Carolina State University Acnowledgment: ECE79-4 slides were adapted
More informationTHIS paper studies the input design problem in system identification.
1534 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 50, NO. 10, OCTOBER 2005 Input Design Via LMIs Admitting Frequency-Wise Model Specifications in Confidence Regions Henrik Jansson Håkan Hjalmarsson, Member,
More informationBeamspace Adaptive Channel Compensation for Sensor Arrays with Faulty Elements
1 2005 Asilomar Conference Beamspace Adaptive Channel Compensation for Sensor Arrays with Faulty Elements Oguz R. Kazanci and Jeffrey L. Krolik Duke University Department of Electrical and Computer Engineering
More informationDESIGN OF PIECEWISE LINEAR POLYPHASE SEQUENCES WITH GOOD CORRELATION PROPERTIES
DESIGN OF PIECEWISE LINEAR POLYPHASE SEQUENCES WITH GOOD CORRELATION PROPERTIES Mojtaba Soltanalian +, Petre Stoica +, Mohammad Mahdi Naghsh and Antonio De Maio + Dept of Information Technology, Uppsala
More informationThe Spectral Density Estimation of Stationary Time Series with Missing Data
The Spectral Density Estimation of Stationary Time Series with Missing Data Jian Huang and Finbarr O Sullivan Department of Statistics University College Cork Ireland Abstract The spectral estimation of
More informationMultirate Digital Signal Processing
Multirate Digital Signal Processing Basic Sampling Rate Alteration Devices Up-sampler - Used to increase the sampling rate by an integer factor Down-sampler - Used to decrease the sampling rate by an integer
More informationAn Efficient and Fast Algorithm for Estimating the Parameters of Sinusoidal Signals
An Efficient and Fast Algorithm for Estimating the Parameters of Sinusoidal Signals Swagata Nandi 1 Debasis Kundu Abstract A computationally efficient algorithm is proposed for estimating the parameters
More informationCONTROL SYSTEMS, ROBOTICS, AND AUTOMATION - Vol. V - Prediction Error Methods - Torsten Söderström
PREDICTIO ERROR METHODS Torsten Söderström Department of Systems and Control, Information Technology, Uppsala University, Uppsala, Sweden Keywords: prediction error method, optimal prediction, identifiability,
More informationII. Nonparametric Spectrum Estimation for Stationary Random Signals - Non-parametric Methods -
II. onparametric Spectrum Estimation for Stationary Random Signals - on-parametric Methods - - [p. 3] Periodogram - [p. 12] Periodogram properties - [p. 23] Modified periodogram - [p. 25] Bartlett s method
More informationExpressions for the covariance matrix of covariance data
Expressions for the covariance matrix of covariance data Torsten Söderström Division of Systems and Control, Department of Information Technology, Uppsala University, P O Box 337, SE-7505 Uppsala, Sweden
More informationPractical Spectral Estimation
Digital Signal Processing/F.G. Meyer Lecture 4 Copyright 2015 François G. Meyer. All Rights Reserved. Practical Spectral Estimation 1 Introduction The goal of spectral estimation is to estimate how the
More informationTwo results in statistical decision theory for detecting signals with unknown distributions and priors in white Gaussian noise.
Two results in statistical decision theory for detecting signals with unknown distributions and priors in white Gaussian noise. Dominique Pastor GET - ENST Bretagne, CNRS UMR 2872 TAMCIC, Technopôle de
More informationVariable and Periodic Signals in Astronomy
Lecture 14: Variability and Periodicity Outline 1 Variable and Periodic Signals in Astronomy 2 Lomb-Scarle diagrams 3 Phase dispersion minimisation 4 Kolmogorov-Smirnov tests 5 Fourier Analysis Christoph
More informationCopyright license. Exchanging Information with the Stars. The goal. Some challenges
Copyright license Exchanging Information with the Stars David G Messerschmitt Department of Electrical Engineering and Computer Sciences University of California at Berkeley messer@eecs.berkeley.edu Talk
More informationExpectation propagation for signal detection in flat-fading channels
Expectation propagation for signal detection in flat-fading channels Yuan Qi MIT Media Lab Cambridge, MA, 02139 USA yuanqi@media.mit.edu Thomas Minka CMU Statistics Department Pittsburgh, PA 15213 USA
More informationOverview of Beamforming
Overview of Beamforming Arye Nehorai Preston M. Green Department of Electrical and Systems Engineering Washington University in St. Louis March 14, 2012 CSSIP Lab 1 Outline Introduction Spatial and temporal
More informationQuantum algorithms for computing short discrete logarithms and factoring RSA integers
Quantum algorithms for computing short discrete logarithms and factoring RSA integers Martin Ekerå, Johan Håstad February, 07 Abstract In this paper we generalize the quantum algorithm for computing short
More informationA6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2013
A6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2013 Lecture 26 Localization/Matched Filtering (continued) Prewhitening Lectures next week: Reading Bases, principal
More informationVariable and Periodic Signals in Astronomy
Lecture 14: Variability and Periodicity Outline 1 Variable and Periodic Signals in Astronomy 2 Lomb-Scarle diagrams 3 Phase dispersion minimisation 4 Kolmogorov-Smirnov tests 5 Fourier Analysis Christoph
More informationMean Likelihood Frequency Estimation
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 48, NO. 7, JULY 2000 1937 Mean Likelihood Frequency Estimation Steven Kay, Fellow, IEEE, and Supratim Saha Abstract Estimation of signals with nonlinear as
More informationEL1820 Modeling of Dynamical Systems
EL1820 Modeling of Dynamical Systems Lecture 10 - System identification as a model building tool Experiment design Examination and prefiltering of data Model structure selection Model validation Lecture
More informationThis is the published version of a paper presented at IEEE International Instrumentation and Measurement Technology Conference (I2MTC), 2013.
http://www.diva-portal.org This is the published version of a paper presented at IEEE International Instrumentation and Measurement Technology Conference (IMTC), 3. Citation for the original published
More informationRobust Capon Beamforming
Robust Capon Beamforming Yi Jiang Petre Stoica Zhisong Wang Jian Li University of Florida Uppsala University University of Florida University of Florida March 11, 2003 ASAP Workshop 2003 1 Outline Standard
More informationTransmit Codes and Receive Filters for Pulse Compression Radar Systems
Transmit Codes and Receive Filters for Pulse Compression Radar Systems Petre Stoica, Jian Li, and ing Xue Abstract Pulse compression radar systems make use of transmit code sequences and receive filters
More informationExploring Granger Causality for Time series via Wald Test on Estimated Models with Guaranteed Stability
Exploring Granger Causality for Time series via Wald Test on Estimated Models with Guaranteed Stability Nuntanut Raksasri Jitkomut Songsiri Department of Electrical Engineering, Faculty of Engineering,
More informationD.S.G. POLLOCK: BRIEF NOTES
BIVARIATE SPECTRAL ANALYSIS Let x(t) and y(t) be two stationary stochastic processes with E{x(t)} = E{y(t)} =. These processes have the following spectral representations: (1) x(t) = y(t) = {cos(ωt)da
More informationTowards Global Design of Orthogonal Filter Banks and Wavelets
Towards Global Design of Orthogonal Filter Banks and Wavelets Jie Yan and Wu-Sheng Lu Department of Electrical and Computer Engineering University of Victoria Victoria, BC, Canada V8W 3P6 jyan@ece.uvic.ca,
More informationEIE6207: Estimation Theory
EIE6207: Estimation Theory Man-Wai MAK Dept. of Electronic and Information Engineering, The Hong Kong Polytechnic University enmwmak@polyu.edu.hk http://www.eie.polyu.edu.hk/ mwmak References: Steven M.
More informationarxiv: v1 [math.st] 1 Dec 2014
HOW TO MONITOR AND MITIGATE STAIR-CASING IN L TREND FILTERING Cristian R. Rojas and Bo Wahlberg Department of Automatic Control and ACCESS Linnaeus Centre School of Electrical Engineering, KTH Royal Institute
More informationFAST RECOVERY ALGORITHMS FOR TIME ENCODED BANDLIMITED SIGNALS. Ernő K. Simonyi
FAST RECOVERY ALGORITHS FOR TIE ENCODED BANDLIITED SIGNALS Aurel A Lazar Dept of Electrical Engineering Columbia University, New York, NY 27, USA e-mail: aurel@eecolumbiaedu Ernő K Simonyi National Council
More informationFalse Alarm Probability based on bootstrap and extreme-value methods for periodogram peaks
False Alarm Probability based on bootstrap and extreme-value methods for periodogram peaks M. Süveges ISDC Data Centre for Astrophysics, University of Geneva, Switzerland Abstract One of the most pertinent
More informationFAST AND ACCURATE DIRECTION-OF-ARRIVAL ESTIMATION FOR A SINGLE SOURCE
Progress In Electromagnetics Research C, Vol. 6, 13 20, 2009 FAST AND ACCURATE DIRECTION-OF-ARRIVAL ESTIMATION FOR A SINGLE SOURCE Y. Wu School of Computer Science and Engineering Wuhan Institute of Technology
More informationMedian Filter Based Realizations of the Robust Time-Frequency Distributions
TIME-FREQUENCY SIGNAL ANALYSIS 547 Median Filter Based Realizations of the Robust Time-Frequency Distributions Igor Djurović, Vladimir Katkovnik, LJubiša Stanković Abstract Recently, somenewefficient tools
More informationChapter 6 THE SAMPLING PROCESS 6.1 Introduction 6.2 Fourier Transform Revisited
Chapter 6 THE SAMPLING PROCESS 6.1 Introduction 6.2 Fourier Transform Revisited Copyright c 2005 Andreas Antoniou Victoria, BC, Canada Email: aantoniou@ieee.org July 14, 2018 Frame # 1 Slide # 1 A. Antoniou
More information9. Model Selection. statistical models. overview of model selection. information criteria. goodness-of-fit measures
FE661 - Statistical Methods for Financial Engineering 9. Model Selection Jitkomut Songsiri statistical models overview of model selection information criteria goodness-of-fit measures 9-1 Statistical models
More informationCOMPARATIVE PERFORMANCE ANALYSIS OF THE ALGORITHMS FOR DETECTING PERIODICALLY EXPRESSED GENES. A Thesis KWADWO SEFA AGYEPONG
COMPARATIVE PERFORMANCE ANALYSIS OF THE ALGORITHMS FOR DETECTING PERIODICALLY EXPRESSED GENES A Thesis by KWADWO SEFA AGYEPONG Submitted to the Office of Graduate Studies of Texas A&M University in partial
More informationTHE estimation of covariance matrices is a crucial component
1 A Subspace Method for Array Covariance Matrix Estimation Mostafa Rahmani and George K. Atia, Member, IEEE, arxiv:1411.0622v1 [cs.na] 20 Oct 2014 Abstract This paper introduces a subspace method for the
More informationAdjusted Empirical Likelihood for Long-memory Time Series Models
Adjusted Empirical Likelihood for Long-memory Time Series Models arxiv:1604.06170v1 [stat.me] 21 Apr 2016 Ramadha D. Piyadi Gamage, Wei Ning and Arjun K. Gupta Department of Mathematics and Statistics
More informationDIGITAL FILTERS Analysis, Design, and Applications by A. Antoniou ERRATA. Page 10, Table 1.1: The upper limit of the summation should be K.
DIGITAL FILTERS Analysis, Design, and Applications by A. Antoniou ERRATA Printing #1 Page vii, line 4 : Replace Geophysicists by Geoscientists. Page 1, Table 1.1: The upper limit of the summation should
More informationAdditional Keplerian Signals in the HARPS data for Gliese 667C from a Bayesian re-analysis
Additional Keplerian Signals in the HARPS data for Gliese 667C from a Bayesian re-analysis Phil Gregory, Samantha Lawler, Brett Gladman Physics and Astronomy Univ. of British Columbia Abstract A re-analysis
More informationOn the exponential convergence of. the Kaczmarz algorithm
On the exponential convergence of the Kaczmarz algorithm Liang Dai and Thomas B. Schön Department of Information Technology, Uppsala University, arxiv:4.407v [cs.sy] 0 Mar 05 75 05 Uppsala, Sweden. E-mail:
More informationAutoregressive Models Fourier Analysis Wavelets
Autoregressive Models Fourier Analysis Wavelets BFR Flood w/10yr smooth Spectrum Annual Max: Day of Water year Flood magnitude vs. timing Jain & Lall, 2000 Blacksmith Fork, Hyrum, UT Analyses of Flood
More informationCo-Prime Arrays and Difference Set Analysis
7 5th European Signal Processing Conference (EUSIPCO Co-Prime Arrays and Difference Set Analysis Usham V. Dias and Seshan Srirangarajan Department of Electrical Engineering Bharti School of Telecommunication
More informationDetection theory. H 0 : x[n] = w[n]
Detection Theory Detection theory A the last topic of the course, we will briefly consider detection theory. The methods are based on estimation theory and attempt to answer questions such as Is a signal
More informationParametric Method Based PSD Estimation using Gaussian Window
International Journal of Engineering Trends and Technology (IJETT) Volume 29 Number 1 - November 215 Parametric Method Based PSD Estimation using Gaussian Window Pragati Sheel 1, Dr. Rajesh Mehra 2, Preeti
More informationFourier Analysis of Signals Using the DFT
Fourier Analysis of Signals Using the DFT ECE 535 Lecture April 29, 23 Overview: Motivation Many applications require analyzing the frequency content of signals Speech processing study resonances of vocal
More informationLOW-COMPLEXITY ROBUST DOA ESTIMATION. P.O.Box 553, SF Tampere, Finland 313 Spl. Independenţei, Bucharest, Romania
LOW-COMPLEXITY ROBUST ESTIMATION Bogdan Dumitrescu 1,2, Cristian Rusu 2, Ioan Tăbuş 1, Jaakko Astola 1 1 Dept. of Signal Processing 2 Dept. of Automatic Control and Computers Tampere University of Technology
More informationIntroduction to system identification
Introduction to system identification Jan Swevers July 2006 0-0 Introduction to system identification 1 Contents of this lecture What is system identification Time vs. frequency domain identification Discrete
More informationEvaluating some Yule-Walker Methods with the Maximum-Likelihood Estimator for the Spectral ARMA Model
TEMA Tend. Mat. Apl. Comput., 9, No. 2 (2008), 175-184. c Uma Publicação da Sociedade Brasileira de Matemática Aplicada e Computacional. Evaluating some Yule-Walker Methods with the Maximum-Likelihood
More informationMaximum Likelihood Methods in Radar Array Signal Processing
Maximum Likelihood Methods in Radar Array Signal Processing A. LEE SWINDLEHURST, MEMBER, IEEE, AND PETRE STOICA, FELLOW, IEEE We consider robust and computationally efficient maximum likelihood algorithms
More informationTest for Nonlinear Input Output Relations in SISO Systems by Preliminary Data Analysis
Aalborg Universitet Test for Nonlinear Input Output Relations in SISO Systems by Preliminary Data Analysis Knudsen, Torben Publication date: 000 Document Version Også kaldet Forlagets PDF Link to publication
More informationOn Identification of Cascade Systems 1
On Identification of Cascade Systems 1 Bo Wahlberg Håkan Hjalmarsson Jonas Mårtensson Automatic Control and ACCESS, School of Electrical Engineering, KTH, SE-100 44 Stockholm, Sweden. (bo.wahlberg@ee.kth.se
More informationLong-Run Covariability
Long-Run Covariability Ulrich K. Müller and Mark W. Watson Princeton University October 2016 Motivation Study the long-run covariability/relationship between economic variables great ratios, long-run Phillips
More informationSignal Modeling, Statistical Inference and Data Mining in Astrophysics
ASTRONOMY 6523 Spring 2013 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Course Approach The philosophy of the course reflects that of the instructor, who takes a dualistic view
More information7. Line Spectra Signal is assumed to consists of sinusoidals as
7. Line Spectra Signal is assumed to consists of sinusoidals as n y(t) = α p e j(ω pt+ϕ p ) + e(t), p=1 (33a) where α p is amplitude, ω p [ π, π] is angular frequency and ϕ p initial phase. By using β
More informationMinimax Design of Complex-Coefficient FIR Filters with Low Group Delay
Minimax Design of Complex-Coefficient FIR Filters with Low Group Delay Wu-Sheng Lu Takao Hinamoto Dept. of Elec. and Comp. Engineering Graduate School of Engineering University of Victoria Hiroshima University
More informationLecture 7 Introduction to Statistical Decision Theory
Lecture 7 Introduction to Statistical Decision Theory I-Hsiang Wang Department of Electrical Engineering National Taiwan University ihwang@ntu.edu.tw December 20, 2016 1 / 55 I-Hsiang Wang IT Lecture 7
More informationADAPTIVE ARRAY DETECTION ALGORITHMS WITH STEERING VECTOR MISMATCH
ADAPTIVE ARRAY DETECTIO ALGORITHMS WITH STEERIG VECTOR MISMATCH LIM Chin Heng, Elias Aboutanios Bernard Mulgrew Institute for Digital Communications School of Engineering & Electronics, University of Edinburgh
More informationSystem Identification Approach Applied to Drift Estimation.
Instrumentation and Measurement Technology Conference - IMTC 007 Warsaw, Poland, May 1-3, 007 System Identification Approach Applied to Drift Estimation. Frans Verbeyst 1,, Rik Pintelon 1, Yves Rolain
More informationIMPROVEMENTS IN MODAL PARAMETER EXTRACTION THROUGH POST-PROCESSING FREQUENCY RESPONSE FUNCTION ESTIMATES
IMPROVEMENTS IN MODAL PARAMETER EXTRACTION THROUGH POST-PROCESSING FREQUENCY RESPONSE FUNCTION ESTIMATES Bere M. Gur Prof. Christopher Niezreci Prof. Peter Avitabile Structural Dynamics and Acoustic Systems
More informationGaussian Estimation under Attack Uncertainty
Gaussian Estimation under Attack Uncertainty Tara Javidi Yonatan Kaspi Himanshu Tyagi Abstract We consider the estimation of a standard Gaussian random variable under an observation attack where an adversary
More information