Numerical Methods in Physics
|
|
- Austin Rose
- 5 years ago
- Views:
Transcription
1 Numerical Methods i Physics Numerische Methode i der Physik, Istructor: Ass. Prof. Dr. Lilia Boeri Room: PH Tel: Address: l.boeri@tugraz.at Room: TDK Semiarraum Time: 8:30-10 a.m. Exercises: Computer Room, PH EG 004 F hpp://itp.tugraz.at/lv/boeri/num_meth/idex.html (Lecture slides, Script, Exercises, etc).
2 TOPICS (this year): Chapter 1: Itroduc8o. Chapter : Numerical methods for the solu8o of liear ihomogeeous systems. Chapter 3: Iterpola\o of poit sets. Chapter 4: Least- Squares Approxima8o. Chapter 5: Numerical solu8o of trascedetal equa8os. Chapter 6: Numerical Itegra\o. Chapter 7: Eigevalues ad Eigevectors of real matrices. Chapter 8: Numerical Methods for the solu8o of ordiary differe8al equa8os: ii8al value problems. Chapter 9: Numerical Methods for the solu\o of ordiary differe\al equa\os: margial value problems.
3 Last week(/10/013) Liear Systems: Direct (LU decomposi8o) vs idirect methods (Gauss- Seidel). Direct methods, itera8ve improvemet of the solu\o (reduce roudoff). Direct methods, special cases: tridiagoal matrices. Idirect methods: itera8ve solu\o for sparse matrices. Idirect methods: Gauss- Seidel itera8o rule. Bad matrices: defii\o ad proper\es. Gauss- Seidel itera\o for bad matrices, storig iforma\o ad efficiecy. Gauss- Seidel method with over ad uder- relaxa8o.
4 Liear Systems Methods of Solu8o (umerical): Âx = b Direct Methods: No methodological error, BUT computa\oally expesive; roudoff errors ca be large. LU decomposi8o. Itera8ve Methods: Simple algorithms; roudoff is easily cotrolled. The solu8o is approximate (truca8o). Gauss- Seidel method.
5 Direct Methods: Âx = b ˆ Ux = y Two- steps procedure: 1) LU decomposi8o (Dooli;le ad Crout): Reformula\o of the Gaussia elimia\o. A real matrix ca always be represeted as the product of two real triagular matrices L ad U, i.e. Â = ˆL ˆ U ) The two auxiliary systems are solved through back ad forward subs\tu\o: ˆ U x = y Back ˆL y = b Forward
6 Gauss- Seidel Method: itera8ve method for liear ihomogeeous sets of equa\os. Advatages: simplicity; the matrix of coefficiets is ot chaged durig the itera\o. Very efficiet for systems with a sparse matrix of coefficiets. x i (t+1) = x i (t) Δx i (t) Δx i (t) = x i (t) + 1 a ii % ' &' j=1( j i) (i =1,..., ) ( a ij x (t) j b i * )* Star\g from a ii\al vector x 0 (star%g vector) oe obtais a sequece of vectors x (t), which coverges to the exact solu\o: lim x (t) x t I prac\ce, the itera\o stops whe a give precisio is reached: xi(t)- xi(t- 1) ε For sparse matrices, oly a few terms survive i the equa\o for the Δx i s.
7 Efficiecy of direct Methods: LU GS Efficiecy of the G- S method compared to LU decomposi\o for a Laplace equa\o (sparse matrix).
8 This week(9/10/013) Least Square Approxima8o: Defii\o of the problem. Sta8s8cal distribu\o of experimetal data: ormal ad Poisso distribu\o. Sta8s8cal proper\es of the fi[g parameters. Model Fuc\os with Liear parameters. How is this implemeted i prac8ce?
9 Least Squares Approximatio Op\mal Fikg of a data set
10 Least- Squares Approxima8o: Basic problem: Data collected from experimets form a discrete set. Physical processes are typically described by mathema\cal expressios, which employ real umbers. These are easily maipulated with the stadard tools of algebra (itegrals, deriva\ves, etc). How to reduce a discrete set to a co8uous fuc8o? Iterpola8o: Fid a aaly\cal expressio (polyomial) which passes exactly through all experimetal poits. Least- squares approxima8o: Fid the best possible aaly\cal formula to approximate the behaviour of the experimetal poits.
11 Example: We wish to measure the resis\vity of a metallic coductor: The data poits (Voltage: V; Curret I) follow Ohm s law, but are affected by a experimetal error. V = R I
12 R Due to experimetal errors, the data do ot lie exactly o top of a straight lie (Ohm s law), but oly approximately.
13 Least- Squares Approxima8o: Mathema8cal Formula8o: Give a set of poits (x i /y i ), we wish to fid a curve f(x) which approximates the poits as closely as possible takig ito accout possible ucertaii\es due to measuremet errors. We also would like to be able to assig differet weights to the poits through suitable weigh\g factors. [ ] χ = g k y k f (x k ;a) k=1 mi y k χ g k > 0 f (x k ;a) a = a 1,..., a q Experimetal values Weighted error sum Weigh8g factors (sta8s8cal) Model fuc8o Model (fi[g) parameters
14 Sta8s8cs: The quality of the results of a LSQ fit depeds o the sta\s\cal quality of the data to be fiped (y k ). The sta\s\cs of the y k s is defied i terms of their expecta\o value (E k ) ad stadard devia\o (σ k ). j E k = 1 N j =1,..., N N ( j) y k j=1 Expecta8o value Number of measuremets
15 Stadard Devia8o (σ k ): measures the spread of the y k s aroud their expecta\o value E k. # % σ k = % % % $ N j=1 ( ( y j) k E k ) & ( ( N 1 ( ( ' 1/ spread
16 Gaussia ormal distribu8o: It is useful to describe experimetal aalogue measuremets. P(x;µ,σ ) = 1 # (x µ) exp% πσ $ σ & ( ' Parameters: µ Media of the distribu\o σ Stadard Devia8o σ Variace Cetral limit theorem: the mea of a radom variable is distributed approximately ormally, irrespec\ve of the form of the origial distribu\o - > very commo!!! Stadard ormal distribu8o: μ=0, σ=0.
17 Normalized Gaussia distribu8o: We will assume that our measured (y) values have a ormal distribu\o, with media = expecta\o value, ad stadard devia\o=exp. stadard devia\o. P(y E) = 1 # (y E) exp% πσ $ σ & ( ' P(y- E) describes the probability with which a specific distace (y- E) betwee the measured value y ad its expecta\o value occurs. The stadard devia\o idicates the spread with which the measured values are distributed aroud their expecta\o value E. Data- sets with a smaller σ are more accurate.
18 P(y E) = 1 # (y E) exp% πσ $ σ & ( ' y = E ±σ Iflec\o poits (P (y)=0).
19 Geometrical Meaig of the Stadard Devia8o: +σ σ +σ σ +3σ 3σ d(y E) d(y E) d(y E) 1 $ (y E) exp& πσ % σ 1 $ (y E) exp& πσ % σ 1 $ (y E) exp& πσ % σ ' ) = ( ' ) = ( ' ) = ( σ σ 3σ
20 Poisso s sta8s8cs: The Poisso s distribu\o is a discrete probability distribu8o that expresses the probability of a give umber of evets occurrig i a fixed iterval of \me ad/or space if these evets occur with a kow average rate ad idepedetly of the \me sice the last evet. Proper8es: P pois (X = k) = f (k;λ) = λ k e λ k! E(X) = λ var(x) = λ Expecta8o Value Variace I a Poisso s distribu\o the expecta\o value ad variace are equal.
21 The Poisso s distribu\o ca be applied to systems with a large umber of possible evets, each of which is rare. Example: Cou\g Experimet (radioac\ve decay): D is a digital couter; the umber of couts has a Poisso distribu\o.
22 Sta8s8cs i the LSQ method: k=1 χ = g k y k f (x k ;a) [ ] mi The weigh8g factors g k deped o the sta\s\cs of the experimetal data sets y k : g k = 1 σ k g k = 1 y k Normal distribu8o (aalogic) Poisso distribu8o (digital)
23 Normal matrix of the LSQ fit: For the weighted error sum χ we ca defie a ormal matrix as: [ N] ij = 1 χ α i α j Its iverse is the covariace matrix, which gives importat iforma\o o the sta\s\cal proper\es of the fikg parameters: C = N 1 Covariace Matrix σ ai = c ii Stadard devia8o of the a i s (fikg parameters) r ij = c ij c ii c jj Correla8o coefficiets for the a i s: ( cij 1) Measure how strogly the i- th ad j- th parameter ifluece each other.
24 Variace ad Stadard Devia8o: V = χ N q Variace q σ V = N q N q Stadard Devia8o Number of fi[g parameters Number of degrees of freedom Quality of the LSQ fit: I case of a ideal model for N>>1, V has approximately a ormal distribu8o with E=1 ad σ=σ V. If V lies sigificatly outside the iterval: [1- σ V,1+σ V ] the fit is bad!
25 Model Fuctios with Liear Parameters
26 Model fuc8os with liear parameters: A model fuc\o with liear parameters has the form: m f (x;a) = a j ϕ j (x) j=1 Here, φ j (x) are arbitrary (liearly idepedet) basis fuc8os. f(x,a) is liear i the fikg parameters, ot i x. I the special case i which oe uses as basis fuc\o a straight lie: φ j (x)=f li (x; a 1,a ) we have the liear regressio. Liear Regressio: f (x;a 1, a ) = a 1 + a x
27 Model fuc8os with liear parameters, LSQ formula: Iser\g the expressio of the model fuc\o with liear parameters i the expressio for the weighted error sum we have: # m & χ = g k % y k a j ϕ j (x k )( k=1 $ % j=1 '( Derivig with respect to the fikg parameters we get: mi χ $ m ' = g k ϕ i (x k )& y k a j ϕ j (x k )) = 0 i =1,..., m a i k=1 %& j=1 () Which ca be recast i liear set of m equa\os: m a j g k ϕ i (x k )ϕ j (x k ) j=1 k=1 = g k y k ϕ i (x k ) Âα = β k=1
28 This ihomogeeous liear set of m equa\os is (i matrix form): Âa = β Â "# α ij $ % α ij = g k ϕ i (x k ) ϕ j(x k ) k=1 β i = g k y k ϕ i (x k ) k=1 Oe ca defie a ormal matrix, give by: ˆN "# ij $ % ij = 1 χ = g k ϕ i (x k ) ϕ j(x k ) = α ij a i a j k=1
29 Stadard Devia8o of the fi[g parameters: We cosider the liear regressio formula (q=): f (x;a 1, a ) = a 1 + a x The least- squares formula gives: [ ] χ = g k y k a bx k k=1 mi a = a opt, b = b opt The ormal matrix of the problem is give by: with: α 11 =! N = α 11 α $ 1 # " α 1 α & % g k α 1 = g k x k k=1, k=1, α = g k x k k=1,
30 The op%mized parameters (a opt ad b opt ) ca be obtaied from:! N # " a b $! & = % # " β 1 β $ & % with: β 1 = g k y k, β = g k x k y k k=1 This gives: k=1 a = a opt = β 1α β α 1 D, b = b opt = β α 11 β 1 α 1 D D is the determiat of the ormal matrix: D = α 11 α α 1
31 The covariace matrix is the iverse of the ormal matrix: C = N 1 = 1 " D $ # α α 1 α 1 α 11 % ' & The diagoal elemets give the stadard devia\os of the fikg parameters: σ a = α D, σ b = α 11 D a = a opt ±σ a, b = b opt ±σ b We will ow show that the same result ca be obtaied usig the Error Propaga8o Rule (EPR).
32 Error Propaga8o Rule: The fikg parameters are fuc\o of the x, y compoets of the data poits: a = a(x 1,..., x ; y 1,..., y ) b = b(x 1,..., x ; y 1,..., y ) We ca the use the stadard formula for the error propaga8o. Give a fuc\o f(x) of variables (x 1,,x ), with errors Δx 1,, Δx the error o the fuc8o is give by: f = f (x 1,..., x ) = f (x) f (x) = f (x)± Δf (x) ( Δf (x)) = )# f & + % ( * + $ x 1 ' # (Δx 1 ) % f $ x & ( ', (Δx ). -.
33 Applyig this to the expressio for the stadard devia\os gives: σ a = σ b = l=1 l=1 (" * a $ )* # x l (" * b $ )* # x l % ' & % ' & " σ (x l ) + $ a # y l " σ (x l ) + $ b # y l % ' & % ' & + σ (y l ) -,- + σ (y l ) -,- If we assume that there is o sta\s\cal error o the x l s, ad that the y l s have a ormal distribu8o, we have: Ad thus: 1 σ (x l ) = 0 g l = σ (y l ) ( σ 1 " a % + a = * $ ' - )* g l # y l &,- ; σ b = l=1 l=1 ( * 1 )* g l " b $ # y l % ' & + -,-
34 To calculate the par8al deriva8ves we have to recall: a = a opt = β 1α β α 1 D, b = b opt = β α 11 β 1 α 1 D α 11 = g k,α 1 = g k x k k=1, k=1,,α = g k x k k=1, α ij y l = 0; β 1 = g k y k, β = g k x k y k k=1 k=1 β 1 y l = g l ; β y l = g l x l Subs\tu\g i the expressios for the stadard devia\os we get: σ a = σ b = l=1 l=1 1 " 1 g l D (α % g l α 1 g l x l ) # $ & ' 1 " 1 g l D ( α % 1g l +α 11 g l x l ) # $ & '
35 We obtai (for example for a): σ a = l=1 1 " 1 g l D (α % g l α 1 g l x l ) # $ & ' 1 ) 1 = g l D g, + l. [(α α 1 x l )] = * - l=1 = 1 g D l [(α α 1 x l )] = l=1 = 1 g D l g k g k' " # x k x k' x l x k x k' + x l x k x k' % & l,k,k' Reamig the idexes of the sum i the last two terms we obtai: σ a = 1 g D l g k g k' x k # $ l,k,k' x k' = x l x k' % & = 1 # g D k x k - k=1 $ - g l g k' l=1 k'=1 x ' g x * k' ), ( l l + l=1 %. &.
36 σ a = 1 ) g D k x k + k=1 * + g l g k' l=1 k'=1 x # g x & k' % ( $ l l ' l=1,. -. D = α 11 α α 1 α 11 = g k,α 1 = g k x k k=1, k=1,,α = g k x k k=1, We thus obtai: σ a = 1 D D g kx k = α D k=1 i.e. the diagoal terms of the covariace matrix give the stadard devia\o for the fikg parameters!!!!
37 I summary (LSQ with liear model parameters): 1) Iput the experimetal data set: x k, y k ad the sta\s\cal weights g k ) Choose a set of basis fuc8os φ j (x): 3) Costruct the auxiliary liear problem: f (x;a) = a j ϕ j (x) Âa = β Â "# α ij $ % α ij = g k ϕ i (x k ) ϕ j(x k ), β i = g k y k ϕ i (x k ) k=1 4) Solve the liear problem (LU decomposi\o); Fid op\mal fikg parameters. a opt = (a opt 1,..., a opt q ) 5) Calculate the covariace matrix (stadard devia\os of the fikg parameters). m j=1 k=1 6) Calculate the value of the op\mal fikg fuc\o o the give data poits: f (x k ;a opt ) 7) Evaluate the weighted error sum. [ ] χ = g k y k f (x k ;a) k=1
38 This week(9/10/013) Least Square Approxima8o: Defii\o of the problem. Sta8s8cal distribu\o of experimetal data: ormal ad Poisso distribu\o. Sta8s8cal proper\es of the fi[g parameters. Model Fuc\os with Liear parameters. How is this implemeted i prac8ce?
Numerical Methods in Physics
Numerical Methods in Physics Numerische Methoden in der Physik, 515.421. Instructor: Ass. Prof. Dr. Lilia Boeri Room: PH 03 090 Tel: +43-316- 873 8191 Email Address: l.boeri@tugraz.at Room: TDK Seminarraum
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationLecture 2: Poisson Sta*s*cs Probability Density Func*ons Expecta*on and Variance Es*mators
Lecture 2: Poisso Sta*s*cs Probability Desity Fuc*os Expecta*o ad Variace Es*mators Biomial Distribu*o: P (k successes i attempts) =! k!( k)! p k s( p s ) k prob of each success Poisso Distributio Note
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2016 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationLecture 3. Properties of Summary Statistics: Sampling Distribution
Lecture 3 Properties of Summary Statistics: Samplig Distributio Mai Theme How ca we use math to justify that our umerical summaries from the sample are good summaries of the populatio? Lecture Summary
More informationApproximations and more PMFs and PDFs
Approximatios ad more PMFs ad PDFs Saad Meimeh 1 Approximatio of biomial with Poisso Cosider the biomial distributio ( b(k,,p = p k (1 p k, k λ: k Assume that is large, ad p is small, but p λ at the limit.
More informationParameter, Statistic and Random Samples
Parameter, Statistic ad Radom Samples A parameter is a umber that describes the populatio. It is a fixed umber, but i practice we do ot kow its value. A statistic is a fuctio of the sample data, i.e.,
More informationLecture 3: August 31
36-705: Itermediate Statistics Fall 018 Lecturer: Siva Balakrisha Lecture 3: August 31 This lecture will be mostly a summary of other useful expoetial tail bouds We will ot prove ay of these i lecture,
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 informationProbability 2 - Notes 10. Lemma. If X is a random variable and g(x) 0 for all x in the support of f X, then P(g(X) 1) E[g(X)].
Probability 2 - Notes 0 Some Useful Iequalities. Lemma. If X is a radom variable ad g(x 0 for all x i the support of f X, the P(g(X E[g(X]. Proof. (cotiuous case P(g(X Corollaries x:g(x f X (xdx x:g(x
More informationThe variance of a sum of independent variables is the sum of their variances, since covariances are zero. Therefore. V (xi )= n n 2 σ2 = σ2.
SAMPLE STATISTICS A radom sample x 1,x,,x from a distributio f(x) is a set of idepedetly ad idetically variables with x i f(x) for all i Their joit pdf is f(x 1,x,,x )=f(x 1 )f(x ) f(x )= f(x i ) The sample
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 informationECONOMETRIC THEORY. MODULE XIII Lecture - 34 Asymptotic Theory and Stochastic Regressors
ECONOMETRIC THEORY MODULE XIII Lecture - 34 Asymptotic Theory ad Stochastic Regressors Dr. Shalabh Departmet of Mathematics ad Statistics Idia Istitute of Techology Kapur Asymptotic theory The asymptotic
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
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 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 informationBHW #13 1/ Cooper. ENGR 323 Probabilistic Analysis Beautiful Homework # 13
BHW # /5 ENGR Probabilistic Aalysis Beautiful Homework # Three differet roads feed ito a particular freeway etrace. Suppose that durig a fixed time period, the umber of cars comig from each road oto the
More informationOutline. Linear regression. Regularization functions. Polynomial curve fitting. Stochastic gradient descent for regression. MLE for regression
REGRESSION 1 Outlie Liear regressio Regularizatio fuctios Polyomial curve fittig Stochastic gradiet descet for regressio MLE for regressio Step-wise forward regressio Regressio methods Statistical techiques
More information17. Joint distributions of extreme order statistics Lehmann 5.1; Ferguson 15
17. Joit distributios of extreme order statistics Lehma 5.1; Ferguso 15 I Example 10., we derived the asymptotic distributio of the maximum from a radom sample from a uiform distributio. We did this usig
More information4. Partial Sums and the Central Limit Theorem
1 of 10 7/16/2009 6:05 AM Virtual Laboratories > 6. Radom Samples > 1 2 3 4 5 6 7 4. Partial Sums ad the Cetral Limit Theorem The cetral limit theorem ad the law of large umbers are the two fudametal theorems
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 information[ ] ( ) ( ) [ ] ( ) 1 [ ] [ ] Sums of Random Variables Y = a 1 X 1 + a 2 X 2 + +a n X n The expected value of Y is:
PROBABILITY FUNCTIONS A radom variable X has a probabilit associated with each of its possible values. The probabilit is termed a discrete probabilit if X ca assume ol discrete values, or X = x, x, x 3,,
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 19 11/17/2008 LAWS OF LARGE NUMBERS II THE STRONG LAW OF LARGE NUMBERS
MASSACHUSTTS INSTITUT OF TCHNOLOGY 6.436J/5.085J Fall 2008 Lecture 9 /7/2008 LAWS OF LARG NUMBRS II Cotets. The strog law of large umbers 2. The Cheroff boud TH STRONG LAW OF LARG NUMBRS While the weak
More informationSequences & Summa,ons
Sequeces & Summa,os Sec,o 2.4 of Rose Sprig 2012 CSCE 235 Itroduc5o to Discrete Structures Course web- page: cse.ul.edu/~cse235 Ques,os: Piazza Outlie Although you are (more or less) familiar with sequeces
More informationProperties and Hypothesis Testing
Chapter 3 Properties ad Hypothesis Testig 3.1 Types of data The regressio techiques developed i previous chapters ca be applied to three differet kids of data. 1. Cross-sectioal data. 2. Time series data.
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 informationThe Poisson Process *
OpeStax-CNX module: m11255 1 The Poisso Process * Do Johso This work is produced by OpeStax-CNX ad licesed uder the Creative Commos Attributio Licese 1.0 Some sigals have o waveform. Cosider the measuremet
More informationWEIGHTED LEAST SQUARES - used to give more emphasis to selected points in the analysis. Recall, in OLS we minimize Q =! % =!
WEIGHTED LEAST SQUARES - used to give more emphasis to selected poits i the aalysis What are eighted least squares?! " i=1 i=1 Recall, i OLS e miimize Q =! % =!(Y - " - " X ) or Q = (Y_ - X "_) (Y_ - X
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 informationExponential Families and Bayesian Inference
Computer Visio Expoetial Families ad Bayesia Iferece Lecture Expoetial Families A expoetial family of distributios is a d-parameter family f(x; havig the followig form: f(x; = h(xe g(t T (x B(, (. where
More information7.1 Convergence of sequences of random variables
Chapter 7 Limit theorems Throughout this sectio we will assume a probability space (Ω, F, P), i which is defied a ifiite sequece of radom variables (X ) ad a radom variable X. The fact that for every ifiite
More informationLecture 33: Bootstrap
Lecture 33: ootstrap Motivatio To evaluate ad compare differet estimators, we eed cosistet estimators of variaces or asymptotic variaces of estimators. This is also importat for hypothesis testig ad cofidece
More informationThis exam contains 19 pages (including this cover page) and 10 questions. A Formulae sheet is provided with the exam.
Probability ad Statistics FS 07 Secod Sessio Exam 09.0.08 Time Limit: 80 Miutes Name: Studet ID: This exam cotais 9 pages (icludig this cover page) ad 0 questios. A Formulae sheet is provided with the
More informationNonlinear regression
oliear regressio How to aalyse data? How to aalyse data? Plot! How to aalyse data? Plot! Huma brai is oe the most powerfull computatioall tools Works differetly tha a computer What if data have o liear
More information1 of 7 7/16/2009 6:06 AM Virtual Laboratories > 6. Radom Samples > 1 2 3 4 5 6 7 6. Order Statistics Defiitios Suppose agai that we have a basic radom experimet, ad that X is a real-valued radom variable
More informationRandom Matrices with Blocks of Intermediate Scale Strongly Correlated Band Matrices
Radom Matrices with Blocks of Itermediate Scale Strogly Correlated Bad Matrices Jiayi Tog Advisor: Dr. Todd Kemp May 30, 07 Departmet of Mathematics Uiversity of Califoria, Sa Diego Cotets Itroductio Notatio
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 informationLast time: Moments of the Poisson distribution from its generating function. Example: Using telescope to measure intensity of an object
6.3 Stochastic Estimatio ad Cotrol, Fall 004 Lecture 7 Last time: Momets of the Poisso distributio from its geeratig fuctio. Gs () e dg µ e ds dg µ ( s) µ ( s) µ ( s) µ e ds dg X µ ds X s dg dg + ds ds
More informationEECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1
EECS564 Estimatio, Filterig, ad Detectio Hwk 2 Sols. Witer 25 4. Let Z be a sigle observatio havig desity fuctio where. p (z) = (2z + ), z (a) Assumig that is a oradom parameter, fid ad plot the maximum
More informationHOMEWORK I: PREREQUISITES FROM MATH 727
HOMEWORK I: PREREQUISITES FROM MATH 727 Questio. Let X, X 2,... be idepedet expoetial radom variables with mea µ. (a) Show that for Z +, we have EX µ!. (b) Show that almost surely, X + + X (c) Fid the
More informationLecture 7: Properties of Random Samples
Lecture 7: Properties of Radom Samples 1 Cotiued From Last Class Theorem 1.1. Let X 1, X,...X be a radom sample from a populatio with mea µ ad variace σ
More informationProbability, Expectation Value and Uncertainty
Chapter 1 Probability, Expectatio Value ad Ucertaity We have see that the physically observable properties of a quatum system are represeted by Hermitea operators (also referred to as observables ) such
More informationLecture 1 Probability and Statistics
Wikipedia: Lecture 1 Probability ad Statistics Bejami Disraeli, British statesma ad literary figure (1804 1881): There are three kids of lies: lies, damed lies, ad statistics. popularized i US by Mark
More informationSignal Processing in Mechatronics
Sigal Processig i Mechatroics Zhu K.P. AIS, UM. Lecture, Brief itroductio to Sigals ad Systems, Review of Liear Algebra ad Sigal Processig Related Mathematics . Brief Itroductio to Sigals What is sigal
More informationThis section is optional.
4 Momet Geeratig Fuctios* This sectio is optioal. The momet geeratig fuctio g : R R of a radom variable X is defied as g(t) = E[e tx ]. Propositio 1. We have g () (0) = E[X ] for = 1, 2,... Proof. Therefore
More informationMath 61CM - Solutions to homework 3
Math 6CM - Solutios to homework 3 Cédric De Groote October 2 th, 208 Problem : Let F be a field, m 0 a fixed oegative iteger ad let V = {a 0 + a x + + a m x m a 0,, a m F} be the vector space cosistig
More informationRandom Variables, Sampling and Estimation
Chapter 1 Radom Variables, Samplig ad Estimatio 1.1 Itroductio This chapter will cover the most importat basic statistical theory you eed i order to uderstad the ecoometric material that will be comig
More informationStatistical Inference (Chapter 10) Statistical inference = learn about a population based on the information provided by a sample.
Statistical Iferece (Chapter 10) Statistical iferece = lear about a populatio based o the iformatio provided by a sample. Populatio: The set of all values of a radom variable X of iterest. Characterized
More informationLinear Regression Models
Liear Regressio Models Dr. Joh Mellor-Crummey Departmet of Computer Sciece Rice Uiversity johmc@cs.rice.edu COMP 528 Lecture 9 15 February 2005 Goals for Today Uderstad how to Use scatter diagrams to ispect
More informationChapter 2 The Monte Carlo Method
Chapter 2 The Mote Carlo Method The Mote Carlo Method stads for a broad class of computatioal algorithms that rely o radom sampligs. It is ofte used i physical ad mathematical problems ad is most useful
More informationChapter 6 Sampling Distributions
Chapter 6 Samplig Distributios 1 I most experimets, we have more tha oe measuremet for ay give variable, each measuremet beig associated with oe radomly selected a member of a populatio. Hece we eed to
More informationResampling Methods. X (1/2), i.e., Pr (X i m) = 1/2. We order the data: X (1) X (2) X (n). Define the sample median: ( n.
Jauary 1, 2019 Resamplig Methods Motivatio We have so may estimators with the property θ θ d N 0, σ 2 We ca also write θ a N θ, σ 2 /, where a meas approximately distributed as Oce we have a cosistet estimator
More informationLecture 1 Probability and Statistics
Wikipedia: Lecture 1 Probability ad Statistics Bejami Disraeli, British statesma ad literary figure (1804 1881): There are three kids of lies: lies, damed lies, ad statistics. popularized i US by Mark
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 informationMATHEMATICAL SCIENCES PAPER-II
MATHEMATICAL SCIENCES PAPER-II. Let {x } ad {y } be two sequeces of real umbers. Prove or disprove each of the statemets :. If {x y } coverges, ad if {y } is coverget, the {x } is coverget.. {x + y } coverges
More information4. Hypothesis testing (Hotelling s T 2 -statistic)
4. Hypothesis testig (Hotellig s T -statistic) Cosider the test of hypothesis H 0 : = 0 H A = 6= 0 4. The Uio-Itersectio Priciple W accept the hypothesis H 0 as valid if ad oly if H 0 (a) : a T = a T 0
More informationSince X n /n P p, we know that X n (n. Xn (n X n ) Using the asymptotic result above to obtain an approximation for fixed n, we obtain
Assigmet 9 Exercise 5.5 Let X biomial, p, where p 0, 1 is ukow. Obtai cofidece itervals for p i two differet ways: a Sice X / p d N0, p1 p], the variace of the limitig distributio depeds oly o p. Use the
More informationStatistics 511 Additional Materials
Cofidece Itervals o mu Statistics 511 Additioal Materials This topic officially moves us from probability to statistics. We begi to discuss makig ifereces about the populatio. Oe way to differetiate probability
More informationTHE KALMAN FILTER RAUL ROJAS
THE KALMAN FILTER RAUL ROJAS Abstract. This paper provides a getle itroductio to the Kalma filter, a umerical method that ca be used for sesor fusio or for calculatio of trajectories. First, we cosider
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 information1 Last time: similar and diagonalizable matrices
Last time: similar ad diagoalizable matrices Let be a positive iteger Suppose A is a matrix, v R, ad λ R Recall that v a eigevector for A with eigevalue λ if v ad Av λv, or equivaletly if v is a ozero
More informationLecture 2: Concentration Bounds
CSE 52: Desig ad Aalysis of Algorithms I Sprig 206 Lecture 2: Cocetratio Bouds Lecturer: Shaya Oveis Ghara March 30th Scribe: Syuzaa Sargsya Disclaimer: These otes have ot bee subjected to the usual scrutiy
More informationLecture 22: Review for Exam 2. 1 Basic Model Assumptions (without Gaussian Noise)
Lecture 22: Review for Exam 2 Basic Model Assumptios (without Gaussia Noise) We model oe cotiuous respose variable Y, as a liear fuctio of p umerical predictors, plus oise: Y = β 0 + β X +... β p X p +
More information7.1 Convergence of sequences of random variables
Chapter 7 Limit Theorems Throughout this sectio we will assume a probability space (, F, P), i which is defied a ifiite sequece of radom variables (X ) ad a radom variable X. The fact that for every ifiite
More information(7 One- and Two-Sample Estimation Problem )
34 Stat Lecture Notes (7 Oe- ad Two-Sample Estimatio Problem ) ( Book*: Chapter 8,pg65) Probability& Statistics for Egieers & Scietists By Walpole, Myers, Myers, Ye Estimatio 1 ) ( ˆ S P i i Poit estimate:
More informationEcon 325 Notes on Point Estimator and Confidence Interval 1 By Hiro Kasahara
Poit Estimator Eco 325 Notes o Poit Estimator ad Cofidece Iterval 1 By Hiro Kasahara Parameter, Estimator, ad Estimate The ormal probability desity fuctio is fully characterized by two costats: populatio
More informationOutput Analysis (2, Chapters 10 &11 Law)
B. Maddah ENMG 6 Simulatio Output Aalysis (, Chapters 10 &11 Law) Comparig alterative system cofiguratio Sice the output of a simulatio is radom, the comparig differet systems via simulatio should be doe
More informationLarge Sample Theory. Convergence. Central Limit Theorems Asymptotic Distribution Delta Method. Convergence in Probability Convergence in Distribution
Large Sample Theory Covergece Covergece i Probability Covergece i Distributio Cetral Limit Theorems Asymptotic Distributio Delta Method Covergece i Probability A sequece of radom scalars {z } = (z 1,z,
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 9 Multicolliearity Dr Shalabh Departmet of Mathematics ad Statistics Idia Istitute of Techology Kapur Multicolliearity diagostics A importat questio that
More informationThe Maximum-Likelihood Decoding Performance of Error-Correcting Codes
The Maximum-Lielihood Decodig Performace of Error-Correctig Codes Hery D. Pfister ECE Departmet Texas A&M Uiversity August 27th, 2007 (rev. 0) November 2st, 203 (rev. ) Performace of Codes. Notatio X,
More informationIE 230 Seat # Name < KEY > Please read these directions. Closed book and notes. 60 minutes.
IE 230 Seat # Name < KEY > Please read these directios. Closed book ad otes. 60 miutes. Covers through the ormal distributio, Sectio 4.7 of Motgomery ad Ruger, fourth editio. Cover page ad four pages of
More informationMeasurement uncertainty of the sound absorption
Measuremet ucertaity of the soud absorptio coefficiet Aa Izewska Buildig Research Istitute, Filtrowa Str., 00-6 Warsaw, Polad a.izewska@itb.pl 6887 The stadard ISO/IEC 705:005 o the competece of testig
More informationConvergence of random variables. (telegram style notes) P.J.C. Spreij
Covergece of radom variables (telegram style otes).j.c. Spreij this versio: September 6, 2005 Itroductio As we kow, radom variables are by defiitio measurable fuctios o some uderlyig measurable space
More informationLecture 19: Convergence
Lecture 19: Covergece Asymptotic approach I statistical aalysis or iferece, a key to the success of fidig a good procedure is beig able to fid some momets ad/or distributios of various statistics. I may
More informationSome Basic Probability Concepts. 2.1 Experiments, Outcomes and Random Variables
Some Basic Probability Cocepts 2. Experimets, Outcomes ad Radom Variables A radom variable is a variable whose value is ukow util it is observed. The value of a radom variable results from a experimet;
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 informationMAT1026 Calculus II Basic Convergence Tests for Series
MAT026 Calculus II Basic Covergece Tests for Series Egi MERMUT 202.03.08 Dokuz Eylül Uiversity Faculty of Sciece Departmet of Mathematics İzmir/TURKEY Cotets Mootoe Covergece Theorem 2 2 Series of Real
More informationDistribution of Random Samples & Limit theorems
STAT/MATH 395 A - PROBABILITY II UW Witer Quarter 2017 Néhémy Lim Distributio of Radom Samples & Limit theorems 1 Distributio of i.i.d. Samples Motivatig example. Assume that the goal of a study is to
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 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 informationEE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course
Sigal-EE Postal Correspodece Course 1 SAMPLE STUDY MATERIAL Electrical Egieerig EE / EEE Postal Correspodece Course GATE, IES & PSUs Sigal System Sigal-EE Postal Correspodece Course CONTENTS 1. SIGNAL
More informationMa 530 Introduction to Power Series
Ma 530 Itroductio to Power Series Please ote that there is material o power series at Visual Calculus. Some of this material was used as part of the presetatio of the topics that follow. What is a Power
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 informationResponse Variable denoted by y it is the variable that is to be predicted measure of the outcome of an experiment also called the dependent variable
Statistics Chapter 4 Correlatio ad Regressio If we have two (or more) variables we are usually iterested i the relatioship betwee the variables. Associatio betwee Variables Two variables are associated
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 informationSummary. Recap ... Last Lecture. Summary. Theorem
Last Lecture Biostatistics 602 - Statistical Iferece Lecture 23 Hyu Mi Kag April 11th, 2013 What is p-value? What is the advatage of p-value compared to hypothesis testig procedure with size α? How ca
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 informationSequences. Notation. Convergence of a Sequence
Sequeces A sequece is essetially just a list. Defiitio (Sequece of Real Numbers). A sequece of real umbers is a fuctio Z (, ) R for some real umber. Do t let the descriptio of the domai cofuse you; it
More informationExercise 4.3 Use the Continuity Theorem to prove the Cramér-Wold Theorem, Theorem. (1) φ a X(1).
Assigmet 7 Exercise 4.3 Use the Cotiuity Theorem to prove the Cramér-Wold Theorem, Theorem 4.12. Hit: a X d a X implies that φ a X (1) φ a X(1). Sketch of solutio: As we poited out i class, the oly tricky
More informationECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015
ECE 8527: Itroductio to Machie Learig ad Patter Recogitio Midterm # 1 Vaishali Ami Fall, 2015 tue39624@temple.edu Problem No. 1: Cosider a two-class discrete distributio problem: ω 1 :{[0,0], [2,0], [2,2],
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 informationEfficient GMM LECTURE 12 GMM II
DECEMBER 1 010 LECTURE 1 II Efficiet The estimator depeds o the choice of the weight matrix A. The efficiet estimator is the oe that has the smallest asymptotic variace amog all estimators defied by differet
More informationLesson 10: Limits and Continuity
www.scimsacademy.com Lesso 10: Limits ad Cotiuity SCIMS Academy 1 Limit of a fuctio The cocept of limit of a fuctio is cetral to all other cocepts i calculus (like cotiuity, derivative, defiite itegrals
More informationComparing Two Populations. Topic 15 - Two Sample Inference I. Comparing Two Means. Comparing Two Pop Means. Background Reading
Topic 15 - Two Sample Iferece I STAT 511 Professor Bruce Craig Comparig Two Populatios Research ofte ivolves the compariso of two or more samples from differet populatios Graphical summaries provide visual
More informationPH 425 Quantum Measurement and Spin Winter SPINS Lab 1
PH 425 Quatum Measuremet ad Spi Witer 23 SPIS Lab Measure the spi projectio S z alog the z-axis This is the experimet that is ready to go whe you start the program, as show below Each atom is measured
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 3 9/11/2013. Large deviations Theory. Cramér s Theorem
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/5.070J Fall 203 Lecture 3 9//203 Large deviatios Theory. Cramér s Theorem Cotet.. Cramér s Theorem. 2. Rate fuctio ad properties. 3. Chage of measure techique.
More informationThe Poisson Distribution
MATH 382 The Poisso Distributio Dr. Neal, WKU Oe of the importat distributios i probabilistic modelig is the Poisso Process X t that couts the umber of occurreces over a period of t uits of time. This
More information