VIII. Coherence and Transfer Function Applications A. Coherence Function Estimates
|
|
- Drusilla Shields
- 6 years ago
- Views:
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 Wb3: SIPE lecture Correlation functions in time & frequency domain Alfred C. Schouten, Dept. of Biomechanical Engineering (BMechE), Fac. 3mE // Delft University
More informationTime 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 informationBasic 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 informationELEG 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 informationX. 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 informationRandomly 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 informationStochastic 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 informationModern 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 informationReview 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 informationSIO 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 informationStatistical 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 informationSinusoidal 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 informationEconometrics 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 informationECE 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 informationStochastic 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 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 information13.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 informationEconometrics 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 informationAn 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 informationINDIAN 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 information2. (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 informationLecture 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 informationNon-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 informationFourier 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 informationModule 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 informationScattering 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 informationWind 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 informationResampling 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 informationChapter 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 informationTSKS01 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 informationThe 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 informationMatrix 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 information6.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 informationLecture 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 informationThe 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 informationComputer 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 informationNONLINEAR 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 informationENSC327 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 informationGATE 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 informationSignal 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 informationChapter 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 informationCommunication 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 informationModule 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 informationVirtual 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 informationGMM, 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 informationLecture 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 informationProf. 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 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 informationSTATISTICAL 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 informationAnalytical 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 informationShort-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 informationEAS 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 informationCalculating 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 informationMATH 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 informationSeparation 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 informationSignal 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 informationHST 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 informationSystem 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 informationIntroduction 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 informationENSC327 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 information2A1H 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 informationA 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 informationElectrical 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 informationNonlinear 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 informationTime 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 informationBIAS 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 information5 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 information3.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 informationMATH 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 informationTime-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 informationUC 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 informationEstimation 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 informationSpatial 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 informationOn 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 informationAA242B: 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 informationHeteroskedasticity. 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 informationEs 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 informationDynamics 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 informationLANGEVIN 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 information4 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 informationIf 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 informationRegression: 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 informationHomework #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 informationEE/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 informationManufacturing 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 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 informationPhase 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 informationEconometrics 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 information16.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 informationthat 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 informationfor 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
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 information2A1H 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 informationFirst 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 informationReliability 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 informationECE 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 information13. 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 informationEL1820 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 informationOn 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 informationSquare 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