Polynomial Interpolation
|
|
- Lilian Cummings
- 5 years ago
- Views:
Transcription
1 Polyomial Iterpolatio 1 Iterpolatig Polyomials the iterpolatio problem 2 Lagrage Iterpolatio a basis of Lagrage polyomials 3 Neville Iterpolatio the value problem Neville s algorithm a Julia fuctio MCS 471 Lecture 7(a) Numerical Aalysis Ja Verschelde, 2 July 2018 Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
2 Polyomial Iterpolatio 1 Iterpolatig Polyomials the iterpolatio problem 2 Lagrage Iterpolatio a basis of Lagrage polyomials 3 Neville Iterpolatio the value problem Neville s algorithm a Julia fuctio Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
3 the iterpolatio problem Ofte we have data collected from some difficult fuctio f (x). With iterpolatio we ca represet the data by a polyomial. Iput: (x i, f i = f (x i )), i = 0, 1,...,, + 1 data poits, x i x j, for all i j, distict values for x. Output: p(x) a polyomial of degree at most so that for all i = 0, 1,..., : p(x i ) = f i. The polyomial p iterpolates the fuctio f (x) at the iterpolatio poits x i, i = 0, 1,...,. We say that p is the iterpolatig polyomial for the fuctio f (x) at x i. Two questios: 1 Is there a uique solutio to the iterpolatio problem? 2 How to efficietly compute the iterpolatig polyomial? Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
4 the coefficiet problem The coefficiet problem asks to compute the coefficiets of p: p(x) = c x + c 1 x c 2 x 2 + c 1 x + c 0 so that p(x i ) = f i for i = 0, 1,...,. Observe that we have + 1 data poits (x i, f i ) give o iput; ad + 1 coefficiets c i, i = 0, 1,..., to compute. Theorem If all iterpolatio poits are mutually distict: x i x j, for all i j, the the polyomial iterpolatio problem has a uique solutio. We prove this by settig up the iterpolatio coditios. Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
5 the iterpolatio coditios For p(x) = c x + + c 2 x 2 + c 1 x + c 0 the coditios p(x i ) = f i, for i = 0, 1,..., lead to a liear system of + 1 equatios: c x 0 + c 1x c 2 x c 1x 0 + c 0 = f 0 c x 1 + c 1x c 2 x c 1x 1 + c 0 = f 1 c x 2 + c 1x c 2 x c 1x 2 + c 0 = f 2. c x 1 + c 1x c 2x c 1x 1 + c 0 = f 1 c x + c 1 x c 2 x 2 + c 1 x + c 0 = f The ukows are the coefficiets c i of the polyomial p. Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
6 the liear system i matrix otatio c x0 + c 1x c 2 x0 2 + c 1x 0 + c 0 = f 0 c x1 + c 1x c 2 x1 2 + c 1x 1 + c 0 = f 1 c x2 + c 1x c 2 x2 2 + c 1x 2 + c 0 = f 2. c x 1 + c 1x c 2x c 1x 1 + c 0 = f 1 c x + c 1 x c 2 x 2 + c 1 x + c 0 = f 1 x 0 x0 2 x 1 0 x0 1 x 1 x1 2 x 1 1 x1 1 x 2 x2 2 x 1 2 x x 1 x 2 1 x x x 2 x 1 x 1 x c 0 c 1 c 2. c 1 c = f 0 f 1 f 2. f 1 f Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
7 the Vadermode matrix The liear system has a uique solutio the determiat of V (x 0, x 1, x 2,..., x 1, x ) = 1 x 0 x0 2 x 1 0 x0 1 x 1 x1 2 x 1 1 x1 1 x 2 x2 2 x 1 2 x x 1 x 2 1 x x x 2 x 1 x 1 x is differet from zero. Defiitio The matrix V (x 0, x 1, x 2,..., x 1, x ) is the Vadermode matrix for the poits x 0, x 1, x 2,..., x 1, x. Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
8 there is a uique solutio Let V = V (x 0, x 1, x 2,..., x 1, x ), the determiat of V is ( ) det(v ) = xi x j. i=0 j=i+1 Example for = 3: det(v ) = (x 0 x 1 )(x 0 x 2 )(x 0 x 3 )(x 1 x 2 )(x 1 x 3 )(x 2 x 3 ). Observe deg(det(v )) = ( + 1)/2 ad det(v ) 0 if x i x j for i j. Theorem The solutio to the iterpolatio problem is uique if x i x j for i j. Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
9 Polyomial Iterpolatio 1 Iterpolatig Polyomials the iterpolatio problem 2 Lagrage Iterpolatio a basis of Lagrage polyomials 3 Neville Iterpolatio the value problem Neville s algorithm a Julia fuctio Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
10 Lagrage iterpolatio Give are + 1 iterpolatio poits (x i, f i ), i = 0, 1,...,, where for all i j: x i x j. The Lagrage iterpolatig polyomial has the form where p(x) = l 0 (x)f 0 + l 1 (x)f (x 1 ) + + l (x)f, l i (x j ) = { 1 if i = j 0 if i j. I this form, we have that p(x i ) = f i, i = 0, 1,...,. Defiitio For + 1 mutually distict poits x i, i = 0, 1,...,, ( x xj the i-th Lagrage polyomial is l i (x) = x i x j j = 0 j i ). Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
11 the Lagrage iterpolatig polyomial The i-th Lagrage polyomial l i (x) is l i (x) = j = 0 j i ( x xj x i x j Example: = 3, i = 1: l 1 (x 0 ) = 0, l 1 (x 1 ) = 1, l 1 (x 2 ) = 0, l 1 (x 3 ) = 0. l 1 (x) = (x x 0)(x x 2 )(x x 3 ) (x 1 x 0 )(x 1 x 2 )(x 1 x 3 ). The solutio of the iterpolatio problem is uique as well, so the form for l 1 (x) is uique. The Lagrage iterpolatig polyomial is p(x) = i=0 j = 0 j i ( x xj x i x j ) f i. This is coveiet if oly the f i s chage while the x i s stays the same. ). Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
12 a exercise Exercise 1: Show that l i (x) = 1. i=0 Cosider the iterpolatio data (x i, f i = 1), i = 0, 1,..., ad remember that the iterpolatio problem has a uique solutio. Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
13 Polyomial Iterpolatio 1 Iterpolatig Polyomials the iterpolatio problem 2 Lagrage Iterpolatio a basis of Lagrage polyomials 3 Neville Iterpolatio the value problem Neville s algorithm a Julia fuctio Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
14 the value problem Iput: (x i, f i = f (x i )), i = 0, 1,...,, + 1 data poits, x i x j, for all i j, distict values for x, x is the value for some x. Output: p(x ) a the value of the iterpolatig polyomial at x. Theorem Let p i = f i, for i = 1, 2,..., ad ( x ) ( x x ) x 0 p 0,1,..., = p 0,..., 1 + p 1,..., x 0 x x x 0 the p 0,1,..., = p(x ). Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
15 Neville iterpolatio Let p i = f i, for i = 0, 1,..., ad ( ) ( ) x x x x0 p 0,1,..., (x) = p 0,..., 1 (x) + p 1,..., (x) x 0 x x x 0 the p 0,1,..., (x i ) = f i, i = 0, 1,...,. ( ) x0 x p 0,1,..., (x 0 ) = p 0,..., = f 0, aaloguous for x. x 0 x For x k, for k > 1 ad k < : ( ) ( ) xk x xk x 0 p 0,1,..., (x k ) = p 0,..., 1 (x k ) + p 1,..., (x k ) x 0 x }{{} x x 0 }{{} f k f k = x k x (x k x 0 ) f k x 0 x = f k Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
16 Polyomial Iterpolatio 1 Iterpolatig Polyomials the iterpolatio problem 2 Lagrage Iterpolatio a basis of Lagrage polyomials 3 Neville Iterpolatio the value problem Neville s algorithm a Julia fuctio Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
17 derivatio of Neville s algorithm For iterpolatio data (x 0, f 0 ), (x 1, f 1 ),..., (x, f ) ad a value x : p i = f i, for i = 0, 1,...,, p i,...,j = p(x ) is the value of the iterpolatig polyomial at x through (x i, f i ), (x i+1, f i+1 ),..., (x j, f j ). The values p i,...,j ca be orgaized i a triagular table. For example, for = 3: x 0 f 0 x 1 f 1 p 0,1 x 2 f 2 p 1,2 p 0,1,2 x 3 f 3 p 2,3 p 1,2,3 p 0,1,2,3 Observe that we may replace f 0 by p 0,1, f 1 by p 1,2, ad f 2 by p 2,3. Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
18 Neville s algorithm For iterpolatio data (x 0, f 0 ), (x 1, f 1 ),..., (x, f ) ad a value x : p i = f i, for i = 0, 1,...,, p i,...,j = p(x ) is the value of the iterpolatig polyomial at x through (x i, f i ), (x i+1, f i+1 ),..., (x j, f j ). for i = 1, 2,..., do for j = 1, 2,..., i do p i j,...,i = (x x i )p i j,...,i 1 (x x i j )p i j+1,...,i x i j x i For efficiet memory usage, relabel p i j,...,i as p[i j], p i j,...,i 1 as p[i j] ad p i j+1,...,i as p[i j + 1]. The cost of Neville s algorithm is O( 2 ). Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
19 Polyomial Iterpolatio 1 Iterpolatig Polyomials the iterpolatio problem 2 Lagrage Iterpolatio a basis of Lagrage polyomials 3 Neville Iterpolatio the value problem Neville s algorithm a Julia fuctio Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
20 a Julia fuctio fuctio eville(x::array{float64,1},f::array{float64,1},xx::float64) # # implemets the iterpolatio algorithm of Neville # # ON ENTRY : # x abscisses, give as a colum vector; # f ordiates, give as a colum vector; # xx poit where to evaluate the iterpolatig # polyomial through (x[i],f[i]). # # ON RETURN : # p last row of Neville s table where p[1] is # the value of the iterpolator at xx. # # EXAMPLE : # x = [32.0, 22.2, 41.6, 10.1, 50.5] # f = [ , , , , ] # xx = 27.5 # p = eville(x,f,xx) # Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
21 Neville s algorithm = legth(x) p = f dx = [0.0 for i=1:] for i=1: dx[i] = xx - x[i] ed for i=2: for j=2:i p[i-j+1] = (dx[i]*p[i-j+1] - dx[i-j+1]*p[i-j+2])/(x[i-j+1]-x[i]) ed ed retur p Numerical Aalysis (MCS 471) Polyomial Iterpolatio L-7(a) 2 July / 21
Find quadratic function which pass through the following points (0,1),(1,1),(2, 3)... 11
Adrew Powuk - http://www.powuk.com- Math 49 (Numerical Aalysis) Iterpolatio... 4. Polyomial iterpolatio (system of equatio)... 4.. Lier iterpolatio... 5... Fid a lie which pass through (,) (,)... 8...
More informationThe Method of Least Squares. To understand least squares fitting of data.
The Method of Least Squares KEY WORDS Curve fittig, least square GOAL To uderstad least squares fittig of data To uderstad the least squares solutio of icosistet systems of liear equatios 1 Motivatio Curve
More informationCALCULATION OF FIBONACCI VECTORS
CALCULATION OF FIBONACCI VECTORS Stuart D. Aderso Departmet of Physics, Ithaca College 953 Daby Road, Ithaca NY 14850, USA email: saderso@ithaca.edu ad Dai Novak Departmet of Mathematics, Ithaca College
More informationLecture 3: Divide and Conquer: Fast Fourier Transform
Lecture 3: Divide ad Coquer: Fast Fourier Trasform Polyomial Operatios vs. Represetatios Divide ad Coquer Algorithm Collapsig Samples / Roots of Uity FFT, IFFT, ad Polyomial Multiplicatio Polyomial operatios
More informationFFTs in Graphics and Vision. The Fast Fourier Transform
FFTs i Graphics ad Visio The Fast Fourier Trasform 1 Outlie The FFT Algorithm Applicatios i 1D Multi-Dimesioal FFTs More Applicatios Real FFTs 2 Computatioal Complexity To compute the movig dot-product
More informationA widely used display of protein shapes is based on the coordinates of the alpha carbons - - C α
Nice plottig of proteis: I A widely used display of protei shapes is based o the coordiates of the alpha carbos - - C α -s. The coordiates of the C α -s are coected by a cotiuous curve that roughly follows
More informationFor a 3 3 diagonal matrix we find. Thus e 1 is a eigenvector corresponding to eigenvalue λ = a 11. Thus matrix A has eigenvalues 2 and 3.
Closed Leotief Model Chapter 6 Eigevalues I a closed Leotief iput-output-model cosumptio ad productio coicide, i.e. V x = x = x Is this possible for the give techology matrix V? This is a special case
More informationThe Discrete Fourier Transform
The Discrete Fourier Trasform Complex Fourier Series Represetatio Recall that a Fourier series has the form a 0 + a k cos(kt) + k=1 b k si(kt) This represetatio seems a bit awkward, sice it ivolves two
More information3.2 Properties of Division 3.3 Zeros of Polynomials 3.4 Complex and Rational Zeros of Polynomials
Math 60 www.timetodare.com 3. Properties of Divisio 3.3 Zeros of Polyomials 3.4 Complex ad Ratioal Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered
More informationLECTURE 8: ORTHOGONALITY (CHAPTER 5 IN THE BOOK)
LECTURE 8: ORTHOGONALITY (CHAPTER 5 IN THE BOOK) Everythig marked by is ot required by the course syllabus I this lecture, all vector spaces is over the real umber R. All vectors i R is viewed as a colum
More informationDr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore
Lecture 0 Coversio Betwee State Space ad Trasfer Fuctio Represetatios i Liear Systems II Dr. Radhakat Padhi Asst. Professor Dept. of Aerospace Egieerig Idia Istitute of Sciece - Bagalore A Alterate First
More informationLinearly Independent Sets, Bases. Review. Remarks. A set of vectors,,, in a vector space is said to be linearly independent if the vector equation
Liearly Idepedet Sets Bases p p c c p Review { v v vp} A set of vectors i a vector space is said to be liearly idepedet if the vector equatio cv + c v + + c has oly the trivial solutio = = { v v vp} The
More informationRecursive Algorithms. Recurrences. Recursive Algorithms Analysis
Recursive Algorithms Recurreces Computer Sciece & Egieerig 35: Discrete Mathematics Christopher M Bourke cbourke@cseuledu A recursive algorithm is oe i which objects are defied i terms of other objects
More informationApply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j.
Eigevalue-Eigevector Istructor: Nam Su Wag eigemcd Ay vector i real Euclidea space of dimesio ca be uiquely epressed as a liear combiatio of liearly idepedet vectors (ie, basis) g j, j,,, α g α g α g α
More informationECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations
ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationMath 155 (Lecture 3)
Math 55 (Lecture 3) September 8, I this lecture, we ll cosider the aswer to oe of the most basic coutig problems i combiatorics Questio How may ways are there to choose a -elemet subset of the set {,,,
More informationApplications in Linear Algebra and Uses of Technology
1 TI-89: Let A 1 4 5 6 7 8 10 Applicatios i Liear Algebra ad Uses of Techology,adB 4 1 1 4 type i: [1,,;4,5,6;7,8,10] press: STO type i: A type i: [4,-1;-1,4] press: STO (1) Row Echelo Form: MATH/matrix
More informationDefinitions and Theorems. where x are the decision variables. c, b, and a are constant coefficients.
Defiitios ad Theorems Remember the scalar form of the liear programmig problem, Miimize, Subject to, f(x) = c i x i a 1i x i = b 1 a mi x i = b m x i 0 i = 1,2,, where x are the decisio variables. c, b,
More information*X203/701* X203/701. APPLIED MATHEMATICS ADVANCED HIGHER Numerical Analysis. Read carefully
X0/70 NATIONAL QUALIFICATIONS 006 MONDAY, MAY.00 PM.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.
More informationACO Comprehensive Exam 9 October 2007 Student code A. 1. Graph Theory
1. Graph Theory Prove that there exist o simple plaar triagulatio T ad two distict adjacet vertices x, y V (T ) such that x ad y are the oly vertices of T of odd degree. Do ot use the Four-Color Theorem.
More information(3) If you replace row i of A by its sum with a multiple of another row, then the determinant is unchanged! Expand across the i th row:
Math 5-4 Tue Feb 4 Cotiue with sectio 36 Determiats The effective way to compute determiats for larger-sized matrices without lots of zeroes is to ot use the defiitio, but rather to use the followig facts,
More informationSPLINE INTERPOLATION AND BEST QUADRATURE FORMULAE. 1. The spline interpolation formula, A spline function S(x), of degree &( = 0)> having the knots
SPLINE INTERPOLATION AND BEST QUADRATURE FORMULAE BY I. J. SCHOENBERG Commuicated by Felix Browder, October 9, 1963 1. The splie iterpolatio formula, A splie fuctio S(x), of degree &( = 0)> havig the kots
More informationBernoulli, Ramanujan, Toeplitz e le matrici triangolari
Due Giori di Algebra Lieare Numerica www.dima.uige.it/ dibeede/gg/home.html Geova, 6 7 Febbraio Beroulli, Ramauja, Toeplitz e le matrici triagolari Carmie Di Fiore, Fracesco Tudisco, Paolo Zellii Speaker:
More informationZeros of Polynomials
Math 160 www.timetodare.com 4.5 4.6 Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered with fidig the solutios of polyomial equatios of ay degree
More informationSection A assesses the Units Numerical Analysis 1 and 2 Section B assesses the Unit Mathematics for Applied Mathematics
X0/70 NATIONAL QUALIFICATIONS 005 MONDAY, MAY.00 PM 4.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.
More informationApplication of Jordan Canonical Form
CHAPTER 6 Applicatio of Jorda Caoical Form Notatios R is the set of real umbers C is the set of complex umbers Q is the set of ratioal umbers Z is the set of itegers N is the set of o-egative itegers Z
More informationPolynomial identity testing and global minimum cut
CHAPTER 6 Polyomial idetity testig ad global miimum cut I this lecture we will cosider two further problems that ca be solved usig probabilistic algorithms. I the first half, we will cosider the problem
More informationLecture 19. Curve fitting I. 1 Introduction. 2 Fitting a constant to measured data
Lecture 9 Curve fittig I Itroductio Suppose we are preseted with eight poits of easured data (x i, y j ). As show i Fig. o the left, we could represet the uderlyig fuctio of which these data are saples
More informationCALCULATING FIBONACCI VECTORS
THE GENERALIZED BINET FORMULA FOR CALCULATING FIBONACCI VECTORS Stuart D Aderso Departmet of Physics, Ithaca College 953 Daby Road, Ithaca NY 14850, USA email: saderso@ithacaedu ad Dai Novak Departmet
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More informationA collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation
Iteratioal Joural of Mathematics Research. ISSN 0976-5840 Volume 9 Number 1 (017) pp. 45-51 Iteratioal Research Publicatio House http://www.irphouse.com A collocatio method for sigular itegral equatios
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationIterative Techniques for Solving Ax b -(3.8). Assume that the system has a unique solution. Let x be the solution. Then x A 1 b.
Iterative Techiques for Solvig Ax b -(8) Cosider solvig liear systems of them form: Ax b where A a ij, x x i, b b i Assume that the system has a uique solutio Let x be the solutio The x A b Jacobi ad Gauss-Seidel
More informationCHAPTER 5. Theory and Solution Using Matrix Techniques
A SERIES OF CLASS NOTES FOR 2005-2006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 3 A COLLECTION OF HANDOUTS ON SYSTEMS OF ORDINARY DIFFERENTIAL
More informationLogit regression Logit regression
Logit regressio Logit regressio models the probability of Y= as the cumulative stadard logistic distributio fuctio, evaluated at z = β 0 + β X: Pr(Y = X) = F(β 0 + β X) F is the cumulative logistic distributio
More informationMEI Casio Tasks for Further Pure
Task Complex Numbers: Roots of Quadratic Equatios. Add a ew Equatio scree: paf 2. Chage the Complex output to a+bi: LpNNNNwd 3. Select Polyomial ad set the Degree to 2: wq 4. Set a=, b=5 ad c=6: l5l6l
More informationChapter 10: Power Series
Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because
More informationFROM GENERALIZED CAUCHY-RIEMANN EQUATIONS TO LINEAR ALGEBRAS. (1) Ê dkij ^ = 0 (* = 1, 2,, (n2- «)),
FROM GENERALIZED CAUCHY-RIEMANN EQUATIONS TO LINEAR ALGEBRAS JAMES A. WARD I a previous paper [l] the author gave a defiitio of aalytic fuctio i liear associative algebras with a idetity. With each such
More informationLøsningsførslag i 4M
Norges tekisk aturviteskapelige uiversitet Istitutt for matematiske fag Side 1 av 6 Løsigsførslag i 4M Oppgave 1 a) A sketch of the graph of the give f o the iterval [ 3, 3) is as follows: The Fourier
More informationLet A(x) and B(x) be two polynomials of degree n 1:
MI-EVY (2011/2012) J. Holub: 4. DFT, FFT ad Patter Matchig p. 2/42 Operatios o polyomials MI-EVY (2011/2012) J. Holub: 4. DFT, FFT ad Patter Matchig p. 4/42 Efficiet Patter Matchig (MI-EVY) 4. DFT, FFT
More informationTheorem: Let A n n. In this case that A does reduce to I, we search for A 1 as the solution matrix X to the matrix equation A X = I i.e.
Theorem: Let A be a square matrix The A has a iverse matrix if ad oly if its reduced row echelo form is the idetity I this case the algorithm illustrated o the previous page will always yield the iverse
More informationLecture 12: September 27
36-705: Itermediate Statistics Fall 207 Lecturer: Siva Balakrisha Lecture 2: September 27 Today we will discuss sufficiecy i more detail ad the begi to discuss some geeral strategies for costructig estimators.
More informationSolution of Linear Constant-Coefficient Difference Equations
ECE 38-9 Solutio of Liear Costat-Coefficiet Differece Equatios Z. Aliyazicioglu Electrical ad Computer Egieerig Departmet Cal Poly Pomoa Solutio of Liear Costat-Coefficiet Differece Equatios Example: Determie
More informationSignal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform
Sigal Processig i Mechatroics Summer semester, 1 Lecture 3, Covolutio, Fourier Series ad Fourier rasform Dr. Zhu K.P. AIS, UM 1 1. Covolutio Covolutio Descriptio of LI Systems he mai premise is that the
More informationCHAPTER I: Vector Spaces
CHAPTER I: Vector Spaces Sectio 1: Itroductio ad Examples This first chapter is largely a review of topics you probably saw i your liear algebra course. So why cover it? (1) Not everyoe remembers everythig
More informationLinear regression. Daniel Hsu (COMS 4771) (y i x T i β)2 2πσ. 2 2σ 2. 1 n. (x T i β y i ) 2. 1 ˆβ arg min. β R n d
Liear regressio Daiel Hsu (COMS 477) Maximum likelihood estimatio Oe of the simplest liear regressio models is the followig: (X, Y ),..., (X, Y ), (X, Y ) are iid radom pairs takig values i R d R, ad Y
More informationAlgorithm Analysis. Algorithms that are equally correct can vary in their utilization of computational resources
Algorithm Aalysis Algorithms that are equally correct ca vary i their utilizatio of computatioal resources time ad memory a slow program it is likely ot to be used a program that demads too much memory
More informationAdvanced Course of Algorithm Design and Analysis
Differet complexity measures Advaced Course of Algorithm Desig ad Aalysis Asymptotic complexity Big-Oh otatio Properties of O otatio Aalysis of simple algorithms A algorithm may may have differet executio
More informationLecture 12: November 13, 2018
Mathematical Toolkit Autum 2018 Lecturer: Madhur Tulsiai Lecture 12: November 13, 2018 1 Radomized polyomial idetity testig We will use our kowledge of coditioal probability to prove the followig lemma,
More informationThe picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled
1 Lecture : Area Area ad distace traveled Approximatig area by rectagles Summatio The area uder a parabola 1.1 Area ad distace Suppose we have the followig iformatio about the velocity of a particle, how
More informationSummary: CORRELATION & LINEAR REGRESSION. GC. Students are advised to refer to lecture notes for the GC operations to obtain scatter diagram.
Key Cocepts: 1) Sketchig of scatter diagram The scatter diagram of bivariate (i.e. cotaiig two variables) data ca be easily obtaied usig GC. Studets are advised to refer to lecture otes for the GC operatios
More informationRecurrence Relations
Recurrece Relatios Aalysis of recursive algorithms, such as: it factorial (it ) { if (==0) retur ; else retur ( * factorial(-)); } Let t be the umber of multiplicatios eeded to calculate factorial(). The
More informationECON 3150/4150, Spring term Lecture 3
Itroductio Fidig the best fit by regressio Residuals ad R-sq Regressio ad causality Summary ad ext step ECON 3150/4150, Sprig term 2014. Lecture 3 Ragar Nymoe Uiversity of Oslo 21 Jauary 2014 1 / 30 Itroductio
More informationx x x 2x x N ( ) p NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS By Newton-Raphson formula
NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS. If g( is cotiuous i [a,b], te uder wat coditio te iterative (or iteratio metod = g( as a uique solutio i [a,b]? '( i [a,b].. Wat
More informationDisjoint set (Union-Find)
CS124 Lecture 7 Fall 2018 Disjoit set (Uio-Fid) For Kruskal s algorithm for the miimum spaig tree problem, we foud that we eeded a data structure for maitaiig a collectio of disjoit sets. That is, we eed
More informationGoodness-of-Fit Tests and Categorical Data Analysis (Devore Chapter Fourteen)
Goodess-of-Fit Tests ad Categorical Data Aalysis (Devore Chapter Fourtee) MATH-252-01: Probability ad Statistics II Sprig 2019 Cotets 1 Chi-Squared Tests with Kow Probabilities 1 1.1 Chi-Squared Testig................
More information2 Geometric interpretation of complex numbers
2 Geometric iterpretatio of complex umbers 2.1 Defiitio I will start fially with a precise defiitio, assumig that such mathematical object as vector space R 2 is well familiar to the studets. Recall that
More informationA brief introduction to linear algebra
CHAPTER 6 A brief itroductio to liear algebra 1. Vector spaces ad liear maps I what follows, fix K 2{Q, R, C}. More geerally, K ca be ay field. 1.1. Vector spaces. Motivated by our ituitio of addig ad
More informationPolynomial Multiplication and Fast Fourier Transform
Polyomial Multiplicatio ad Fast Fourier Trasform Com S 477/577 Notes Ya-Bi Jia Sep 19, 2017 I this lecture we will describe the famous algorithm of fast Fourier trasform FFT, which has revolutioized digital
More informationMon Feb matrix inverses. Announcements: Warm-up Exercise:
Math 225-4 Week 6 otes We will ot ecessarily fiish the material from a give day's otes o that day We may also add or subtract some material as the week progresses, but these otes represet a i-depth outlie
More informationA Pseudo Spline Methods for Solving an Initial Value Problem of Ordinary Differential Equation
Joural of Matematics ad Statistics 4 (: 7-, 008 ISSN 549-3644 008 Sciece Publicatios A Pseudo Splie Metods for Solvig a Iitial Value Problem of Ordiary Differetial Equatio B.S. Ogudare ad G.E. Okeca Departmet
More informationIntensive Algorithms Lecture 11. DFT and DP. Lecturer: Daniel A. Spielman February 20, f(n) O(g(n) log c g(n)).
Itesive Algorithms Lecture 11 DFT ad DP Lecturer: Daiel A. Spielma February 20, 2018 11.1 Itroductio The purpose of this lecture is to lear how use the Discrete Fourier Trasform to save space i Dyamic
More information(bilinearity), a(u, v) M u V v V (continuity), a(v, v) m v 2 V (coercivity).
Precoditioed fiite elemets method Let V be a Hilbert space, (, ) V a ier product o V ad V the correspodig iduced orm. Let a be a coercive, cotiuous, biliear form o V, that is, a : V V R ad there exist
More informationEngineering Analysis ( & ) Lec(7) CH 2 Higher Order Linear ODEs
Philadelphia Uiversit/Facult of Egieerig Commuicatio ad Electroics Egieerig Egieerig Aalsis (6500 & 6300) Higher Order Liear ODEs Istructor: Eg. Nada khatib Email: khatib@philadelphia.edu.jo Higher order
More informationM A T H F A L L CORRECTION. Algebra I 1 4 / 1 0 / U N I V E R S I T Y O F T O R O N T O
M A T H 2 4 0 F A L L 2 0 1 4 HOMEWORK ASSIGNMENT #4 CORRECTION Algebra I 1 4 / 1 0 / 2 0 1 4 U N I V E R S I T Y O F T O R O N T O P r o f e s s o r : D r o r B a r - N a t a Correctio Homework Assigmet
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationChimica Inorganica 3
himica Iorgaica Irreducible Represetatios ad haracter Tables Rather tha usig geometrical operatios, it is ofte much more coveiet to employ a ew set of group elemets which are matrices ad to make the rule
More informationSeptember 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1
September 0 s (Edecel) Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright
More informationPhysics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.
Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;
More informationLinear Classifiers III
Uiversität Potsdam Istitut für Iformatik Lehrstuhl Maschielles Lere Liear Classifiers III Blaie Nelso, Tobias Scheffer Cotets Classificatio Problem Bayesia Classifier Decisio Liear Classifiers, MAP Models
More informationGlobal polynomial interpolants can have extremely undesirable properties if the sample points are not chosen correctly.
Ruge Pheomeo Global polyomial iterpolats ca have extremely udesirable properties if the sample poits are ot chose correctly. Example: Iterpolate the fuctio f(x) = over the iterval [, ] at the equally spaced
More information10-701/ Machine Learning Mid-term Exam Solution
0-70/5-78 Machie Learig Mid-term Exam Solutio Your Name: Your Adrew ID: True or False (Give oe setece explaatio) (20%). (F) For a cotiuous radom variable x ad its probability distributio fuctio p(x), it
More informationPAPER : IIT-JAM 2010
MATHEMATICS-MA (CODE A) Q.-Q.5: Oly oe optio is correct for each questio. Each questio carries (+6) marks for correct aswer ad ( ) marks for icorrect aswer.. Which of the followig coditios does NOT esure
More information15.081J/6.251J Introduction to Mathematical Programming. Lecture 21: Primal Barrier Interior Point Algorithm
508J/65J Itroductio to Mathematical Programmig Lecture : Primal Barrier Iterior Poit Algorithm Outlie Barrier Methods Slide The Cetral Path 3 Approximatig the Cetral Path 4 The Primal Barrier Algorithm
More informationWe will conclude the chapter with the study a few methods and techniques which are useful
Chapter : Coordiate geometry: I this chapter we will lear about the mai priciples of graphig i a dimesioal (D) Cartesia system of coordiates. We will focus o drawig lies ad the characteristics of the graphs
More informationR is a scalar defined as follows:
Math 8. Notes o Dot Product, Cross Product, Plaes, Area, ad Volumes This lecture focuses primarily o the dot product ad its may applicatios, especially i the measuremet of agles ad scalar projectio ad
More informationFall 2011, EE123 Digital Signal Processing
Lecture 5 Miki Lustig, UCB September 14, 211 Miki Lustig, UCB Motivatios for Discrete Fourier Trasform Sampled represetatio i time ad frequecy umerical Fourier aalysis requires a Fourier represetatio that
More informationSolutions. Number of Problems: 4. None. Use only the prepared sheets for your solutions. Additional paper is available from the supervisors.
Quiz November 4th, 23 Sigals & Systems (5-575-) P. Reist & Prof. R. D Adrea Solutios Exam Duratio: 4 miutes Number of Problems: 4 Permitted aids: Noe. Use oly the prepared sheets for your solutios. Additioal
More informationTMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods
TMA4205 Numerical Liear Algebra The Poisso problem i R 2 : diagoalizatio methods September 3, 2007 c Eiar M Røquist Departmet of Mathematical Scieces NTNU, N-749 Trodheim, Norway All rights reserved A
More informationChapter 4 : Laplace Transform
4. Itroductio Laplace trasform is a alterative to solve the differetial equatio by the complex frequecy domai ( s = σ + jω), istead of the usual time domai. The DE ca be easily trasformed ito a algebraic
More informationINTRODUCTION TO MATRIX ALGEBRA. a 11 a a 1n a 21 a a 2n...
INTRODUCTION TO MATRIX ALGEBRA DEFINITIONOFAMATRIXANDAVECTOR DefiitioofamatrixAmatrixisarectagulararrayofumbersarrageditorowsad colums It is writte as a a 2 a a 2 a 22 a 2 () a m a m2 a m Theabovearrayiscalledamby(m
More informationInverse Matrix. A meaning that matrix B is an inverse of matrix A.
Iverse Matrix Two square matrices A ad B of dimesios are called iverses to oe aother if the followig holds, AB BA I (11) The otio is dual but we ofte write 1 B A meaig that matrix B is a iverse of matrix
More informationEstimation for Complete Data
Estimatio for Complete Data complete data: there is o loss of iformatio durig study. complete idividual complete data= grouped data A complete idividual data is the oe i which the complete iformatio of
More informationMachine Learning Regression I Hamid R. Rabiee [Slides are based on Bishop Book] Spring
Machie Learig Regressio I Hamid R. Rabiee [Slides are based o Bishop Book] Sprig 015 http://ce.sharif.edu/courses/93-94//ce717-1 Liear Regressio Liear regressio: ivolves a respose variable ad a sigle predictor
More informationTable 12.1: Contingency table. Feature b. 1 N 11 N 12 N 1b 2 N 21 N 22 N 2b. ... a N a1 N a2 N ab
Sectio 12 Tests of idepedece ad homogeeity I this lecture we will cosider a situatio whe our observatios are classified by two differet features ad we would like to test if these features are idepedet
More informationDeterminants of order 2 and 3 were defined in Chapter 2 by the formulae (5.1)
5. Determiats 5.. Itroductio 5.2. Motivatio for the Choice of Axioms for a Determiat Fuctios 5.3. A Set of Axioms for a Determiat Fuctio 5.4. The Determiat of a Diagoal Matrix 5.5. The Determiat of a Upper
More informationDivide & Conquer. Divide-and-conquer algorithms. Conventional product of polynomials. Conventional product of polynomials.
Divide-ad-coquer algorithms Divide & Coquer Strategy: Divide the problem ito smaller subproblems of the same type of problem Solve the subproblems recursively Combie the aswers to solve the origial problem
More informationCS583 Lecture 02. Jana Kosecka. some materials here are based on E. Demaine, D. Luebke slides
CS583 Lecture 02 Jaa Kosecka some materials here are based o E. Demaie, D. Luebke slides Previously Sample algorithms Exact ruig time, pseudo-code Approximate ruig time Worst case aalysis Best case aalysis
More informationx c the remainder is Pc ().
Algebra, Polyomial ad Ratioal Fuctios Page 1 K.Paulk Notes Chapter 3, Sectio 3.1 to 3.4 Summary Sectio Theorem Notes 3.1 Zeros of a Fuctio Set the fuctio to zero ad solve for x. The fuctio is zero at these
More informationMatrix Representation of Data in Experiment
Matrix Represetatio of Data i Experimet Cosider a very simple model for resposes y ij : y ij i ij, i 1,; j 1,,..., (ote that for simplicity we are assumig the two () groups are of equal sample size ) Y
More informationSignal Processing. Lecture 02: Discrete Time Signals and Systems. Ahmet Taha Koru, Ph. D. Yildiz Technical University.
Sigal Processig Lecture 02: Discrete Time Sigals ad Systems Ahmet Taha Koru, Ph. D. Yildiz Techical Uiversity 2017-2018 Fall ATK (YTU) Sigal Processig 2017-2018 Fall 1 / 51 Discrete Time Sigals Discrete
More informationME NUMERICAL METHODS Fall 2007
ME - 310 NUMERICAL METHODS Fall 2007 Group 02 Istructor: Prof. Dr. Eres Söylemez (Rm C205, email:eres@metu.edu.tr ) Class Hours ad Room: Moday 13:40-15:30 Rm: B101 Wedesday 12:40-13:30 Rm: B103 Course
More informationPolynomials. Computer Programming for Engineers (2014 Spring)
Computer Programmig for Egieers (014 Sprig) Polyomials Hyougshick Kim Departmet of Computer Sciece ad Egieerig College of Iformatio ad Commuicatio Egieerig Sugkyukwa Uiversity Polyomials A th order polyomial
More informationTHE REPRESENTATION OF THE REMAINDER IN CLASSICAL BERNSTEIN APPROXIMATION FORMULA
Global Joural of Advaced Research o Classical ad Moder Geometries ISSN: 2284-5569, Vol.6, 2017, Issue 2, pp.119-125 THE REPRESENTATION OF THE REMAINDER IN CLASSICAL BERNSTEIN APPROXIMATION FORMULA DAN
More informationMachine Learning for Data Science (CS 4786)
Machie Learig for Data Sciece CS 4786) Lecture & 3: Pricipal Compoet Aalysis The text i black outlies high level ideas. The text i blue provides simple mathematical details to derive or get to the algorithm
More informationCLOSED FORM FORMULA FOR THE NUMBER OF RESTRICTED COMPOSITIONS
Submitted to the Bulleti of the Australia Mathematical Society doi:10.1017/s... CLOSED FORM FORMULA FOR THE NUMBER OF RESTRICTED COMPOSITIONS GAŠPER JAKLIČ, VITO VITRIH ad EMIL ŽAGAR Abstract I this paper,
More informationUnit 4: Polynomial and Rational Functions
48 Uit 4: Polyomial ad Ratioal Fuctios Polyomial Fuctios A polyomial fuctio y px ( ) is a fuctio of the form p( x) ax + a x + a x +... + ax + ax+ a 1 1 1 0 where a, a 1,..., a, a1, a0are real costats ad
More informationTAMS24: Notations and Formulas
TAMS4: Notatios ad Formulas Basic otatios ad defiitios X: radom variable stokastiska variabel Mea Vätevärde: µ = X = by Xiagfeg Yag kpx k, if X is discrete, xf Xxdx, if X is cotiuous Variace Varias: =
More information