Gaussian Quadrature We use the following formula with weights w i and - PowerPoint PPT Presentation

Dec 06, 2022 •331 likes •384 views

Gaussian Quadrature We use the following formula with weights w i and nodes x i 1 n f ( x ) dx w i f ( x i ) 1 i =1 The nodes and weights are of the form: n (# of points) x i w i 2 0.5773502691896257 1.0000000000000000

Gaussian Quadrature We use the following formula with weights w i and nodes x i � 1 n � f ( x ) dx ≈ w i f ( x i ) − 1 i =1 The nodes and weights are of the form: n (# of points) x i w i 2 0.5773502691896257 1.0000000000000000 -0.5773502691896257 1.0000000000000000 3 0.7745966692414834 0.5555555555555556 0 0.8888888888888888 -0.7745966692414834 0.5555555555555556 4 0.8611363115940525 0.3478548451374544 0.3399810435848563 0.6521451548625460 -0.3399810435848563 0.6521451548625460 -0.8611363115940525 0.3478548451374544 · · · · · · · · · These can be found in most Numerical analysis texts. September 24, 2019 1 / 4
Matlab implementation We approximate an integral over [ a , b ] as follows: � b n � f dx = w i f (˜ ˜ x i ) a i =1 b − a x i = ( b + a ) + x i ( b − a ) where ˜ w i = w i and ˜ 2 2 The function lgwt computes ˜ w i and ˜ x i for any n on an interval [ a , b ]. [x,w]=lgwt(4,0,pi); We can then specify the function and approximate the integral >> f=sin(x); >> approx_int = sum(w.*f) approx_int = 1.999984228457722 September 24, 2019 2 / 4
Comparison to Composite Trapezoidal rule Gaussian Quadrature is far superior! n Trapeziodal Error Gaussian Error 2 0.429203673205103 0.064180425348863 4 0.103881102062960 1.577154227838662e-05 8 0.025768398054449 4.440892098500626e-15 16 0.006429656227660 4.440892098500626e-16 September 24, 2019 3 / 4
All is not lost for Composite methods! � 2 π e sin( x ) dx 0 n Trapeziodal Error Gaussian Error 2 1.67174121383325e+00 1.52774484383558e+00 4 3.43969188091968e-02 2.02753056344147e-01 8 1.25168893738703e-06 1.222935441335338e-03 16 6.21724893790088e-15 4.724131308364576e-09 September 24, 2019 4 / 4
They perfom quite well for periodic integrals September 24, 2019 4 / 4

Recommend

Adaptive Quadrature for Nystrm Yuchen Su 12/6/17 Quadrature Quadrature Quadrature

Adaptive Quadrature for Nystrm Yuchen Su 12/6/17 Quadrature Quadrature Quadrature Quadrature Quadrature Quadrature Composite Newton-Cotes Quadrature Composite Gaussian Error Estimation When should we be satisfied with our Quadrature?

792 views • 66 slides

Gaussian Quadrature September 25, 2011 Interpolation Approximation of integrals Approximation

Interpolation Gaussian Quadrature September 25, 2011 Interpolation Approximation of integrals Approximation of integrals by quadrature Many definite integrals cannot be computed in closed form, and must be approximated numerically.

197 views • 16 slides

Numerical integration (a.k.a quadrature formulas or quadrature rules) Quadrature rules are used

Numerical integration (a.k.a quadrature formulas or quadrature rules) Quadrature rules are used to approximate integrals of functions that we are not able to compute exactly. Given g : [ c , d ] R , the most common quadrature rules look like

310 views • 19 slides

4. Numerical Quadrature Where analytical abilities end . . . 4. Numerical Quadrature Numerical

Preliminary Remarks Simple and Composite Rules Extrapolation Other Approaches Applications of Quadrature 4. Numerical Quadrature Where analytical abilities end . . . 4. Numerical Quadrature Numerical Programming I (for CSE), Hans-Joachim

477 views • 32 slides

From numerical quadrature to Pad approximation Claude Brezinski University of Lille - France

From numerical quadrature to Pad approximation Claude Brezinski University of Lille - France Numerical quadrature Interpolatory quadratures Gaussian quadratures Error estimates Pad approximation of power series

658 views • 35 slides

4. Numerical Quadrature: Where analytical abilities end . . . 4. Numerical Quadrature: Where

Preliminary Remarks Simple and Composite Rules Extrapolation Other Approaches Applications of Quadrature 4. Numerical Quadrature: Where analytical abilities end . . . 4. Numerical Quadrature: Where analytical abilities end . . . Numerical

463 views • 32 slides

Quadrature CS3220 - Summer 2008 Jonathan Kaldor What is Quadrature? Numerical evaluation /

Quadrature CS3220 - Summer 2008 Jonathan Kaldor What is Quadrature? Numerical evaluation / approximation of definite integrals ab f(x) d x What is Quadrature? Recall from calculus that we can compute the integral of a function

710 views • 40 slides

Lesson 2 Quadrature rules 1 Quadrature is the computation of integrals Z b f ( x )

Lesson 2 Quadrature rules 1 Quadrature is the computation of integrals Z b f ( x ) d x a Only in very special cases can integrals be computed exactly, otherwise, we must compute them numerically (i.e., approximately) The

626 views • 50 slides

Gaussian Filter The Gaussian filter 1 2 1 A Gaussian kernel gives less 1 2 4 2 weight to

Gaussian Filter The Gaussian filter 1 2 1 A Gaussian kernel gives less 1 2 4 2 weight to pixels further from 16 the center of the window 1 2 1 This kernel is an approximation of a Gaussian function Gaussian filtering versus mean

405 views • 11 slides

Quadrature Domains in Complex Variables Alan Legg Department of Mathematical Sciences, IPFW

Quadrature Domains in Complex Variables Alan Legg Department of Mathematical Sciences, IPFW October 8, 2016 Quadrature Domains Quadrature domains (QDs) are a special type of domain, traditionally in C , where integrating some class of

896 views • 86 slides

Fast Bayesian automatic Fast Bayesian automatic adaptive quadrature adaptive quadrature Gh.

Fast Bayesian automatic Fast Bayesian automatic adaptive quadrature adaptive quadrature Gh. Adam, S. Adam LIT-JINR Dubna, Russia & IFIN-HH Bucharest, Romania Invited Lecture at the RO-LCG 2014 conference, November 3-5, 2014 Outline

551 views • 52 slides

Method of Finite Elements I PhD Candidate - Charilaos Mylonas HIL H33.1 Quadrature and Boundary

Method of Finite Elements I PhD Candidate - Charilaos Mylonas HIL H33.1 Quadrature and Boundary Conditions , 26 March, 2018 Institute of Structural Engineering Method of Finite Elements I 1 Quadrature Boundary conditions Outline Quadrature

2.03k views • 27 slides

Non-Gaussian likelihoods for Gaussian Processes Alan Saul Outline Motivation Non-Gaussian

Non-Gaussian likelihoods for Gaussian Processes Alan Saul Outline Motivation Non-Gaussian posteriors Approximate methods Laplace approximation Variational bayes Expectation propagation Comparisons GP regression - recap so far Model the

1.41k views • 113 slides

about quadrature Nick Trefethen, March 2016 With thanks to Andr Weideman for #2, #5, #8, #9

Ten things you should know about quadrature Nick Trefethen, March 2016 With thanks to Andr Weideman for #2, #5, #8, #9 Nick Hale for #4, #5, #7, #8 Nick Higham for #5, #8 Folkmar Bornemann for #6. 1/26 A lifetime ago 1. Gauss quadrature

479 views • 27 slides

Quadrature on the positive real line with quasi and pseudo orthogonal rational functions Adhemar

Gauss-type Quadrature ORF with a parameter Computation Quadrature on the positive real line with quasi and pseudo orthogonal rational functions Adhemar Bultheel (joint work with P. Gonz alez-Vera, E. Hendriksen, O. Nj astad) Department

843 views • 80 slides

Quadrature Generation Techniques in CMOS Relaxation Oscillators S. Aniruddhan Indian Institute

Quadrature Generation Techniques in CMOS Relaxation Oscillators S. Aniruddhan Indian Institute of Technology Madras Chennai, India ISCAS 2012 Outline Introduction & Motivation Quadrature Relaxation Oscillators (QRXO)

269 views • 22 slides

Numerical Integration = Quadrature 1. Numerical integration plays key roles for computing, and

Numerical Integration = Quadrature 1. Numerical integration plays key roles for computing, and some well-known functions are definied as integrals. 2. Numerical integration, also known as quadrature : Given a sample of n points for a function f ,

651 views • 8 slides

Multi-rate Signal Processing 6. Quadrature Mirror Filter (QMF) Bank Electrical & Computer

6 Quadrature Mirror Filter (QMF) Bank Appendix: Detailed Derivations Multi-rate Signal Processing 6. Quadrature Mirror Filter (QMF) Bank Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides

531 views • 38 slides

Multivariate Gaussian Mean vector: Covariance matrix: 2 1 Conditioning a Gaussian Joint

Readings: K&F: 6.1, 6.2, 6.3, 14.1, 14.2, 14.3, 14.4, Kalman Filters Gaussian MNs Graphical Models 10708 Carlos Guestrin Carnegie Mellon University December 1 st , 2008 1 Multivariate Gaussian Mean vector: Covariance matrix: 2 1

257 views • 11 slides

Quadrature Domains and Equilibrium on the Sphere Alan Legg Department of Mathematical Sciences,

Quadrature Domains and Equilibrium on the Sphere Alan Legg Department of Mathematical Sciences, Purdue Fort Wayne MWAA October 7, 2018 Alan Legg Quadrature Domains and Equilibrium on the Sphere This talk will represent work with P. Dragnev,

840 views • 81 slides

1 Gaussian Fun Facts Well add to these as we go along! First, consider a Gaussian random

Statistical Modeling and Analysis of Neural Data (NEU 560) Princeton University, Spring 2018 Jonathan Pillow Lecture 11 notes: MAP inference and regularization Thurs, 3.15 1 Gaussian Fun Facts Well add to these as we go along! First,

259 views • 4 slides

Lecture 3 Capacity of Multiuser Gaussian Channels The Gaussian uplink: 6.1 The fading

Lecture 3 Capacity of Multiuser Gaussian Channels The Gaussian uplink: 6.1 The fading Gaussian uplink: 6.3 (parts) The Gaussian downlink: 6.2 The fading Gaussian downlink: 6.4 (parts) Mikael Skoglund, Theoretical Foundations of

347 views • 10 slides

Outline n Univariate Gaussian n Multivariate Gaussian n Law of Total Probability n

Gaussians Pieter Abbeel UC Berkeley EECS Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics Outline n Univariate Gaussian n Multivariate Gaussian n Law of Total Probability n Conditioning (Bayes rule)

580 views • 10 slides

The Gaussian Distribution Chris Williams School of Informatics, University of Edinburgh October

The Gaussian Distribution Chris Williams School of Informatics, University of Edinburgh October 2007 1 / 19 Overview Probability density functions Univariate Gaussian Multivariate Gaussian Mahalanobis distance Properties of Gaussian

940 views • 19 slides