Stability For a stable numerical scheme, the errors in the initial condition will not grow unboundedly with time.
|
|
- Ruth Sherman
- 5 years ago
- Views:
Transcription
1 .3.5. Stablty Aalyss Readg: Taehll et al. Secto 3.6. Stablty For a stable umercal scheme, the errors the tal codto wll ot grow uboudedly wth tme. I ths secto, we dscuss the methods for determg the stablty of F.D. schemes. Ths s very mportat whe desgg a F.D. scheme ad for uderstadg t behavor. There are several methods: Eergy method vo Neuma method Matrx method (for systems of equatos) Dscrete perturbato method (wll ot dscuss) Note: Stablty refers to the F.D. Does ot volve B.C. or I.C. Refers to tme-matchg problems oly The eergy method Read Durra Secto. (hadout, see web l)
2 Ths method s used much less ofte tha the vo Neuma method. It's attractve because t wors for olear problem ad problems wthout perod B.C. The ey s to show that a postve defte quatty le ( u ) s bouded for all. u u We llustrate ths method usg the upstream-forward scheme for wave (advecto) equato + c = 0 : t x u u u u + c = t x + 0 Let µ = c t/ x u = ( µ ) u + µ u + Squarg both sdes ad summg over all grd pots: ( u ) = [( µ ) ( u ) + ( µ ) µ u u + µ ( u ) ] () + Assumg perodc B.C. ( u ) = ( u )
3 ad usg the Schwarz equalty (whch says that for two vectors U ad V, U V U V ), uu ( u) ( u ) = ( u). If µ ( µ ) 0, all coeffcets of RHS terms () are postve, ad we have ( u ) [( µ ) + ( µ ) µ + µ ] ( u ) = ( u ), +.e., the L orm at + s o greater tha that at, therefore the scheme s stable! The codto µ ( µ ) 0 gves µ = c t/ x whch s the stablty codto for ths scheme. As we dscussed earler, t says that waves caot propagate more tha oe grd terval durg oe t order to mata stablty. vo Neuma method Read Taehll et al, Secto I a sese, we have already used ths method to fd stablty of the -D dffuso equato whe we were provg the covergece of the soluto usg FTCS scheme. We foud the 3
4 u + = [ M(t) ] + u 0 where M(t) = amplfcato factor. If M(t) by some measure (M ca be a matrx or a complex umber), the u + u ad the soluto caot grow tme the scheme s stable. Essetally, vo Neuma method expads the F.D.E. a Fourer seres, fds the amplfcato factor ad determes uder what codto the factor s less tha or equal to for stablty. Assumptos:. The equato has to be Lear wth costat coeffcets.. It s assumed that the soluto s perodc. Wth ths method, the depedet varable s decomposed to a complex (or a real) Fourer seres: uxyzt (,,, ) = Ulm (,, )exp[ x ( + ly+ mz ωt)] () lm,, where U s the complex ampltude ad ω = ω R + ω Ι s the complex frequecy. I fact ω R gves the wave propagato speed ad ω Ι gves the growth ad decayg rate. e = e = e e ωt [ ω + ω ] t ω t ω t R I R I Rt e ω I e ω t - phase fucto of Fourer compoets - growth or decay rate 4
5 If ω I >0, the soluto wll grow expoetally tme. Example : -D dffuso equato wth FTCS scheme. u u = µ ( u u + u ) (3) + + K t where µ =. ( x) Let's exame a sgle wave : u = U e ( x ωt) + ( x ωt) ω t = = ( x ωt) ± x ± u U e e u U e e Substtute the above to (3) Ue ( e ) = µ Ue [ e + e ] ( x ωt) ω t ( x ωt) x x ( x ωt ) ω t Ue [ e µ (cos x )] = 0 [ + 4µ s ( /)] = 0 ( x ωt) ω t Ue e x 5
6 For o-trval soluto, we requre e µ ω t + 4 s ( x/) = 0 e µ ω t = 4 s ( x/) Here, t e ω s actually the amplfcato factor, the same as the M dscussed earler. λ t e ω = u + /u - the amplfcato factor. For stablty we requre λ µ 4 s ( x/) µ / as before! I practce, whe µ=/, the soluto (amplfcato factor) swtches betwee ad + for x waves ( s ( x/ ) = ) very other step, whch s urealstc. The stadard requremet s therefore µ /4. Therefore, do ot avely th dffuso terms a umercal model does ot cause umercal stablty! Whe tegrated stably, the dffuso term CFD models teds to stablze the soluto by llg off/dampg small scale waves, but whe stablty codto s ot met, t tem tself wll cause problem! 6
7 Read Pele (984) secto 0.. (hadout, see web l) Implct Methods Read Taehll et al, secod part of secto So far, we have dealt wth oly explct schemes whch have the form of + u = f( u, u,...). Wth these schemes, the future state at each grd pot s oly depedet o the curret ad past tme levels, therefore the soluto ca be obtaed drectly or explctly. c.f., explct fuctos such as y = x. Implct scheme volves varables of the future tme level at more tha oe grd pot (ofte resultg from fte dfferece of varable(s) at the future tme level). Mathematcally t ca be expressed as + + u = f( u, u, u,...). Ths s aalogous to mplct fuctos such as x = s(x). As oe ca mage, mplct schemes are more dffcult to solve. Usually matrx verso s volved. We wll frst loo at the stablty property of a mplct scheme. Example. Cosder the -D dffuso equato u t = K u xx aga. It s approxmated by the followg F.D. scheme: 7
8 + + tu = K[ xxu + ( ) xxu ] (4) δ αδ α δ (Note: shorthad otatos for F.D. are used. See Appedx. e.g., δ + ( )/ u u u t + t = ). Whe α = 0, the scheme s explct, ad s the FTCS scheme dscussed earler. Whe α = /, t s mplct ad s called Cra-Ncolso scheme. For other values of α, t s a geeral mplct scheme. We ca show that 4 u ( x) τ = K K t α O x t x (show f for yourself!). 4 ( ) We ca see that whe α = /, t s d-order accurate tme ad space. Otherwse, t's frst-order tme whch s expected for u-cetered tme-dfferecg scheme (whe α=/, the rght had sde s a averaged betwee the curret ad future tme levels vald at +/. Relatve to ths RHS, the LHS tme dfferece becomes cetered tme. We ow that the smplest cetered dfferece scheme s secod-order accurate). Whe [ ] = 0, the scheme becomes fourth-order space. Let's perform stablty aalyss o (4) usg vo Neuma method. u = U e = U e e U λ e (5) ( x ωt) ω t j x j x j 8
9 Hereλ t e ω Substtute (5) to (4). Note that we are ow usg j as the grd pot dex. j x + j x + x x Ue ( λ λ ) = µ Ue αλ + ( α) λ ( e + e ) Dvdg Ue j x λ o both sdes, ad rearragg λ = µ x αλ + α 4 s ( / )[ ( )] λ = 4( αµ ) s ( x/) 4 s ( x/) + αµ Loo at several cases: Case I: α = 0, λ µ = 4 s ( x/) µ / as before. The scheme s codtoally stable. Case II: α = / (Cra-Ncolso) x λ = µ s ( /) + µ s ( x/) for all values of µ, therefore the scheme s absolutely or ucodtoally stable. Case III: α =, the tme dfferece s bacward, relatve to the RHS terms. 9
10 λ =, aga for all values of µ, + µ 4 s ( x/) therefore the scheme s also absolutely stable. However, ths scheme s oly frst-order accurate tme, as dscussed earler (cosstet wth the tme dfferece scheme beg u-cetered). I geeral, whe 0 α < /, t s requred that µ /( - 4α ), therefore the scheme s codtoally stable. Whe / α, the scheme s ucodtoally stable (t s sometmes referred to as the forward-based scheme). I the ARPS, the mplct dffuso scheme s a opto for treatg the vertcal turbulet mxg terms. Ths treatmet s ecessary order to remove the severe stablty costrat from these terms whe vertcal mxg s strog sde the plaetary boudary layer (PBL). The latter occurs whe the PBL s covectvely ustably ad the o-local PBL mxg s voed wth the Su ad Chag (986) parameterzato. Parameter alfcoef arps.put correspods to -α here (see hadout). Fally, we ote that for mult-tme level schemes, there s usually multple solutos for the amplfcato factor λ. some of them mght represet spurous computatoal modes due to the use of extra (artfcal) tal codtos. The expresso of λ ca be too complcated so that a graphc plottg s eeded to uderstad ts depedecy o wave umber. λ has to be o greater tha for all possble waves. The shortest wave resolvable o a grd has a wavelegth of, ad the logest s L, where L s the doma wdth. Trdagoal Solver 0
11 -D mplct method ofte leads to trdagoal systems of lear algebrac equatos. (I the ARPS, ths appears twce oce whe soud waves are treated mplctly the vertcal drecto ad oce whe the vertcal turbulece mxg s treated mplctly). For example, Eq. (4) ca be rewrtte as u u u u + u+ = α K + t x f (6) where f cotas other ow tmes at tme level +. It ca be rearraged to + tk ( ) u = α u u + u+ + tf + u x tk + tk + tk + α u + ( + α u ) α u + ) = tf + u x x x. Let A = C = αµ, B = ( + αµ ), ad d = tf + u, where tk + µ = u, we have x Au + Bu + Cu = d (7) for =,,, N-, assumg the boudares are at =0 ad N. If we have Drchlet boudary codtos,.e., u at =0 ad N are ow, the for =, the equato becomes B u + Cu = d Au (8)
12 ad for = N, the equato s A u + B u = d C u. (9) N N N N N N N For =, 3,., N-, the equato remas of the form Eq.(7). If we wrte the equatos (7-9) a matrx form, we have B C u D A B C u D..... A B C u = D..... AN BN CN un DN A B u D N N N N (0) where D = t f u Au +, 0 0 N N N N N D = t f u C u +, D = tf u for =,, N-.
13 If we have Neuma boudary codtos,.e., we ow the gradet of u at the boudares whch dscretzed form are u u0 = L ad un un = R. Plug these relatos to Eq.(7) for = ad =N-, we obta equatos smlar to (8) ad (9): ( A + B) u + Cu = d + AL () N N ( N N ) N N N A u + B + C u = d C R. () I ths case, the fal coeffcets (0) are dfferet for the frst ad last equato. Sce each except for the frst ad last row of the coeffcet matrx, oly three elemets are o-zero ad the o-zero elemets of the matrx are alged alog the dagoal axs, ths system s called trdagoal system of equatos. It ca be solved effcetly usg Thomas Algorthm. The procedure cossts of two parts. Frst, Eq.(7) s mapulated to the followg form: C ' u D' C ' u D'..... C ' u = D'..... C ' N un D' N u D' N N (3) whch the subdagoal coeffcets A are elmated ad the dagoal coeffcets are ormalzed. For the frst equato 3
14 C D C' =, D' =. (4a) B B For the geeral equatos: C D AD' C' =, D' = B AC B AC ' '. (4b) Equatos (4) represet a forward sweep step (see fgure below). It s followed by a bacward substtuto step that fds soluto u from (3). The soluto s: u = D' N N u = D' u C' for from N to. + (5) 4
15 Note both (4) ad (5) volve reducto, the algorthm s heretly o-parallelzable. Fortuately, for multdmesoal problems, multple systems of equatos ofte eed to be solved, ad oe ca explot parallelsm alog other dmesos (e.g., j stead of drecto). Read Secto of Taehll et al ad Appedx A Stablty Aalyss for Systems of Equatos Whe we are dealg wth a system of equatos, we ca also apply the vo Neuma method to fd the stablty property of a gve F.D. scheme. As wth sgle equatos, vo Neuma ca oly be used for lear systems of equatos. For olear systems, learzato has to be performed frst. Wthout gog to detals, we pot out that a system of lear equatos ca be expressed a matrx form le u u + [ A] = 0 t x The equato s frst dscretzed usg certa F.D. scheme, u ca be wrtte terms of a dscrete Fourer seres ad the wave compoet s the substtuted to the dscrete equato to obta somethg le: U + = [ M( t, x)] U (7) where U s the ampltude vector for wave at tme level, ad [M] s called the amplfcato matrx. The scheme s stable whe the maxmum absolute egevalue of [M] s o greater tha. (6) 5
16 Why the maxmum absolute egevalue? Because as you saw earler ( Chapter ) that a system of equato le (7) ca be trasformed to a system of decoupled equatos, ad the egevalues of [M] become the amplfcato factors for each of the ew depedet varables, v (the elemet of vectorv ),.e., we ca obta from (7) V [ N] V, + = where [N] s a dagoal matrx wth the egevalues of [M] as ts dagoal elemets. {[T] - [M] [T] = [N] [ T] U + = [ T] [ M][ T][ T] U V + = [ N] V where + V = [ N] U. Therefore ( v ) + = λ ( v ) where v s a elemet of vectorv.} Sce the system s stable oly whe all depedet varables rema bouded, the absolute value of the maxmum egevalue has to be o greater tha. Read secto 3.6. of Taehl et al. 6
17 Appedx Shorthad otatos for dscrete/fte dfferece operators ad dscretzato dettes Notatos: A A x j+ x+ Aj x x = δ x = + x A A δ x = A δ x = x Idettes: A j+ / j / j+ j x A x A x j j δ A δ A δ A x x x x x Aδ B δ ( A B) Bδ A x x x x δx δx( /) A A A x x
2.3. Quantitative Properties of Finite Difference Schemes. Reading: Tannehill et al. Sections and
.3. Quattatve Propertes of Fte Dfferece Schemes.3.1. Cosstecy, Covergece ad Stablty of F.D. schemes Readg: Taehll et al. Sectos 3.3.3 ad 3.3.4. Three mportat propertes of F.D. schemes: Cosstecy A F.D.
More information2.3. Quantitative Properties of Finite Difference Schemes. Reading: Tannehill et al. Sections and
.3. Quattatve Propertes of Fte Dfferece Schemes.3.. Cosstecy, Covergece ad Stablty of F.D. schemes Readg: Taehll et al. Sectos 3.3.3 ad 3.3.4. Three mportat propertes of F.D. schemes: Cosstecy A F.D. represetato
More informationLecture 07: Poles and Zeros
Lecture 07: Poles ad Zeros Defto of poles ad zeros The trasfer fucto provdes a bass for determg mportat system respose characterstcs wthout solvg the complete dfferetal equato. As defed, the trasfer fucto
More informationBeam Warming Second-Order Upwind Method
Beam Warmg Secod-Order Upwd Method Petr Valeta Jauary 6, 015 Ths documet s a part of the assessmet work for the subject 1DRP Dfferetal Equatos o Computer lectured o FNSPE CTU Prague. Abstract Ths documet
More informationChapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II
CEE49b Chapter - Free Vbrato of Mult-Degree-of-Freedom Systems - II We ca obta a approxmate soluto to the fudametal atural frequecy through a approxmate formula developed usg eergy prcples by Lord Raylegh
More informationbest estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best
Error Aalyss Preamble Wheever a measuremet s made, the result followg from that measuremet s always subject to ucertaty The ucertaty ca be reduced by makg several measuremets of the same quatty or by mprovg
More informationX X X E[ ] E X E X. is the ()m n where the ( i,)th. j element is the mean of the ( i,)th., then
Secto 5 Vectors of Radom Varables Whe workg wth several radom varables,,..., to arrage them vector form x, t s ofte coveet We ca the make use of matrx algebra to help us orgaze ad mapulate large umbers
More informationFunctions of Random Variables
Fuctos of Radom Varables Chapter Fve Fuctos of Radom Varables 5. Itroducto A geeral egeerg aalyss model s show Fg. 5.. The model output (respose) cotas the performaces of a system or product, such as weght,
More informationd dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin
Learzato of the Swg Equato We wll cover sectos.5.-.6 ad begg of Secto 3.3 these otes. 1. Sgle mache-fte bus case Cosder a sgle mache coected to a fte bus, as show Fg. 1 below. E y1 V=1./_ Fg. 1 The admttace
More informationENGI 4421 Propagation of Error Page 8-01
ENGI 441 Propagato of Error Page 8-01 Propagato of Error [Navd Chapter 3; ot Devore] Ay realstc measuremet procedure cotas error. Ay calculatos based o that measuremet wll therefore also cota a error.
More informationEnd of Finite Volume Methods Cartesian grids. Solution of the Navier-Stokes Equations. REVIEW Lecture 17: Higher order (interpolation) schemes
REVIEW Lecture 17: Numercal Flud Mechacs Sprg 2015 Lecture 18 Ed of Fte Volume Methods Cartesa grds Hgher order (terpolato) schemes Soluto of the Naver-Stokes Equatos Dscretzato of the covectve ad vscous
More informationCHAPTER 4 RADICAL EXPRESSIONS
6 CHAPTER RADICAL EXPRESSIONS. The th Root of a Real Number A real umber a s called the th root of a real umber b f Thus, for example: s a square root of sce. s also a square root of sce ( ). s a cube
More informationCubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem
Joural of Amerca Scece ;6( Cubc Nopolyomal Sple Approach to the Soluto of a Secod Order Two-Pot Boudary Value Problem W.K. Zahra, F.A. Abd El-Salam, A.A. El-Sabbagh ad Z.A. ZAk * Departmet of Egeerg athematcs
More informationNumerical Analysis Formulae Booklet
Numercal Aalyss Formulae Booklet. Iteratve Scemes for Systems of Lear Algebrac Equatos:.... Taylor Seres... 3. Fte Dfferece Approxmatos... 3 4. Egevalues ad Egevectors of Matrces.... 3 5. Vector ad Matrx
More informationMOLECULAR VIBRATIONS
MOLECULAR VIBRATIONS Here we wsh to vestgate molecular vbratos ad draw a smlarty betwee the theory of molecular vbratos ad Hückel theory. 1. Smple Harmoc Oscllator Recall that the eergy of a oe-dmesoal
More informationUNIT 2 SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS
Numercal Computg -I UNIT SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS Structure Page Nos..0 Itroducto 6. Objectves 7. Ital Approxmato to a Root 7. Bsecto Method 8.. Error Aalyss 9.4 Regula Fals Method
More informationThird handout: On the Gini Index
Thrd hadout: O the dex Corrado, a tala statstca, proposed (, 9, 96) to measure absolute equalt va the mea dfferece whch s defed as ( / ) where refers to the total umber of dvduals socet. Assume that. The
More information1 Lyapunov Stability Theory
Lyapuov Stablty heory I ths secto we cosder proofs of stablty of equlbra of autoomous systems. hs s stadard theory for olear systems, ad oe of the most mportat tools the aalyss of olear systems. It may
More informationEECE 301 Signals & Systems
EECE 01 Sgals & Systems Prof. Mark Fowler Note Set #9 Computg D-T Covoluto Readg Assgmet: Secto. of Kame ad Heck 1/ Course Flow Dagram The arrows here show coceptual flow betwee deas. Note the parallel
More informationPGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation
PGE 30: Formulato ad Soluto Geosystems Egeerg Dr. Balhoff Iterpolato Numercal Methods wth MATLAB, Recktewald, Chapter 0 ad Numercal Methods for Egeers, Chapra ad Caale, 5 th Ed., Part Fve, Chapter 8 ad
More informationA Primer on Summation Notation George H Olson, Ph. D. Doctoral Program in Educational Leadership Appalachian State University Spring 2010
Summato Operator A Prmer o Summato otato George H Olso Ph D Doctoral Program Educatoal Leadershp Appalacha State Uversty Sprg 00 The summato operator ( ) {Greek letter captal sgma} s a structo to sum over
More information2006 Jamie Trahan, Autar Kaw, Kevin Martin University of South Florida United States of America
SOLUTION OF SYSTEMS OF SIMULTANEOUS LINEAR EQUATIONS Gauss-Sedel Method 006 Jame Traha, Autar Kaw, Kev Mart Uversty of South Florda Uted States of Amerca kaw@eg.usf.edu Itroducto Ths worksheet demostrates
More informationUNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS
UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS Exam: ECON430 Statstcs Date of exam: Frday, December 8, 07 Grades are gve: Jauary 4, 08 Tme for exam: 0900 am 00 oo The problem set covers 5 pages Resources allowed:
More information5 Short Proofs of Simplified Stirling s Approximation
5 Short Proofs of Smplfed Strlg s Approxmato Ofr Gorodetsky, drtymaths.wordpress.com Jue, 20 0 Itroducto Strlg s approxmato s the followg (somewhat surprsg) approxmato of the factoral,, usg elemetary fuctos:
More informationMultiple Choice Test. Chapter Adequacy of Models for Regression
Multple Choce Test Chapter 06.0 Adequac of Models for Regresso. For a lear regresso model to be cosdered adequate, the percetage of scaled resduals that eed to be the rage [-,] s greater tha or equal to
More informationApplying the condition for equilibrium to this equilibrium, we get (1) n i i =, r G and 5 i
CHEMICAL EQUILIBRIA The Thermodyamc Equlbrum Costat Cosder a reversble reacto of the type 1 A 1 + 2 A 2 + W m A m + m+1 A m+1 + Assgg postve values to the stochometrc coeffcets o the rght had sde ad egatve
More informationLaboratory I.10 It All Adds Up
Laboratory I. It All Adds Up Goals The studet wll work wth Rema sums ad evaluate them usg Derve. The studet wll see applcatos of tegrals as accumulatos of chages. The studet wll revew curve fttg sklls.
More informationmeans the first term, a2 means the term, etc. Infinite Sequences: follow the same pattern forever.
9.4 Sequeces ad Seres Pre Calculus 9.4 SEQUENCES AND SERIES Learg Targets:. Wrte the terms of a explctly defed sequece.. Wrte the terms of a recursvely defed sequece. 3. Determe whether a sequece s arthmetc,
More informationLecture 12 APPROXIMATION OF FIRST ORDER DERIVATIVES
FDM: Appromato of Frst Order Dervatves Lecture APPROXIMATION OF FIRST ORDER DERIVATIVES. INTRODUCTION Covectve term coservato equatos volve frst order dervatves. The smplest possble approach for dscretzato
More informationCan we take the Mysticism Out of the Pearson Coefficient of Linear Correlation?
Ca we tae the Mstcsm Out of the Pearso Coeffcet of Lear Correlato? Itroducto As the ttle of ths tutoral dcates, our purpose s to egeder a clear uderstadg of the Pearso coeffcet of lear correlato studets
More informationChapter 9 Jordan Block Matrices
Chapter 9 Jorda Block atrces I ths chapter we wll solve the followg problem. Gve a lear operator T fd a bass R of F such that the matrx R (T) s as smple as possble. f course smple s a matter of taste.
More informationCHAPTER VI Statistical Analysis of Experimental Data
Chapter VI Statstcal Aalyss of Expermetal Data CHAPTER VI Statstcal Aalyss of Expermetal Data Measuremets do ot lead to a uque value. Ths s a result of the multtude of errors (maly radom errors) that ca
More informationIntroduction to local (nonparametric) density estimation. methods
Itroducto to local (oparametrc) desty estmato methods A slecture by Yu Lu for ECE 66 Sprg 014 1. Itroducto Ths slecture troduces two local desty estmato methods whch are Parze desty estmato ad k-earest
More informationNewton s Power Flow algorithm
Power Egeerg - Egll Beedt Hresso ewto s Power Flow algorthm Power Egeerg - Egll Beedt Hresso The ewto s Method of Power Flow 2 Calculatos. For the referece bus #, we set : V = p.u. ad δ = 0 For all other
More informationTESTS BASED ON MAXIMUM LIKELIHOOD
ESE 5 Toy E. Smth. The Basc Example. TESTS BASED ON MAXIMUM LIKELIHOOD To llustrate the propertes of maxmum lkelhood estmates ad tests, we cosder the smplest possble case of estmatg the mea of the ormal
More informationThe number of observed cases The number of parameters. ith case of the dichotomous dependent variable. the ith case of the jth parameter
LOGISTIC REGRESSION Notato Model Logstc regresso regresses a dchotomous depedet varable o a set of depedet varables. Several methods are mplemeted for selectg the depedet varables. The followg otato s
More informationu(x, t) = u 0 (x ct). This Riemann invariant u is constant along characteristics λ with x = x 0 +ct (u(x, t) = u 0 (x 0 )):
x, t, h x The Frst-Order Wave Eqato The frst-order wave advecto eqato s c > 0 t + c x = 0, x, t = 0 = 0x. The solto propagates the tal data 0 to the rght wth speed c: x, t = 0 x ct. Ths Rema varat s costat
More informationLecture Note to Rice Chapter 8
ECON 430 HG revsed Nov 06 Lecture Note to Rce Chapter 8 Radom matrces Let Y, =,,, m, =,,, be radom varables (r.v. s). The matrx Y Y Y Y Y Y Y Y Y Y = m m m s called a radom matrx ( wth a ot m-dmesoal dstrbuto,
More informationRademacher Complexity. Examples
Algorthmc Foudatos of Learg Lecture 3 Rademacher Complexty. Examples Lecturer: Patrck Rebesch Verso: October 16th 018 3.1 Itroducto I the last lecture we troduced the oto of Rademacher complexty ad showed
More informationSTATISTICAL PROPERTIES OF LEAST SQUARES ESTIMATORS. x, where. = y - ˆ " 1
STATISTICAL PROPERTIES OF LEAST SQUARES ESTIMATORS Recall Assumpto E(Y x) η 0 + η x (lear codtoal mea fucto) Data (x, y ), (x 2, y 2 ),, (x, y ) Least squares estmator ˆ E (Y x) ˆ " 0 + ˆ " x, where ˆ
More informationBlock-Based Compact Thermal Modeling of Semiconductor Integrated Circuits
Block-Based Compact hermal Modelg of Semcoductor Itegrated Crcuts Master s hess Defese Caddate: Jg Ba Commttee Members: Dr. Mg-Cheg Cheg Dr. Daqg Hou Dr. Robert Schllg July 27, 2009 Outle Itroducto Backgroud
More informationInitial-Value Problems for ODEs. numerical errors (round-off and truncation errors) Consider a perturbed system: dz dt
Ital-Value Problems or ODEs d GIVEN: t t,, a FIND: t or atb umercal errors (roud-o ad trucato errors) Cosder a perturbed sstem: dz t, z t, at b z a a Does z(t) (t)? () (uqueess) a uque soluto (t) exsts
More informationCS286.2 Lecture 4: Dinur s Proof of the PCP Theorem
CS86. Lecture 4: Dur s Proof of the PCP Theorem Scrbe: Thom Bohdaowcz Prevously, we have prove a weak verso of the PCP theorem: NP PCP 1,1/ (r = poly, q = O(1)). Wth ths result we have the desred costat
More informationESS Line Fitting
ESS 5 014 17. Le Fttg A very commo problem data aalyss s lookg for relatoshpetwee dfferet parameters ad fttg les or surfaces to data. The smplest example s fttg a straght le ad we wll dscuss that here
More informationCHAPTER 3 POSTERIOR DISTRIBUTIONS
CHAPTER 3 POSTERIOR DISTRIBUTIONS If scece caot measure the degree of probablt volved, so much the worse for scece. The practcal ma wll stck to hs apprecatve methods utl t does, or wll accept the results
More informationC-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory
ROAD MAP... AE301 Aerodyamcs I UNIT C: 2-D Arfols C-1: Aerodyamcs of Arfols 1 C-2: Aerodyamcs of Arfols 2 C-3: Pael Methods C-4: Th Arfol Theory AE301 Aerodyamcs I Ut C-3: Lst of Subects Problem Solutos?
More informationF. Inequalities. HKAL Pure Mathematics. 進佳數學團隊 Dr. Herbert Lam 林康榮博士. [Solution] Example Basic properties
進佳數學團隊 Dr. Herbert Lam 林康榮博士 HKAL Pure Mathematcs F. Ieualtes. Basc propertes Theorem Let a, b, c be real umbers. () If a b ad b c, the a c. () If a b ad c 0, the ac bc, but f a b ad c 0, the ac bc. Theorem
More informationAssignment 7/MATH 247/Winter, 2010 Due: Friday, March 19. Powers of a square matrix
Assgmet 7/MATH 47/Wter, 00 Due: Frday, March 9 Powers o a square matrx Gve a square matrx A, ts powers A or large, or eve arbtrary, teger expoets ca be calculated by dagoalzg A -- that s possble (!) Namely,
More informationMATH 247/Winter Notes on the adjoint and on normal operators.
MATH 47/Wter 00 Notes o the adjot ad o ormal operators I these otes, V s a fte dmesoal er product space over, wth gve er * product uv, T, S, T, are lear operators o V U, W are subspaces of V Whe we say
More information1 Convergence of the Arnoldi method for eigenvalue problems
Lecture otes umercal lear algebra Arold method covergece Covergece of the Arold method for egevalue problems Recall that, uless t breaks dow, k steps of the Arold method geerates a orthogoal bass of a
More informationECE 595, Section 10 Numerical Simulations Lecture 19: FEM for Electronic Transport. Prof. Peter Bermel February 22, 2013
ECE 595, Secto 0 Numercal Smulatos Lecture 9: FEM for Electroc Trasport Prof. Peter Bermel February, 03 Outle Recap from Wedesday Physcs-based devce modelg Electroc trasport theory FEM electroc trasport
More information( ) = ( ) ( ) Chapter 13 Asymptotic Theory and Stochastic Regressors. Stochastic regressors model
Chapter 3 Asmptotc Theor ad Stochastc Regressors The ature of eplaator varable s assumed to be o-stochastc or fed repeated samples a regresso aalss Such a assumpto s approprate for those epermets whch
More informationWe have already referred to a certain reaction, which takes place at high temperature after rich combustion.
ME 41 Day 13 Topcs Chemcal Equlbrum - Theory Chemcal Equlbrum Example #1 Equlbrum Costats Chemcal Equlbrum Example #2 Chemcal Equlbrum of Hot Bured Gas 1. Chemcal Equlbrum We have already referred to a
More informationAssignment 5/MATH 247/Winter Due: Friday, February 19 in class (!) (answers will be posted right after class)
Assgmet 5/MATH 7/Wter 00 Due: Frday, February 9 class (!) (aswers wll be posted rght after class) As usual, there are peces of text, before the questos [], [], themselves. Recall: For the quadratc form
More informationLogistic regression (continued)
STAT562 page 138 Logstc regresso (cotued) Suppose we ow cosder more complex models to descrbe the relatoshp betwee a categorcal respose varable (Y) that takes o two (2) possble outcomes ad a set of p explaatory
More informationx y exp λ'. x exp λ 2. x exp 1.
egecosmcd Egevalue-egevector of the secod dervatve operator d /d hs leads to Fourer seres (se, cose, Legedre, Bessel, Chebyshev, etc hs s a eample of a systematc way of geeratg a set of mutually orthogoal
More informationInvestigation of Partially Conditional RP Model with Response Error. Ed Stanek
Partally Codtoal Radom Permutato Model 7- vestgato of Partally Codtoal RP Model wth Respose Error TRODUCTO Ed Staek We explore the predctor that wll result a smple radom sample wth respose error whe a
More informationANALYSIS ON THE NATURE OF THE BASIC EQUATIONS IN SYNERGETIC INTER-REPRESENTATION NETWORK
Far East Joural of Appled Mathematcs Volume, Number, 2008, Pages Ths paper s avalable ole at http://www.pphm.com 2008 Pushpa Publshg House ANALYSIS ON THE NATURE OF THE ASI EQUATIONS IN SYNERGETI INTER-REPRESENTATION
More informationFeature Selection: Part 2. 1 Greedy Algorithms (continued from the last lecture)
CSE 546: Mache Learg Lecture 6 Feature Selecto: Part 2 Istructor: Sham Kakade Greedy Algorthms (cotued from the last lecture) There are varety of greedy algorthms ad umerous amg covetos for these algorthms.
More informationLecture 3 Probability review (cont d)
STATS 00: Itroducto to Statstcal Iferece Autum 06 Lecture 3 Probablty revew (cot d) 3. Jot dstrbutos If radom varables X,..., X k are depedet, the ther dstrbuto may be specfed by specfyg the dvdual dstrbuto
More informationPart 4b Asymptotic Results for MRR2 using PRESS. Recall that the PRESS statistic is a special type of cross validation procedure (see Allen (1971))
art 4b Asymptotc Results for MRR usg RESS Recall that the RESS statstc s a specal type of cross valdato procedure (see Alle (97)) partcular to the regresso problem ad volves fdg Y $,, the estmate at the
More information{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:
Chapter 4 Exercses Samplg Theory Exercse (Smple radom samplg: Let there be two correlated radom varables X ad A sample of sze s draw from a populato by smple radom samplg wthout replacemet The observed
More informationEconometric Methods. Review of Estimation
Ecoometrc Methods Revew of Estmato Estmatg the populato mea Radom samplg Pot ad terval estmators Lear estmators Ubased estmators Lear Ubased Estmators (LUEs) Effcecy (mmum varace) ad Best Lear Ubased Estmators
More informationDerivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations
Dervato of -Pot Block Method Formula for Solvg Frst Order Stff Ordary Dfferetal Equatos Kharul Hamd Kharul Auar, Kharl Iskadar Othma, Zara Bb Ibrahm Abstract Dervato of pot block method formula wth costat
More informationBounds on the expected entropy and KL-divergence of sampled multinomial distributions. Brandon C. Roy
Bouds o the expected etropy ad KL-dvergece of sampled multomal dstrbutos Brado C. Roy bcroy@meda.mt.edu Orgal: May 18, 2011 Revsed: Jue 6, 2011 Abstract Iformato theoretc quattes calculated from a sampled
More informationPTAS for Bin-Packing
CS 663: Patter Matchg Algorthms Scrbe: Che Jag /9/00. Itroducto PTAS for B-Packg The B-Packg problem s NP-hard. If we use approxmato algorthms, the B-Packg problem could be solved polyomal tme. For example,
More informationDiscrete Mathematics and Probability Theory Fall 2016 Seshia and Walrand DIS 10b
CS 70 Dscrete Mathematcs ad Probablty Theory Fall 206 Sesha ad Walrad DIS 0b. Wll I Get My Package? Seaky delvery guy of some compay s out delverg packages to customers. Not oly does he had a radom package
More informationL5 Polynomial / Spline Curves
L5 Polyomal / Sple Curves Cotets Coc sectos Polyomal Curves Hermte Curves Bezer Curves B-Sples No-Uform Ratoal B-Sples (NURBS) Mapulato ad Represetato of Curves Types of Curve Equatos Implct: Descrbe a
More informationAnalysis of System Performance IN2072 Chapter 5 Analysis of Non Markov Systems
Char for Network Archtectures ad Servces Prof. Carle Departmet of Computer Scece U Müche Aalyss of System Performace IN2072 Chapter 5 Aalyss of No Markov Systems Dr. Alexader Kle Prof. Dr.-Ig. Georg Carle
More informationFourth Order Four-Stage Diagonally Implicit Runge-Kutta Method for Linear Ordinary Differential Equations ABSTRACT INTRODUCTION
Malasa Joural of Mathematcal Sceces (): 95-05 (00) Fourth Order Four-Stage Dagoall Implct Ruge-Kutta Method for Lear Ordar Dfferetal Equatos Nur Izzat Che Jawas, Fudzah Ismal, Mohamed Sulema, 3 Azm Jaafar
More informationLyapunov Stability. Aleksandr Mikhailovich Lyapunov [1] 1 Autonomous Systems. R into Nonlinear Systems in Mechanical Engineering Lesson 5
Joh vo Neuma Perre de Fermat Joseph Fourer etc 858 Nolear Systems Mechacal Egeerg Lesso 5 Aleksadr Mkhalovch Lyapuov [] Bor: 6 Jue 857 Yaroslavl Russa Ded: 3 November 98 Bega hs educato at home Graduated
More informationAhmed Elgamal. MDOF Systems & Modal Analysis
DOF Systems & odal Aalyss odal Aalyss (hese otes cover sectos from Ch. 0, Dyamcs of Structures, Al Chopra, Pretce Hall, 995). Refereces Dyamcs of Structures, Al K. Chopra, Pretce Hall, New Jersey, ISBN
More informationUNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS
UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS Postpoed exam: ECON430 Statstcs Date of exam: Jauary 0, 0 Tme for exam: 09:00 a.m. :00 oo The problem set covers 5 pages Resources allowed: All wrtte ad prted
More informationHomework 1: Solutions Sid Banerjee Problem 1: (Practice with Asymptotic Notation) ORIE 4520: Stochastics at Scale Fall 2015
Fall 05 Homework : Solutos Problem : (Practce wth Asymptotc Notato) A essetal requremet for uderstadg scalg behavor s comfort wth asymptotc (or bg-o ) otato. I ths problem, you wll prove some basc facts
More informationBayes (Naïve or not) Classifiers: Generative Approach
Logstc regresso Bayes (Naïve or ot) Classfers: Geeratve Approach What do we mea by Geeratve approach: Lear p(y), p(x y) ad the apply bayes rule to compute p(y x) for makg predctos Ths s essetally makg
More informationChapter 5 Properties of a Random Sample
Lecture 6 o BST 63: Statstcal Theory I Ku Zhag, /0/008 Revew for the prevous lecture Cocepts: t-dstrbuto, F-dstrbuto Theorems: Dstrbutos of sample mea ad sample varace, relatoshp betwee sample mea ad sample
More informationFractional Order Finite Difference Scheme For Soil Moisture Diffusion Equation And Its Applications
IOS Joural of Mathematcs (IOS-JM e-iss: 78-578. Volume 5, Issue 4 (Ja. - Feb. 3, PP -8 www.osrourals.org Fractoal Order Fte Dfferece Scheme For Sol Mosture Dffuso quato Ad Its Applcatos S.M.Jogdad, K.C.Takale,
More informationQR Factorization and Singular Value Decomposition COS 323
QR Factorzato ad Sgular Value Decomposto COS 33 Why Yet Aother Method? How do we solve least-squares wthout currg codto-squarg effect of ormal equatos (A T A A T b) whe A s sgular, fat, or otherwse poorly-specfed?
More informationLecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model
Lecture 7. Cofdece Itervals ad Hypothess Tests the Smple CLR Model I lecture 6 we troduced the Classcal Lear Regresso (CLR) model that s the radom expermet of whch the data Y,,, K, are the outcomes. The
More informationGeneral Method for Calculating Chemical Equilibrium Composition
AE 6766/Setzma Sprg 004 Geeral Metod for Calculatg Cemcal Equlbrum Composto For gve tal codtos (e.g., for gve reactats, coose te speces to be cluded te products. As a example, for combusto of ydroge wt
More informationChapter 4 (Part 1): Non-Parametric Classification (Sections ) Pattern Classification 4.3) Announcements
Aoucemets No-Parametrc Desty Estmato Techques HW assged Most of ths lecture was o the blacboard. These sldes cover the same materal as preseted DHS Bometrcs CSE 90-a Lecture 7 CSE90a Fall 06 CSE90a Fall
More informationLecture 16: Backpropogation Algorithm Neural Networks with smooth activation functions
CO-511: Learg Theory prg 2017 Lecturer: Ro Lv Lecture 16: Bacpropogato Algorthm Dsclamer: These otes have ot bee subected to the usual scruty reserved for formal publcatos. They may be dstrbuted outsde
More informationPoint Estimation: definition of estimators
Pot Estmato: defto of estmators Pot estmator: ay fucto W (X,..., X ) of a data sample. The exercse of pot estmato s to use partcular fuctos of the data order to estmate certa ukow populato parameters.
More informationEstimation of Stress- Strength Reliability model using finite mixture of exponential distributions
Iteratoal Joural of Computatoal Egeerg Research Vol, 0 Issue, Estmato of Stress- Stregth Relablty model usg fte mxture of expoetal dstrbutos K.Sadhya, T.S.Umamaheswar Departmet of Mathematcs, Lal Bhadur
More informationDKA method for single variable holomorphic functions
DKA method for sgle varable holomorphc fuctos TOSHIAKI ITOH Itegrated Arts ad Natural Sceces The Uversty of Toushma -, Mamhosama, Toushma, 770-8502 JAPAN Abstract: - Durad-Kerer-Aberth (DKA method for
More informationECE606: Solid State Devices Lecture 13 Solutions of the Continuity Eqs. Analytical & Numerical
ECE66: Sold State Devces Lecture 13 Solutos of the Cotuty Eqs. Aalytcal & Numercal Gerhard Klmeck gekco@purdue.edu Outle Aalytcal Solutos to the Cotuty Equatos 1) Example problems ) Summary Numercal Solutos
More informationFor combinatorial problems we might need to generate all permutations, combinations, or subsets of a set.
Addtoal Decrease ad Coquer Algorthms For combatoral problems we mght eed to geerate all permutatos, combatos, or subsets of a set. Geeratg Permutatos If we have a set f elemets: { a 1, a 2, a 3, a } the
More information( q Modal Analysis. Eigenvectors = Mode Shapes? Eigenproblem (cont) = x x 2 u 2. u 1. x 1 (4.55) vector and M and K are matrices.
4.3 - Modal Aalyss Physcal coordates are ot always the easest to work Egevectors provde a coveet trasformato to modal coordates Modal coordates are lear combato of physcal coordates Say we have physcal
More informationOverview of the weighting constants and the points where we evaluate the function for The Gaussian quadrature Project two
Overvew of the weghtg costats ad the pots where we evaluate the fucto for The Gaussa quadrature Project two By Ashraf Marzouk ChE 505 Fall 005 Departmet of Mechacal Egeerg Uversty of Teessee Koxvlle, TN
More informationDATE: 21 September, 1999 TO: Jim Russell FROM: Peter Tkacik RE: Analysis of wide ply tube winding as compared to Konva Kore CC: Larry McMillan
M E M O R A N D U M DATE: 1 September, 1999 TO: Jm Russell FROM: Peter Tkack RE: Aalyss of wde ply tube wdg as compared to Kova Kore CC: Larry McMlla The goal of ths report s to aalyze the spral tube wdg
More information8.1 Hashing Algorithms
CS787: Advaced Algorthms Scrbe: Mayak Maheshwar, Chrs Hrchs Lecturer: Shuch Chawla Topc: Hashg ad NP-Completeess Date: September 21 2007 Prevously we looked at applcatos of radomzed algorthms, ad bega
More informationNumerical Simulations of the Complex Modied Korteweg-de Vries Equation. Thiab R. Taha. The University of Georgia. Abstract
Numercal Smulatos of the Complex Moded Korteweg-de Vres Equato Thab R. Taha Computer Scece Departmet The Uversty of Georga Athes, GA 002 USA Tel 0-542-2911 e-mal thab@cs.uga.edu Abstract I ths paper mplemetatos
More informationGeneralized Linear Regression with Regularization
Geeralze Lear Regresso wth Regularzato Zoya Bylsk March 3, 05 BASIC REGRESSION PROBLEM Note: I the followg otes I wll make explct what s a vector a what s a scalar usg vec t or otato, to avo cofuso betwee
More informationMu Sequences/Series Solutions National Convention 2014
Mu Sequeces/Seres Solutos Natoal Coveto 04 C 6 E A 6C A 6 B B 7 A D 7 D C 7 A B 8 A B 8 A C 8 E 4 B 9 B 4 E 9 B 4 C 9 E C 0 A A 0 D B 0 C C Usg basc propertes of arthmetc sequeces, we fd a ad bm m We eed
More informationf f... f 1 n n (ii) Median : It is the value of the middle-most observation(s).
CHAPTER STATISTICS Pots to Remember :. Facts or fgures, collected wth a defte pupose, are called Data.. Statstcs s the area of study dealg wth the collecto, presetato, aalyss ad terpretato of data.. The
More informationInvestigating Cellular Automata
Researcher: Taylor Dupuy Advsor: Aaro Wootto Semester: Fall 4 Ivestgatg Cellular Automata A Overvew of Cellular Automata: Cellular Automata are smple computer programs that geerate rows of black ad whte
More informationECONOMETRIC THEORY. MODULE VIII Lecture - 26 Heteroskedasticity
ECONOMETRIC THEORY MODULE VIII Lecture - 6 Heteroskedastcty Dr. Shalabh Departmet of Mathematcs ad Statstcs Ida Isttute of Techology Kapur . Breusch Paga test Ths test ca be appled whe the replcated data
More informationTransforms that are commonly used are separable
Trasforms s Trasforms that are commoly used are separable Eamples: Two-dmesoal DFT DCT DST adamard We ca the use -D trasforms computg the D separable trasforms: Take -D trasform of the rows > rows ( )
More informationThe Selection Problem - Variable Size Decrease/Conquer (Practice with algorithm analysis)
We have covered: Selecto, Iserto, Mergesort, Bubblesort, Heapsort Next: Selecto the Qucksort The Selecto Problem - Varable Sze Decrease/Coquer (Practce wth algorthm aalyss) Cosder the problem of fdg the
More information1 Mixed Quantum State. 2 Density Matrix. CS Density Matrices, von Neumann Entropy 3/7/07 Spring 2007 Lecture 13. ψ = α x x. ρ = p i ψ i ψ i.
CS 94- Desty Matrces, vo Neuma Etropy 3/7/07 Sprg 007 Lecture 3 I ths lecture, we wll dscuss the bascs of quatum formato theory I partcular, we wll dscuss mxed quatum states, desty matrces, vo Neuma etropy
More information