Solving the advection PDE in explicit FTCS, Lax, Implicit FTCS and Crank-Nicolson methods for constant and varying speed.
|
|
- Kenneth McBride
- 6 years ago
- Views:
Transcription
1 Solvig te advectio PDE i explicit FTCS, Lax, Implicit FTCS ad Crak-Nicolso metods for costat ad varyig speed. Accuracy, stability ad software aimatio Report submitted for fulfillmet of te Requiremets for MAE 294 Masters degree project Supervisor: Dr Doald Dabdub, UCI. Writte by Nasser M. Abbasi. Masters degree cadidate studet. Mecaical egieerig departmet Uiversity of Califoria, Irvie Cotets 1 Itroductio Backward differece (Upwid) Forward differece (dowwid) Ceter differece Numerical scemes Explicit Metods FTCS Dowwid Upwid LAX Lax-Wedroff Leap-frog Implicit Metods Implicit FTCS Wedrof Crak-Nicolso Stability aalysis Stability aalysis for FTCS Stability aalysis of te dowwid metod Stability aalysis of te upwid metod Stability aalysis of Lax Stability of Lax-Wedroff
2 3.6 Stability aalysis of te Implicit FTCS Solutio Results ad Output Case Case Case Case case case case case CPU compariso tables Accuracy compariso tables Coclusio 24 6 Appedix Plots case case case case case case case case Source code Refereces 48 2
3 1 Itroductio Te goal of tis project is to aalyze ad compare differet umerical metods for solvig te first order advectio PDE. Followig te aalytical aalysis for stability of te umerical sceme, aimatio were doe to visually illustrate ad cofirm tese results. Tis was carried for differet parameters. Te aimatio was programmed i Matematica ad saved to aimated gif files wic was te loaded ito te HTML versio of tis report located ere. Fortra 95 was used for te computatio part, wile Matematica was used for te aimatio ad grapics part. Te above lik cotais all te supportig material for te project, icludig te Fortra program (i source ad widows executable format) used to carry te mai computatio, ad te Matematica program used to do te aimatio ad te Uix bas file used to process te computatio for differet parameters. Te specific PDE example used for te aalysis ad aimatio was te oe provided by Professor Doald Dabdub for te fial exam for is MAE 185 course (Numerical metods for mecaical egieers) i sprig Tis PDE is described below: Solve umerically c t + u c x = 0 Were c(x, t) is te cocetratio of a give material as a fuctio of time ad space. Te above is solved for te followig 2 cases 1. u (te advectio speed, or te speed at wic te mass is beig trasported) is a costat value give as (2 ft/mi. 2. u is a fuctio of time defied as u(t) = t 20 ft/mi Te problem parameters are: t 0 0 x L Were L = 100 feet.. Iitial coditios are { ( ( 1 + cos π x 30 )) 5 25 x 25 c (x, 0) = F (x) = 0 oterwise Te boudary coditios are c (0, t) = 0 c(l, t) = 0 Tis PDE is a example of a IBVP (Iitial ad Boudary Value Problem). Differet umerical metods are used to solve te above PDE. Te metods are compared for stability usig Vo Neuma stability aalysis. Te umerical metods are also compared for accuracy. Tis was doe by comparig te umerical solutio to te kow aalytical solutio at eac time step. Te compariso was doe by computig te root mea square error (RMSE) betwee te umerical ad te aalytical solutio at eac time step. Te metod wit te least RMSE at te ed of te simulatio is cosidered te most accurate. Te above PDE as a kow aalytical solutio wic is C (x, t) = F (x ut) 3
4 Te above aalytical solutio idicates tat te iitial cocetratio will move from left to rigt wit te advectio speed u. Te formulatio of eac umerical metod is sow below. is used to represet x, te space betwee 2 space grid poit, or te space step size, ad τ is used to represet t, te time step. Te space lie as N grid poits. Te spacig was fixed at 0.01 ft for all te metods ad for all te test cases, wile τ was caged. Tis made comparig te differet metods simpler. Te followig diagram illustrates te discretizatio used. Sould we cosider te lower left ad te lower rigt grid poits above as part of te iitial coditios, or part of te boudary coditios? Stability of eac metod is derived. Stability is importat, sice by te Lax-Rictmyer equivalece teorem 1, stability implies covergece of te solutio. Covergece of te umerical solutios implies tat as te step size becomes smaller, te umerical solutio coverges to te aalytical solutio. Explicit ad implicit umerical metods are used. We solvig for te future value of te solutio at a sigle ode i terms of oly past values, te metod is called a explicit metod. I oter words, we te oly ukow is te future value of te solutio at a sigle ode, ad everytig else o te rigt ad side of te fiite differece equatio is a solutio derived at earlier time step, te metod is explicit. A implicit metod is oe i wic te fiite differece equatio cotais te solutio at a at future time at more ta oe ode. I oter words, future solutio are beig solved for at more ta oe ode i terms of te solutio at earlier time. Implicit metods terefor are usually solved by matrix metods by solvig Ax = b were b represets preset preset kow solutio values, ad x are te ukow future solutio values, ad A is te coefficiet matrix wic will usually be block diagoal (or tri diagoal) i sape. I te derivatios below, te otatio of Ci is used to idicate te solutio at time step ad at space ode i. Hece C (x i, t ) is writte as Ci. Tis otatio seems to be more clear ta te C i, otatio. Differet fiite differece scemes for solvig a PDE are obtaied by usig differet metods of approximatig te derivative terms i te PDE. Tis will be illustrated usig te space derivative c x. Tis derivative ca be approximated i oe of te followig 3 ways (all at time step ) 1.1 Backward differece (Upwid) c x C i Ci 1 1 Rictmyer ad Morto p45): Give a properly posed liear iitial value problem ad a fiite differece approximatio to it tat satisfies te cosistecy coditio, stability is te ecessary ad sufficiet coditio for covergece. 4
5 1.2 Forward differece (dowwid) c x C i+1 C i 1.3 Ceter differece c x C i+1 C i 1 2 Te followig are te derivatio of a umber of metods for solvig te advectio PDE obtaied by usig te above defiitios for te derivative we applied to bot space ad time. 5
6 2 Numerical scemes 2.1 Explicit Metods FTCS Wit FTCS, te forward time derivative, ad te cetered space derivative are used. Hece te advectio PDE ca be writte as Solvig for C +1 i results i C +1 i C +1 i Ci τ ( C = u i+1 Ci 1 ) 2 = Ci uτ ( C 2 i+1 Ci 1) Tis metod will be sow to be ucoditioally ustable Dowwid Here, te forward time derivative for C t is used ad also te forward space derivative for C x. Tis results i ) C +1 i C +1 i Ci τ ( C = u i+1 Ci = C i uτ ( C i+1 Ci ) Tis metod will be sow to be ucoditioally ustable as well Upwid Here, te forward time derivative for C t is used, ad te backward derivative for C x is used. Tis results i ( C i Ci 1 ) C +1 i C +1 i Tis will be sow to be stable if uτ LAX Ci τ = u = C i uτ Lookig at te FTCS eq (1) above, ad sow below agai Te term C i C +1 i ( C i C i 1) = Ci uτ ( C 2 i+1 Ci 1) above is replaced by its average value C i+1 +C i 1 2 to obtai te LAX metod (0) (1) Tis metod will be sow to be stable if uτ 1 Ci +1 = 1 ( C 2 i+1 + Ci 1) uτ ( C 2 i+1 Ci 1) (4) 6
7 2.1.5 Lax-Wedroff By usig te secod-order fiite differece sceme for te time derivative, te metod of Lax-Wedroff metod is obtaied Leap-frog C +1 i = Ci uτ ( C 2 i+1 Ci 1) u 2 τ 2 ( + C 2 2 i+1 + Ci 1 2Ci ) I tis metod, te cetered derivative is used for bot time ad space. Tis results i C +1 i Ci 1 2τ ( C = u i+1 Ci 1 ) 2 Tis metod requires a special startig procedure due to te term Ci 1. Aoter sceme suc as Lax ca be used to kick start tis metod. 2.2 Implicit Metods Implicit FTCS Give te explicit FTCS derived above C +1 i Ci τ ( C = u i+1 Ci 1 ) 2 Te above is modified it by evaluatig te space ceter derivative at time step + 1 istead of at time step, tis results i ( Ci +1 Ci C +1 i+1 = u ) C+1 i 1 (5A) τ 2 Hece C +1 i Writig it i matrix form, first lettig α = uτ 2 results i α 1 α α 1 α α 1 α 0. + uτ 2 C+1 i+1 uτ 2 C+1 i 1 = C i (5B) Were N is te umber of space grid poits. Te above is writte as Ax = b C +1 0 C +1 1 C +1 2 C C +1 N 1 = Solvig for x, wic represets te solutio at time step + 1 or at time t = ( + 1) τ. b represets te curret solutio at time step, ad A is te matrix of te coefficiets sow above. Due to te form of te A matrix, (Called tri diagoal, or Block diagoal), a algoritm tat takes advatages of tis form is used. Tis is called te Tomas algoritm. Tis greatly speeds up te solutio. If we ad used a geeral algoritm to solve tis system suc as te Gauss elimiatio metod, it would ave bee muc slower, C 0 C 1 C 2 C 3. C N 1 7
8 makig te implicit metod ot practical to use. (Some tests o te same data sowed te Tomas algoritm to be 50 times faster ta Gaussia elimiatio) Wedrof Tis metod uses ceter differece for te derivative aroud te space step ( i + 1 2) ad te time step ( + 1 2) τ Tis leads to te followig sceme ( 1 uτ ) ( Ci uτ ) C +1 i+1 = ( 1 + uτ ) ( Ci + 1 uτ ) C i+1 Tis ca also be solved usig similar matrix metod to tat used for te implicit FTCS. Tis metod is ot used i tis report Crak-Nicolso By takig te average of te explicit FTCS ad te implicit FTCS formulatios (sow agai below), te C-N sceme is derived C +1 i C +1 i Ci τ Ci τ ( C = u i+1 Ci 1 ) 2 ( C +1 i+1 = u ) C+1 i 1 2 Takig te average of te above results i Ci +1 Ci = u ( C i+1 C ) ( i 1 u C +1 i+1 ) C+1 i 1 τ as C +1 i + uτ 4 C+1 i+1 uτ 4 C+1 i 1 = C i uτ 4 C i+1 + uτ 4 C i 1 Now te system Ax = b is setup to solve for future values as follows. Let α = uτ 4, te system ca be writte α 1 α α 1 α α 1 α α C +1 0 C +1 1 C +1 2 C C +1 N 1 Tomas algoritm is used to solve te above system for C +1 i. C 0 C1 αc 2 + αc 0 C = 2 αc 3 + αc 1 C3 αc 4 + αc 2. C N 1 8
9 3 Stability aalysis A umerical solutio is stable if te eergy cotet remai below some limitig value o matter ow log te solutio is itegrated. I essece, tis meas tat te solutio does ot blow up after some time. Tis ca be called BIBO stability (Bouded I Bouded Out). Hece oe way to aalyze te stability of te umerical solutio is to determie a expressio tat relates te amplitude of te solutio betwee 2 time steps, ad to determie if tis ratio remai less ta or equal to a uity as more ad more time steps are take. Tis type of aalysis is called Vo Neuma stability aalysis for umerical metods. Te aalysis is based of fidig a expressio for te magificatio factor of te wave amplitude at eac step. Te solutio will be stable if tis magificatio factor is less ta oe. Let te magificatio factor be ξ. Te umerical sceme is stable iff ζ 1 Te Courat Friedrics Lewy (CFL) criteria for stability says tat ζ 1 uτ 1 Were u,, ad τ are as defied above: u is te wave speed, = x ad τ = t. Te umber uτ is also called te courat umber. Some umerical metods will be sow to be ucoditioally ustable (suc as explicit FTCS ad te explicit upwid). Tis meas tat eve if courat umber was 1, te umerical solutios will evetually become ustable. Some explicit metods suc as LAX, are coditioally stable if te courat umber was 1. Implicit metods are ucoditioally stable, ece courat umber is ot used for tese metods. However, tis does ot mea oe ca take as large step as oe wats wit te implicit metods, sice accuracy will be affected eve if te solutio remai stable. Hece, te best umerical sceme is oe i wic te largest step size ca be take, wit te least amout of iaccuracy i te umerical solutio wile remaiig stable. For umerical sceme tat are coditioally stable, it ca be see from te CFL coditio tat for a fixed speed u ad fixed, te maximum time step tat ca be take is give by τ max u It ca be immediately see from above, tat for te case we te advectio speed is varyig ad is a fuctio of time suc as te case we u (t) = t 20 implyig tat te speed is icreasig wit time, te we usig a fixed time step τ it will evetually become larger ta u ad te umerical sceme will be ustable. Tis is because as u (t) is becomig larger ad larger, wile is fixed, te term u will become smaller ad smaller. Hece to keep te courat umber uτ 1, te time step take must remai less ta u, ece usig a fixed time step wit icreasig u will evetually lead to istability. Tis will affect te explicit metods tat are coditioally stable suc as te LAX metod, sice te Lax metod is explicit ad depeds o satisfyig te CFL all te time for its stability. Implicit metods are stable for ay time step. I te followig we derive te details of te stability aalysis ad use Vo Neuma aalysis to derive a expressio for te amplificatio factor ζ for differet umerical scemes. So to summarize: 1. Explicit FTCS is ucoditioally ustable. 2. Explicit LAX is stable if uτ 3. Implicit FTCS ad C-R are stable for all τ 1, or i oter words, τ max u 9
10 3.1 Stability aalysis for FTCS Usig Vo Neuma metod, te followig trial solutio to te PDE is assumed c (x, t) = A (t) e jkx were j = 1 ad k is te wave umber ad A is te amplitude of te wave, as a fuctio of time. Hece te solutio at time step ad at x = x i = i is writte as A e jki (2) Substitute tis trial solutio (2) ito te (1) results i A +1 e jki = A e jki uτ ( A e jk(i+1) A e jk(i 1)) (3) 2 Let ξ be te ratio of te amplitude of te wave at time step + 1 relative to tat at time step. ece Divide (3) by A results i Divide te above by e jki ξ = A+1 A ξe jki = e jki uτ (e jk(i+1) e jk(i 1)) 2 ξ = 1 uτ (e jk e jk) 2 = 1 uτ j si (k) Hece ( uτ ) 2 ξ = 1 + si (k) Tis implies tat ξ 1 regardless of te time step τ selected or te space step, ece FTCS is ucoditioally ustable. For a fixed speed u, te istability ca be delayed by makig τ smaller, but it could ot be preveted. Evetually tis umerical solutio will blow up. Tis will be illustrated below i a aimatio. See case 3 ad 4 as examples. Te istability ca be delayed by makig τ smaller for a fixed, or by makig larger for a fixed τ. 3.2 Stability aalysis of te dowwid metod Ci +1 = Ci uτ ( C i+1 Ci ) Substitute te trial solutio A e jki ito te above A +1 e jki = A e jki uτ (A e jk(i+1) A e jki) ξ = 1 uτ ( ) e jk 1 = 1 + uτ uτ ejk = 1 + uτ uτ (cos (k) + j si (k)) = 1 + uτ (1 cos k) j uτ 10 si k
11 Let uτ = λ Hece ξ = 1 + λ (1 cos k) jλ si k ξ 2 = (1 + λ (1 cos k)) 2 + (λ si k) 2 = 1 + 2λ (1 cos k) + λ 2 (1 cos k) 2 + λ 2 si 2 k = 1 + 2λ (1 cos k) + λ 2 ( 1 2 cos k + cos 2 k ) + λ 2 si 2 k = 1 + 2λ 2λ cos k + λ 2 2λ 2 cos k + λ 2 cos 2 k + λ 2 si 2 k = 1 + 2λ 2λ cos k + 2λ 2 2λ 2 cos k = 1 + 2λ (1 + λ) (1 cos k) Hece for stability it is required tat 1 + 2λ (1 + λ) (1 cos k) 1 or 2λ (1 + λ) (1 cos k) 0 sice λ = uτ, a positive quatity, te te above coditio ca ot be satisfied. Hece te dowwid metod is ucoditioally ustable. 3.3 Stability aalysis of te upwid metod Ci +1 = Ci uτ Substitute te trial solutio A e jki ito te above ( C i C i 1) Let uτ Hece Hece = λ A +1 e jki = A e jki uτ (A e jki A e jk(i 1)) ξ = 1 uτ (1 e jk) = 1 uτ + uτ e jk = 1 uτ + uτ (cos (k) j si (k)) = 1 uτ (1 cos k) j uτ ξ = 1 λ (1 cos k) jλ si k si k ξ 2 = (1 λ (1 cos k)) 2 + (λ si k) 2 = 1 2λ (1 cos k) + λ 2 (1 cos k) 2 + λ 2 si 2 k = 1 2λ + 2λ cos k + λ 2 ( 1 + cos 2 k 2 cos k ) + λ 2 si 2 k = 1 2λ + 2λ cos k + λ 2 + λ 2 cos 2 k 2λ 2 cos k + λ 2 si 2 k = 1 2λ + 2λ cos k + 2λ 2 2λ 2 cos k = 1 2λ (1 λ) (1 cos k) 11
12 Hece for stability it is required tat or 1 2λ (1 λ) (1 cos k) 1 2λ (1 λ) (1 cos k) 0 Wic will be true oly if (1 λ) 0 or λ 1 ece tis implies uτ 1 Hece te upwid metod is stable if te CFL coditio is satisfied. Tis will be see as te same stability coditio for te Lax metod below. 3.4 Stability aalysis of Lax Replace te trial fuctio from (2) i Lax formulatio i (4) ad obtai A +1 e jki = 1 2 ( A e jk(i+1) + A e jk(i 1)) uτ (A e jk(i+1) A e jk(i 1)) 2 Divide by A e jki, te magificatio factor ζ is obtaied Hece ζ = 1 ( e jk + e jk) uτ (e jk e jk) 2 2 = cos (k) j uτ si (k) ζ = cos 2 (k) + ( uτ ) 2 si 2 (k) Sice cos 2 (k) 1 ad si 2 (k) 1, te it is see tat ζ 1 if uτ 1 Hece te followig is te coditio for stability τ u As metioed earlier, tis is called te CFL coditio. Te Lax metod is stable for τ u owever, a modified versio of tis metod is more accurate, wic is te Lax-Wedroff metod. 3.5 Stability of Lax-Wedroff Tis is te same as te Lax metod. Te metod is stable if τ u 3.6 Stability aalysis of te Implicit FTCS Replace te trial fuctio from (2) i (5B) results i Divide by A e jki A +1 e jki + uτ 2 A+1 e jk(i+1) uτ 2 A+1 e jk(i 1) = A e jki 12
13 Hece Hece tis sows tat te ξ + uτ 2 ξejk uτ 2 ξe jk = 1 ( ξ 1 + uτ 2 ejk uτ 2 e jk) = 1 ( ξ 1 + j uτ ) si (k) = 1 ξ = 1 ξ = 1 + j uτ si (k) = 1 j uτ si (k) 1 + uτ si (k) 1 + ( ) uτ 2 si 2 (k) 1 + uτ si (k) < 1 Implicit FTCS metod is ucoditioally stable. Tis property is commo to all implicit metods. Eve toug te implicit FTCS is stable, it is ot very accurate. See case 8 below for a example. 13
14 4 Solutio Results ad Output For te Fortra implemetatio, te followig metods are implemeted. Te explicit FTCS, Explicit Lax, Implicit FTCS, ad Implicit Crak-Nicolso. For eac metod, te followig was geerated 1. CPU time used for te ru. 2. sap sot of te solutio at t = 0, t = 15, ad t = 30 miutes. 3. RMSE betwee te umerical solutio ad te aalytical solutio. 4. Aimatio of te umerical solutio. Te aimatio was doe by takig sapsots of te solutio at regular itervals i Fortra. Tese were saved to disk. Te Matematica was used to geerate te aimatio ad te plots. To compare te stability ad accuracy of te metods, te time step was caged (icreased) ad a ew ru was made. 8 differet values of time steps are used. So tere are 8 tests cases. Tese 8 test cases were ru for bot fixed speed (u = 2 ft/mi) ad for u = t 20 ft/mi. Tis table below summarizes tese cases. Te appedix cotais all te plots. Te aimatios are added as HTML liks. 4.1 Case 1 τ = sec, = 0.1 ft Speed Metod CPU time (sec) RMSE Aimatio (2D) plots U=2 Explicit FTCS Explicit LAX Implicit FTCS C-R U=t/20 Explicit FTCS Explicit LAX Implicit FTCS C-R Note te followig: Te explicit FTCS remaied stable trougout te ru due to te small time step. All oter metods were stable as well durig te ru. For te CPU for te varyig u case, otice tat for te implicit metods tis value is larger ta te CPU for te same metods but we u is fixed. Tis is due to te fact tat te matrix A is o loger costat, ad must be recomputed at eac time step before callig Tomas algoritm to solve Ax = b system. Also otice tat te CPU time for te implicit metods is larger ta te explicit metods. Tis is due to te extra computatioal cost i solvig Ax = b. Eve we usig Tomas algoritm, tis is still more expesive ta te explicit metods we umber of time steps is large. 14
15 4.2 Case 2 τ = sec, = 0.1 ft Speed Metod CPU time (sec) RMSE Aimatio (2D) plots U=2 Explicit FTCS Explicit LAX Implicit FTCS C-R U=t/20 Explicit FTCS Explicit LAX Implicit FTCS C-R Te explicit FTCS is stable for most of te ru, ear te ed it is startig to be become ustable. Notice tat aroud 26 miutes tat bubbles are startig to sow up i te umerical solutio dowstream. Tis is a caracteristic of ow tis metod becomes ustable. Tis will be more clear i te ext test cases we te time step is made larger. For te varyig speed case, te explicit metod usig te same time step remaied stable durig te wole 30 miutes. Tis is because te average speed was less ta 2 ft/mi, ece te mass did ot ave to travel as log a distace as wit fixed speed of u = 2, ad so te istability did ot sow up. Matematically tis ca be explaied by lookig at te term uτ, ece for smaller u, te courat umber is smaller. Notice also te RMSE is smaller for variable speed compared to fixed speed. Agai tis is related to te smaller average speed makig te courat umber smaller. 4.3 Case 3 I tis case, we sligtly make te time step loger ta before. We start to see te istability of FTCS. τ = sec, = 0.1 ft, uτ = for fixed u Speed Metod CPU time (sec) RMSE Aimatio (2D) plots U=2 Explicit FTCS Explicit LAX Implicit FTCS C-R U=t/20 Explicit FTCS Explicit LAX Implicit FTCS C-R For explicit FTCS, Te solutio ow startig to sow istability at 25 miutes. Lax remaied stable sice CFL is satisfied. Explicit FTCS is becomig less accurate as well. Explicit Lax is most accurate at tis time step. 4.4 Case 4 I tis case, we sligtly make te time step eve loger ta before. Now FTCS becomes more ustable. τ = sec, = 0.1 ft, uτ =
16 Speed Metod CPU time (sec) RMSE Aimatio (2D) plots U=2 Explicit FTCS Explicit LAX Implicit FTCS C-R U=t/20 Explicit FTCS Explicit LAX Implicit FTCS 4, C-R FTCS Istability starts at aroud 20 miutes. LAX remaied stable sice CFL is satisfied. Lax remaied te most accurate at tis time step. It accuracy actually improved as te time step became larger. 4.5 case 5 Agai te time step is made loger ta before. Now te explicit FTCS is completely ustable. τ = sec, = 0.1 ft For te case of fixed U, we ave uτ = = 0.9 1, wile for varyig U, te maximum value will be at te ed of te ru, wic is 30/20 = 1.5 ft/mi., ece te CFL coditio is cagig, wit a value of = at te ed of te ru wic is still 1 Speed Metod CPU time (sec) RMSE Aimatio (2D) plots U=2 Explicit FTCS 0.73 blows up Explicit LAX Implicit FTCS C-R U=t/20 Explicit FTCS 0.28 blow up Explicit LAX Implicit FTCS C-R For te varyig speed case, te explicit FTCS remaied stable for te duratio of te ru as compared to te case wit te fixed speed. Tis is because te average wave speed is less ta wit te fixed wave speed case. Te magificatio factor depeds o te speed of te wave. ( uτ ) 2 ξ = 1 + si (k) Wit te varyig speed case, te coefficiet uτ was smaller durig te wole ru, sice te maximum speed u attaied will be 1.5 ft/mi. as compared to 2 ft/mi. i te fixed u case. We see ta te smaller te speed u te smaller te magificatio (wit everytig else beig fixed). If we ave ru te simulatio a little loger for te varyig speed case, we will see te istability wit explicit FTCS. Tis below is a diagram sowig 2 rus usig te explicit FTCS bot wit u = t 20 ft/mi, oe was ru for 30 miutes, ad te secod for 53 miutes. Te ru to 30 miutes sowed o istability wile te ru for 53 miutes sowed te istability. Tis sow te explicit FTCS will evetually become ustable. 16
17 Tis is a aimatio of te above 17
18 4.6 case 6 I tis case, te time step is icreased so tat uτ is just above te CFL coditio. Notice ow tat te Explicit LAX metod become ustable as expected. Te oter implicit metods remai stable. te explicit FTCS metod ow is completely ustable. Te implicit FTCS metod is startig to become less accurate. τ = sec, = 0.1 ft, uτ = = > 1 Speed Metod CPU time (sec) RMSE Aimatio (2D) plots U=2 Explicit FTCS 0.7 blows up N/A blows up Explicit LAX Implicit FTCS C-R U=t/20 Explicit FTCS 0.28 blows up N/A blows up Explicit LAX Implicit FTCS C-R Notice tat explicit LAX takes muc less CPU ta ay oter metod. 4.7 case 7 τ = 0.06 sec, = 0.1 ft, uτ = = 1. 2 > 1 Speed Metod CPU time (sec) RMSE Aimatio (2D) plots U=2 Explicit FTCS 0.65 blows up N/A blows up Explicit LAX 0.9 blows up Implicit FTCS C-R U=t/20 Explicit FTCS blows up N/A blows up Explicit LAX Implicit FTCS C-R Notice tat te CPU for te implicit metod we speed is fixed is ow iger ta te CPU for te explicit metods. Tis ca be explaied as follows: sice te time step ow is larger ta before, te umber of times to solve Ax = b as bee reduced. Tis made te implicit metods faster. Tis implies tat Usig a relatively large time step, implicit metods become faster ta te explicit metods. 18
19 4.8 case 8 τ = 0.07 sec, = 0.1 ft, uτ = = 1. 4 > 1 Speed Metod CPU time (sec) RMSE Aimatio (2D) plots U=2 Explicit FTCS 0.5 blows up N/A blows up Explicit LAX 0.89 blows up Implicit FTCS C-R U=t/20 Explicit FTCS blows up N/A blows up Explicit LAX Implicit FTCS C-R
20 4.9 CPU compariso tables As expected, CPU time usage will be less as te time step is icreased. Tere is a aomaly cased oticed were te CPU time icreased for te Lax metod we te time step is icreased from to 0.06, Tis eeds furter ivestigatio. Tis table below summarizes te CPU time i secods used by eac metod for te case of costat speed as time step is icreased. τ sec Explicit F T CS Explicit LAX Implicit F T CS C R Tis is te plot of te above table Tis table below summarizes te CPU time i secods used by eac metod for te case of varyig speed as time step is icreased. τ sec Explicit F T CS Explicit LAX Implicit F T CS C R Tis is te plot of te above table 20
21 Tis plot below compares te CPU time for eac metod we te speed is costat vs. we te speed was cagig wit time. 21
22 22
23 4.10 Accuracy compariso tables Tis table below summarizes te RMS error from eac umerical metod as a fuctio of cagig te time step size. Tis is for case of costat speed. time step Explicit F T CS Explicit LAX Implicit F T CS C R blows up blows up blows up blows up blows up blows up Notice tat te Lax metod became more accurate we te time step was icreased from to 0.04 secods, te it starts to become less accurate as time step is icreased. Tis is couter ituitive to wat oe ca expect. It will be iterestig to ivestigate tis furter to obtai a matematical explaatio for tis strage peomea. Te accuracy of te implicit FTCS, ad C-R also icreased sligtly as te time step became larger from to , te te implicit FTCS became worst i terms of accuracy as te time step icreased. C-R metod accuracy did ot deteriorate as muc wit icreasig te time step. Tis sows te C-R sceme to be more robust. Tis table below summarizes te RMS error from eac umerical metod as a fuctio of cagig te time step size. Tis is for case of cagig speed. time step Explicit F T CS Explicit LAX Implicit F T CS C R blows up blows up blows up blows up Te effect of avig te speed defied as µ = t 20is to delay istability for te explicit metods as time step is icreased. Notice also ere te case were te Lax metod became more accurate as te time step is icreased from to
24 5 Coclusio 4 differet umerical fiite differece scemes are examied for CPU time, stability ad accuracy i solvig te advectio PDE for costat speed ad for a speed wic is a fuctio of time. For accuracy, a iterestig result is observed. Te Lax sceme is te most accurate for Courat umber close to uity. Tis meas as te time step is icreased, te Lax become more accurate of te 4 metods. But beyod te CFL coditio, Bot explicit metods (FTCS ad Lax) became less accurate. Explicit FTCS became ustable sooer ta Lax, wile te implicit metods remaied stable. Te implicit FTCS was less accurate ta te C-R metod. Tis implies tat oe sould use te Lax metod if oe ca be satisfied wit a time step suc tat te courat umber is close to a uit. For stability, Crak-Nicolso was te most stable of all metods. Stability by itself is ot sufficiet coditio to use to select a umerical sceme. It must also be accurate. Te C-R metod as bot tese properties for te rage of te time steps cosidered. But as metioed above, tere is a rage of time steps i wic te Lax metod is more accurate ta all te oter metods. For CPU usage, te explicit metods used less CPU time we te time step was small, up to sec. Tis ca be explaied as follows: for small step size, te umber of time to solve Ax = b is large. Hece te implicit metods will be slower. As te time step is icreased to te rage of sec ad over, te implicit metods actually became more CPU efficiet due to te fact tat te umber of times to solve Ax = b is less because te umber of steps is less. I coclusio, te selectio of a fiite differece sceme depeds o may factors. Stability ad accuracy beig te most importat. Te time step size plays a critical rule. For Courat umber close to a uity, te Lax metod is te most attractive. For larger time steps, te C-R metod sould be cosidered. 24
25 6 Appedix 6.1 Plots case 1 25
26 26
27 6.1.2 case 2 27
28 28
29 6.1.3 case 3 29
30 30
31 6.1.4 case 4 31
32 32
33 6.1.5 case 5 33
34 34
35 6.1.6 case 6 35
36 36
37 6.1.7 case 7 37
38 38
39 6.1.8 case 8 39
40 40
41 6.2 Source code 1!******************************************* 2!* 3!* Solve te advectio PDE usig Explicit FTCS, 4!* Explicit Lax, Implicit FTCS, ad implicit Crak-Nicolso 5!* metods for costat ad varyig speed. 6!* 7!* Solve dc/dt = -u dc/dx 8!* u = t/20 ft/miute 9!* ad 10!* u costat 11!* 12!* Compiler used: gu 95 (g95) o Cygwi. Gcc !* Date: Jue !* 15!* by Nasser Abbasi 16!******************************************* PROGRAM advectio 19 IMPLICIT NONE REAL :: DT,DX,max_ru_time,legt,sapsot_, & 22 first_limit,secod_limit 23 INTEGER :: N,SNAPSHOTS 24 caracter(10) :: cmd_arg! to read time step from commad lie INTEGER :: metod! 1=FTCS, 2=LAX, 3=Implicit FTCS, 4=C-R 27 INTEGER :: mode! 1=Fixed wid speed, 2=speed fuctio of time REAL :: t_start, t_ed, cpu_time_used,ed_lie(1002) 30 INTEGER :: ALL_DATA_FILE_ID 31 PARAMETER(ALL_DATA_FILE_ID=900) 32 33! Iitialize data. All metods will use te same 34! parameters to make comparig tem easier 35 36! read t from commad lie. 37 CALL getarg(1,cmd_arg) 38 cmd_arg=trim(cmd_arg) 39 prit *,'= ', cmd_arg 40 read(cmd_arg,*)dt! i time, i miutes prit *,'Dt=',DT N = 1000! umber of grid poits i space 45 legt = 100! legt of space solutio i feet first_limit = 0.25*legt 48 secod_limit = 0.35*legt DX = legt/n! i space, i feets max_ru_time = 30.0! ow log to ru for i miutes 53 SNAPSHOTS = 200! umber of sapsots per ru. Used for aimatio sapsot_ = max_ru_time / SNAPSHOTS! time betwee eac sap sot prit *,'DT=',DT,' miutes, DX=',DX,' feets' 58 prit *,'takig sapsots every ', sapsot_,' miutes' DO mode=1,2 61 prit*,'=======> processig mode ',mode 62 DO metod=1,4! No eumeratio data types i Fotra CALL CPU_TIME(t_start)! get curret CPU time 65 CALL process(mode,metod,n,dt,dx,max_ru_time,sapsot_,& 66 first_limit,secod_limit) 67 CALL CPU_TIME(t_ed)! get curret CPU time cpu_time_used = t_ed - t_start WRITE(*,FMT='(A,I2,A,F12.5)') 'CPU TIME used for metod', metod, ' = ', cpu_time_used 72! Now record test case parameters i last lie 41
42 73 ed_lie=0 74 ed_lie(1)=cpu_time_used 75 ed_lie(2)=dt 76 ed_lie(3)=dx 77 ed_lie(4)=mode 78 ed_lie(5)=metod WRITE(UNIT=ALL_DATA_FILE_ID,FMT=*) ed_lie 81 CLOSE(ALL_DATA_FILE_ID) END DO 84 END DO END PROGRAM advectio 87!************************************ 88!* 89!* 90!************************************ 91 SUBROUTINE process(mode,metod,n,dt,dx,max_ru_time,sapsot_,& 92 first_limit,secod_limit) 93 IMPLICIT NONE INTEGER, INTENT(IN) :: mode,metod,n 96 REAL, INTENT(IN) :: DT,DX,max_ru_time,sapsot_,& 97 first_limit,secod_limit INTEGER :: I 100 LOGICAL :: sap_sot_at_15_take 101 INTEGER :: ALL_DATA_FILE_ID 102 PARAMETER(ALL_DATA_FILE_ID=900) 103 REAL :: sap_curret_time 104 REAL :: curret_time 105 REAL :: C(N)! curret solutio 106 REAL :: CNEW(N)! future solutio 107 REAL :: CEXACT(N)! curret exact solutio 108 REAL :: curret_first_limit 109 REAL :: A(N,N),aa(N),b(2:N),cc(N-1),CTEMP(N)! for C-R ad implicit FTCS 110 REAL :: K,speed 111 REAL :: error,rms! root mea square error betwee curret ad iitial sol curret_time = sap_curret_time = CALL iitialize_solutio(c,n,dx,first_limit,secod_limit) 117 CEXACT = C 118 curret_first_limit = first_limit CALL pre_loop_iitializatio(mode,metod,curret_time,k, & DT,DX,N,C,ALL_DATA_FILE_ID, & 122 A,aa,b,cc ) sap_sot_at_15_take=.false DO WHILE(curret_time < max_ru_time) IF( sap_curret_time sapsot_ ) THEN 129 sap_curret_time = WRITE( UNIT=ALL_DATA_FILE_ID, FMT=*) curret_time, error, C 131 END IF SELECT CASE(metod) CASE( 1:2 ) IF(metod==1) THEN! ftcs 138 IF(mode==2)THEN 139 K = speed(mode,curret_time)*dt/(2.*dx) 140 ENDIF DO I = 2,N CNEW(I) = C(I) - K * ( C(I+1) - C(I-1) ) 144 END DO 145 ELSE!lax 146 IF(mode == 2) THEN 147 K = speed(mode,curret_time)*dt/(dx) 148 ENDIF 42
43 DO I = 2,N CNEW(I) = C(I) - K/2. * ( C(I+1) - C(I-1) ) + & 152 (K**2.)/2 * ( C(I+1) +C(I-1)-2.*C(I) ) 153 END DO 154 END IF CNEW(1) = C(1) 157 CNEW(N) = C(N)! Boudary coditios 158 C=CNEW CASE( 3 )! implicit ftcs IF( mode == 2) THEN! oly eed to update Matrix for varyig U 163 K = speed(mode,curret_time)*dt/(2.*dx) CALL iit_a_matrix(a,k,n) CALL iit_diagoal_vectors(n,a,cc,aa,b) 167 END IF CALL solve_tomas_algoritm(n,aa,b,cc,c,cnew) 170 C = CNEW CASE( 4 )! C-R IF(mode == 2) THEN!oly eed to update A if U cages 175 K = speed(mode,curret_time)*dt/(4*dx)! C-R 176 CALL iit_a_matrix(a,k,n) 177 CALL iit_diagoal_vectors(n,a,cc,aa,b) 178 END IF CTEMP(1) = C(1) 181 CTEMP(N) = C(N) DO I=2,N CTEMP(I)=C(I)+K*C(I-1)-K*C(I+1) 185 END DO CALL solve_tomas_algoritm(n,aa,b,cc,ctemp,c) END SELECT IF( curret_time 15.0.AND. (.NOT. sap_sot_at_15_take)) THEN 192 sap_sot_at_15_take =.TRUE. 193 CALL take_oe_sap_sot(mode,metod,15,n,c,dx) 194 END IF curret_time = curret_time + DT 197 curret_first_limit = curret_first_limit + speed(mode,curret_time)*dt 198 CALL get_curret_exact_solutio(cexact,n,curret_first_limit,dx) 199 error = RMS(CEXACT,C,N) sap_curret_time = sap_curret_time + DT END DO CALL take_oe_sap_sot(mode,metod,30,n,c,dx) END SUBROUTINE process 208!************************************ 209!* 210!* 211!************************************ 212 SUBROUTINE pre_loop_iitializatio(mode,metod,curret_time,k,& 213 DT,DX,N,C,ALL_DATA_FILE_ID,& 214 A,aa,b,cc) 215 IMPLICIT NONE INTEGER, INTENT(IN) 218 REAL, INTENT(IN) :: mode,metod,n,all_data_file_id :: C(N),DT,DX,curret_time 219 REAL, INTENT(OUT) :: K,A(N,N),aa(N),b(2:N),cc(N-1) 220 REAL :: speed SELECT CASE(metod) 223 CASE( 1 )! FTCS
44 225 K = speed(mode,curret_time)*dt/(2.*dx) IF(mode==1) THEN 228 OPEN(UNIT=ALL_DATA_FILE_ID, file='expall.txt')! all time sots 229 CALL prit_to_file(c,'exp0.txt',n,dx) 230 ELSE 231 OPEN(UNIT=ALL_DATA_FILE_ID, file='exp_extraall.txt')! all time sots 232 CALL prit_to_file(c,'exp_extra0.txt',n,dx) 233 END IF CASE( 2 )! Lax K = speed(mode,curret_time)*dt/(dx) IF(mode==1) THEN 240 OPEN(UNIT=ALL_DATA_FILE_ID, file='laxall.txt')! all time sots 241 CALL prit_to_file(c,'lax0.txt',n,dx) 242 ELSE 243 OPEN(UNIT=ALL_DATA_FILE_ID, file='lax_extraall.txt')! all time sots 244 CALL prit_to_file(c,'lax_extra0.txt',n,dx) 245 END IF CASE( 3 )! Implicit FTCS K = speed(mode,curret_time)*dt/(2.*dx) CALL iit_a_matrix(a,k,n) 252 CALL iit_diagoal_vectors(n,a,cc,aa,b) IF(mode==1) THEN 255 OPEN(UNIT=ALL_DATA_FILE_ID, file='impall.txt')! all time sots 256 CALL prit_to_file(c,'imp0.txt',n,dx) 257 ELSE 258 OPEN(UNIT=ALL_DATA_FILE_ID, file='imp_extraall.txt')! all time sots 259 CALL prit_to_file(c,'imp_extra0.txt',n,dx) 260 END IF CASE( 4 )! C-R K = speed(mode,curret_time)*dt/(4*dx)! C-R CALL iit_a_matrix(a,k,n) 267 CALL iit_diagoal_vectors(n,a,cc,aa,b) IF(mode==1) THEN 270 OPEN(UNIT=ALL_DATA_FILE_ID, file='crall.txt')! all time sots 271 CALL prit_to_file(c,'cr0.txt',n,dx) 272 ELSE 273 OPEN(UNIT=ALL_DATA_FILE_ID, file='cr_extraall.txt')! all time sots 274 CALL prit_to_file(c,'cr_extra0.txt',n,dx) 275 END IF 276 END SELECT WRITE( UNIT=ALL_DATA_FILE_ID, FMT=*) curret_time,0, C END SUBROUTINE pre_loop_iitializatio 281!************************************ 282!* 283!* 284!************************************ 285 SUBROUTINE iit_diagoal_vectors(n,a,cc,aa,b) 286 IMPLICIT NONE INTEGER, INTENT(IN) ::N 289 REAL, INTENT(IN) ::A(N,N) 290 REAL, INTENT(OUT) ::aa(n),b(2:n),cc(n-1) INTEGER ::I,J J=2 295 DO I=1,N cc(i)=a(i,j) 297 J=J END DO 299 cc(1)=
45 301 DO I=1,N 302 aa(i)=a(i,i) 303 END DO J=1 306 DO I=2,N 307 b(i)=a(i,j) 308 J=J END DO END SUBROUTINE iit_diagoal_vectors 312!************************************ 313!* 314!* 315!************************************ 316 SUBROUTINE iitialize_solutio(c,n,dx,first_limit,secod_limit) 317 IMPLICIT NONE INTEGER, INTENT(IN) :: N 320 REAL, INTENT(IN) :: DX,first_limit,secod_limit 321 REAL, INTENT(INOUT) :: C(0:N-1) INTEGER :: I 324 REAL :: x, PI,av,R PARAMETER( PI = ACOS(-1.) ) x = av = (secod_limit+first_limit)/ R = av - first_limit C = DO I=0,N IF( x first_limit.and. x secod_limit ) THEN 337 C(I) = 1 + COS( PI * (x-av)/r ) 338 END IF x = x + DX 341 END DO END SUBROUTINE iitialize_solutio 344!************************************ 345!* 346!* 347!************************************ 348 SUBROUTINE prit_to_file(c,file_ame,n,dx) 349 IMPLICIT NONE REAL, INTENT(IN) :: C(N),DX 353 INTEGER, INTENT(IN) :: N CHARACTER* (*), INTENT(IN) :: file_ame INTEGER :: I 358 INTEGER :: FILE_ID 359 PARAMETER(FILE_ID=999) 360 REAL :: curret_positio OPEN(UNIT=FILE_ID, file=file_ame) curret_positio = 0; 365 DO I=1,N WRITE( UNIT=FILE_ID, FMT=* ) curret_positio,'\t', C(I) 368 curret_positio = curret_positio + DX END DO CLOSE(FILE_ID) END SUBROUTINE prit_to_file !************************************!* 45
46 377!* 378!************************************ 379 SUBROUTINE iit_a_matrix(a,k,n) 380 IMPLICIT NONE INTEGER, INTENT(IN) ::N 383 REAL, INTENT(IN) ::K 384 REAL, INTENT(OUT) ::A(N,N) INTEGER ::I DO I = 2,N A(I,I-1) = -K 390 A(I,I) = A(I,I+1) = K 392 END DO A(1,1) = A(N,N) = END SUBROUTINE iit_a_matrix 398!************************************ 399!* 400!* 401!************************************ 402 SUBROUTINE solve_tomas_algoritm(n,aa,b,c,old_c,ew_c) 403 IMPLICIT NONE REAL, INTENT(IN) :: aa(n),b(2:n),c(n-1),old_c(n) 406 INTEGER, INTENT(IN) :: N 407 REAL, INTENT(INOUT) :: ew_c(n) INTEGER :: I 410 REAL :: alpa(n),beta(2:n),g(n) alpa(1) = aa(1) 413 DO I=2,N 414 beta(i)=b(i)/alpa(i-1) 415 alpa(i)=aa(i)-beta(i)*c(i-1) 416 END DO g(1)=old_c(1) 419 DO I=2,N 420 g(i)=old_c(i)-beta(i)*g(i-1) 421 END DO ew_c(n)=g(n)/alpa(n) 424 DO I=N-1,1, ew_c(i)=(g(i)-c(i)*ew_c(i+1))/alpa(i) 426 END DO END SUBROUTINE solve_tomas_algoritm 429!************************************ 430!* 431!* 432!************************************ 433 REAL FUNCTION speed(mode,time) 434 IMPLICIT NONE INTEGER, INTENT(IN) :: MODE 437 REAL, INTENT(IN) :: time IF( MODE == 1 ) THEN 440 speed= ELSE 442 speed=time/ END IF END FUNCTION speed 446!************************************ 447!* 448!* 449!************************************ 450 SUBROUTINE take_oe_sap_sot(mode,metod,time,n,c,dx) 451 IMPLICIT NONE
47 453 INTEGER, INTENT(IN) ::TIME,mode,metod,N 454 REAL, INTENT(IN) ::C(N),DX IF(TIME==15) THEN 457 SELECT CASE(metod) 458 CASE(1) 459 IF(mode==1) THEN 460 CALL prit_to_file(c,'exp15.txt',n,dx) 461 ELSE 462 CALL prit_to_file(c,'exp_extra15.txt',n,dx) 463 END IF 464 CASE(2) 465 IF(mode==1) THEN 466 CALL prit_to_file(c,'lax15.txt',n,dx) 467 ELSE 468 CALL prit_to_file(c,'lax_extra15.txt',n,dx) 469 ENDIF 470 CASE(3) 471 IF(mode==1) THEN 472 CALL prit_to_file(c,'imp15.txt',n,dx) 473 ELSE 474 CALL prit_to_file(c,'imp_extra15.txt',n,dx) 475 END IF 476 CASE(4) 477 IF(mode==1) THEN 478 CALL prit_to_file(c,'cr15.txt',n,dx) 479 ELSE 480 CALL prit_to_file(c,'cr_extra15.txt',n,dx) 481 END IF 482 END SELECT 483 ELSE 484 SELECT CASE(metod) 485 CASE(1) 486 IF(mode==1) THEN 487 CALL prit_to_file(c,'exp30.txt',n,dx) 488 ELSE 489 CALL prit_to_file(c,'exp_extra30.txt',n,dx) 490 END IF 491 CASE(2) 492 IF(mode==1) THEN 493 CALL prit_to_file(c,'lax30.txt',n,dx) 494 ELSE 495 CALL prit_to_file(c,'lax_extra30.txt',n,dx) 496 ENDIF 497 CASE(3) 498 IF(mode==1) THEN 499 CALL prit_to_file(c,'imp30.txt',n,dx) 500 ELSE 501 CALL prit_to_file(c,'imp_extra30.txt',n,dx) 502 END IF 503 CASE(4) 504 IF(mode==1) THEN 505 CALL prit_to_file(c,'cr30.txt',n,dx) 506 ELSE 507 CALL prit_to_file(c,'cr_extra30.txt',n,dx) 508 END IF 509 END SELECT 510 END IF 511 END SUBROUTINE take_oe_sap_sot 512!************************************ 513!* 514!* 515!************************************ 516 REAL FUNCTION RMS(CEXACT,C,N) 517 IMPLICIT NONE REAL, INTENT(IN) :: CEXACT(N),C(N) 520 INTEGER, INTENT(IN) :: N INTEGER :: I RMS= DO I=1,N 526 RMS = RMS+(CEXACT(I)-C(I))**2 527 END DO
48 529 RMS = RMS/N 530 RMS = SQRT(RMS) 531 END FUNCTION RMS 532!************************************ 533!* !************************************ 536 SUBROUTINE get_curret_exact_solutio(cexact,n,curret_first_limit,dx) 537 IMPLICIT NONE 538 REAL, INTENT(IN) :: curret_first_limit,dx 539 REAL, INTENT(OUT) :: CEXACT(0:N-1) 540 INTEGER, INTENT(IN) :: N INTEGER :: I 543 REAL :: first_limit 544 REAL :: secod_limit 545 REAL :: av,r,sift,x,pi PARAMETER( PI = ACOS(-1.) ) first_limit = secod_limit = sift = curret_first_limit - FIRST_LIMIT 553 first_limit = curret_first_limit 554 secod_limit = secod_limit + sift av = (secod_limit+first_limit)/ R = av - first_limit CEXACT = x = DO I = 0,N IF( x first_limit.and. x secod_limit ) THEN 564 CEXACT(I) = 1 + COS( PI * (x -av)/r ) 565 END IF x = x + DX 568 END DO 569 END SUBROUTINE get_curret_exact_solutio 7 Refereces 1. Numerical Metods for pysics. Secod editio. Alejadro Garcia 2. Applied Numerical Metods for Egieers. Terrece Akal. 3. Computatioal Teciques for fluid dyamics. Secod editio. C.A.J.Fletcer 48
Mixed and Implicit Schemes Implicit Schemes. Exercise: Verify that ρ is unimodular: ρ = 1.
Mixed ad Implicit Schemes 3..4 The leapfrog scheme is stable for the oscillatio equatio ad ustable for the frictio equatio. The Euler forward scheme is stable for the frictio equatio but ustable for the
More informationBuilding a Dynamic Two Dimensional Heat Transfer Model part #1
Buildig a Dyamic Two Dimesioal Heat Trasfer Model part #1 - Tis is te first alf of a tutorial wic sows ow to build a basic dyamic eat coductio model of a square plate. Te same priciple could be used to
More informationKernel Density Estimation. Let X be a random variable with continuous distribution F (x) and density f(x) = d
Kerel Desity Estimatio Let X be a radom variable wit cotiuous distributio F (x) ad desity f(x) = d dx F (x). Te goal is to estimate f(x). Wile F (x) ca be estimated by te EDF ˆF (x), we caot set ˆf(x)
More informationLimits of sequences. Contents 1. Introduction 2 2. Some notation for sequences The behaviour of infinite sequences 3
Limits of sequeces I this uit, we recall what is meat by a simple sequece, ad itroduce ifiite sequeces. We explai what it meas for two sequeces to be the same, ad what is meat by the -th term of a sequece.
More information1 Estimating sensitivities
Copyright c 27 by Karl Sigma 1 Estimatig sesitivities Whe estimatig the Greeks, such as the, the geeral problem ivolves a radom variable Y = Y (α) (such as a discouted payoff) that depeds o a parameter
More informationStatistics for Economics & Business
Statistics for Ecoomics & Busiess Cofidece Iterval Estimatio Learig Objectives I this chapter, you lear: To costruct ad iterpret cofidece iterval estimates for the mea ad the proportio How to determie
More information43. A 000 par value 5-year bod with 8.0% semiaual coupos was bought to yield 7.5% covertible semiaually. Determie the amout of premium amortized i the 6 th coupo paymet. (A).00 (B).08 (C).5 (D).5 (E).34
More informationHopscotch and Explicit difference method for solving Black-Scholes PDE
Mälardale iversity Fiacial Egieerig Program Aalytical Fiace Semiar Report Hopscotch ad Explicit differece method for solvig Blac-Scholes PDE Istructor: Ja Röma Team members: A Gog HaiLog Zhao Hog Cui 0
More informationCombining imperfect data, and an introduction to data assimilation Ross Bannister, NCEO, September 2010
Combiig imperfect data, ad a itroductio to data assimilatio Ross Baister, NCEO, September 00 rbaister@readigacuk The probability desity fuctio (PDF prob that x lies betwee x ad x + dx p (x restrictio o
More informationChapter 8. Confidence Interval Estimation. Copyright 2015, 2012, 2009 Pearson Education, Inc. Chapter 8, Slide 1
Chapter 8 Cofidece Iterval Estimatio Copyright 2015, 2012, 2009 Pearso Educatio, Ic. Chapter 8, Slide 1 Learig Objectives I this chapter, you lear: To costruct ad iterpret cofidece iterval estimates for
More informationOverlapping Generations
Eco. 53a all 996 C. Sims. troductio Overlappig Geeratios We wat to study how asset markets allow idividuals, motivated by the eed to provide icome for their retiremet years, to fiace capital accumulatio
More informationThe material in this chapter is motivated by Experiment 9.
Chapter 5 Optimal Auctios The material i this chapter is motivated by Experimet 9. We wish to aalyze the decisio of a seller who sets a reserve price whe auctioig off a item to a group of bidders. We begi
More informationCalculation of the Annual Equivalent Rate (AER)
Appedix to Code of Coduct for the Advertisig of Iterest Bearig Accouts. (31/1/0) Calculatio of the Aual Equivalet Rate (AER) a) The most geeral case of the calculatio is the rate of iterest which, if applied
More informationpoint estimator a random variable (like P or X) whose values are used to estimate a population parameter
Estimatio We have oted that the pollig problem which attempts to estimate the proportio p of Successes i some populatio ad the measuremet problem which attempts to estimate the mea value µ of some quatity
More informationSequences and Series
Sequeces ad Series Matt Rosezweig Cotets Sequeces ad Series. Sequeces.................................................. Series....................................................3 Rudi Chapter 3 Exercises........................................
More information11.7 (TAYLOR SERIES) NAME: SOLUTIONS 31 July 2018
.7 (TAYLOR SERIES NAME: SOLUTIONS 3 July 08 TAYLOR SERIES ( The power series T(x f ( (c (x c is called the Taylor Series for f(x cetered at x c. If c 0, this is called a Maclauri series. ( The N-th partial
More informationAn Empirical Study of the Behaviour of the Sample Kurtosis in Samples from Symmetric Stable Distributions
A Empirical Study of the Behaviour of the Sample Kurtosis i Samples from Symmetric Stable Distributios J. Marti va Zyl Departmet of Actuarial Sciece ad Mathematical Statistics, Uiversity of the Free State,
More informationWe learned: $100 cash today is preferred over $100 a year from now
Recap from Last Week Time Value of Moey We leared: $ cash today is preferred over $ a year from ow there is time value of moey i the form of willigess of baks, busiesses, ad people to pay iterest for its
More informationOutline. Plotting discrete-time signals. Sampling Process. Discrete-Time Signal Representations Important D-T Signals Digital Signals
Outlie Discrete-Time Sigals-Itroductio. Plottig discrete-time sigals. Samplig Process. Discrete-Time Sigal Represetatios Importat D-T Sigals Digital Sigals Discrete-Time Sigals-Itroductio The time variable
More informationCAPITAL PROJECT SCREENING AND SELECTION
CAPITAL PROJECT SCREEIG AD SELECTIO Before studyig the three measures of ivestmet attractiveess, we will review a simple method that is commoly used to scree capital ivestmets. Oe of the primary cocers
More informationUsing Math to Understand Our World Project 5 Building Up Savings And Debt
Usig Math to Uderstad Our World Project 5 Buildig Up Savigs Ad Debt Note: You will have to had i aswers to all umbered questios i the Project Descriptio See the What to Had I sheet for additioal materials
More informationWhen you click on Unit V in your course, you will see a TO DO LIST to assist you in starting your course.
UNIT V STUDY GUIDE Percet Notatio Course Learig Outcomes for Uit V Upo completio of this uit, studets should be able to: 1. Write three kids of otatio for a percet. 2. Covert betwee percet otatio ad decimal
More informationResearch Article The Probability That a Measurement Falls within a Range of n Standard Deviations from an Estimate of the Mean
Iteratioal Scholarly Research Network ISRN Applied Mathematics Volume 0, Article ID 70806, 8 pages doi:0.540/0/70806 Research Article The Probability That a Measuremet Falls withi a Rage of Stadard Deviatios
More informationAPPLICATION OF GEOMETRIC SEQUENCES AND SERIES: COMPOUND INTEREST AND ANNUITIES
APPLICATION OF GEOMETRIC SEQUENCES AND SERIES: COMPOUND INTEREST AND ANNUITIES Example: Brado s Problem Brado, who is ow sixtee, would like to be a poker champio some day. At the age of twety-oe, he would
More informationEXERCISE - BINOMIAL THEOREM
BINOMIAL THOEREM / EXERCISE - BINOMIAL THEOREM LEVEL I SUBJECTIVE QUESTIONS. Expad the followig expressios ad fid the umber of term i the expasio of the expressios. (a) (x + y) 99 (b) ( + a) 9 + ( a) 9
More informationModels of Asset Pricing
APPENDIX 1 TO CHAPTER 4 Models of Asset Pricig I this appedix, we first examie why diversificatio, the holdig of may risky assets i a portfolio, reduces the overall risk a ivestor faces. The we will see
More informationModels of Asset Pricing
APPENDIX 1 TO CHAPTER4 Models of Asset Pricig I this appedix, we first examie why diversificatio, the holdig of may risky assets i a portfolio, reduces the overall risk a ivestor faces. The we will see
More information5. Best Unbiased Estimators
Best Ubiased Estimators http://www.math.uah.edu/stat/poit/ubiased.xhtml 1 of 7 7/16/2009 6:13 AM Virtual Laboratories > 7. Poit Estimatio > 1 2 3 4 5 6 5. Best Ubiased Estimators Basic Theory Cosider agai
More informationMinhyun Yoo, Darae Jeong, Seungsuk Seo, and Junseok Kim
Hoam Mathematical J. 37 (15), No. 4, pp. 441 455 http://dx.doi.org/1.5831/hmj.15.37.4.441 A COMPARISON STUDY OF EXPLICIT AND IMPLICIT NUMERICAL METHODS FOR THE EQUITY-LINKED SECURITIES Mihyu Yoo, Darae
More informationTwitter: @Owe134866 www.mathsfreeresourcelibrary.com Prior Kowledge Check 1) State whether each variable is qualitative or quatitative: a) Car colour Qualitative b) Miles travelled by a cyclist c) Favourite
More informationFOUNDATION ACTED COURSE (FAC)
FOUNDATION ACTED COURSE (FAC) What is the Foudatio ActEd Course (FAC)? FAC is desiged to help studets improve their mathematical skills i preparatio for the Core Techical subjects. It is a referece documet
More informationModels of Asset Pricing
4 Appedix 1 to Chapter Models of Asset Pricig I this appedix, we first examie why diversificatio, the holdig of may risky assets i a portfolio, reduces the overall risk a ivestor faces. The we will see
More informationThe Time Value of Money in Financial Management
The Time Value of Moey i Fiacial Maagemet Muteau Irea Ovidius Uiversity of Costata irea.muteau@yahoo.com Bacula Mariaa Traia Theoretical High School, Costata baculamariaa@yahoo.com Abstract The Time Value
More informationAUTOMATIC GENERATION OF FUZZY PAYOFF MATRIX IN GAME THEORY
AUTOMATIC GENERATION OF FUZZY PAYOFF MATRIX IN GAME THEORY Dr. Farha I. D. Al Ai * ad Dr. Muhaed Alfarras ** * College of Egieerig ** College of Coputer Egieerig ad scieces Gulf Uiversity * Dr.farha@gulfuiversity.et;
More informationSubject CT5 Contingencies Core Technical. Syllabus. for the 2011 Examinations. The Faculty of Actuaries and Institute of Actuaries.
Subject CT5 Cotigecies Core Techical Syllabus for the 2011 Examiatios 1 Jue 2010 The Faculty of Actuaries ad Istitute of Actuaries Aim The aim of the Cotigecies subject is to provide a groudig i the mathematical
More informationOnline appendices from The xva Challenge by Jon Gregory. APPENDIX 10A: Exposure and swaption analogy.
APPENDIX 10A: Exposure ad swaptio aalogy. Sorese ad Bollier (1994), effectively calculate the CVA of a swap positio ad show this ca be writte as: CVA swap = LGD V swaptio (t; t i, T) PD(t i 1, t i ). i=1
More informationA random variable is a variable whose value is a numerical outcome of a random phenomenon.
The Practice of Statistics, d ed ates, Moore, ad Stares Itroductio We are ofte more iterested i the umber of times a give outcome ca occur tha i the possible outcomes themselves For example, if we toss
More informationAppendix 1 to Chapter 5
Appedix 1 to Chapter 5 Models of Asset Pricig I Chapter 4, we saw that the retur o a asset (such as a bod) measures how much we gai from holdig that asset. Whe we make a decisio to buy a asset, we are
More informationEVEN NUMBERED EXERCISES IN CHAPTER 4
Joh Riley 7 July EVEN NUMBERED EXERCISES IN CHAPTER 4 SECTION 4 Exercise 4-: Cost Fuctio of a Cobb-Douglas firm What is the cost fuctio of a firm with a Cobb-Douglas productio fuctio? Rather tha miimie
More informationEstimating Proportions with Confidence
Aoucemets: Discussio today is review for midterm, o credit. You may atted more tha oe discussio sectio. Brig sheets of otes ad calculator to midterm. We will provide Scatro form. Homework: (Due Wed Chapter
More informationRafa l Kulik and Marc Raimondo. University of Ottawa and University of Sydney. Supplementary material
Statistica Siica 009: Supplemet 1 L p -WAVELET REGRESSION WITH CORRELATED ERRORS AND INVERSE PROBLEMS Rafa l Kulik ad Marc Raimodo Uiversity of Ottawa ad Uiversity of Sydey Supplemetary material This ote
More informationMaximum Empirical Likelihood Estimation (MELE)
Maximum Empirical Likelihood Estimatio (MELE Natha Smooha Abstract Estimatio of Stadard Liear Model - Maximum Empirical Likelihood Estimator: Combiatio of the idea of imum likelihood method of momets,
More informationof Asset Pricing R e = expected return
Appedix 1 to Chapter 5 Models of Asset Pricig EXPECTED RETURN I Chapter 4, we saw that the retur o a asset (such as a bod) measures how much we gai from holdig that asset. Whe we make a decisio to buy
More informationChapter 8: Estimation of Mean & Proportion. Introduction
Chapter 8: Estimatio of Mea & Proportio 8.1 Estimatio, Poit Estimate, ad Iterval Estimate 8.2 Estimatio of a Populatio Mea: σ Kow 8.3 Estimatio of a Populatio Mea: σ Not Kow 8.4 Estimatio of a Populatio
More informationInstitute of Actuaries of India Subject CT5 General Insurance, Life and Health Contingencies
Istitute of Actuaries of Idia Subject CT5 Geeral Isurace, Life ad Health Cotigecies For 2017 Examiatios Aim The aim of the Cotigecies subject is to provide a groudig i the mathematical techiques which
More informationTopic-7. Large Sample Estimation
Topic-7 Large Sample Estimatio TYPES OF INFERENCE Ò Estimatio: É Estimatig or predictig the value of the parameter É What is (are) the most likely values of m or p? Ò Hypothesis Testig: É Decidig about
More informationOnline appendices from Counterparty Risk and Credit Value Adjustment a continuing challenge for global financial markets by Jon Gregory
Olie appedices from Couterparty Risk ad Credit Value Adjustmet a APPENDIX 8A: Formulas for EE, PFE ad EPE for a ormal distributio Cosider a ormal distributio with mea (expected future value) ad stadard
More informationii. Interval estimation:
1 Types of estimatio: i. Poit estimatio: Example (1) Cosider the sample observatios 17,3,5,1,18,6,16,10 X 8 X i i1 8 17 3 5 118 6 16 10 8 116 8 14.5 14.5 is a poit estimate for usig the estimator X ad
More informationToday: Finish Chapter 9 (Sections 9.6 to 9.8 and 9.9 Lesson 3)
Today: Fiish Chapter 9 (Sectios 9.6 to 9.8 ad 9.9 Lesso 3) ANNOUNCEMENTS: Quiz #7 begis after class today, eds Moday at 3pm. Quiz #8 will begi ext Friday ad ed at 10am Moday (day of fial). There will be
More informationx satisfying all regularity conditions. Then
AMS570.01 Practice Midterm Exam Sprig, 018 Name: ID: Sigature: Istructio: This is a close book exam. You are allowed oe-page 8x11 formula sheet (-sided). No cellphoe or calculator or computer is allowed.
More informationCHAPTER 2 PRICING OF BONDS
CHAPTER 2 PRICING OF BONDS CHAPTER SUARY This chapter will focus o the time value of moey ad how to calculate the price of a bod. Whe pricig a bod it is ecessary to estimate the expected cash flows ad
More informationMonetary Economics: Problem Set #5 Solutions
Moetary Ecoomics oblem Set #5 Moetary Ecoomics: oblem Set #5 Solutios This problem set is marked out of 1 poits. The weight give to each part is idicated below. Please cotact me asap if you have ay questios.
More informationSampling Distributions and Estimation
Cotets 40 Samplig Distributios ad Estimatio 40.1 Samplig Distributios 40. Iterval Estimatio for the Variace 13 Learig outcomes You will lear about the distributios which are created whe a populatio is
More informationof Asset Pricing APPENDIX 1 TO CHAPTER EXPECTED RETURN APPLICATION Expected Return
APPENDIX 1 TO CHAPTER 5 Models of Asset Pricig I Chapter 4, we saw that the retur o a asset (such as a bod) measures how much we gai from holdig that asset. Whe we make a decisio to buy a asset, we are
More informationThe Limit of a Sequence (Brief Summary) 1
The Limit of a Sequece (Brief Summary). Defiitio. A real umber L is a it of a sequece of real umbers if every ope iterval cotaiig L cotais all but a fiite umber of terms of the sequece. 2. Claim. A sequece
More informationSTRAND: FINANCE. Unit 3 Loans and Mortgages TEXT. Contents. Section. 3.1 Annual Percentage Rate (APR) 3.2 APR for Repayment of Loans
CMM Subject Support Strad: FINANCE Uit 3 Loas ad Mortgages: Text m e p STRAND: FINANCE Uit 3 Loas ad Mortgages TEXT Cotets Sectio 3.1 Aual Percetage Rate (APR) 3.2 APR for Repaymet of Loas 3.3 Credit Purchases
More informationFaculdade de Economia da Universidade de Coimbra
Faculdade de Ecoomia da Uiversidade de Coimbra Grupo de Estudos Moetários e Fiaceiros (GEMF) Av. Dias da Silva, 65 300-5 COIMBRA, PORTUGAL gemf@fe.uc.pt http://www.uc.pt/feuc/gemf PEDRO GODINHO Estimatig
More informationMA Lesson 11 Section 1.3. Solving Applied Problems with Linear Equations of one Variable
MA 15200 Lesso 11 Sectio 1. I Solvig Applied Problems with Liear Equatios of oe Variable 1. After readig the problem, let a variable represet the ukow (or oe of the ukows). Represet ay other ukow usig
More informationMATH : EXAM 2 REVIEW. A = P 1 + AP R ) ny
MATH 1030-008: EXAM 2 REVIEW Origially, I was havig you all memorize the basic compoud iterest formula. I ow wat you to memorize the geeral compoud iterest formula. This formula, whe = 1, is the same as
More informationINTERVAL GAMES. and player 2 selects 1, then player 2 would give player 1 a payoff of, 1) = 0.
INTERVAL GAMES ANTHONY MENDES Let I ad I 2 be itervals of real umbers. A iterval game is played i this way: player secretly selects x I ad player 2 secretly ad idepedetly selects y I 2. After x ad y are
More informationChapter 5: Sequences and Series
Chapter 5: Sequeces ad Series 1. Sequeces 2. Arithmetic ad Geometric Sequeces 3. Summatio Notatio 4. Arithmetic Series 5. Geometric Series 6. Mortgage Paymets LESSON 1 SEQUENCES I Commo Core Algebra I,
More informationIntroduction to Probability and Statistics Chapter 7
Itroductio to Probability ad Statistics Chapter 7 Ammar M. Sarha, asarha@mathstat.dal.ca Departmet of Mathematics ad Statistics, Dalhousie Uiversity Fall Semester 008 Chapter 7 Statistical Itervals Based
More informationInferential Statistics and Probability a Holistic Approach. Inference Process. Inference Process. Chapter 8 Slides. Maurice Geraghty,
Iferetial Statistics ad Probability a Holistic Approach Chapter 8 Poit Estimatio ad Cofidece Itervals This Course Material by Maurice Geraghty is licesed uder a Creative Commos Attributio-ShareAlike 4.0
More informationProblem Set 1a - Oligopoly
Advaced Idustrial Ecoomics Sprig 2014 Joha Steek 6 may 2014 Problem Set 1a - Oligopoly 1 Table of Cotets 2 Price Competitio... 3 2.1 Courot Oligopoly with Homogeous Goods ad Differet Costs... 3 2.2 Bertrad
More informationCAPITAL ASSET PRICING MODEL
CAPITAL ASSET PRICING MODEL RETURN. Retur i respect of a observatio is give by the followig formula R = (P P 0 ) + D P 0 Where R = Retur from the ivestmet durig this period P 0 = Curret market price P
More informationJournal of Statistical Software
JSS Joural of Statistical Software Jue 2007, Volume 19, Issue 6. http://www.jstatsoft.org/ Ratioal Arithmetic Mathematica Fuctios to Evaluate the Oe-sided Oe-sample K-S Cumulative Samplig Distributio J.
More information0.1 Valuation Formula:
0. Valuatio Formula: 0.. Case of Geeral Trees: q = er S S S 3 S q = er S S 4 S 5 S 4 q 3 = er S 3 S 6 S 7 S 6 Therefore, f (3) = e r [q 3 f (7) + ( q 3 ) f (6)] f () = e r [q f (5) + ( q ) f (4)] = f ()
More informationParametric Density Estimation: Maximum Likelihood Estimation
Parametric Desity stimatio: Maimum Likelihood stimatio C6 Today Itroductio to desity estimatio Maimum Likelihood stimatio Itroducto Bayesia Decisio Theory i previous lectures tells us how to desig a optimal
More information0.07. i PV Qa Q Q i n. Chapter 3, Section 2
Chapter 3, Sectio 2 1. (S13HW) Calculate the preset value for a auity that pays 500 at the ed of each year for 20 years. You are give that the aual iterest rate is 7%. 20 1 v 1 1.07 PV Qa Q 500 5297.01
More information1 + r. k=1. (1 + r) k = A r 1
Perpetual auity pays a fixed sum periodically forever. Suppose a amout A is paid at the ed of each period, ad suppose the per-period iterest rate is r. The the preset value of the perpetual auity is A
More informationDESCRIPTION OF MATHEMATICAL MODELS USED IN RATING ACTIVITIES
July 2014, Frakfurt am Mai. DESCRIPTION OF MATHEMATICAL MODELS USED IN RATING ACTIVITIES This documet outlies priciples ad key assumptios uderlyig the ratig models ad methodologies of Ratig-Agetur Expert
More informationA New Approach to Obtain an Optimal Solution for the Assignment Problem
Iteratioal Joural of Sciece ad Research (IJSR) ISSN (Olie): 231-7064 Idex Copericus Value (2013): 6.14 Impact Factor (2015): 6.31 A New Approach to Obtai a Optimal Solutio for the Assigmet Problem A. Seethalakshmy
More informationFurther Pure 1 Revision Topic 5: Sums of Series
The OCR syllabus says that cadidates should: Further Pure Revisio Topic 5: Sums of Series Cadidates should be able to: (a) use the stadard results for Σr, Σr, Σr to fid related sums; (b) use the method
More informationA point estimate is the value of a statistic that estimates the value of a parameter.
Chapter 9 Estimatig the Value of a Parameter Chapter 9.1 Estimatig a Populatio Proportio Objective A : Poit Estimate A poit estimate is the value of a statistic that estimates the value of a parameter.
More informationBinomial Model. Stock Price Dynamics. The Key Idea Riskless Hedge
Biomial Model Stock Price Dyamics The value of a optio at maturity depeds o the price of the uderlyig stock at maturity. The value of the optio today depeds o the expected value of the optio at maturity
More information4.5 Generalized likelihood ratio test
4.5 Geeralized likelihood ratio test A assumptio that is used i the Athlete Biological Passport is that haemoglobi varies equally i all athletes. We wish to test this assumptio o a sample of k athletes.
More informationCreditRisk + Download document from CSFB web site:
CreditRis + Dowload documet from CSFB web site: http://www.csfb.com/creditris/ Features of CreditRis+ pplies a actuarial sciece framewor to the derivatio of the loss distributio of a bod/loa portfolio.
More information1 ECON4415: International Economics Problem Set 4 - Solutions
ECON445: Iteratioal Ecoomics Problem Set 4 - Solutios. I Moopolistic competitio. Moopolistic competitio is a market form where May rms producig di eret varieties. Each rm has moopoly power over its ow
More informationOptimizing of the Investment Structure of the Telecommunication Sector Company
Iteratioal Joural of Ecoomics ad Busiess Admiistratio Vol. 1, No. 2, 2015, pp. 59-70 http://www.aisciece.org/joural/ijeba Optimizig of the Ivestmet Structure of the Telecommuicatio Sector Compay P. N.
More informationEconomic Computation and Economic Cybernetics Studies and Research, Issue 2/2016, Vol. 50
Ecoomic Computatio ad Ecoomic Cyberetics Studies ad Research, Issue 2/216, Vol. 5 Kyoug-Sook Moo Departmet of Mathematical Fiace Gacho Uiversity, Gyeoggi-Do, Korea Yuu Jeog Departmet of Mathematics Korea
More informationMath 312, Intro. to Real Analysis: Homework #4 Solutions
Math 3, Itro. to Real Aalysis: Homework #4 Solutios Stephe G. Simpso Moday, March, 009 The assigmet cosists of Exercises 0.6, 0.8, 0.0,.,.3,.6,.0,.,. i the Ross textbook. Each problem couts 0 poits. 0.6.
More informationsetting up the business in sage
3 settig up the busiess i sage Chapter itroductio Settig up a computer accoutig program for a busiess or other orgaisatio will take some time, but as log as the correct data is etered i the correct format
More informationAnomaly Correction by Optimal Trading Frequency
Aomaly Correctio by Optimal Tradig Frequecy Yiqiao Yi Columbia Uiversity September 9, 206 Abstract Uder the assumptio that security prices follow radom walk, we look at price versus differet movig averages.
More informationChapter 4 - Consumer. Household Demand and Supply. Solving the max-utility problem. Working out consumer responses. The response function
Almost essetial Cosumer: Optimisatio Chapter 4 - Cosumer Osa 2: Household ad supply Cosumer: Welfare Useful, but optioal Firm: Optimisatio Household Demad ad Supply MICROECONOMICS Priciples ad Aalysis
More informationMulti-Criteria Flow-Shop Scheduling Optimization
Multi-Criteria Flow-Shop Schedulig Optimizatio A Seior Project Submitted I Partial Fulfillmet Of the Requiremets for the Degree of Bachelor of Sciece i Idustrial Egieerig Preseted to: The Faculty of Califoria
More informationr i = a i + b i f b i = Cov[r i, f] The only parameters to be estimated for this model are a i 's, b i 's, σe 2 i
The iformatio required by the mea-variace approach is substatial whe the umber of assets is large; there are mea values, variaces, ad )/2 covariaces - a total of 2 + )/2 parameters. Sigle-factor model:
More informationA New Constructive Proof of Graham's Theorem and More New Classes of Functionally Complete Functions
A New Costructive Proof of Graham's Theorem ad More New Classes of Fuctioally Complete Fuctios Azhou Yag, Ph.D. Zhu-qi Lu, Ph.D. Abstract A -valued two-variable truth fuctio is called fuctioally complete,
More informationACTUARIAL RESEARCH CLEARING HOUSE 1990 VOL. 2 INTEREST, AMORTIZATION AND SIMPLICITY. by Thomas M. Zavist, A.S.A.
ACTUARIAL RESEARCH CLEARING HOUSE 1990 VOL. INTEREST, AMORTIZATION AND SIMPLICITY by Thomas M. Zavist, A.S.A. 37 Iterest m Amortizatio ad Simplicity Cosider simple iterest for a momet. Suppose you have
More informationLecture 9: The law of large numbers and central limit theorem
Lecture 9: The law of large umbers ad cetral limit theorem Theorem.4 Let X,X 2,... be idepedet radom variables with fiite expectatios. (i) (The SLLN). If there is a costat p [,2] such that E X i p i i=
More information18.S096 Problem Set 5 Fall 2013 Volatility Modeling Due Date: 10/29/2013
18.S096 Problem Set 5 Fall 2013 Volatility Modelig Due Date: 10/29/2013 1. Sample Estimators of Diffusio Process Volatility ad Drift Let {X t } be the price of a fiacial security that follows a geometric
More informationPositivity Preserving Schemes for Black-Scholes Equation
Research Joural of Fiace ad Accoutig IN -97 (Paper) IN -7 (Olie) Vol., No.7, 5 Positivity Preservig chemes for Black-choles Equatio Mohammad Mehdizadeh Khalsaraei (Correspodig author) Faculty of Mathematical
More informationSubject CT1 Financial Mathematics Core Technical Syllabus
Subject CT1 Fiacial Mathematics Core Techical Syllabus for the 2018 exams 1 Jue 2017 Subject CT1 Fiacial Mathematics Core Techical Aim The aim of the Fiacial Mathematics subject is to provide a groudig
More informationEstimating Forward Looking Distribution with the Ross Recovery Theorem
roceedigs of the Asia acific Idustrial Egieerig & Maagemet Systems Coferece 5 Estimatig Forward Lookig Distributio with the Ross Recovery Theorem Takuya Kiriu Graduate School of Sciece ad Techology Keio
More informationSolutions to Problem Sheet 1
Solutios to Problem Sheet ) Use Theorem.4 to prove that p log for all real x 3. This is a versio of Theorem.4 with the iteger N replaced by the real x. Hit Give x 3 let N = [x], the largest iteger x. The,
More informationThese characteristics are expressed in terms of statistical properties which are estimated from the sample data.
0. Key Statistical Measures of Data Four pricipal features which characterize a set of observatios o a radom variable are: (i) the cetral tedecy or the value aroud which all other values are buched, (ii)
More informationA Hybrid Finite Difference Method for Valuing American Puts
Proceedigs of the World Cogress o Egieerig 29 Vol II A Hybrid Fiite Differece Method for Valuig America Puts Ji Zhag SogPig Zhu Abstract This paper presets a umerical scheme that avoids iteratios to solve
More informationad covexity Defie Macaulay duratio D Mod = r 1 = ( CF i i k (1 + r k) i ) (1.) (1 + r k) C = ( r ) = 1 ( CF i i(i + 1) (1 + r k) i+ k ) ( ( i k ) CF i
Fixed Icome Basics Cotets Duratio ad Covexity Bod Duratios ar Rate, Spot Rate, ad Forward Rate Flat Forward Iterpolatio Forward rice/yield, Carry, Roll-Dow Example Duratio ad Covexity For a series of cash
More information5 Statistical Inference
5 Statistical Iferece 5.1 Trasitio from Probability Theory to Statistical Iferece 1. We have ow more or less fiished the probability sectio of the course - we ow tur attetio to statistical iferece. I statistical
More informationThis article is part of a series providing
feature Bryce Millard ad Adrew Machi Characteristics of public sector workers SUMMARY This article presets aalysis of public sector employmet, ad makes comparisos with the private sector, usig data from
More informationForecasting bad debt losses using clustering algorithms and Markov chains
Forecastig bad debt losses usig clusterig algorithms ad Markov chais Robert J. Till Experia Ltd Lambert House Talbot Street Nottigham NG1 5HF {Robert.Till@uk.experia.com} Abstract Beig able to make accurate
More information