VIII. Coherence and Transfer Function Applications A. Coherence Function Estimates

Size: px
Start display at page:

Download "VIII. Coherence and Transfer Function Applications A. Coherence Function Estimates"

Transcription

1 VIII. Coherence and Transfer Function Applications A. Coherence Function Estimates Consider the application of these ideas to the specific problem of atmospheric turbulence measurements outlined in Figure VIII.. acceleration Output v(t) Vertical C.. Input u(t) Vertical gust velocity Figure VIII.. Airplane flying through atmospheric turbulence, where u(t) measured vertical gust velocity and v(t) measured vertical acceleration (acceleration in g's). Consider the spectral energy density and coherence functions (see Figures VIII. & VIII.3) for T 0 minute (600 sec) records of u(t) and v(t) in with a fundamental frequency band f 0.05 cps over the frequency range 0. to 4.0 cps. What are the degrees of freedom for these calculations? What is the zero coherence confidence level? Coherence Analysis High Coherence ( ) in the range 0.3 < f <.0 cps Low f < 0.3 cps coherence low due to maneer loads induced by pilot, rather than atmospheric turbulence! i.e. contributions to y(t) from other f >.0 cps low coherence due to reduced energy in response relative to measurement noise. Response energy decreases at higher frequencies due to decreasing plane response to turbulence and reduced turbulence energy. The ordinary coherence squared, defined as, (VIII.) where 0 is a statistical estimator related to the frequency response function as described next. 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown

2 Figure VIII. Power spectra of gust velocity and response acceleration. (Data from studies funded by the NASA Langley Research Center, Hampton, Virginia; Contract NAS -8538). Figure VIII.3 Coherence function between gust velocity and response acceleration. 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown

3 ********************************************************************* EXAMPLES - The Extremes ********************************************************************* Extreme : Consider the ideal noise free, constant parameter system in Figure VIII.4, for which H and H. (VIII.) Figure VIII.4 In this model, a noise-free input produces a noise-free output. Thus the coherence squared is H (VIII.3) H Conclusion: The coherence between input and output in a noise-free system is ALWAYS equal to Extreme : Consider a "real world" situation in which the output v(t) is completely unrelated to the input u(t). Thus 0. Normally falls in the range 0 for one or more of the following reasons: (a) There is extraneous noise in the measurement of the input and output; (b) The system relating u(t) and v(t) is nonlinear; and/or (c) The output v(t) is due to u(t) as well as other unmeasured inputs. In these Anormal situations@ for a linear system (at each frequency f); can be thought of as the fractional part of the mean square output v(t ) due to u(t). Conversely the quantity - is the fractional part of the mean square output error due to sources other than u(t); and thus represents the level of noise in the system. To explore these ideas, consider the linear process in Figure VIII.5.in which () noise n(t) is added to the noise-free input u(t) at the input; and () noise m(t) is added to the noise-free output v(t) at the output. 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 3

4 Figure VIII.5 Schematic of linear system in which distinctly different noises are added to the noise-free input and output respectively. Under these circumstances, the respective measured input and output are u (t) v (t) u(t) + n(t) v(t) + m(t) And the respective auto- and cross-spectral density functions for the system are u nn mm um un vm nv nu mv nm, (VIII.4) (VIII.5a) where H (VIII..5b) and H. (VIII..5c) Note that all possible correlations among signals and noise are allowed in this general case. So lets consider these special cases of interest. Case. No input noise [i.e. n(t) 0]; Output noise is uncorrelated (See Figure VIII.6) Figure VIII.6 Noise-free input and uncorrelated output noise. Thus Eqs. (VIII.5) are reduced to the following set of equations (VIII.6, 7, 8 and 9): ' u ( t ) u ( t ) (VIII.6a) ' v ( t ) v ( t ) + m ( t ) (VIII.6b) where ( f ) 0 vm u (VIII.7a) 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 4

5 v + mm (VIII.7b) H H u (VIII.8a) H (VIII.8b) H (VIII.8) (VIII.9a) (VIII.9b) Thus the noise-free output spectrum can be computed using only the measured input and output ( i.e. u and! ). Therefore the spectrum of the noise at the output mm - mm can also be computed according to. This case is important because in reality () Inputs can be defined; and () Input noise can be minimized to simulate the conditions of the case. However we have no control over output noise because of () non-linearities and/or () unmeasured inputs! The ordinary coherence squared function for this case is (VIII.0) [ + mm ] u + mm, (VIII.) which compares to the previous noise-free case in which 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 5

6 Thus for > 0 mm. However because of (VIII.9) u < (VIII.) So we can use measurements to compute the Acoherent output spectrum@ or Further from above (t) mm and we can compute the spectrum of the output noise (VIII.3) - (- mm v ) (VIII.4) ************************************************************** Case. No output noise [m(t) 0]; Input noise is uncorrelated (See Figure VIII.7) Figure VIII.7. Uncorrelated noise is added only to the input. (t) u(t) + n(t), (VIII.5a) with un 0. (t) v(t) (VIII.5b) u + nn v H (VIII.6a) 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 6

7 H * [ ] H (VIII.6b) u Alternatively it can be shown from the definitions of From above > (VIII.8) H v H * H S u that. (VIII.7) * H and can be computed from u (t), (t). Also since and > v H > nn u - (VIII.9) (VIII.0) The spectra of both the noise-free input and the noise at the input can be computed. However, do not apply this model if output noise is expected. The ordinary coherence squared in this case is v [ + nn] (VIII.a) + nn (VIII.b) Thus, like case, the ordinary coherence squared function can be used to determine the spectrum of the noise-free input according to (VIII.) AND the input noise spectrum according to (- ) nn (VIII.3) 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 7

8 The Effects of eneral Extraneous Noise Under ideal conditions, the measured spectral density function u + nn (VIII.4a) v + mm (VIII.4b) u H H (VIII.5) u H v - nn (VIII.6) Note H u v but need to know the noise. H f( nn ), - mm - nn not f( mm ), H f( mm, mn ) The coherence squared, based on measurements, is always less than the true coherence as shown below û ˆ (VIII.7) [ + nn ][ + mm] + nn + mm nn mm + < Set output noise mm 0. B. Frequency Response Function Estimates Consider a single-input output linear system u(t) v(t) 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 8

9 v(t) 0 h( τ ) u(t -τ ) dτ h( τ ) 0 for τ < 0 It follows that i.e. physically realizable v(t) v(t + τ ) Taking expected values of both sides (zero Similarly R R 0 h( α ) h( β ) u(t - β ) u(t + τ -α ) dα d β u, v ) input/output auto-correlation function (t) 0 ( τ ) h( α) h( β ) R 0 h( α) R Fourier Transforms of the above in terms of one-sided-spectra ( τ -α) dα ( τ + β -α ) dα d β H _ H (VIII.8a) H > H (VIII.8b) Consider the ratio of the two different ways of computing the squared gain factor: H H (VIII.9) Different estimates of the squared gain factor: () Ĥ ) (VIII.30a) which is a biased estimate except at 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 9

10 () Ĥ (VIII.30b) which is a biased estimate if there is input noise n(t) However, if m(t) is the only output noise, then the ratio of () and () is an unbiased estimate. Ĥ Ĥ (VIII.30c) The sample frequency response function is Ĥ Ĥ e +i ˆ θ (VIII.3) the bias error sources in this case are: () Inherent in the estimation process; () Nonlinear and/or time varying system parameters (i.e. non-stationarity); (3) Bias in the power- and cross-spectral density estimates; (4) Extraneous noise measured at input (incoherent noise does not cause bias); (5) Contributions to the measured output from other unmeasured inputs correlated with the measured input; (6) Bias due to propagation time delays. Remarks: () In general there is bias in the estimation process. That is } E E{ E{ } H } E{ Ĥ (VIII.3) E{ Ĥ} H, thus Ĥ is a biased estimator as shown below. 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 0

11 In practice this source of error can be made relatively small with a good choice of degrees of freedom. () The linearity assumption is often violated if the operating range is sufficiently wide. Nevertheless it should be noted that the application of equation (a) above to nonlinear systems represents the "best" linear approximation (in the least squares sense) of the frequency response function of the system under the specified input and output conditions. (3) Note: The bias of spectral density estimates can be large in and around spectral peaks. However, these errors can be suppressed by choosing f s small enough to define spectral peaks. (4) In the case, where extraneous noise at input is uncorrelated with input or output signal (i.e. Case ), (see Figure VIII.8). Figure VIII.8 Extraneous noise at the input is uncorrelated with input or output signal. Consider input signal u(t) + n(t) (t) for the case in which n(t) passes through system, but u(t) E{ u } (VIII.33) E{ } E{H } H (VIII.34) E{ } + nn + nn does. If all other bias errors are ignored, then with the normalized bias error E{Ĥ } nn ε b - + nn -- (VIII.35) H mm + nn downward bias Thus the estimated frequency response function, Ĥ, is biased downwards by an amount eb equal to the noise to signal ratio at input e.g. N/S 0.0 or (5) When noise m(t) at output is correlated with measured input u(t) (Figure VIII.9) 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown

12 Figure VIII.9 Noise at the output m(t) is correlated with the measured input u(t); u(t) is unmeasured input. u(t) - measured input; u'(t) - an unmeasured input; v(t) - output; and v'(t) - measured output, and from the definition of second edition of Bendat and Piersol) that Ĥ above, it can be shown (Eq. 7.3 of section 7... of the u u (VIII.36) If all other bias errors are ignored, then E{H} H - u v - u u v H u [- ] [- u u u ] v (VIII.37) Note, further that if u and are uncorrelated, then u u u 0 (VIII.38) E{Ĥ} H which demonstrates that unmeasured uncorrelated noise at the input causes no bias. 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown

13 Random errors (see Figure VIII.0) If H is unknown Figure VIII.0 A system with noise z(t) at measured output that is uncorrelated with input or output. Let u(t) - input v(t) - output z(t) - noise uncorrelated with input or output (i.e. residual extraneous noise at the output related only to the measurement only) v'(t) measured output z(t) v (t) - v(t) (VIII.39) The Finite Fourier Transform for z(t) is z(t) v (t) - 0 h( τ ) u(t -τ ) dτ Z(f, T) V (f,t) - H U(f, T) + [- Ĥ U(f, T) + Ĥ U(f,T)] where Ĥ v - a biased estimator or Z(f, T) Ẑ(f,T) + U(f,T) [Ĥ(f,T) - H], (VIII.40) where Ẑ (f,t) V (f,t) - Ĥ(f,T) U(f, T) (VIII.4) with Ẑ best estimate of noise It follows that + Ĥ - H zz ẑẑ (VIII.4) 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 3

14 If (VIII.4) is smoothed so that all the terms, including the ẑ ẑ estimate, have n f s T degrees of freedom (DOF), then it will be distributed according to n zz χ n (VIII.43) zz The two terms on the RHS of (VIII.4) are statistically independent, thus their DOF will be additive according to χ n χ n n + n-n χ n n where DOF for term (a) of (VIII.4) n is reduced by number of constraints on From (VIII.40) and (VIII.4) Ẑ (f,t) F(Ĥ), ẑ ẑ. where Ẑ(f, T) is defined in terms of Ĥ (not H), which implies independent constraints on Ẑ ; one each for the real and imaginary part. Hence ẑẑ will have less DOF than zz, thus n zz χ n zz (VIII.44a) n Ĥ - H zz y n χ (VIII.44b) ẑẑ v v u v but we know [- ] so (VIII.44a) becomes n [- ] v v u v χ n- zz y n- (VIII.45) 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 4

15 F,n n y / n y / n n n n- F,n- rˆ Ĥ - H F n- rˆ (n- ) Ĥ - H (VIII.46) [- ˆu v ] v v [- ˆ ],n-; v n ft; n- α (VIII.47) _ ( ft-) ft For a ratio of estimators with Χ n distributions we have a new distribution From (VIII.46), if H is unknown, then the confidence intervals for H with confidence coefficient -α, can be determined by a quantity rˆ such that the variance Var[ Ĥ] rˆ H is a ratio of χ distributed estimates so the appropriate distribution function is v F,n-;α input spectrum 00 α percentage point of an F-distribution with n and n n- degrees of freedom, f n f measured output spectrum s T T f Ts s ˆ v sample estimate of the ordinary coherence between u(t) and v (t) Approximate (-α) confidence intervals [ Ĥ - rˆ H Ĥ + rˆ ] [ ˆ φ - ˆ φ φ ˆ φ + ˆ] φ (VIII.48) (See B&P for Tables) ˆ - φ sin rˆ Ĥ 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 5

16 C. Estimating Tidal Transfer Function Error Assume that the observed tide v'(t) is the outcome of a linear process, characterized by the true transfer function H and corresponding weighting function h(t), plus some uncorrelated random noise at the output m(t) as shown below. u(t) h(t) /H v(t) + v'(t) (VIII.49a) m(t) The input u(t) in this case is the noise-free astronomical forcing, which is specified in terms of the equilibrium tide. The measured response of this system v'(t) is the sum of the noise-free tidal response v(t) and the noise at the output m(t) according to (t) v(t) + m(t) (VIII.49b) The tidal response v(t) is described here in terms of the harmonic constants, (amplitude H n and phase K n respectively) for each tidal constituent of the equilibrium tide (with it s unique frequency f n ). Here we seek an estimate of the uncertainty in the estimates of H n, K n in the presence of noise. In the frequency domain (VIII.49b) becomes v + where the 's are one-sided spectral density functions. We have shown that since mm v(t) 0 h( β ) u(t- β ) d β, H (VIII.50) A corresponding relation is u ' where Ĥ is a transfer function estimate. Ĥ v, (VIII.5) 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 6

17 Here is the complex one-sided cross-spectral density function, C + i Q, (VIII.5) where C and Q are the co- and quadrature-spectrum respectively. In polar notation, Eq. (VIII.50) becomes where and Thus and -iθ -i Φ e (C tan H e + Q [Q ) / C / ] - θ. H θ Φ., (VIII.53) Since is noise-free, the uncertainty in computing estimates of, θ uy and H, Φ is equivalent. There are several references on estimating uncertainties for H and Φ. From the Appendix B of Munk and Cartwright (966), transfer function estimates true transfer function T and a random noise according to Here / Ĥ e -iν mm u v ' / Ĥ are written in terms of the T + ε H e -i ε is a complex random variable, whose (a) real and + ε. imaginary parts have approximate normal probability distributions with zero mean and (b) whose variance in the fundamental frequency band f is given by / [ mm f ]/ [ f ] 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 7 p σ, where p f e T; in which f e is the effective bandwidth used in the spectral computation and T is the length of series. (Note that [ mm f]/p is the mean square error in estimating a mean square value in the presence of band limited white noise).

18 We normalizeσ / by the magnitude of the true transfer function H / and find ' '' σ σ H ( f ) [ mm f ]/ [ f ] They compute the probability density function and the 95% confidence limits for the normalized variables Ĥ /H and β Φˆ - Φ in terms of σ. (See Figure VIII.). u v p, Figure VIII.The upper diagrams show the probability distributions of (left) ρ sample admittance/true admittance, and (right) θ sample phase - true phase, for stated values of the noise parameter σ. Lower diagrams show the 95% confidence limits of the same quantities and also the mean r.m.s. values of ρ, all plotted against σ. (From Munk and Cartwright, 966). 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 8

19 D. Summary of Error Relations Error Relations (a) Sample Mean --- is unbiased Random Error: In general Var [u] T T -T when τ» T; T τ - R T 0 T ( τ ) dτ ; Var [u] - R ( τ ) dτ Normalized Mean Square Error: Band-limited white noise (BT > 5) ε u σ u BT µ where B and τ o is the integral time scale τ o (b) Sample Mean Square --- is unbiased Var [u ] T T T -T - τ - (R T (R ( τ ) + µ R ( τ ) + µ R u u ( τ )) dτ ( τ )) dτ ; whenτ Normalized Mean Square Error: Band limited white noise (BT 5) max» T Var u [ ^] R (O) BT τ o ε BT T 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 9

20 (c) Sample Cross-Covariance --- is unbiased Var[ Rˆ ( τ )] T T T α - [R -T T [R - ( α) R ( α) R ( α) + R ( α) + R ( α + τ ) R ( α + τ ) R vu vu ( α -τ )] dα ( α -τ )] dα ; for τ» T Normalized Mean Square Error: Band limited white noise BT > 5 ε BT + R (0) R (0), R ( τ ) (d) where Sample Spectral Density Var[ ] b[ ft " -8π - τ R ( τ ) e f ] 4 -iπ f τ dτ " Normalized mean square error E{ - } ε (-α) confidence interval ft 4 f " n χ n; α / n < χ Similar for cross-spectral density. n;-α / n f s T (e) Sample Coherence ˆ / + ˆ - w ln tan h - 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown 0

21 (-α) confidence interval [tanh w - (n- ) - - σ w α / z ) < - tanh (w - (n- ) + σ w zα / ] Frequency Response: Ĥ (a biased ratio of Χ distributed variables) rˆ F n- (-α) confidence interval [- ˆ,n-; α ] [( Ĥ - rˆ) H ( Ĥ + rˆ)] [ ˆ θ - ˆ θ θ ˆ θ + ˆ], θ where ˆ - θ sin rˆ Ĥ 7 February 008 Chapter VIII Coherence & Transfer Apps 8Wendell S. Brown

System Identification & Parameter Estimation

System Identification & Parameter Estimation System Identification & Parameter Estimation Wb3: SIPE lecture Correlation functions in time & frequency domain Alfred C. Schouten, Dept. of Biomechanical Engineering (BMechE), Fac. 3mE // Delft University

More information

Time series models in the Frequency domain. The power spectrum, Spectral analysis

Time series models in the Frequency domain. The power spectrum, Spectral analysis ime series models in the Frequency domain he power spectrum, Spectral analysis Relationship between the periodogram and the autocorrelations = + = ( ) ( ˆ α ˆ ) β I Yt cos t + Yt sin t t= t= ( ( ) ) cosλ

More information

Basic Descriptions and Properties

Basic Descriptions and Properties CHAPTER 1 Basic Descriptions and Properties This first chapter gives basic descriptions and properties of deterministic data and random data to provide a physical understanding for later material in this

More information

ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process

ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process Department of Electrical Engineering University of Arkansas ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process Dr. Jingxian Wu wuj@uark.edu OUTLINE 2 Definition of stochastic process (random

More information

X. Cross Spectral Analysis

X. Cross Spectral Analysis X. Cross Spectral Analysis Cross Spectrum We have already dealt with the crosscovariance function (ccvf): ˆ 1 C (k) = X(i)Y(i + k) N N k i= 0 The Fourier transform of the ccvf is called the cross spectrum.

More information

Randomly Modulated Periodic Signals

Randomly Modulated Periodic Signals Randomly Modulated Periodic Signals Melvin J. Hinich Applied Research Laboratories University of Texas at Austin hinich@mail.la.utexas.edu www.la.utexas.edu/~hinich Rotating Cylinder Data Fluid Nonlinearity

More information

Stochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno

Stochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno Stochastic Processes M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno 1 Outline Stochastic (random) processes. Autocorrelation. Crosscorrelation. Spectral density function.

More information

Modern Navigation. Thomas Herring

Modern Navigation. Thomas Herring 12.215 Modern Navigation Thomas Herring Estimation methods Review of last class Restrict to basically linear estimation problems (also non-linear problems that are nearly linear) Restrict to parametric,

More information

Review of Doppler Spread The response to exp[2πift] is ĥ(f, t) exp[2πift]. ĥ(f, t) = β j exp[ 2πifτ j (t)] = exp[2πid j t 2πifτ o j ]

Review of Doppler Spread The response to exp[2πift] is ĥ(f, t) exp[2πift]. ĥ(f, t) = β j exp[ 2πifτ j (t)] = exp[2πid j t 2πifτ o j ] Review of Doppler Spread The response to exp[2πift] is ĥ(f, t) exp[2πift]. ĥ(f, t) = β exp[ 2πifτ (t)] = exp[2πid t 2πifτ o ] Define D = max D min D ; The fading at f is ĥ(f, t) = 1 T coh = 2D exp[2πi(d

More information

SIO 221B, Rudnick adapted from Davis 1. 1 x lim. N x 2 n = 1 N. { x} 1 N. N x = 1 N. N x = 1 ( N N x ) x = 0 (3) = 1 x N 2

SIO 221B, Rudnick adapted from Davis 1. 1 x lim. N x 2 n = 1 N. { x} 1 N. N x = 1 N. N x = 1 ( N N x ) x = 0 (3) = 1 x N 2 SIO B, Rudnick adapted from Davis VII. Sampling errors We do not have access to the true statistics, so we must compute sample statistics. By this we mean that the number of realizations we average over

More information

Statistical errors in the estimation of time-averaged acoustic energy density using the two-microphone method

Statistical errors in the estimation of time-averaged acoustic energy density using the two-microphone method Statistical errors in the estimation of time-averaged acoustic energy density using the two-microphone method Short title: Statistical errors in acoustic energy density Justin Ghan, Ben Cazzolato, and

More information

Sinusoidal Modeling. Yannis Stylianou SPCC University of Crete, Computer Science Dept., Greece,

Sinusoidal Modeling. Yannis Stylianou SPCC University of Crete, Computer Science Dept., Greece, Sinusoidal Modeling Yannis Stylianou University of Crete, Computer Science Dept., Greece, yannis@csd.uoc.gr SPCC 2016 1 Speech Production 2 Modulators 3 Sinusoidal Modeling Sinusoidal Models Voiced Speech

More information

Econometrics II - EXAM Outline Solutions All questions have 25pts Answer each question in separate sheets

Econometrics II - EXAM Outline Solutions All questions have 25pts Answer each question in separate sheets Econometrics II - EXAM Outline Solutions All questions hae 5pts Answer each question in separate sheets. Consider the two linear simultaneous equations G with two exogeneous ariables K, y γ + y γ + x δ

More information

ECE 636: Systems identification

ECE 636: Systems identification ECE 636: Systems identification Lectures 7 8 onparametric identification (continued) Important distributions: chi square, t distribution, F distribution Sampling distributions ib i Sample mean If the variance

More information

Stochastic Dynamics of SDOF Systems (cont.).

Stochastic Dynamics of SDOF Systems (cont.). Outline of Stochastic Dynamics of SDOF Systems (cont.). Weakly Stationary Response Processes. Equivalent White Noise Approximations. Gaussian Response Processes as Conditional Normal Distributions. Stochastic

More information

IMPROVEMENTS IN MODAL PARAMETER EXTRACTION THROUGH POST-PROCESSING FREQUENCY RESPONSE FUNCTION ESTIMATES

IMPROVEMENTS 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 information

13.42 READING 6: SPECTRUM OF A RANDOM PROCESS 1. STATIONARY AND ERGODIC RANDOM PROCESSES

13.42 READING 6: SPECTRUM OF A RANDOM PROCESS 1. STATIONARY AND ERGODIC RANDOM PROCESSES 13.42 READING 6: SPECTRUM OF A RANDOM PROCESS SPRING 24 c A. H. TECHET & M.S. TRIANTAFYLLOU 1. STATIONARY AND ERGODIC RANDOM PROCESSES Given the random process y(ζ, t) we assume that the expected value

More information

Econometrics Summary Algebraic and Statistical Preliminaries

Econometrics Summary Algebraic and Statistical Preliminaries Econometrics Summary Algebraic and Statistical Preliminaries Elasticity: The point elasticity of Y with respect to L is given by α = ( Y/ L)/(Y/L). The arc elasticity is given by ( Y/ L)/(Y/L), when L

More information

An Introduction to Bayesian Linear Regression

An Introduction to Bayesian Linear Regression An Introduction to Bayesian Linear Regression APPM 5720: Bayesian Computation Fall 2018 A SIMPLE LINEAR MODEL Suppose that we observe explanatory variables x 1, x 2,..., x n and dependent variables y 1,

More information

INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -12 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.

INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -12 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -12 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Data Extension & Forecasting Moving

More information

2. (a) What is gaussian random variable? Develop an equation for guassian distribution

2. (a) What is gaussian random variable? Develop an equation for guassian distribution Code No: R059210401 Set No. 1 II B.Tech I Semester Supplementary Examinations, February 2007 PROBABILITY THEORY AND STOCHASTIC PROCESS ( Common to Electronics & Communication Engineering, Electronics &

More information

Lecture 2. Fading Channel

Lecture 2. Fading Channel 1 Lecture 2. Fading Channel Characteristics of Fading Channels Modeling of Fading Channels Discrete-time Input/Output Model 2 Radio Propagation in Free Space Speed: c = 299,792,458 m/s Isotropic Received

More information

Non-parametric identification

Non-parametric identification Non-parametric Non-parametric Transient Step-response using Spectral Transient Correlation Frequency function estimate Spectral System Identification, SSY230 Non-parametric 1 Non-parametric Transient Step-response

More information

Fourier Analysis Linear transformations and lters. 3. Fourier Analysis. Alex Sheremet. April 11, 2007

Fourier Analysis Linear transformations and lters. 3. Fourier Analysis. Alex Sheremet. April 11, 2007 Stochastic processes review 3. Data Analysis Techniques in Oceanography OCP668 April, 27 Stochastic processes review Denition Fixed ζ = ζ : Function X (t) = X (t, ζ). Fixed t = t: Random Variable X (ζ)

More information

Module 29.3: nag tsa spectral Time Series Spectral Analysis. Contents

Module 29.3: nag tsa spectral Time Series Spectral Analysis. Contents Time Series Analysis Module Contents Module 29.3: nag tsa spectral Time Series Spectral Analysis nag tsa spectral calculates the smoothed sample spectrum of a univariate and bivariate time series. Contents

More information

Scattering of Electromagnetic Radiation. References:

Scattering of Electromagnetic Radiation. References: Scattering of Electromagnetic Radiation References: Plasma Diagnostics: Chapter by Kunze Methods of experimental physics, 9a, chapter by Alan Desilva and George Goldenbaum, Edited by Loveberg and Griem.

More information

Wind Loads. 2 nd Part Conversion of wind speeds to wind loads, aerodynamical coefficients, wind pressure and suction, and practical design methods

Wind Loads. 2 nd Part Conversion of wind speeds to wind loads, aerodynamical coefficients, wind pressure and suction, and practical design methods Wind Loads nd Part Conversion of wind speeds to wind loads, aerodynamical coefficients, wind pressure and suction, and practical design methods March 007 Júlíus Sólnes, Professor of civil and environmental

More information

Resampling techniques for statistical modeling

Resampling techniques for statistical modeling Resampling techniques for statistical modeling Gianluca Bontempi Département d Informatique Boulevard de Triomphe - CP 212 http://www.ulb.ac.be/di Resampling techniques p.1/33 Beyond the empirical error

More information

Chapter 6: Nonparametric Time- and Frequency-Domain Methods. Problems presented by Uwe

Chapter 6: Nonparametric Time- and Frequency-Domain Methods. Problems presented by Uwe System Identification written by L. Ljung, Prentice Hall PTR, 1999 Chapter 6: Nonparametric Time- and Frequency-Domain Methods Problems presented by Uwe System Identification Problems Chapter 6 p. 1/33

More information

TSKS01 Digital Communication Lecture 1

TSKS01 Digital Communication Lecture 1 TSKS01 Digital Communication Lecture 1 Introduction, Repetition, and Noise Modeling Emil Björnson Department of Electrical Engineering (ISY) Division of Communication Systems Emil Björnson Course Director

More information

The statistics of ocean-acoustic ambient noise

The statistics of ocean-acoustic ambient noise The statistics of ocean-acoustic ambient noise Nicholas C. Makris Naval Research Laboratory, Washington, D.C. 0375, USA Abstract With the assumption that the ocean-acoustic ambient noise field is a random

More information

Matrix power converters: spectra and stability

Matrix power converters: spectra and stability Matrix power converters: spectra and stability Stephen Cox School of Mathematical Sciences, University of Nottingham supported by EPSRC grant number EP/E018580/1 Making It Real Seminar, Bristol 2009 Stephen

More information

6.435, System Identification

6.435, System Identification System Identification 6.435 SET 3 Nonparametric Identification Munther A. Dahleh 1 Nonparametric Methods for System ID Time domain methods Impulse response Step response Correlation analysis / time Frequency

More information

Lecture 1: Introduction to System Modeling and Control. Introduction Basic Definitions Different Model Types System Identification

Lecture 1: Introduction to System Modeling and Control. Introduction Basic Definitions Different Model Types System Identification Lecture 1: Introduction to System Modeling and Control Introduction Basic Definitions Different Model Types System Identification What is Mathematical Model? A set of mathematical equations (e.g., differential

More information

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture 6-5 November, 2015

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture 6-5 November, 2015 The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture 6-5 November, 015 Institute of Structural Engineering Method of Finite Elements II 1 Introduction

More information

Computer Science, Informatik 4 Communication and Distributed Systems. Simulation. Discrete-Event System Simulation. Dr.

Computer Science, Informatik 4 Communication and Distributed Systems. Simulation. Discrete-Event System Simulation. Dr. Simulation Discrete-Event System Simulation Chapter 0 Output Analysis for a Single Model Purpose Objective: Estimate system performance via simulation If θ is the system performance, the precision of the

More information

NONLINEAR INTEGRAL MINIMUM VARIANCE-LIKE CONTROL WITH APPLICATION TO AN AIRCRAFT SYSTEM

NONLINEAR INTEGRAL MINIMUM VARIANCE-LIKE CONTROL WITH APPLICATION TO AN AIRCRAFT SYSTEM NONLINEAR INTEGRAL MINIMUM VARIANCE-LIKE CONTROL WITH APPLICATION TO AN AIRCRAFT SYSTEM D.G. Dimogianopoulos, J.D. Hios and S.D. Fassois DEPARTMENT OF MECHANICAL & AERONAUTICAL ENGINEERING GR-26500 PATRAS,

More information

ENSC327 Communications Systems 2: Fourier Representations. Jie Liang School of Engineering Science Simon Fraser University

ENSC327 Communications Systems 2: Fourier Representations. Jie Liang School of Engineering Science Simon Fraser University ENSC327 Communications Systems 2: Fourier Representations Jie Liang School of Engineering Science Simon Fraser University 1 Outline Chap 2.1 2.5: Signal Classifications Fourier Transform Dirac Delta Function

More information

GATE EE Topic wise Questions SIGNALS & SYSTEMS

GATE EE Topic wise Questions SIGNALS & SYSTEMS www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)

More information

Signal Design for Band-Limited Channels

Signal Design for Band-Limited Channels Wireless Information Transmission System Lab. Signal Design for Band-Limited Channels Institute of Communications Engineering National Sun Yat-sen University Introduction We consider the problem of signal

More information

Chapter 11. Output Analysis for a Single Model Prof. Dr. Mesut Güneş Ch. 11 Output Analysis for a Single Model

Chapter 11. Output Analysis for a Single Model Prof. Dr. Mesut Güneş Ch. 11 Output Analysis for a Single Model Chapter Output Analysis for a Single Model. Contents Types of Simulation Stochastic Nature of Output Data Measures of Performance Output Analysis for Terminating Simulations Output Analysis for Steady-state

More information

Communication Systems Lecture 21, 22. Dong In Kim School of Information & Comm. Eng. Sungkyunkwan University

Communication Systems Lecture 21, 22. Dong In Kim School of Information & Comm. Eng. Sungkyunkwan University Communication Systems Lecture 1, Dong In Kim School of Information & Comm. Eng. Sungkyunkwan University 1 Outline Linear Systems with WSS Inputs Noise White noise, Gaussian noise, White Gaussian noise

More information

Module 4. Signal Representation and Baseband Processing. Version 2 ECE IIT, Kharagpur

Module 4. Signal Representation and Baseband Processing. Version 2 ECE IIT, Kharagpur Module Signal Representation and Baseband Processing Version ECE II, Kharagpur Lesson 8 Response of Linear System to Random Processes Version ECE II, Kharagpur After reading this lesson, you will learn

More information

Virtual Array Processing for Active Radar and Sonar Sensing

Virtual Array Processing for Active Radar and Sonar Sensing SCHARF AND PEZESHKI: VIRTUAL ARRAY PROCESSING FOR ACTIVE SENSING Virtual Array Processing for Active Radar and Sonar Sensing Louis L. Scharf and Ali Pezeshki Abstract In this paper, we describe how an

More information

GMM, HAC estimators, & Standard Errors for Business Cycle Statistics

GMM, HAC estimators, & Standard Errors for Business Cycle Statistics GMM, HAC estimators, & Standard Errors for Business Cycle Statistics Wouter J. Den Haan London School of Economics c Wouter J. Den Haan Overview Generic GMM problem Estimation Heteroskedastic and Autocorrelation

More information

Lecture 11: Spectral Analysis

Lecture 11: Spectral Analysis Lecture 11: Spectral Analysis Methods For Estimating The Spectrum Walid Sharabati Purdue University Latest Update October 27, 2016 Professor Sharabati (Purdue University) Time Series Analysis October 27,

More information

Prof. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis

Prof. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis Introduction to Time Series Analysis 1 Contents: I. Basics of Time Series Analysis... 4 I.1 Stationarity... 5 I.2 Autocorrelation Function... 9 I.3 Partial Autocorrelation Function (PACF)... 14 I.4 Transformation

More information

A6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2013

A6523 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 information

STATISTICAL ORBIT DETERMINATION

STATISTICAL ORBIT DETERMINATION STATISTICAL ORBIT DETERMINATION Satellite Tracking Example of SNC and DMC ASEN 5070 LECTURE 6 4.08.011 1 We will develop a simple state noise compensation (SNC) algorithm. This algorithm adds process noise

More information

Analytical Methods. Session 3: Statistics II. UCL Department of Civil, Environmental & Geomatic Engineering. Analytical Methods.

Analytical Methods. Session 3: Statistics II. UCL Department of Civil, Environmental & Geomatic Engineering. Analytical Methods. Analytical Methods Session 3: Statistics II More statistical tests Quite a few more questions that we might want to ask about data that we have. Is one value significantly different to the rest, or to

More information

Short-Term Wave Analysis

Short-Term Wave Analysis Chapter 3 Short-Term Wave Analysis 3.1 Introduction In analysis of wave data, it is important to distinguish between shortterm and long-term wave analysis. Short-term analysis refers to analysis of waves

More information

EAS 305 Random Processes Viewgraph 1 of 10. Random Processes

EAS 305 Random Processes Viewgraph 1 of 10. Random Processes EAS 305 Random Processes Viewgraph 1 of 10 Definitions: Random Processes A random process is a family of random variables indexed by a parameter t T, where T is called the index set λ i Experiment outcome

More information

Calculating Accurate Waveforms for LIGO and LISA Data Analysis

Calculating Accurate Waveforms for LIGO and LISA Data Analysis Calculating Accurate Waveforms for LIGO and LISA Data Analysis Lee Lindblom Theoretical Astrophysics, Caltech HEPL-KIPAC Seminar, Stanford 17 November 2009 Results from the Caltech/Cornell Numerical Relativity

More information

MATH 251 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam

MATH 251 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam MATH 51 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam A collection of previous exams could be found at the coordinator s web: http://www.math.psu.edu/tseng/class/m51samples.html

More information

Separation of a Signal of Interest from a Seasonal Effect in Geophysical Data: I. El Niño/La Niña Phenomenon

Separation of a Signal of Interest from a Seasonal Effect in Geophysical Data: I. El Niño/La Niña Phenomenon International Journal of Geosciences, 2011, 2, **-** Published Online November 2011 (http://www.scirp.org/journal/ijg) Separation of a Signal of Interest from a Seasonal Effect in Geophysical Data: I.

More information

Signal processing Frequency analysis

Signal processing Frequency analysis Signal processing Frequency analysis Jean-Hugh Thomas (jean-hugh.thomas@univ-lemans.r) Fourier series and Fourier transorm (h30 lecture+h30 practical work) 2 Sampling (h30+h30) 3 Power spectrum estimation

More information

HST 583 FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA ANALYSIS AND ACQUISITION A REVIEW OF STATISTICS FOR FMRI DATA ANALYSIS

HST 583 FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA ANALYSIS AND ACQUISITION A REVIEW OF STATISTICS FOR FMRI DATA ANALYSIS HST 583 FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA ANALYSIS AND ACQUISITION A REVIEW OF STATISTICS FOR FMRI DATA ANALYSIS EMERY N. BROWN AND CHRIS LONG NEUROSCIENCE STATISTICS RESEARCH LABORATORY DEPARTMENT

More information

System Modeling and Identification CHBE 702 Korea University Prof. Dae Ryook Yang

System Modeling and Identification CHBE 702 Korea University Prof. Dae Ryook Yang System Modeling and Identification CHBE 702 Korea University Prof. Dae Ryook Yang 1-1 Course Description Emphases Delivering concepts and Practice Programming Identification Methods using Matlab Class

More information

Introduction to Probability and Stochastic Processes I

Introduction to Probability and Stochastic Processes I Introduction to Probability and Stochastic Processes I Lecture 3 Henrik Vie Christensen vie@control.auc.dk Department of Control Engineering Institute of Electronic Systems Aalborg University Denmark Slides

More information

ENSC327 Communications Systems 2: Fourier Representations. School of Engineering Science Simon Fraser University

ENSC327 Communications Systems 2: Fourier Representations. School of Engineering Science Simon Fraser University ENSC37 Communications Systems : Fourier Representations School o Engineering Science Simon Fraser University Outline Chap..5: Signal Classiications Fourier Transorm Dirac Delta Function Unit Impulse Fourier

More information

2A1H Time-Frequency Analysis II

2A1H Time-Frequency Analysis II 2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 209 For any corrections see the course page DW Murray at www.robots.ox.ac.uk/ dwm/courses/2tf. (a) A signal g(t) with period

More information

A note on the uniform perturbation index 1

A note on the uniform perturbation index 1 Rostock. Math. Kolloq. 52, 33 46 (998) Subject Classification (AMS) 65L5 M. Arnold A note on the uniform perturbation index ABSTRACT. For a given differential-algebraic equation (DAE) the perturbation

More information

Electrical Engineering Written PhD Qualifier Exam Spring 2014

Electrical Engineering Written PhD Qualifier Exam Spring 2014 Electrical Engineering Written PhD Qualifier Exam Spring 2014 Friday, February 7 th 2014 Please do not write your name on this page or any other page you submit with your work. Instead use the student

More information

Nonlinear Control Systems

Nonlinear Control Systems Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 3. Fundamental properties IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Example Consider the system ẋ = f

More information

Time Series Analysis. Solutions to problems in Chapter 7 IMM

Time Series Analysis. Solutions to problems in Chapter 7 IMM Time Series Analysis Solutions to problems in Chapter 7 I Solution 7.1 Question 1. As we get by subtraction kx t = 1 + B + B +...B k 1 )ǫ t θb) = 1 + B +... + B k 1 ), and BθB) = B + B +... + B k ), 1

More information

BIAS ON THE WIND-SPEED GUST CAUSED BY QUANTIZATION AND DISJUNCT SAMPLING OF THE CUP-ANEMOMETER SIGNAL

BIAS ON THE WIND-SPEED GUST CAUSED BY QUANTIZATION AND DISJUNCT SAMPLING OF THE CUP-ANEMOMETER SIGNAL BIAS ON THE WIND-SPEED GUST CAUSED BY QUANTIZATION AND DISJUNCT SAMPLING OF THE CUP-ANEMOMETER SIGNAL Presentation of Main Results L. Kristensen & Ole Frost Hansen January 16, 7 We are here concerned with

More information

5 Analog carrier modulation with noise

5 Analog carrier modulation with noise 5 Analog carrier modulation with noise 5. Noisy receiver model Assume that the modulated signal x(t) is passed through an additive White Gaussian noise channel. A noisy receiver model is illustrated in

More information

3.1.1 Lightcurve, colour-colour and hardness intensity diagram

3.1.1 Lightcurve, colour-colour and hardness intensity diagram Chapter 3 X ray data analysis methods 3.1 Data Analysis Procedure The analysis and reduction procedure of astronomical data can be broadly classified into two categories - (1) count rate variations as

More information

MATH 23 Exam 2 Review Solutions

MATH 23 Exam 2 Review Solutions MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution

More information

Time-modulation of electron-capture decay factor detected at GSI, Darmstadt

Time-modulation of electron-capture decay factor detected at GSI, Darmstadt Time-modulation of electron-capture decay factor detected at GSI, Darmstadt Byung Kyu Park Department of Physics University of California, Berkeley Physics 250 March 20, 2008 Byung Kyu Park (UC Berkeley)

More information

UC Berkeley Math 10B, Spring 2015: Midterm 2 Prof. Sturmfels, April 9, SOLUTIONS

UC Berkeley Math 10B, Spring 2015: Midterm 2 Prof. Sturmfels, April 9, SOLUTIONS UC Berkeley Math 10B, Spring 2015: Midterm 2 Prof. Sturmfels, April 9, SOLUTIONS 1. (5 points) You are a pollster for the 2016 presidential elections. You ask 0 random people whether they would vote for

More information

Estimation of Parameters

Estimation of Parameters CHAPTER Probability, Statistics, and Reliability for Engineers and Scientists FUNDAMENTALS OF STATISTICAL ANALYSIS Second Edition A. J. Clark School of Engineering Department of Civil and Environmental

More information

Spatial Regression. 3. Review - OLS and 2SLS. Luc Anselin. Copyright 2017 by Luc Anselin, All Rights Reserved

Spatial Regression. 3. Review - OLS and 2SLS. Luc Anselin.   Copyright 2017 by Luc Anselin, All Rights Reserved Spatial Regression 3. Review - OLS and 2SLS Luc Anselin http://spatial.uchicago.edu OLS estimation (recap) non-spatial regression diagnostics endogeneity - IV and 2SLS OLS Estimation (recap) Linear Regression

More information

On Input Design for System Identification

On Input Design for System Identification On Input Design for System Identification Input Design Using Markov Chains CHIARA BRIGHENTI Masters Degree Project Stockholm, Sweden March 2009 XR-EE-RT 2009:002 Abstract When system identification methods

More information

AA242B: MECHANICAL VIBRATIONS

AA242B: MECHANICAL VIBRATIONS AA242B: MECHANICAL VIBRATIONS 1 / 50 AA242B: MECHANICAL VIBRATIONS Undamped Vibrations of n-dof Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations:

More information

Heteroskedasticity. Part VII. Heteroskedasticity

Heteroskedasticity. Part VII. Heteroskedasticity Part VII Heteroskedasticity As of Oct 15, 2015 1 Heteroskedasticity Consequences Heteroskedasticity-robust inference Testing for Heteroskedasticity Weighted Least Squares (WLS) Feasible generalized Least

More information

Es e j4φ +4N n. 16 KE s /N 0. σ 2ˆφ4 1 γ s. p(φ e )= exp 1 ( 2πσ φ b cos N 2 φ e 0

Es e j4φ +4N n. 16 KE s /N 0. σ 2ˆφ4 1 γ s. p(φ e )= exp 1 ( 2πσ φ b cos N 2 φ e 0 Problem 6.15 : he received signal-plus-noise vector at the output of the matched filter may be represented as (see (5-2-63) for example) : r n = E s e j(θn φ) + N n where θ n =0,π/2,π,3π/2 for QPSK, and

More information

Dynamics of structures

Dynamics of structures Dynamics of structures 2.Vibrations: single degree of freedom system Arnaud Deraemaeker (aderaema@ulb.ac.be) 1 One degree of freedom systems in real life 2 1 Reduction of a system to a one dof system Example

More information

LANGEVIN THEORY OF BROWNIAN MOTION. Contents. 1 Langevin theory. 1 Langevin theory 1. 2 The Ornstein-Uhlenbeck process 8

LANGEVIN THEORY OF BROWNIAN MOTION. Contents. 1 Langevin theory. 1 Langevin theory 1. 2 The Ornstein-Uhlenbeck process 8 Contents LANGEVIN THEORY OF BROWNIAN MOTION 1 Langevin theory 1 2 The Ornstein-Uhlenbeck process 8 1 Langevin theory Einstein (as well as Smoluchowski) was well aware that the theory of Brownian motion

More information

4 Classical Coherence Theory

4 Classical Coherence Theory This chapter is based largely on Wolf, Introduction to the theory of coherence and polarization of light [? ]. Until now, we have not been concerned with the nature of the light field itself. Instead,

More information

If we want to analyze experimental or simulated data we might encounter the following tasks:

If we want to analyze experimental or simulated data we might encounter the following tasks: Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction

More information

Regression: Lecture 2

Regression: Lecture 2 Regression: Lecture 2 Niels Richard Hansen April 26, 2012 Contents 1 Linear regression and least squares estimation 1 1.1 Distributional results................................ 3 2 Non-linear effects and

More information

Homework #2 Due Monday, April 18, 2012

Homework #2 Due Monday, April 18, 2012 12.540 Homework #2 Due Monday, April 18, 2012 Matlab solution codes are given in HW02_2012.m This code uses cells and contains the solutions to all the questions. Question 1: Non-linear estimation problem

More information

EE/PEP 345. Modeling and Simulation. Spring Class 11

EE/PEP 345. Modeling and Simulation. Spring Class 11 EE/PEP 345 Modeling and Simulation Class 11 11-1 Output Analysis for a Single Model Performance measures System being simulated Output Output analysis Stochastic character Types of simulations Output analysis

More information

Manufacturing Equipment Control

Manufacturing Equipment Control QUESTION 1 An electric drive spindle has the following parameters: J m = 2 1 3 kg m 2, R a = 8 Ω, K t =.5 N m/a, K v =.5 V/(rad/s), K a = 2, J s = 4 1 2 kg m 2, and K s =.3. Ignore electrical dynamics

More information

System Identification Approach Applied to Drift Estimation.

System 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 information

Phase Synchronization

Phase Synchronization Phase Synchronization Lecture by: Zhibin Guo Notes by: Xiang Fan May 10, 2016 1 Introduction For any mode or fluctuation, we always have where S(x, t) is phase. If a mode amplitude satisfies ϕ k = ϕ k

More information

Econometrics II - EXAM Answer each question in separate sheets in three hours

Econometrics II - EXAM Answer each question in separate sheets in three hours Econometrics II - EXAM Answer each question in separate sheets in three hours. Let u and u be jointly Gaussian and independent of z in all the equations. a Investigate the identification of the following

More information

16.584: Random (Stochastic) Processes

16.584: Random (Stochastic) Processes 1 16.584: Random (Stochastic) Processes X(t): X : RV : Continuous function of the independent variable t (time, space etc.) Random process : Collection of X(t, ζ) : Indexed on another independent variable

More information

that efficiently utilizes the total available channel bandwidth W.

that efficiently utilizes the total available channel bandwidth W. Signal Design for Band-Limited Channels Wireless Information Transmission System Lab. Institute of Communications Engineering g National Sun Yat-sen University Introduction We consider the problem of signal

More information

for valid PSD. PART B (Answer all five units, 5 X 10 = 50 Marks) UNIT I

for valid PSD. PART B (Answer all five units, 5 X 10 = 50 Marks) UNIT I Code: 15A04304 R15 B.Tech II Year I Semester (R15) Regular Examinations November/December 016 PROBABILITY THEY & STOCHASTIC PROCESSES (Electronics and Communication Engineering) Time: 3 hours Max. Marks:

More information

ω 0 = 2π/T 0 is called the fundamental angular frequency and ω 2 = 2ω 0 is called the

ω 0 = 2π/T 0 is called the fundamental angular frequency and ω 2 = 2ω 0 is called the he ime-frequency Concept []. Review of Fourier Series Consider the following set of time functions {3A sin t, A sin t}. We can represent these functions in different ways by plotting the amplitude versus

More information

2A1H Time-Frequency Analysis II Bugs/queries to HT 2011 For hints and answers visit dwm/courses/2tf

2A1H Time-Frequency Analysis II Bugs/queries to HT 2011 For hints and answers visit   dwm/courses/2tf Time-Frequency Analysis II (HT 20) 2AH 2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 20 For hints and answers visit www.robots.ox.ac.uk/ dwm/courses/2tf David Murray. A periodic

More information

First Year Physics: Prelims CP1. Classical Mechanics: Prof. Neville Harnew. Problem Set III : Projectiles, rocket motion and motion in E & B fields

First Year Physics: Prelims CP1. Classical Mechanics: Prof. Neville Harnew. Problem Set III : Projectiles, rocket motion and motion in E & B fields HT017 First Year Physics: Prelims CP1 Classical Mechanics: Prof Neville Harnew Problem Set III : Projectiles, rocket motion and motion in E & B fields Questions 1-10 are standard examples Questions 11-1

More information

Reliability Theory of Dynamic Loaded Structures (cont.) Calculation of Out-Crossing Frequencies Approximations to the Failure Probability.

Reliability Theory of Dynamic Loaded Structures (cont.) Calculation of Out-Crossing Frequencies Approximations to the Failure Probability. Outline of Reliability Theory of Dynamic Loaded Structures (cont.) Calculation of Out-Crossing Frequencies Approximations to the Failure Probability. Poisson Approximation. Upper Bound Solution. Approximation

More information

ECE 636: Systems identification

ECE 636: Systems identification ECE 636: Systems identification Lectures 9 0 Linear regression Coherence Φ ( ) xy ω γ xy ( ω) = 0 γ Φ ( ω) Φ xy ( ω) ( ω) xx o noise in the input, uncorrelated output noise Φ zz Φ ( ω) = Φ xy xx ( ω )

More information

13. Power Spectrum. For a deterministic signal x(t), the spectrum is well defined: If represents its Fourier transform, i.e., if.

13. Power Spectrum. For a deterministic signal x(t), the spectrum is well defined: If represents its Fourier transform, i.e., if. For a deterministic signal x(t), the spectrum is well defined: If represents its Fourier transform, i.e., if jt X ( ) = xte ( ) dt, (3-) then X ( ) represents its energy spectrum. his follows from Parseval

More information

EL1820 Modeling of Dynamical Systems

EL1820 Modeling of Dynamical Systems EL1820 Modeling of Dynamical Systems Lecture 9 - Parameter estimation in linear models Model structures Parameter estimation via prediction error minimization Properties of the estimate: bias and variance

More information

On dealing with spatially correlated residuals in remote sensing and GIS

On dealing with spatially correlated residuals in remote sensing and GIS On dealing with spatially correlated residuals in remote sensing and GIS Nicholas A. S. Hamm 1, Peter M. Atkinson and Edward J. Milton 3 School of Geography University of Southampton Southampton SO17 3AT

More information

Square Root Raised Cosine Filter

Square Root Raised Cosine Filter Wireless Information Transmission System Lab. Square Root Raised Cosine Filter Institute of Communications Engineering National Sun Yat-sen University Introduction We consider the problem of signal design

More information