Evaluation of maintenance policies for equipment subject to quality shifts and failures

Size: px
Start display at page:

Download "Evaluation of maintenance policies for equipment subject to quality shifts and failures"

Transcription

1 Evaluaion of mainenance policies for equipmen subjec o qualiy shifs and failures Sofia Panagioidou, George Tagaras To cie his version: Sofia Panagioidou, George Tagaras. Evaluaion of mainenance policies for equipmen subjec o qualiy shifs and failures. Inernaional Journal of Producion Research, Taylor Francis, 00, (0), pp.-. <0.00/000000>. <hal-00> HAL Id: hal-00 hps://hal.archives-ouveres.fr/hal-00 Submied on Sep 00 HAL is a muli-disciplinary open access archive for he deposi and disseminaion of scienific research documens, wheher hey are published or no. The documens may come from eaching and research insiuions in France or abroad, or from public or privae research ceners. L archive ouvere pluridisciplinaire HAL, es desinée au dépô e à la diffusion de documens scienifiques de niveau recherche, publiés ou non, émanan des éablissemens d enseignemen e de recherche français ou érangers, des laboraoires publics ou privés.

2 Inernaional Journal of Producion Research Evaluaion of mainenance policies for equipmen subjec o qualiy shifs and failures Journal: Inernaional Journal of Producion Research Manuscrip ID: TPRS-00-IJPR-00.R Manuscrip Type: Original Manuscrip Dae Submied by he Auhor: -Aug-00 Complee Lis of Auhors: Panagioidou, Sofia; Arisole Universiy of Thessaloniki, Deparmen of Mechanical Engineering Tagaras, George; Arisole Universiy of Thessaloniki, Deparmen of Mechanical Engineering Keywords: STOCHASTIC MODELS, PREVENTIVE MAINTENANCE, MAINTENANCE PLANNING, PROCESS MONITORING Keywords (user):

3 Page of Inernaional Journal of Producion Research Evaluaion of Mainenance Policies for Equipmen Subjec o Qualiy Shifs and Failures SOFIA PANAGIOTIDOU (span@auh.gr) and GEORGE TAGARAS * (agaras@auh.gr) Arisole Universiy of Thessaloniki Deparmen of Mechanical Engineering Thessaloniki, Greece el.: fax.: Revised Augus 00 * Corresponding auhor

4 Inernaional Journal of Producion Research Page of Evaluaion of Mainenance Policies for Equipmen Subjec o Qualiy Shifs and Failures Absrac We develop an economic model for he opimizaion of mainenance procedures in a producion process wih wo qualiy saes. In addiion o deerioraing wih age, he equipmen may experience a jump o an ou-of-conrol sae (qualiy shif), which is characerized by lower producion revenues and higher endency o failure. The imes o qualiy shif and failure are allowed o be generally disribued random variables. We consider wo ypes of mainenance: minimal mainenance (MM) ha upgrades he qualiy sae of he equipmen wihou affecing is age and perfec prevenive mainenance (PM) ha fully upgrades he equipmen o he as-good-as-new condiion. We derive he expression for he expeced profi per ime uni and we invesigae, hrough a large number of numerical examples, he ype of he opimal soluion. I is concluded ha in pracically every case he opimal mainenance policy is an exreme one: i eiher calls for immediae MM as soon as a qualiy shif occurs (acive policy) or i allows operaion in he ou-of-conrol sae unil he ime of a scheduled PM acion (passive policy).

5 Page of Inernaional Journal of Producion Research Inroducion Equipmen mainenance is a very imporan operaion in almos every producion sysem. An appropriae prevenive mainenance (PM) policy no only reduces he probabiliy of equipmen failure bu also improves he working condiion of he equipmen resuling in lower producion coss and/or higher produc qualiy. The failure rae of he equipmen is ypically assumed o increase wih ime or usage (afer he iniial infan moraliy period) and consequenly age-based PM models have been widely sudied in he lieraure. These models usually assume ha he operaing condiion of he equipmen remains sable hroughou producion and no deerioraion mechanism exiss oher han complee failure. Someimes, hough, he equipmen may deeriorae o a less desirable working condiion before failing alogeher. This inferior condiion may be associaed wih boh higher operaing or qualiy coss and increased failure raes. For example, consider he cuing process of meal frames using power saws. From ime o ime, he cuing disc loses is balance and consequenly is abiliy o produce perfecly fla bur-free cus on he meal pars. In addiion o he lower qualiy of he produced meal pars, he misbalance leads o higher cenral axle faigue and higher probabiliy of breaking down (failure). Similarly, misbalance is a common malfuncion in many elecric and elecronic devices where a fan is used for he cooling process of he device. When his happens, he consequences are poor cooling, higher elecriciy consumpion, and higher proneness o failure of he elecric moor due o overheaing. In cases like he above condiion-monioring can provide useful informaion regarding he operaing sae of he equipmen and consequenly can lead o PM acions ha proec he equipmen more effecively agains failures. Several condiion-based PM models have been developed o address he problem of properly mainaining

6 Inernaional Journal of Producion Research Page of equipmen wih muliple operaing saes; see he surveys by Pierskalla and Voelker () and Valdez-Flores and Feldman (). However, he grea majoriy of hese models focus on he invesigaion of replacemen or perfec PM policies ha resore he equipmen o an as-good-as-new sae, while relaively few conribuions have been made o he field of imperfec PM. Imperfec PM is assumed o improve he working condiion of he equipmen bu no necessarily o an as-good-as-new sae. Various imperfec PM models have been suggesed in he lieraure, some of which are applicable o muli-sae producion processes/equipmen. See he survey by Pham and Wang () for more deails. Here we are mainly ineresed in cases where an imperfec PM upgrades he operaing sae of he equipmen by reducing is failure rae wihou affecing is age. A special case of his ype of imperfec PM is minimal mainenance (MM), which improves he equipmen condiion by one sae only. The objecive of his paper is o sudy a mainenance policy, including boh perfec and imperfec mainenance acions, which is appropriae for a producion process (equipmen) wih wo operaing saes, namely an in-conrol sae and an ou-ofconrol sae in Saisical Process Conrol erminology, and a failure sae. The wo disinc operaing saes are characerized by differen operaing and qualiy-relaed coss and by differen failure raes; he in-conrol sae has generally lower operaing/qualiy cos per ime uni and lower failure rae for he same equipmen age compared o he ou-of-conrol sae. Producion processes wih muliple operaing saes (qualiy saes) are ypically encounered in he conex of saisical process conrol, where a common pracice agains qualiy deerioraion is o bring he process back o is in-conrol sae afer he deecion of a qualiy shif o an ou-of-conrol sae. In his paper, we combine qualiy

7 Page of Inernaional Journal of Producion Research adjusmens o upgrade he process qualiy sae wih convenional mainenance acions o deal wih failures. More specifically, he proposed mainenance policy comprises he following ypes of mainenance acions: Perfec correcive mainenance (CM) upon failure: i resores he failed equipmen o he as-good-as-new sae. Perfec prevenive mainenance (PM) a some criical age: i resores he working equipmen o he as-good-as-new sae. Minimal mainenance (MM) applied only when he process is ou-ofconrol: i upgrades he equipmen from he ou-of-conrol o he in-conrol sae wihou affecing he equipmen age. The main conribuion of he presen paper is ha by allowing for generally disribued imes o qualiy shifs and o failures (wih non-decreasing failure raes) i grealy expands he model s realism and applicabiliy o acual producion sysems compared o exising models ha require a leas one of he above ime disribuions o be exponenial. An addiional conribuion lies in he resuls; i is documened ha in pracically all cases i suffices o consider and compare only exreme qualiy mainenance policies (acive or passive), which are easy o undersand and implemen. Secion presens a brief lieraure review while Secion describes he problem in deail and inroduces he necessary noaion. In secion we develop he mahemaical model while in secion we discuss he form of he opimal policy. Secion provides numerical examples and a sysemaic discussion of he effec of several model parameers on he opimal mainenance policy. The las secion summarizes he basic resuls and suggess direcions for fuure research.

8 Inernaional Journal of Producion Research Page of Lieraure Review The academic lieraure is replee wih models of mainenance policies comprising various ypes of acions (imperfec mainenance, minimal repair ec.). In his brief lieraure review we resric our aenion o models for producion processes (equipmen) wih muliple disinc operaing saes and a leas one failure sae, which assume perfec failure resoraion and include he noion of imperfec PM. More specifically we consider imperfec PM ha upgrades he operaing sae of he equipmen by reducing is failure rae wihou affecing is age. One of he mos imporan disinguishing feaures among hese models is he ype of he failure mechanism. We herefore sar he review wih models based on purely Markovian deerioraion mechanisms (exponenial disribuion of he ime o failure) and nex we proceed o he presenaion of more general non-markovian models. Purely Markovian deerioraion The pioneering work of Derman () concerns a repair-replacemen policy for a muli-sae Markovian deerioraing producion process. The model allows several alernaive mainenance decisions a every operaing sae of he process, which do no necessarily resore he equipmen o an as-good-as-new sae bu hey affec he sae ransiion probabiliies. In addiion, Derman () considers muliple inoperaive saes where replacemen is he only feasible acion. The objecive is o find he mainenance policy ha maximizes he expeced ime beween replacemens. I is shown ha his kind of problem can be expressed hrough a linear programming formulaion. Özekici and Günlük () sudy a similar problem bu wih deerminisic mainenance effecs (every mainenance acion leads o a cerain process sae) and a more general cos funcion. The objecive is o selec he mos appropriae mainenance

9 Page of Inernaional Journal of Producion Research acion for each process sae so as o minimize he expeced cos. Özekici and Günlük () provide srucural properies of he opimal mainenance policy under various cos srucures. A slighly differen approach o modeling of imperfec mainenance (repair) for producion processes wih muliple operaing saes and an absorbing failure sae has been proposed by Chiang and Yuan (00). The process is periodically inspeced and: a) nohing is done if he equipmen is found o be wihin a group of good saes, or b) he equipmen is repaired o a beer sae if found o be wihin a group of inermediae saes, or c) he equipmen is replaced if found o be wihin a group of inferior saes including he failure sae. The opimal policy is derived by minimizing he expeced long-run cos rae. A Markovian model which explicily combines qualiy conrol wih mainenance procedures has been developed by Tagaras (); i concerns a producion process subjec o qualiy shifs and failures, assuming muliple ou-of-conrol saes. The objecive is o opimize boh qualiy conrol schemes and PM procedures. In addiion o CM (resoraion) following a failure and periodic PM, qualiy adjusmens are carried ou whenever he process is found o operae in an ou-of-conrol sae. Non-Markovian deerioraion Mousafa e al. (00) sudy a muli-sae semi-markovian deerioraing sysem, allowing boh replacemen and minimal mainenance. They assume ha minimal mainenance upgrades he equipmen condiion by one sae. They use a conrol limi policy wih wo hreshold saes in a way similar o Chiang and Yuan (00), bu hey show ha his ype of policy is no always opimal.

10 Inernaional Journal of Producion Research Page of Makis and Fung (, ) sudy he inegraed problem of deermining he opimal qualiy conrol schedule and he opimal producion quaniy in a producion process wih wo qualiy saes and a single failure sae. The model allows for periodic prevenive replacemen as well. In boh models he ime o qualiy shif is assumed o be exponenially disribued while he ime o failure is assumed o be generally disribued. As soon as a qualiy shif is deeced he process is resored o he in-conrol sae wih he same equipmen age. I is worh noing ha he failure ime disribuion is assumed o be independen of he acual qualiy sae of he process and consequenly resoraions of qualiy shifs do no consiue PM acions agains failures. The producion process sudied in his paper is similar o ha of Makis and Fung (wo qualiy saes and a single failure sae) bu we focus on mainenance procedures raher han on he deerminaion of producion quaniies. In addiion, we consider general (no necessarily Markovian) deerioraion mechanisms no only for failure bu for qualiy shifs as well. Furhermore, we allow he failure ime disribuion o depend no only on he equipmen age bu also on is sae. Before concluding his brief lieraure review i should be noed, for he sae of compleeness, ha here exis a number of papers combining saisical process conrol (SPC) wih prevenive mainenance procedures. Alhough hese papers share some degree of similariy wih our work, hey approach he problem mosly from a qualiyoriened poin of view, since hey only consider qualiy deerioraion mechanisms. One of he earlies models in his field is ha of Rahim and Banerjee (). Their objecive is he inegraed opimizaion of he SPC parameers and he PM ime in producion processes subjec o qualiy shifs, which resul in an inferior-qualiy, ye operaing, sae. Oher relevan models have been developed more recenly by Cassady e al.

11 Page of Inernaional Journal of Producion Research (000), Lee and Rahim (00) and Linderman e al. (00). Ben Daya and Rahim (000) sudy a similar problem incorporaing he noion of imperfec prevenive mainenance. However, hese papers do no ake ino accoun he possibiliy of a complee failure ha would enforce an immediae cease of operaion, which is a ypical elemen of pracically all mainenance problems and models.. Problem definiion, assumpions and noaion We consider a producion process ha may operae in one of wo possible qualiy saes; in-conrol sae or sae 0 and ou-of-conrol sae or sae. Regardless of he acual qualiy sae of he process, he equipmen may suffer a failure a any ime, resuling in complee soppage of operaion. The failure rae of he process in boh qualiy saes is a non-decreasing funcion of he equipmen age. The failure rae in he ou-of-conrol sae is assumed o be higher han ha of he in-conrol sae for he same equipmen age. Apar from ha, he wo qualiy saes also differ in erms of producion revenues; he ou-of-conrol sae is assumed o be less profiable han he in-conrol sae. Equal failure raes or equal producion revenues are easily reaed as special cases. Qualiy shifs from he in-conrol o he ou-of-conrol sae may occur a any ime bu hey are immediaely observed and consequenly he acual qualiy sae of he process is always known wih accuracy. Noe ha he assumpion of coninuous and accurae knowledge of he acual process qualiy sae is no characerisic of ypical SPC problems. However, i is ofen realisic since here are many cases where he qualiy shif is direcly observable or even self announced, e.g. hrough some disinc noise. Besides, even when a qualiy shif has o be deeced hrough inspecion, he inspecion process is someimes performed on a coninuous basis and is very accurae; e.g. inspecion by means of on - line sensors.

12 Inernaional Journal of Producion Research Page 0 of Since he failure rae increases wih he age of he equipmen due o physical deerioraion even in he in-conrol sae, i is reasonable o prevenively mainain he equipmen when i reaches some criical age ( m0 ), beyond which he probabiliy of failure is unaccepably high. PM upgrades he equipmen in he as-good-as-new condiion and oally renews he process. Alhough he frequency of failures can be radically decreased hrough PM, he equipmen will ineviably fail occasionally and producion will be inerruped. CM is implemened as soon as a failure occurs and he equipmen is resored again o he as-good-as-new condiion. From he above descripion i is clear ha he producion-mainenance process consiss of a series of independen and sochasically idenical cycles. Each cycle begins wih he process in he as-good-as-new condiion (in-conrol sae and zero equipmen age) and erminaes eiher wih a PM a m0 or wih a CM following a failure, whichever occurs firs. Wihin each cycle he process may shif o he ou-of-conrol sae leading o decreased producion revenues and/or increased failure rae. To avoid producion under hese poor condiions i is possible o sop he process and bring he equipmen back o he in-conrol sae wihou affecing he equipmen age. In oher words, an MM can ake place o improve he qualiy sae of he process. I is worh noing ha such MM acions are usually assumed o be performed immediaely afer he process is deeced o operae a or above some hreshold sae. Our model is more flexible han ha, since he process may inenionally be allowed o operae for some ime in he ou-of-conrol sae before MM is performed. In general, MM is no necessarily preferable o operaion in he ou-of-conrol sae. Wheher MM is worhwhile or no depends on he radeoff beween is cos on one hand and he lower

13 Page of Inernaional Journal of Producion Research revenues and higher failure rae associaed wih ou-of-conrol operaion on he oher. We consider he following hree alernaive siuaions/policies:. There is a criical age of he equipmen, m (0< m < m0 ), beyond which operaion in he ou-of-conrol sae is uneconomical due o he unaccepably high failure rae. Tha is, if a qualiy shif occurs a ime > m hen MM is performed immediaely and he process coninues is operaion in he in-conrol sae. If a qualiy shif occurs a ime < m, he process is allowed o coninue in he ou-of-conrol sae unil m and i is only hen resored o he in-conrol sae (unless a failure occurs before m ).. Operaion in he ou-of-conrol sae is so cosly ha he process is no allowed o operae a all in his sae. In such cases MM is always performed as soon as a qualiy shif occurs ( m =0). Borrowing from he auomaic conrol erminology, a policy wih m =0 will be called Acive Qualiy Mainenance policy, or simply AQM.. The cos of MM is oo high relaive o is benefis and consequenly he equipmen is no resored o he in-conrol sae earlier han m0 ( m = m0 ). A policy wih m = m0 will be called Passive Qualiy Mainenance policy, or simply PQM. Noe ha when m < m0 (cases and above) he process may shif o he ou-ofconrol sae more han once in each producion cycle and as a resul MM will be performed several imes in a cycle. To summarize, he proposed mainenance policy is characerized by wo criical imes (equipmen ages) m and m0 ( 0 m m0 ) and he objecive is o find he

14 Inernaional Journal of Producion Research Page of opimal values of m0 and m ha maximize he expeced profi per ime uni. The noaion ha will be used o develop he opimizaion model is presened below: f() densiy funcion of he ime o qualiy shif F() cumulaive disribuion funcion of he ime o qualiy shif; F() = F( ) h( ) φ i () = f () F() (qualiy shif rae) densiy funcion of he ime of failure (equipmen age) if he process is in sae i (i=0,) a =0; noe ha he densiy funcion of he ime o failure if a qualiy shif occurs a ime s is φ () Φ ( s ) for > s Φ i () cumulaive disribuion funcion of he ime of failure in sae i; hi = i i m0 m Z Z P Z M R i W W P W M i i Φ () = Φ φ () Φ () (failure rae in sae i) equipmen age; =0 a he beginning of each cycle scheduled prevenive mainenance ime in he in-conrol sae scheduled minimal mainenance ime if a qualiy shif occurs a < m expeced ime o perform correcive mainenance expeced ime o perform prevenive mainenance expeced ime o perform minimal mainenance expeced ne revenue per ime uni of operaion in sae i cos of correcive mainenance cos of prevenive mainenance cos of minimal mainenance 0

15 Page of Inernaional Journal of Producion Research n Expeced number of MM acions in a cycle E(T) expeced cycle lengh E(P) expeced cycle profi EPT expeced profi per ime uni Noe ha R 0 and R ake ino accoun he poenial cos of low qualiy iems and he operaing cos. Consequenly i is reasonable o assume ha R R 0.. Model developmen In order o formulae he EPT funcion we firs develop he expressions for he expeced cycle lengh E(T) and he expeced cycle profi E(P). Each cycle begins wih zero equipmen age (=0) and erminaes eiher wih repair afer failure before m0 or wih prevenive mainenance a m0. In boh cases he process may never shif o he ou-ofconrol sae, may shif jus once or may shif several imes. Thus, he oal duraion of a cycle consiss of he following sub-periods: operaing ime in he in-conrol sae (random variable T 0 ) operaing ime in he ou-of-conrol sae, if a qualiy shif occurs prior o failure and before m (random variable T ) repair ime, if a failure occurs before m0 prevenive mainenance ime, if he equipmen reaches age m0 wihou failure ime for minimal mainenance acions. Similarly, he expeced profi per cycle consiss of he following componens: producion ne revenue in he in-conrol sae

16 Inernaional Journal of Producion Research Page of producion ne revenue in he ou-of-conrol sae, if a qualiy shif occurs prior o failure and before m repair cos, if a failure occurs before m0 prevenive mainenance cos, if he equipmen reaches age m0 wihou failure cos of minimal mainenance acions. Operaing ime in he in-conrol sae Once he equipmen reaches age m (eiher in he in-conrol or in he ou-of-conrol sae) i is no allowed o operae in he ou-of-conrol sae in he remainder of he cycle, alhough qualiy shifs may sill occur (bu will be immediaely removed hrough MM). Consequenly, he expeced in-conrol period during he cycle can be divided ino wo pars; one before m and one afer m. Regarding he firs par, he in-conrol period lass unil m only if neiher a failure nor a qualiy shif occurs by ha ime; oherwise, i lass unil some ime < m. Thus, he expeced duraion of he in-conrol period before m is m. Ε(Τ before ) = Φ F + φ F d+ f Φ d 0 m m 0 m m Noing ha φ ( ) d= -dφ ( ) and inegraing F( ) dφ ( ) 0 0 m 0 by pars yields 0 m Ε(Τ before ) = Φ F d. () 0 m 0 0 The in-conrol period afer m lass eiher unil m0 or unil an equipmen failure before m0. Thus, he expeced lengh of he in-conrol period afer m provided ha he equipmen survives unil m is m

17 Page of Inernaional Journal of Producion Research m 0 Φ0 m0 φ0 0 m ( m0 m) ( m) Φ0 m Φ 0 m Ε(Τ afer ) = + - d. Again, using φ ( ) d -dφ ( ) 0 0 = and inegraing by pars yields m m 0 m 0 Φ Ε(Τ0 afer m) = d. () Φ 0 m The probabiliy ha he equipmen will survive unil m wihou a failure is Φ = +. () m m p Φ m 0 m F m f Φ0 d Φ 0 The firs erm in he righ hand side of () expresses he probabiliy ha neiher a failure nor a qualiy shif occurs before m (he process remains in he in-conrol sae unil m ), while he second erm expresses he probabiliy ha a qualiy shif occurs prior o m, ye no failure occurs unil m (he operaing ime in he ou-of-conrol sae is sricly posiive). The oal expeced ime ha he process operaes in he in-conrol sae in a cycle is given by m E T = Ε(Τ before ) + p Ε(Τ afer ), 0 0 m 0 m which, combining () hrough () and simplifying resuls in = + + m m 0 m m 0 Φ m Φ0 0 0 m 0 0 Φ 0 m 0 Φ m 0 m. () E T Φ F d F Φ d f Φ d d Operaing ime in he ou-of-conrol sae The process spends some ime in he ou-of-conrol sae if and only if a qualiy shif occurs prior o m. Then, he operaing ime in he ou-of-conrol sae lass unil m if no failure occurs by ha ime or unil an equipmen failure before m. Thus, he

18 Inernaional Journal of Producion Research Page of expeced period ha he process operaes in he ou-of-conrol sae in a cycle is given by: ( - )φ Φ Ε(Τ ) f Φ d d - f ()Φ d m m m m = 0 + ( m ) 0 Φ 0 Φ 0. () Inegraing by pars he second inegral of he firs erm in he righ hand side of () yields ( ) Φ m m Ε(Τ ) f( Φ ) 0( ) d d Φ 0 =. () Probabiliy of prevenive mainenance and probabiliy of failure Prevenive mainenance is only performed whenever he equipmen reaches age m0 wihou a failure. The probabiliy ha he equipmen will reach age m wihou a failure, p m, is given by (), while he probabiliy ha he equipmen will no fail in he ime inerval from m o m0, provided ha i survived unil m, is Φ0( m0) Φ0( m). The fac ha he process may shif (once or more) o he ou-of-conrol sae afer m has no effec on he failure probabiliy of he process since he equipmen is immediaely resored o he in-conrol sae. Thus, he probabiliy of prevenive mainenance in a cycle is p p Φ F f Φ d m Φ0 m0 Φ m Φ0 m0 PM = = m 0 m m + 0 Φ0 m Φ 0 Φ0 m = Φ Φ Φ F + f Φ d. () m 0 m0 m 0( m0) ( m) 0 Φ0 m Φ 0 The probabiliy of failure in a cycle is simply -p PM.

19 Page of Inernaional Journal of Producion Research Expeced number of minimal mainenance acions The earlies ime ha a minimal mainenance can ake place is m ; his happens if he equipmen reaches age m in he ou-of-conrol sae. Therefore, he expeced number of minimal mainenance acions by m, denoed n, is equal o he probabiliy of ha even Φ =. () m m n f( Φ ) 0( ) d Φ 0 The equipmen is also minimally mainained every ime a qualiy shif occurs afer m. The expeced number of minimal mainenance acions afer m, provided ha he equipmen survived unil m, is denoed n and compued by he following lemma, he exensive proof of which is presened in he Appendix. Lemma The expeced number of minimal mainenance acions afer m in a cycle, provided ha he equipmen survives unil m, is n 0 m0 = h d. () m Φ Φ 0 m The expeced oal number of minimal mainenance acions hroughou a cycle is n= n + p n. Combining equaions (), () and () and simplifying we obain m Φ Φ Φ n= f Φ d+ F Φ h d+ f Φ d h d m m 0 m m0 m m 0 0 ( m) 0 0 Φ Φ 0 m 0 Φ m 0 m. Expeced profi per ime uni The oal expeced cycle lengh is given by Ε(Τ) = Ε(Τ ) + E(T ) + Z p + Z - p + Z n, 0 P PM PM M while he oal expeced profi per cycle is given by (0)

20 Inernaional Journal of Producion Research Page of Ε(P) = R Ε(Τ ) + R E(T ) - W p - W - p W n. 0 0 P PM PM M Finally, since he process is a renewal reward process, he expeced profi per ime uni can be expressed as he raio of he expeced profi per cycle o he expeced cycle lengh: EPT (, ) m0 m E(P) =. E(T) Noe ha he model holds for any value of m in he inerval [0, m0 ] including he special cases m =0 (immediae MM afer a qualiy shif; acive qualiy mainenance policy) and m = m0 (he equipmen is allowed o operae in he ou-of-conrol sae unil m0 wihou any MM inervenion; passive qualiy mainenance policy). In he laer case, if no failure occurs and he equipmen succeeds o survive unil m0 operaing in he ou-of-conrol sae, he scheduled PM a m0 is combined wih an MM, wih oal cos W M + W P and oal ime Z M + Z P, so ha he nex producion cycle begins wih he process in a perfec condiion (equipmen as-good-as-new and in-conrol).. Wha is (usually) he ype of he opimal policy? The mainenance model developed in he previous secion allows m o be anywhere in he range [0, m0 ]. For given m0, he model essenially compares he cos of MM o upgrade he process qualiy o he in-conrol sae agains he benefis of such an acion and reurns he opimal value of m. However, he mainenance managemen of mos acual producion sysems eiher adops an Acive Qualiy Mainenance (AQM) policy, whereby he equipmen is never allowed o operae in he ou-of-conrol sae ( m =0) or i employs a Passive Qualiy Mainenance (PQM) policy, whereby he equipmen is reaed he same regardless of is qualiy sae ( m = m0 ). The choice beween AQM and PQM is ypically made empirically; if he ou-of-conrol sae is considered

21 Page of Inernaional Journal of Producion Research unaccepable hen AQM is adoped, while if process inerrupions are considered cosly and undesirable hen PQM is preferred. To invesigae he acual ype of he opimal policy and evaluae he qualiy of AQM and PQM policies we solved a large number of numerical examples wih differen problem parameers. Specifically, we solved 0 examples wih all process parameers randomly seleced from a wide range of values. In addiion, we have allowed boh Weibull and Gamma disribuions o describe he process failure and qualiy shif mechanisms. The densiy funcions of hese disribuions are presened in Table. Noe ha ensuring ha he failure rae when operaing in he ou-of-conrol sae wih equipmen age is a leas as large as he failure rae in he in-conrol sae wih he same equipmen age, requires c 0 =c and λ0 λ for boh disribuion ypes. The allowable values of all process parameers are shown in Table. [Inser Tables and abou here] The opimizaion is performed by means of an exhausive search over all possible m0 and m values, under he consrain 0 m m0, o ensure ha he global opimum of EPT is obained. For compuaional simpliciy we resric our numerical invesigaion o ineger values of m0 and m wih iniial values equal o zero. The opimizaion algorihm compues he expeced profi for all m0 and m values and he search procedure sops as soon as m0 and m reach some hreshold values beyond which heir effec on he expeced profi is insignifican. The convergence of EPT is assured by he fac ha he probabiliy of failure will evenually reach uniy for large values of m0. The compuaional ime required for finding he opimal soluion ypically varies beween 0 and 0 minues on a Penium IV. GHz personal compuer, depending on he parameers of he example.

22 Inernaional Journal of Producion Research Page 0 of The main findings are summarized below: In all 0 cases, he opimal policy is eiher AQM ( m =0) or PQM ( m = m0 ). In he vas majoriy of he cases examined ( ou of 0) he opimal policy is AQM. Among he cases where AQM is no opimal, he percenage loss ha would resul from using he bes possible AQM insead of he opimal PQM is.% on average, ranging beween 0.0% and.%. Among he cases where PQM is no opimal, he percenage loss ha would resul from using he bes possible PQM is.% on average, ranging beween 0.0% and.%. When he failure imes follow Gamma disribuions, he opimal soluion usually dicaes ha he equipmen should no undergo PM a all ( m0 ). This can be explained by he fac ha he failure rae of a Gamma disribuion is sabilized as he equipmen age grows large (similar o he memoryless exponenial disribuion) and i may be beer o coninue operaion from some poin on raher han mainain he equipmen. I is imporan o menion here ha alhough he opimal soluion in all numerical examples solved in he course of his research (including he 0 examples of his secion and he examples of he nex secion) is obained for an exreme value of m (eiher 0 or m0 ), in general he opimal soluion is no necessarily unique. In some cases he expeced profi funcion EPT is so fla wihin a range of m values (from m =0 o some criical m value) ha pracically all m in his range can be considered opimal. In hese cases i may be preferable o choose he larges m in he opimal range so as o minimize downimes due o MM acions.

23 Page of Inernaional Journal of Producion Research Also noe ha our exensive numerical invesigaion of he behaviour of EPT has shown ha an exreme value of m (eiher m =0 or m = m0 ) is always opimal, even when here are muliple opimal soluions wih 0< m < m0. Alhough a formal proof of such a propery, if i is indeed rue in general, remains elusive, he pracical implicaion is ha i is sufficien o search for he opimal soluion only beween he wo exreme policies AQM and PQM. In his way he compuaional requiremens are significanly reduced since i suffices o opimize wo single-variable funcions insead of searching over a wo-dimensional decision space.. The effec of process parameers on he opimal policy To invesigae sysemaically he effec of he process parameers on he opimal policy we have solved problems, differing subsanially in key model feaures such as he relaive coss of operaion and mainenance aciviies and he equipmen proneness o qualiy shifs and failures. Specifically, we express he qualiy shif mechanism and he failure mechanism in boh qualiy saes (i=0,) by Weibull disribuions of he equipmen age. The parameers λ, λ, R and W P are examined a levels each as shown in Table, while he pair W M, Z M is examined a hree levels as follows: (a) W M =0.W P and Z M =0.Z P (b) W M =0.W P and Z M =0.Z P (c) W M =0.W P and Z M =0.Z P The remaining parameers are se equal o he following values in all cases: c 0 =c =, c=., λ 0 =0.00, R 0 =00, Z=Z P =.0, W=00. The cases are numbered a, b, c o a, b, c where a, b or c indicaes he W M, Z M combinaion. [Inser Table abou here]

24 Inernaional Journal of Producion Research Page of Table shows he opimal criical equipmen ages m, m0 and he corresponding maximum expeced profi per ime uni, EPT, for he cases. Noe ha, exacly as in he 0 examples of he previous secion, in all cases he opimal policy is eiher AQM ( cases) or PQM ( cases). In of he cases where a passive qualiy mainenance (PQM) policy is opimal, i is opimal o only use correcive mainenance upon failure and never resor o PM or MM acions ( m = m0 = ). The las columns of Table presen he opimal AQM and he opimal PQM policies along wih he percenage losses when using each policy. Obviously, when one of hese policies is opimal is respecive percenage loss is zero. [Inser Table abou here] The effecs of he parameers on he ype of he opimal policy as well as on he savings ha can be achieved by using he proposed model as opposed o blindly following he AQM or he PQM policy are summarized as follows: Large R ( R R 0 ) and/or small λ ( λ λ 0 ) implies ha sae (ou-ofconrol) is no ha inferior o sae 0 (in-conrol) in erms of profi and/or failure rae. Consequenly i may be more economical o allow operaion in sae raher han upgrade he process qualiy o sae 0. Thus, PQM ends o ouperform AQM in such cases. The percenage loss associaed wih he use of an AQM policy when PQM is opimal increases as R increases and/or as λ decreases. This is because increasing R or decreasing λ increases he expeced profi per ime uni for m = m0 (PQM policy), while he expeced profi per ime uni for m = 0 remains unalered (AQM policy); consequenly he wo soluions diverge. 0

25 Page of Inernaional Journal of Producion Research Combinaions b (W M =0.W P, Z M =0.Z P ) and c (W M =0.W P, Z M =0.Z P ) are characerized by a higher direc or indirec (downime) cos of MM han combinaion a (W M =0.W P, Z M =0.Z P ). Consequenly, a PQM policy is more likely o be opimal in combinaions b and c. The percenage loss due o incorrecly adoping an AQM policy (while PQM is opimal) increases when moving from combinaion a of W M and Z M o combinaion b or c, while he percenage loss associaed wih incorrec adopion of a PQM policy decreases. This is because he increased cos of MM in combinaions b and c ends o have a greaer negaive impac on he soluions wih m = 0 han on hose wih m = m0. Large W P and/or large λ signify more expensive or more frequen MM acions and consequenly a PQM policy is more likely o be opimal. The percenage loss associaed wih incorrec use of an AQM policy increases as W P and/or λ increase, due o he negaive effec of boh hese parameers on he oal cos of MM acions. However, he percenage loss associaed wih incorrec use of a PQM policy can eiher increase or decrease as W P and/or λ increase. We have also sudied he isolaed effec of W P, keeping W M consan; in all cases he ype of he opimal policy remained unalered. Finally, he effec of failure ime variabiliy has also been invesigaed using c0= c= insead of c0 = c = in all cases and modifying he values of λ 0 and λ, so as o keep he mean imes o failure unalered for boh qualiy saes. This invesigaion has no revealed any sysemaic effec of he failure ime variabiliy on he

26 Inernaional Journal of Producion Research Page of ype of he opimal policy nor on he percenage loss associaed wih consisenly adoping AQM or PQM.. Conclusions In his paper we have sudied mainenance procedures in a producion process subjec o boh qualiy shifs and failures. The equipmen deerioraes coninuously due o he ageing process and a he same ime i may experience a jump o an inferior qualiy sae (ou-of-conrol sae) upon he occurrence of an assignable cause. Transiions o he ouof-conrol sae have a dual impac on he process; hey resul in lower producion qualiy implying lower producion revenues bu hey also increase he failure rae of he process. The proposed mainenance policy is a combinaion of an age-based prevenive mainenance policy a some criical age m0, wih addiional minimal mainenance acions, which upgrade he process qualiy o he in-conrol sae. Such MM acions (qualiy adjusmens) are commonly used in qualiy conrol o eliminae he negaive effecs of assignable causes. In our model resoring he process o he in-conrol sae afer he occurrence of an assignable cause no only improves producion qualiy bu also decreases he failure rae of he process. In conras o ypical qualiy conrol models, he MM considered here is no necessarily implemened immediaely afer qualiy shifs bu may inenionally be posponed unil some criical equipmen age. Our invesigaion showed ha in pracically every case he opimal mainenance policy eiher calls for immediae MM as soon as a qualiy shif occurs (acive qualiy mainenance, AQM) or allows operaion in he ou-of-conrol sae and an MM is only implemened along wih prevenive mainenance a predeermined imes (passive qualiy mainenance, PQM). Thus, i suffices o consider only he opimal AQM and

27 Page of Inernaional Journal of Producion Research PQM policies, compare hem and use he one ha is mos effecive in each paricular case. Neverheless, using he wrong exreme policy may resul in significan loss. Our numerical invesigaion has shown ha such losses can be as high as abou 0% of he opimal expeced profi. There are some ineresing exensions of he proposed mainenance model ha are worh sudying, such as he case of incomplee informaion abou he qualiy sae of he process. In addiion, full inegraion of qualiy conrol procedures wih equipmen mainenance in deerioraing producion processes under general (non-resricive) assumpions would be of grea pracical ineres. Appendix: Proof of Lemma We firs deermine he probabiliy ha exacly n (n=,,..) qualiy shifs occur in he inerval ( m, m0 ) provided ha he equipmen survives unil ime/age m. We sar wih he derivaion for n= and hen generalize for higher values of n. The process shifs exacly once from he in-conrol o he ou-of-conrol sae (n=) if eiher one of he following wo scenarios maerializes: a) The process shifs o he ou-of-conrol sae a ime ( m < < m0 ) prior o failure, an MM is immediaely implemened upgrading he process qualiy o he in-conrol sae and operaion coninues wihou inervenion (neiher a failure nor a qualiy shif occurs) unil m0. b) The process shifs o he ou-of-conrol sae a ime ( m < < m0 ) prior o failure, an MM is immediaely implemened upgrading he process qualiy o he in-conrol sae and operaion coninues unil he occurrence of a failure a ime before m0. No oher qualiy shifs occur beween and.

28 Inernaional Journal of Producion Research Page of Under hese scenarios he probabiliy of a qualiy shif a prior o failure provided ha he equipmen survives unil m is p = f Φ. (A) m 0 0 qs d F m m Φ0 m The probabiliy ha neiher a failure nor a qualiy shif occurs in (, m0 ) is p Φ F = 0 m0 m0 a Φ0 F, (A) while he probabiliy of failure in (, m0 ) provided ha he equipmen operaes coninuously in he in-conrol sae is p φ F m0 = 0 b d Φ 0 F. (A) Combining (A) hrough (A) he probabiliy ha he process shifs exacly once from he in-conrol o he ou-of-conrol sae during he inerval ( m, m0 ) is given by P = p p + p qs a b f Φ φ F f Φ Φ F = + m 0 m 0 m m0 m0 dd d F m m Φ0 m Φ 0 F F m m Φ0 m Φ0 F. Using h( ) f( ) F( ) = and simplifying yields m 0 m0 m 0 P = h( ) φ0( ) F( ) dd+ Φ0( m0) F( m0) h( ) d F( m) Φ0( m) m m. Reversing he order of inegraion in he double inegral resuls in m 0 m0 P = φ0( ) F( ) h( ) dd Φ0( m0) F( m0) h( ) d F( m) Φ0( m) + m m m (A)

29 Page of Inernaional Journal of Producion Research = m 0 m where H( ) + φ0 F H H m d Φ0 m0 F m0 H m0 H m = h x dx. 0 Φ ( ) F m 0 m Similarly, he process shifs exacly wice (n=) from he in-conrol o he ou-ofconrol sae during he inerval ( m, m0 ) whenever a qualiy shif occurs a ime ( m < < m0 ) prior o failure (probabiliy p qs ), an MM is immediaely implemened upgrading he process qualiy o he in-conrol sae and hen he process shifs exacly once more (a ime/age > ) o he ou-of-conrol sae during he res of he cycle (, m0 ). Since he probabiliy of ha las even is analogous o P of (A) bu refers o he inerval (, m0 ), i is denoed P ( ). Thus, properly adaping (A) o he inerval (, m0 ) leads o he following expression for he probabiliy P of he even n=: = we ge f Φ f Φ φ F P = p P = dd d + m 0 m0 m qs ( ) F m m Φ0 m F Φ0 Φ 0 F m 0 f Φ f Φ Φ F m m0 m0 F Φ F Φ Φ F m m 0 m 0 0 d d m 0 m0 m 0 m 0 m0 + h h φ F dd d Φ F h h d d 0 0 m0 m0 m m Φ ( ) F m 0 m Reversing he order of inegraion in boh inegrals and hen using b a b H H a H b H a h( ) H( ) H( a) d = = a.

30 Inernaional Journal of Producion Research Page of = P = m 0 m0 + φ F h h d d d Φ F h h d d 0 0 m0 m0 m m m m m Φ ( ) F m 0 m m 0 m 0 + φ F h H H d d Φ F h H H d 0 m 0 m0 m0 m m m m = m 0 m Φ ( ) F m 0 m Φ ( ) H H m H m0 H m φ0( ) F( ) d+ Φ0( m0) F( m0) F m 0 m Exending he preceding analysis o n> we arrive (by inducion) a he following general expression for he probabiliy P n ha he process shifs exacly n imes from he in-conrol o he ou-of-conrol sae during he inerval ( m, m0 ): = and m 0 m P n = m 0 m n Φ ( ) m 0 m H H m H m0 H m φ0( ) F( ) d+ Φ0( m0) F( m0) n! n! F The expeced number of MM acions in he inerval ( m, m0 ) is = npn n= n Φ ( ) n n H H m H m0 H m 0 + 0( m0) ( m0) n= n! n= n! φ F n d Φ F n F Using successively he facs ha e H = F n= n= m 0 m. n n x x n = x = xe n! n! ( ) x n..

31 Page of Inernaional Journal of Producion Research he expression for n becomes n = m 0 m H( ) H( m) H( m 0) H( m) 0 m + 0 m0 m0 m0 m φ F H H e d Φ F H H e Φ ( ) F m 0 m = + Φ m 0 φ0( ) H( ) H( m) d Φ0( m0) H( m0) H( m) 0( m). m Finally, using φ ( ) d -dφ ( ) expression of Lemma : 0 0 = and inegraing by pars resuls in he simple n m 0 m 0 Φ = h( ) d. Φ 0 m References Ben-Daya, M. and Rahim, M.A., Effec of Mainenance on he Economic Design of X - Conrol Char. European Journal of Operaional Research, 000, 0, -. Cassady, C.R., Bowden, R.O., Liew, L. and Pohl, E.A., Combining Prevenive Mainenance and Saisical Process Conrol: a Preliminary Invesigaion. IIE Transacions, 000,, -. Chiang, G. H. and Yuan, J., Opimal mainenance policy for a Markovian sysem under periodic inspecion. Reliabiliy Engineering and Sysem Safey, 00,, -. Derman, C., Opimal replacemen and mainenance under Markovian deerioraion wih probabiliy bounds on failure. Managemen Science,,, -. Lee, B.H. and Rahim, M.A., An inegraed economic design model for qualiy conrol, replacemen, and mainenance. Qualiy Engineering, 00,, -.

32 Inernaional Journal of Producion Research Page 0 of Linderman, K., McKone-Swee, K.E. and Anderson, J.C., An inegraed sysems approach o process conrol and mainenance. European Journal of Operaional Research, 00,, -0. Makis, V. and Fung, J., Opimal prevenive replacemen, lo sizing and inspecion policy for a deerioraing producion sysem. Journal of Qualiy in Mainenance Engineering,,, -. Makis, V. and Fung, J., An EMQ model wih inspecions and random machine failures. Journal of he Operaional Research Sociey,,, -. Mousafa, M. S., Abdel Maksoud, E. Y. and Sadek, S., Opimal major and minimal mainenance policies for deerioraing sysems. Reliabiliy Engineering and Sysem Safey, 00,, -. Özekici, S. and Günlük, N. O., Mainenance of a device wih age-dependen exponenial failures. Naval Research Logisics,,, -. Pham, H. and Wang, H., Imperfec mainenance. European Journal of Operaional Research,,, -. Pierskalla, W. P. and Voelker, J. A., A survey of mainenance models: The conrol and surveillance of deerioraing sysems. Naval Research Logisics Quarerly,,, -. Rahim, M.A. and Banerjee, P.K., A generalized model for he economic design of X - conrol chars for producion sysems wih increasing failure rae and early replacemen. Naval Research Logisics,, 0, -0. Tagaras, G., An inegraed cos model for he join opimizaion of process conrol and mainenance. Journal of he Operaional Research Sociey,,, -.

33 Page of Inernaional Journal of Producion Research Valdez-Flores, C. and Feldman, R.M., A survey of Prevenive Mainenance models for sochasically deerioraing single-uni sysems. Naval Research Logisics Quarerly,,, -.

34 Inernaional Journal of Producion Research Page of Table. Densiy funcions for qualiy shif and failure mechanisms Disribuion Qualiy shif mechanism Failure mechanism in sae i (i = 0,) Weibull Gamma c c λ f = λc e, >0, c>0, λ>0 >0, c i >0, λ i >0 c c λ f = λ e Γ c, >0, c>0, λ>0 c c i i λi i i i φ = λ c e, c i c i λ i φ = λ e Γ c, i i i >0, c i >0, λ i >0 0

35 Page of Inernaional Journal of Producion Research Table. Range of process parameer values Process parameer Range of values R R 0 - R 0 W W P W M 0.W - W 0.W P - W P Z = Z P Z M φ i : Weibull (mean µ i ) f: Weibull 0.Z - Z c i. - (0.) µ 0-0 µ 0 - µ 0 c. - (0.) (mean µ) µ - 0 φ i : Gamma (mean µ i ) f: Gamma c i - - µ 0-0 µ 0 - µ 0 c (mean µ) µ - 0

36 Inernaional Journal of Producion Research Page of Table. Parameer values for he numerical invesigaion λ λ R W P

37 Page of Inernaional Journal of Producion Research Se Table. Opimal soluions and evaluaion of he AQM and PQM policies Parameers λ λ R W P MM Opimum AQM Loss PQM Loss Policy m m0 EPT m0 (%) ( m =0) m = m0 (%) a b c a b.. 0 c a b c a b c a b c a b.. 0 c a b c a.. 0 b c a b c a b c a b c a b c a b c a b.0. 0 c a b c a b.. 0 c

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

Inventory Control of Perishable Items in a Two-Echelon Supply Chain

Inventory Control of Perishable Items in a Two-Echelon Supply Chain Journal of Indusrial Engineering, Universiy of ehran, Special Issue,, PP. 69-77 69 Invenory Conrol of Perishable Iems in a wo-echelon Supply Chain Fariborz Jolai *, Elmira Gheisariha and Farnaz Nojavan

More information

Evaluation of Mean Time to System Failure of a Repairable 3-out-of-4 System with Online Preventive Maintenance

Evaluation of Mean Time to System Failure of a Repairable 3-out-of-4 System with Online Preventive Maintenance American Journal of Applied Mahemaics and Saisics, 0, Vol., No., 9- Available online a hp://pubs.sciepub.com/ajams/// Science and Educaion Publishing DOI:0.69/ajams--- Evaluaion of Mean Time o Sysem Failure

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

Probabilistic Models for Reliability Analysis of a System with Three Consecutive Stages of Deterioration

Probabilistic Models for Reliability Analysis of a System with Three Consecutive Stages of Deterioration Yusuf I., Gaawa R.I. Volume, December 206 Probabilisic Models for Reliabiliy Analysis of a Sysem wih Three Consecuive Sages of Deerioraion Ibrahim Yusuf Deparmen of Mahemaical Sciences, Bayero Universiy,

More information

Maintenance Models. Prof. Robert C. Leachman IEOR 130, Methods of Manufacturing Improvement Spring, 2011

Maintenance Models. Prof. Robert C. Leachman IEOR 130, Methods of Manufacturing Improvement Spring, 2011 Mainenance Models Prof Rober C Leachman IEOR 3, Mehods of Manufacuring Improvemen Spring, Inroducion The mainenance of complex equipmen ofen accouns for a large porion of he coss associaed wih ha equipmen

More information

Errata (1 st Edition)

Errata (1 st Edition) P Sandborn, os Analysis of Elecronic Sysems, s Ediion, orld Scienific, Singapore, 03 Erraa ( s Ediion) S K 05D Page 8 Equaion (7) should be, E 05D E Nu e S K he L appearing in he equaion in he book does

More information

CENTRALIZED VERSUS DECENTRALIZED PRODUCTION PLANNING IN SUPPLY CHAINS

CENTRALIZED VERSUS DECENTRALIZED PRODUCTION PLANNING IN SUPPLY CHAINS CENRALIZED VERSUS DECENRALIZED PRODUCION PLANNING IN SUPPLY CHAINS Georges SAHARIDIS* a, Yves DALLERY* a, Fikri KARAESMEN* b * a Ecole Cenrale Paris Deparmen of Indusial Engineering (LGI), +3343388, saharidis,dallery@lgi.ecp.fr

More information

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

International Journal of Computer Science Trends and Technology (IJCST) Volume 3 Issue 6, Nov-Dec 2015

International Journal of Computer Science Trends and Technology (IJCST) Volume 3 Issue 6, Nov-Dec 2015 Inernaional Journal of Compuer Science Trends and Technology (IJCST) Volume Issue 6, Nov-Dec 05 RESEARCH ARTICLE OPEN ACCESS An EPQ Model for Two-Parameer Weibully Deerioraed Iems wih Exponenial Demand

More information

On Multicomponent System Reliability with Microshocks - Microdamages Type of Components Interaction

On Multicomponent System Reliability with Microshocks - Microdamages Type of Components Interaction On Mulicomponen Sysem Reliabiliy wih Microshocks - Microdamages Type of Componens Ineracion Jerzy K. Filus, and Lidia Z. Filus Absrac Consider a wo componen parallel sysem. The defined new sochasic dependences

More information

On a Discrete-In-Time Order Level Inventory Model for Items with Random Deterioration

On a Discrete-In-Time Order Level Inventory Model for Items with Random Deterioration Journal of Agriculure and Life Sciences Vol., No. ; June 4 On a Discree-In-Time Order Level Invenory Model for Iems wih Random Deerioraion Dr Biswaranjan Mandal Associae Professor of Mahemaics Acharya

More information

Applying Genetic Algorithms for Inventory Lot-Sizing Problem with Supplier Selection under Storage Capacity Constraints

Applying Genetic Algorithms for Inventory Lot-Sizing Problem with Supplier Selection under Storage Capacity Constraints IJCSI Inernaional Journal of Compuer Science Issues, Vol 9, Issue 1, No 1, January 2012 wwwijcsiorg 18 Applying Geneic Algorihms for Invenory Lo-Sizing Problem wih Supplier Selecion under Sorage Capaciy

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

Solutions Problem Set 3 Macro II (14.452)

Solutions Problem Set 3 Macro II (14.452) Soluions Problem Se 3 Macro II (14.452) Francisco A. Gallego 04/27/2005 1 Q heory of invesmen in coninuous ime and no uncerainy Consider he in nie horizon model of a rm facing adjusmen coss o invesmen.

More information

A Study of Inventory System with Ramp Type Demand Rate and Shortage in The Light Of Inflation I

A Study of Inventory System with Ramp Type Demand Rate and Shortage in The Light Of Inflation I Inernaional Journal of Mahemaics rends and echnology Volume 7 Number Jan 5 A Sudy of Invenory Sysem wih Ramp ype emand Rae and Shorage in he Ligh Of Inflaion I Sangeea Gupa, R.K. Srivasava, A.K. Singh

More information

5. Stochastic processes (1)

5. Stochastic processes (1) Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

This document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC

This document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC his documen was generaed a 1:4 PM, 9/1/13 Copyrigh 213 Richard. Woodward 4. End poins and ransversaliy condiions AGEC 637-213 F z d Recall from Lecure 3 ha a ypical opimal conrol problem is o maimize (,,

More information

An Inventory Model for Time Dependent Weibull Deterioration with Partial Backlogging

An Inventory Model for Time Dependent Weibull Deterioration with Partial Backlogging American Journal of Operaional Research 0, (): -5 OI: 0.593/j.ajor.000.0 An Invenory Model for Time ependen Weibull eerioraion wih Parial Backlogging Umakana Mishra,, Chaianya Kumar Tripahy eparmen of

More information

1. An introduction to dynamic optimization -- Optimal Control and Dynamic Programming AGEC

1. An introduction to dynamic optimization -- Optimal Control and Dynamic Programming AGEC This documen was generaed a :45 PM 8/8/04 Copyrigh 04 Richard T. Woodward. An inroducion o dynamic opimizaion -- Opimal Conrol and Dynamic Programming AGEC 637-04 I. Overview of opimizaion Opimizaion is

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based

More information

Seminar 4: Hotelling 2

Seminar 4: Hotelling 2 Seminar 4: Hoelling 2 November 3, 211 1 Exercise Par 1 Iso-elasic demand A non renewable resource of a known sock S can be exraced a zero cos. Demand for he resource is of he form: D(p ) = p ε ε > A a

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

Shiva Akhtarian MSc Student, Department of Computer Engineering and Information Technology, Payame Noor University, Iran

Shiva Akhtarian MSc Student, Department of Computer Engineering and Information Technology, Payame Noor University, Iran Curren Trends in Technology and Science ISSN : 79-055 8hSASTech 04 Symposium on Advances in Science & Technology-Commission-IV Mashhad, Iran A New for Sofware Reliabiliy Evaluaion Based on NHPP wih Imperfec

More information

Mathematical Theory and Modeling ISSN (Paper) ISSN (Online) Vol 3, No.3, 2013

Mathematical Theory and Modeling ISSN (Paper) ISSN (Online) Vol 3, No.3, 2013 Mahemaical Theory and Modeling ISSN -580 (Paper) ISSN 5-05 (Online) Vol, No., 0 www.iise.org The ffec of Inverse Transformaion on he Uni Mean and Consan Variance Assumpions of a Muliplicaive rror Model

More information

A Dynamic Model of Economic Fluctuations

A Dynamic Model of Economic Fluctuations CHAPTER 15 A Dynamic Model of Economic Flucuaions Modified for ECON 2204 by Bob Murphy 2016 Worh Publishers, all righs reserved IN THIS CHAPTER, OU WILL LEARN: how o incorporae dynamics ino he AD-AS model

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

An Improved Adaptive CUSUM Control Chart for Monitoring Process Mean

An Improved Adaptive CUSUM Control Chart for Monitoring Process Mean An Improved Adapive CUSUM Conrol Char for Monioring Process Mean Jun Du School of Managemen Tianjin Universiy Tianjin 37, China Zhang Wu, Roger J. Jiao School of Mechanical and Aerospace Engineering Nanyang

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8) I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression

More information

Christos Papadimitriou & Luca Trevisan November 22, 2016

Christos Papadimitriou & Luca Trevisan November 22, 2016 U.C. Bereley CS170: Algorihms Handou LN-11-22 Chrisos Papadimiriou & Luca Trevisan November 22, 2016 Sreaming algorihms In his lecure and he nex one we sudy memory-efficien algorihms ha process a sream

More information

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology

More information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

More information

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy,

More information

) were both constant and we brought them from under the integral.

) were both constant and we brought them from under the integral. YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha

More information

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin ACE 56 Fall 005 Lecure 4: Simple Linear Regression Model: Specificaion and Esimaion by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Simple Regression: Economic and Saisical Model

More information

1. An introduction to dynamic optimization -- Optimal Control and Dynamic Programming AGEC

1. An introduction to dynamic optimization -- Optimal Control and Dynamic Programming AGEC This documen was generaed a :37 PM, 1/11/018 Copyrigh 018 Richard T. Woodward 1. An inroducion o dynamic opimiaion -- Opimal Conrol and Dynamic Programming AGEC 64-018 I. Overview of opimiaion Opimiaion

More information

3.1 More on model selection

3.1 More on model selection 3. More on Model selecion 3. Comparing models AIC, BIC, Adjused R squared. 3. Over Fiing problem. 3.3 Sample spliing. 3. More on model selecion crieria Ofen afer model fiing you are lef wih a handful of

More information

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070

More information

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN Inernaional Journal of Scienific & Engineering Research, Volume 4, Issue 10, Ocober-2013 900 FUZZY MEAN RESIDUAL LIFE ORDERING OF FUZZY RANDOM VARIABLES J. EARNEST LAZARUS PIRIYAKUMAR 1, A. YAMUNA 2 1.

More information

An introduction to the theory of SDDP algorithm

An introduction to the theory of SDDP algorithm An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

Energy Storage Benchmark Problems

Energy Storage Benchmark Problems Energy Sorage Benchmark Problems Daniel F. Salas 1,3, Warren B. Powell 2,3 1 Deparmen of Chemical & Biological Engineering 2 Deparmen of Operaions Research & Financial Engineering 3 Princeon Laboraory

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems. di ernardo, M. (995). A purely adapive conroller o synchronize and conrol chaoic sysems. hps://doi.org/.6/375-96(96)8-x Early version, also known as pre-prin Link o published version (if available):.6/375-96(96)8-x

More information

A First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18

A First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18 A Firs ourse on Kineics and Reacion Engineering lass 19 on Uni 18 Par I - hemical Reacions Par II - hemical Reacion Kineics Where We re Going Par III - hemical Reacion Engineering A. Ideal Reacors B. Perfecly

More information

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,

More information

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves Rapid Terminaion Evaluaion for Recursive Subdivision of Bezier Curves Thomas F. Hain School of Compuer and Informaion Sciences, Universiy of Souh Alabama, Mobile, AL, U.S.A. Absrac Bézier curve flaening

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A

More information

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017 Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =

More information

Introduction to Probability and Statistics Slides 4 Chapter 4

Introduction to Probability and Statistics Slides 4 Chapter 4 Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random

More information

CHAPTER 2 Signals And Spectra

CHAPTER 2 Signals And Spectra CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

Methodology. -ratios are biased and that the appropriate critical values have to be increased by an amount. that depends on the sample size.

Methodology. -ratios are biased and that the appropriate critical values have to be increased by an amount. that depends on the sample size. Mehodology. Uni Roo Tess A ime series is inegraed when i has a mean revering propery and a finie variance. I is only emporarily ou of equilibrium and is called saionary in I(0). However a ime series ha

More information

Math 333 Problem Set #2 Solution 14 February 2003

Math 333 Problem Set #2 Solution 14 February 2003 Mah 333 Problem Se #2 Soluion 14 February 2003 A1. Solve he iniial value problem dy dx = x2 + e 3x ; 2y 4 y(0) = 1. Soluion: This is separable; we wrie 2y 4 dy = x 2 + e x dx and inegrae o ge The iniial

More information

Single-Pass-Based Heuristic Algorithms for Group Flexible Flow-shop Scheduling Problems

Single-Pass-Based Heuristic Algorithms for Group Flexible Flow-shop Scheduling Problems Single-Pass-Based Heurisic Algorihms for Group Flexible Flow-shop Scheduling Problems PEI-YING HUANG, TZUNG-PEI HONG 2 and CHENG-YAN KAO, 3 Deparmen of Compuer Science and Informaion Engineering Naional

More information

Deteriorating Inventory Model When Demand Depends on Advertisement and Stock Display

Deteriorating Inventory Model When Demand Depends on Advertisement and Stock Display Inernaional Journal of Operaions Research Inernaional Journal of Operaions Research Vol. 6, No. 2, 33 44 (29) Deerioraing Invenory Model When Demand Depends on Adverisemen and Sock Display Nia H. Shah,

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

ACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H.

ACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H. ACE 56 Fall 5 Lecure 8: The Simple Linear Regression Model: R, Reporing he Resuls and Predicion by Professor Sco H. Irwin Required Readings: Griffihs, Hill and Judge. "Explaining Variaion in he Dependen

More information

Deteriorating Inventory Model with Time. Dependent Demand and Partial Backlogging

Deteriorating Inventory Model with Time. Dependent Demand and Partial Backlogging Applied Mahemaical Sciences, Vol. 4, 00, no. 7, 36-369 Deerioraing Invenory Model wih Time Dependen Demand and Parial Backlogging Vinod Kumar Mishra Deparmen of Compuer Science & Engineering Kumaon Engineering

More information

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

More information

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

More information

Random Walk with Anti-Correlated Steps

Random Walk with Anti-Correlated Steps Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

Comparing Means: t-tests for One Sample & Two Related Samples

Comparing Means: t-tests for One Sample & Two Related Samples Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion

More information

Solutions for Assignment 2

Solutions for Assignment 2 Faculy of rs and Science Universiy of Torono CSC 358 - Inroducion o Compuer Neworks, Winer 218 Soluions for ssignmen 2 Quesion 1 (2 Poins): Go-ack n RQ In his quesion, we review how Go-ack n RQ can be

More information

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

More information

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15. SMT Calculus Tes Soluions February 5,. Le f() = and le g() =. Compue f ()g (). Answer: 5 Soluion: We noe ha f () = and g () = 6. Then f ()g () =. Plugging in = we ge f ()g () = 6 = 3 5 = 5.. There is a

More information

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17 EES 16A Designing Informaion Devices and Sysems I Spring 019 Lecure Noes Noe 17 17.1 apaciive ouchscreen In he las noe, we saw ha a capacior consiss of wo pieces on conducive maerial separaed by a nonconducive

More information

Solutions to Odd Number Exercises in Chapter 6

Solutions to Odd Number Exercises in Chapter 6 1 Soluions o Odd Number Exercises in 6.1 R y eˆ 1.7151 y 6.3 From eˆ ( T K) ˆ R 1 1 SST SST SST (1 R ) 55.36(1.7911) we have, ˆ 6.414 T K ( ) 6.5 y ye ye y e 1 1 Consider he erms e and xe b b x e y e b

More information

DEPARTMENT OF STATISTICS

DEPARTMENT OF STATISTICS A Tes for Mulivariae ARCH Effecs R. Sco Hacker and Abdulnasser Haemi-J 004: DEPARTMENT OF STATISTICS S-0 07 LUND SWEDEN A Tes for Mulivariae ARCH Effecs R. Sco Hacker Jönköping Inernaional Business School

More information

A Hop Constrained Min-Sum Arborescence with Outage Costs

A Hop Constrained Min-Sum Arborescence with Outage Costs A Hop Consrained Min-Sum Arborescence wih Ouage Coss Rakesh Kawara Minnesoa Sae Universiy, Mankao, MN 56001 Email: Kawara@mnsu.edu Absrac The hop consrained min-sum arborescence wih ouage coss problem

More information

Block Diagram of a DCS in 411

Block Diagram of a DCS in 411 Informaion source Forma A/D From oher sources Pulse modu. Muliplex Bandpass modu. X M h: channel impulse response m i g i s i Digial inpu Digial oupu iming and synchronizaion Digial baseband/ bandpass

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION DOI: 0.038/NCLIMATE893 Temporal resoluion and DICE * Supplemenal Informaion Alex L. Maren and Sephen C. Newbold Naional Cener for Environmenal Economics, US Environmenal Proecion

More information

This document was generated at 7:34 PM, 07/27/09 Copyright 2009 Richard T. Woodward

This document was generated at 7:34 PM, 07/27/09 Copyright 2009 Richard T. Woodward his documen was generaed a 7:34 PM, 07/27/09 Copyrigh 2009 Richard. Woodward 15. Bang-bang and mos rapid approach problems AGEC 637 - Summer 2009 here are some problems for which he opimal pah does no

More information

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

More information

Stochastic models and their distributions

Stochastic models and their distributions Sochasic models and heir disribuions Couning cusomers Suppose ha n cusomers arrive a a grocery a imes, say T 1,, T n, each of which akes any real number in he inerval (, ) equally likely The values T 1,,

More information

Reliability of Technical Systems

Reliability of Technical Systems eliabiliy of Technical Sysems Main Topics Inroducion, Key erms, framing he problem eliabiliy parameers: Failure ae, Failure Probabiliy, Availabiliy, ec. Some imporan reliabiliy disribuions Componen reliabiliy

More information

control properties under both Gaussian and burst noise conditions. In the ~isappointing in comparison with convolutional code systems designed

control properties under both Gaussian and burst noise conditions. In the ~isappointing in comparison with convolutional code systems designed 535 SOFT-DECSON THRESHOLD DECODNG OF CONVOLUTONAL CODES R.M.F. Goodman*, B.Sc., Ph.D. W.H. Ng*, M.S.E.E. Sunnnary Exising majoriy-decision hreshold decoders have so far been limied o his paper a new mehod

More information

Macroeconomic Theory Ph.D. Qualifying Examination Fall 2005 ANSWER EACH PART IN A SEPARATE BLUE BOOK. PART ONE: ANSWER IN BOOK 1 WEIGHT 1/3

Macroeconomic Theory Ph.D. Qualifying Examination Fall 2005 ANSWER EACH PART IN A SEPARATE BLUE BOOK. PART ONE: ANSWER IN BOOK 1 WEIGHT 1/3 Macroeconomic Theory Ph.D. Qualifying Examinaion Fall 2005 Comprehensive Examinaion UCLA Dep. of Economics You have 4 hours o complee he exam. There are hree pars o he exam. Answer all pars. Each par has

More information

MATHEMATICAL DESCRIPTION OF THEORETICAL METHODS OF RESERVE ECONOMY OF CONSIGNMENT STORES

MATHEMATICAL DESCRIPTION OF THEORETICAL METHODS OF RESERVE ECONOMY OF CONSIGNMENT STORES MAHEMAICAL DESCIPION OF HEOEICAL MEHODS OF ESEVE ECONOMY OF CONSIGNMEN SOES Péer elek, József Cselényi, György Demeer Universiy of Miskolc, Deparmen of Maerials Handling and Logisics Absrac: Opimizaion

More information

Expert Advice for Amateurs

Expert Advice for Amateurs Exper Advice for Amaeurs Ernes K. Lai Online Appendix - Exisence of Equilibria The analysis in his secion is performed under more general payoff funcions. Wihou aking an explici form, he payoffs of he

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information