Aalborg Universitet. Published in: IEEE Transactions on Smart Grid. DOI (link to publication from Publisher): /TSG.2017.

Size: px
Start display at page:

Download "Aalborg Universitet. Published in: IEEE Transactions on Smart Grid. DOI (link to publication from Publisher): /TSG.2017."

Transcription

1 Aalborg Unvertet A Data-Drven Stochatc Reactve Power Optmzaton Conderng Uncertante n Actve Dtrbuton Networ and Decompoton Method Dng, ao; Yang, Qngrun; Yang, Yongheng; L, Cheng; Be, Zhaohong; Blaaberg, Frede Publhed n: IEEE ranacton on Smart Grd DOI (ln to publcaton from Publher): 0.09/SG Publcaton date: 07 Document Veron Accepted author manucrpt, peer revewed veron Ln to publcaton from Aalborg Unverty Ctaton for publhed veron (APA): Dng,., Yang, Q., Yang, Y., L, C., Be, Z., & Blaaberg, F. (07). A Data-Drven Stochatc Reactve Power Optmzaton Conderng Uncertante n Actve Dtrbuton Networ and Decompoton Method. IEEE ranacton on Smart Grd, PP(99), -. DOI: 0.09/SG General rght Copyrght and moral rght for the publcaton made acceble n the publc portal are retaned by the author and/or other copyrght owner and t a condton of acceng publcaton that uer recogne and abde by the legal requrement aocated wth thee rght.? Uer may download and prnt one copy of any publcaton from the publc portal for the purpoe of prvate tudy or reearch.? You may not further dtrbute the materal or ue t for any proft-mang actvty or commercal gan? You may freely dtrbute the URL dentfyng the publcaton n the publc portal? ae down polcy If you beleve that th document breache copyrght pleae contact u at vbn@aub.aau.d provdng detal, and we wll remove acce to the wor mmedately and nvetgate your clam. Downloaded from vbn.aau.d on: november 30, 07

2 A Data-Drven Stochatc Reactve Power Optmzaton Conderng Uncertante n Actve Dtrbuton Networ and Decompoton Method ao Dng, Member, IEEE, Qngrun Yang, Student Member, IEEE, Yongheng Yang, Member, IEEE, Cheng L, Student Member, IEEE, Zhaohong Be, Senor Member, IEEE, Frede Blaaberg, Fellow, IEEE Abtract o addre the uncertan output of dtrbuted generator (DG) for reactve power optmzaton n actve dtrbuton networ, the tochatc programg model wdely ued. he model employed to fnd an optmal control trategy wth mum expected networ lo whle atfyng all the phycal contrant. heren, the probablty dtrbuton of uncertante n the tochatc model alway pre-defned by the htorcal data. However, the emprcal dtrbuton can be baed due to a lmted amount of htorcal data and thu reult n a uboptmal control decon. herefore, n th paper, a data-drven modelng approach ntroduced to aume that the probablty dtrbuton from the htorcal data uncertan wthn a confdence et. Furthermore, a data-drven tochatc programg model formulated a a two-tage problem, where the frt-tage varable fnd the optmal control for dcrete reactve power compenaton equpment under the wort probablty dtrbuton of the econd tage recoure. he econd-tage varable are aduted to uncertan probablty dtrbuton. In partcular, th two-tage problem ha a pecal tructure o that the econd-tage problem can be drectly decompoed nto everal mall-cale ub-problem, whch can be handled n parallel wthout the nformaton of dual problem. Numercal tudy on two dtrbuton ytem ha been performed. Comparon wth the two-tage tochatc and robut approache demontrate the effectvene of the propoal. Index erm Stochatc optmzaton; reactve power optmzaton; column-and-contrant generaton algorthm; actve dtrbuton networ; dtrbuted generaton NOMENCLAURE Indce and Set,, Index for bue t Index for tme perod B Set of bue E Set of branche h wor wa upported n part by Natonal Key Reearch and Development Program of Chna (06YFB090904), n part by Natonal Natural Scence Foundaton of Chna (Grant ), n part by Chna Potdoctoral Scence Foundaton (05M580847), n part by Natural Scence Ba Reearch Plan n Shaanx Provnce of Chna (06JQ505), the proect of State Key Laboratory of Electrcal Inulaton and n part by Power Equpment n X an Jaotong Unverty (EIPE630).. Dng (e-mal: tdng5@mal.xtu.edu.cn), Q. Yang, C. L and Z. Be are wth the State Key Laboratory of Electrcal Inulaton and Power Equpment, Department of Electrcal Engneerng, X an Jaotong Unverty, X an, Shaanx, 70049, Chna; Y. Yang and F. Blaaberg are wth the Department of Energy echnology, Aalborg Unverty, Aalborg DK-90, Denmar (e-mal: yoy@et.aau.d; fbl@et.aau.d). Θ Set of branche wth tranformer Set of bue for reactve power compenator D Set of bue for hunt capactor/reactor () Set of all parent of bu () Set of all chldren of bu Confdence et of the probablty dtrbuton Y Feable regon of contnuou varable under -th cenaro Parameter M A large number me horzon N w Cardnalty of Θ N c Cardnalty of D N r Cardnalty of D N Number of cenaro K Number of obervaton for uncertan parameter n Number of tap rato at tranformer branch (, ) r, x Retance/reactance of branch (, ) b, Shunt uceptance from to ground C /C Upper/lower bound of hunt capactor/reactor capacty at bu W / W Upper/lower bound of tranformer rato lmt at branch (, ) c, Specfed operatonal tme for hunt capactor/reactor at bu w, Specfed operatonal tme for tranformer (, ) Step ze of hunt capactor/reactor at bu w, ap rato on -th level of the tranformer (, ) U /U Upper/lower bound of voltage magntude at bu I l Current capacty lmt of branch (, ) Q c, / Q c, Upper/lower bound of reactve power compenaton for contnuou reactve power compenator at bu u Uncertan parameter u Uncertan parameter under -th cenaro Number of auxlary bnary varable,0,...,, p 0 H, G Probablty from the htorcal data A parameter that can control the ze of the confdence et ung -norm to control the ze of the confdence et ung nf-norm to control the ze of the confdence et Confdence level Varable Actve/reactve power flow from bu to

3 U Voltage magntude of bu P, Q Inected actve/reactve power of bu l Squared branch current at branch (, ) w ap rato of the tranformer branch (, ) C Value of hunt capactor/reactor at bu o Optmal 0- decon on -th level of the tranformer (, ) ρ Optmal tep of hunt capactor/reactor at bu Q c, Value of reactve power compenaton for contnuou reactve power compenator at bu v Squared voltage magntude of bu z Dcrete decon varable y Contnuou decon varable y Contnuou decon varable under -th cenaro,0,..., Auxlary bnary varable to expre the nteger, varable ρ by bnary code D I. INRODUCION ISRIBUED networ, characterzed by ther motly radal topology, are featured wth heavly fluctuatng load, whch may lead to large power loe and voltage drop near the end of feeder, adverely affectng ndutral manufacture and daly lve. o mprove the power qualty, reactve power optmzaton, ervng for tertary voltage control (VC), am to mze the total tranmon loe and mprove the voltage profle by controllng reactve power compenator and tranformer tap rato over everal perod, whle atfyng pecfc phycal and operatng contrant. Generally, the controlled equpment can be clafed a contnuou and dcrete controllable devce. he dcrete controllable devce are controlled va wtchng on/off and they hould not be aduted qute frequently due to ther ervce lfetme and extng manufacture technque. hu, the total number of operatng tme of dcrete controllable devce lmted, whch lead to the development of the dynamc reactve power optmzaton (DRPO) model [-3]. h model actually a large-cale mxed-nteger nonlnear programg and everal technque ncludng ntellgent earche and tandard branch-and-bound/cut method were propoed to olve th complex model [4]-[6]. Wth a propoal of a two-tage mult-perod mxed-nteger convex model, [7] analyzed the tradeoff between r mtgaton and nvetment cot mzaton. In [8], a voltage ecurty contraned mult-perod optmal reactve power flow model wa propoed baed on the generalzed Bender decompoton method wth an optmal condton decompoton approach to olve t. However, the ze of data are a a reult of large-cale mxed-nteger nonlnear programg problem wth mult-perod, ncreang the computatonal burden and tme. Recently, the conc relaxaton technque wa tuded n dtrbuton networ, whch gve a ound oluton whle gnfcantly mprovng the computatonal performance [9]-[]. For ntance, n [], the econd-order cone to relax the non-convex power flow equaton were propoed n order to obtan a mxed nteger econd order coned programg model, after whch a entvty-baed relaxaton and decompoton method wa ntroduced to further mprove the computaton. After deterng the total ze of the dtrbuted energy torage (DES, e.g., battere) and optmal locaton for the DES, [3] appled the econd order cone programg relaxaton to obtan the globally optmal oluton and avod the problem of NP-hardne. Furthermore, [4] dealt wth a ont problem of reactve power optmzaton and networ reconfguraton to mze power loe and mprove the voltage profle, the orgnal non-convex model of whch wa converted nto a mxed nteger econd order cone programg model ung the econd-order cone relaxaton, the bg-m method and the pecewe lnearzaton technque. Neverthele, an ncreang number of dtrbuted generator (DG) ncludng wnd power and photovoltac (PV) cog nto dtrbuted networ nowaday. he dtrbuted networ ntegrated wth DG, termed a actve dtrbuted networ, are facng crtcal techncal challenge to tradtonal operaton due to the tochatc nature of DG, whch may reult n uncertan output, and thu everer voltage volaton. o cope wth the uncertan output of DG n the optmzaton operaton n actve dtrbuton networ, tochatc programg [5-7], chance-contraned baed tochatc programg [8-0] and robut optmzaton [-3] have been extenvely explored. For example, a mult-cenaro framewor for optmal power flow under the wort wnd cenaro and tranmon N- contngency to properly addre the uncertan wnd power generaton wa propoed n [4]. A tochatc mult-obectve framewor for dtrbuton feeder reconfguraton wa employed n [5], frtly convertng t nto pecfc detertc cenaro among random cenaro of wnd/load forecat varaton and then mplementng mult-obectve formulaton for each detertc cenaro n the frt tage. In [6], a chance-contraned programg for optmal power flow under uncertanty conderng nonlnear model wth multple uncertan nput wa tuded, where a bac-mappng approach and lnear approxmaton of nonlnear model equaton were performed. Furthermore, [7] converted the chance-contraned tochatc programg formulaton nto a lnear detertc problem and a decompoton-baed method to olve the day-ahead chedulng problem. Although lnearzed model enable to mprove computatonal effcency, the accuracy of lnearzaton hould be enured. Generally, tochatc programg method cannot cover all the poble realzaton of uncertante. In order to addre th problem, robut optmzaton wa propoed to mmunze agant the oluton wthn a gven uncertanty et. A preented n [8], a two-tage robut reactve power optmzaton to coordnate the dcrete and contnuou reactve power compenator wa et up, whle hedgng agant any poble realzaton wthn uncertan wnd power output. A mxed-nteger two-tage robut optmzaton formulaton and a decompoton algorthm n a mater-lave tructure to acheve mum networ loe were dcued n [9], conderng the wort condton over uncertanty et. Although the robut optmzaton can protect the ytem agant a pre-defned uncertanty et, t alway gve a more conervatve oluton than the tochatc approach. In practce, htorcal data of DG output may be avalable at ISO/RO. herefore, t poble to derve a more effcent oluton that robut whle le conervatve, whch ncorporate the uperorty of both tochatc and robut approache. Accordng to the htorcal data, a confdence et contructed for the probablty dtrbuton of the uncertante to fnd an optmal oluton under the wort probablty dtrbuton [30]-[35]. herefore, a data-drven two-tage tochatc dynamc reactve power optmzaton model developed n th wor to coordnate the dcrete and contnuou controllable

4 3 devce, whle addreng the uncertan DG output. he contrbuton of the paper are ummarzed a follow: ) It the frt tme to et up a data-drven tochatc programg model n the dtrbuton networ, where the econd order cone programg relaxaton utlzed to relax the nonconvex feable regon caued by the branch flow equaton. Furthermore, the dynamc reactve power optmzaton can be termed a a large-cale mxed-nteger econd order cone programg model. ) It found that the propoed model ha a pecal tructure n the econd-tage b-level model, where the feable regon of the uncertanty et dont wth the operatng regon. A a reult, a new column-and-contrant generaton algorthm propoed to decompoe the b-level problem nto everal mall-cale ub-problem to be handled n parallel, whch doe not requre the dualty nformaton a the tradtonal method. he ret of the paper organzed a follow: Secton II preent a general dynamc reactve power optmzaton baed on econd order cone programg relaxaton for actve dtrbuton networ. In Secton III, a data-drven tochatc reactve power optmzaton model propoed wth the conderaton of uncertan DG output. Furthermore, a new dualty-free baed column-and-contrant generaton algorthm preented to olve the propoed reactve power optmzaton model n Secton IV. In Secton V, numercal reult obtaned on a 33-bu ytem demontrate the effectvene of the propoal, whch alo compared wth the two tradtonal approache. Fnally, concluon are drawn n Secton VI. II. REACIVE POWER OPIMIZAION MODEL IN ACIVE DISRIBUION NEWORKS A. Formulaton of Reactve Power Optmzaton Model Dtrbuton networ, dfferent from tranmon networ, have the property that the topology radal, o t very common to utlze the branch flow formulaton for decrbng the power flow n dtrbuton networ [], [8], [36]. P H H r l, B Q G G xl + b, U B U U r H xg r x l,, E \ () U U r H xg r x l,, w H G lu,, E where, E \ denote, E, but,. he frt and econd equaton decrbe the actve and reactve power balance at each bu; the thrd and fourth equaton decrbe the voltage drop at each lne and tranformer; the lat equaton decrbe the relatonhp among voltage, current and power. he reactve power optmzaton problem eentally am to mze total power loe by controllng the reactve power compenator and tranformer tap rato over a gven number of tme horzon whle atfyng varou phycal contrant. Here, the reactve power compenator can be clafed a contnuou adutment equpment uch a DG output, and dcrete adutment equpment ncludng capactor ban. It common that the electrc devce ncludng tranformer tap rato and wtched capactor ban cannot be aduted very frequently due to the lmtaton of ther ervce lfetme and extng manufacture technque. herefore, the mum allowable operatonal tme hould be condered n the model and the reactve power optmzaton model can be exactly wrtten a follow Q t, t, ot c r l t () t, E.t, PDG, t PL, t H t H t r l t B, t,..., (3) U tc t+ Qc, t QL, t, B, G t G t xl t+ b, U t t,..., (4) n 0 U t U t r H t x G t r x l t, o t U t U t, w,, E /, t,..., (5),, t r H t x G t r t x t l t,,,,..., (6) H t G t l t U t E (7) n o, t,, 0 (8) U U t U, B (9),, 0 l t I E (0) Qc, Qc, t Qc,, \ D () n t 0 C t C t, D () C C t C, D (3) t t c,, D t o t o t, (4) w,, (5) t Z o, t 0,,,, (6), D (7) where () am to mze total networ lo over tme perod; (3)-(4) denote the power balance at each bu; (5)-(7) how the Ohm' law for each branch, ncludng (6) for tranformer branch; (8) how a choce contrant by whch only one trap rato level choen; (9)-(0) are contrant for voltage

5 4 magntude and branch current; () the contrant for the contnuou reactve power compenator; ()-(3) are the contrant for dcrete reactve power compenator; (4)-(5) are retrcton that the total allowable operatonal tme by dcrete adutment equpment hould be lmted. However, the model ()-(7) a mxed nteger nonlnear nonconvex programg whch very dffcult to olve. However, the non-convexty come only from the nonlnear power flow contrant. o addre th ue, the em-defnte programg (SDP) and econd order cone programg (SOCP) were propoed to convexfy the feable regon encloed by the power flow contrant [9]-[0]. It wa hown n [9]-[0] that SOCP and SDP relaxaton method are equvalent for the radal networ, but the computatonal tme from the former one much le than the latter one. h becaue both SOCP and SDP are olved by the tandard prmal-dual nteror pont method, but SOCP ha much better wore-cae complexty than SDP [37]. heoretcally, the complexty of SOCP O(n 3 ) wherea O(n 4 ) of SDP. Here, n the number of varable. hu, for a large power ytem wth numerou varable, SOCP would perform much fater than SDP and thu elected n th wor. B. SOCP Relaxaton for Reactve Power Optmzaton Model At frt, let U t v t (4)-(7), (9) wll become v t C t + Qc, t QL, t for B and then contrant +, G t G t xl t b v t, B, t,..., (8) v t v t r H t x G t r x l t, n 0, E /, t,..., (9) o t v t u t r H t x G t r t x t l t, w,, t,..., (0),,,, H t G t l t v t E () U v t U, B () he contrant n () a nonlnear equalty, reultng n the nonconvex problem. o addre th ue, the econd order cone relaxaton performed by relaxng the quadratc equalty nto nequalty, yeldng H t G t l t v t l t u t,, E (3) After th relaxaton, the orgnal reactve power optmzaton model wll lead to be a mxed nteger econd order cone programg model, but not a tandard mxed nteger econd order cone programg model nce there are tll many blnear term n the above model, and we can mplfy them by reformulaton n the appendx, leadng to (A0)-(A). Subequently, the reactve power optmzaton model n can be mathematcally formulated a a general problem a a y (4).t. zy, Y Az b, z 0, (5) y Cy f, Q y q c,,..., y d n Y (6) Dy g Gz, Ey u where C, Q, q, c, D, g, G, E and d are matrx/vector form wth repect to the orgnal model. III. DAA-DRIVEN SOCHASIC REACIVE POWER OPIMIZAION CONSIDERING UNCERAINIES In the lat ecton, the reactve power optmzaton model only conducted under a gven load demand curve over multple tme perod. However, to addre the uncertan generaton output of the dtrbuted generator (.e., u n (6)), the tochatc programg employed to coordnate the dcrete and contnuou reactve power compenator. Specfcally, the dcrete decon varable (.e., z n (4)) hould be detered before the uncertanty revealed nce uch equpment hould not be aduted qute frequently, wherea the contnuou decon varable (.e., y n (4)) can be flexble wth the revealed uncertanty. h framewor gve a two-tage framewor and for the N cenaro of uncertante from dcretzng the gven N u and the corre- ). he obectve functon probablty dtrbuton, uch that pondng probablty ( p,, N p u,, mze the total expected networ lo. hen, the general data-drven tochatc reactve power optmzaton model formulated a.t. zy, Y N p a y (7) Az b, z 0, (8),,,..., y Cy f Q y q c y d n Y (9) Dy g Gz, Ey u Due to the lmted nformaton from the htorcal data, the probablty dtrbuton of uncertante cannot be exactly detered by the data. A a reult, we allow the probablty dtrbuton of uncertante to be arbtrary wthn a pre-defned confdence et contructed from the htorcal data. hu, the propoed data-drven tochatc reactve power optmzaton model am to fnd the optmal oluton under the wort-cae probablty dtrbuton, uch that.t. z p yy N p a y (30) Az b, z 0, (3),,,..., y Cy f Q y q c y d n Y (3) Dy g Gz, Ey u

6 5 In [35], two popular confdence et baed on norm- and norm-nf were preented for, whch can be expreed a N N N p R p p0 R p 0, p p (33) N N R 0 R p p 0, N p p p p (34) Suppong N cenaro from K obervaton, we have the followng relatonhp between the number of htorcal data and : K / 0 Ne K 0 Ne N Pr p p (35) Pr p p (36) It can be found that the rght-hand de of (35)-(36) actually the confdence level of the confdence et. hen, the relatonhp between confdence level (.e., the rght-hand de of (35)-(36)) and the value of gven by N N ln (37) K N ln K (38) Furthermore, (37) and (38) how that wth the ncreae of the number of htorcal data,.e., M, the etmated probablty dtrbuton wll be cloer to t true dtrbuton. hat mean, wll become maller untl to zero. Moreover, for the ame, maller than. IV. COLUMN-AND-CONSRAIN GENERAION ALGORIHM he propoed data-drven tochatc reactve power optmzaton model can be cat a a two-tage optmzaton problem whch generally can be olved by the Bender decompoton method or tandard column-and-contrant generaton method (C&CG). hee method are mplemented n a mater-ubproblem framewor: ub-problem (SP) am to fnd the crtcal cenaro of the uncertan et for a gven frt-tage decon varable that provde an upper bound; then new varable and contrant are added to the mater problem (MP) to obtan a lower bound. he MP and SP are olved teratvely and the method top untl the gap between the upper and lower bound maller than a pre-et convergence tolerance. A. C&CG-Sub-problem For a gven pecfc frt-tage varable n the -th teraton a z, we can et up a econd-tage b-level - model from (30)-(3) to fnd the wort-cae cenaro, yeldng p y Y N p a y (39),,,..., y Cy f Q y q c y d n.t. Y Dy g Gz, Ey u,..., N (40) It can be oberved that the model (39)-(40) ha ome pecal properte: () the ub-feable regon (Y,,Y,,Y N ) are eparable; () the decon varable p are all nonnegatve; () the feable regon of and Y are abolutely dont. For the frt and econd properte that the ub-feable regon (Y,,Y,,Y N ) are eparable and the decon varable p are all nonnegatve, we can exchange the ummaton operator and operator, o the econd-tage - problem can be reformulated a p N p a y (4) y Y,,,..., y Cy f Q y q c y d n.t. Y Dy g Gz, Ey u,..., N (4) For convenence, let become.t. h p N y Y arg h ph a y and the above model y Y (43) a y (44),,,..., y Cy f Q y q c y d n.t. Y Dy g Gz, Ey u,..., N (45) Accordng to the property (), the feable regon Y for varable y and the feable regon for varable p are abolutely dont. hat mean, the feable regon of upper-level model doen t affect the lower-level model and for any gven value p, the optmal oluton y unque. A a reult, the blevel model can be olved by equentally olvng upper-level and lower-level model, repectvely. Moreover, the frt property tell that the ub-feable regon (Y,,Y,,Y N ) are eparable, o lower-level model of the b-level model can be further decompoed nto N ndependent optmzaton model. h gve the fact that the b-level model can be decoupled by the followng tructure: For each u, t generate a econd order cone programg model, uch that arg a y (46) h y Y,,,..., y Cy f Q y q c y d n.t. Y (47) Dy g Gz, Ey u It can be oberved that the above econd order cone programg model are N mall model, comparng to the orgnal model (43)-(45), nce (46)-(47) only contan varable y for each model wherea (43)-(45) contan (y,,y,,y N ) multaneouly n one model. Moreover, the N mall model can be handled n parallel. After obtanng the optmal oluton ( h,, h above N mall model, we have p N hp N ) for the (48)

7 6 hu, we can ee that the orgnal b-level model can be olved by N mall econd order cone programg model that can be handled n parallel and one mall lnear programg. When the SP olved, an optmal value Q z and the wort-cae probablty p are obtaned, whch n fact gve an upper bound for the orgnal model. hen, a et of extra varable y,+ and aocated contrant are generated and added nto mater problem by fxng the optmal probablty p from the above model n (48). If the SP feable, we can create varable, y and agn the followng contrant to C&CG-mater problem, whch called optmalty cut. N p a y (49), Cy f, Q y q c y d,,..., n (50),,, N (5),, Dy g Gz, Ey u,,..., where a dummy contnuou varable. If the SP nfeable, t poble to create varable x and agn the followng contrant to C&CG-mater problem, whch called feablty cut. Cy f, Q y q c y d,,..., n (5),,, N (53),, Dy g Gz, Ey u,,..., B. C&CG-Mater Problem he MP am to relax the orgnal optmzaton model and provde a lower bound. After K teraton have been preceded, the mater problem can be decrbed a follow: (54).t. N z Az b, z 0, (55) p a y,,,... K (56), Cy f, Q y q c y d,,..., n,,,,,,... K,,..., N (57),, Dy g Gz, Ey u,,,... K,,..., N (58) he above MP a tandard mxed nteger econd order cone programg model that can be ealy handled by the off-the-hell commercal olver, uch a MOSEK, CPLEX, GUROBI, etc.. Solvng the MP gve the optmal dcrete varable z and optmal contnuou varable (y,,,y, ) that are generated n SP for the uncertanty et. he SP and the MP are olved teratvely untl the gven convergence crtera atfed and thu the global optmal oluton obtaned. A preented n [4], the column-and-contrant generaton algorthm can be converged n fnte teraton where all poble realzaton of are needed to be enumerated [4]. However, t hould be noted that the propoed column-and-contrant generaton algorthm a lttle dfferent from that n [4], where the propoed model ha a pecal tructure, o that there no need a [4] to dualze the nner - b-level model nto a ngle level model to olve t. Wth repect to the pecal tructure, we propoe a novel C&CG decompoton method wth the nformaton of dualty. For the gven convergence error, the mplementaton tep of the propoed algorthm are gven n able I. ABLE I PROCEDURE OF C&CG ALGORIHM C&CG algorthm. Set LB=, UB=, =0;. Solve the mater problem (54)-(58). Derve an optmal oluton z and (, y,,,y,+ ) for =,,N. hen, update the Update UB={UB, lower bound LB c y ; 3. Fx z and olve the ubproblem (46)-(47) n parallel, yeldng an optmal value Q( z ) and wort-cae probablty p. 4 c z Q z }. If UB-LB<, terate; ele, go to tep 5. 5 Generate varable ( y,,, y,,, y N, ). Add the new varable and contrant to mater problem accordng to (49)-(53). Update =+, and go to tep. V. NUMERICAL ANALYSIS A. et Sytem and Data Collecton In th ecton, a 33-bu dtrbuton networ that plotted n Fg. analyzed to verfy the propoed method. We conder the tep of tap rato (R) of the tranformer n the ubtaton 0.0 and the range [0.94,.06]. wo wtchable capactor/reactor (SCR) are connected to bue {#3, #9} whoe capacty are both [-0.0, +0.0] MVar, where the tep are 0.00 and MVar. he mum operatng tme over 4 hour for SCR are 8 and 6, repectvely. Bede, fve DG are ntalled at bue {#9, #5, #8, #3, #33} wth the capacty beng 0. MW, 0. MW, 0.3 MW, 0.3 MW and 0.3 MW repectvely. he forecated load demand and DG generaton factor over 4 hour are depcted n Fg.., where t aumed that the uncertan DG output follow a multvarate normal dtrbuton wth the varance equvalent to /5 of the mean value (a..a., forecated value). We randomly generate 000 ample by Monte Carlo mulaton to mulate the et of the htorcal data. ang the detered by (37) for example, the relatonhp among, N and hown n Fg. 3. h reveal that for the gven number of ample, wth the ncreae of the number of cenaro N and confdence level, become larger and the uncertanty et wll become larger a well. It obvou that the ze of uncertanty et wll affect the optmal oluton, o n the followng tudy, we wll chooe dfferent N and to how the mpact of uncertanty et on the reactve power optmzaton model. he computatonal ta were performed on a.0 GHz peronal computer wth 4 GB RAM, and the propoed method wa programmed n MALAB where the mxed nteger econd order cone programg were olved ung CPLEX.5.

8 7 Load/Generaton Factor Fg.. opology of 33-bu ytem Load Factor DG Factor at bue #9, #5, #8 DG Factor at bue #3, # me Perod (hour) Fg.. Load/Generaton factor over 4 hour N Fg. 3. he relatonhp among, N and. B. Reult and Comparon on A 33-Bu et Sytem he propoed method compared wth the tradtonal two popular method, two-tage tochatc and robut optmzaton approache, denoted by S and R repectvely. o compare wth the tradtonal tochatc programg model wth the detertc multvarate normal dtrbuton, we olve the tradtonal model and fx the frt-tage decon varable. hen, we randomly chooe 0000 dfferent probablte from the uncertan et and olve the econd-tage problem for each gven probablty, where t found that the oluton wth the mum networ lo erved a the wort-cae cenaro for the tochatc approache, denoted a Swt. Furthermore, the comparon of the three method preented n able II. he reult how that the two-tage robut optmzaton method yeld the hghet networ lo (.8545 MW) and the two-tage tochatc optmzaton method arrve at the lowet networ lo (.07 MW). he two-tage robut 9 SCR DG 0 optmzaton method optmze the optmal oluton under the wort-cae for all the poble realzaton, whch lead to the larget optmal oluton. he two-tage tochatc optmzaton method neglect the uncertanty of probablty of each cenaro, whch lead to the mallet optmal oluton. Moreover, he two method alway yeld the ame oluton for dfferent. he propoed method under both and gve a mld optmal oluton and can be termed a a budget that can control the ze of uncertanty et and further affect the optmal oluton. Moreover, the networ lo from the wort cae of tochatc approach (.e, Swt) conderng uncertan probablty dtrbuton about 0%~30% larger than the tradtonal two-tage tochatc programg. Increang confdence level lead to a larger uncertanty et, o that the wort-cae oluton wll become larger. Comparng the propoed method wth the tradtonal tochatc approach, t can be oberved that the networ lo from the propoed method under both uncertan et nf and larger than the tradtonal tochatc programg, whle t maller than that from the wort cae of tochatc approach. In partcular, a maller confdence level lead to a larger gap between Swt and the propoed method. Bede, the networ lo by the propoed method under dfferent uncertanty et gve dfferent value, but for the ame confdence level, the optmal oluton very cloe and the optmal oluton under lghtly maller than that under nf. Fnally, the dcrete control acton by the three method are tuded and compared. ae the frt SCR for llutraton and Fg. 4 depct that four and eght operatng tme of SCR are obtaned by robut and tochatc optmzaton method, wherea the propoed method operated between 4 and 8 tme. Here, we only chooe =0.5 and =0.99 for comparon due to the lmted pace. It oberve that wth the ncreae of, the optmal control acton over 4 hour cloer to that of robut optmzaton method. h becaue the ncreae of wll enlarge the uncertanty et, whch cloer to the uncertanty et of robut optmzaton approach. able II. Comparon of networ lo by three method under dfferent nf Networ Lo (MW) R S Swt

9 8 Locaton of SCR Robut Optmzaton Method Stochatc Optmzaton Method Propoed Method wth =0.99 Propoed Method wth = me Perod/Hour Fg. 4. Comparon of SCR on the three method C. Comparon of Computatonal Performance Between the Propoed Method and radtonal Approache he comparon of computatonal performance among the three approache hown n able III, preentng the teraton (Iter.) and computatonal tme (me) of each method. For the two-tage tochatc prograg model, t need to olve a large-cale mxed nteger econd order cone programg, whch actually a ngle-level model that can be drectly handled by the off-the-helf olver. However, the computatonal tme ncreae gnfcantly wth the ncreae of the number of cenaro. For the two-tage robut optmzaton model, the number of teraton only 3, where a new wort-cae cenaro dentfed at each teraton. It very tme-conug becaue olvng the nner b-level - problem need to tae dual and furthermore to olve a large-cale mxed nteger econd order cone programg model. In contrat, the propoed data-drven tochatc programg model ha a very pecal tructure, n whch the feable regon of econd-tage problem dont wth the uncertanty et, o a new column-and-contrant generaton algorthm propoed to decompoe the SP nto N mall-cale SP that can be olved n parallel. Meanwhle, the SP are everal econd order cone programg model, dfferent from the robut optmzaton model where the SP a large-cale mxed nteger econd order cone programg. he computatonal tme can be further reduced gnfcantly. It oberved from able III that the propoed method much fater than the two-tage robut optmzaton method. Moreover, when N mall, two-tage tochatc programg model a lttle fater than the propoed method, but wth ncreang N, the two-tage tochatc programg model become gnfcantly lower due to the large number of varable and contrant from the cenaro, wherea the computatonal tme of the propoed model ncreae only lghtly than to the decompoton method. herefore, the propoed method perform fater than the two-tage tochatc programg model epecally for the cae wth a large number of cenaro. Another tet ytem from a 3-bu tet ytem wth 0 DG and fve wtchable capactor/reactor (SCR) connected to bu, 35, 54, 76, and 08, whch hown n Fg. 5. he detaled nformaton can be avalable from [8]. he comparon of computatonal performance among the three approache hown n able IV, where t can be oberved that the robut optmzaton need x teraton for convergence by ue of column-and-contrant generaton algorthm and the total tme about 473. he computatonal tme of the tochatc optmzaton wll ncreae gnfcantly wth ncreang the number of cenaro. h becaue the tochatc optmzaton model contan N et of decon varable and contrant. Large N wll gnfcantly ncreae the number of total decon varable and contrant and thu need more computatonal tme. A for the propoed method, the computatonal peed more than 0 tme fater than the robut and tochatc optmzaton model when N large. Snce the ncreae of N wll enlarge the uncertanty et. herefore, t need more teraton for convergence and the total computatonal tme wll ncreae a well. Fnally, t hould be mentoned that the mum gap of conc relaxaton for any tet ytem maller than 0-4 MW, uggetng that the econd order cone programg relaxaton alway exact to the orgnal nonconvex model. able III. Comparon of computatonal effcency by three method on 33-bu tet ytem Propoed Robut Stochatc N Iter. me () Iter. me () Iter. me () Subtaton Fg. 5. opology of 3-bu ytem able IV. Comparon of computatonal effcency by three method on 3-bu tet ytem Propoed Robut Stochatc N Iter. me () Iter. me () Iter. me ()

10 9 VI. CONCLUSIONS h wor propoe a data-drven tochatc reactve power optmzaton model to addre uncertan dtrbuted generator ntegrated nto actve dtrbuton networ. Accordng to the htorcal data, the propoed method contruct a confdence et for the probablty dtrbuton of the uncertante and am to fnd an optmal oluton under the wort probablty dtrbuton. Furthermore, conc relaxaton employed to utlze to relax the feable regon encloed by power flow equaton. It noted that the propoed model ha a pecal tructure, o that a new column-and-contrant generaton algorthm propoed to decompoe the econd-tage b-level nner problem nto everal mall-cale ubproblem that can be handled n parallel. he comparon wth the tradtonal two-tage tochatc and robut approache on two tet ytem how that the propoed model can acheve better optmal oluton and computatonal performance than tradtonal method. APPENDIX A dcued n Secton II, the orgnal reactve power optmzaton model not a tandard mxed nteger econd order cone programg model nce there are tll many blnear term n the above model. Now, we can mplfy them o a to contruct a tandard mxed nteger econd order cone programg model. () Reformulaton for contrant n ()-(4) and (8) he dcrete reactve power compenator n ()-(4) and (8) are nonnegatve nteger, rather than 0- bnary varable. For the tandard mxed nteger programg model, t expected to formulate the model wth bnary varable. herefore, we hould reformulate each nteger varable t nto a combnaton of 0- bnary varable. Snce any nteger number ha a unque bnary code, the bnary code of t can be expreed by the combnaton of bnary varable t, t,..., t a,0,, 0 t t t t,0,..., (A) Accordng to the bound contrant n () and (3) that C C t C C t C t, we can derve 0,0,, t t... t C C (A) Snce t, t,..., t 0,, the mum value,0,, hould be. herefore, the mum value of hould be C C C C log log (A3) v t C t become Accordng to (A) and (), v t C t v t C t 0,0,, C v t t v t v t... t v t (A4), (8) derve For convenence, let, t, tv t 0,0, C t t..., t + Q c, t QL, t +, G t G t xl t b v t Furthermore, t tv t,, mean of the bg-m approach, uch that, t,..., (A5) can be lnearzed by,,, M t t M t M t t v t M t gve,,, For (A6), nce, t 0,,, t 0 0, t 0 0, t v t 0, t gve M t M, we can fnd that D,,.., (A6) M t v t M, D,,, t 0 ; t v t. herefore, (A6) equvalent to t tv t,, and ()-(4) and (8) can be expreed a (A5) wth addtonal contrant (A6). Meanwhle, tang (A) nto (4) lead to, t, t c,, D t 0 (A7) () Reformulaton for contrant n (0) Smlar to the method for lnearzng blnear term o t v t can be alo lnear- v tc t, the blnear term, zed ung the bg-m approach. Let h, t o, tv t the contrant (0) contanng blnear term and o t v t wll become n h, t v t r H t xg t r t x tl t, 0 w,,,,, Mo, t h, Mo, t,,, t,..., (A8) M o t h t v t M o t,,.., n (A9) Accordng to the above reformulaton and relaxaton n ecton II, the reactve power optmzaton model can be cat a a tandard 0- mxed nteger econd order cone programg a follow:,,

11 0 Q t, t, ot c r l t (A0) t, E.t (3), (8), (0)-(), (5)-(6), (9), ()-(3), (A5)-(A9) (A) t0,, o t0,, Q c t Contnou (A) REFERENCES [] S. Cor, P. Marannno, N. Lognore, G. Morechn et al., Coordnaton between the reactve power chedulng functon and the herarchcal voltage control of the EHV ENEL ytem, IEEE ran. Power Syt., vol. 0, no., pp , 995. []. Dng, Q. Guo, H. Sun, et al., "A quadratc robut optmzaton model for automatc voltage control on wnd farm de," IEEE Power and Energy Socety General Meetng, Vancouver, Canada, 03, pp. -5. [3]. Dng, Q. Guo, H. Sun, et al., A robut two-level coordnated tatc voltage ecurty regon for centrally ntegrated wnd farm, IEEE ran. Smart Grd, vol. PP, no. 99, pp. -, 05. [4] J. Zhao, L. Jud n Z. Da, et al., Voltage tablty contraned dynamc optmal reactve power flow bed on branch-bound and prmal-dual nteror pont method, Int. J. Electr. Power Energy Syt., vol. 7, no.3, pp , 05. [5] B. Zhao, C. X. Guo, Y. J. Cao, A multagent-baed partcle warm optmzaton approach for optmal reactve power dpatch, IEEE ran. Power Syt., vol. 0, no., pp , 005. [6]. Dng, R. Bo, F. L, et al., A b-level branch and bound method for economc dpatch wth dont prohbted zone conderng networ loe, IEEE ran. Power Syt., vol. 30, no. 6, pp , 05. [7] J. Lopez, J. Contrera and J. R. S Mantovan, Reactve power plannng under condtonal-value-at-r aement ung chance-contraned optmzaton, IE Gener. ranm. Dtrb., vol. 9, no. 3, pp. 3-40, 04. [8] A. Rabee and M. Parnan, Voltage ecurty contraned mult-perod optmal reactve power flow ung bender and optmalty condton decompoton, IEEE ran. Power Syt., vol. 8, no., 03. [9] M. Farvar, and S. H. Low, Branch Flow Model: Relaxaton and Convexfcaton-Part I, IEEE ran. Power Syt., vol. 8, no. 3, pp , 03. [0] M. Farvar, and S. H. Low, Branch Flow Model: Relaxaton and Convexfcaton-Part II, IEEE ran. Power Syt., vol. 8, no. 3, pp , 03. [] B. Mohamadreza, H. M. Reza and G. Mehrdad, Second-order cone programg for optmal power flow n VSC-type-AC-DC grd, IEEE ran. Power Syt., vol. 8, no. 4, pp , 03. []. Dng, S. Lu, Z. Wu and Z. Be, Sentvty-baed relaxaton and decompoton method to dynamc reactve power optmzaton conderng DG n actve dtrbuton networ, IE Gener. ranm. Dtrb.,vol. PP, no. 99, pp. -, 06. [3] Q. L, R. Ayyanar and V. Vttal, Convex optmzaton for DES plannng and operaton n radal dtrbuton ytem wth hgh penetraton of photovoltac reource, IEEE ran. Sutan. Energy, vol. 7, no.3, 06. [4] Z. an, W. Wu, B. Zhang and A. Boe, Mxed-nteger econd-order cone programg model for VAR optmzaton and networ reconfguraton n actve dtrbuton, IE Gener. ranm. Dtrb., vol. 0, no. 8, pp , 06. [5] J. Brge and F. Louveaux, Introducton to tochatc programg, New Yor, NY,USA: Sprnger-Verlag,997 [6] A. Papavalou, S. Oren and R. O Nell, Reerve requrement for wnd power ntegraton: a cenaro-baed tochatc programg framewor, IEEE ran. Power Syt., vol. 55, no. 4, pp.97-06, 0. [7] C. Sahn, M. Shahdehpour and I. Ermen, Allocaton of hourly reerve veru demand repone for ecurty-contraned chedulng of tochatc wnd energy, IEEE ran. Sutan. Energy, vol. 4, no., pp.9-8, 03. [8] A. Charne, W. W. Cooper, Chance-contraned programg, Manage. Sc., vol. 6, no., pp , 959. [9] L. Roald, S. Mra,. Kraue and G. Anderon, Correctve control to handle forecat uncertanty: a chance contraned optmal power flow, IEEE ran. Power Syt., vol. PP, no. 99, pp. -, 06. [0] H. Heaz and H. M. Rad, Energy torage plannng n actve dtrbuton grd: a chance-contraned optmzaton wth non-parametrc probablty functon, IEEE ran. Smart Grd, vol. PP, no. 99, pp. -, 06. [] J. Lopez, D. Pozo, J. Contrera and J. Mantovan, A mult-obectve regret robut VAr plannng model, IEEE ran. Power Syt., vol. PP, no. 99, pp. -, 06. [] R. A. Jabr, Adutable robut OPF wth renewable energy ource, IEEE ran. Power Syt., vol. 8, no. 4, pp , 03. [3] L. Zhhuan, Y. Ln and X. Duan, Non-doated ortng genetc algorthm-ii for robut mult-obectve optmal reactve power dpatch, IE Gener. ranm. Dtrb., vol. 4, no. 9, pp , 00. [4] X. Fang, F. L, Y. We, R. Azm and Y. Xu, Reactve power plannng under hgh penetraton of wnd energy ung Bender decompoton, IE Gener. ranm. Dtrb., vol. 9, no. 4, pp , 05. [5]. Nnam, A. Kavoufard and J. Aghae, Scenaro-baed mult-obectve dtrbuton feeder reconfguraton conderng wnd power ung adaptve modfed partcle warm optmzaton, IE Renewable Power Gener., vol. 6, no. 4,pp , 0. [6] H. Zhang and P. L. Chance contraned programg for optmal power flow under uncertanty, IEEE ran. Power Syt., vol. 6, no. 4, pp , 0. [7] H. Wu, M. Shahdehpour, Z. L and W. an, Chance-contraned day-ahead chedulng n tochatc power ytem operaton, IEEE ran. Power Syt., vol. 9, no. 4,pp ,04. [8]. Dng, S. Lu, W. Yuan Z. Be and B. Zeng, A two-tage robut reactve power optmzaton conderng uncertan wnd power ntegraton n actve dtrbuton networ, IEEE ran. Sutan. Energy, vol. 7, no., pp. 30-3,06. [9] H. Haghghat and B. Zeng, Dtrbuton ytem reconfguraton under uncertan load and renewable generaton, IEEE ran. Power Syt., vol. 3, no. 4, pp , 06. [30] Z. Wang, Q. Ban, H. Xn and D. Gan, A dtrbutonally robut co-ordnated reerve chedulng model conderng CVaR-baed wnd power reerve requrement, IEEE ran. Sutan. Energy, vol. 7, no., pp , 06. [3] P. Xong, P. Jruttaroen and C. Sngh, A dtrbutonally robut optmzaton model for unt commtment conderng uncertan wnd power generaton, IEEE ran. Power Syt., vol. PP, no. 99, pp.-, 06. [3] Y. Zhang, S. Shen and J. L. Matheu, Dtrbutonally robut chance-contraned optmal power flow wth uncertan renewable and uncertan reerve provded by load, IEEE ran. Power Syt., vol. PP, no. 99, pp. -, 06. [33] W. We, F. Lu and S. Me, Dtrbutonally robut co-optmzaton of energy and reerve dpatch, IEEE ran. Sutan. Energy, vol. 7, no., pp , 06. [34] Q. Ban, H. Xn, Z. Wang, D. Gan and K. P. Wong, Dtrbutonally robut oluton to the reerve chedulng problem wth partal nformaton of wnd power, IEEE ran. Power Syt., vol. 30, no. 5, pp. 8-83, 05. [35] C. Zhao and Y. Guan, Data-drven tochatc unt commtment for ntegratng wnd generaton, IEEE ran. Power Syt., vol. 3, no. 4, pp , 06. [36] N. L, L. Chen, S. H. Low, Exact convex relaxaton of OPF for radal networ ung branch flow model, IEEE hrd Internatonal Conference on Smart Grd Communcaton (SmartGrdComm), pp. 7-, 0. [37] M. S. Lobo, L. Vandenberghe, S. Boyd and et al, Applcaton of econd-order cone programg, Lnear algebra and t applcaton, vol. 84, no., pp. 93-8, 998. ao Dng (S 3 M 5) receved the B.S.E.E. and M.S.E.E. degree from Southeat Unverty, Nanng, Chna, n 009 and 0, repectvely, and the Ph.D. degree from nghua Unverty, Beng, Chna, n 05. Durng 03~04, he wa a vtng cholar wth the Department of Electrcal Engneerng and Computer Scence, he Unverty of enneee, Knoxvlle (UK), N, USA. He receved the excellent mater and doctoral dertaton from Southeat Unverty and nghua Unverty, repectvely, and outtandng graduate award of Beng Cty. Dr. Dng ha publhed more than 60 techncal paper and authored by Sprnger hee recognzng outtandng Ph.D. reearch around the world and acro the phycal cence Power Sytem Operaton wth Large Scale Stochatc Wnd Power Integraton. He currently an aocate profeor n the State Key Laboratory of Electrcal Inulaton and Power Equpment, the School of Electrcal Engneerng,

12 X an Jaotong Unverty. H current reearch nteret nclude electrcty maret, power ytem economc and optmzaton method, and power ytem plannng and relablty evaluaton. QngrunYang receved the B.S. degree from the School of Electrcal Engnee rng, X an Jaotong Unverty, X an, Chna, n 07. He currently purun g the M.S. degree at X an Jaotong Unverty. H maor reearch nteret n clude power ytem optmzaton and demand repone. Yongheng Yang (S - M 5) receved the B.Eng. degree n 009 from Northwetern Polytechncal Unverty, Chna and the Ph.D. degree n 04 from Aalborg Unverty, Denmar. He wa a potgraduate wth Southeat Unverty, Chna, from 009 to 0. In 03, he wa a Vtng Scholar wth exa A&M Unverty, USA. Snce 04, he ha been wth the Department of Energy echnology, Aalborg Unverty, where currently he an Atant Profeor. H reearch nteret are focued on grd ntegraton of renewable energy ytem, power converter degn, analy and control, harmonc dentfcaton and mtgaton, and relablty n power electronc. Dr. Yang ha publhed more than 80 techncal paper and co-authored a boo Perodc Control of Power Electronc Converter (London, UK: IE, 07). Dr. Yang a Member of the IEEE Power Electronc Socety (PELS) Student and Young Profeonal Commttee, where he erve a the Global Strategy Char and reponble for the IEEE PELS Student and Young Profeonal Actvte. He erved a a Guet Aocate Edtor of IEEE JOURNAL OF EMERGING AND SELECED OPICS IN POWER ELECRONICS, and ha alo been nvted a a Guet Edtor of Appled Scence. He an actve revewer for relevant top-ter ournal. Cheng L receved the B.S. degree from the School of Electrcal Engneerng, X an Jaotong Unverty, X an, Chna, n 06. He currently purung the M.S. degree at X an Jaotong Unverty. H maor reearch nteret nclude power ytem optmzaton and renewable energy ntegraton. Zhaohong Be (M 98 SM ) receved the B.S. and M.S. degree from the Electrc Power Department of Shandong Unverty, Jnan, Chna, n 99 and 994, repectvely, and the Ph.D. degree from X an Jaotong Unverty, X an, Chna, n 998. Currently, he a Profeor n the State Key Laboratory of Electrcal Inulaton and Power Equpment and the School of Electrcal Engneerng, X an Jaotong Unverty. Her man nteret are power ytem plannng and relablty evaluaton, a well a the ntegraton of the renewable energy. Frede Blaaberg (S 86 - M 88 - SM 97 - F 03) wa wth ABB-Scanda, Rander, Denmar, from 987 to 988. From 988 to 99, he wa a Ph.D. Student wth Aalborg Unverty, Aalborg, Denmar. He became an Atant Profeor n 99, an Aocate Profeor n 996, and a Full Profeor of power electronc and drve n 998 at Aalborg Unverty. H current reearch nteret nclude power electronc and t applcaton uch a n wnd turbne, PV ytem, relablty, harmonc and adutable peed drve. Prof. Blaaberg ha receved 7 IEEE Prze Paper Award, the IEEE Power Electronc Socety (PELS) Dtnguhed Servce Award n 009, the EPE-PEMC Councl Award n 00, the IEEE Wllam E. Newell Power Electronc Award n 04 and the Vllum Kann Ramuen Reearch Award n 04. He wa an Edtor-n-Chef of the IEEE RANSACIONS ON POWER ELECRONICS from 006 to 0. Prof. Blaaberg wa noated n 04 and 05 by homon Reuter to be between the mot 50 cted reearcher n Engneerng n the world.

Additional File 1 - Detailed explanation of the expression level CPD

Additional File 1 - Detailed explanation of the expression level CPD Addtonal Fle - Detaled explanaton of the expreon level CPD A mentoned n the man text, the man CPD for the uterng model cont of two ndvdual factor: P( level gen P( level gen P ( level gen 2 (.).. CPD factor

More information

Distributed Control for the Parallel DC Linked Modular Shunt Active Power Filters under Distorted Utility Voltage Condition

Distributed Control for the Parallel DC Linked Modular Shunt Active Power Filters under Distorted Utility Voltage Condition Dtrbted Control for the Parallel DC Lnked Modlar Shnt Actve Power Flter nder Dtorted Utlty Voltage Condton Reearch Stdent: Adl Salman Spervor: Dr. Malabka Ba School of Electrcal and Electronc Engneerng

More information

Team. Outline. Statistics and Art: Sampling, Response Error, Mixed Models, Missing Data, and Inference

Team. Outline. Statistics and Art: Sampling, Response Error, Mixed Models, Missing Data, and Inference Team Stattc and Art: Samplng, Repone Error, Mxed Model, Mng Data, and nference Ed Stanek Unverty of Maachuett- Amhert, USA 9/5/8 9/5/8 Outlne. Example: Doe-repone Model n Toxcology. ow to Predct Realzed

More information

The multivariate Gaussian probability density function for random vector X (X 1,,X ) T. diagonal term of, denoted

The multivariate Gaussian probability density function for random vector X (X 1,,X ) T. diagonal term of, denoted Appendx Proof of heorem he multvarate Gauan probablty denty functon for random vector X (X,,X ) px exp / / x x mean and varance equal to the th dagonal term of, denoted he margnal dtrbuton of X Gauan wth

More information

Specification -- Assumptions of the Simple Classical Linear Regression Model (CLRM) 1. Introduction

Specification -- Assumptions of the Simple Classical Linear Regression Model (CLRM) 1. Introduction ECONOMICS 35* -- NOTE ECON 35* -- NOTE Specfcaton -- Aumpton of the Smple Clacal Lnear Regreon Model (CLRM). Introducton CLRM tand for the Clacal Lnear Regreon Model. The CLRM alo known a the tandard lnear

More information

Improvements on Waring s Problem

Improvements on Waring s Problem Improvement on Warng Problem L An-Png Bejng, PR Chna apl@nacom Abtract By a new recurve algorthm for the auxlary equaton, n th paper, we wll gve ome mprovement for Warng problem Keyword: Warng Problem,

More information

Small signal analysis

Small signal analysis Small gnal analy. ntroducton Let u conder the crcut hown n Fg., where the nonlnear retor decrbed by the equaton g v havng graphcal repreentaton hown n Fg.. ( G (t G v(t v Fg. Fg. a D current ource wherea

More information

Start Point and Trajectory Analysis for the Minimal Time System Design Algorithm

Start Point and Trajectory Analysis for the Minimal Time System Design Algorithm Start Pont and Trajectory Analy for the Mnmal Tme Sytem Degn Algorthm ALEXANDER ZEMLIAK, PEDRO MIRANDA Department of Phyc and Mathematc Puebla Autonomou Unverty Av San Claudo /n, Puebla, 757 MEXICO Abtract:

More information

Harmonic oscillator approximation

Harmonic oscillator approximation armonc ocllator approxmaton armonc ocllator approxmaton Euaton to be olved We are fndng a mnmum of the functon under the retrcton where W P, P,..., P, Q, Q,..., Q P, P,..., P, Q, Q,..., Q lnwgner functon

More information

Solution Methods for Time-indexed MIP Models for Chemical Production Scheduling

Solution Methods for Time-indexed MIP Models for Chemical Production Scheduling Ian Davd Lockhart Bogle and Mchael Farweather (Edtor), Proceedng of the 22nd European Sympoum on Computer Aded Proce Engneerng, 17-2 June 212, London. 212 Elever B.V. All rght reerved. Soluton Method for

More information

Introduction to Interfacial Segregation. Xiaozhe Zhang 10/02/2015

Introduction to Interfacial Segregation. Xiaozhe Zhang 10/02/2015 Introducton to Interfacal Segregaton Xaozhe Zhang 10/02/2015 Interfacal egregaton Segregaton n materal refer to the enrchment of a materal conttuent at a free urface or an nternal nterface of a materal.

More information

728. Mechanical and electrical elements in reduction of vibrations

728. Mechanical and electrical elements in reduction of vibrations 78. Mechancal and electrcal element n reducton of vbraton Katarzyna BIAŁAS The Slean Unverty of Technology, Faculty of Mechancal Engneerng Inttute of Engneerng Procee Automaton and Integrated Manufacturng

More information

Estimation of Finite Population Total under PPS Sampling in Presence of Extra Auxiliary Information

Estimation of Finite Population Total under PPS Sampling in Presence of Extra Auxiliary Information Internatonal Journal of Stattc and Analy. ISSN 2248-9959 Volume 6, Number 1 (2016), pp. 9-16 Reearch Inda Publcaton http://www.rpublcaton.com Etmaton of Fnte Populaton Total under PPS Samplng n Preence

More information

Discrete Simultaneous Perturbation Stochastic Approximation on Loss Function with Noisy Measurements

Discrete Simultaneous Perturbation Stochastic Approximation on Loss Function with Noisy Measurements 0 Amercan Control Conference on O'Farrell Street San Francco CA USA June 9 - July 0 0 Dcrete Smultaneou Perturbaton Stochatc Approxmaton on Lo Functon wth Noy Meaurement Q Wang and Jame C Spall Abtract

More information

Method Of Fundamental Solutions For Modeling Electromagnetic Wave Scattering Problems

Method Of Fundamental Solutions For Modeling Electromagnetic Wave Scattering Problems Internatonal Workhop on MehFree Method 003 1 Method Of Fundamental Soluton For Modelng lectromagnetc Wave Scatterng Problem Der-Lang Young (1) and Jhh-We Ruan (1) Abtract: In th paper we attempt to contruct

More information

Chapter 6 The Effect of the GPS Systematic Errors on Deformation Parameters

Chapter 6 The Effect of the GPS Systematic Errors on Deformation Parameters Chapter 6 The Effect of the GPS Sytematc Error on Deformaton Parameter 6.. General Beutler et al., (988) dd the frt comprehenve tudy on the GPS ytematc error. Baed on a geometrc approach and aumng a unform

More information

Chapter 11. Supplemental Text Material. The method of steepest ascent can be derived as follows. Suppose that we have fit a firstorder

Chapter 11. Supplemental Text Material. The method of steepest ascent can be derived as follows. Suppose that we have fit a firstorder S-. The Method of Steepet cent Chapter. Supplemental Text Materal The method of teepet acent can be derved a follow. Suppoe that we have ft a frtorder model y = β + β x and we wh to ue th model to determne

More information

On the SO 2 Problem in Thermal Power Plants. 2.Two-steps chemical absorption modeling

On the SO 2 Problem in Thermal Power Plants. 2.Two-steps chemical absorption modeling Internatonal Journal of Engneerng Reearch ISSN:39-689)(onlne),347-53(prnt) Volume No4, Iue No, pp : 557-56 Oct 5 On the SO Problem n Thermal Power Plant Two-tep chemcal aborpton modelng hr Boyadjev, P

More information

Confidence intervals for the difference and the ratio of Lognormal means with bounded parameters

Confidence intervals for the difference and the ratio of Lognormal means with bounded parameters Songklanakarn J. Sc. Technol. 37 () 3-40 Mar.-Apr. 05 http://www.jt.pu.ac.th Orgnal Artcle Confdence nterval for the dfference and the rato of Lognormal mean wth bounded parameter Sa-aat Nwtpong* Department

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

The Study of Teaching-learning-based Optimization Algorithm

The Study of Teaching-learning-based Optimization Algorithm Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute

More information

ENTROPY BOUNDS USING ARITHMETIC- GEOMETRIC-HARMONIC MEAN INEQUALITY. Guru Nanak Dev University Amritsar, , INDIA

ENTROPY BOUNDS USING ARITHMETIC- GEOMETRIC-HARMONIC MEAN INEQUALITY. Guru Nanak Dev University Amritsar, , INDIA Internatonal Journal of Pure and Appled Mathematc Volume 89 No. 5 2013, 719-730 ISSN: 1311-8080 prnted veron; ISSN: 1314-3395 on-lne veron url: http://.jpam.eu do: http://dx.do.org/10.12732/jpam.v895.8

More information

arxiv: v1 [cs.gt] 15 Jan 2019

arxiv: v1 [cs.gt] 15 Jan 2019 Model and algorthm for tme-content rk-aware Markov game Wenje Huang, Pham Vet Ha and Wllam B. Hakell January 16, 2019 arxv:1901.04882v1 [c.gt] 15 Jan 2019 Abtract In th paper, we propoe a model for non-cooperatve

More information

Variable Structure Control ~ Basics

Variable Structure Control ~ Basics Varable Structure Control ~ Bac Harry G. Kwatny Department of Mechancal Engneerng & Mechanc Drexel Unverty Outlne A prelmnary example VS ytem, ldng mode, reachng Bac of dcontnuou ytem Example: underea

More information

Pythagorean triples. Leen Noordzij.

Pythagorean triples. Leen Noordzij. Pythagorean trple. Leen Noordz Dr.l.noordz@leennoordz.nl www.leennoordz.me Content A Roadmap for generatng Pythagorean Trple.... Pythagorean Trple.... 3 Dcuon Concluon.... 5 A Roadmap for generatng Pythagorean

More information

EEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming

EEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming EEL 6266 Power System Operaton and Control Chapter 3 Economc Dspatch Usng Dynamc Programmng Pecewse Lnear Cost Functons Common practce many utltes prefer to represent ther generator cost functons as sngle-

More information

MULTIPLE REGRESSION ANALYSIS For the Case of Two Regressors

MULTIPLE REGRESSION ANALYSIS For the Case of Two Regressors MULTIPLE REGRESSION ANALYSIS For the Cae of Two Regreor In the followng note, leat-quare etmaton developed for multple regreon problem wth two eplanator varable, here called regreor (uch a n the Fat Food

More information

Experience with Automatic Generation Control (AGC) Dynamic Simulation in PSS E

Experience with Automatic Generation Control (AGC) Dynamic Simulation in PSS E Semens Industry, Inc. Power Technology Issue 113 Experence wth Automatc Generaton Control (AGC) Dynamc Smulaton n PSS E Lu Wang, Ph.D. Staff Software Engneer lu_wang@semens.com Dngguo Chen, Ph.D. Staff

More information

A Simple Heuristic for Reducing the Number of Scenarios in Two-stage Stochastic Programming

A Simple Heuristic for Reducing the Number of Scenarios in Two-stage Stochastic Programming A Smle Heurtc for Reducng the Number of Scenaro n wo-tage Stochatc Programmng Ramumar aruah Marano Martn and gnaco E. Gromann * Deartment of Chemcal Engneerng Carnege Mellon Unverty Pttburgh PA 5 U.S.A.

More information

A Weighted UTASTAR Method for the Multiple Criteria Decision Making with Interval Numbers

A Weighted UTASTAR Method for the Multiple Criteria Decision Making with Interval Numbers 3rd Internatonal Conference on Management Scence and Management Innovaton MSMI 2016) A Weghted UTASTAR Method for the Multple Crtera Decon Makng wth Interval Number Wen-Tao Xong Jng Cheng School of Mathematc

More information

Research Article Runge-Kutta Type Methods for Directly Solving Special Fourth-Order Ordinary Differential Equations

Research Article Runge-Kutta Type Methods for Directly Solving Special Fourth-Order Ordinary Differential Equations Hndaw Publhng Corporaton Mathematcal Problem n Engneerng Volume 205, Artcle ID 893763, page http://dx.do.org/0.55/205/893763 Reearch Artcle Runge-Kutta Type Method for Drectly Solvng Specal Fourth-Order

More information

KEY POINTS FOR NUMERICAL SIMULATION OF INCLINATION OF BUILDINGS ON LIQUEFIABLE SOIL LAYERS

KEY POINTS FOR NUMERICAL SIMULATION OF INCLINATION OF BUILDINGS ON LIQUEFIABLE SOIL LAYERS KY POINTS FOR NUMRICAL SIMULATION OF INCLINATION OF BUILDINGS ON LIQUFIABL SOIL LAYRS Jn Xu 1, Xaomng Yuan, Jany Zhang 3,Fanchao Meng 1 1 Student, Dept. of Geotechncal ngneerng, Inttute of ngneerng Mechanc,

More information

Root Locus Techniques

Root Locus Techniques Root Locu Technque ELEC 32 Cloed-Loop Control The control nput u t ynthezed baed on the a pror knowledge of the ytem plant, the reference nput r t, and the error gnal, e t The control ytem meaure the output,

More information

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

More information

Improvements on Waring s Problem

Improvements on Waring s Problem Imrovement on Warng Problem L An-Png Bejng 85, PR Chna al@nacom Abtract By a new recurve algorthm for the auxlary equaton, n th aer, we wll gve ome mrovement for Warng roblem Keyword: Warng Problem, Hardy-Lttlewood

More information

Optimal Dual-Connectivity Traffic Offloading in Energy-Harvesting Small-Cell Networks

Optimal Dual-Connectivity Traffic Offloading in Energy-Harvesting Small-Cell Networks Th artcle ha been accepted for publcaton n a future ue of th journal, but ha not been fully edted. Content may change pror to fnal publcaton. Ctaton nformaton: DOI.9/TGCN.28.285843, IEEE Tranacton on Green

More information

Comparative Study on Electromagnetic and Electromechanical Transient Model for Grid-connected Photovoltaic Power System

Comparative Study on Electromagnetic and Electromechanical Transient Model for Grid-connected Photovoltaic Power System Energy and Power Engneerng, 13, 5, 47-5 do:1.436/epe.13.54b48 Publhed Onlne July 13 (http://www.crp.org/journal/epe) Comparatve Study on and Tranent Model for Grd-connected Photovoltac Power Sytem Man

More information

A NUMERICAL MODELING OF MAGNETIC FIELD PERTURBATED BY THE PRESENCE OF SCHIP S HULL

A NUMERICAL MODELING OF MAGNETIC FIELD PERTURBATED BY THE PRESENCE OF SCHIP S HULL A NUMERCAL MODELNG OF MAGNETC FELD PERTURBATED BY THE PRESENCE OF SCHP S HULL M. Dennah* Z. Abd** * Laboratory Electromagnetc Sytem EMP BP b Ben-Aknoun 606 Alger Algera ** Electronc nttute USTHB Alger

More information

A Preliminary Study on Material Utilization of Stiffened Cylindrical Shells

A Preliminary Study on Material Utilization of Stiffened Cylindrical Shells Reearch Journal of Appled Scence, Engneerng and echnology 6(5): 757-763, 03 ISSN: 040-7459; e-issn: 040-7467 Maxwell Scentfc Organzaton, 03 Submtted: December 8, 0 Accepted: February 08, 03 Publhed: Augut

More information

Information Acquisition in Global Games of Regime Change (Online Appendix)

Information Acquisition in Global Games of Regime Change (Online Appendix) Informaton Acquton n Global Game of Regme Change (Onlne Appendx) Mchal Szkup and Iabel Trevno Augut 4, 05 Introducton Th appendx contan the proof of all the ntermedate reult that have been omtted from

More information

ELG3336: Op Amp-based Active Filters

ELG3336: Op Amp-based Active Filters ELG6: Op Amp-baed Actve Flter Advantage: educed ze and weght, and thereore paratc. Increaed relablty and mproved perormance. Smpler degn than or pave lter and can realze a wder range o uncton a well a

More information

Chapter - 2. Distribution System Power Flow Analysis

Chapter - 2. Distribution System Power Flow Analysis Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load

More information

Foresighted Resource Reciprocation Strategies in P2P Networks

Foresighted Resource Reciprocation Strategies in P2P Networks Foreghted Reource Recprocaton Stratege n PP Networ Hyunggon Par and Mhaela van der Schaar Electrcal Engneerng Department Unverty of Calforna Lo Angele (UCLA) Emal: {hgpar mhaela@ee.ucla.edu Abtract We

More information

A Hybrid Evolution Algorithm with Application Based on Chaos Genetic Algorithm and Particle Swarm Optimization

A Hybrid Evolution Algorithm with Application Based on Chaos Genetic Algorithm and Particle Swarm Optimization Natonal Conference on Informaton Technology and Computer Scence (CITCS ) A Hybrd Evoluton Algorthm wth Applcaton Baed on Chao Genetc Algorthm and Partcle Swarm Optmzaton Fu Yu School of Computer & Informaton

More information

Modeling of Wave Behavior of Substrate Noise Coupling for Mixed-Signal IC Design

Modeling of Wave Behavior of Substrate Noise Coupling for Mixed-Signal IC Design Modelng of Wave Behavor of Subtrate Noe Couplng for Mxed-Sgnal IC Degn Georgo Veron, Y-Chang Lu, and Robert W. Dutton Center for Integrated Sytem, Stanford Unverty, Stanford, CA 9435 yorgo@gloworm.tanford.edu

More information

MODELLING OF STOCHASTIC PARAMETERS FOR CONTROL OF CITY ELECTRIC TRANSPORT SYSTEMS USING EVOLUTIONARY ALGORITHM

MODELLING OF STOCHASTIC PARAMETERS FOR CONTROL OF CITY ELECTRIC TRANSPORT SYSTEMS USING EVOLUTIONARY ALGORITHM MODELLING OF STOCHASTIC PARAMETERS FOR CONTROL OF CITY ELECTRIC TRANSPORT SYSTEMS USING EVOLUTIONARY ALGORITHM Mkhal Gorobetz, Anatoly Levchenkov Inttute of Indutral Electronc and Electrotechnc, Rga Techncal

More information

Robust Capacitated Facility Location Problem: Optimization Model and Solution Algorithms

Robust Capacitated Facility Location Problem: Optimization Model and Solution Algorithms 222222222214 Journal of Uncertan Sytem Vol.7, No.1, pp.22-35, 2013 Onlne at: www.u.org.uk Robut Capactated Faclty Locaton Problem: Optmzaton Model and Soluton Algorthm Ragheb Rahmanan *, Mohammad Sad-Mehrabad,

More information

Performance Bounds for P2P Streaming System with Transcoding

Performance Bounds for P2P Streaming System with Transcoding Appl. Math. Inf. Sc. 7, No. 6, 2477-2484 2013 2477 Appled Mathematc & Informaton Scence An Internatonal Journal http://dx.do.org/10.12785/am/070641 Performance Bound for P2P Streamng Sytem wth Trancodng

More information

A Robust Method for Calculating the Correlation Coefficient

A Robust Method for Calculating the Correlation Coefficient A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal

More information

Chapter 2 A Class of Robust Solution for Linear Bilevel Programming

Chapter 2 A Class of Robust Solution for Linear Bilevel Programming Chapter 2 A Class of Robust Soluton for Lnear Blevel Programmng Bo Lu, Bo L and Yan L Abstract Under the way of the centralzed decson-makng, the lnear b-level programmng (BLP) whose coeffcents are supposed

More information

A Quadratic Constraint Total Least-squares Algorithm for Hyperbolic Location

A Quadratic Constraint Total Least-squares Algorithm for Hyperbolic Location I. J. Communcaton, Network and Sytem Scence, 8,, 5-6 Publhed Onlne May 8 n ScRe (http://www.srpublhng.org/journal/jcn/). A Quadratc Contrant otal Leat-quare Algorthm for Hyperbolc Locaton Ka YANG, Janpng

More information

Adaptive Memory Programming for the Robust Capacitated International Sourcing Problem

Adaptive Memory Programming for the Robust Capacitated International Sourcing Problem Adaptve Memory Programmng for the Robut Capactated Internatonal Sourcng Problem Joé Lu González Velarde Centro de Stema de Manufactura, ITESM Monterrey, Méxco Lugonzal@campu.mty.tem.mx Rafael Martí Departamento

More information

Operational comparison of a new FACTS controller (RHFC) with other FACTS devices considering modified steady-state model

Operational comparison of a new FACTS controller (RHFC) with other FACTS devices considering modified steady-state model Rohanak REZAEIPOUR, Ahad KAZEMI, Mlad TAYEI 3 Iran Unverty of cence and Technology (,, 3 Operatonal comparon of a new FAT controller (RHF wth other FAT devce conderng modfed teadytate model Abtract. In

More information

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan Wnter 2008 CS567 Stochastc Lnear/Integer Programmng Guest Lecturer: Xu, Huan Class 2: More Modelng Examples 1 Capacty Expanson Capacty expanson models optmal choces of the tmng and levels of nvestments

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

m = 4 n = 9 W 1 N 1 x 1 R D 4 s x i

m = 4 n = 9 W 1 N 1 x 1 R D 4 s x i GREEDY WIRE-SIZING IS LINEAR TIME Chr C. N. Chu D. F. Wong cnchu@c.utexa.edu wong@c.utexa.edu Department of Computer Scence, Unverty of Texa at Autn, Autn, T 787. ABSTRACT In nterconnect optmzaton by wre-zng,

More information

Wind - Induced Vibration Control of Long - Span Bridges by Multiple Tuned Mass Dampers

Wind - Induced Vibration Control of Long - Span Bridges by Multiple Tuned Mass Dampers Tamkang Journal of Scence and Engneerng, Vol. 3, o., pp. -3 (000) Wnd - Induced Vbraton Control of Long - Span Brdge by Multple Tuned Ma Damper Yuh-Y Ln, Ch-Mng Cheng and Davd Sun Department of Cvl Engneerng

More information

The lower and upper bounds on Perron root of nonnegative irreducible matrices

The lower and upper bounds on Perron root of nonnegative irreducible matrices Journal of Computatonal Appled Mathematcs 217 (2008) 259 267 wwwelsevercom/locate/cam The lower upper bounds on Perron root of nonnegatve rreducble matrces Guang-Xn Huang a,, Feng Yn b,keguo a a College

More information

Hybrid State Feedback Controller Design of Networked Switched Control Systems with Packet Dropout

Hybrid State Feedback Controller Design of Networked Switched Control Systems with Packet Dropout Amercan Control Conference Marrott Waterfront, Baltmore, MD, USA June 3-July, WeB6.6 Hybrd State Feedbac Controller Degn of etwored Swtched Control Sytem wth Pacet Dropout Dan Ma, Georg M. Dmrov and Jun

More information

A Computational Method for Solving Two Point Boundary Value Problems of Order Four

A Computational Method for Solving Two Point Boundary Value Problems of Order Four Yoge Gupta et al, Int. J. Comp. Tec. Appl., Vol (5), - ISSN:9-09 A Computatonal Metod for Solvng Two Pont Boundary Value Problem of Order Four Yoge Gupta Department of Matematc Unted College of Engg and

More information

Operating conditions of a mine fan under conditions of variable resistance

Operating conditions of a mine fan under conditions of variable resistance Paper No. 11 ISMS 216 Operatng condtons of a mne fan under condtons of varable resstance Zhang Ynghua a, Chen L a, b, Huang Zhan a, *, Gao Yukun a a State Key Laboratory of Hgh-Effcent Mnng and Safety

More information

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS HCMC Unversty of Pedagogy Thong Nguyen Huu et al. A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS Thong Nguyen Huu and Hao Tran Van Department of mathematcs-nformaton,

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

Seismic Reliability Analysis and Topology Optimization of Lifeline Networks

Seismic Reliability Analysis and Topology Optimization of Lifeline Networks The 4 th World Conference on Earthquake Engneerng October 2-7, 2008, Beng, Chna Semc Relablty Analy and Topology Optmzaton of Lfelne Network ABSTRACT: Je L and We Lu 2 Profeor, Dept. of Buldng Engneerng,

More information

An Interactive Optimisation Tool for Allocation Problems

An Interactive Optimisation Tool for Allocation Problems An Interactve Optmsaton ool for Allocaton Problems Fredr Bonäs, Joam Westerlund and apo Westerlund Process Desgn Laboratory, Faculty of echnology, Åbo Aadem Unversty, uru 20500, Fnland hs paper presents

More information

j=0 s t t+1 + q t are vectors of length equal to the number of assets (c t+1 ) q t +1 + d i t+1 (1) (c t+1 ) R t+1 1= E t β u0 (c t+1 ) R u 0 (c t )

j=0 s t t+1 + q t are vectors of length equal to the number of assets (c t+1 ) q t +1 + d i t+1 (1) (c t+1 ) R t+1 1= E t β u0 (c t+1 ) R u 0 (c t ) 1 Aet Prce: overvew Euler equaton C-CAPM equty premum puzzle and rk free rate puzzle Law of One Prce / No Arbtrage Hanen-Jagannathan bound reoluton of equty premum puzzle Euler equaton agent problem X

More information

Two Approaches to Proving. Goldbach s Conjecture

Two Approaches to Proving. Goldbach s Conjecture Two Approache to Provng Goldbach Conecture By Bernard Farley Adved By Charle Parry May 3 rd 5 A Bref Introducton to Goldbach Conecture In 74 Goldbach made h mot famou contrbuton n mathematc wth the conecture

More information

Resonant FCS Predictive Control of Power Converter in Stationary Reference Frame

Resonant FCS Predictive Control of Power Converter in Stationary Reference Frame Preprnt of the 9th World Congre The Internatonal Federaton of Automatc Control Cape Town, South Afrca. Augut -9, Reonant FCS Predctve Control of Power Converter n Statonary Reference Frame Lupng Wang K

More information

This appendix presents the derivations and proofs omitted from the main text.

This appendix presents the derivations and proofs omitted from the main text. Onlne Appendx A Appendx: Omtted Dervaton and Proof Th appendx preent the dervaton and proof omtted from the man text A Omtted dervaton n Secton Mot of the analy provded n the man text Here, we formally

More information

Integrated approach in solving parallel machine scheduling and location (ScheLoc) problem

Integrated approach in solving parallel machine scheduling and location (ScheLoc) problem Internatonal Journal of Industral Engneerng Computatons 7 (2016) 573 584 Contents lsts avalable at GrowngScence Internatonal Journal of Industral Engneerng Computatons homepage: www.growngscence.com/ec

More information

A New Inverse Reliability Analysis Method Using MPP-Based Dimension Reduction Method (DRM)

A New Inverse Reliability Analysis Method Using MPP-Based Dimension Reduction Method (DRM) roceedng of the ASME 007 Internatonal Degn Engneerng Techncal Conference & Computer and Informaton n Engneerng Conference IDETC/CIE 007 September 4-7, 007, La Vega, eada, USA DETC007-35098 A ew Inere Relablty

More information

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem H.K. Pathak et. al. / (IJCSE) Internatonal Journal on Computer Scence and Engneerng Speedng up Computaton of Scalar Multplcaton n Ellptc Curve Cryptosystem H. K. Pathak Manju Sangh S.o.S n Computer scence

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

bounds compared to SB and SBB bounds as the former two have an index parameter, while the latter two

bounds compared to SB and SBB bounds as the former two have an index parameter, while the latter two 1 Queung Procee n GPS and PGPS wth LRD Traffc Input Xang Yu, Ian L-Jn Thng, Yumng Jang and Chunmng Qao Department of Computer Scence and Engneerng State Unverty of New York at Buffalo Department of Electrcal

More information

Research Article Green s Theorem for Sign Data

Research Article Green s Theorem for Sign Data Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of

More information

ELE B7 Power Systems Engineering. Power Flow- Introduction

ELE B7 Power Systems Engineering. Power Flow- Introduction ELE B7 Power Systems Engneerng Power Flow- Introducton Introducton to Load Flow Analyss The power flow s the backbone of the power system operaton, analyss and desgn. It s necessary for plannng, operaton,

More information

Calculating Jacobian coefficients of primitive constraints with respect to Euler parameters

Calculating Jacobian coefficients of primitive constraints with respect to Euler parameters Calculatng Jacoban coeffcent of prmtve contrant wth repect to Euler parameter Yong Lu, Ha-Chuan Song, Jun-Ha Yong To cte th veron: Yong Lu, Ha-Chuan Song, Jun-Ha Yong. Calculatng Jacoban coeffcent of prmtve

More information

Estimation: Part 2. Chapter GREG estimation

Estimation: Part 2. Chapter GREG estimation Chapter 9 Estmaton: Part 2 9. GREG estmaton In Chapter 8, we have seen that the regresson estmator s an effcent estmator when there s a lnear relatonshp between y and x. In ths chapter, we generalzed the

More information

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem Appled Mathematcal Scences Vol 5 0 no 65 3 33 Interactve B-Level Mult-Objectve Integer Non-lnear Programmng Problem O E Emam Department of Informaton Systems aculty of Computer Scence and nformaton Helwan

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

A Two-Stage Modeling and Solution Framework for Multisite Midterm Planning under Demand Uncertainty

A Two-Stage Modeling and Solution Framework for Multisite Midterm Planning under Demand Uncertainty Ind. Eng. Chem. Re. 2000, 39, 3799-3813 3799 A Two-Stage Modelng and Soluton Framework for Multte Mdterm Plannng under Demand Uncertanty Anhuman Gupta and Cota D. Marana* Department of Chemcal Engneerng,

More information

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010 Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton

More information

M. Mechee, 1,2 N. Senu, 3 F. Ismail, 3 B. Nikouravan, 4 and Z. Siri Introduction

M. Mechee, 1,2 N. Senu, 3 F. Ismail, 3 B. Nikouravan, 4 and Z. Siri Introduction Hndaw Publhng Corporaton Mathematcal Problem n Engneerng Volume 23, Artcle ID 795397, 7 page http://dx.do.org/.55/23/795397 Reearch Artcle A Three-Stage Ffth-Order Runge-Kutta Method for Drectly Solvng

More information

Application of Nonbinary LDPC Codes for Communication over Fading Channels Using Higher Order Modulations

Application of Nonbinary LDPC Codes for Communication over Fading Channels Using Higher Order Modulations Applcaton of Nonbnary LDPC Codes for Communcaton over Fadng Channels Usng Hgher Order Modulatons Rong-Hu Peng and Rong-Rong Chen Department of Electrcal and Computer Engneerng Unversty of Utah Ths work

More information

Deep Reinforcement Learning with Experience Replay Based on SARSA

Deep Reinforcement Learning with Experience Replay Based on SARSA Deep Renforcement Learnng wth Experence Replay Baed on SARSA Dongbn Zhao, Hatao Wang, Kun Shao and Yuanheng Zhu Key Laboratory of Management and Control for Complex Sytem Inttute of Automaton Chnee Academy

More information

Assortment Optimization under MNL

Assortment Optimization under MNL Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.

More information

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 ISSN 1549-644 1 Scence Publcatons do:1.844/jmssp.1.4.8 Publshed Onlne 9 (1) 1 (http://www.thescpub.com/jmss.toc) A MODIFIED METHOD FOR SOLVING SYSTEM OF

More information

Computing Correlated Equilibria in Multi-Player Games

Computing Correlated Equilibria in Multi-Player Games Computng Correlated Equlbra n Mult-Player Games Chrstos H. Papadmtrou Presented by Zhanxang Huang December 7th, 2005 1 The Author Dr. Chrstos H. Papadmtrou CS professor at UC Berkley (taught at Harvard,

More information

PROBABILITY-CONSISTENT SCENARIO EARTHQUAKE AND ITS APPLICATION IN ESTIMATION OF GROUND MOTIONS

PROBABILITY-CONSISTENT SCENARIO EARTHQUAKE AND ITS APPLICATION IN ESTIMATION OF GROUND MOTIONS PROBABILITY-COSISTET SCEARIO EARTHQUAKE AD ITS APPLICATIO I ESTIATIO OF GROUD OTIOS Q-feng LUO SUARY Th paper preent a new defnton of probablty-content cenaro earthquae PCSE and an evaluaton method of

More information

Stability Analysis of Inverter for Renewable Energy

Stability Analysis of Inverter for Renewable Energy vance n Power an Energy Sytem Stablty naly of nverter for Renewable Energy TOORU MOR, JUNCH R Electrcal Engneerng an Electronc ogakun Unverty 1-24-2 Nh-hnjuku, Shnjuku-ku, Tokyo 163-8677 JPN cm1134@n.kogakun.ac.jp

More information

BOUNDARY ELEMENT METHODS FOR VIBRATION PROBLEMS. Ashok D. Belegundu Professor of Mechanical Engineering Penn State University

BOUNDARY ELEMENT METHODS FOR VIBRATION PROBLEMS. Ashok D. Belegundu Professor of Mechanical Engineering Penn State University BOUNDARY ELEMENT METHODS FOR VIBRATION PROBLEMS by Aho D. Belegundu Profeor of Mechancal Engneerng Penn State Unverty ahobelegundu@yahoo.com ASEE Fello, Summer 3 Colleague at NASA Goddard: Danel S. Kaufman

More information

Scattering of two identical particles in the center-of. of-mass frame. (b)

Scattering of two identical particles in the center-of. of-mass frame. (b) Lecture # November 5 Scatterng of two dentcal partcle Relatvtc Quantum Mechanc: The Klen-Gordon equaton Interpretaton of the Klen-Gordon equaton The Drac equaton Drac repreentaton for the matrce α and

More information

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k) ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of

More information

STOCHASTIC BEHAVIOUR OF COMMUNICATION SUBSYSTEM OF COMMUNICATION SATELLITE

STOCHASTIC BEHAVIOUR OF COMMUNICATION SUBSYSTEM OF COMMUNICATION SATELLITE IJS 4 () July Sharma & al ehavour of Subytem of ommuncaton Satellte SOHSI HVIOU O OMMUNIION SUSYSM O OMMUNIION SLLI SK Mttal eepankar Sharma & Neelam Sharma 3 S he author n th paper have dcued the tochatc

More information

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there

More information

Power law and dimension of the maximum value for belief distribution with the max Deng entropy

Power law and dimension of the maximum value for belief distribution with the max Deng entropy Power law and dmenson of the maxmum value for belef dstrbuton wth the max Deng entropy Bngy Kang a, a College of Informaton Engneerng, Northwest A&F Unversty, Yanglng, Shaanx, 712100, Chna. Abstract Deng

More information

MODELLING OF TRANSIENT HEAT TRANSPORT IN TWO-LAYERED CRYSTALLINE SOLID FILMS USING THE INTERVAL LATTICE BOLTZMANN METHOD

MODELLING OF TRANSIENT HEAT TRANSPORT IN TWO-LAYERED CRYSTALLINE SOLID FILMS USING THE INTERVAL LATTICE BOLTZMANN METHOD Journal o Appled Mathematc and Computatonal Mechanc 7, 6(4), 57-65 www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.4.6 e-issn 353-588 MODELLING OF TRANSIENT HEAT TRANSPORT IN TWO-LAYERED CRYSTALLINE SOLID

More information

Numerical Methods for Solving Turbulent Flows by Using Parallel Technologies

Numerical Methods for Solving Turbulent Flows by Using Parallel Technologies Journal of Computer and Communcaton, 0,, -5 do:0.46/cc.0.00 Publhed Onlne February 0 (http://www.crp.org/ournal/cc) Numercal Method for Solvng urbulent Flow by Ung Parallel echnologe Albek Iakhov Department

More information

SCENARIO SELECTION PROBLEM IN NORDIC POWER MARKETS. Juha Ojala. Systems Analysis Laboratory, Helsinki University of Technology

SCENARIO SELECTION PROBLEM IN NORDIC POWER MARKETS. Juha Ojala. Systems Analysis Laboratory, Helsinki University of Technology 2.2.2003 SCENARIO SELECTION ROBLEM IN NORDIC OER MARKETS MAT-2.08 INDEENDENT RESEARCH ROJECT IN ALIED MATHEMATICS Juha Oala oala@cc.hut.f Sytem Analy Laboratory, Heln Unverty of Technology ortfolo Management

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information