Optimal state estimation for numerical weather prediction using reduced order models
|
|
- Maurice Harrell
- 6 years ago
- Views:
Transcription
1 Optimal state estimation for numerical weather prediction using reduced order models A.S. Lawless 1, C. Boess 1, N.K. Nichols 1 and A. Bunse-Gerstner 2 1 School of Mathematical and Physical Sciences, University of Reading, U.K. 2 ZeTeM, Universtaet Bremen, Germany.
2 Outline 1 4D-Var and Gauss-Newton 2 Model reduction in data assimilation 3 Extension to unstable systems 4 Numerical results 5 Conclusions
3 4-Dimensional Variational Data Assimilation (4D-Var) Aim To find the best estimate of the true state of a system consistent with observations distributed in time and with the system dynamics. x x b t 0 Time t N
4 4-Dimensional Variational Data Assimilation (4D-Var) Aim To find the best estimate of the true state of a system consistent with observations distributed in time and with the system dynamics. x x a x b t 0 Time t N
5 Nonlinear least squares problem We consider a dynamical system x i+1 = M i (x i ) y i = H(x i ) + η i with M i Model dynamics x i State vector O( ) y i Observations at time t i O( ) η i Observational noise Observation operator H i
6 Nonlinear least squares problem Then the data assimilation problem is to minimize J [x 0 ] = 1 2 (x 0 x b ) T B 1 0 (x 0 x b )+ 1 2 N i=0 subject to the dynamical system, where x b A priori (background) estimate B 0 Background error covariance matrix Observation error covariance matrix R i (H i [x i ] y i ) T R 1 i (H i [x i ] y i )
7 This can be written as a standard nonlinear least squares problem J [x] = 1 2 f(x) 2 2 with f(x 0 ) = B (x 0 x b ) R (H 0 [x 0 ] y 0 ). R 1 2 N (H N [x N ] y N )
8 Gauss-Newton method (Incremental 4D-Var) For the nonlinear least squares problem J [x] = 1 2 f(x) 2 2 the Gauss-Newton method is equivalent to the following iteration Start with iterate x (0) For k = 0,...,K 1 Solve min δx J(x (k) )δx + f(x (k) ) 2 2 Update x (k+1) = x (k) + δx
9 For the data assimilation problem the inner cost function is J (k) [δx (k) 0 ] = 1 2 (δx(k) 0 [x b x 0 (k) ]) T B 1 0 (δx(k) 0 [x b x 0 (k) ]) N i=0 (H i δx (k) i d (k) i ) T R 1 (H i δx (k) i i d (k) i ) with d i = y i H i [x i ], subject to the linear dynamical system δx i+1 = M i δx i d i = H i δx i.
10 In practice the inner problem can only be solved approximately. Common approximations include 1 Truncate iterative solution, so that exact minimum not found. 2 Use approximate linear model. 3 Solve the problem in a restricted space (usually a lower spatial resolution). Gratton, Lawless and Nichols (SIOPT, 2007) proved conditions for convergence of Gauss-Newton under approximations 1 & 2. In this work we consider approximation 3.
11 Low resolution G-N in data assimilation We introduce linear restriction operators U T i R r n such that δˆx i = U T i δx i R r, r < n. We also define linear prolongation operators V i R n r, with U T i V i = I r and V i U T i a projection.
12 Low resolution G-N in data assimilation We introduce linear restriction operators U T i R r n such that δˆx i = U T i δx i R r, r < n. We also define linear prolongation operators V i R n r, with U T i V i = I r and V i U T i a projection. We can then define a restricted dynamical system in R r δˆx i+1 = ˆM i δˆx i, ˆd i = Ĥ i δˆx i, where and V i ˆM i U T i approximates M i Ĥ i U T i approximates H i.
13 We solve the minimization problem in R r ˆ J (k) [δˆx (k) 0 ] = 1 2 (δˆx(k) 0 U T 0 [x b x 0 (k) ]) T ˆB 1 0 (δˆx(k) 0 U T 0 [x b x (k) 0 ]) + 1 N (Ĥ i δˆx (k) i d (k) i ) T R 1 i (Ĥ i δˆx (k) i 2 i=0 d (k) i ), subject to the restricted system, and then update the G-N iterate using δx (k) 0 = V 0 δˆx (k) 0
14 In practice the restricted system is chosen to be a low resolution version of the full system. This does not take into account the dynamics or observations. Can we do better using ideas of model reduction?
15 Model reduction Given a time-invariant linear system δx i+1 d i = Mδx i + Gu i = Hδx i where u i N(0,W), find projection matrices U,V with U T V = I r and r << n, such that the reduced order system δˆx i+1 ˆd i = U T MV i δˆx i + U T Gu i = HV i δˆx i approximates the full system with δx i Vδˆx i.
16 We seek matrices U,V such that we minimize { ] lim [ˆd [ˆd T ] } E i d i R 1 i d i i over all inputs of normalized unit length.
17 We seek matrices U,V such that we minimize { ] lim [ˆd [ˆd T ] } E i d i R 1 i d i i over all inputs of normalized unit length. Bernstein et al. (1986) derive necessary conditions for a minimum, but it is not practicable to obtain the optimal matrices. The method of balanced truncation finds an approximation solution, with error bounded in terms of the Hankel singular values.
18 Balanced truncation Balanced truncation removes states which are least affected by inputs and which have least effect on outputs. There are 2 steps 1 Balancing - Transform the system into one in which these states are the same. 2 Truncation - Truncate states related to the smallest singular values of the transformed covariance matrices (Hankel singular values).
19 Step 1: The balancing transformation simultaneously diagonalizes the state covariance matrices P and Q associated with the inputs and outputs. This requires the solution of the Stein equations P = MPM T + GG T, Q = M T QM + H T R 1 H. The transformation Ψ, given by the matrix of eigenvectors of PQ, is used to transform the system into balanced form.
20 Step 1: The balancing transformation simultaneously diagonalizes the state covariance matrices P and Q associated with the inputs and outputs. This requires the solution of the Stein equations P = MPM T + GG T, Q = M T QM + H T R 1 H. The transformation Ψ, given by the matrix of eigenvectors of PQ, is used to transform the system into balanced form. Step 2: The balanced system is truncated using U T = [I r,0]ψ 1, V = Ψ [ Ir 0 ].
21 Application to data assimilation In data assimilation we have the (time-invariant) stochastic dynamical system δx i+1 d i = Mδx i = Hδx i with stochastic initial condition δx 0 N(0,B 0 ). To apply model reduction we need to write this in the form δx i+1 d i = Mδx i + Gu i = Hδx i
22 We write our original system in the form δx i+1 d i = Mδx i + Gũ i = Hδx i with and For 4D-Var we have G = 0. G = [G I n ], ũ i = [ ui [ W 0 ũ i N(0, W), W = 0 B 0 v i ] ].
23 Reduced order cost function The reduced order minimization problem is then min ˆ J (k) [δˆx (k) 0 ] = 1 2 (δˆx(k) 0 U T [x b x 0 (k) ]) T (U T B 0 U) 1 (δˆx (k) 0 U T [x b x (k) 0 ]) + 1 N (HVδˆx (k) i d (k) i ) T R 1 (HVδˆx (k) i d (k) i ), 2 i=0 subject to and set δˆx i+1 = U T MVδˆx i δx (k) 0 = Vδˆx (k) 0.
24 Extension to unstable systems Previously we have shown the benefit of using reduced order models in data assimilation when the system matrix is stable (Lawless et. al, 2008, Mon. Wea. Rev.; Lawless et. al, 2008, Int. J. Numer. Methods in Fluids.) But How can we choose the projection operators U T and V when the system matrix is unstable?
25 We consider 3 methods: 1 Choose U T as a low resolution spatial restriction. 2 Balanced truncation, with standard splitting for unstable systems. 3 New α-bounded balanced truncation for unstable systems.
26 Standard extension of balanced truncation Divide the system S into its stable and unstable parts S = S + + S Apply the usual balanced truncation method to the stable part S + (S + ) red Combine the reduced stable part with the unchanged unstable part Ŝ = (S + ) red + S
27 Standard extension of balanced truncation Divide the system S into its stable and unstable parts S = S + + S Apply the usual balanced truncation method to the stable part S + (S + ) red Combine the reduced stable part with the unchanged unstable part Ŝ = (S + ) red + S Problem: Only works if the number of unstable poles is small.
28 New approach for unstable systems I Boess, Bunse-Gerstner and Nichols (2010) introduce an extension of balanced truncation for α-bounded systems, i.e. systems whose eigenvalues lie in a circle of radius α > 1. 1 Find α R + such that all eigenvalues of the system matrix M lie inside an α-circle around the origin. 2 Choose a reduction order r n. 3 Apply model reduction to the α-scaled system δx (α) i+1 = M α δx (α) i + G α u i d (α) i = H α δx (α) i with M α := M/α, G α := G/ α, H α := H/ α.
29 New approach for unstable systems II This yields projection matrices U α,v α R n r and a low order system: δˆx (α) i+1 = ˆM α δˆx (α) i + G α u i ˆd (α) i = Ĥ α δˆx (α) i with ˆM α := U T αm α V α, Ĝ α := U T α G α and Ĥ α := H α V α. 4 The scaling α is then reversed to give the reduced system δˆx i+1 = ˆMδˆx i + Gu i ˆd i = Ĥδˆx i with ˆM = α ˆM α, G = α G α and Ĥ = αĥ α.
30 Reduced order cost function The reduced order minimization problem is now min ˆ J (k) [δˆx (k) 0 ] = 1 2 (δˆx(k) 0 U T α [x b x 0 (k) ]) T (U T αb 0 U α ) 1 (δˆx (k) 0 U T α [x b x (k) 0 ]) + 1 N (HV α δˆx (k) i d (k) i ) T R 1 (HV α δˆx (k) i d (k) i ), 2 i=0 subject to the reduced system. We then set δx (k) 0 = V α δˆx (k) 0.
31 Numerical results - 2D Eady model The Eady model is a 2D x z model describing the evolution of a buoyancy b on the boundaries and a potential vorticity q in the interior. We have the evolution equations ( t + z x ) b = ψ x, on z = ±1, x [0, X], 2 and ( t + z ) q = 0, in z [ 1 x 2, 1 ], x [0, X], 2 where the streamfunction ψ satisfies 2 ψ x ψ z 2 = q, with boundary conditions ψ z = b, in z [ 1 2, 1 ], x [0, X], 2 on z = ±1, x [0, X]. 2
32 Methodology Define the true solution at the initial time. Run the model to generate synthetic observations at times t 0,...,t 5 (12h window). Define a linear least squares problem using these observations and a background field. Solve the reduced order least squares problem using low resolution model; reduced order system using standard balanced truncation; reduced order system using new method. Lift solution to full resolution and compare with solution to full linear least squares problem.
33 Experiment 1: Zero PV When the interior potential voriticity is zero, the only evolution is on the upper and lower boundaries. We consider 20 horizontal grid points (system size 40). True state is the eigenvector associated with the largest eigenvalue of the model. Zero background state with inverse Laplacian covariance matrix of variance 1. Observation error covariance matrix is identity. Observations taken of buoyancy on lower boundary only. Reduced problem solved with size 20.
34 0.2 Eady40, lower buoyancy, α= 1.12: Solution plot 0.15 x 0 x 0 (low), order x 0 (std), order 20 x 0 (alpha), order solution vector components 10 0 Eady40, lower buoyancy, α= 1.12: Error plot error e 1 =x 0 x 0 (low), order 20 e 2 =x 0 x 0 (std), order 20 e 3 =x 0 x 0 (alpha), order vector components Figure: Lower boundary solution (top) and error (bottom).
35 0.2 x 0 Eady40, upper buoyancy, α= 1.12: Solution plot x 0 (low), order 20 x 0 (std), order 20 x 0 (alpha), order solution vector components 10 0 Eady40, upper buoyancy, α= 1.12: Error plot error e 1 =x 0 x 0 (low), order 20 e 2 =x 0 x 0 (std), order 20 e 3 =x 0 x 0 (alpha), order vector components Figure: Upper boundary solution (top) and error (bottom).
36 Eigenvalues of α reduced M, α = 1.12 Eigenvalues of low resolution M eigs M eigs M 0.6 eigs M eigs M α 0.6 low unit circle unit circle imaginary axis imaginary axis real axis real axis Eigenvalues of standard reduced M eigs M eigs M std unit circle 0.4 imaginary axis real axis Figure: Comaprison of eigenvalues with full system for α-reduced (top left), standard reduced (top right) and low resolution (bottom).
37 Experiment 2: Non-zero PV We now consider a case when there is evolution of the PV in the interior. The full system size is now 1040 (80 horizontal grid points). True solution is non-modal growing mode. Observations are taken at every second grid point. Background error variances are 10 5 on the boundaries and 1 in the interior. Reduced problem solved with size 520.
38 2e 006 2e Contour plot of true solution x e e 006 6e 006 8e 006 2e 006 4e e 006 6e Contour plot of lifted low res solution, k=520 x e e 006 4e 006 6e e e 006 6e 006 4e Figure: True solution (top) and low resolution (bottom).
39 2e 006 2e 007 2e Contour plot of true solution x e e 006 6e 006 8e 006 2e 006 4e e 006 6e Contour plot of lifted mod red solution using standard BT, k=520 x e 007 2e 007 4e 007 6e 007 4e 007 2e 007 4e 007 2e 007 2e Figure: True solution (top) and standard BT (bottom).
40 2e 006 2e 006 2e 006 2e Contour plot of true solution x e e 006 6e 006 8e 006 2e 006 4e e 006 6e Contour plot of lifted mod red solution using α BT, α = 1.13, k=520 x e e 006 6e 006 8e 006 2e 006 4e e 006 6e Figure: True solution (top) and α BT (bottom).
41 Conclusions In many applications dynamical systems may be unstable over finite time intervals.
42 Conclusions In many applications dynamical systems may be unstable over finite time intervals. A new method of α-bounded balanced truncation has been incorporated into 4D-Var.
43 Conclusions In many applications dynamical systems may be unstable over finite time intervals. A new method of α-bounded balanced truncation has been incorporated into 4D-Var. Numerical results show a more accurate solution than using low resolution models or standard model reduction.
44 Conclusions In many applications dynamical systems may be unstable over finite time intervals. A new method of α-bounded balanced truncation has been incorporated into 4D-Var. Numerical results show a more accurate solution than using low resolution models or standard model reduction. In part this can be explained by preservation of the model eigenvalues.
45 Remaining issues In order to be implemented in real systems the method still requires some further work:
46 Remaining issues In order to be implemented in real systems the method still requires some further work: Extension to time-varying systems.
47 Remaining issues In order to be implemented in real systems the method still requires some further work: Extension to time-varying systems. Efficient reduction methods for large-scale systems.
48 Remaining issues In order to be implemented in real systems the method still requires some further work: Extension to time-varying systems. Efficient reduction methods for large-scale systems. Update of reduced order system within iterative process.
49 Remaining issues In order to be implemented in real systems the method still requires some further work: Extension to time-varying systems. Efficient reduction methods for large-scale systems. Update of reduced order system within iterative process. Nevertheless results indicate that this is research that is worth pursuing.
50 References I Boess, C., Nichols, N.K. and Bunse-Gerstner, A. (2010), Model reduction for discrete unstable control systems using a balanced truncation approach. Preprint MPS , Department of Mathematics, University of Reading. Gratton, S., Lawless, A.S. and Nichols, N.K. (2007), Approximate Gauss-Newton methods for nonlinear least squares problems. SIAM J. on Optimization, 18, Lawless, A.S., Nichols, N.K., Boess, C. and Bunse-Gerstner, A. (2008), Using model reduction methods within incremental 4D-Var. Mon. Wea. Rev., 136,
51 References II Lawless, A.S., Nichols, N.K., Boess, C. and Bunse-Gerstner, A. (2008), Approximate Gauss-Newton methods for optimal state estimation using reduced order models. Int. J. Numer. Methods in Fluids, 56, Boess, C., Lawless, A.S., Nichols, N.K. and Bunse-Gerstner, A. (2011), State estimation using model order reduction for unstable systems. Comput. Fluids, 46,
Using Model Reduction Methods within Incremental Four-Dimensional Variational Data Assimilation
APRIL 2008 L A W L E S S E T A L. 1511 Using Model Reduction Methods within Incremental Four-Dimensional Variational Data Assimilation A. S. LAWLESS AND N. K. NICHOLS Department of Mathematics, University
More informationOptimal Control and Estimation:
Optimal Control and Estimation: An ICFD Retrospective N.K. Nichols University of Reading ICFD 25th Anniversary With thanks to: CEGB / National Power SSEB BP Exploration British Gas Institute of Hydrology
More informationParallel Algorithms for Four-Dimensional Variational Data Assimilation
Parallel Algorithms for Four-Dimensional Variational Data Assimilation Mie Fisher ECMWF October 24, 2011 Mie Fisher (ECMWF) Parallel 4D-Var October 24, 2011 1 / 37 Brief Introduction to 4D-Var Four-Dimensional
More informationMathematical Concepts of Data Assimilation
Mathematical Concepts of Data Assimilation N.K. Nichols 1 Introduction Environmental systems can be realistically described by mathematical and numerical models of the system dynamics. These models can
More informationA comparison of the model error and state formulations of weak-constraint 4D-Var
A comparison of the model error and state formulations of weak-constraint 4D-Var A. El-Said, A.S. Lawless and N.K. Nichols School of Mathematical and Physical Sciences University of Reading Work supported
More informationRelationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF
Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF Eugenia Kalnay and Shu-Chih Yang with Alberto Carrasi, Matteo Corazza and Takemasa Miyoshi 4th EnKF Workshop, April 2010 Relationship
More informationDepartment of Mathematics and Statistics
School of Mathematical and Physical Sciences Department of Mathematics and Statistics Preprint MPS-01-17 15 August 01 Data assimilation with correlated observation errors: analysis accuracy with approximate
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. OPTIM. Vol. 18, No. 1, pp. 106 13 c 007 Society for Industrial and Applied Mathematics APPROXIMATE GAUSS NEWTON METHODS FOR NONLINEAR LEAST SQUARES PROBLEMS S. GRATTON, A. S. LAWLESS, AND N. K.
More informationWeak Constraints 4D-Var
Weak Constraints 4D-Var Yannick Trémolet ECMWF Training Course - Data Assimilation May 1, 2012 Yannick Trémolet Weak Constraints 4D-Var May 1, 2012 1 / 30 Outline 1 Introduction 2 The Maximum Likelihood
More informationRelationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF
Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF Eugenia Kalnay and Shu-Chih Yang with Alberto Carrasi, Matteo Corazza and Takemasa Miyoshi ECODYC10, Dresden 28 January 2010 Relationship
More informationModel Reduction for Unstable Systems
Model Reduction for Unstable Systems Klajdi Sinani Virginia Tech klajdi@vt.edu Advisor: Serkan Gugercin October 22, 2015 (VT) SIAM October 22, 2015 1 / 26 Overview 1 Introduction 2 Interpolatory Model
More informationParameter-Dependent Eigencomputations and MEMS Applications
Parameter-Dependent Eigencomputations and MEMS Applications David Bindel UC Berkeley, CS Division Parameter-Dependent Eigencomputations and MEMS Applications p.1/36 Outline Some applications Mathematical
More informationSubspace-based Identification
of Infinite-dimensional Multivariable Systems from Frequency-response Data Department of Electrical and Electronics Engineering Anadolu University, Eskişehir, Turkey October 12, 2008 Outline 1 2 3 4 Noise-free
More informationNonlinear error dynamics for cycled data assimilation methods
Nonlinear error dynamics for cycled data assimilation methods A J F Moodey 1, A S Lawless 1,2, P J van Leeuwen 2, R W E Potthast 1,3 1 Department of Mathematics and Statistics, University of Reading, UK.
More informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
More informationSolutions Preliminary Examination in Numerical Analysis January, 2017
Solutions Preliminary Examination in Numerical Analysis January, 07 Root Finding The roots are -,0, a) First consider x 0 > Let x n+ = + ε and x n = + δ with δ > 0 The iteration gives 0 < ε δ < 3, which
More informationNumerical Weather prediction at the European Centre for Medium-Range Weather Forecasts (2)
Numerical Weather prediction at the European Centre for Medium-Range Weather Forecasts (2) Time series curves 500hPa geopotential Correlation coefficent of forecast anomaly N Hemisphere Lat 20.0 to 90.0
More informationPractical Aspects of Ensemble-based Kalman Filters
Practical Aspects of Ensemble-based Kalman Filters Lars Nerger Alfred Wegener Institute for Polar and Marine Research Bremerhaven, Germany and Bremen Supercomputing Competence Center BremHLR Bremen, Germany
More informationVariational data assimilation
Background and methods NCEO, Dept. of Meteorology, Univ. of Reading 710 March 2018, Univ. of Reading Bayes' Theorem Bayes' Theorem p(x y) = posterior distribution = p(x) p(y x) p(y) prior distribution
More informationVariational assimilation Practical considerations. Amos S. Lawless
Variational assimilation Practical considerations Amos S. Lawless a.s.lawless@reading.ac.uk 4D-Var problem ] [ ] [ 2 2 min i i i n i i T i i i b T b h h J y R y B,,, n i i i i f subject to Minimization
More informationSuperparameterization and Dynamic Stochastic Superresolution (DSS) for Filtering Sparse Geophysical Flows
SP and DSS for Filtering Sparse Geophysical Flows Superparameterization and Dynamic Stochastic Superresolution (DSS) for Filtering Sparse Geophysical Flows Presented by Nan Chen, Michal Branicki and Chenyue
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More information, b = 0. (2) 1 2 The eigenvectors of A corresponding to the eigenvalues λ 1 = 1, λ 2 = 3 are
Quadratic forms We consider the quadratic function f : R 2 R defined by f(x) = 2 xt Ax b T x with x = (x, x 2 ) T, () where A R 2 2 is symmetric and b R 2. We will see that, depending on the eigenvalues
More informationQuantifying observation error correlations in remotely sensed data
Quantifying observation error correlations in remotely sensed data Conference or Workshop Item Published Version Presentation slides Stewart, L., Cameron, J., Dance, S. L., English, S., Eyre, J. and Nichols,
More informationThe ECMWF Hybrid 4D-Var and Ensemble of Data Assimilations
The Hybrid 4D-Var and Ensemble of Data Assimilations Lars Isaksen, Massimo Bonavita and Elias Holm Data Assimilation Section lars.isaksen@ecmwf.int Acknowledgements to: Mike Fisher and Marta Janiskova
More information5.1 2D example 59 Figure 5.1: Parabolic velocity field in a straight two-dimensional pipe. Figure 5.2: Concentration on the input boundary of the pipe. The vertical axis corresponds to r 2 -coordinate,
More informationModel reduction for linear systems by balancing
Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,
More informationBackground and observation error covariances Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience
Background and observation error covariances Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience Sarah Dance School of Mathematical and Physical Sciences, University of Reading
More informationThe Important State Coordinates of a Nonlinear System
The Important State Coordinates of a Nonlinear System Arthur J. Krener 1 University of California, Davis, CA and Naval Postgraduate School, Monterey, CA ajkrener@ucdavis.edu Summary. We offer an alternative
More informationOptimal Resolution of State Space in Chaotic Hyperbolic 2-dimensional Maps
Optimal Resolution of State Space in Chaotic Hyperbolic -dimensional Maps Gable M. Wadsworth and Daniel T. Kovari (Dated: May, ) Noise is a factor in every physical system. To deal with this physical constraint
More informationCoupled atmosphere-ocean data assimilation in the presence of model error. Alison Fowler and Amos Lawless (University of Reading)
Coupled atmosphere-ocean data assimilation in the presence of model error Alison Fowler and Amos Lawless (University of Reading) Introduction Coupled DA methods are being developed to initialise forecasts
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationData Assimilation for Weather Forecasting: Reducing the Curse of Dimensionality
Data Assimilation for Weather Forecasting: Reducing the Curse of Dimensionality Philippe Toint (with S. Gratton, S. Gürol, M. Rincon-Camacho, E. Simon and J. Tshimanga) University of Namur, Belgium Leverhulme
More informationRegularization Parameter Estimation for Least Squares: A Newton method using the χ 2 -distribution
Regularization Parameter Estimation for Least Squares: A Newton method using the χ 2 -distribution Rosemary Renaut, Jodi Mead Arizona State and Boise State September 2007 Renaut and Mead (ASU/Boise) Scalar
More informationProblem set 5 solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 24: MODEL REDUCTION Problem set 5 solutions Problem 5. For each of the stetements below, state
More informationThe Impact of Background Error on Incomplete Observations for 4D-Var Data Assimilation with the FSU GSM
The Impact of Background Error on Incomplete Observations for 4D-Var Data Assimilation with the FSU GSM I. Michael Navon 1, Dacian N. Daescu 2, and Zhuo Liu 1 1 School of Computational Science and Information
More informationThe University of Reading
The University of Reading Radial Velocity Assimilation and Experiments with a Simple Shallow Water Model S.J. Rennie 2 and S.L. Dance 1,2 NUMERICAL ANALYSIS REPORT 1/2008 1 Department of Mathematics 2
More informationDOA Estimation using MUSIC and Root MUSIC Methods
DOA Estimation using MUSIC and Root MUSIC Methods EE602 Statistical signal Processing 4/13/2009 Presented By: Chhavipreet Singh(Y515) Siddharth Sahoo(Y5827447) 2 Table of Contents 1 Introduction... 3 2
More informationCHAPTER 5. Basic Iterative Methods
Basic Iterative Methods CHAPTER 5 Solve Ax = f where A is large and sparse (and nonsingular. Let A be split as A = M N in which M is nonsingular, and solving systems of the form Mz = r is much easier than
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More informationData Assimilation: Finding the Initial Conditions in Large Dynamical Systems. Eric Kostelich Data Mining Seminar, Feb. 6, 2006
Data Assimilation: Finding the Initial Conditions in Large Dynamical Systems Eric Kostelich Data Mining Seminar, Feb. 6, 2006 kostelich@asu.edu Co-Workers Istvan Szunyogh, Gyorgyi Gyarmati, Ed Ott, Brian
More informationA Local Ensemble Kalman Filter for Atmospheric Data Assimilation
Tellus 000, 000 000 (0000) Printed 1 April 2004 (Tellus LATEX style file v2.2) A Local Ensemble Kalman Filter for Atmospheric Data Assimilation By EDWARD OTT 1, BRIAN R. HUNT 2, ISTVAN SZUNYOGH 3, ALEKSEY
More information1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0
Numerical Analysis 1 1. Nonlinear Equations This lecture note excerpted parts from Michael Heath and Max Gunzburger. Given function f, we seek value x for which where f : D R n R n is nonlinear. f(x) =
More informationMethods for sparse analysis of high-dimensional data, II
Methods for sparse analysis of high-dimensional data, II Rachel Ward May 26, 2011 High dimensional data with low-dimensional structure 300 by 300 pixel images = 90, 000 dimensions 2 / 55 High dimensional
More informationBasic Math for
Basic Math for 16-720 August 23, 2002 1 Linear Algebra 1.1 Vectors and Matrices First, a reminder of a few basic notations, definitions, and terminology: Unless indicated otherwise, vectors are always
More informationOn Dominant Poles and Model Reduction of Second Order Time-Delay Systems
On Dominant Poles and Model Reduction of Second Order Time-Delay Systems Maryam Saadvandi Joint work with: Prof. Karl Meerbergen and Dr. Elias Jarlebring Department of Computer Science, KULeuven ModRed
More informationRecent Developments in Numerical Methods for 4d-Var
Recent Developments in Numerical Methods for 4d-Var Mike Fisher Slide 1 Recent Developments Numerical Methods 4d-Var Slide 2 Outline Non-orthogonal wavelets on the sphere: - Motivation: Covariance Modelling
More informationMASS, STIFFNESS AND DAMPING IDENTIFICATION OF A TWO-STORY BUILDING MODEL
COMPDYN 2 3 rd ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, M. Fragiadakis, V. Plevris (eds.) Corfu, Greece, 25-28 May 2 MASS,
More informationNon-linear least squares
Non-linear least squares Concept of non-linear least squares We have extensively studied linear least squares or linear regression. We see that there is a unique regression line that can be determined
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 19: More on Arnoldi Iteration; Lanczos Iteration Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical Analysis I 1 / 17 Outline 1
More informationIntroduction to Quantum Mechanics Physics Thursday February 21, Problem # 1 (10pts) We are given the operator U(m, n) defined by
Department of Physics Introduction to Quantum Mechanics Physics 5701 Temple University Z.-E. Meziani Thursday February 1, 017 Problem # 1 10pts We are given the operator Um, n defined by Ûm, n φ m >< φ
More informationPractical Optimization: Basic Multidimensional Gradient Methods
Practical Optimization: Basic Multidimensional Gradient Methods László Kozma Lkozma@cis.hut.fi Helsinki University of Technology S-88.4221 Postgraduate Seminar on Signal Processing 22. 10. 2008 Contents
More informationNeural Network Training
Neural Network Training Sargur Srihari Topics in Network Training 0. Neural network parameters Probabilistic problem formulation Specifying the activation and error functions for Regression Binary classification
More informationDATA ASSIMILATION FOR FLOOD FORECASTING
DATA ASSIMILATION FOR FLOOD FORECASTING Arnold Heemin Delft University of Technology 09/16/14 1 Data assimilation is the incorporation of measurement into a numerical model to improve the model results
More informationGaussian Filtering Strategies for Nonlinear Systems
Gaussian Filtering Strategies for Nonlinear Systems Canonical Nonlinear Filtering Problem ~u m+1 = ~ f (~u m )+~ m+1 ~v m+1 = ~g(~u m+1 )+~ o m+1 I ~ f and ~g are nonlinear & deterministic I Noise/Errors
More informationIterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations
Iterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations Klajdi Sinani Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University
More informationQuantifying observation error correlations in remotely sensed data
Quantifying observation error correlations in remotely sensed data Conference or Workshop Item Published Version Presentation slides Stewart, L., Cameron, J., Dance, S. L., English, S., Eyre, J. and Nichols,
More informationSecond-order cone programming
Outline Second-order cone programming, PhD Lehigh University Department of Industrial and Systems Engineering February 10, 2009 Outline 1 Basic properties Spectral decomposition The cone of squares The
More informationReal Eigenvalue Extraction and the Distance to Uncontrollability
Real Eigenvalue Extraction and the Distance to Uncontrollability Emre Mengi Computer Science Department Courant Institute of Mathematical Sciences New York University mengi@cs.nyu.edu May 22nd, 2006 Emre
More informationRevision of TR-09-25: A Hybrid Variational/Ensemble Filter Approach to Data Assimilation
Revision of TR-9-25: A Hybrid Variational/Ensemble ilter Approach to Data Assimilation Adrian Sandu 1 and Haiyan Cheng 1 Computational Science Laboratory Department of Computer Science Virginia Polytechnic
More informationA HYBRID ENSEMBLE KALMAN FILTER / 3D-VARIATIONAL ANALYSIS SCHEME
A HYBRID ENSEMBLE KALMAN FILTER / 3D-VARIATIONAL ANALYSIS SCHEME Thomas M. Hamill and Chris Snyder National Center for Atmospheric Research, Boulder, Colorado 1. INTRODUCTION Given the chaotic nature of
More informationComparisons between 4DEnVar and 4DVar on the Met Office global model
Comparisons between 4DEnVar and 4DVar on the Met Office global model David Fairbairn University of Surrey/Met Office 26 th June 2013 Joint project by David Fairbairn, Stephen Pring, Andrew Lorenc, Neill
More informationCSE 554 Lecture 7: Alignment
CSE 554 Lecture 7: Alignment Fall 2012 CSE554 Alignment Slide 1 Review Fairing (smoothing) Relocating vertices to achieve a smoother appearance Method: centroid averaging Simplification Reducing vertex
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More informationApproximating Optimal State Estimation
Approximating Optimal State Estimation Brian F. Farrell and Petros J. Ioannou Harvard University Cambridge MA 238, U.S.A farrell@deas.harvard.edu, pji@cc.uoa.gr ABSTRACT Minimizing forecast error requires
More informationModel order reduction of electrical circuits with nonlinear elements
Model order reduction of electrical circuits with nonlinear elements Andreas Steinbrecher and Tatjana Stykel 1 Introduction The efficient and robust numerical simulation of electrical circuits plays a
More informationBALANCING-RELATED MODEL REDUCTION FOR DATA-SPARSE SYSTEMS
BALANCING-RELATED Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Computational Methods with Applications Harrachov, 19 25 August 2007
More informationMACHINE LEARNING. Methods for feature extraction and reduction of dimensionality: Probabilistic PCA and kernel PCA
1 MACHINE LEARNING Methods for feature extraction and reduction of dimensionality: Probabilistic PCA and kernel PCA 2 Practicals Next Week Next Week, Practical Session on Computer Takes Place in Room GR
More informationSTA141C: Big Data & High Performance Statistical Computing
STA141C: Big Data & High Performance Statistical Computing Numerical Linear Algebra Background Cho-Jui Hsieh UC Davis May 15, 2018 Linear Algebra Background Vectors A vector has a direction and a magnitude
More informationVariational Data Assimilation Current Status
Variational Data Assimilation Current Status Eĺıas Valur Hólm with contributions from Mike Fisher and Yannick Trémolet ECMWF NORDITA Solar and stellar dynamos and cycles October 2009 Eĺıas Valur Hólm (ECMWF)
More informationDevelopment of wavelet methodology for weather Data Assimilation
Development of wavelet methodology for weather Data Assimilation Aimé Fournier Thomas Auligné Mesoscale and Microscale Meteorology Division National Center for Atmospheric Research Newton Institute Mathematical
More informationBackground Error Covariance Modelling
Background Error Covariance Modelling Mike Fisher Slide 1 Outline Diagnosing the Statistics of Background Error using Ensembles of Analyses Modelling the Statistics in Spectral Space - Relaxing constraints
More informationScientific Computing WS 2018/2019. Lecture 9. Jürgen Fuhrmann Lecture 9 Slide 1
Scientific Computing WS 2018/2019 Lecture 9 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 9 Slide 1 Lecture 9 Slide 2 Simple iteration with preconditioning Idea: Aû = b iterative scheme û = û
More informationPart 3: Trust-region methods for unconstrained optimization. Nick Gould (RAL)
Part 3: Trust-region methods for unconstrained optimization Nick Gould (RAL) minimize x IR n f(x) MSc course on nonlinear optimization UNCONSTRAINED MINIMIZATION minimize x IR n f(x) where the objective
More informationInverse Singular Value Problems
Chapter 8 Inverse Singular Value Problems IEP versus ISVP Existence question A continuous approach An iterative method for the IEP An iterative method for the ISVP 139 140 Lecture 8 IEP versus ISVP Inverse
More informationFloquet Theory for Internal Gravity Waves in a Density-Stratified Fluid. Yuanxun Bill Bao Senior Supervisor: Professor David J. Muraki August 3, 2012
Floquet Theory for Internal Gravity Waves in a Density-Stratified Fluid Yuanxun Bill Bao Senior Supervisor: Professor David J. Muraki August 3, 212 Density-Stratified Fluid Dynamics Density-Stratified
More informationIndefinite and physics-based preconditioning
Indefinite and physics-based preconditioning Jed Brown VAW, ETH Zürich 2009-01-29 Newton iteration Standard form of a nonlinear system F (u) 0 Iteration Solve: Update: J(ũ)u F (ũ) ũ + ũ + u Example (p-bratu)
More informationEcon 204 Supplement to Section 3.6 Diagonalization and Quadratic Forms. 1 Diagonalization and Change of Basis
Econ 204 Supplement to Section 3.6 Diagonalization and Quadratic Forms De La Fuente notes that, if an n n matrix has n distinct eigenvalues, it can be diagonalized. In this supplement, we will provide
More informationLeast squares: introduction to the network adjustment
Least squares: introduction to the network adjustment Experimental evidence and consequences Observations of the same quantity that have been performed at the highest possible accuracy provide different
More informationGFD II: Balance Dynamics ATM S 542
GFD II: Balance Dynamics ATM S 542 DARGAN M. W. FRIERSON UNIVERSITY OF WASHINGTON, DEPARTMENT OF ATMOSPHERIC SCIENCES WEEK 8 SLIDES Eady Model Growth rates (imaginary part of frequency) Stable for large
More informationMachine vision, spring 2018 Summary 4
Machine vision Summary # 4 The mask for Laplacian is given L = 4 (6) Another Laplacian mask that gives more importance to the center element is given by L = 8 (7) Note that the sum of the elements in the
More informationCorrelations of control variables in variational data assimilation
Quarterly Journal of the Royal Meteorological Society Q. J. R. Meteorol. Soc. 37: 62 63, April 2 A Correlations of control variables in variational data assimilation D. Katz a,a.s.lawless* a,n.k.nichols
More informationComparing Local Ensemble Transform Kalman Filter with 4D-Var in a Quasi-geostrophic model
Comparing Local Ensemble Transform Kalman Filter with 4D-Var in a Quasi-geostrophic model Shu-Chih Yang 1,2, Eugenia Kalnay 1, Matteo Corazza 3, Alberto Carrassi 4 and Takemasa Miyoshi 5 1 University of
More informationSpectral transforms. Contents. ARPEGE-Climat Version 5.1. September Introduction 2
Spectral transforms ARPEGE-Climat Version 5.1 September 2008 Contents 1 Introduction 2 2 Spectral representation 2 2.1 Spherical harmonics....................... 2 2.2 Collocation grid.........................
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationRepresentation of inhomogeneous, non-separable covariances by sparse wavelet-transformed matrices
Representation of inhomogeneous, non-separable covariances by sparse wavelet-transformed matrices Andreas Rhodin, Harald Anlauf German Weather Service (DWD) Workshop on Flow-dependent aspects of data assimilation,
More informationNumerical solutions of nonlinear systems of equations
Numerical solutions of nonlinear systems of equations Tsung-Ming Huang Department of Mathematics National Taiwan Normal University, Taiwan E-mail: min@math.ntnu.edu.tw August 28, 2011 Outline 1 Fixed points
More informationMath Models of OR: The Revised Simplex Method
Math Models of OR: The Revised Simplex Method John E. Mitchell Department of Mathematical Sciences RPI, Troy, NY 12180 USA September 2018 Mitchell The Revised Simplex Method 1 / 25 Motivation Outline 1
More informationReduced-dimension Models in Nonlinear Finite Element Dynamics of Continuous Media
Reduced-dimension Models in Nonlinear Finite Element Dynamics of Continuous Media Petr Krysl, Sanjay Lall, and Jerrold E. Marsden, California Institute of Technology, Pasadena, CA 91125. pkrysl@cs.caltech.edu,
More informationHigh-resolution Parametric Subspace Methods
High-resolution Parametric Subspace Methods The first parametric subspace-based method was the Pisarenko method,, which was further modified, leading to the MUltiple SIgnal Classification (MUSIC) method.
More information4. 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 information7 EQUATIONS OF MOTION FOR AN INVISCID FLUID
7 EQUATIONS OF MOTION FOR AN INISCID FLUID iscosity is a measure of the thickness of a fluid, and its resistance to shearing motions. Honey is difficult to stir because of its high viscosity, whereas water
More informationThe 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 informationDiagnosis of observation, background and analysis-error statistics in observation space
Q. J. R. Meteorol. Soc. (2005), 131, pp. 3385 3396 doi: 10.1256/qj.05.108 Diagnosis of observation, background and analysis-error statistics in observation space By G. DESROZIERS, L. BERRE, B. CHAPNIK
More informationFundamentals of Data Assimila1on
014 GSI Community Tutorial NCAR Foothills Campus, Boulder, CO July 14-16, 014 Fundamentals of Data Assimila1on Milija Zupanski Cooperative Institute for Research in the Atmosphere Colorado State University
More informationRobust Multivariable Control
Lecture 2 Anders Helmersson anders.helmersson@liu.se ISY/Reglerteknik Linköpings universitet Today s topics Today s topics Norms Today s topics Norms Representation of dynamic systems Today s topics Norms
More informationOptimal Interpolation ( 5.4) We now generalize the least squares method to obtain the OI equations for vectors of observations and background fields.
Optimal Interpolation ( 5.4) We now generalize the least squares method to obtain the OI equations for vectors of observations and background fields. Optimal Interpolation ( 5.4) We now generalize the
More information1 Number Systems and Errors 1
Contents 1 Number Systems and Errors 1 1.1 Introduction................................ 1 1.2 Number Representation and Base of Numbers............. 1 1.2.1 Normalized Floating-point Representation...........
More informationDacian N. Daescu 1, I. Michael Navon 2
6.4 REDUCED-ORDER OBSERVATION SENSITIVITY IN 4D-VAR DATA ASSIMILATION Dacian N. Daescu 1, I. Michael Navon 2 1 Portland State University, Portland, Oregon, 2 Florida State University, Tallahassee, Florida
More information