u t + f(u) x = 0, (12.1) f(u) x dx = 0. u(x, t)dx = f(u(a)) f(u(b)).
|
|
- Jack Townsend
- 6 years ago
- Views:
Transcription
1 12 Fiite Volume Methos Whe solvig a PDE umerically, how o we eal with iscotiuous iitial ata? The Fiite Volume metho has particular stregth i this area. It is commoly use for hyperbolic PDEs whose solutios ca spotaeously evelop iscotiuities as they evolve i time. These solutios are ofte calle shock waves. Coservatio Laws Cosier the coservatio law u t + f(u) x =, (12.1) where u is a (spatially) oe-imesioal coserve quatity, a f(u) is the flux of u. The cotiuous itegral formulatio of (12.1) states that b t a t b a u(x, t)x + b a f(u) x x =. u(x, t)x may be thought of as the time evolutio of the total mass of u across the omai [a, b], a is epeet oly o the flux through the bouaries, sice t b a u(x, t)x = f(u(a)) f(u(b)). This fact is a importat iea utilize by fiite volume methos, which geerally cosier the evolutio of u ot at a give poit, but istea i volume-average regios. For example, let {x i } be a gri of equally space poits with spacig x, a let C i be the i-th volume (subiterval) efie by (x i 1/2, x i+1/2 ). We are itereste i the evolutio of the volume average of u over this iterval, U i = 1 x C i u(x, t )x, where {t } is the time iscretizatio. The evolutio of these volume-average quatities will epe oly o the flux through the cell eges, so that u(x, t)x = f(u(x i 1/2, t)) f(u(x i+1/2, t)). (12.2) t C i 1
2 2 Lab 12. Fiite Volume Methos U i +1 t +1 F i 1/2 F i+1/2 t U i 1 U i U i+1 Figure 12.1: A schematic of the fluxes for the fiite volume metho as iicate by (12.3). We ca the costruct a time-steppig metho where i U i x (the total mass of the system) is coserve from oe time step to the ext. Let F i 1/2 = 1 t t +1 t t +1 t f(u(x i 1/2, t)) t. The [ ] u(x, t) x = u(x, t +1 ) u(x, t ) x, t C i C i ( ) = t Fi 1/2 F i+1/2. Thus, by itegratig (12.2) i time, we may approximate the evolutio of the cell ( volume ) averages with the metho where Ui = 1 x C i u(x, t ) t. U +1 i = Ui t ( ) Fi+1/2 x F i 1/2. (12.3) This formulatio guaratees the coservatio properties that are so esirable for coservatio laws, if the time-average fluxes Fi 1/2 ca be iscretize i a atural way. The key cotributio of fiite volume methos is the computatio of Fi 1/2. For a truly oliear f(u) this ca be rather complicate a messy, a typically will ivolve solvig what is usually referre to as the Riema problem for the coservatio law. The itereste stuet ca look at [Le22] for a very thorough itrouctio a iscussio o the subject. We will cosier the liear problem i oe imesio. The aalog to higher imesios is obtaie by cosierig the eigevector ecompositio of ay liear system. Noliear equatios complicate thigs further. The liear avectio equatio a upwiig The simplest coservatio law escribes the avectio or trasport of a quatity. The PDE is give by u t + au x =, (12.4) a escribes the motio of a cocetratio of some costituet u by a costat velocity oeimesioal wi a >.
3 3 I higher imesios this is a importat problem i may fiels, for example the trasport of chemicals i the atmosphere a oceas, proper mixig of various properties i metallurgy, a the passig of iformatio alog a etwork. Note that wheever u(x, t) is a solutio of the avectio equatio, the u(x at, t ) (for ay fixe t ) is also a solutio. Thus, if u(x, ) = u (x) the the solutio for all time ca be represete by u(x, t) = u (x at). This is a importat property of (12.4), a gives a ew meaig to the term avectio: this equatio merely takes the iitial coitios a passively trasports them with velocity a. For this equatio the computatio of the flux appears straightforwar: F i 1/2 = au i 1/2 where the U i 1/2 refers to the time average of U i 1/2 over the iterval t to t +1. Let us etermie how to approximate this time average. Note from Figure 12.2 that whe a > the flux that etermies U +1 i will be epeet o the value of Ui 1. Thus, oe possibility is to approximate the flux by F i 1/2 = au i 1. Usig this approximatio of the flux together with the flux ifferecig formula (12.3) yiels the first orer upwi metho, give by U +1 i = Ui a t ( U x i Ui 1 ). Aother way to erive the upwi metho is to istea suppose that what we wat to o is recostruct u(x) at each time step isie each cell (x i 1/2, x i+1/2 )from the mea values i that cell a its surrouig eighbors. This recostructe ũ(x) is the efie piecewise for each cell i. The solutio at the ext time step ca be fou as ũ +1 (x) = ũ (x a t) which allows us to etermie the fluxes F i 1/2 oce we have settle o a metho for etermiig ũ (x) i each cell. The simplest approach is ũ (x) = U i for x (x i 1/2, x i+1/2 ) This leas to fluxes give by F i 1/2 = a t = a t t+1 t t = au i 1. The followig coe solves the problem ũ (x i 1/2, t) t, (12.5) ũ (x i 1/2 at) t, (12.6) u t + au x =, < x < 1, u(x, t) = f(x), u(, t) = u(1, t), (12.7) where a ( ) (x.3) 2 f(x) = exp + χ.5 (.6,.7) x.6 χ (.6,.7.) = 1.6 < x <.7 x.7 Notice that this PDE has perioic bouary coitios. Essetially we are evolvig the sigal arou the uit circle. This allows us to evolve the sigal much further to test our umerical methos,
4 4 Lab 12. Fiite Volume Methos sice we oly have to iscretize the iterval [, 1] istea of a much larger omai. To see how to implemet the bouary coitios, cosier a gri = x < x 1 <... < x N 1 < x N = 1 of evely space poits. Sice u(x) is perioic the u(x N ) = u(x ), so it is sufficiet to track x,..., x N 1. import umpy as p from matplotlib import pyplot as plt from math import floor ef upwi(u, a, xmi, xmax, t_fial, t): """ Solve the avectio equatio with perioic bouary coitios o the iterval [xmi, xmax] usig the upwi fiite volume scheme. Use u as the iitial coitios. a is the costat from the PDE. Use the size of u as the umber of oes i the spatial imesio. Let t be the umber of spaces i the time imesio (this is the same as the umber of steps if you o ot iclue the iitial state). Plot a show the compute solutio alog with the exact solutio. """ t = float(t_fial) / t # Sice we are oig perioic bouary coitios, # we ee to ivie by u.size istea of (u.size - 1). x = float(xmax - xmi) / u.size lamba_ = a * t / x u = u.copy() for j i rage(t): # The Upwi metho. The p.roll fuctio helps us # accout for the perioic bouary coitios. u -= lamba_ * (u - p.roll(u, 1)) # Get the x values for the plots. x = p.lispace(xmi, xmax, u.size+1)[:-1] # Plot the compute solutio. plt.plot(x, u, label='upwi Metho') # Fi the exact solutio a plot it. istace = a * t_fial roll = it((istace - floor(istace)) * u.size) plt.plot(x, p.roll(u, roll), label='exact solutio') # Show the plot with the lege. plt.lege(loc='best') plt.show() # Defie the iitial coitios. # Leave off the last poit sice we're usig perioic # bouary coitios. x = 3 t = x * 3 // 2 x = p.lispace(., 1., x+1)[:-1]
5 5 u = p.exp(-(x -.3)**2 /.5) arr = (.6 < x) & (x <.7 ) u[arr] += 1. # Ru the simulatio. upwi(u, 1.2,, 1, 1.2, t) Try ruig the previous coe block with x set to 3, 6, 12, a 24. You will otice that the umerical solutio iffuses with time. It iffuses especially fast at the poits of iscotiuity. Piecewise liear recostructio a slope limiters The upwi metho is formally oly first orer, a actually oes relatively poorly i terms of actually trasportig the iitial ata with velocity a. You ca otice from the example coe that the upwi metho has errors that are iffusive meaig that the iitial ata is iffuse as time evolves, losig the peaks a fie etails. This is because the error for the upwi metho is o the orer of the seco erivative of u which is of a iffusive ature. To get a improve metho, cosier a better recostructio isie each cell, i.e. ũ (x) = U i + m i (x x i ) for x (x i 1/2, x i+1/2 ) (12.8) where the slope of this liear recostructio m i is etermie as a fuctio of the eighborig cell averages at time a Ui itself. The the flux is give by F i 1/2 = a t = a t ( = a t+1 t t ũ (x i 1/2 at) t, Ui 1 + m i (x i 1/2 at x i ), ) Ui 1 + m i 1 ( x a t). 2 (12.9) Oe of the most atural approaches is to just estimate the slope epeig o the cell i a a eighborig cell i + 1 or i 1. This leas to two popular methos, the Lax-Weroff metho a the Beam-Warmig metho (that really is the ame). The Lax-Weroff metho has a slope chose as m i = U i+1 U i x. (12.1) which it turs out is formally seco-orer accurate. It turs out though that the errors for this metho are ispersive, meaig that ear very steep graiets, the metho will geerate very rapi oscillatios (ue to the thir erivative of u ot beig approximate accurately). Aother way to cosier how these errors arise is to otice from Figure 12.2 that if the piecewise liear recostructio is avocate by some positive wi a the there will be places where the iscotiuous ature of the recostructio will itrouce spurious maxima or miima ito the solutio. These become the spurious waves see i simulatios usig the Lax-Weroff metho. A solutio to this ilemma betwee balacig the iffusive a ispersive errors comes from costructig slopes m i that esure o such o-mootoic trasport takes place. The basic iea is to costrai the slope so that the recostructe piecewise liear fuctio ũ (x) will ot geerate
6 6 Lab 12. Fiite Volume Methos U i+1 U i+2 U i U i 2 U i 1 Figure 12.2: The piecewise liear recostructio for the upwi a Lax-Weroff methos. The soli lies represet the simplest recostructio of the cell averages leaig to the upwi metho, a the ashe lies are those whose slope is obtaie via the Lax-Weroff metho. Note that the LW metho itrouces a spurious maximum at i + 3/2 (the cell ege betwee U i+1 a U i+2 ) a the miimum at i 3/2 will be uphysical exaggerate. The upwi metho avois this ifficulties, but clearly loses a sigificat amout of the available iformatio. This provies the motivatio for the slope limiters. uphysical extremal values whe it is avocate by some fiite wi a. The Mimo limiter chooses the slope as ( U m i = mimo i Ui 1, U i+1 U ) i (12.11) x x where a if a < b a ab > mimo(a, b) = b if b < a a ab > if ab <. (12.12)
7 Aalytic solutio Upwi Lax Weroff Mi mo Figure 12.3: Solutios of (12.7) at time t = 1.2 usig various methos. Here the avectio coefficiet is a = 1.2, a there are N = 1 subitervals i space, 15 subitervals i time. Problem 1. Implemet the Lax Weroff metho a use it to solve (12.7). For N = 3, 6, 12, 24, plot the aalytic solutio a the Lax-Weroff solutio. (You shoul have 4 separate plots, each with 3 graphs.) You shoul be able to tell that the Lax Weroff metho approximates the smooth portio of the sigal much better, as it oes ot struggle with iffusio. Ufortuately, it has some ifficulty with the iscotiuous portio, where uphysical oscillatios are see. Recall that we saw somethig similar i the waves lab whe there were iscotiuous iitial coitios. Hit: Use equatios 12.9 a Problem 2. Implemet the Mimo metho a use it to solve (12.7). For N = 3, 6, 12, 24, plot the aaytic solutio a the Mimo solutio. (You shoul have 4 separate plots, each with 2 graphs.) Be sure to vectorize the mimo operatio. Hit: Use equatios 12.9 a 12.3.
Lecture #3. Math tools covered today
Toay s Program:. Review of previous lecture. QM free particle a particle i a bo. 3. Priciple of spectral ecompositio. 4. Fourth Postulate Math tools covere toay Lecture #3. Lear how to solve separable
More informationAP Calculus BC Review Chapter 12 (Sequences and Series), Part Two. n n th derivative of f at x = 5 is given by = x = approximates ( 6)
AP Calculus BC Review Chapter (Sequeces a Series), Part Two Thigs to Kow a Be Able to Do Uersta the meaig of a power series cetere at either or a arbitrary a Uersta raii a itervals of covergece, a kow
More information3. Calculus with distributions
6 RODICA D. COSTIN 3.1. Limits of istributios. 3. Calculus with istributios Defiitio 4. A sequece of istributios {u } coverges to the istributio u (all efie o the same space of test fuctios) if (φ, u )
More informationRepresenting Functions as Power Series. 3 n ...
Math Fall 7 Lab Represetig Fuctios as Power Series I. Itrouctio I sectio.8 we leare the series c c c c c... () is calle a power series. It is a uctio o whose omai is the set o all or which it coverges.
More information6.451 Principles of Digital Communication II Wednesday, March 9, 2005 MIT, Spring 2005 Handout #12. Problem Set 5 Solutions
6.51 Priciples of Digital Commuicatio II Weesay, March 9, 2005 MIT, Sprig 2005 Haout #12 Problem Set 5 Solutios Problem 5.1 (Eucliea ivisio algorithm). (a) For the set F[x] of polyomials over ay fiel F,
More informationMechatronics II Laboratory Exercise 5 Second Order Response
Mechatroics II Laboratory Exercise 5 Seco Orer Respose Theoretical Backgrou Seco orer ifferetial equatios approximate the yamic respose of may systems. The respose of a geeric seco orer system ca be see
More informationClassical Electrodynamics
A First Look at Quatum Physics Classical Electroyamics Chapter Itrouctio a Survey Classical Electroyamics Prof. Y. F. Che Cotets A First Look at Quatum Physics. Coulomb s law a electric fiel. Electric
More informationAnalytic Number Theory Solutions
Aalytic Number Theory Solutios Sea Li Corell Uiversity sl6@corell.eu Ja. 03 Itrouctio This ocumet is a work-i-progress solutio maual for Tom Apostol s Itrouctio to Aalytic Number Theory. The solutios were
More information1 Approximating Integrals using Taylor Polynomials
Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................
More informationLecture 6 Testing Nonlinear Restrictions 1. The previous lectures prepare us for the tests of nonlinear restrictions of the form:
Eco 75 Lecture 6 Testig Noliear Restrictios The previous lectures prepare us for the tests of oliear restrictios of the form: H 0 : h( 0 ) = 0 versus H : h( 0 ) 6= 0: () I this lecture, we cosier Wal,
More information1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.
.3 Covergece Theorems of Fourier Series I this sectio, we preset the covergece of Fourier series. A ifiite sum is, by defiitio, a limit of partial sums, that is, a cos( kx) b si( kx) lim a cos( kx) b si(
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 informationSome Nonlinear Equations with Double Solutions: Soliton and Chaos
Some Noliear Equatios with Double Solutios: Solito a Chaos Yi-Fag Chag Departmet of Physics, Yua Uiversity, Kumig, 659, Chia (E-mail: yifagchag@hotmail.com) Abstract The fuametal characteristics of solito
More informationThe structure of Fourier series
The structure of Fourier series Valery P Dmitriyev Lomoosov Uiversity, Russia Date: February 3, 2011) Fourier series is costructe basig o the iea to moel the elemetary oscillatio 1, +1) by the expoetial
More informationd dx where k is a spring constant
Vorlesug IX Harmoic Oscillator 1 Basic efiitios a properties a classical mechaics Oscillator is efie as a particle subject to a liear force fiel The force F ca be epresse i terms of potetial fuctio V F
More informationBENDING FREQUENCIES OF BEAMS, RODS, AND PIPES Revision S
BENDING FREQUENCIES OF BEAMS, RODS, AND PIPES Revisio S By Tom Irvie Email: tom@vibratioata.com November, Itrouctio The fuametal frequecies for typical beam cofiguratios are give i Table. Higher frequecies
More informationRIEMANN ZEROS AND AN EXPONENTIAL POTENTIAL
RIEMANN ZEROS AND AN EXPONENTIAL POTENTIAL Jose Javier Garcia Moreta Grauate stuet of Physics at the UPV/EHU (Uiversity of Basque coutry) I Soli State Physics Ares: Practicates Aa y Grijalba 5 G P.O 644
More informationFINALTERM EXAMINATION Fall 9 Calculus & Aalytical Geometry-I Questio No: ( Mars: ) - Please choose oe Let f ( x) is a fuctio such that as x approaches a real umber a, either from left or right-had-side,
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 informationNotes on iteration and Newton s method. Iteration
Notes o iteratio ad Newto s method Iteratio Iteratio meas doig somethig over ad over. I our cotet, a iteratio is a sequece of umbers, vectors, fuctios, etc. geerated by a iteratio rule of the type 1 f
More informationFourier Series and the Wave Equation
Fourier Series ad the Wave Equatio We start with the oe-dimesioal wave equatio u u =, x u(, t) = u(, t) =, ux (,) = f( x), u ( x,) = This represets a vibratig strig, where u is the displacemet of the strig
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 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 informationMATH 10550, EXAM 3 SOLUTIONS
MATH 155, EXAM 3 SOLUTIONS 1. I fidig a approximate solutio to the equatio x 3 +x 4 = usig Newto s method with iitial approximatio x 1 = 1, what is x? Solutio. Recall that x +1 = x f(x ) f (x ). Hece,
More information4.1 Sigma Notation and Riemann Sums
0 the itegral. Sigma Notatio ad Riema Sums Oe strategy for calculatig the area of a regio is to cut the regio ito simple shapes, calculate the area of each simple shape, ad the add these smaller areas
More informationSparsification using Regular and Weighted. Graphs
Sparsificatio usig Regular a Weighte 1 Graphs Aly El Gamal ECE Departmet a Cooriate Sciece Laboratory Uiversity of Illiois at Urbaa-Champaig Abstract We review the state of the art results o spectral approximatio
More informationChapter 2 Transformations and Expectations
Chapter Trasformatios a Epectatios Chapter Distributios of Fuctios of a Raom Variable Problem: Let be a raom variable with cf F ( ) If we efie ay fuctio of, say g( ) g( ) is also a raom variable whose
More informationRIEMANN ZEROS AND A EXPONENTIAL POTENTIAL
RIEMANN ZEROS AND A EXPONENTIAL POTENTIAL Jose Javier Garcia Moreta Grauate stuet of Physics at the UPV/EHU (Uiversity of Basque coutry) I Soli State Physics Ares: Practicates Aa y Grijalba 5 G P.O 644
More informationRiemann Sums y = f (x)
Riema Sums Recall that we have previously discussed the area problem I its simplest form we ca state it this way: The Area Problem Let f be a cotiuous, o-egative fuctio o the closed iterval [a, b] Fid
More informationTHE LEGENDRE POLYNOMIALS AND THEIR PROPERTIES. r If one now thinks of obtaining the potential of a distributed mass, the solution becomes-
THE LEGENDRE OLYNOMIALS AND THEIR ROERTIES The gravitatioal potetial ψ at a poit A at istace r from a poit mass locate at B ca be represete by the solutio of the Laplace equatio i spherical cooriates.
More informationMost text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t
Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said
More informationFinite Dierence Schemes
MATH-459 Numerical Methods for Coservatio Laws by Prof. Ja S. Hesthave Solutio set 2: Fiite Dierece Schemes Exercise 2. Cosistecy A method is cosistet if its local trucatio error T k satises T k (x, t)
More informationMA131 - Analysis 1. Workbook 3 Sequences II
MA3 - Aalysis Workbook 3 Sequeces II Autum 2004 Cotets 2.8 Coverget Sequeces........................ 2.9 Algebra of Limits......................... 2 2.0 Further Useful Results........................
More information6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer.
6 Itegers Modulo I Example 2.3(e), we have defied the cogruece of two itegers a,b with respect to a modulus. Let us recall that a b (mod ) meas a b. We have proved that cogruece is a equivalece relatio
More informationWe are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n
Review of Power Series, Power Series Solutios A power series i x - a is a ifiite series of the form c (x a) =c +c (x a)+(x a) +... We also call this a power series cetered at a. Ex. (x+) is cetered at
More informationKinetics of Complex Reactions
Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet
More informationTaylor expansion: Show that the TE of f(x)= sin(x) around. sin(x) = x - + 3! 5! L 7 & 8: MHD/ZAH
Taylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. A ay poit i the eighbourhood of 0, the fuctio ƒ() ca be represeted by a power series of the followig form: X 0 f(a) f() f() ( ) f( ) ( )
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 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 information(average number of points per unit length). Note that Equation (9B1) does not depend on the
EE603 Class Notes 9/25/203 Joh Stesby Appeix 9-B: Raom Poisso Poits As iscusse i Chapter, let (t,t 2 ) eote the umber of Poisso raom poits i the iterval (t, t 2 ]. The quatity (t, t 2 ) is a o-egative-iteger-value
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 informationMathematics 1 Outcome 1a. Pascall s Triangle and the Binomial Theorem (8 pers) Cumulative total = 8 periods. Lesson, Outline, Approach etc.
prouce for by Tom Strag Pascall s Triagle a the Biomial Theorem (8 pers) Mathematics 1 Outcome 1a Lesso, Outlie, Approach etc. Nelso MIA - AH M1 1 Itrouctio to Pascal s Triagle via routes alog a set of
More informationMath 257: Finite difference methods
Math 257: Fiite differece methods 1 Fiite Differeces Remember the defiitio of a derivative f f(x + ) f(x) (x) = lim 0 Also recall Taylor s formula: (1) f(x + ) = f(x) + f (x) + 2 f (x) + 3 f (3) (x) +...
More information4.3 Growth Rates of Solutions to Recurrences
4.3. GROWTH RATES OF SOLUTIONS TO RECURRENCES 81 4.3 Growth Rates of Solutios to Recurreces 4.3.1 Divide ad Coquer Algorithms Oe of the most basic ad powerful algorithmic techiques is divide ad coquer.
More informationMoment closure for biochemical networks
Momet closure for biochemical etworks João Hespaha Departmet of Electrical a Computer Egieerig Uiversity of Califoria, Sata Barbara 9-9 email: hespaha@ece.ucsb.eu Abstract Momet closure is a techique use
More informationALGEBRA HW 7 CLAY SHONKWILER
ALGEBRA HW 7 CLAY SHONKWILER Prove, or isprove a salvage: If K is a fiel, a f(x) K[x] has o roots, the K[x]/(f(x)) is a fiel. Couter-example: Cosier the fiel K = Q a the polyomial f(x) = x 4 + 3x 2 + 2.
More informationSigma notation. 2.1 Introduction
Sigma otatio. Itroductio We use sigma otatio to idicate the summatio process whe we have several (or ifiitely may) terms to add up. You may have see sigma otatio i earlier courses. It is used to idicate
More informationLecture 2: Monte Carlo Simulation
STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?
More informationSequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet
More informationCS321. Numerical Analysis and Computing
CS Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 September 8 5 What is the Root May physical system ca
More informationk=1 s k (x) (3) and that the corresponding infinite series may also converge; moreover, if it converges, then it defines a function S through its sum
0. L Hôpital s rule You alreay kow from Lecture 0 that ay sequece {s k } iuces a sequece of fiite sums {S } through S = s k, a that if s k 0 as k the {S } may coverge to the it k= S = s s s 3 s 4 = s k.
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam will cover.-.9. This sheet has three sectios. The first sectio will remid you about techiques ad formulas that you should kow. The secod gives a umber of practice questios for you
More informationL 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5!
aylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. At ay poit i the eighbourhood of =0, the fuctio ca be represeted as a power series of the followig form: X 0 f(a) f() ƒ() f()= ( ) f( ) (
More informationInfinite Sequences and Series
Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet
More informationThe Phi Power Series
The Phi Power Series I did this work i about 0 years while poderig the relatioship betwee the golde mea ad the Madelbrot set. I have fially decided to make it available from my blog at http://semresearch.wordpress.com/.
More informationPolynomial Functions and Their Graphs
Polyomial Fuctios ad Their Graphs I this sectio we begi the study of fuctios defied by polyomial expressios. Polyomial ad ratioal fuctios are the most commo fuctios used to model data, ad are used extesively
More informationTitle. Author(s)Cho, Yonggeun; Jin, Bum Ja. CitationJournal of Mathematical Analysis and Applications, 3. Issue Date Doc URL.
Title Blow-up of viscous heat-couctig compressible flow Author(s)Cho, Yoggeu; Ji, Bum Ja CitatioJoural of Mathematical Aalysis a Applicatios, 3 Issue Date 26-8-15 Doc URL http://hl.hale.et/2115/1442 Type
More informationDefinition 2 (Eigenvalue Expansion). We say a d-regular graph is a λ eigenvalue expander if
Expaer Graphs Graph Theory (Fall 011) Rutgers Uiversity Swastik Kopparty Throughout these otes G is a -regular graph 1 The Spectrum Let A G be the ajacecy matrix of G Let λ 1 λ λ be the eigevalues of A
More informationTHE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.
THE SOLUTION OF NONLINEAR EQUATIONS f( ) = 0. Noliear Equatio Solvers Bracketig. Graphical. Aalytical Ope Methods Bisectio False Positio (Regula-Falsi) Fied poit iteratio Newto Raphso Secat The root of
More informationQuadratic Functions. Before we start looking at polynomials, we should know some common terminology.
Quadratic Fuctios I this sectio we begi the study of fuctios defied by polyomial expressios. Polyomial ad ratioal fuctios are the most commo fuctios used to model data, ad are used extesively i mathematical
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationA COMPUTATIONAL STUDY UPON THE BURR 2-DIMENSIONAL DISTRIBUTION
TOME VI (year 8), FASCICULE 1, (ISSN 1584 665) A COMPUTATIONAL STUDY UPON THE BURR -DIMENSIONAL DISTRIBUTION MAKSAY Ştefa, BISTRIAN Diaa Alia Uiversity Politehica Timisoara, Faculty of Egieerig Hueoara
More informationLinear Regression Demystified
Liear Regressio Demystified Liear regressio is a importat subject i statistics. I elemetary statistics courses, formulae related to liear regressio are ofte stated without derivatio. This ote iteds to
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam 4 will cover.-., 0. ad 0.. Note that eve though. was tested i exam, questios from that sectios may also be o this exam. For practice problems o., refer to the last review. This
More informationA sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece,, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet as
More informationMA131 - Analysis 1. Workbook 2 Sequences I
MA3 - Aalysis Workbook 2 Sequeces I Autum 203 Cotets 2 Sequeces I 2. Itroductio.............................. 2.2 Icreasig ad Decreasig Sequeces................ 2 2.3 Bouded Sequeces..........................
More informationInhomogeneous Poisson process
Chapter 22 Ihomogeeous Poisso process We coclue our stuy of Poisso processes with the case of o-statioary rates. Let us cosier a arrival rate, λ(t), that with time, but oe that is still Markovia. That
More informationCastiel, Supernatural, Season 6, Episode 18
13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio
More informationSequences I. Chapter Introduction
Chapter 2 Sequeces I 2. Itroductio A sequece is a list of umbers i a defiite order so that we kow which umber is i the first place, which umber is i the secod place ad, for ay atural umber, we kow which
More informationCS537. Numerical Analysis and Computing
CS57 Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 Jauary 9 9 What is the Root May physical system ca be
More informationCALCULUS AB SECTION I, Part A Time 60 minutes Number of questions 30 A CALCULATOR MAY NOT BE USED ON THIS PART OF THE EXAM.
AP Calculus AB Portfolio Project Multiple Choice Practice Name: CALCULUS AB SECTION I, Part A Time 60 miutes Number of questios 30 A CALCULATOR MAY NOT BE USED ON THIS PART OF THE EXAM. Directios: Solve
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 informationA note on equiangular tight frames
A ote o equiagular tight frames Thomas Strohmer Departmet of Mathematics, Uiversity of Califoria, Davis, CA 9566, USA Abstract We settle a cojecture of Joseph Rees about the existece a costructio of certai
More informationRational Functions. Rational Function. Example. The degree of q x. If n d then f x is an improper rational function. x x. We have three forms
Ratioal Fuctios We have three forms R Ratioal Fuctio a a1 a p p 0 b b1 b0 q q p a a1 a0 q b b1 b0 The egree of p The egree of q is is If the f is a improper ratioal fuctio Compare forms Epae a a a b b1
More information... and realizing that as n goes to infinity the two integrals should be equal. This yields the Wallis result-
INFINITE PRODUTS Oe defies a ifiite product as- F F F... F x [ F ] Takig the atural logarithm of each side oe has- l[ F x] l F l F l F l F... So that the iitial ifiite product will coverge oly if the sum
More informationTime-Domain Representations of LTI Systems
2.1 Itroductio Objectives: 1. Impulse resposes of LTI systems 2. Liear costat-coefficiets differetial or differece equatios of LTI systems 3. Bloc diagram represetatios of LTI systems 4. State-variable
More informationMachine Learning for Data Science (CS 4786)
Machie Learig for Data Sciece CS 4786) Lecture 9: Pricipal Compoet Aalysis The text i black outlies mai ideas to retai from the lecture. The text i blue give a deeper uderstadig of how we derive or get
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 informationHOMEWORK #10 SOLUTIONS
Math 33 - Aalysis I Sprig 29 HOMEWORK # SOLUTIONS () Prove that the fuctio f(x) = x 3 is (Riema) itegrable o [, ] ad show that x 3 dx = 4. (Without usig formulae for itegratio that you leart i previous
More information( a) ( ) 1 ( ) 2 ( ) ( ) 3 3 ( ) =!
.8,.9: Taylor ad Maclauri Series.8. Although we were able to fid power series represetatios for a limited group of fuctios i the previous sectio, it is ot immediately obvious whether ay give fuctio has
More informationIP Reference guide for integer programming formulations.
IP Referece guide for iteger programmig formulatios. by James B. Orli for 15.053 ad 15.058 This documet is iteded as a compact (or relatively compact) guide to the formulatio of iteger programs. For more
More informationMath 312 Lecture Notes One Dimensional Maps
Math 312 Lecture Notes Oe Dimesioal Maps Warre Weckesser Departmet of Mathematics Colgate Uiversity 21-23 February 25 A Example We begi with the simplest model of populatio growth. Suppose, for example,
More informationDiscrete Mathematics for CS Spring 2005 Clancy/Wagner Notes 21. Some Important Distributions
CS 70 Discrete Mathematics for CS Sprig 2005 Clacy/Wager Notes 21 Some Importat Distributios Questio: A biased coi with Heads probability p is tossed repeatedly util the first Head appears. What is the
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 information18.440, March 9, Stirling s formula
Stirlig s formula 8.44, March 9, 9 The factorial fuctio! is importat i evaluatig biomial, hypergeometric, ad other probabilities. If is ot too large,! ca be computed directly, by calculators or computers.
More informationAddition: Property Name Property Description Examples. a+b = b+a. a+(b+c) = (a+b)+c
Notes for March 31 Fields: A field is a set of umbers with two (biary) operatios (usually called additio [+] ad multiplicatio [ ]) such that the followig properties hold: Additio: Name Descriptio Commutativity
More informationMath 113, Calculus II Winter 2007 Final Exam Solutions
Math, Calculus II Witer 7 Fial Exam Solutios (5 poits) Use the limit defiitio of the defiite itegral ad the sum formulas to compute x x + dx The check your aswer usig the Evaluatio Theorem Solutio: I this
More informationCAMI Education linked to CAPS: Mathematics. Grade The main topics in the FET Mathematics Curriculum NUMBER
- 1 - CAMI Eucatio like to CAPS: Grae 1 The mai topics i the FET Curriculum NUMBER TOPIC 1 Fuctios Number patters, sequeces a series 3 Fiace, growth a ecay 4 Algebra 5 Differetial Calculus 6 Probability
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 information1. Do the following sequences converge or diverge? If convergent, give the limit. Explicitly show your reasoning. 2n + 1 n ( 1) n+1.
Solutio: APPM 36 Review #3 Summer 4. Do the followig sequeces coverge or iverge? If coverget, give the limit. Eplicitly show your reasoig. a a = si b a = { } + + + 6 c a = e Solutio: a Note si a so, si
More informationADVANCED DIGITAL SIGNAL PROCESSING
ADVANCED DIGITAL SIGNAL PROCESSING PROF. S. C. CHAN (email : sccha@eee.hku.hk, Rm. CYC-702) DISCRETE-TIME SIGNALS AND SYSTEMS MULTI-DIMENSIONAL SIGNALS AND SYSTEMS RANDOM PROCESSES AND APPLICATIONS ADAPTIVE
More informationlecture 3: Interpolation Error Bounds
6 lecture 3: Iterpolatio Error Bouds.6 Covergece Theory for Polyomial Iterpolatio Iterpolatio ca be used to geerate low-degree polyomials that approimate a complicated fuctio over the iterval [a, b]. Oe
More informationME 375 FINAL EXAM Friday, May 6, 2005
ME 375 FINAL EXAM Friay, May 6, 005 Divisio: Kig 11:30 / Cuigham :30 (circle oe) Name: Istructios (1) This is a close book examiatio, but you are allowe three 8.5 11 crib sheets. () You have two hours
More informationAreas and Distances. We can easily find areas of certain geometric figures using well-known formulas:
Areas ad Distaces We ca easily fid areas of certai geometric figures usig well-kow formulas: However, it is t easy to fid the area of a regio with curved sides: METHOD: To evaluate the area of the regio
More informationIntegrable Functions. { f n } is called a determining sequence for f. If f is integrable with respect to, then f d does exist as a finite real number
MATH 532 Itegrable Fuctios Dr. Neal, WKU We ow shall defie what it meas for a measurable fuctio to be itegrable, show that all itegral properties of simple fuctios still hold, ad the give some coditios
More information! " * (x,t) " (x,t) dx =! #(x,t) dx = 1 all space
Chapter-4 Formalism 4- Schroiger Equatio Durig the early ays of i evelopmet of QM Schroiger a Heiseberg le the charge. Schroiger evelope a QM theory Schroiger Picture base o his famous equato. Heiseberg
More informationx 2 x x x x x + x x +2 x
Math 5440: Notes o particle radom walk Aaro Fogelso September 6, 005 Derivatio of the diusio equatio: Imagie that there is a distributio of particles spread alog the x-axis ad that the particles udergo
More informationNUMERICAL METHODS COURSEWORK INFORMAL NOTES ON NUMERICAL INTEGRATION COURSEWORK
NUMERICAL METHODS COURSEWORK INFORMAL NOTES ON NUMERICAL INTEGRATION COURSEWORK For this piece of coursework studets must use the methods for umerical itegratio they meet i the Numerical Methods module
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 informationOptimization Methods MIT 2.098/6.255/ Final exam
Optimizatio Methods MIT 2.098/6.255/15.093 Fial exam Date Give: December 19th, 2006 P1. [30 pts] Classify the followig statemets as true or false. All aswers must be well-justified, either through a short
More information