Demand Side Management in Smart Grids using a Repeated Game Framework

Size: px
Start display at page:

Download "Demand Side Management in Smart Grids using a Repeated Game Framework"

Transcription

1 Demand Sde Managemen n Smar Grds usng a Repeaed Game Framework Lnq Song, Yuanzhang Xao and Mhaela van der Schaar, Fellow, IEEE Absrac Demand sde managemen (DSM) s a key soluon for reducng he peak-me power consumpon n smar grds. To provde ncenves for consumers o shf her consumpon o off-peak mes, he uly company charges consumers dfferenal prcng for usng power a dfferen mes of he day. Consumers ake no accoun hese dfferenal prces when decdng when and how much power o consume daly. Imporanly, whle consumers enjoy lower bllng coss when shfng her power usage o off-peak mes, hey also ncur dscomfor coss due o he alerng of her power consumpon paerns. Exsng works propose saonary sraeges for he myopc consumers o mnmze her shor-erm bllng and dscomfor coss. In conras, we model he neracon emergng among self-neresed, foresghed consumers as a repeaed energy schedulng game and prove ha he saonary sraeges are subopmal n erms of long-erm oal bllng and dscomfor coss. Subsequenly, we propose a novel framework for deermnng opmal nonsaonary DSM sraeges, n whch consumers can choose dfferen daly power consumpon paerns dependng on her preferences, rounes, and needs. As a drec consequence of he nonsaonary DSM polcy, dfferen subses of consumers are allowed o use power n peak mes a a low prce. The subse of consumers ha are seleced daly o have her jon dscomfor and bllng coss mnmzed s deermned based on he consumers power consumpon preferences as well as on he pas hsory of whch consumers have shfed her usage prevously. Imporanly, we show ha he proposed sraeges are ncenve-compable. Smulaons confrm ha, gven he same peak-o-average rao, he proposed sraegy can reduce he oal cos (bllng and dscomfor coss) by up o 50% compared o exsng DSM sraeges. Keywords Smar Grds; Demand Sde Managemen; Crcal Peak Prcng; Consumer Dscomfor; Repeaed Games; Incenve Desgn. I. ITRODUCTIO Smar grds am o provde a more relable, envronmenally frendly and economcally effcen power sysem [1][2]. The uly company sells elecrcy o consumers, who are equpped wh smar meers. Smar meers exchange nformaon beween consumers and he uly company, and schedule he household energy consumpon for consumers. The nformaon gahered hrough smar meers can be used by he uly company o adjus he elecrcy prces. Demand Sde Managemen (DSM) s a key mechansm o make smar grds more effcen and cos-effecve [1][2]. DSM refers o he programs adoped by uly companes o drecly or L. Song, Y. Xao and M. van der Schaar are wh Deparmen of Elecrcal Engneerng, UCLA. Emal: songlnq@ucla.edu, {yxao, mhaela}@ee.ucla.edu.

2 2 ndrecly nfluence he consumers power consumpon behavor n order o reduce he Peak-o- Average Rao (PAR) of he oal load n he smar grd sysem. A hgher PAR resuls n much hgher operaon coss and possbly ouages of he sysem. DSM ams o ncenvze consumers o shf her peak-me power consumpon o off-peak mes, hereby resulng n sgnfcan PAR reducons n he power sysem. A. Relaed Work Drec Load Conrol (DLC) and Smar Prcng (SP) are wo popular exsng approaches for mplemenng DSM. DLC refers o he program n whch he uly company can remoely manage a fracon of consumers applances o shf her peak-me power usage o off-peak mes [3]. Alernavely, SP [6]-[18] provdes an economc ncenve for consumers o volunarly manage her power usage. Examples are Real-Tme Prcng (RTP) [6], Tme-Of-Use Prcng (TOU) [14], Crcal Peak Prcng (CPP) [16]-[18], ec. However, he above works [3][6]-[8][16]- [18] do no consder he consumers dscomfor coss whch s nduced by alerng her power consumpon paerns. Some recen works consdered consumers dscomfor coss [4][5][9]-[15][20][23] and amed o jonly mnmze he consumers bllng and dscomfor coss (referred o subsequenly as he oal cos). These works can be classfed no wo caegores, dependng on he deployed consumer model. The frs caegory assumed ha he consumers are prce-akng (.e., hey do no consder how her consumpon wll affec he prces). By assumng ha consumers are prceakers, he decson makng of a sngle foresghed consumer s formulaed as a sochasc conrol problem amng o mnmze s long-erm oal cos n [9]-[11]. Alernavely, n [12][13], mulple myopc consumers am o mnmze her curren oal coss and her decsons are formulaed as sac opmzaon problems among cooperave users for whch dsrbued algorhms are proposed o fnd he opmal prces. The second caegory assumed ha he consumers are myopc and prce-ancpang (.e., hey ake no accoun how her consumpon wll affec he prces). In hs case, each consumer s power usage affecs he oher consumers bllng coss. These works [6]-[8][15] modeled he neracons emergng among myopc consumers as one-sho games and suded he ash

3 3 DLC or SP Prce-akng or prce-ancpang consumers TABLE I. COMPARISO WITH EXISTIG WORKS Myopc or foresghed consumers Sngle consumer or mulple consumers Model Consumer Dscomfor [3] DLC - Myopc Mulple Opmzaon o [4][5] DLC - Myopc Mulple Opmzaon Yes [6]-[8] SP Prce-ancpang Myopc Mulple One-sho game o [9]-[11] SP Prce-akng Foresghed Sngle Sochasc conrol Yes [12][13] SP Prce-akng Myopc Mulple Opmzaon Yes [14] SP Prce-akng Myopc Sngle One-sho game Yes [15] SP Prce-ancpang Myopc Mulple One-sho game Yes [16] SP Prce-akng Myopc Mulple Opmzaon o Our Work SP Prce-ancpang Foresghed Mulple Repeaed game Yes equlbrum (E) of he emergng game. In hs paper, we also model he consumers as prceancpang. However, n our model, he consumers nerac wh each oher repeaedly and are foresghed, hereby engagng n a repeaed game. I s well-known ha n one-sho games wh myopc prce-ancpang consumers, he sysem performance (.e., he oal cos) a he equlbrum can be much worse han he opmal performance acheved by myopc prce-akng consumers [25]. One way o acheve he same performance for myopc prce-ancpang consumers as n he case of myopc prce-akng consumers s o model he sysem as a repeaed game, and o ncenvze he prce-ancpang consumers o choose (n equlbrum) he opmal consumpon adoped by myopc prce-akng consumers [24]. Imporanly, n hs paper we go one sep furher, and adop a novel repeaed game framework whch can sgnfcanly ouperform he opmal performance obaned by he myopc prce-akng consumers [12][13]. All he exsng works consderng mulple consumers [3]-[8][12]-[16] assumed ha he consumers are myopc and mnmze her curren coss. The opmal DSM sraeges n hese works are saonary,.e., all consumers adop fxed daly/weekly power consumpon paerns as long as he sysem parameers (e.g., he consumers desred power consumpon paerns) do no change. However, as we wll show laer n he paper, he saonary DSM sraeges are subopmal n erms of he long-erm oal cos. To mnmze he oal cos, some consumers are requred o shf her peak-me power usage o he off-peak mes whle he remanng consumers can use energy when desred. By deployng hs opmal sraegy, he consumers who shf her peak-me consumpon ncur dscomfor coss, bu hs leads o a reducon of he peak-me prce and of he bllng cos of all he consumers. Imporanly, our proposed nonsaonary DSM

4 4 sraegy can acheve he opmal oal cos whle ensurng farness among consumers by recommendng dfferen subses of consumers (referred o as he acve se) o shf her peakme consumpon each day. The acve se s deermned by he consumers preferences and he pas selecon of acve ses. A dealed comparson of our work and exsng works s hghlghed n Table I. B. Our Conrbuons In hs paper, snce he consumers say n he sysem for a long me and nerac wh each oher repeaedly, we formulae he consumers neracons as a repeaed game. The repeaed naure of he neracon provdes ncenves for prce-ancpang consumers o cooperae (as shown n [24]). Alhough he proposed framework can mprove he performance of he saonary DSM sraeges dscussed for any SP scheme, we focus on he CPP scheme, whch has been wdely used for resdenal consumers and s shown o work well n praccal scenaros [17]-[19]. CPP defnes peak days n a year or peak mes n a day, and charges hgher prces durng hese peak hours f CPP evens, such as sysem load warnng, exreme weaher condons, and sysem emergences, occur [18]. Based on he repeaed game model, we propose an opmal nonsaonary 1 DSM mechansm ha mnmzes he oal cos and ouperforms he opmal saonary DSM sraegy. In addon, he proposed sraegy s Incenve-Compable (IC), namely he self-neresed consumers wll fnd n her self-neres o follow he recommended sraegy. Each day, he DSM sraegy selecs an acve se of consumers based on her preferences of wheher or when o shf and on he pas hsory of consumpon paern shfs. These consumers sacrfce her curren dscomfor coss o mnmze he oal bllng cos. In reurn, hey wll enjoy n he fuure lower bllng coss whou ncurrng dscomfor coss when oher consumers are chosen n he acve se. In hs way, he proposed sraegy mnmzes he long-erm oal cos whle ensurng farness among he consumers. In summary, he man conrbuons of our work are as follows: 1 Recall ha n a saonary DSM mechansm, he consumers adop fxed daly/weekly power consumpon paerns as long as he sysem parameers do no change. In conras, n a nonsaonary DSM mechansm, he consumers may adop dfferen daly/weekly power consumpon paerns (e.g., a consumer may shf s peak-me consumpon oday bu no shf omorrow) even f he sysem parameers reman he same.

5 5 A Repeaed Game Framework: A repeaed game framework s proposed o model he neracons among foresghed prce-ancpang consumers over me. Jon Bllng and Dscomfor Coss Mnmzaon: The proposed DSM mechansm consders no only he consumers bllng coss bu also her dscomfor coss. Opmal onsaonary DSM mechansm: We analycally prove ha n he energy schedulng game wh dscomfor coss, nonsaonary DSM mechansms can ouperform saonary mechansms, and propose an opmal nonsaonary DSM mechansm ha can be easly mplemened. IC sraeges: The DSM mechansm s IC, meanng ha he self-neresed consumers have no ncenve o devae from he recommended sraegy. Consumer heerogeney: Our framework can model dfferen ypes of consumers wh dfferen dscomfor coss. Moreover, dfferen consumers may have dfferen preferences on how and when her consumpon paerns are shfed. We sudy he mpac of dfferen ypes of consumers on he performance of he proposed mechansm. The res of hs paper s organzed as follows. Secon II models he repeaed energy schedulng game and formulaes he DSM mechansm desgn problem. Secon III formally nroduces he proposed algorhm for consrucng he DSM sraegy and dscusses how o mplemen he proposed algorhm n he smar grd sysem. Secon IV provdes smulaon resuls o valdae he performance of he proposed algorhm. Secon V concludes he paper. II. SYSTEM MODEL A. Energy Schedulng Game A smar grd sysem consss of a uly company and mulple consumers, as shown n Fg.1. The DSM s mplemened hrough smar meers on he consumer sde and a DSM cener n he uly company [2]. We denoe he se of consumers by {1,2,..., }. Tme s dvded no perods 0,1,2, L. We assume ha each perod s dvded no H me slos wh equal lengh and denoe he se of me slos by H {1,2,..., H}. oe ha we use perod here o denoe

6 6 Informaon Flow Power Flow Smar Meer da (% ), a da (% ), a 1 1 a1 p l p,, 1 Uly Company a Smar Meer Consumer 1... Fg. 1. Smar Grd Sysem Model. Consumer each sage of he neracon among consumers and use me slo o denoe he dscree me o schedule power usage whn a perod. In hs paper, we consder a perod o be one day as n [6]- [9][12]-[16], and each slo can be one or mulple hours (e.g., H 24, H 12 ). For each consumer n perod, s acon s he power consumpon paern n ha perod, whch s he vecor of power consumpon a each me slo and denoed by a ( a,, a ), wh a A beng he power consumpon and A beng he se of power,,1 H, h consumpon a each me slo. The power consumpon a h, a me slo h consss of non-shfable load and shfable load. The non-shfable load, such as lghng, cookng, wachng TV, s no conrollable by smar meers, whle he shfable load, such as dsh and clohes washng, heang and coolng sysems, can be conrolled by he smar meers [6]-[8]. We denoe he non-shfable and shfable loads a me slo h H by sasfes a h,, h, h b h, 0 and b s and a s, 0, respecvely. Thus, consumer s power consumpon h b h, h,. oe ha here he shfable load s no only lmed o he deferrable load, bu also ncludes smar hermosas amed a conrollng heang/coolng sysems, whch can shf he energy usage ahead of me and have been shown [19] o sgnfcanly reduce he peak load n pracce, ec. We denoe by H a A he daly oal power consumpon for resdenal consumer, h1 h, where A s eher a consan or slowly varyng as n [6]-[8][13]-[16]. We denoe by a ( a, a, La 1 2 ) A he power consumpon profle of all consumers, where 1 A A. H

7 7 The oal load a me slo h, denoed by l a h, 1 h, s he sum of all consumers power consumpon. We assume ha he desred power consumpon paern n each perod for consumer s a [ a, a, La ] A, whch refers o s preferred daly power consumpon paern [22]. H,1,2, H The correspondng desred power consumpon profle of all consumers and oal load are denoed by a [ a, a, La ] A and 1 2 l a h 1 h,. Defne h argmax lh hh as he peak me of he day, he lengh of whch changes accordng o how H s se. Emprcal sudes show ha, compared wh ndusral and commercal consumers, resdenal consumers have very smlar peak-me shfable loads [19], mplyng haa b s, for each consumer. h, h, h Remark: Frs, dfferen consumers may have dfferen desred power consumpon paerns a. For example, mos consumers may have her hghes power consumpon n he evenng, whle some consumers may have hghes power consumpon n he mornng. Second, each consumer can have s ndvdual peak consumpon n hours oher han h ; h s deermned based on he aggregae power consumpon of all consumers, nsead of he power consumpon of ndvdual consumers. Thrd, n he model currenly presened, we le he desred daly power consumpon paern a o be he same for each day. One may argue ha he desred daly power consumpon should be dfferen n weekdays and n weekends. In hs case, we can easly exend he model such ha one perod s one week and a s he desred weekly power consumpon paern, where he desred daly power consumpon paerns are dfferen each day. The cos of a consumer consss of bllng and dscomfor coss. The bllng cos s he power consumed mulpled by he un prce and he dscomfor cos s he consumer s dscomfor caused by reschedulng s daly power consumpon from s desred power consumpon paern. We denoe c : A a R o be he cos funcon of consumer : H h1 h h, c( a ) p ( a ) a d ( a ), (1)

8 8 where p : h A a R denoes he prce a me slo h, a denoes he daly bllng H p ( ) a h1 h h, H cos, d : A a R denoes he dscomfor cos of consumer. Whou loss of generaly, he oal cos s he sum, nsead of he weghed sum, of bllng and dscomfor coss, because any wegh he consumer pus on he dscomfor cos can be absorbed n he expresson of he dscomfor cos funcon d. We wll dscuss laer he prce p and dscomfor cos d n deal. h To model he neracons among consumers, we frs formalze he one-sho energy schedulng game G {, { A },{ c } }, where, 1 1 { A } and 1 { } denoe he se of consumers, c 1 ses of acons and cos funcons for each consumer, respecvely. ex, we formalze he repeaed game model. In each perod, consumer deermnes s power consumpon paern a based on s hsory, whch s a collecon of all s pas power consumpon paerns and he pas prces made publc o he consumers. The hsory of consumer up o perod can be wren as: he nal hsory s defned as { a, p, a, p,..., a, p } ) H H ( A P for 0 and 0. The correspondng publc hsory s defned as { p, p,..., p } ( P ) for 0 and he nal hsory s defned as H 0. The publc sraegy of consumer s defned as a mappng from publc hsory o curren acon, denoed by H : U( P ) 0 H 0 a A, where ( P ) [24]. Due o realzaon equvalence prncple [24, Lemma 7.1.2], we lose nohng, n erms of he achevable operang pons, by resrcng o publc sraeges, compared o sraeges usng he enre hsory. Gven he sraegy profle of all consumers, denoed by (,, L ), consumer s 1 2 average long-erm cos s dscouned by a facor : C ( ) (1 ) c ( ( )) 0, (2) where c ( ( )) s he cos of consumer n perod. The dscoun facor [0,1) represens how consumers dscoun her moneary coss n he fuure, e.g. due o neres raes/nflaon. A smaller means ha he consumers dscoun her fuure coss more. Snce he neres

9 9 rae/nflaon s usually he same for all he consumers, we assume o be equal for all consumers as n [24][25]. The correspondng long-erm dscouned dscomfor cos s denoed by D ( ). Hence, he nfnely repeaed game model can be wren as RG H {, U ( ),{ },{ C ( )} } P, where, H ( ) 0 U P, { } and 1 { C ( )} 1 denoe he se of consumers, se of publc hsores, ses of sraeges and ses of cos funcons, respecvely. B. Crcal Peak Prcng Scheme Recall ha n Crcal Peak Prcng (CPP) scheme, he uly company charges a hgher prce n he crcal peak me when CPP evens occur (such as sysem load warnng, exreme weaher condons, and sysem emergences), amng a reducng he peak-me load of he sysem [16]- [18]. The maxmum number of crcal days and of crcal hours whn a day s predefned. We model he CPP prcng scheme wh a sngle crcal peak me and only consder he CPP evens rggered by he oal load n he sysem. The me-varyng prce funcon p ( ) h a s defned as: where p H p Lo p, 0 l l p ( ) p ( l ) h h h p l l Lo h h, H h h a, (3) are he peak prce and off-peak prce of he prcng model and l h s he hreshold of he oal load. When l h l h, he hgher prce wll no be rggered and p wll be adoped. Lo When l h l, he CPP even occurs and he hgher prce h H p wll be adoped. To provde ncenves for consumers o shf her peak-me loads, he hreshold l h s se o be below he peak load and above he off-peak load, namely, l l, l l, h H \{ h}. (4) h h h h Gven peak load reducon goal, namely h l, we furher se m l l / s h h h, where m s he smalles number of consumers needed o shf her peak-me consumpon such ha he peak-me prce s low. We denoe he prces whn a day by p ( p, p, Lp ). 1 2 H

10 10 C. Consumer Dscomfor Cos We use a dscomfor cos funcon o model he consumers dscomfor from reschedulng her power consumpon paerns. Many papers defne he dscomfor cos funcon based on he mposed change n her consumpon paerns,.e., he dsance beween consumer s desred demand and acual consumpon [14][20]-[23]. As n [14][20], we use a lnear weghed funcon o model he dscomfor cos: where k h,, H k ( a a ), a a h h, h, h, ( ) 1 0, a a d a, (5) R are parameers of he dscomfor cos funcon. Consumer s dscomfor cos d ( a ) s an ncreasng funcon of he dsance beween he rescheduled and he desred power consumpon paerns. The above dscomfor cos funcon s able o capure several mporan consumer preferences n erms of consumpon paern shfng. Frs, capures he consumer s wllngness o shf: a consumer wh a large s less wllng o shf. Second, k h, capures he consumer s preference on how o shf. A larger k, h ndcaes ha he consumer does no wan o reduce s peak-me consumpon; a larger k h, for h h ndcaes ha he consumer does no wan o shf s peakme consumpon o me slo h. We are neresed n wo coss ha consumer can acheve: consumer s mnmum cos achevable by any power consumpon profle, denoed by c =mn c ( a )= p A, and consumer aa Lo s mnmum cos achevable by any power consumpon profle n whch consumer shfs all s peak-me shfable load, denoed by c% mn c ( a ) p A d ( a% ) a A, a Lo h, bh,, where a% argmn d ( a ). Clearly, we always have c c%, because c% s he mnmum cos H a A, a b h, h, achevable when consraned o consumer s ceran power consumpon paerns (whch cause dscomfor coss). The cos c s mporan because s he mnmum cos ha consumer could

11 11 possbly acheve. The cos c% s also mporan because represens he mnmum cos ha consumer could possbly acheve when shfs all of s peak-me consumpon (.e. when s seleced o be n he acve se). Based on he relaonshp beween bllng and dscomfor coss, we can classfy he consumers no hree classes. Consumers wh low dscomfor coss: Ths class ncludes consumers whose dscomfor coss are low compared wh ph p Lo. In hs case, he ncrease n he prce maers much more o consumers han he dscomfor coss. For hese consumers, he dscomfor coss can be reasonably gnored and only he bllng coss need o be consdered. Consumers wh medum dscomfor coss: Ths class ncludes consumers who have medum dscomfor coss wh respec o ph p Lo. In hs case, boh bllng and dscomfor coss need o be consdered. We say a consumer has medum dscomfor cos f ( p p ) l H Lo h m d ( a% ) and [( p p )] a. (6) H Lo h, The frs nequaly mples ha he dscomfor cos s no oo large, such ha he consumers are wllng o shf her peak-me consumpon f by dong so her bllng coss can be grealy reduced. The second nequaly mples ha he dscomfor cos s no oo small, such ha each consumer does care abou s own dscomfor cos and s no wllng o shf s peak-me consumpon every day. Consumers wh hgh dscomfor coss: Ths class ncludes consumers who have hgh dscomfor coss parameers compared wh ph p Lo. In hs case, dscomfor coss are hgh when alerng her power consumpon paerns and hence hey should be recommended no o change her power consumpon paerns. In a praccal sysem, he above hree classes of consumers coexs. However, snce he recommendaons o consumers wh low and hgh dscomfor coss can be deermned based on her cos characerscs (namely, consumers wh low dscomfor coss always shf and choose a%, whle consumers wh hgh dscomfor coss canno shf her power consumpon paern and

12 12 choose a ), we only need o consder he DSM sraegy for consumers wh medum dscomfor coss. Hence, n Secon III, we focus on dervng he opmal DSM sraegy for consumers wh medum dscomfor coss. The opmal DSM sraegy derved for consumers wh medum dscomfor coss, along wh he predeermned recommended sraeges for consumers wh hgh and low dscomfor coss, consues he DSM sraegy used n he praccal sysem wh hree classes of consumers. D. Problem Formulaon In hs subsecon, we consder he sysem wh self-neresed consumers and formulae he opmal IC DSM mechansm desgn problem. The desgner s he benevolen uly company ha ams a mnmzng he oal cos (maxmzng he socal welfare) n he smar grd sysem. However, excep for he oal cos, mananng farness among all he consumers s essenal [19]. Hence, he mechansm wll ensure ha he average dscomfor cos of consumer s no greaer han a maxmal value Therefore, he DSM mechansm Desgn Problem (DDP) can be formulaed as π h, h, H (DDP): mnmze C ( ) subjec o a b, T A a, T h1 h, D( ) D, s IC,max. D,max. The uly company wll solve hs problem, hen recommend he consumers wh he opmal soluon. The abovemenoned energy schedulng game model can be exended or revsed n dfferen ways o accommodae varous sysems. Frs, he CPP scheme may be only appled n he summer and wner seasons, when he sysem load s hgh. In hs case, he DSM mechansm can updae he sysem parameers and begn o solve he DDP problem agan, a he begnnng of a new season. Second, he daly power consumpon A for a specfc consumer can vary for dfferen days durng a week, due o her dfferen weekly rounes, e.g. hey do he laundry on Frdays ec.

13 13 In hs case, we can se he perod n he model o be a week, nsead of a day, and se he daly oal power consumpon and desred power consumpon paern as and H mod7 a A h1 h, a a mod7, where { } 6 A and 6 { a } 0 0 are he daly power consumpon and desred power consumpon paerns of a week and are dfferen for dfferen days of a week. III. OPTIMAL STRATEGIES In hs secon, we solve he DDP problem defned n he prevous secon. We frs dscuss a benchmark case - he performance of he one-sho energy schedulng game. Subsequenly, we characerze he Pareo effcen operang pons of he repeaed energy schedulng game and propose our nonsaonary algorhm ha acheves he opmal soluon of he DDP problem. Fnally, we descrbe he mplemenaon of he proposed DSM n he smar grd sysem. A. Benchmark The One-sho Energy Schedulng Game In he unque E of he one-sho energy schedulng game, he consumer chooses s desred power consumpon paern. We sae hs formally n he followng heorem. Theorem 1 (ash Equlbrum of he One-sho Game): The one-sho energy schedulng game has a unque E, n whch each consumer chooses s desred power usage as a a,. (7) Proof: The dea of he proof s o show ha a s he domnan sraegy for consumer. The complee proof s gven n Appendx A. The resul of Theorem 1 s he well-known ragedy of commons. Each consumer ams o mnmze s ndvdual bllng and dscomfor coss and hus, wll myopcally fnd n s selfneres o mnmze s ndvdual cos by sckng o s desred power consumpon paern. Ths resuls n a hgh prce and low socal welfare. Accordng o Theorem 1, he oal cos a he unque E s C 1 E,, where C c E ( a ). In he followng secon we wll quanfy he neffcency, of he E and propose a novel DSM mechansm ha can acheve he opmal socal welfare.

14 14 B. Pareo-opmal Regon of he Repeaed Energy Schedulng Game We formally characerze he achevable operang pons of he repeaed energy schedulng game. I s well known ha as long as s suffcenly close o 1, he achevable regon of he repeaed energy schedulng game s he convex hull of he one-sho energy schedulng game [24, Lemma 3.7.1]. Hence, we can wre he achevable regon of repeaed energy schedulng game as H C conv{ ca ( ) a A, a A, a b }, where conv{ X } s he convex hull of X and c c L c 1 2 h1 h, h, h, ca ( ) ( ( a), ( a), ( a )) s he cos profle of he consumers. We wll prove ha he Pareoopmal regon of he repeaed energy schedulng game (.e., he Pareo boundary of he se of achevable cos profles), denoed by B, s a face of C ha s a polyope wh dmenson 1. Ths means ha we can analycally express he Pareo-opmal regon. Theorem 2 formally characerzes he Pareo-opmal regon B analycally. Theorem 2: The Pareo-opmal regon of he repeaed energy schedulng game s B { C ( C, C, L, C ) ( C c %)/( c c %) mc, c %}. (8) In addon, he saonary DSM sraeges can only acheve he exreme pons 2 of B. Proof: See Appendx B. Theorem 2 does no only characerze he Pareo-opmal regon (.e., par of a hyperplane) of he repeaed energy schedulng game, bu also proves ha by usng saonary DSM sraeges, we canno acheve any pons on he Pareo-opmal regon oher han he exreme pons. oe, however, ha, he exreme pons can only be acheved by he acon profles n whch m consumers shf her peak consumpon and ncur cos c, whle he oher consumers do no shf and ncur cos c%. I s clear ha he exreme pons are no desrable operang pons, because an exreme pon can be acheved only when a fxed se of m consumers shf her peak-me consumpon all he me, whch s unfar for hese m consumers because hey ncur hgh dscomfor coss all he me. Hence, he desred operang pons le n he neror of he Pareoopmal regon, whch can be acheved by nonsaonary sraeges n he repeaed energy schedulng game accordng o Theorem 2. 2 An exreme pon of a convex se s he pon ha s no he convex combnaon of any oher pons n hs se. In our case, snce B s par of a hyperplane, he exreme pons wll be he verces of B.

15 15 By addng IC consrans and he consrans on he maxmum dscomfor coss, he feasble Pareo-opmal regon can be wren as B C { C ( C, C, L, C ) ( C c% )/( c c% ) mc, c%, C C }, (9) mn{, }, C c% D.,max,max where C C C,max E, C. onsaonary DSM Mechansm Gven he Pareo-opmal regon, we can hen reformulae he DDP problem as a lnear programmng problem: mnmze subjec o C C B. (10) The soluon of (10) s an exreme pon of B, denoed by C C C C (,, L ). C 1 2 Theorem 3: The feasble Pareo-opmal regon B s achevable f he dscoun facor C sasfes C Proof: See Appendx C (11) m 1 Gven a desred operang pon n B C, we can use he onsaonary DSM (-DSM) algorhm, descrbed n Table II, o consruc he DSM sraegy. In perod, he -DSM algorhm chooses he acve se I () I conssng of m ou of consumers o reschedule her power consumpon paerns, where I s he se of all possble ndex combnaon ha conanng m consumers ou of. The choce of whch m consumers are seleced depends on how far hey are from her arge cos and hs s measured by ndex g () and he m consumers who have he larges g () wll aler her power consumpon paerns. Theorem 4: When he dscoun facor sasfes (11), he -DSM algorhm s IC and can acheve he opmal operang pon C C C C (,, L ). Proof: See Appendx D. 1 2

16 16 Theorem 3 and 4 sae ha when he dscoun facor sasfes (11), he opmal nonsaonary DSM mechansm can be consruced by he -DSM algorhm. In pracce, he daly neres rae s low enough o make condon (11) be sasfed and he proposed DSM mechansm can accommodae a large number of consumers. D. DSM Implemenaon TABLE II. OSTATIOARY DSM (-DSM) ALGORITHM Inpu: Targe average cos vecorc ( C1, C2, LC), 0. Oupu: Opmal sraegy 1: Se g () ( C c% ) ( c c% ). j j j j j 2: Repea 3: If p ( ) p,, hen h H 4: Recommend acon a o all consumer. 5: Else 6: Fnd he acve se I () {,, L } of m consumers who have he m larges ndces gj (). 1 2 m 7: Recommend acon a% o I (), and a o I (). 8: Observe consumers acon a. 9: If all consumers follow he recommendaon, hen 10: Updae g( 1) [ g() ( 1) 1 { I()} ]/ for all. 11: Broadcas p () p for all h. 12: Else h 13: Broadcas p () ph and p () p, h h. 14: End f 15: End f 16: 1 17: End Repea h Lo h In hs subsecon, we descrbe how o mplemen he DSM mechansm n he smar grd sysem based on he proposed -DSM algorhm. The publc nformaon, such as prce nformaon and oal load, s known o every consumer and he prvae nformaon of consumer, such as dscomfor cos, power consumpon, s only known o s own. There are wo phases n mplemenng he DSM mechansm: he nalzaon phase and he run me phase. As n [2], he uly company recommends he sraegy one day n advance. 3 Consumers can schedule her power usage and exchange nformaon wh uly company hrough smar meers for he upcomng day so as o mnmze her bllng and dscomfor coss based on he recommendaon. Lo 3 The DSM can also be mplemened n real-me or T-day-ahead, as long as here s enough me for he consumers o respond o he recommendaons.

17 17 TABLE III. IFORMATIO EXCHAGE OF IMPLEMETATIO Informaon exchange Publc nformaon (known o all consumers) Consumer s prvae nformaon Uly company o consumer Consumer o uly company Inalzaon pp,,, ( ), H plo lh m l h d ( a ), a 1: p, p, l, l, C H Lo h h Run me 4:, 7: p Inalzaon 2: d ( a% )( or c, c% ), a Run me 56: a Fg.2. The Implemenaon of DSM. The deals of publc and prvae nformaon, he wo phases and he correspondng nformaon exchange of he wo phases are shown n Fg. 2 and Table III. Inalzaon: 1 The uly company broadcass he CPP parameers: p Lo, p and l H h. 2 Consumer calculaes s cos parameers c% and c, and sends hem o he uly company. 3 The uly company calculaes he arge cos vecor C ( C, C, L C ) by solvng DDP 1 2 problem, and hen sends he arge cos vecor o all consumers. Run me: In each perod: 4 The DSM cener calculaes he acve se I() and hen makes he recommendaon () o each consumer. 5 Each consumer chooses s energy schedulng sraegy. oe ha he consumers can devae from he recommended sraegy. However, snce he sraegy s IC, s n her self-neres o comply wh he recommendaon. 6 The smar meer repors he power consumpon paerns a o DSM cener. 7 The uly company calculaes he prce based on he oal load and hen broadcass he prces. Go o nex perod and repea 4-7. The consumers dscomfor cos funcons may vary over me. For example, hey may have dfferen desred power consumpon paerns n dfferen seasons/monhs, or hey may prefer o shf he peak-me consumpon o dfferen off-peak mes (.e., dfferen k h, ) due o her change of schedules (e.g., due o a new job), or hey smply change her wllngness (.e., dfferen ) o shf. To accommodae for such occasonal changes (e.g., once a monh), he uly company can

18 18 run he nalzaon phase agan n order o calculae he new opmal arge cos vecor, and hen run he -DSM algorhm based on he updaed arge cos vecor. IV. UMERICAL RESULTS In hs secon, we compare he performance of our proposed DSM mechansm wh hose obaned usng he one-sho energy schedulng games wh myopc prce-ancpang consumers [6]-[8][15], jon opmzaon wh myopc prce-akng consumers [12][13], as well as sochasc conrol mehods wh a sngle foresghed consumer [9]-[11]. Then we sudy he mpac of dfferen consumer preferences. In addon, we sudy he sysem performance when he percenage of shfable load and he lengh of he peak me vary. Throughou hs secon, we use he followng sysem parameers by defaul unless we change some of hem explcly. We consder he scenaro ha H 24 me slos and se he dscoun facor of he consumers o be 0.995, whch ensures ha (11) s sasfed. The prcng scheme ses he peak prce and off-peak prce o be p 0.8 $/kwh and p 0.1 $/kwh 4. Accordng H Lo o he uly company s PAR goal (namely he percenage of reducon n peak-me load, wren as l l / l h h h ), he hreshold l h n (3) wll be se o an approprae value o conrol he parameer m. We smulae boh he scenaro wh heerogeneous consumers wh parameers shown n Table IV and Fg. 3 and he scenaro wh homogeneous consumers wh he same parameers as Type 1 consumers descrbed n Table IV and Fg. 3. In hs expermen, he shfable load of each consumer s se o be 40% of he consumer s oal load. A. Comparson wh Exsng Mechansms In hs subsecon we compare our proposed -DSM algorhm wh he exsng ones [6][7]. TABLE V. COMPARISOS OF DIFFERET MECHAISMS TABLE IV. PARAMETERS OF THREE TYPES OF COSUMERS A (kwh) kh, ( h 1 o 14) / ($) D,max ($) kh, ( h 15 o 24) ($/kwh) Type / Type / Type / Works Algorhm Sraeges [6]-[8] H OG-DSM mn{ ( ), ( )} [15] p 1 h a a h h d a a [12][13] JO-DSM mn [ ph( a ) ah, d( a)] a H n1 h1 [9]-[11] SC-DSM mn{ pa, ( 1) d ( a )} Our work a H h1 h h -DSM mn{( 1) c ( ( ))} 0 4 Accordng o [17][19], he peak prce s ofen a leas 6 mes hgher han he off-peak prce.

19 19 Toal cos Performance gan Fg.3. The Desred Power Consumpon Paerns of Type 1, 2, 3 Consumers. TABLE VI. COMPARISO OF TOTAL COSTS ACHIEVED BY DIFFERET ALGORITHMS umber of consumers (Homogeneous, PAR<2.280) umber of consumers (Heerogeneous, PAR<2.359 ) OG-DSM JO-DSM SC-DSM DSM Over OG-DSM 38% 39% 39% 39% 40% 49% 49% 50% 50% 50% Over JO-DSM 27% 28% 28% 28% 28% 38% 39% 39% 40% 40% Over SC-DSM 34% 35% 35% 35% 35% 45% 46% 46% 47% 47% We dvde he exsng algorhms no hree caegores as shown n Table V. The One-sho Game based saonary DSM (OG-DSM) algorhms wh myopc prceancpang consumers [6]-[8][15] calculae he E and operae a E of he one-sho energy schedulng game, whch s characerzed n Theorem 1. The Jon Opmzaon (JO-DSM) algorhms wh myopc prce-akng consumers [12][13] assume ha he obeden consumers jonly mnmze he oal cos of he sysem. In hs case, he opmal performance of saonary DSM mechansm can be acheved by approprae prcng schemes. The Sngle-consumer Sochasc Conrol (SC-DSM) mehods [9]-[11] ry o use sochasc conrol mehods o mnmze he oal cos of a sngle consumer. In hs case, he uly company ses he prce as ph p H and p p, h h, and he consumer buys energy n advance accordng o s scheduled power h Lo consumpon paern a. We assume ha renewable energy s avalable wh probably 5 0.8, n whch case he consumer can reschedule s power consumpon paern o he desred paern whou sufferng he dscomfor cos snce he energy supply s abundan. The renewable energy s no avalable wh probably 1 0.2, n whch case he consumer mus comply wh s scheduled power consumpon paern and wll ncur dscomfor cos d ( a ). 5 Ths probably comes from he uncerany of renewable energy generaon (wheher s wndy n wnd energy generaon, wheher s shny n solar energy generaon, ec.).

20 20 Gven he same PAR goal, he comparson of oal coss usng hese four algorhms s shown n Table VI. We can see ha when he number of consumers ncreases, he -DSM algorhm sgnfcanly ouperforms oher hree algorhms. The cos reducons compared o OG-DSM, JO- DSM and SC-DSM are 40%, 28% and 35% n homogeneous case and 50%, 40% and 47% n heerogeneous case, respecvely. oe ha our algorhm, whch s IC, can sgnfcanly ouperform he JO-DSM algorhm, even hough s no IC. B. Impac of Dscomfor Cos TABLE VII. COMPARISO WITH BILLIG COST MIIMIZATIO ALGORITHM Toal Cos/ umber of consumers (Homogeneous) umber of consumers (Heerogeneous) Performance gan DSM Bllng cos mnmzaon Performance gan 35% 36% 36% 36% 36% 48% 48% 48% 48% 48% Some works [6]-[8][16] consder mnmzng he bllng cos only, whou akng no accoun he dscomfor cos; whle we consder he problem of jonly mnmzng bllng and dscomfor coss as n [9]-[15]. We compare he resuls of he bllng cos mnmzaon algorhm wh our proposed algorhm n Table VII. In he smulaon, we assume ha he PAR reducon goal s 10%. For he bllng cos mnmzaon algorhm, snce here are numerous opmal soluons whch acheve he mnmal bllng cos, we choose a far soluon where all consumers shf he same amoun of peak-me consumpon o off-peak mes, and hen calculae he bllng and dscomfor coss. By comparng he resuls n Table VII, we can see ha he performance of our proposed algorhm sgnfcanly ouperforms he bllng cos mnmzaon algorhm, wh around 36% cos reducon for he homogeneous scenaro and 48% cos reducon for he heerogeneous scenaro, respecvely. In fac, he bllng cos mnmzaon algorhm does no consder he mpac of consumers behavor on dscomfor coss, resulng n a hgher dscomfor cos han our proposed -DSM algorhm. The -DSM algorhm nduces he consumers o cooperae wh each oher o reduce her long-erm dscomfor coss and he consumers wh hgher dscomfor coss benef more hrough cooperaon n he heerogeneous scenaro. Thus, he performance gan of he -DSM algorhm over he bllng cos mnmzaon polcy n he heerogeneous scenaro s hgher han ha n he homogeneous scenaro.

21 21 C. Impac of Consumer Preferences Recall from Secon III ha n order o acheve he socal opmum,.e., he soluon o he DDP problem, he sysem runs he -DSM algorhm whch requres only a subse of he consumers o reduce her peak-me power usage. However, some consumers may have ceran preferences of when o shf her consumpon. Hence we smulae he scenaro of 3, m 1, o show he mpac of consumer preferences. The consumers ypes are Type 1, 2 and 3. We compare four cases: Case 1: all consumers can aler her power consumpon paerns. Case 2: Type 1 consumer canno aler s power consumpon paern on Mondays; Type 2 consumer canno aler s power consumpon paern on Tuesdays; Type 3 consumer canno aler s power consumpon paern on Wednesdays. Case 3: Type 1 and 3 consumers canno aler her power consumpon paerns on Mondays; Type 2 consumer canno aler s power consumpon paern on Tuesdays. Case 4: Type 1, 2 and 3 consumers canno aler her power consumpon paerns on Mondays. oe ha he four cases represen dfferen levels of heerogeney n consumer preferences, whch resul n dfferen levels of flexbly n he power consumpon schedulng. In Case 1, no consumer has specfc preferences on how o shf he power consumpon paern. Hence, he flexbly of schedulng s he hghes. In Case 2, each consumer has dsnc preferences on how o shf power consumpon. Snce her preferences are dfferen, he flexbly of schedulng s sll hgh. In Case 3, Type 1 and 3 consumers have smlar preferences (.e., boh canno shf on Fg.4. Comparson wh Dfferen Consumer Preferences.

22 22 Toal Cos JO-DSM SC-DSM -DSM TABLE VIII. IMPACT OF PERCETAGE OF SHIFTABLE LOAD PAR Percenage of shfable load (Homogeneous) Percenage of shfable load (Heerogeneous) Reducon goal 20% 30% 40% 50% 60% 20% 30% 40% 50% 60% 10.0% % % % % % % % % Mondays), whch makes he schedulng less flexble (.e. we can only aler Type 2 consumer s consumpon on Mondays). In Case 4, all consumers have he same preference (.e., all consumers canno shf on Mondays). Hence, he schedulng s he leas flexble. We can see from Fg. 4 ha he average dscouned daly cos ncreases from Case 1 o Case 4, as expeced. Ths mples ha he sysem performance depends on he heerogeney of consumer preferences. Regardng he mpac of he dscoun facor on he convergence rae, we can see ha here s no sgnfcan performance loss beween he case when 0.90 and he case when D. Impac of Percenage of Shfable Load In order o evaluae he sysem performance wh dfferen percenages of he shfable load, we show he oal cos wh dfferen PAR reducon goals, as shown n Table VIII. The OG-DSM canno mee he PAR reducon goals, so we compare he oher hree algorhms. We se he number of consumers o be 100. Smulaon resuls show ha gven he consran mposed by he PAR reducon goal, our algorhm can acheve around 8% and 9% reducon n he oal cos when he percenage of shfable load vares from 30% o 60%, n he homogeneous and heerogeneous scenaros, respecvely; whle hese reducons are 7% and 10% wh JO-DSM algorhm and 12% and 9% wh SC-DSM algorhm. Ths mples ha he performance of our algorhm s no sgnfcanly affeced by consumers beng homogeneous or heerogeneous, whle he performances of he oher wo algorhms sgnfcanly depend on he consumer heerogeney. A rade-off beween he opmal oal cos and he PAR can be observed: when a hgher hreshold l h (.e. a smaller m ) s chosen, fewer consumers are requred o change her power

23 23 TABLE IX. IMPACT OF THE LEGTH OF THE PEAK TIME Toal Cos umber of consumers (Homogeneous) umber of consumers (Heerogeneous) H= JO-DSM H= H= H= SC-DSM H= H= H= DSM H= H= consumpon paerns, resulng n a hgher PAR bu a lower dscomfor cos. As a resul, hs rade-off should be consdered when choosng he desgn parameer l h (or m ). E. Impac of he Lengh of he Peak me In hs subsecon, we evaluae he sysem performance n erms of he oal cos n Table IX, when he lengh of he peak me vares. We se he lengh of he peak me o be 1 hour ( H 24 ), 2 hours (H 12 ) and 4 hours (H 6 ). We need o change he desred power consumpon paern accordng o he lengh of he peak me (for example, he desred power consumpon level n me slo 1 when H 12 s he sum of he desred power consumpon levels n me slos 1 and 2 when H 24 ). We fx he PAR reducon goal a 10%. The OG- DSM canno mee he PAR reducon goal. We compare he oher hree algorhms. Frs, we can see ha our algorhm acheves a lower oal cos han he oher wo under all lenghs of he peak me wh boh homogeneous and heerogeneous consumers. Second, we can observe ha he oal cos rses when he lengh of he peak me ncreases. The ncrease n he oal cos comes from he ncreased dscomfor cos, because when he lengh of he peak me ncreases, he consumers need o shf more energy o off-peak mes n order o oban a low peak-me prce. Hence, s mporan o see how he performances of he algorhms vary wh he lengh of he peak me. Under our algorhm, here s a slgh ncrease n he oal cos,.e., 3% (homogeneous) and 2% (heerogeneous), respecvely, when he lengh of he peak me vares from 1 hour o 4 hours. In conras, he ncreases n he oal cos are sgnfcan usng he oher wo algorhms, namely 40% (homogeneous) and 38% (heerogeneous) for JO-DSM, and 50% (homogeneous) and 50% (heerogeneous) for SC-DSM. Clearly, our algorhm s less sensve o he lengh of he peak

24 24 me, yeldng hgh performance for dfferen lenghs of he peak me as compared o he oher wo algorhms, whose performances degrade a lo when he lengh of he peak me ncreases. V. COCLUSIOS In hs paper, we proposed a nonsaonary DSM mechansm, whch explos he repeaed neracons of he consumers over me. We rgorously prove ha he proposed DSM mechansm can acheve he socal opmum n erms of he oal cos, and ouperform exsng saonary DSM sraeges. Moreover, he proposed mechansm s IC, meanng ha each self-neresed consumer volunarly follows he power consumpon paerns recommended by he opmal DSM mechansm. Smulaon resuls valdae our analycal resuls on he DSM mechansm desgn and demonsrae up o 50% performance gans compared wh exsng mechansms, especally when here are a large number of heerogeneous consumers n he sysems. In addon, compared o he exsng mechansms, he performance of he proposed mechansm s much less sensve and hus much more robus o sysem parameers such as he consumer heerogeney and he lengh of he peak me. APPEDIX A. Proof of Theorem 1 Gven he oher consumers acon a, we calculae he bes response of consumer. We denoe l h a., j, j jh, (a) If l a l, he bes response of consumer s obvously a h, h, h, snce p A d ( a ) p A. Lo Lo (b) If l a l, hen h, h, h a A, a a, we have H p A ( p p ) a d ( a ) p A [ k ( a a ) ] Lo H Lo h, Lo h1 h, h, h, H h1 pa d ( a) h h, where he frs nequaly s due o (6). Therefore, he unque E s a a. H, (12)

25 25 B. Proof of Theorem 2 We frs prove ha -DSM can acheve he operang pons n B. Gven any subse I of m consumer, le I choose acon a% and I () choose acon a. Then s cos s C [ I] c% 1 c 1. I s easy o see ha he cos profle { I} { I} ( C [ I], C [ I],, C [ I]) 1 2 L s n B. Snce he oher cos profles n B are convex combnaons of {( C1[ I], C2[ I], L, C [ I])} I, all he cos profles n B can be acheved [24, Lemma 3.7.1]. ex, we prove B o be he Pareo-opmal regon by showng ha % % a, ( c ( a ) c )/( c c ) m and equaly holds only by choosng he acon profles 1 descrbed above, whch proves he second saemen of he heorem. To show hese, we analyze he soluon of he followng problem: ( C1, C2, L, C ) C mnmze ( C c% )/( c c% ) 1 (13) (a) Suppose p ( a ) h p H, hen he opmal acon s a and he correspondng opmal cos vecor s ca ( ) ( C, C, L, C ) due o Theorem 1. However, 1, E 2, E E, ( C c% )/( c c% ) ( p p ) a / d ( a% ) 1 E, 1 H Lo h, ( p p ) a /max{ d ( a% )} ( p p ) l /max{ d ( a% )} m H Lo 1 h, H Lo h where he las nequaly s due o (6).. (14) (b) Suppose p ( a ) p and he opmal acon s a. Obvously, h Lo 1 a h, l h, oherwse a leas one consumer can reduce s dscomfor cos by shfng peak-me power whle keepng he oal bllng cos and oher consumers dscomfor coss unchanged. ex, we show ha eher a a or a b. h, h, h, h, Suppose for and j, we have b a a and b a a. Whou loss of h, h, h, jh, jh, jh, generaly, we assume ( k mn k )/ d ( a% ) ( k mn k )/ d ( a% ). Then we compare he h, h, jh, jh, j j hh hh coss of and j usng acon a, where a a V, h, h, a argmn c ( aˆ ), h h h, a ˆA, aˆ a V h, h,,

26 26 a a V, jh, jh, a argmn c ( aˆ ), h h jh, a ˆA, aˆ jh, a V jh, j, wh V mn{ a b, a a } h, h, jh, jh,, and a a, n j,. Thus, we have: n n ( c ( a) c% ) ( c ( a) c% ) j j d ( a ) d ( a ) j j ( c c% ) ( c c% ) d ( a% ) d ( a% ) j j j j [ ( k mn k )( a a )]/ d ( a% ) [ ( k mn k )( a a )]/ d ( a% ) h, hh h, h, h, j jh, hh jh, jh, jh, j j.(15) [ ( k mn k )( a ( a ))]/ d ( a% ) [ ( k mn k )( a ( a ))]/ d ( a% ) h, hh h, h, h, j jh, hh jh, jh, jh, j j ( c ( a) c% ) ( c ( a) c% ) j j ( c c% ) ( c c% ) j j We noce ha eher a a% or a a,.e., eher c( a ) c or c ( a ) % c. By repeang j j j hs procedure, we can fnally ge a subse I, where acon a. Therefore, he soluon of problem (13) s n he form of I chooses acon a%, and I chooses C [ I] and he Pareo-opmal regon can be wren as (8). In addon, any oher acon profle wll resul n a % % ( c ( ) c )/( c c ) m. In oher words, he oher cos profles n he se B canno be 1 acheved by saonary sraeges. C. Proof of Theorem 3 Based on -DSM, n me perod, c %, c I () and c c, I (). Then, C (1 ) [ c 1 c 1 ] 0 { I()} { I()} where 1{} g s he ndcaor funcon. We denoe g g C c% g (,, L, g ), g 1 2 c c% % (16), g C c% c c% G { ( %)/ %), C B }. Then we noce ha and g g C c ( c c C C g (1 ) 0 { I ()} 1 and C B s equvalen o g G C C. We denoe by g () (1 ) 1 and g (0) g he connuaon cos a me { I ( )} and 0, respecvely. Then we use a backward nducon mehod o show ha connuaon cos a

27 27 any gven me can be decomposed of he curren cos and he connuaon cos a me 1. We call he vecor g() G C a feasble vecor and hrough he decomposon we wll show ha when he dscoun facor sasfes (11), he connuaon cos a any me 1 s also a feasble vecor, namely g( 1) show hs. C For me perod, suppose g() G C, hen G. Snce he orgnal g(0) G C, we use mahemacal nducon o Or equvalenly, We need o show g () (1 ) g ( 1) g ( 1) { I()} 1. (17) [ g () (1 )]/, I () g ()/, I () g ( 1) m and 0 g ( 1) g 1. (18). The former one s obvous due o hypohess. For he laer one, 0 g ( 1) g s obvous for I (), and for I (), we also need Ths can be smplfed as g () (1 ). (19) 0 g, I () 1 g () 1 g, I () 1 g () The frs erm s sasfed due o hypohess. The second erm requres max mn max{1 g ()} g () G I() I I() C Whou loss of generaly, we sor g () n an decreasng order, namely,. (20). (21) I s easy o calculae: g () g () L g () L g (). (22) 1 2 m max{1 g ()} 1 g (), (23) I() max{} I( ) and mnmax{1 g ()} (1 g ()) I() I I() m. (24)

28 28 Thus, he wors case g () n (21) s Thus we requre m1 m1 ( m g ) ( m g ) 1 1 g () ( g, L, g,, L, ). (25) 1 m1 ( m 1) ( m 1) m 1 Snce ( m g ) m ( m 1) 1 ( m ) 1 ( m 1) m1 g 1. (26), B s achevable, when (11) s sasfed. C D. Proof of Theorem 4 Obvously, C C C C (,, L ) s n B. Accordng o he proof of Theorem 3, C can be 1 2 C decomposed of curren cos plus connuaon cos: ( ) % ( %) () (1 )[ % ] ( ). (27) C c c c g c 1 c 1 C { I()} { I() } 1 When he dscoun facor sasfes (11), he opmal operang pon C can be acheved by decomposng no curren cos plus connuaon cos n each perod. We also noce ha he bes sraegy of I() n (21) s o choose I() {1,2, L m},.e., he ndexed wh he m larges g (). Therefore, he proposed -DSM can acheve opmal performance C C n he sysem. To show -DSM s IC, we need o show ha gven he operang pon ( C, C, LC 1 2 ) B C and an arbrary me perod, he connuaon cos s he mnmum cos achevable for consumer. By decomposon, we have ( ) (1 )[ % ] ( ). (28) C c1 c 1 C { I( )} { I( ) } 1 Obvously for consumer I (), has no ncenve o devae from sraegy For consumer c a (, ( )) I (), f devaes from sraegy n perod., he mnmum cos can acheve s, snce he bes response s o choose a accordng o Theorem 1. Accordng o one- sho devaon prncple [24], when consumer devaes, he oher consumers wll play E sraegy from me perod 1. So he cos from perod s C ( ) (1 ) c ( a, ( )) C ( ) C. (29) E 1 E,

Non-stationary Demand Side Management Method for Smart Grids

Non-stationary Demand Side Management Method for Smart Grids on-saonary Demand Sde Managemen Meod for Smar Grds Lnq Song Yuanzang Xao Maela van der Scaar Elecrcal Engneerng Deparmen UCLA Emal: songlnq@uclaedu ABSTRACT Demand sde managemen (DSM) s a key soluon for

More information

Graduate Macroeconomics 2 Problem set 5. - Solutions

Graduate Macroeconomics 2 Problem set 5. - Solutions Graduae Macroeconomcs 2 Problem se. - Soluons Queson 1 To answer hs queson we need he frms frs order condons and he equaon ha deermnes he number of frms n equlbrum. The frms frs order condons are: F K

More information

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION INTERNATIONAL TRADE T. J. KEHOE UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 27 EXAMINATION Please answer wo of he hree quesons. You can consul class noes, workng papers, and arcles whle you are workng on he

More information

Solution in semi infinite diffusion couples (error function analysis)

Solution in semi infinite diffusion couples (error function analysis) Soluon n sem nfne dffuson couples (error funcon analyss) Le us consder now he sem nfne dffuson couple of wo blocks wh concenraon of and I means ha, n a A- bnary sysem, s bondng beween wo blocks made of

More information

Political Economy of Institutions and Development: Problem Set 2 Due Date: Thursday, March 15, 2019.

Political Economy of Institutions and Development: Problem Set 2 Due Date: Thursday, March 15, 2019. Polcal Economy of Insuons and Developmen: 14.773 Problem Se 2 Due Dae: Thursday, March 15, 2019. Please answer Quesons 1, 2 and 3. Queson 1 Consder an nfne-horzon dynamc game beween wo groups, an ele and

More information

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005 Dynamc Team Decson Theory EECS 558 Proec Shruvandana Sharma and Davd Shuman December 0, 005 Oulne Inroducon o Team Decson Theory Decomposon of he Dynamc Team Decson Problem Equvalence of Sac and Dynamc

More information

Lecture 11 SVM cont

Lecture 11 SVM cont Lecure SVM con. 0 008 Wha we have done so far We have esalshed ha we wan o fnd a lnear decson oundary whose margn s he larges We know how o measure he margn of a lnear decson oundary Tha s: he mnmum geomerc

More information

Reactive Methods to Solve the Berth AllocationProblem with Stochastic Arrival and Handling Times

Reactive Methods to Solve the Berth AllocationProblem with Stochastic Arrival and Handling Times Reacve Mehods o Solve he Berh AllocaonProblem wh Sochasc Arrval and Handlng Tmes Nsh Umang* Mchel Berlare* * TRANSP-OR, Ecole Polyechnque Fédérale de Lausanne Frs Workshop on Large Scale Opmzaon November

More information

Online Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading

Online Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading Onlne Supplemen for Dynamc Mul-Technology Producon-Invenory Problem wh Emssons Tradng by We Zhang Zhongsheng Hua Yu Xa and Baofeng Huo Proof of Lemma For any ( qr ) Θ s easy o verfy ha he lnear programmng

More information

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS R&RATA # Vol.) 8, March FURTHER AALYSIS OF COFIDECE ITERVALS FOR LARGE CLIET/SERVER COMPUTER ETWORKS Vyacheslav Abramov School of Mahemacal Scences, Monash Unversy, Buldng 8, Level 4, Clayon Campus, Wellngon

More information

On One Analytic Method of. Constructing Program Controls

On One Analytic Method of. Constructing Program Controls Appled Mahemacal Scences, Vol. 9, 05, no. 8, 409-407 HIKARI Ld, www.m-hkar.com hp://dx.do.org/0.988/ams.05.54349 On One Analyc Mehod of Consrucng Program Conrols A. N. Kvko, S. V. Chsyakov and Yu. E. Balyna

More information

TSS = SST + SSE An orthogonal partition of the total SS

TSS = SST + SSE An orthogonal partition of the total SS ANOVA: Topc 4. Orhogonal conrass [ST&D p. 183] H 0 : µ 1 = µ =... = µ H 1 : The mean of a leas one reamen group s dfferen To es hs hypohess, a basc ANOVA allocaes he varaon among reamen means (SST) equally

More information

Variants of Pegasos. December 11, 2009

Variants of Pegasos. December 11, 2009 Inroducon Varans of Pegasos SooWoong Ryu bshboy@sanford.edu December, 009 Youngsoo Cho yc344@sanford.edu Developng a new SVM algorhm s ongong research opc. Among many exng SVM algorhms, we wll focus on

More information

CS286.2 Lecture 14: Quantum de Finetti Theorems II

CS286.2 Lecture 14: Quantum de Finetti Theorems II CS286.2 Lecure 14: Quanum de Fne Theorems II Scrbe: Mara Okounkova 1 Saemen of he heorem Recall he las saemen of he quanum de Fne heorem from he prevous lecure. Theorem 1 Quanum de Fne). Le ρ Dens C 2

More information

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!) i+1,q - [(! ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL The frs hng o es n wo-way ANOVA: Is here neracon? "No neracon" means: The man effecs model would f. Ths n urn means: In he neracon plo (wh A on he horzonal

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm H ( q, p, ) = q p L( q, q, ) H p = q H q = p H = L Equvalen o Lagrangan formalsm Smpler, bu

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm Hqp (,,) = qp Lqq (,,) H p = q H q = p H L = Equvalen o Lagrangan formalsm Smpler, bu wce as

More information

Robustness Experiments with Two Variance Components

Robustness Experiments with Two Variance Components Naonal Insue of Sandards and Technology (NIST) Informaon Technology Laboraory (ITL) Sascal Engneerng Dvson (SED) Robusness Expermens wh Two Varance Componens by Ana Ivelsse Avlés avles@ns.gov Conference

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Lnear Response Theory: The connecon beween QFT and expermens 3.1. Basc conceps and deas Q: ow do we measure he conducvy of a meal? A: we frs nroduce a weak elecrc feld E, and hen measure

More information

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4 CS434a/54a: Paern Recognon Prof. Olga Veksler Lecure 4 Oulne Normal Random Varable Properes Dscrmnan funcons Why Normal Random Varables? Analycally racable Works well when observaon comes form a corruped

More information

EEL 6266 Power System Operation and Control. Chapter 5 Unit Commitment

EEL 6266 Power System Operation and Control. Chapter 5 Unit Commitment EEL 6266 Power Sysem Operaon and Conrol Chaper 5 Un Commmen Dynamc programmng chef advanage over enumeraon schemes s he reducon n he dmensonaly of he problem n a src prory order scheme, here are only N

More information

( ) () we define the interaction representation by the unitary transformation () = ()

( ) () we define the interaction representation by the unitary transformation () = () Hgher Order Perurbaon Theory Mchael Fowler 3/7/6 The neracon Represenaon Recall ha n he frs par of hs course sequence, we dscussed he chrödnger and Hesenberg represenaons of quanum mechancs here n he chrödnger

More information

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany Herarchcal Markov Normal Mxure models wh Applcaons o Fnancal Asse Reurns Appendx: Proofs of Theorems and Condonal Poseror Dsrbuons John Geweke a and Gann Amsano b a Deparmens of Economcs and Sascs, Unversy

More information

Comb Filters. Comb Filters

Comb Filters. Comb Filters The smple flers dscussed so far are characered eher by a sngle passband and/or a sngle sopband There are applcaons where flers wh mulple passbands and sopbands are requred Thecomb fler s an example of

More information

Math 128b Project. Jude Yuen

Math 128b Project. Jude Yuen Mah 8b Proec Jude Yuen . Inroducon Le { Z } be a sequence of observed ndependen vecor varables. If he elemens of Z have a on normal dsrbuon hen { Z } has a mean vecor Z and a varancecovarance marx z. Geomercally

More information

Tight results for Next Fit and Worst Fit with resource augmentation

Tight results for Next Fit and Worst Fit with resource augmentation Tgh resuls for Nex F and Wors F wh resource augmenaon Joan Boyar Leah Epsen Asaf Levn Asrac I s well known ha he wo smple algorhms for he classc n packng prolem, NF and WF oh have an approxmaon rao of

More information

Efficient Asynchronous Channel Hopping Design for Cognitive Radio Networks

Efficient Asynchronous Channel Hopping Design for Cognitive Radio Networks Effcen Asynchronous Channel Hoppng Desgn for Cognve Rado Neworks Chh-Mn Chao, Chen-Yu Hsu, and Yun-ng Lng Absrac In a cognve rado nework (CRN), a necessary condon for nodes o communcae wh each oher s ha

More information

Notes on the stability of dynamic systems and the use of Eigen Values.

Notes on the stability of dynamic systems and the use of Eigen Values. Noes on he sabl of dnamc ssems and he use of Egen Values. Source: Macro II course noes, Dr. Davd Bessler s Tme Seres course noes, zarads (999) Ineremporal Macroeconomcs chaper 4 & Techncal ppend, and Hamlon

More information

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim Korean J. Mah. 19 (2011), No. 3, pp. 263 272 GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS Youngwoo Ahn and Kae Km Absrac. In he paper [1], an explc correspondence beween ceran

More information

Dual Approximate Dynamic Programming for Large Scale Hydro Valleys

Dual Approximate Dynamic Programming for Large Scale Hydro Valleys Dual Approxmae Dynamc Programmng for Large Scale Hydro Valleys Perre Carpener and Jean-Phlppe Chanceler 1 ENSTA ParsTech and ENPC ParsTech CMM Workshop, January 2016 1 Jon work wh J.-C. Alas, suppored

More information

Modeling and Solving of Multi-Product Inventory Lot-Sizing with Supplier Selection under Quantity Discounts

Modeling and Solving of Multi-Product Inventory Lot-Sizing with Supplier Selection under Quantity Discounts nernaonal ournal of Appled Engneerng Research SSN 0973-4562 Volume 13, Number 10 (2018) pp. 8708-8713 Modelng and Solvng of Mul-Produc nvenory Lo-Szng wh Suppler Selecon under Quany Dscouns Naapa anchanaruangrong

More information

2 Aggregate demand in partial equilibrium static framework

2 Aggregate demand in partial equilibrium static framework Unversy of Mnnesoa 8107 Macroeconomc Theory, Sprng 2012, Mn 1 Fabrzo Perr Lecure 1. Aggregaon 1 Inroducon Probably so far n he macro sequence you have deal drecly wh represenave consumers and represenave

More information

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model Probablsc Model for Tme-seres Daa: Hdden Markov Model Hrosh Mamsuka Bonformacs Cener Kyoo Unversy Oulne Three Problems for probablsc models n machne learnng. Compung lkelhood 2. Learnng 3. Parsng (predcon

More information

Optimal environmental charges under imperfect compliance

Optimal environmental charges under imperfect compliance ISSN 1 746-7233, England, UK World Journal of Modellng and Smulaon Vol. 4 (28) No. 2, pp. 131-139 Opmal envronmenal charges under mperfec complance Dajn Lu 1, Ya Wang 2 Tazhou Insue of Scence and Technology,

More information

Robust and Accurate Cancer Classification with Gene Expression Profiling

Robust and Accurate Cancer Classification with Gene Expression Profiling Robus and Accurae Cancer Classfcaon wh Gene Expresson Proflng (Compuaonal ysems Bology, 2005) Auhor: Hafeng L, Keshu Zhang, ao Jang Oulne Background LDA (lnear dscrmnan analyss) and small sample sze problem

More information

An introduction to Support Vector Machine

An introduction to Support Vector Machine An nroducon o Suppor Vecor Machne 報告者 : 黃立德 References: Smon Haykn, "Neural Neworks: a comprehensve foundaon, second edon, 999, Chaper 2,6 Nello Chrsann, John Shawe-Tayer, An Inroducon o Suppor Vecor Machnes,

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 0 Canoncal Transformaons (Chaper 9) Wha We Dd Las Tme Hamlon s Prncple n he Hamlonan formalsm Dervaon was smple δi δ Addonal end-pon consrans pq H( q, p, ) d 0 δ q ( ) δq ( ) δ

More information

Lecture 6: Learning for Control (Generalised Linear Regression)

Lecture 6: Learning for Control (Generalised Linear Regression) Lecure 6: Learnng for Conrol (Generalsed Lnear Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure 6: RLSC - Prof. Sehu Vjayakumar Lnear Regresson

More information

2/20/2013. EE 101 Midterm 2 Review

2/20/2013. EE 101 Midterm 2 Review //3 EE Mderm eew //3 Volage-mplfer Model The npu ressance s he equalen ressance see when lookng no he npu ermnals of he amplfer. o s he oupu ressance. I causes he oupu olage o decrease as he load ressance

More information

Advanced Macroeconomics II: Exchange economy

Advanced Macroeconomics II: Exchange economy Advanced Macroeconomcs II: Exchange economy Krzyszof Makarsk 1 Smple deermnsc dynamc model. 1.1 Inroducon Inroducon Smple deermnsc dynamc model. Defnons of equlbrum: Arrow-Debreu Sequenal Recursve Equvalence

More information

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5 TPG460 Reservor Smulaon 08 page of 5 DISCRETIZATIO OF THE FOW EQUATIOS As we already have seen, fne dfference appromaons of he paral dervaves appearng n he flow equaons may be obaned from Taylor seres

More information

The Dynamic Programming Models for Inventory Control System with Time-varying Demand

The Dynamic Programming Models for Inventory Control System with Time-varying Demand The Dynamc Programmng Models for Invenory Conrol Sysem wh Tme-varyng Demand Truong Hong Trnh (Correspondng auhor) The Unversy of Danang, Unversy of Economcs, Venam Tel: 84-236-352-5459 E-mal: rnh.h@due.edu.vn

More information

2 Aggregate demand in partial equilibrium static framework

2 Aggregate demand in partial equilibrium static framework Unversy of Mnnesoa 8107 Macroeconomc Theory, Sprng 2009, Mn 1 Fabrzo Perr Lecure 1. Aggregaon 1 Inroducon Probably so far n he macro sequence you have deal drecly wh represenave consumers and represenave

More information

SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β

SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β SARAJEVO JOURNAL OF MATHEMATICS Vol.3 (15) (2007), 137 143 SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β M. A. K. BAIG AND RAYEES AHMAD DAR Absrac. In hs paper, we propose

More information

Epistemic Game Theory: Online Appendix

Epistemic Game Theory: Online Appendix Epsemc Game Theory: Onlne Appendx Edde Dekel Lucano Pomao Marcano Snscalch July 18, 2014 Prelmnares Fx a fne ype srucure T I, S, T, β I and a probably µ S T. Le T µ I, S, T µ, βµ I be a ype srucure ha

More information

Polymerization Technology Laboratory Course

Polymerization Technology Laboratory Course Prakkum Polymer Scence/Polymersaonsechnk Versuch Resdence Tme Dsrbuon Polymerzaon Technology Laboraory Course Resdence Tme Dsrbuon of Chemcal Reacors If molecules or elemens of a flud are akng dfferen

More information

II. Light is a Ray (Geometrical Optics)

II. Light is a Ray (Geometrical Optics) II Lgh s a Ray (Geomercal Opcs) IIB Reflecon and Refracon Hero s Prncple of Leas Dsance Law of Reflecon Hero of Aleandra, who lved n he 2 nd cenury BC, posulaed he followng prncple: Prncple of Leas Dsance:

More information

A Systematic Framework for Dynamically Optimizing Multi-User Wireless Video Transmission

A Systematic Framework for Dynamically Optimizing Multi-User Wireless Video Transmission A Sysemac Framework for Dynamcally Opmzng ul-user Wreless Vdeo Transmsson Fangwen Fu, haela van der Schaar Elecrcal Engneerng Deparmen, UCLA {fwfu, mhaela}@ee.ucla.edu Absrac In hs paper, we formulae he

More information

5th International Conference on Advanced Design and Manufacturing Engineering (ICADME 2015)

5th International Conference on Advanced Design and Manufacturing Engineering (ICADME 2015) 5h Inernaonal onference on Advanced Desgn and Manufacurng Engneerng (IADME 5 The Falure Rae Expermenal Sudy of Specal N Machne Tool hunshan He, a, *, La Pan,b and Bng Hu 3,c,,3 ollege of Mechancal and

More information

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Journal of Appled Mahemacs and Compuaonal Mechancs 3, (), 45-5 HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Sansław Kukla, Urszula Sedlecka Insue of Mahemacs,

More information

Bundling with Customer Self-Selection: A Simple Approach to Bundling Low Marginal Cost Goods On-Line Appendix

Bundling with Customer Self-Selection: A Simple Approach to Bundling Low Marginal Cost Goods On-Line Appendix Bundlng wh Cusomer Self-Selecon: A Smple Approach o Bundlng Low Margnal Cos Goods On-Lne Appendx Lorn M. H Unversy of Pennsylvana, Wharon School 57 Jon M. Hunsman Hall Phladelpha, PA 94 lh@wharon.upenn.edu

More information

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes. umercal negraon of he dffuson equaon (I) Fne dfference mehod. Spaal screaon. Inernal nodes. R L V For hermal conducon le s dscree he spaal doman no small fne spans, =,,: Balance of parcles for an nernal

More information

Lecture VI Regression

Lecture VI Regression Lecure VI Regresson (Lnear Mehods for Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure VI: MLSC - Dr. Sehu Vjayakumar Lnear Regresson Model M

More information

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation Global Journal of Pure and Appled Mahemacs. ISSN 973-768 Volume 4, Number 6 (8), pp. 89-87 Research Inda Publcaons hp://www.rpublcaon.com Exsence and Unqueness Resuls for Random Impulsve Inegro-Dfferenal

More information

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model BGC1: Survval and even hsory analyss Oslo, March-May 212 Monday May 7h and Tuesday May 8h The addve regresson model Ørnulf Borgan Deparmen of Mahemacs Unversy of Oslo Oulne of program: Recapulaon Counng

More information

Chapter 6: AC Circuits

Chapter 6: AC Circuits Chaper 6: AC Crcus Chaper 6: Oulne Phasors and he AC Seady Sae AC Crcus A sable, lnear crcu operang n he seady sae wh snusodal excaon (.e., snusodal seady sae. Complee response forced response naural response.

More information

10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current :

10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current : . A. IUITS Synopss : GOWTH OF UNT IN IUIT : d. When swch S s closed a =; = d. A me, curren = e 3. The consan / has dmensons of me and s called he nducve me consan ( τ ) of he crcu. 4. = τ; =.63, n one

More information

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s Ordnary Dfferenal Equaons n Neuroscence wh Malab eamples. Am - Gan undersandng of how o se up and solve ODE s Am Undersand how o se up an solve a smple eample of he Hebb rule n D Our goal a end of class

More information

Cubic Bezier Homotopy Function for Solving Exponential Equations

Cubic Bezier Homotopy Function for Solving Exponential Equations Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 omoopy Funcon for Solvng Eponenal Equaons S. S. Raml *,,. Mohamad Nor,a, N. S. Saharzan,b and M.

More information

( ) [ ] MAP Decision Rule

( ) [ ] MAP Decision Rule Announcemens Bayes Decson Theory wh Normal Dsrbuons HW0 due oday HW o be assgned soon Proec descrpon posed Bomercs CSE 90 Lecure 4 CSE90, Sprng 04 CSE90, Sprng 04 Key Probables 4 ω class label X feaure

More information

12d Model. Civil and Surveying Software. Drainage Analysis Module Detention/Retention Basins. Owen Thornton BE (Mech), 12d Model Programmer

12d Model. Civil and Surveying Software. Drainage Analysis Module Detention/Retention Basins. Owen Thornton BE (Mech), 12d Model Programmer d Model Cvl and Surveyng Soware Dranage Analyss Module Deenon/Reenon Basns Owen Thornon BE (Mech), d Model Programmer owen.hornon@d.com 4 January 007 Revsed: 04 Aprl 007 9 February 008 (8Cp) Ths documen

More information

FTCS Solution to the Heat Equation

FTCS Solution to the Heat Equation FTCS Soluon o he Hea Equaon ME 448/548 Noes Gerald Reckenwald Porland Sae Unversy Deparmen of Mechancal Engneerng gerry@pdxedu ME 448/548: FTCS Soluon o he Hea Equaon Overvew Use he forward fne d erence

More information

ISSN MIT Publications

ISSN MIT Publications MIT Inernaonal Journal of Elecrcal and Insrumenaon Engneerng Vol. 1, No. 2, Aug 2011, pp 93-98 93 ISSN 2230-7656 MIT Publcaons A New Approach for Solvng Economc Load Dspach Problem Ansh Ahmad Dep. of Elecrcal

More information

Volatility Interpolation

Volatility Interpolation Volaly Inerpolaon Prelmnary Verson March 00 Jesper Andreasen and Bran Huge Danse Mares, Copenhagen wan.daddy@danseban.com brno@danseban.com Elecronc copy avalable a: hp://ssrn.com/absrac=69497 Inro Local

More information

e-journal Reliability: Theory& Applications No 2 (Vol.2) Vyacheslav Abramov

e-journal Reliability: Theory& Applications No 2 (Vol.2) Vyacheslav Abramov June 7 e-ournal Relably: Theory& Applcaons No (Vol. CONFIDENCE INTERVALS ASSOCIATED WITH PERFORMANCE ANALYSIS OF SYMMETRIC LARGE CLOSED CLIENT/SERVER COMPUTER NETWORKS Absrac Vyacheslav Abramov School

More information

Computational results on new staff scheduling benchmark instances

Computational results on new staff scheduling benchmark instances TECHNICAL REPORT Compuaonal resuls on new saff shedulng enhmark nsanes Tm Curos Rong Qu ASAP Researh Group Shool of Compuer Sene Unersy of Nongham NG8 1BB Nongham UK Frs pulshed onlne: 19-Sep-2014 las

More information

National Exams December 2015 NOTES: 04-BS-13, Biology. 3 hours duration

National Exams December 2015 NOTES: 04-BS-13, Biology. 3 hours duration Naonal Exams December 205 04-BS-3 Bology 3 hours duraon NOTES: f doub exss as o he nerpreaon of any queson he canddae s urged o subm wh he answer paper a clear saemen of any assumpons made 2 Ths s a CLOSED

More information

CHAPTER 10: LINEAR DISCRIMINATION

CHAPTER 10: LINEAR DISCRIMINATION CHAPER : LINEAR DISCRIMINAION Dscrmnan-based Classfcaon 3 In classfcaon h K classes (C,C,, C k ) We defned dscrmnan funcon g j (), j=,,,k hen gven an es eample, e chose (predced) s class label as C f g

More information

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction ECOOMICS 35* -- OTE 9 ECO 35* -- OTE 9 F-Tess and Analyss of Varance (AOVA n he Smple Lnear Regresson Model Inroducon The smple lnear regresson model s gven by he followng populaon regresson equaon, or

More information

SUPPLEMENT TO EFFICIENT DYNAMIC MECHANISMS IN ENVIRONMENTS WITH INTERDEPENDENT VALUATIONS: THE ROLE OF CONTINGENT TRANSFERS

SUPPLEMENT TO EFFICIENT DYNAMIC MECHANISMS IN ENVIRONMENTS WITH INTERDEPENDENT VALUATIONS: THE ROLE OF CONTINGENT TRANSFERS SUPPLEMENT TO EFFICIENT DYNAMIC MECHANISMS IN ENVIRONMENTS WITH INTERDEPENDENT VALUATIONS: THE ROLE OF CONTINGENT TRANSFERS HENG LIU In hs onlne appendx, we dscuss budge balance and surplus exracon n dynamc

More information

Genetic Algorithm in Parameter Estimation of Nonlinear Dynamic Systems

Genetic Algorithm in Parameter Estimation of Nonlinear Dynamic Systems Genec Algorhm n Parameer Esmaon of Nonlnear Dynamc Sysems E. Paeraks manos@egnaa.ee.auh.gr V. Perds perds@vergna.eng.auh.gr Ah. ehagas kehagas@egnaa.ee.auh.gr hp://skron.conrol.ee.auh.gr/kehagas/ndex.hm

More information

(,,, ) (,,, ). In addition, there are three other consumers, -2, -1, and 0. Consumer -2 has the utility function

(,,, ) (,,, ). In addition, there are three other consumers, -2, -1, and 0. Consumer -2 has the utility function MACROECONOMIC THEORY T J KEHOE ECON 87 SPRING 5 PROBLEM SET # Conder an overlappng generaon economy le ha n queon 5 on problem e n whch conumer lve for perod The uly funcon of he conumer born n perod,

More information

Bandlimited channel. Intersymbol interference (ISI) This non-ideal communication channel is also called dispersive channel

Bandlimited channel. Intersymbol interference (ISI) This non-ideal communication channel is also called dispersive channel Inersymol nererence ISI ISI s a sgnal-dependen orm o nererence ha arses ecause o devaons n he requency response o a channel rom he deal channel. Example: Bandlmed channel Tme Doman Bandlmed channel Frequency

More information

2.1 Constitutive Theory

2.1 Constitutive Theory Secon.. Consuve Theory.. Consuve Equaons Governng Equaons The equaons governng he behavour of maerals are (n he spaal form) dρ v & ρ + ρdv v = + ρ = Conservaon of Mass (..a) d x σ j dv dvσ + b = ρ v& +

More information

. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue.

. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue. Lnear Algebra Lecure # Noes We connue wh he dscusson of egenvalues, egenvecors, and dagonalzably of marces We wan o know, n parcular wha condons wll assure ha a marx can be dagonalzed and wha he obsrucons

More information

Department of Economics University of Toronto

Department of Economics University of Toronto Deparmen of Economcs Unversy of Torono ECO408F M.A. Economercs Lecure Noes on Heeroskedascy Heeroskedascy o Ths lecure nvolves lookng a modfcaons we need o make o deal wh he regresson model when some of

More information

Time-interval analysis of β decay. V. Horvat and J. C. Hardy

Time-interval analysis of β decay. V. Horvat and J. C. Hardy Tme-nerval analyss of β decay V. Horva and J. C. Hardy Work on he even analyss of β decay [1] connued and resuled n he developmen of a novel mehod of bea-decay me-nerval analyss ha produces hghly accurae

More information

Let s treat the problem of the response of a system to an applied external force. Again,

Let s treat the problem of the response of a system to an applied external force. Again, Page 33 QUANTUM LNEAR RESPONSE FUNCTON Le s rea he problem of he response of a sysem o an appled exernal force. Agan, H() H f () A H + V () Exernal agen acng on nernal varable Hamlonan for equlbrum sysem

More information

Learning Objectives. Self Organization Map. Hamming Distance(1/5) Introduction. Hamming Distance(3/5) Hamming Distance(2/5) 15/04/2015

Learning Objectives. Self Organization Map. Hamming Distance(1/5) Introduction. Hamming Distance(3/5) Hamming Distance(2/5) 15/04/2015 /4/ Learnng Objecves Self Organzaon Map Learnng whou Exaples. Inroducon. MAXNET 3. Cluserng 4. Feaure Map. Self-organzng Feaure Map 6. Concluson 38 Inroducon. Learnng whou exaples. Daa are npu o he syse

More information

Dynamic Team Decision Theory

Dynamic Team Decision Theory Dynamc Team Decson Theory EECS 558 Proec Repor Shruvandana Sharma and Davd Shuman December, 005 I. Inroducon Whle he sochasc conrol problem feaures one decson maker acng over me, many complex conrolled

More information

Single-loop System Reliability-Based Design & Topology Optimization (SRBDO/SRBTO): A Matrix-based System Reliability (MSR) Method

Single-loop System Reliability-Based Design & Topology Optimization (SRBDO/SRBTO): A Matrix-based System Reliability (MSR) Method 10 h US Naonal Congress on Compuaonal Mechancs Columbus, Oho 16-19, 2009 Sngle-loop Sysem Relably-Based Desgn & Topology Opmzaon (SRBDO/SRBTO): A Marx-based Sysem Relably (MSR) Mehod Tam Nguyen, Junho

More information

EP2200 Queuing theory and teletraffic systems. 3rd lecture Markov chains Birth-death process - Poisson process. Viktoria Fodor KTH EES

EP2200 Queuing theory and teletraffic systems. 3rd lecture Markov chains Birth-death process - Poisson process. Viktoria Fodor KTH EES EP Queung heory and eleraffc sysems 3rd lecure Marov chans Brh-deah rocess - Posson rocess Vora Fodor KTH EES Oulne for oday Marov rocesses Connuous-me Marov-chans Grah and marx reresenaon Transen and

More information

MANY real-world applications (e.g. production

MANY real-world applications (e.g. production Barebones Parcle Swarm for Ineger Programmng Problems Mahamed G. H. Omran, Andres Engelbrech and Ayed Salman Absrac The performance of wo recen varans of Parcle Swarm Opmzaon (PSO) when appled o Ineger

More information

APOC #232 Capacity Planning for Fault-Tolerant All-Optical Network

APOC #232 Capacity Planning for Fault-Tolerant All-Optical Network APOC #232 Capacy Plannng for Faul-Toleran All-Opcal Nework Mchael Kwok-Shng Ho and Kwok-wa Cheung Deparmen of Informaon ngneerng The Chnese Unversy of Hong Kong Shan, N.T., Hong Kong SAR, Chna -mal: kwcheung@e.cuhk.edu.hk

More information

January Examinations 2012

January Examinations 2012 Page of 5 EC79 January Examnaons No. of Pages: 5 No. of Quesons: 8 Subjec ECONOMICS (POSTGRADUATE) Tle of Paper EC79 QUANTITATIVE METHODS FOR BUSINESS AND FINANCE Tme Allowed Two Hours ( hours) Insrucons

More information

Survival Analysis and Reliability. A Note on the Mean Residual Life Function of a Parallel System

Survival Analysis and Reliability. A Note on the Mean Residual Life Function of a Parallel System Communcaons n Sascs Theory and Mehods, 34: 475 484, 2005 Copyrgh Taylor & Francs, Inc. ISSN: 0361-0926 prn/1532-415x onlne DOI: 10.1081/STA-200047430 Survval Analyss and Relably A Noe on he Mean Resdual

More information

Clustering (Bishop ch 9)

Clustering (Bishop ch 9) Cluserng (Bshop ch 9) Reference: Daa Mnng by Margare Dunham (a slde source) 1 Cluserng Cluserng s unsupervsed learnng, here are no class labels Wan o fnd groups of smlar nsances Ofen use a dsance measure

More information

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS THE PREICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS INTROUCTION The wo dmensonal paral dfferenal equaons of second order can be used for he smulaon of compeve envronmen n busness The arcle presens he

More information

Bernoulli process with 282 ky periodicity is detected in the R-N reversals of the earth s magnetic field

Bernoulli process with 282 ky periodicity is detected in the R-N reversals of the earth s magnetic field Submed o: Suden Essay Awards n Magnecs Bernoull process wh 8 ky perodcy s deeced n he R-N reversals of he earh s magnec feld Jozsef Gara Deparmen of Earh Scences Florda Inernaonal Unversy Unversy Park,

More information

DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL

DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL Sco Wsdom, John Hershey 2, Jonahan Le Roux 2, and Shnj Waanabe 2 Deparmen o Elecrcal Engneerng, Unversy o Washngon, Seale, WA, USA

More information

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations.

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations. Soluons o Ordnary Derenal Equaons An ordnary derenal equaon has only one ndependen varable. A sysem o ordnary derenal equaons consss o several derenal equaons each wh he same ndependen varable. An eample

More information

Online Appendix for. Strategic safety stocks in supply chains with evolving forecasts

Online Appendix for. Strategic safety stocks in supply chains with evolving forecasts Onlne Appendx for Sraegc safey socs n supply chans wh evolvng forecass Tor Schoenmeyr Sephen C. Graves Opsolar, Inc. 332 Hunwood Avenue Hayward, CA 94544 A. P. Sloan School of Managemen Massachuses Insue

More information

Introduction to Boosting

Introduction to Boosting Inroducon o Boosng Cynha Rudn PACM, Prnceon Unversy Advsors Ingrd Daubeches and Rober Schapre Say you have a daabase of news arcles, +, +, -, -, +, +, -, -, +, +, -, -, +, +, -, + where arcles are labeled

More information

A Tour of Modeling Techniques

A Tour of Modeling Techniques A Tour of Modelng Technques John Hooker Carnege Mellon Unversy EWO Semnar February 8 Slde Oulne Med neger lnear (MILP) modelng Dsuncve modelng Knapsack modelng Consran programmng models Inegraed Models

More information

FI 3103 Quantum Physics

FI 3103 Quantum Physics /9/4 FI 33 Quanum Physcs Aleander A. Iskandar Physcs of Magnesm and Phooncs Research Grou Insu Teknolog Bandung Basc Conces n Quanum Physcs Probably and Eecaon Value Hesenberg Uncerany Prncle Wave Funcon

More information

Testing a new idea to solve the P = NP problem with mathematical induction

Testing a new idea to solve the P = NP problem with mathematical induction Tesng a new dea o solve he P = NP problem wh mahemacal nducon Bacground P and NP are wo classes (ses) of languages n Compuer Scence An open problem s wheher P = NP Ths paper ess a new dea o compare he

More information

EXECUTION COSTS IN FINANCIAL MARKETS WITH SEVERAL INSTITUTIONAL INVESTORS

EXECUTION COSTS IN FINANCIAL MARKETS WITH SEVERAL INSTITUTIONAL INVESTORS EXECUION COSS IN FINANCIAL MARKES WIH SEVERAL INSIUIONAL INVESORS Somayeh Moazen, Yuyng L, Kae Larson Cheron School of Compuer Scence Unversy of Waerloo, Waerloo, ON, Canada emal: {smoazen, yuyng, klarson}@uwaerlooca

More information

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC CH.3. COMPATIBILITY EQUATIONS Connuum Mechancs Course (MMC) - ETSECCPB - UPC Overvew Compably Condons Compably Equaons of a Poenal Vecor Feld Compably Condons for Infnesmal Srans Inegraon of he Infnesmal

More information

Li An-Ping. Beijing , P.R.China

Li An-Ping. Beijing , P.R.China A New Type of Cpher: DICING_csb L An-Png Bejng 100085, P.R.Chna apl0001@sna.com Absrac: In hs paper, we wll propose a new ype of cpher named DICING_csb, whch s derved from our prevous sream cpher DICING.

More information

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas)

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas) Lecure 8: The Lalace Transform (See Secons 88- and 47 n Boas) Recall ha our bg-cure goal s he analyss of he dfferenal equaon, ax bx cx F, where we emloy varous exansons for he drvng funcon F deendng on

More information