Chapter 7. Solution of Ordinary Differential Equations
|
|
- Catherine Phillips
- 6 years ago
- Views:
Transcription
1 Ordar Dfferetal Equatos - 4 Chapter 7. Soluto of Ordar Dfferetal Equatos 7.. Itroducto The damc behavor of ma relevat sstems ad materals ca be descrbed wth ordar dfferetal equatos ODEs. I ths chapter, we provde a troducto to the techques for umercal soluto of ODEs. We beg wth a sgle, frst-order ODE tal value problem. We the et the process to hgh-order ODEs, sstems of ODEs ad boudar value problems. The techques descrbed ths chapter work for both lear ad olear ODEs. However, we pot out that for lear ODEs ad some o-lear ODEs there ma est more elegat aaltcal solutos. I ths book, we do ot vestgate the aaltcal solutos but lmt ourselves to umercal methods that provde vald solutos. 7.. Ital Value Problems Regardless of whether oe ts to solve a ordar dfferetal equato ODE wth aaltcal or umercal techques, the problem must frst be properl posed. Let us tall lmt ourselves to a sgle, frst-order ODE, volvg a sgle depet varable ad a sgle depet varable. Dfferetal equatos that volve more tha oe depet varable are called partal dfferetal equatos PDEs ad are ot cosdered ths book. Dfferetal equatos that volve more tha oe depet varable costtute sstems of ODEs ad addressed later ths chapter. The soluto to a ODE s a fucto,. The most geeral formulato of ths ODE s f,, 7. where we have voked the shorthad otato, d d, for the dervatve of wth respect to. Ofte, t s possble through algebrac mapulato of equato 7. to solate the dervatve o the LHS, f, 7.
2 Ordar Dfferetal Equatos - 5 Ths s perhaps the most famlar form of a ODE, but t s ot the most geeral. ODEs are frequetl categorzed as lear ad olear ODEs. It s mportat to remember that the leart of the ODE s defed b the leart of ol ad specfcall ot of. The ukow must be operated o eclusvel b lear operators order for the ODE to be cosdered lear. Recall that the dfferetal operator s a lear operator. The geeral form of a lear frst order ODE s a b c 7.3 The coeffcets, a, b ad c ca be olear fuctos of. Regardless of whether oe s solvg a ODE of the form gve 7., 7. or 7.3, the problem s ot properl posed ad does ot have a uque soluto utl a tal codto s suppled. The tal codto provdes a value of the fucto at oe pot. 7.4 The tal value problem IVP s the combato of the ODE ad the tal codto Euler Method It s structve to beg our stud of umercal techques for solvg ODEs wth a frst order method. I practce, we wll ever use a frst-order method to solve a ODE, due to ts low accurac. Nevertheless, the frst-order method plas a mportat role as a structoal tool. The Euler method reles o a Talor seres trucated after the lear term. f df f h O h 3. d whch ca be rewrtte the omeclature of ODEs as 7.5 where the dscretzato s gve b. Equato 7.5 s the Euler method. Oe begs at the tal codto,, as gve the problem statemet. Oe choose a dscretzato,, whch s a crucal decso. Oe the evaluates the dervatve at usg the ODE equato 7., 7. or 7.3, whch also was gve the problem statemet. The Euler method equato 7.5 ca the be used drectl to estmate the value of the ukow fucto at
3 Ordar Dfferetal Equatos - 6. Ths process ca be repeated movg further dow the -as for as log as terest remas. Before we troduce specfc eamples, two pots of dscusso are ecessar. Frst, the choce of dscretzato or step-sze s crucal. It s true that the smaller the step-sze the more accurate the soluto. However, we must balace the desre for accurac wth the eed for computatoal effcec. If our choce of step-sze s so small that we eed a tractable umber of steps, the that does t help us. If the step sze s too bg, the method wll fal. I fact, for the smple methods preseted ths book, the ol reaso a umercal route for solvg a ODE wll fal s because the step-sze was too large. The reaso ca be see Fgure 7.. I the top fgure, t s clear that as the step-sze creases, the error creases. The error s quatfed the mddle fgure. For some fuctos t s possble that the error leads to values of the fucto that are ot permtted, such as egatve values of fuctos that requre a square root. The secod pot of terest s the determato of whe to stop the Euler teratos. The fal value of s smpl chose b the user. Frequetl, the stoppg pot s determed whe oe has reached a stead state. For eample the bottom part of Fgure 7., the fucto clearl goes to zero. Cotug to solve the ODE out to large values of tpcall serves o useful purpose. If oe does t kow how log t wll take to reach a stead state, the oe ma have to guess the fal value. If t s too short, the oe ca smpl et the Euler procedure usg the fal pot of the prevous process as the tal codto for the f relatve error f aaltcal soluto Euler d =. Euler d=. Euler d =.5 tal codto Euler d =. Euler d=. Euler d = aaltcal soluto Euler d =. tal codto Fgure 7. Applcato of Euler method to =- wth tal codto ==. Top: soluto up to =. Mddle: error up to =. Bottom: soluto up to =.
4 Ordar Dfferetal Equatos - 7 cotuato. A MATLAB code whch mplemets the Euler method for a sgle frst-order ODE s provded later ths chapter Classcal Fourth-Order Ruge-Kutta Method The obvous soluto to accurac of the Euler method s to crease the order of the method. There are umerous flavors of ODE solvers of all orders. Here we smpl preset the Classcal Fourth-Order Ruge-Kutta RK4 Method, whch s a tred ad true soluto techque. The RK4 method proceeds as does the Euler method b startg at the tal codto ad followg a estmate of the average slope over the terval,, out to the et dscretzato pot. 7.6 I the Euler method, the estmate of the average slope was, smpl the value of the slope at the begg of the terval. The RK4 method uses a hgher order appromato for the average slope over the terval, amel where k k k3 k k k k k 3 4 ' ' ' ',, k, k, k Wthout a rgorous dervato, we ca see that the Ruge-Kutta method mproves the qualt of the estmate of the slope over the terval b evaluatg the slope at the begg, mddle twce ad of the terval. The RK4 method s a fourth-order method, so the error decreases as the stepsze to the fourth power. I Fgure 7., we preset a applcato of the RK4 method. I the left fgure we show that eve for coarse dscretzato, the RK4 method ca be ver accurate. I the rght fgure we quatf the relatve error. Note that the error s ow show o a log as
5 Ordar Dfferetal Equatos - 8.E+..E- aaltcal soluto RK4 d =. RK4 d=.5.e- Euler d =. tal codto relatve error f.8.6.e-3 RK4 d =..E-4 RK4 d=.5 Euler d =..E-5.4.E-6..E-7.E Fgure 7.. Applcato of RK4 method to =- wth tal codto ==. Left: soluto up to =. Rght: error up to =. A MATLAB code whch mplemets the classcal fourth-order Ruge-Kutta method for a sgle frst-order ODE s provded later ths chapter Applcato to Sstems of Ordar Dfferetal Equatos The eteso of the Euler or Ruge-Kutta method to sstems of ODEs s ver straghtforward for a tal value problem. Let s suppose that we have ODEs, each of whch ca be wrtte the form of equato 7. as d f,,, 3,..., d d f,,, 3,..., d... d f,,, 3,..., d 7.9 wth tal codtos:, for to 7. The fact that all of the codtos are gve at the same value of s what makes ths problem a IVP. Note carefull, that there are ODE s, thus fuctos, f, ad ukow fuctos,,
6 Ordar Dfferetal Equatos - 9 but there s ol oe depet varable. Ths s what makes ths sstem a sstem of ODEs rather tha PDEs. The techques ths chapter wll ot solve PDEs. The Euler method for a sstem of of ODEs ca be wrtte as 3 * f,,,,..., for = to 7. There s o dfferece betwee ths equato ad the equato for the sgle sstem ODE usg Euler s method equato 7.5. I ths case, remember that the subscrpt attached to the ad the f deotes dfferet fuctos. The subscrpt attached to the varable deotes steps or teratos. The Euler algorthm for a sstem of ODEs smpl requres oe to evaluate all f at terato from the ODEs the problem statemet equato 7.9, the compute all at terato + from equato 7. ad repeat the process. The eteso of RK4 to sstems of ODE s s ust as smple as the eteso of Euler s method. However, because RK4 has a lttle more sophstcato, the eteso looks more complcated, whe t s reall ot. The RK4 equato for a sstem of equatos s gve b 6 k k k k h,, 3, 4, for = to 7. where k k k k,, 3, 4, f, m f, m k m 3, m, m f, m k f, k, m for = to 7.3 I equato 7.3, the braces m represet the set of all m from m = to., all of whch are evaluated at. The mult-ode RK4 algorthm requres that oe frst evaluate all k,. Wth k, had, oe ca compute the argumets requred for k,. Oce all of the k, are kow, oe ca et compute the argumets requred for k 3,. Oce all of the k 3, are kow, oe ca et compute the argumets requred for k 4,. Oce all the k,, k,, k 3, ad k 4, are kow, oe ca use the RK4. method equato 7. to evaluate the value of the ukow fuctos at the et step, MATLAB codes whch mplemet the Euler method ad the classcal fourth-order Ruge- Kutta method for a sstem of frst-order ODEs are provded later ths chapter.
7 Ordar Dfferetal Equatos Hgher-Order ODEs Regardless of whether oe s pursug a aaltcal or umercal soluto to a hgher tha frst order ODE, there s a smple trck to covertg a th -order ODE to a sstem of frst-order ODEs. You make a substtuto that trasforms the th-order dfferetal equato to frstorder dfferetal equatos. Cosder a th -order ODE of the form,,, d d d d d d f d d 7.4 wth the followg tal codtos:,,, d d d d d d 7.5 The coverso s a three step process. The frst step s defg ew varables, whch alwas have the followg form: 3 d d d d d d 7.6 The secod step s wrtg oe frst-order ODE for each of these varables. The frst - ODEs have a trval form. The last ODE s smpl the orgal hgher-order ODE, equato 7.4, wth the ew varables, equato 7.6 substtuted.
8 Ordar Dfferetal Equatos - d d d 3 d 7.6 d d d f,, d, 3 The thrd ad fal step of the coverso process s to rewrte the tal codtos, equato 7.5, terms of the ew varables, 7.7,, 3, After the coverso equatos 7.6 ad 7.7 have costtute a sstem of frst-order ODEs wth tal codtos. Ths IVP ca be solved usg the methods of the prevous secto Boudar Value Problems Whe oe vokes Fck s law of dffuso a materal balace, or Fourer s law of heat coducto a eerg balace, or Newto s law of vscost a mometum balace, oe ca up wth a stead-state model that s a secod order ODE, where the depet varable s posto. What dstgushes these problems from other hgher-order ODEs s that the codtos are ot gve at the same pot. From a mathematcal pot of vew, we ca wrte ths equato as d d f,, d d 7.8 wth the followg boudar codtos: ad f f 7.9 Istead of beg gve the value of the fucto ad ts dervatve at a sgle pot,, we ow have a costrat o the fucto at two pots, ad f. For eample, we kow the temperature at ether sde of a metal rod,where our thermocouples are attached, but we kow othg about the
9 Ordar Dfferetal Equatos - dervatve of the temperature. The boudar codtos equato 7.9 defe what s called the Boudar Value Problem BVP. As preseted, the Euler or RK4 ol work for IVPs, ot BVPs. The soluto les creatg a umercal method that combes our ablt to umercall solve olear algebrac equatos d wth ODEs. If we had bee gve a tal value of the frst dervatve,, we d would have a IVP ad we could solve the problem easl, usg the techques Secto 7.6 breakg up the secod-order ODE to two frst-order ODEs ad Secto 7.5 solvg a sstem of frst-order ODEs a IVP. Therefore, we pla to our stregths. We wll guess a value of. We wll solve the sstem of ODEs from to f. We wll compare the calculated value of f wth the boudar codto, f. If the match wth a tolerace, we made a good guess. Probabl, the do t match ad we must make a ew guess for ad terate aga. What should gude our guess for? Wh, I have a ust the thg! We should use the Newto-Raphso method wth umercal dervatves for solvg a sgle olear algebrac equato. The ukow s. The algebrac equato s f g 7. f I other words, we must choose the correct tal codto o the dervatve to satsf the gve boudar codto o the fucto. We uderstad that sce f s arrved at va the soluto of the set of ODEs, that chagg the value of wll chage the value of f. Ever evaluato of the fucto gve equato 7. requres a soluto of the sstem of ODEs. A MATLAB code whch mplemets the soluto of a BVP for a sstem of ODEs usg the Newto-Raphso method ad the classcal fourth-order Ruge-Kutta method s provded later ths chapter Subroute Codes I ths secto, we provde routes for mplemetg the varous umercal ODE solver methods descrbed above. Note that these codes correspod to the theor ad otato eactl as lad out ths book. These codes do ot cota etesve error checkg, whch would complcate the codg ad defeat ther purpose as learg tools. That sad, these codes work ad ca be used to solve problems. As before, o the course webste, two etrel equvalet versos of ths code are provded ad are ttled code.m ad code_short.m. The short verso s preseted here. The loger verso, cotag structos ad servg more as a learg tool, s ot preseted here. The umercal mechacs of the two versos of the code are detcal.
10 Ordar Dfferetal Equatos - 3 Code 7.. Euler Method ODE euler_short fucto [,]=euler_short,o,f,o; d = f-o/; = zeros+,; for = ::+ = o + -*d; = zeros+,; = o; for = :: dd = fukeval,; + = + d*dd; close all; plot = ; f plot == plot,,'k-o', label '', label '' ; fd = fope'euler_out.tt','w'; fprtffd,' \'; fprtffd,'%3.5e %3.5e \', [,]'; fclosefd; fucto dd = fukeval,; dd = -.*; A eample of usg euler_short s gve below.» [,]=euler_short,,,; Ths program geerates outputs three forms. Frst, the ad vectors are stored memor ad ca be drectl accessed. Secod, the program geerates a plot of vs.. Thrd, the program geerates a output fle, euler_out.tt, that cotas ad vectors tabulated form. Code 7.. Fourth-Order Ruge-Kutta Method ODE rk4_short fucto [,]=rk4_short,o,f,o; d = f-o/; = zeros+,; for = ::+ = o + -*d; = zeros+,; = o; for = :: = ; = ; k = fukeval,; = +.5*d;
11 = +.5*d*k; k = fukeval,; 3 = +.5*d; 3 = +.5*d*k; k3 = fukeval3,3; 4 = + d; 4 = + d*k3; k4 = fukeval4,4; dd =./6.*k +.*k +.*k3 + k4; + = + d*dd; close all; plot = ; f plot == plot,,'k-o', label '', label '' ; fd = fope'rk4_out.tt','w'; fprtffd,' \'; fprtffd,'%3.5e %3.5e \', [,]'; fclosefd; fucto dd = fukeval,; dd = -.* A eample of usg rk4_short s gve below.» [,]=rk4_short,,,; Ordar Dfferetal Equatos - 4 Ths commad geerates outputs three forms. Frst, the ad vectors are stored memor ad ca be drectl accessed. Secod, the program geerates a plot of vs.. Thrd, the program geerates a output fle, rk4_out.tt, that cotas ad vectors tabulated form. Code 7.3. Euler Method ODEs euler_short fucto [,]=euler,o,f,o; d = f-o/; = zeros+,; for = ::+ = o + -*d; m=maszeo; = zeros+,m;,:m = o:m; dd = zerosm,; for = :: dd = fukeval,,:m; +,:m =,:m + d*dd:m; close all; plot = ; f plot ==
12 Ordar Dfferetal Equatos - 5 for = ::m color_de = get_plot_color; plot :,:,,color_de; hold o; hold off; label '' ; label '' ; leg tstr[:m]'; fd = fope'euler_out.tt','w'; fprtffd,'... m \'; for = ::+ fprtffd,'%3.5e ', ; for = ::m fprtffd,'%3.5e ',,; fprtffd,' \'; fclosefd; fucto dd = fukeval,; dd = -.* -.* -.5*3; dd = -.* - 4.* -.5*3; dd3 = -.5* -.4* -.*3; A eample of usg euler_short s gve below.» [,]=euler,,,[,,]; Note that the fourth put argumet, o, s ow a vector. Ths program geerates outputs three forms. Frst, the vector ad matr are stored memor ad ca be drectl accessed. Secod, the program geerates a plot of vs.. Thrd, the program geerates a output fle, euler_out.tt, that cotas ad tabulated form. Also ote that ths code calls a acllar fucto for plottg, get_plot_color, whch s reproduced below. The fucto assgs a dfferet color to each fucto so that the curves are dstgushable o the graph. Ths fucto ca be cluded at the bottom of the fle euler_short.m. fucto color_de = get_plot_color; f == color_de = 'k-'; elsef == color_de = 'r-'; elsef == 3 color_de = 'b-'; elsef == 4 color_de = 'g-'; elsef == 5
13 Ordar Dfferetal Equatos - 6 color_de = 'm-'; elsef == 6 color_de = 'k:'; elsef == 7 color_de = 'r:'; elsef == 8 color_de = 'b:'; elsef == 9 color_de = 'g:'; elsef == color_de = 'm:'; else color_de = 'k-'; Code 7.4. Classcal Fourth-Order Ruge-Kutta Method ODEs rk4_short fucto [,]=rk4_short,o,f,o; d = f-o/; = zeros+,; for = ::+ = o + -*d; m=maszeo; = zeros+,m;,:m = o:m; dd = zeros,m; temp = zeros,m; k = zeros,m; k = zeros,m; k3 = zeros,m; k4 = zeros,m; for = :: = ; temp:m =,:m; k:m = fukeval,temp; = +.5*d; temp:m =,:m +.5*d*k:m; k:m = fukeval,temp; 3 = +.5*d; temp:m =,:m +.5*d*k:m; k3:m = fukeval3,temp; 4 = + d; temp:m =,:m + d*k3:m; k4:m = fukeval4,temp; dd:m =./6.*k:m +.*k:m +.*k3:m + k4:m; +,:m =,:m + d*dd:m; close all; plot = ; f plot == for = ::m
14 Ordar Dfferetal Equatos - 7 color_de = get_plot_color; plot :,:,,color_de; hold o; hold off; label '' ; label '' ; leg tstr[:m]'; fd = fope'rk4_out.tt','w'; fprtffd,'... m \'; for = ::+ fprtffd,'%3.5e ', ; for = ::m fprtffd,'%3.5e ',,; fprtffd,' \'; fclosefd; fucto dd = fukeval,; dd = -.* -.* -.5*3; dd = -.* - 4.* -.5*3; dd3 = -.5* -.4* -.*3; A eample of usg rk4_short s gve below.» [,]=rk4,,,[,,]; Note that the fourth put argumet, o, s ow a vector. Ths program geerates outputs three forms. Frst, the vector ad matr are stored memor ad ca be drectl accessed. Secod, the program geerates a plot of vs.. Thrd, the program geerates a output fle, euler_out.tt, that cotas ad tabulated form. Also ote that ths code calls a acllar fucto for plottg, get_plot_color, whch s reproduced above. The fucto assgs a dfferet color to each fucto so that the curves are dstgushable o the graph. Ths fucto ca be cluded at the bottom of the fle rk4_short.m. Code 7.5. Newto-Raphso/Ruge-Kutta Method ODEs BVP rd & rk4 No ew codes were used to solve the BVP problem. Istead the put fles for the Newto- Raphso wth umercal dervatves rd.m ad the classcal fourth-order Ruge-Kutta method for a sstem of ODEs rk4.m were modfed. I the fle rk4.m, we etered the ODEs fucto dd = fukeval,; dd = ; dd = -.* -.5*;
15 Ordar Dfferetal Equatos - 8 I the fle rd.m, the put fucto was specfed as fucto f = fukeval o = ; o_ = ; o_ = ; f =.; f =.; = ; [,]=rk4,o,f,[o_,o_]; f_calc = +,; f = f_calc-f; I ths put fle, we call the Ruge-Kutta route. We provde the gve boudar codtos. We allow the tal codto for the secod fucto to var wth each terato, sce t s the ukow the Newto-Raphso procedure. A eample of usg solvg ths BVP s gve below.» [,err] = rd. cout = old =.e- f = -.45e+ df = 5.85e- ew =.836e+ err =.e+ cout = old =.836e+ f = e-3 df = 5.85e- ew =.836e+ err =.374e-3 =.83 err =.37e-3 The problem was solved oe terato because the ODEs were lear. It took a secod terato to detf covergece. The tal value of the dervatve requred to solve the boudar codto s.83. The rk4.m program geerates a plot of vs. ad geerates a output fle, rk4_out.tt, that cotas ad tabulated form. However, be wared that ths plot ad table of data correspod to the last call of the Ruge-Kutta route b the Newto-Raphso code, whch would have bee to evaluate the umercal dervatve. Therefore, ths plot does ot correspod eactl to the soluto. To geerate the soluto, ru the Ruge-Kutta code wth the approprate tal codto.» [,]=rk4,,,[,.83]; The resultg plot s show Fgure 7.3. It demostrates that the fal value of s deed the specfed value of.
16 Ordar Dfferetal Equatos Problems Problems are located o course webste Fgure 7.3. Soluto to boudar value problem. d d subect to the d d boudar codtos ad. The two curves correspod to black ad the frst dervatve red.
UNIT 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 information( ) 2 2. Multi-Layer Refraction Problem Rafael Espericueta, Bakersfield College, November, 2006
Mult-Layer Refracto Problem Rafael Espercueta, Bakersfeld College, November, 006 Lght travels at dfferet speeds through dfferet meda, but refracts at layer boudares order to traverse the least-tme path.
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 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 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 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 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 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 information7.0 Equality Contraints: Lagrange Multipliers
Systes Optzato 7.0 Equalty Cotrats: Lagrage Multplers Cosder the zato of a o-lear fucto subject to equalty costrats: g f() R ( ) 0 ( ) (7.) where the g ( ) are possbly also olear fuctos, ad < otherwse
More information1 0, x? x x. 1 Root finding. 1.1 Introduction. Solve[x^2-1 0,x] {{x -1},{x 1}} Plot[x^2-1,{x,-2,2}] 3
Adrew Powuk - http://www.powuk.com- Math 49 (Numercal Aalyss) Root fdg. Itroducto f ( ),?,? Solve[^-,] {{-},{}} Plot[^-,{,-,}] Cubc equato https://e.wkpeda.org/wk/cubc_fucto Quartc equato https://e.wkpeda.org/wk/quartc_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 informationSupervised learning: Linear regression Logistic regression
CS 57 Itroducto to AI Lecture 4 Supervsed learg: Lear regresso Logstc regresso Mlos Hauskrecht mlos@cs.ptt.edu 539 Seott Square CS 57 Itro to AI Data: D { D D.. D D Supervsed learg d a set of eamples s
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 information1 Onto functions and bijections Applications to Counting
1 Oto fuctos ad bectos Applcatos to Coutg Now we move o to a ew topc. Defto 1.1 (Surecto. A fucto f : A B s sad to be surectve or oto f for each b B there s some a A so that f(a B. What are examples of
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 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 information0/1 INTEGER PROGRAMMING AND SEMIDEFINTE PROGRAMMING
CONVEX OPIMIZAION AND INERIOR POIN MEHODS FINAL PROJEC / INEGER PROGRAMMING AND SEMIDEFINE PROGRAMMING b Luca Buch ad Natala Vktorova CONENS:.Itroducto.Formulato.Applcato to Kapsack Problem 4.Cuttg Plaes
More informationNumerical Solutions of Second Order Initial Value Problems of Bratu-Type equation Using Higher Ordered Rungu-Kutta Method
Iteratoal Joural of Scetfc ad Research Publcatos Volume 7 Issue October 7 87 ISSN 5-353 Numercal Solutos of Secod Order Ital Value Problems of Bratu-Tpe equato Usg Hgher Ordered Rugu-Kutta Method Habtamu
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 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 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 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 informationApplied Fitting Theory VII. Building Virtual Particles
Appled Fttg heory II Paul Avery CBX 98 38 Jue 8, 998 Apr. 7, 999 (rev.) Buldg rtual Partcles I Statemet of the problem I may physcs aalyses we ecouter the problem of mergg a set of partcles to a sgle partcle
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 informationCS 1675 Introduction to Machine Learning Lecture 12 Support vector machines
CS 675 Itroducto to Mache Learg Lecture Support vector maches Mlos Hauskrecht mlos@cs.ptt.edu 539 Seott Square Mdterm eam October 9, 7 I-class eam Closed book Stud materal: Lecture otes Correspodg chapters
More informationG S Power Flow Solution
G S Power Flow Soluto P Q I y y * 0 1, Y y Y 0 y Y Y 1, P Q ( k) ( k) * ( k 1) 1, Y Y PQ buses * 1 P Q Y ( k1) *( k) ( k) Q Im[ Y ] 1 P buses & Slack bus ( k 1) *( k) ( k) Y 1 P Re[ ] Slack bus 17 Calculato
More informationLecture Notes 2. The ability to manipulate matrices is critical in economics.
Lecture Notes. Revew of Matrces he ablt to mapulate matrces s crtcal ecoomcs.. Matr a rectagular arra of umbers, parameters, or varables placed rows ad colums. Matrces are assocated wth lear equatos. lemets
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 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 informationChapter 5. Curve fitting
Chapter 5 Curve ttg Assgmet please use ecell Gve the data elow use least squares regresso to t a a straght le a power equato c a saturato-growthrate equato ad d a paraola. Fd the r value ad justy whch
More informationMaximum Likelihood Estimation
Marquette Uverst Maxmum Lkelhood Estmato Dael B. Rowe, Ph.D. Professor Departmet of Mathematcs, Statstcs, ad Computer Scece Coprght 08 b Marquette Uverst Maxmum Lkelhood Estmato We have bee sag that ~
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 informationEvaluating Polynomials
Uverst of Nebraska - Lcol DgtalCommos@Uverst of Nebraska - Lcol MAT Exam Expostor Papers Math the Mddle Isttute Partershp 7-7 Evaluatg Polomals Thomas J. Harrgto Uverst of Nebraska-Lcol Follow ths ad addtoal
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 information2/20/2013. Topics. Power Flow Part 1 Text: Power Transmission. Power Transmission. Power Transmission. Power Transmission
/0/0 Topcs Power Flow Part Text: 0-0. Power Trassso Revsted Power Flow Equatos Power Flow Proble Stateet ECEGR 45 Power Systes Power Trassso Power Trassso Recall that for a short trassso le, the power
More informationSolving Constrained Flow-Shop Scheduling. Problems with Three Machines
It J Cotemp Math Sceces, Vol 5, 2010, o 19, 921-929 Solvg Costraed Flow-Shop Schedulg Problems wth Three Maches P Pada ad P Rajedra Departmet of Mathematcs, School of Advaced Sceces, VIT Uversty, Vellore-632
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 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 informationAlgorithms behind the Correlation Setting Window
Algorths behd the Correlato Settg Wdow Itroducto I ths report detaled forato about the correlato settg pop up wdow s gve. See Fgure. Ths wdow s obtaed b clckg o the rado butto labelled Kow dep the a scree
More informationComparison of Dual to Ratio-Cum-Product Estimators of Population Mean
Research Joural of Mathematcal ad Statstcal Sceces ISS 30 6047 Vol. 1(), 5-1, ovember (013) Res. J. Mathematcal ad Statstcal Sc. Comparso of Dual to Rato-Cum-Product Estmators of Populato Mea Abstract
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 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 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 informationRegression and the LMS Algorithm
CSE 556: Itroducto to Neural Netorks Regresso ad the LMS Algorthm CSE 556: Regresso 1 Problem statemet CSE 556: Regresso Lear regresso th oe varable Gve a set of N pars of data {, d }, appromate d b a
More informationBinary classification: Support Vector Machines
CS 57 Itroducto to AI Lecture 6 Bar classfcato: Support Vector Maches Mlos Hauskrecht mlos@cs.ptt.edu 539 Seott Square CS 57 Itro to AI Supervsed learg Data: D { D, D,.., D} a set of eamples D, (,,,,,
More informationGiven a table of data poins of an unknown or complicated function f : we want to find a (simpler) function p s.t. px (
Iterpolato 1 Iterpolato Gve a table of data pos of a ukow or complcated fucto f : y 0 1 2 y y y y 0 1 2 we wat to fd a (smpler) fucto p s.t. p ( ) = y for = 0... p s sad to terpolate the table or terpolate
More informationMechanics of Materials CIVL 3322 / MECH 3322
Mechacs of Materals CVL / MECH Cetrods ad Momet of erta Calculatos Cetrods = A = = = A = = Cetrod ad Momet of erta Calculatos z= z A = = Parallel As Theorem f ou kow the momet of erta about a cetrodal
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 informationIdeal multigrades with trigonometric coefficients
Ideal multgrades wth trgoometrc coeffcets Zarathustra Brady December 13, 010 1 The problem A (, k) multgrade s defed as a par of dstct sets of tegers such that (a 1,..., a ; b 1,..., b ) a j = =1 for all
More informationSimple Linear Regression
Statstcal Methods I (EST 75) Page 139 Smple Lear Regresso Smple regresso applcatos are used to ft a model descrbg a lear relatoshp betwee two varables. The aspects of least squares regresso ad correlato
More informationChapter 2 Supplemental Text Material
-. Models for the Data ad the t-test Chapter upplemetal Text Materal The model preseted the text, equato (-3) s more properl called a meas model. ce the mea s a locato parameter, ths tpe of model s also
More informationLecture 5: Interpolation. Polynomial interpolation Rational approximation
Lecture 5: Iterpolato olyomal terpolato Ratoal appromato Coeffcets of the polyomal Iterpolato: Sometme we kow the values of a fucto f for a fte set of pots. Yet we wat to evaluate f for other values perhaps
More informationMMJ 1113 FINITE ELEMENT METHOD Introduction to PART I
MMJ FINITE EEMENT METHOD Cotut requremets Assume that the fuctos appearg uder the tegral the elemet equatos cota up to (r) th order To esure covergece N must satsf Compatblt requremet the fuctos must have
More informationSTK4011 and STK9011 Autumn 2016
STK4 ad STK9 Autum 6 Pot estmato Covers (most of the followg materal from chapter 7: Secto 7.: pages 3-3 Secto 7..: pages 3-33 Secto 7..: pages 35-3 Secto 7..3: pages 34-35 Secto 7.3.: pages 33-33 Secto
More information1. A real number x is represented approximately by , and we are told that the relative error is 0.1 %. What is x? Note: There are two answers.
PROBLEMS A real umber s represeted appromately by 63, ad we are told that the relatve error s % What s? Note: There are two aswers Ht : Recall that % relatve error s What s the relatve error volved roudg
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 informationSupport vector machines II
CS 75 Mache Learg Lecture Support vector maches II Mlos Hauskrecht mlos@cs.ptt.edu 539 Seott Square Learl separable classes Learl separable classes: here s a hperplae that separates trag staces th o error
More informationJohns Hopkins University Department of Biostatistics Math Review for Introductory Courses
Johs Hopks Uverst Departmet of Bostatstcs Math Revew for Itroductor Courses Ratoale Bostatstcs courses wll rel o some fudametal mathematcal relatoshps, fuctos ad otato. The purpose of ths Math Revew s
More informationAnalysis of Lagrange Interpolation Formula
P IJISET - Iteratoal Joural of Iovatve Scece, Egeerg & Techology, Vol. Issue, December 4. www.jset.com ISS 348 7968 Aalyss of Lagrage Iterpolato Formula Vjay Dahya PDepartmet of MathematcsMaharaja Surajmal
More informationHomework Assignment Number Eight Solutions
D Keer MSE 0 Dept o Materals Scece & Egeerg Uversty o Teessee Kovlle Homework Assgmet Number Eght Solutos Problem Fd the soluto to the ollowg system o olear algebrac equatos ear () Soluto: s Sce ths s
More informationCorrelation and Regression Analysis
Chapter V Correlato ad Regresso Aalss R. 5.. So far we have cosdered ol uvarate dstrbutos. Ma a tme, however, we come across problems whch volve two or more varables. Ths wll be the subject matter of the
More informationLecture 1: Introduction to Regression
Lecture : Itroducto to Regresso A Eample: Eplag State Homcde Rates What kds of varables mght we use to epla/predct state homcde rates? Let s cosder just oe predctor for ow: povert Igore omtted varables,
More informationJohns Hopkins University Department of Biostatistics Math Review for Introductory Courses
Johs Hopks Uverst Departmet of Bostatstcs Math Revew for Itroductor Courses Ratoale Bostatstcs courses wll rel o some fudametal mathematcal relatoshps, fuctos ad otato. The purpose of ths Math Revew s
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 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 informationhp calculators HP 30S Statistics Averages and Standard Deviations Average and Standard Deviation Practice Finding Averages and Standard Deviations
HP 30S Statstcs Averages ad Stadard Devatos Average ad Stadard Devato Practce Fdg Averages ad Stadard Devatos HP 30S Statstcs Averages ad Stadard Devatos Average ad stadard devato The HP 30S provdes several
More informationENGI 4430 Numerical Integration Page 5-01
ENGI 443 Numercal Itegrato Page 5-5. Numercal Itegrato I some o our prevous work, (most otaly the evaluato o arc legth), t has ee dcult or mpossle to d the dete tegral. Varous symolc algera ad calculus
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 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 informationLecture Notes Forecasting the process of estimating or predicting unknown situations
Lecture Notes. Ecoomc Forecastg. Forecastg the process of estmatg or predctg ukow stuatos Eample usuall ecoomsts predct future ecoomc varables Forecastg apples to a varet of data () tme seres data predctg
More informationPHYS Look over. examples 2, 3, 4, 6, 7, 8,9, 10 and 11. How To Make Physics Pay PHYS Look over. Examples: 1, 4, 5, 6, 7, 8, 9, 10,
PHYS Look over Chapter 9 Sectos - Eamples:, 4, 5, 6, 7, 8, 9, 0, PHYS Look over Chapter 7 Sectos -8 8, 0 eamples, 3, 4, 6, 7, 8,9, 0 ad How To ake Phscs Pa We wll ow look at a wa of calculatg where the
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 informationObjectives of Multiple Regression
Obectves of Multple Regresso Establsh the lear equato that best predcts values of a depedet varable Y usg more tha oe eplaator varable from a large set of potetal predctors {,,... k }. Fd that subset of
More informationSummary of the lecture in Biostatistics
Summary of the lecture Bostatstcs Probablty Desty Fucto For a cotuos radom varable, a probablty desty fucto s a fucto such that: 0 dx a b) b a dx A probablty desty fucto provdes a smple descrpto of the
More informationOn Modified Interval Symmetric Single-Step Procedure ISS2-5D for the Simultaneous Inclusion of Polynomial Zeros
It. Joural of Math. Aalyss, Vol. 7, 2013, o. 20, 983-988 HIKARI Ltd, www.m-hkar.com O Modfed Iterval Symmetrc Sgle-Step Procedure ISS2-5D for the Smultaeous Icluso of Polyomal Zeros 1 Nora Jamalud, 1 Masor
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 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 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 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 informationOrdinary Least Squares Regression. Simple Regression. Algebra and Assumptions.
Ordary Least Squares egresso. Smple egresso. Algebra ad Assumptos. I ths part of the course we are gog to study a techque for aalysg the lear relatoshp betwee two varables Y ad X. We have pars of observatos
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 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 informationDescriptive Statistics
Page Techcal Math II Descrptve Statstcs Descrptve Statstcs Descrptve statstcs s the body of methods used to represet ad summarze sets of data. A descrpto of how a set of measuremets (for eample, people
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 informationå 1 13 Practice Final Examination Solutions - = CS109 Dec 5, 2018
Chrs Pech Fal Practce CS09 Dec 5, 08 Practce Fal Examato Solutos. Aswer: 4/5 8/7. There are multle ways to obta ths aswer; here are two: The frst commo method s to sum over all ossbltes for the rak of
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 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 informationPower Flow S + Buses with either or both Generator Load S G1 S G2 S G3 S D3 S D1 S D4 S D5. S Dk. Injection S G1
ower Flow uses wth ether or both Geerator Load G G G D D 4 5 D4 D5 ecto G Net Comple ower ecto - D D ecto s egatve sg at load bus = _ G D mlarl Curret ecto = G _ D At each bus coservato of comple power
More information2SLS Estimates ECON In this case, begin with the assumption that E[ i
SLS Estmates ECON 3033 Bll Evas Fall 05 Two-Stage Least Squares (SLS Cosder a stadard lear bvarate regresso model y 0 x. I ths case, beg wth the assumto that E[ x] 0 whch meas that OLS estmates of wll
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 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 informationMultiple Regression. More than 2 variables! Grade on Final. Multiple Regression 11/21/2012. Exam 2 Grades. Exam 2 Re-grades
STAT 101 Dr. Kar Lock Morga 11/20/12 Exam 2 Grades Multple Regresso SECTIONS 9.2, 10.1, 10.2 Multple explaatory varables (10.1) Parttog varablty R 2, ANOVA (9.2) Codtos resdual plot (10.2) Trasformatos
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 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 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 informationX ε ) = 0, or equivalently, lim
Revew for the prevous lecture Cocepts: order statstcs Theorems: Dstrbutos of order statstcs Examples: How to get the dstrbuto of order statstcs Chapter 5 Propertes of a Radom Sample Secto 55 Covergece
More informationStability For a stable numerical scheme, the errors in the initial condition will not grow unboundedly with time.
.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
More informationMean is only appropriate for interval or ratio scales, not ordinal or nominal.
Mea Same as ordary average Sum all the data values ad dvde by the sample sze. x = ( x + x +... + x Usg summato otato, we wrte ths as x = x = x = = ) x Mea s oly approprate for terval or rato scales, ot
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 informationLecture 1: Introduction to Regression
Lecture : Itroducto to Regresso A Eample: Eplag State Homcde Rates What kds of varables mght we use to epla/predct state homcde rates? Let s cosder just oe predctor for ow: povert Igore omtted varables,
More information