Systems Approaches to Estimation Problems in Thin Film Processing

Size: px
Start display at page:

Download "Systems Approaches to Estimation Problems in Thin Film Processing"

Transcription

1 Systems Approaches to Estimation Problems in Thin Film Processing Tyrone Vincent October 20, 2008 T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

2 Acknowledgements Matt Hilt, Numerica Corporation Renee Spinhirne, Lockheed-Martin Bharat Joshi, University of North Carolina, Charlotte Lin Simpson, NREL Kameshwar Poolla, UCB Costas Spanos, UCB T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

3 The Problem of Interest State Estimation of Dynamic Systems x 3 u x 2 y x 1 goal: determine current internal state from observable measurements T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

4 The Problem of Interest State Estimation of Dynamic Systems x 3 u(t) y(t) u t x 2 t y x 1 goal: determine current internal state from observable measurements T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

5 The Problem of Interest State Estimation of Dynamic Systems x 3 u(t) y(t) u t x 2 t y x 1 goal: determine current internal state from observable measurements T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

6 The Problem of Interest State Estimation of Dynamic Systems x 3 u(t) y(t) u t x 2 t y x 1 goal: determine current internal state from observable measurements T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

7 The Problem of Interest State Estimation of Dynamic Systems x 3 u(t) y(t) u t x 2 t y x 1 goal: determine current internal state from observable measurements T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

8 The Problem of Interest State Estimation of Dynamic Systems x 3 u(t) y(t) u t x 2 t y x 1 goal: determine current internal state from observable measurements T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

9 The Problem of Interest Examples given input/output data, estimate processing conditions given in-situ reflectance measurement, estimate deposition rate given partial measurement of wafer parameter, estimate parameter value across the wafer Systems Questions Experiment Design Sensor Selection Measurement Sequencing T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

10 The Problem of Interest Examples given input/output data, estimate processing conditions given in-situ reflectance measurement, estimate deposition rate given partial measurement of wafer parameter, estimate parameter value across the wafer Systems Questions Experiment Design Sensor Selection Measurement Sequencing T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

11 Outline 1 Enhancement of Estimation of Relative Deposition Rate 2 Observability Enhancement for Optical Measurements 3 Measurement Sequencing in Semiconductor Manufacturing 4 Conclusion T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

12 Enhancement of Estimation of Relative Deposition Rate Outline 1 Enhancement of Estimation of Relative Deposition Rate 2 Observability Enhancement for Optical Measurements 3 Measurement Sequencing in Semiconductor Manufacturing 4 Conclusion T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

13 Enhancement of Estimation of Relative Deposition Rate Multi-Zone Co-Evaporation Deposition on moving substrate. Effusion sources create metal plume via resistive heating. Different materials can be deposited in each zone Composition measured at exit T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

14 Enhancement of Estimation of Relative Deposition Rate Deposition Model T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

15 Enhancement of Estimation of Relative Deposition Rate Multi-Zone Model with Transport Delay Compensation T set 1 L therm T 1 N T set 2 L therm T 2 N r 1 L plume + r 2 L plume + y With constant set-points, individual deposition rates are not observable. T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

16 Enhancement of Estimation of Relative Deposition Rate Perturbation Experiment idea: add small perturbations to source set-points output perturbations encode relative sensitivity of deposition rate to temperature perturbations. perturbations chosen to meet manufacturing tolerances T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

17 Enhancement of Estimation of Relative Deposition Rate Effusion Source Model effusion rate assuming supersaturation: r e = κ exp(a) ( exp b ) (τ melt ) c+0.5 τ melt deposition rate to good approximation ( r = δη exp b τ melt ) T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

18 Enhancement of Estimation of Relative Deposition Rate Rate Sensitivity biased measurement τ of melt temperature ατ + β = τ melt deposition rate sensitivity to temperature dr dτ = αb τ 2 melt r T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

19 Enhancement of Estimation of Relative Deposition Rate The advantage of ratios ratio of sensitivities dr 1 dτ dr 1 dτ ( τ 2 ) 2 = melt α1 b 1 r 1 τmelt 1 α 2 b 2 r 2 same metal implies b 1 b 2 same type of thermocouple implies α 1 α 2 ratios cancel (most) unknowns dr 1 dτ dr 1 dτ ( τ 2 ) 2 = melt r1 τmelt 1 r 2 melt temperature ratio τ 2 melt τ 1 melt is fairly well known ( ). T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

20 Enhancement of Estimation of Relative Deposition Rate Manufacturing Tolerances During processing, perturbations are undesired in Total amount deposited Cu/(In+Ga) (controls microstructure) Ga/(In+Ga) (controls bandgap) T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

21 Enhancement of Estimation of Relative Deposition Rate Example Experiment Data from production chamber courtesy Global Solar Energy, Tucson, AZ Independent perturbation applied to sources in two zones depositing the same material. Experiment increased output standard deviation from 20 to 40 Angstroms. Temperature Change (K) Relative Thickness (A) Perturbation Sequences Composition Variation Zone 1 Zone Time (s) T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

22 Enhancement of Estimation of Relative Deposition Rate Experiment Analysis T set 1 L therm T 1 dr 1 dτ T set 2 L therm T 2 dr 2 dτ r 1 L plume + r 2 L plume + y Find gains to best match small scale input/output behavior T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

23 Enhancement of Estimation of Relative Deposition Rate Relative Deposition Rate Estimate filtered output model fit 100 Ratio estimated from shutter experiment: 3.5 Ratio estimated from system identification experiment: 2.4 Magnitude Time (s) T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

24 Enhancement of Estimation of Relative Deposition Rate Block diagram in terms of input variance assumption: input has characteristics of i.i.d random sequence with covariance R. T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

25 Enhancement of Estimation of Relative Deposition Rate Estimate Quality Analysis s N( s, S 1 ) g = s 1 s s s 1 3 p g (g) g T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

26 Enhancement of Estimation of Relative Deposition Rate Estimate Quality Analysis s N( s, S 1 ) g = s 1 s s s 1 3 p g (g) g T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

27 Enhancement of Estimation of Relative Deposition Rate Estimate Quality Result T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

28 Enhancement of Estimation of Relative Deposition Rate Convex approximations for ratio distributions E(S) = { s (s s) S(s s) δ 2} s s 1 3 p g (g) g T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

29 Enhancement of Estimation of Relative Deposition Rate Convex approximations for ratio distributions 2 φ max (S) = { max subject to ratio(s) s E(S) 1.5 s min s s 1 3 p g (g) 2 1 φ min g T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

30 Enhancement of Estimation of Relative Deposition Rate Convex approximations for ratio distributions 2 φ min (S) = { min subject to ratio(s) s E(S) 1.5 s min s s max s 1 3 p g (g) 2 1 φ min φ max g T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

31 Enhancement of Estimation of Relative Deposition Rate Convex approximations for ratio distributions φ(s) = φ max (S) φ min (S) s min s s max s 1 3 p g (g) 2 1 φ min φ max g T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

32 Enhancement of Estimation of Relative Deposition Rate Convex approximations for ratio distributions ρ = s s min 0.5 s max s 1 3 p g (g) 2 1 φ min φ max g T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

33 Enhancement of Estimation of Relative Deposition Rate Properties of ratio metrics Input and Noise Covariances Estimate Covariance φ for ratio estimate is a convex function of the parameters of the input covariance Input Covariance and Disturbance Covariance Performance Variable Covariance φ for ratio of performance variable is a quasi-linear function of the parameters of the input covariance. optimal input covariance found by solving convex optimization problem. T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

34 Enhancement of Estimation of Relative Deposition Rate Example Two zone/two metal process goal: optimize ratio estimate for indium T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

35 Enhancement of Estimation of Relative Deposition Rate Example Case 1: Indium perturbations only, no correlation R = T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

36 Enhancement of Estimation of Relative Deposition Rate Example Case 2: Indium perturbations only R = T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

37 Enhancement of Estimation of Relative Deposition Rate Example Case 2: Indium and Gallium perturbations R = T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

38 Outline Observability Enhancement for Optical Measurements 1 Enhancement of Estimation of Relative Deposition Rate 2 Observability Enhancement for Optical Measurements 3 Measurement Sequencing in Semiconductor Manufacturing 4 Conclusion T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

39 n the next section, we describe how the optical data is currently used to obtain lm Observability Enhancement for Optical Measurements ness. in-situ Reflectometry Current Practice any optical measurements of lm thickness used in the semiconductor industry depend in lm interference eects. This behavior is depicted in Figure 3.4. In this example, ht wave traveling through air strikes a wafer made up of a stack of parallel thin lms. incident light wave is partially reected and partially transmitted at the boundary monitor of wafer surface during deposition via reflectance een two media with dierent refractive indices ni. Reections from lower in the stack ine and interfere with the primary reection o the top of the stack. The result is a ted light wave whose intensity varies as a function of the incident light wavelength, the of incidence, media refractive indices, the polarization state, and most importantly, hicknesses of the lms in the stack Maximum singular value Time (s) Output Time (s) (a) Single Wavelength Observability Max and Min Singular Values Minimum singular value Maximum singular value Time (s) Output Figure 3.4: Thin lm interference he following is an abstraction of the thickness measurement problem which contains the Time (s) r of many of the measurement techniques used today. A static nonlinear model for the T. Vincent (IMPACT) Systems Approaches to Estimation (b) Dual Wavelength October 20, / 45

40 Observability Enhancement for Optical Measurements System block diagram Disturbance Power Flow Throttle PECVD Thickness Composition r(t, C) + Voltage Gain Noise ẋ = f(x, u) + v y = g(x, λ) + n T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

41 Observability Enhancement for Optical Measurements Parameter extraction techniques Least Squares Moving Horizon Filter Extended Kalman Filter 12 Output Time (s) T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

42 Observability Enhancement for Optical Measurements Parameter extraction techniques Least Squares Moving Horizon Filter Extended Kalman Filter 12 Output Time (s) T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

43 Observability Enhancement for Optical Measurements Parameter extraction techniques Least Squares Moving Horizon Filter Extended Kalman Filter 12 Output Time (s) T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

44 Observability Enhancement for Optical Measurements Observability 12 thickness, deposition rate, gains Output Time (s) observability rank of output mapping T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

45 Observability Enhancement for Optical Measurements Wavelength Selection Time (s) (a) Single Wavelength Observability Max and Min Singular Values 10 8 Minimum singular value Observability Max and Min Singular Values Minimum singular value Maximum singular value Time (s) Output 10 1 Maximum singular value Time (s) Output Time (s) (a) Single Wavelength Time (s) (b) Dual Wavelength Observability Max and Min Singular Values T. Vincent (IMPACT) Systems Approaches 3 to Estimation October 20, / 45 Figure 3.12: Maximum and minimum singular values of obse

46 Measurement Sequencing in Semiconductor Manufacturing Outline 1 Enhancement of Estimation of Relative Deposition Rate 2 Observability Enhancement for Optical Measurements 3 Measurement Sequencing in Semiconductor Manufacturing 4 Conclusion T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

47 Measurement Sequencing in Semiconductor Manufacturing Starting point: Hierarchical Spatial Variability Model Jason P. Cain and Costas J. Spanos, Electrical linewidth metrology for systematic CD variation characterization and causal analysis Kun Qian and Costas J. Spanos, A Comprehensive Model of Process Variability for Statistical Timing Optimization L gate (f, d, k) = L 0 + L a (k) + L b (f, d, k) + L c (f, k) +L d (d, k) + L e (d) + L f (f, d, k) L gate (f) (e) (d) (c) (b) (a) pos T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

48 Measurement Sequencing in Semiconductor Manufacturing Terms in the Model L 0 - nominal value L a (k) - random variable N(0, σ a ) L b (f, d, k) = φ 1 (f, d)θ 1 (k) - second order polynomial L c (f, k) - i.i.d random sequence N(0, σ c ) L d (f, k) = φ 2 (f)θ 2 (k) - second order polynomial L e (d) - known layout dependent function L f (f, d, k) - i.i.d random sequence N(0, σ f ) T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

49 Measurement Sequencing in Semiconductor Manufacturing Are Lengths Correlated? Want to know: L gate (f, d, k) for all f, d, k. Does measurement of L gate (f 1, d 1, k) help me predict L gate (f 2, d 2, k)? Definition The covariance of random variables x 1 and x 2 is Result cov (x 1, x 2 ) = E [(x 1 E(x 1 ))(x 2 E(x 2 ))] If x 1, x 2 are Gaussian, and cov (x 1, x 2 ) = 0, p x1 (x 1 x 2 ) = p x1 (x 1 ). Measurement of x 1 does not help prediction of x 2 when cov (x 1, x 2 ) = 0. T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

50 Measurement Sequencing in Semiconductor Manufacturing Are Lengths Correlated? Want to know: L gate (f, d, k) for all f, d, k. Does measurement of L gate (f 1, d 1, k) help me predict L gate (f 2, d 2, k)? Definition The covariance of random variables x 1 and x 2 is Result cov (x 1, x 2 ) = E [(x 1 E(x 1 ))(x 2 E(x 2 ))] If x 1, x 2 are Gaussian, and cov (x 1, x 2 ) = 0, p x1 (x 1 x 2 ) = p x1 (x 1 ). Measurement of x 1 does not help prediction of x 2 when cov (x 1, x 2 ) = 0. T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

51 Measurement Sequencing in Semiconductor Manufacturing Are Lengths Correlated? Want to know: L gate (f, d, k) for all f, d, k. Does measurement of L gate (f 1, d 1, k) help me predict L gate (f 2, d 2, k)? Definition The covariance of random variables x 1 and x 2 is Result cov (x 1, x 2 ) = E [(x 1 E(x 1 ))(x 2 E(x 2 ))] If x 1, x 2 are Gaussian, and cov (x 1, x 2 ) = 0, p x1 (x 1 x 2 ) = p x1 (x 1 ). Measurement of x 1 does not help prediction of x 2 when cov (x 1, x 2 ) = 0. T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

52 Measurement Sequencing in Semiconductor Manufacturing Covariance Let L 1 = L gate (f 1, d 1, k), L 2 = L gate (f 2, d 2, k) cov (L 1, L 2 ) = var (L a ) (+ var (L c )) + φ 1 (f 1, d 1 )φ 1 (f 2, d 2 ) var (θ 1 ) + φ 2 (f 1 )φ 2 (f 2 ) var (θ 2 ) Prediction possible if L a, θ 1, and θ 2 vary with k. T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

53 Measurement Sequencing in Semiconductor Manufacturing Illustration of Reduction in Uncertainty L gate (f) (e) (d) (c) (b) (a) T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45 pos

54 Measurement Sequencing in Semiconductor Manufacturing Illustration of Reduction in Uncertainty L gate (f) (e) (d) (c) (b) (a) T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45 pos

55 Measurement Sequencing in Semiconductor Manufacturing Illustration of Reduction in Uncertainty L gate (f) (e) (d) (c) (b) (a) T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45 pos

56 Measurement Sequencing in Semiconductor Manufacturing Dynamic Model What is the variation from wafer to wafer? L gate (b) (a) pos T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

57 Measurement Sequencing in Semiconductor Manufacturing Dynamic Model What is the variation from wafer to wafer? L gate (b) (a) pos T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

58 Measurement Sequencing in Semiconductor Manufacturing Dynamic Model What is the variation from wafer to wafer? L gate (b) (a) pos T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

59 Measurement Sequencing in Semiconductor Manufacturing Dynamic Model What is the variation from wafer to wafer? L gate (b) (a) pos T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

60 Measurement Sequencing in Semiconductor Manufacturing Time Correlation R theta (τ) θ k+1 = Aθ k + w k τ T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

61 Measurement Sequencing in Semiconductor Manufacturing Choosing a Sampling Sequence Dynamic Model, x = [ θ1 θ 2 ]. x k+1 = Ax k + w k y k = C(f, d) k x k + e k Uncertainty Prediction - Kalman Filter T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

62 Measurement Sequencing in Semiconductor Manufacturing Sampling Sequence Selection Uncertainty for wafer k after measurement L gate (f) (e) (d) (c) (b) (a) pos T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

63 Measurement Sequencing in Semiconductor Manufacturing Sampling Sequence Selection Uncertainty for wafer k + 1 before measurement L gate (f) (e) (d) (c) (b) (a) pos T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

64 Outline Conclusion 1 Enhancement of Estimation of Relative Deposition Rate 2 Observability Enhancement for Optical Measurements 3 Measurement Sequencing in Semiconductor Manufacturing 4 Conclusion T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

65 Conclusion Conclusion Three illustrations of analysis of measurement accuracy Tools: uncertainty modeling and propagation optimization linking measurements over time using dynamic models Dynamic modeling of wafer processing may lead to more efficient measurement sequencing Supported in part by NSF grant ECS and DOE Programs DE-ZDO and DE-FG36-08GO T. Vincent (IMPACT) Systems Approaches to Estimation October 20, / 45

Subject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)

Subject: Optimal Control Assignment-1 (Related to Lecture notes 1-10) Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must

More information

ik () uk () Today s menu Last lecture Some definitions Repeatability of sensing elements

ik () uk () Today s menu Last lecture Some definitions Repeatability of sensing elements Last lecture Overview of the elements of measurement systems. Sensing elements. Signal conditioning elements. Signal processing elements. Data presentation elements. Static characteristics of measurement

More information

Sensors and Metrology. Outline

Sensors and Metrology. Outline Sensors and Metrology A Survey 1 Outline General Issues & the SIA Roadmap Post-Process Sensing (SEM/AFM, placement) In-Process (or potential in-process) Sensors temperature (pyrometry, thermocouples, acoustic

More information

VARIANCE COMPUTATION OF MODAL PARAMETER ES- TIMATES FROM UPC SUBSPACE IDENTIFICATION

VARIANCE COMPUTATION OF MODAL PARAMETER ES- TIMATES FROM UPC SUBSPACE IDENTIFICATION VARIANCE COMPUTATION OF MODAL PARAMETER ES- TIMATES FROM UPC SUBSPACE IDENTIFICATION Michael Döhler 1, Palle Andersen 2, Laurent Mevel 1 1 Inria/IFSTTAR, I4S, Rennes, France, {michaeldoehler, laurentmevel}@inriafr

More information

Lecture 2: Thin Films. Thin Films. Calculating Thin Film Stack Properties. Jones Matrices for Thin Film Stacks. Mueller Matrices for Thin Film Stacks

Lecture 2: Thin Films. Thin Films. Calculating Thin Film Stack Properties. Jones Matrices for Thin Film Stacks. Mueller Matrices for Thin Film Stacks Lecture 2: Thin Films Outline Thin Films 2 Calculating Thin Film Stack Properties 3 Jones Matrices for Thin Film Stacks 4 Mueller Matrices for Thin Film Stacks 5 Mueller Matrix for Dielectrica 6 Mueller

More information

Swing Curves. T h e L i t h o g r a p h y T u t o r (Summer 1994) Chris A. Mack, FINLE Technologies, Austin, Texas

Swing Curves. T h e L i t h o g r a p h y T u t o r (Summer 1994) Chris A. Mack, FINLE Technologies, Austin, Texas T h e L i t h o g r a p h y T u t o r (Summer 1994) Swing Curves Chris A. Mack, FINLE Technologies, Austin, Texas In the last edition of this column, we saw that exposing a photoresist involves the propagation

More information

Calculating Thin Film Stack Properties. Polarization Properties of Thin Films

Calculating Thin Film Stack Properties. Polarization Properties of Thin Films Lecture 6: Thin Films Outline 1 Thin Films 2 Calculating Thin Film Stack Properties 3 Polarization Properties of Thin Films 4 Anti-Reflection Coatings 5 Interference Filters Christoph U. Keller, Utrecht

More information

Kalman Filter. Predict: Update: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q

Kalman Filter. Predict: Update: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q Kalman Filter Kalman Filter Predict: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q Update: K = P k k 1 Hk T (H k P k k 1 Hk T + R) 1 x k k = x k k 1 + K(z k H k x k k 1 ) P k k =(I

More information

GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL OF ELECTRICAL AND COMPUTER ENGINEERING Final Examination - Fall 2015 EE 4601: Communication Systems

GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL OF ELECTRICAL AND COMPUTER ENGINEERING Final Examination - Fall 2015 EE 4601: Communication Systems GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL OF ELECTRICAL AND COMPUTER ENGINEERING Final Examination - Fall 2015 EE 4601: Communication Systems Aids Allowed: 2 8 1/2 X11 crib sheets, calculator DATE: Tuesday

More information

Calculating Thin Film Stack Properties

Calculating Thin Film Stack Properties Lecture 5: Thin Films Outline 1 Thin Films 2 Calculating Thin Film Stack Properties 3 Fabry-Perot Tunable Filter 4 Anti-Reflection Coatings 5 Interference Filters Christoph U. Keller, Leiden University,

More information

Stochastic Models, Estimation and Control Peter S. Maybeck Volumes 1, 2 & 3 Tables of Contents

Stochastic Models, Estimation and Control Peter S. Maybeck Volumes 1, 2 & 3 Tables of Contents Navtech Part #s Volume 1 #1277 Volume 2 #1278 Volume 3 #1279 3 Volume Set #1280 Stochastic Models, Estimation and Control Peter S. Maybeck Volumes 1, 2 & 3 Tables of Contents Volume 1 Preface Contents

More information

Sensors and Metrology

Sensors and Metrology Sensors and Metrology A Survey 1 Outline General Issues & the SIA Roadmap Post-Process Sensing (SEM/AFM, placement) In-Process (or potential in-process) Sensors temperature (pyrometry, thermocouples, acoustic

More information

Robot Localization and Kalman Filters

Robot Localization and Kalman Filters Robot Localization and Kalman Filters Rudy Negenborn rudy@negenborn.net August 26, 2003 Outline Robot Localization Probabilistic Localization Kalman Filters Kalman Localization Kalman Localization with

More information

COS Lecture 16 Autonomous Robot Navigation

COS Lecture 16 Autonomous Robot Navigation COS 495 - Lecture 16 Autonomous Robot Navigation Instructor: Chris Clark Semester: Fall 011 1 Figures courtesy of Siegwart & Nourbakhsh Control Structure Prior Knowledge Operator Commands Localization

More information

Chapter 3. Design for manufacture of Fabry-Perot cavity sensors

Chapter 3. Design for manufacture of Fabry-Perot cavity sensors Chapter 3 Design for manufacture of Fabry-Perot cavity sensors When Fabry-Perot cavity sensors are manufactured, the thickness of each layer must be tightly controlled to achieve the target performance

More information

Electrodynamics Qualifier Examination

Electrodynamics Qualifier Examination Electrodynamics Qualifier Examination August 15, 2007 General Instructions: In all cases, be sure to state your system of units. Show all your work, write only on one side of the designated paper, and

More information

Observability and state estimation

Observability and state estimation EE263 Autumn 2015 S Boyd and S Lall Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time observability

More information

Lecture 19 Observability and state estimation

Lecture 19 Observability and state estimation EE263 Autumn 2007-08 Stephen Boyd Lecture 19 Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time

More information

Convergence of Square Root Ensemble Kalman Filters in the Large Ensemble Limit

Convergence of Square Root Ensemble Kalman Filters in the Large Ensemble Limit Convergence of Square Root Ensemble Kalman Filters in the Large Ensemble Limit Evan Kwiatkowski, Jan Mandel University of Colorado Denver December 11, 2014 OUTLINE 2 Data Assimilation Bayesian Estimation

More information

Nonlinear Control Systems

Nonlinear Control Systems Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 5. Input-Output Stability DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Input-Output Stability y = Hu H denotes

More information

Across-wafer CD Uniformity Enhancement through Control of Multi-zone PEB Profiles

Across-wafer CD Uniformity Enhancement through Control of Multi-zone PEB Profiles Across-wafer CD Uniformity Enhancement through Control of Multi-zone PEB Profiles Qiaolin Zhang *a, Paul Friedberg b, Cherry Tang c Bhanwar Singh c, Kameshwar Poolla a, Costas J. Spanos b a Dept of Mechanical

More information

3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE

3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE 3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE 3.0 INTRODUCTION The purpose of this chapter is to introduce estimators shortly. More elaborated courses on System Identification, which are given

More information

Probabilistic Graphical Models

Probabilistic Graphical Models Probabilistic Graphical Models Brown University CSCI 2950-P, Spring 2013 Prof. Erik Sudderth Lecture 12: Gaussian Belief Propagation, State Space Models and Kalman Filters Guest Kalman Filter Lecture by

More information

1 Kalman Filter Introduction

1 Kalman Filter Introduction 1 Kalman Filter Introduction You should first read Chapter 1 of Stochastic models, estimation, and control: Volume 1 by Peter S. Maybec (available here). 1.1 Explanation of Equations (1-3) and (1-4) Equation

More information

Question 1. (Marks 16)

Question 1. (Marks 16) 5 Question 1. (Marks 16) Consider the circuit shown in the figure, where C 1 = 6.00µF, C 2 = 3.00µF, and V = 20.0V. Capacitor C 1 is first charged by closing switch S 1. Switch S 1 is then opened, and

More information

ELE539A: Optimization of Communication Systems Lecture 15: Semidefinite Programming, Detection and Estimation Applications

ELE539A: Optimization of Communication Systems Lecture 15: Semidefinite Programming, Detection and Estimation Applications ELE539A: Optimization of Communication Systems Lecture 15: Semidefinite Programming, Detection and Estimation Applications Professor M. Chiang Electrical Engineering Department, Princeton University March

More information

E190Q Lecture 11 Autonomous Robot Navigation

E190Q Lecture 11 Autonomous Robot Navigation E190Q Lecture 11 Autonomous Robot Navigation Instructor: Chris Clark Semester: Spring 013 1 Figures courtesy of Siegwart & Nourbakhsh Control Structures Planning Based Control Prior Knowledge Operator

More information

Overview of Beamforming

Overview of Beamforming Overview of Beamforming Arye Nehorai Preston M. Green Department of Electrical and Systems Engineering Washington University in St. Louis March 14, 2012 CSSIP Lab 1 Outline Introduction Spatial and temporal

More information

CDS 101/110a: Lecture 10-1 Robust Performance

CDS 101/110a: Lecture 10-1 Robust Performance CDS 11/11a: Lecture 1-1 Robust Performance Richard M. Murray 1 December 28 Goals: Describe how to represent uncertainty in process dynamics Describe how to analyze a system in the presence of uncertainty

More information

Lecture Note 12: Kalman Filter

Lecture Note 12: Kalman Filter ECE 645: Estimation Theory Spring 2015 Instructor: Prof. Stanley H. Chan Lecture Note 12: Kalman Filter LaTeX prepared by Stylianos Chatzidakis) May 4, 2015 This lecture note is based on ECE 645Spring

More information

Phys 531 Lecture 27 6 December 2005

Phys 531 Lecture 27 6 December 2005 Phys 531 Lecture 27 6 December 2005 Final Review Last time: introduction to quantum field theory Like QM, but field is quantum variable rather than x, p for particle Understand photons, noise, weird quantum

More information

Aircraft Turbofan Engine Health Estimation Using Constrained Kalman Filtering. Funded by the NASA Aviation Safety Program June 16, 2003

Aircraft Turbofan Engine Health Estimation Using Constrained Kalman Filtering. Funded by the NASA Aviation Safety Program June 16, 2003 Aircraft Turbofan Engine Health Estimation Using Constrained Kalman Filtering Dan Simon Electrical Engineering Dept. Cleveland State University Cleveland, Ohio Donald L. Simon US Army Research Laboratory

More information

State Estimation and Motion Tracking for Spatially Diverse VLC Networks

State Estimation and Motion Tracking for Spatially Diverse VLC Networks State Estimation and Motion Tracking for Spatially Diverse VLC Networks GLOBECOM Optical Wireless Communications Workshop December 3, 2012 Anaheim, CA Michael Rahaim mrahaim@bu.edu Gregary Prince gbprince@bu.edu

More information

Sensor Tasking and Control

Sensor Tasking and Control Sensor Tasking and Control Sensing Networking Leonidas Guibas Stanford University Computation CS428 Sensor systems are about sensing, after all... System State Continuous and Discrete Variables The quantities

More information

Multiuser Capacity Analysis of WDM in Nonlinear Fiber Optics

Multiuser Capacity Analysis of WDM in Nonlinear Fiber Optics Multiuser Capacity Analysis of WDM in Nonlinear Fiber Optics Mohammad H. Taghavi N., George C. Papen, and Paul H. Siegel Dept. of ECE, UCSD, La Jolla, CA 92093 Email: {mtaghavi, gpapen, psiegel}@ucsd.edu

More information

Dynamic measurement: application of system identification in metrology

Dynamic measurement: application of system identification in metrology 1 / 25 Dynamic measurement: application of system identification in metrology Ivan Markovsky Dynamic measurement takes into account the dynamical properties of the sensor 2 / 25 model of sensor as dynamical

More information

4. DATA ASSIMILATION FUNDAMENTALS

4. DATA ASSIMILATION FUNDAMENTALS 4. DATA ASSIMILATION FUNDAMENTALS... [the atmosphere] "is a chaotic system in which errors introduced into the system can grow with time... As a consequence, data assimilation is a struggle between chaotic

More information

Lecture 5: Polarization. Polarized Light in the Universe. Descriptions of Polarized Light. Polarizers. Retarders. Outline

Lecture 5: Polarization. Polarized Light in the Universe. Descriptions of Polarized Light. Polarizers. Retarders. Outline Lecture 5: Polarization Outline 1 Polarized Light in the Universe 2 Descriptions of Polarized Light 3 Polarizers 4 Retarders Christoph U. Keller, Leiden University, keller@strw.leidenuniv.nl ATI 2016,

More information

VASE. J.A. Woollam Co., Inc. Ellipsometry Solutions

VASE. J.A. Woollam Co., Inc. Ellipsometry Solutions VASE J.A. Woollam Co., Inc. Ellipsometry Solutions Accurate Capabilities The VASE is our most accurate and versatile ellipsometer for research on all types of materials: semiconductors, dielectrics, polymers,

More information

Nonlinear Identification of Backlash in Robot Transmissions

Nonlinear Identification of Backlash in Robot Transmissions Nonlinear Identification of Backlash in Robot Transmissions G. Hovland, S. Hanssen, S. Moberg, T. Brogårdh, S. Gunnarsson, M. Isaksson ABB Corporate Research, Control Systems Group, Switzerland ABB Automation

More information

Development and Applications of Acoustic Charge Transport

Development and Applications of Acoustic Charge Transport Development and Applications of Acoustic Charge Transport CID #478 Physics 222 Term Paper #1 Abstract Acoustic charge transport is accomplished through the use of a surface acoustic wave passing through

More information

Introduction to Reflectometry and Small Angle Scattering under Grazing Incidence

Introduction to Reflectometry and Small Angle Scattering under Grazing Incidence Mitglied der Helmholtz-Gemeinschaft Introduction to Reflectometry and Small Angle Scattering under Grazing Incidence E. Kentzinger Jülich Center for Neutron Science and Peter Grünberg Institut Jülich Research

More information

Propagation losses in optical fibers

Propagation losses in optical fibers Chapter Dielectric Waveguides and Optical Fibers 1-Fev-017 Propagation losses in optical fibers Charles Kao, Nobel Laureate (009) Courtesy of the Chinese University of Hong Kong S.O. Kasap, Optoelectronics

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

Autonomous Mobile Robot Design

Autonomous Mobile Robot Design Autonomous Mobile Robot Design Topic: Extended Kalman Filter Dr. Kostas Alexis (CSE) These slides relied on the lectures from C. Stachniss, J. Sturm and the book Probabilistic Robotics from Thurn et al.

More information

POLARIZATION OF LIGHT

POLARIZATION OF LIGHT POLARIZATION OF LIGHT OVERALL GOALS The Polarization of Light lab strongly emphasizes connecting mathematical formalism with measurable results. It is not your job to understand every aspect of the theory,

More information

FEEDBACK CONTROL OF GROWTH RATE AND SURFACE ROUGHNESS IN THIN FILM GROWTH. Yiming Lou and Panagiotis D. Christofides

FEEDBACK CONTROL OF GROWTH RATE AND SURFACE ROUGHNESS IN THIN FILM GROWTH. Yiming Lou and Panagiotis D. Christofides FEEDBACK CONTROL OF GROWTH RATE AND SURFACE ROUGHNESS IN THIN FILM GROWTH Yiming Lou and Panagiotis D. Christofides Department of Chemical Engineering University of California, Los Angeles IEEE 2003 Conference

More information

State Observers and the Kalman filter

State Observers and the Kalman filter Modelling and Control of Dynamic Systems State Observers and the Kalman filter Prof. Oreste S. Bursi University of Trento Page 1 Feedback System State variable feedback system: Control feedback law:u =

More information

Organization. I MCMC discussion. I project talks. I Lecture.

Organization. I MCMC discussion. I project talks. I Lecture. Organization I MCMC discussion I project talks. I Lecture. Content I Uncertainty Propagation Overview I Forward-Backward with an Ensemble I Model Reduction (Intro) Uncertainty Propagation in Causal Systems

More information

Statistics 910, #15 1. Kalman Filter

Statistics 910, #15 1. Kalman Filter Statistics 910, #15 1 Overview 1. Summary of Kalman filter 2. Derivations 3. ARMA likelihoods 4. Recursions for the variance Kalman Filter Summary of Kalman filter Simplifications To make the derivations

More information

A Class of Impedance Tomography based Sensors for Semiconductor Manufacturing

A Class of Impedance Tomography based Sensors for Semiconductor Manufacturing A Class of Impedance Tomography based Sensors for Semiconductor Manufacturing Michiel V.P. Krüger, Kameshwar Poolla, and Costas J. Spanos Department of Mechanical Engineering, Department of Electrical

More information

Lab #3 Transparent Conductors

Lab #3 Transparent Conductors Lab #3 Transparent Conductors R.J. Ellingson and M.J. Heben Sept. 23, 2014 PHYS 4580, 6/7280 Today at 4:30 pm, R1 Room 1010 Impact of Optical Loss in Window Layer in PV Cells O Quantum Efficiency, aka

More information

Improving Yield for High Pin Count Wafer Probing Applications

Improving Yield for High Pin Count Wafer Probing Applications John R. Goulding Martin Elzingre Larry Hendler Technical Marketing Manager 200mm Product Manager Consultant Electroglas, Inc. Electroglas, Inc. Electroglas, Inc. South West Test Workshop June 11-14 th

More information

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

VIII. Coherence and Transfer Function Applications A. Coherence Function Estimates 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

More information

Cross Directional Control

Cross Directional Control Cross Directional Control Graham C. Goodwin Day 4: Lecture 4 16th September 2004 International Summer School Grenoble, France 1. Introduction In this lecture we describe a practical application of receding

More information

Simulations of an Interference Birefringent Thin-Film Filter Used as a Narrow-Band Polarizer

Simulations of an Interference Birefringent Thin-Film Filter Used as a Narrow-Band Polarizer imulations of an Interference Birefringent Thin-Film Filter Used as a Narrow-Band olarizer Hugo uertas de Araújo and ebastião Gomes dos antos Filho University of ão aulo LI / I / EU Av. rof. Luciano Gualberto

More information

Fitting Linear Statistical Models to Data by Least Squares II: Weighted

Fitting Linear Statistical Models to Data by Least Squares II: Weighted Fitting Linear Statistical Models to Data by Least Squares II: Weighted Brian R. Hunt and C. David Levermore University of Maryland, College Park Math 420: Mathematical Modeling April 21, 2014 version

More information

Timing Recovery at Low SNR Cramer-Rao bound, and outperforming the PLL

Timing Recovery at Low SNR Cramer-Rao bound, and outperforming the PLL T F T I G E O R G A I N S T I T U T E O H E O F E A L P R O G R ESS S A N D 1 8 8 5 S E R V L O G Y I C E E C H N O Timing Recovery at Low SNR Cramer-Rao bound, and outperforming the PLL Aravind R. Nayak

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

Dr Ian R. Manchester

Dr Ian R. Manchester Week Content Notes 1 Introduction 2 Frequency Domain Modelling 3 Transient Performance and the s-plane 4 Block Diagrams 5 Feedback System Characteristics Assign 1 Due 6 Root Locus 7 Root Locus 2 Assign

More information

ECE 636: Systems identification

ECE 636: Systems identification ECE 636: Systems identification Lectures 3 4 Random variables/signals (continued) Random/stochastic vectors Random signals and linear systems Random signals in the frequency domain υ ε x S z + y Experimental

More information

ENGR352 Problem Set 02

ENGR352 Problem Set 02 engr352/engr352p02 September 13, 2018) ENGR352 Problem Set 02 Transfer function of an estimator 1. Using Eq. (1.1.4-27) from the text, find the correct value of r ss (the result given in the text is incorrect).

More information

Statistics Part I Introduction. Joe Nahas University of Notre Dame

Statistics Part I Introduction. Joe Nahas University of Notre Dame Statistics Part I Introduction Joe Nahas University of Notre Dame A Very Simple Example: A Pair of Die A pair of six sided die Values for each die: 1, 2, 3, 4, 5, 6. Values for the pair: 2, 3, 4, 5, 6,

More information

Introduction to Dynamic Light Scattering with Applications. Onofrio Annunziata Department of Chemistry Texas Christian University Fort Worth, TX, USA

Introduction to Dynamic Light Scattering with Applications. Onofrio Annunziata Department of Chemistry Texas Christian University Fort Worth, TX, USA Introduction to Dynamic Light Scattering with Applications Onofrio Annunziata Department of Chemistry Texas Christian University Fort Worth, TX, USA Outline Introduction to dynamic light scattering Particle

More information

Uniformly Uniformly-ergodic Markov chains and BSDEs

Uniformly Uniformly-ergodic Markov chains and BSDEs Uniformly Uniformly-ergodic Markov chains and BSDEs Samuel N. Cohen Mathematical Institute, University of Oxford (Based on joint work with Ying Hu, Robert Elliott, Lukas Szpruch) Centre Henri Lebesgue,

More information

LARGE-SCALE TRAFFIC STATE ESTIMATION

LARGE-SCALE TRAFFIC STATE ESTIMATION Hans van Lint, Yufei Yuan & Friso Scholten A localized deterministic Ensemble Kalman Filter LARGE-SCALE TRAFFIC STATE ESTIMATION CONTENTS Intro: need for large-scale traffic state estimation Some Kalman

More information

Optical Measurements of Critical Dimensions at Several Stages of the Mask Fabrication Process

Optical Measurements of Critical Dimensions at Several Stages of the Mask Fabrication Process Optical Measurements of Critical Dimensions at Several Stages of the Mask Fabrication Process John C. Lam, Alexander Gray n&k Technology, Inc., Santa Clara, CA ABSTRACT Critical dimension (CD) metrology

More information

Chapter 2 Process Variability. Overview. 2.1 Sources and Types of Variations

Chapter 2 Process Variability. Overview. 2.1 Sources and Types of Variations Chapter 2 Process Variability Overview Parameter variability has always been an issue in integrated circuits. However, comparing with the size of devices, it is relatively increasing with technology evolution,

More information

1 Name: Student number: DEPARTMENT OF PHYSICS AND PHYSICAL OCEANOGRAPHY MEMORIAL UNIVERSITY OF NEWFOUNDLAND. Fall :00-11:00

1 Name: Student number: DEPARTMENT OF PHYSICS AND PHYSICAL OCEANOGRAPHY MEMORIAL UNIVERSITY OF NEWFOUNDLAND. Fall :00-11:00 1 Name: DEPARTMENT OF PHYSICS AND PHYSICAL OCEANOGRAPHY MEMORIAL UNIVERSITY OF NEWFOUNDLAND Final Exam Physics 3000 December 11, 2012 Fall 2012 9:00-11:00 INSTRUCTIONS: 1. Answer all seven (7) questions.

More information

2.830J / 6.780J / ESD.63J Control of Manufacturing Processes (SMA 6303)

2.830J / 6.780J / ESD.63J Control of Manufacturing Processes (SMA 6303) MIT OpenCourseWare http://ocw.mit.edu 2.830J / 6.780J / ESD.63J Control of Manufacturing Processes (SMA 6303) Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.

More information

Variation-Resistant Dynamic Power Optimization for VLSI Circuits

Variation-Resistant Dynamic Power Optimization for VLSI Circuits Process-Variation Variation-Resistant Dynamic Power Optimization for VLSI Circuits Fei Hu Department of ECE Auburn University, AL 36849 Ph.D. Dissertation Committee: Dr. Vishwani D. Agrawal Dr. Foster

More information

Techniques for directly measuring the absorbance of photoresists at EUV wavelengths

Techniques for directly measuring the absorbance of photoresists at EUV wavelengths Techniques for directly measuring the absorbance of photoresists at EUV wavelengths Manish Chandhok, a Heidi Cao, a Wang Yueh, a Eric Gullikson, b Robert Brainard, c Stewart Robertson c a Intel Corporation,

More information

Expectation. DS GA 1002 Statistical and Mathematical Models. Carlos Fernandez-Granda

Expectation. DS GA 1002 Statistical and Mathematical Models.   Carlos Fernandez-Granda Expectation DS GA 1002 Statistical and Mathematical Models http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall16 Carlos Fernandez-Granda Aim Describe random variables with a few numbers: mean, variance,

More information

Aircraft Turbofan Engine Health Estimation Using Constrained Kalman Filtering

Aircraft Turbofan Engine Health Estimation Using Constrained Kalman Filtering Cleveland State University EngagedScholarship@CSU Electrical Engineering & Computer Science Faculty Publications Electrical Engineering & Computer Science Department 6-16-2003 Aircraft Turbofan Engine

More information

Example 1 linear elastic structure; forces f 1 ;:::;f 100 induce deections d 1 ;:::;d f i F max i, several hundred other constraints: max load p

Example 1 linear elastic structure; forces f 1 ;:::;f 100 induce deections d 1 ;:::;d f i F max i, several hundred other constraints: max load p Convex Optimization in Electrical Engineering Stephen Boyd and Lieven Vandenberghe EE Dept., Stanford ISL EE370 Seminar, 2/9/96 1 Main idea many electrical engineering design problems can be cast as convex

More information

Photolithography II ( Part 1 )

Photolithography II ( Part 1 ) 1 Photolithography II ( Part 1 ) Chapter 14 : Semiconductor Manufacturing Technology by M. Quirk & J. Serda Bjørn-Ove Fimland, Department of Electronics and Telecommunication, Norwegian University of Science

More information

Sequential Monte Carlo Methods (for DSGE Models)

Sequential Monte Carlo Methods (for DSGE Models) Sequential Monte Carlo Methods (for DSGE Models) Frank Schorfheide University of Pennsylvania, PIER, CEPR, and NBER October 23, 2017 Some References These lectures use material from our joint work: Tempered

More information

2-D SENSOR POSITION PERTURBATION ANALYSIS: EQUIVALENCE TO AWGN ON ARRAY OUTPUTS. Volkan Cevher, James H. McClellan

2-D SENSOR POSITION PERTURBATION ANALYSIS: EQUIVALENCE TO AWGN ON ARRAY OUTPUTS. Volkan Cevher, James H. McClellan 2-D SENSOR POSITION PERTURBATION ANALYSIS: EQUIVALENCE TO AWGN ON ARRAY OUTPUTS Volkan Cevher, James H McClellan Georgia Institute of Technology Atlanta, GA 30332-0250 cevher@ieeeorg, jimmcclellan@ecegatechedu

More information

9. Model Selection. statistical models. overview of model selection. information criteria. goodness-of-fit measures

9. Model Selection. statistical models. overview of model selection. information criteria. goodness-of-fit measures FE661 - Statistical Methods for Financial Engineering 9. Model Selection Jitkomut Songsiri statistical models overview of model selection information criteria goodness-of-fit measures 9-1 Statistical models

More 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

Chapter 7. Solar Cell

Chapter 7. Solar Cell Chapter 7 Solar Cell 7.0 Introduction Solar cells are useful for both space and terrestrial application. Solar cells furnish the long duration power supply for satellites. It converts sunlight directly

More information

Time-Varying Parameters

Time-Varying Parameters Kalman Filter and state-space models: time-varying parameter models; models with unobservable variables; basic tool: Kalman filter; implementation is task-specific. y t = x t β t + e t (1) β t = µ + Fβ

More information

The Kalman Filter. Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience. Sarah Dance

The Kalman Filter. Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience. Sarah Dance The Kalman Filter Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience Sarah Dance School of Mathematical and Physical Sciences, University of Reading s.l.dance@reading.ac.uk July

More information

NUMERICAL SIMULATION OF PULSE LASER ABLATION 1

NUMERICAL SIMULATION OF PULSE LASER ABLATION 1 NUMERICAL SIMULATION OF PULSE LASER ABLATION 1 Pritamkumar Dake Mechanical engineering department SRES College of engineering Kopargaon, Maharastra, India Email: 1 Pritam.dake@gmail.com Abstract Pulsed

More information

On the Assessment of a Bayesian Validation Methodology for Data Reduction Models Relevant to Shock Tube Experiments

On the Assessment of a Bayesian Validation Methodology for Data Reduction Models Relevant to Shock Tube Experiments ICES REPORT -36 November 2 On the Assessment of a Bayesian Validation Methodology for Data Reduction Models Relevant to Shock Tube Experiments by M. Panesi, K. Miki, S. Prudhomme, and A. Brandis The Institute

More information

Statistics & Data Sciences: First Year Prelim Exam May 2018

Statistics & Data Sciences: First Year Prelim Exam May 2018 Statistics & Data Sciences: First Year Prelim Exam May 2018 Instructions: 1. Do not turn this page until instructed to do so. 2. Start each new question on a new sheet of paper. 3. This is a closed book

More information

Advanced Adaptive Control for Unintended System Behavior

Advanced Adaptive Control for Unintended System Behavior Advanced Adaptive Control for Unintended System Behavior Dr. Chengyu Cao Mechanical Engineering University of Connecticut ccao@engr.uconn.edu jtang@engr.uconn.edu Outline Part I: Challenges: Unintended

More information

ẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)

ẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7) EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and

More information

Q. 1 Q. 25 carry one mark each.

Q. 1 Q. 25 carry one mark each. GATE 5 SET- ELECTRONICS AND COMMUNICATION ENGINEERING - EC Q. Q. 5 carry one mark each. Q. The bilateral Laplace transform of a function is if a t b f() t = otherwise (A) a b s (B) s e ( a b) s (C) e as

More information

Modelling and System Identification Microexam 1

Modelling and System Identification Microexam 1 Modelling and System Identification Microexam 1 Prof Dr Moritz Diehl, IMTEK, Universität Freiburg, and OPTEC/ESAT-STADIUS, KU Leuven December 4, 013, 16:00-17:00, Freiburg achname: Vorname: Matrikelnummer:

More information

Feedback Linearization based Arc Length Control for Gas Metal Arc Welding

Feedback Linearization based Arc Length Control for Gas Metal Arc Welding 5 American Control Conference June 8-, 5. Portland, OR, USA FrA5.6 Feedback Linearization based Arc Length Control for Gas Metal Arc Welding Jesper S. Thomsen Abstract In this paper a feedback linearization

More information

3 Gramians and Balanced Realizations

3 Gramians and Balanced Realizations 3 Gramians and Balanced Realizations In this lecture, we use an optimization approach to find suitable realizations for truncation and singular perturbation of G. It turns out that the recommended realizations

More information

Modeling and Analysis of Dynamic Systems

Modeling and Analysis of Dynamic Systems Modeling and Analysis of Dynamic Systems by Dr. Guillaume Ducard Fall 2016 Institute for Dynamic Systems and Control ETH Zurich, Switzerland based on script from: Prof. Dr. Lino Guzzella 1/33 Outline 1

More information

5 Transfer function modelling

5 Transfer function modelling MSc Further Time Series Analysis 5 Transfer function modelling 5.1 The model Consider the construction of a model for a time series (Y t ) whose values are influenced by the earlier values of a series

More information

Data assimilation with and without a model

Data assimilation with and without a model Data assimilation with and without a model Tyrus Berry George Mason University NJIT Feb. 28, 2017 Postdoc supported by NSF This work is in collaboration with: Tim Sauer, GMU Franz Hamilton, Postdoc, NCSU

More information

UCSD ECE153 Handout #40 Prof. Young-Han Kim Thursday, May 29, Homework Set #8 Due: Thursday, June 5, 2011

UCSD ECE153 Handout #40 Prof. Young-Han Kim Thursday, May 29, Homework Set #8 Due: Thursday, June 5, 2011 UCSD ECE53 Handout #40 Prof. Young-Han Kim Thursday, May 9, 04 Homework Set #8 Due: Thursday, June 5, 0. Discrete-time Wiener process. Let Z n, n 0 be a discrete time white Gaussian noise (WGN) process,

More information

Estimation of the Capacity of Multipath Infrared Channels

Estimation of the Capacity of Multipath Infrared Channels Estimation of the Capacity of Multipath Infrared Channels Jeffrey B. Carruthers Department of Electrical and Computer Engineering Boston University jbc@bu.edu Sachin Padma Department of Electrical and

More information

Bayes Filter Reminder. Kalman Filter Localization. Properties of Gaussians. Gaussians. Prediction. Correction. σ 2. Univariate. 1 2πσ e.

Bayes Filter Reminder. Kalman Filter Localization. Properties of Gaussians. Gaussians. Prediction. Correction. σ 2. Univariate. 1 2πσ e. Kalman Filter Localization Bayes Filter Reminder Prediction Correction Gaussians p(x) ~ N(µ,σ 2 ) : Properties of Gaussians Univariate p(x) = 1 1 2πσ e 2 (x µ) 2 σ 2 µ Univariate -σ σ Multivariate µ Multivariate

More information

Robot Localisation. Henrik I. Christensen. January 12, 2007

Robot Localisation. Henrik I. Christensen. January 12, 2007 Robot Henrik I. Robotics and Intelligent Machines @ GT College of Computing Georgia Institute of Technology Atlanta, GA hic@cc.gatech.edu January 12, 2007 The Robot Structure Outline 1 2 3 4 Sum of 5 6

More information

Slide a window along the input arc sequence S. Least-squares estimate. σ 2. σ Estimate 1. Statistically test the difference between θ 1 and θ 2

Slide a window along the input arc sequence S. Least-squares estimate. σ 2. σ Estimate 1. Statistically test the difference between θ 1 and θ 2 Corner Detection 2D Image Features Corners are important two dimensional features. Two dimensional image features are interesting local structures. They include junctions of dierent types Slide 3 They

More information