On the Convexity of the System Loss Function
|
|
- Ellen Booker
- 5 years ago
- Views:
Transcription
1 O the Covety of the System Loss Fucto Sebasta de la orre, Member, IEEE ad Fracsco D. Galaa, Fellow, IEEE Abstract We show that the system fucto a power etwork s bouded below by ay umber of supportg hyperplaes the space of geeralzed jectos. A supportg hyperplae s defed by the lear aylor seres epaso of the system fucto aroud a gve operatg pot. he supportg hyper-plae property s vald provded that the epaso pot s suffcetly ear the flat-voltage profle (F). We have assessed epermetally the rage of valdty of ths assumpto, called the rage of hyper-plae support (RHS), showg that for typcal etworks, the RHS s broad, partcularly whe the bus voltages are cotrolled ear per ut. he supportg hyper-plae model was also tested as part of a ecoomc dspatch wth trasmsso es to demostrate that ths lear model provdes the same results as whe the es are treated as a o-lear fucto. Ide erms System, load flow, lear appromato, supportg hyper-plaes, covety, formulae. I. OMECLAURE Defto: Flat oltage rofle (F) s the operatg codto uder whch all comple bus voltages are equal to some arbtrary level, e + jf rectagular coordates or δ polar coordates, where e cos( ) f s ( ) = δ ad = δ. Wthout of geeralty, we assume that bus s the referece wth δ = or, equvaletly, wth f =. hs the mples that δ = ad that f =. arameters: umber of etwork buses; umber of load flow equatos ad ukows = -; G Real part of the etwork admttace matr; G ' Sub-matr of G wthout the -th row ad colum; L Block dagoal matr formed wth G ad G' ; B Imagary part of the etwork admttace matr; B ser Compoet of B due to le seres admttaces; B sh Compoet of B due to shut admttaces; ector of s of dmeso ; ector or matr of s; -dmesoal vector of s wth at posto ; hs work was supported by the atural Sceces ad Egeerg Research Coucl (SERC), Caada, by the Fods ature et techologes, Quebec, ad by the Europea Uo project FEDER-CICY FD ad grat FU-A99 of the Mstero de Educacó, Cultura y Deporte of Spa. F. D. Galaa s wth the Departmet of Electrcal ad Computg Egeerg at McGll Uversty, Motreal, Québec, Caada. S. de la orre, was wth the E..S.I. Idustrales, Uversdad de Castlla- La Macha, Cudad Real, Spa, ad s curretly wth the E..S.I. elecomucaoes, Uversdad de Malaga, Spa her e-mals are: galaa@ece.mcgll.ca ad storre@uma.es. Ω Set of odes coected to ode. arables: ector of comple odal voltages; e Real part of ; f Imagary part of ; f ' f ecludg the referece bus ; ector of magtudes squared of ; ector of magtudes of ; δ ector of phase agles of ecludg the referece bus ; -dmesoal vector of rectagular coordate voltage e compoets, ' ; f alue of uder F; e Real part of uder F; Real power jecto at bus ; Q Reactve power jecto at bus ; ector of real power jectos at all buses; ' ecludg the slack bus; Q ector of reactve power jectos at all buses; S Set of buses wth specfed real power jecto;.e., all buses ecept for the slack. Set of buses wth specfed reactve power jecto; S Q S Set of buses where the voltage magtude s specfed; z ector of geeralzed jectos of dmeso comprsed of all S, all Q S Q ad all S ; z() Fucto of characterzg the geeralzed jectos; System ; I ( ) Fucto of characterzg the system. II. IRODUCIO problems such as ut commtmet, resource schedulg, or etwork epaso, the power etwork s ofte modeled a appromate maer or s altogether eglected. hs s doe order to reduce the umber of varables ad the overall computatoal complety, but also to make the problem lear thereby permttg the use of powerful med teger lear programmg tools []. oday, such tools ca relably solve problems wth very large umbers of costrats ad varables, both teger ad cotuous, whle o-lear med-teger programmg tools are stll urelable ad wth lmted applcatos. he smplest etwork models completely gore the trasmsso ad descrbe the etwork as a sgle ode. he et level of complety models the system trasmsso,, the power balace equato as a appromate eplct polyomal fucto of the geerato levels, g,
2 m g = d + ( g) () = Loss appromato formulae were etesvely studed as early as the 96 s, begg wth the B-coeffcets approach [] whch the system fucto was epressed as a costat coeffcet quadratc fucto of the power geerato levels, g. I the early 98 s, takg advatage of specfc aalytc propertes of the load flow equatos rectagular coordates [3], more eact eplct appromato formulae up to thrd order were derved based o aylor seres epasos of terms of z aroud a specfed load flow operatg pot [4]. A smlar result to [4] was publshed recetly wth applcatos to voltage stablty [5]. Other authors have eamed eplct appromatos based o polyomals the geerato levels [6, 7]. It s also possble to model the of every trasmsso le a power etwork by a quadratc appromato of the correspodg le phase agle dfferece. hs leads to a varato of the DC load flow model of the form, = B' δ [ ] = B δ + q( δ) f f f f f whch cludes the le power flows, f, ther lmts () f, ad the appromato of the dvdual le real trasmsso es, q(δ), [8, 9]. otwthstadg the avalablty of the above metoed appromato formulae, whe solvg problems that do ot volve teger decsos, such as a optmal power flow, our eperece t s questoable whether such appromatos offer ay sgfcat advatage over a olear optmzato method that accouts for the full olear load flow model. I problems volvg teger decsos, such as ut commtmet, oe dsadvatage of the more accurate quadratc or hgher order appromato formulae s that order to make use of med-teger lear programmg solvers, t s ecessary to further decompose each o-lear term to pece-wse lear compoets. hs step requres the troducto of ew cotuous ad possbly teger varables, therefore addg to the overall modelg ad computatoal complety. O the other had, the use of lear appromato formulae, although readly compatble wth med-teger lear programmg tools, s subject to greater appromato accuraces. Oe drawback of all appromato formulae, whether lear or hgher order, s that the errors troduced by the appromato comparso wth the eact model caot be easly quatfed or estmated. he ma cotrbuto of ths paper s a theoretcal result that had bee prevously hypotheszed [], but ot rgorously prove or tested, amely to demostrate that, uder certa relatvely weak assumptos, the o-lear system fucto,, s cove the geeralzed power flow jectos, z. As such, a frst order aylor seres epaso ca be show to defe a lower boud appromato of the system lear z, kow as a supportg hyper-plae. Such represetatos are ot appromatos the usual learzato sese, sce they also costtute a lower boud o the o-lear fucto. hs property leads to two mportat results: () As each such equalty s a ecessary codto, ew hyper-plaes ca be added wthout valdatg the estg oes; () he error betwee the o-lear fucto ad the descrbg supportg hyper-plaes s greater tha or equal to zero. hus, the olear system behavour ca be approached wth arbtrary accuracy by creasg the umber of lear ecessary codtos. It s also mportat to ote that the learzed model based o the well-kow cremetal trasmsso (IL) coeffcets s a specal case of a supportg hyperplae. Aother cotrbuto of ths paper s to quatfy the rage of valdty of the covety property of the fucto sde of whch the system es ca be modeled by a set of lear equaltes. umercal results o systems up to 8 buses suggest that the covety property s vald over a wde rage of operatg codtos, cludg typcal operatg pots. he system boudg theorem has also bee appled to solve a y ecoomc dspatch problem to llustrate the accuracy of the hyper-plae appromato. III. SYSEM LOSS FUCIO SUORIG HYER-LAES A. o-lear Load Flow Model I load flow models the decso varables ad goverg equatos are as follows, = ( ); S Q = Q ( ); S Q = ( ); S = = ( ) ( ) = = B. Geeralzed Load Flow Ijectos he vector of geeralzed load flow jectos, z R (3), s defed by the set of, Q ad of equato (3). As show Apped A, Result, the relato betwee z ad s deoted by, z= z( ) J( ) (4) C. System Loss Sestvty ector he system sestvty vector aroud a gve operatg pot, deoted by β, s defed by the sestvtes of the system wth respect to the geeralzed load flow jectos, z, that s, ( z ) z ( ) ( ) ( ) β = = z (5) = J ( ) L
3 3 where we have used the deftos of z() equato (4) ad of () Result 5 of Apped A. ote that the well-kow cremetal trasmsso (IL) coeffcets are the compoets of β correspodg to the real power jectos. he sestvty vector, β, characterzes the frst order aylor seres epaso of the fucto aroud the operatg pot, that s, ( ) = + β z z +ε ( ) ( ) where the quatty ε represets the error ths lear appromato due to hgher order terms. Oe of the propertes of the rectagular load flow formulato derved Apped B s that, ( ) = β z ( ) (7) whch whe substtuted to (6) reduces the frst order epaso of the system to, = β z + ε (8) Apped B also shows that whe s suffcetly ear the F, the appromato error ε s o-egatve for ay feasble z, that s, for ay z for whch there ests a correspodg load flow soluto,. hese results allow us to formally eucate the followg theorem. D. System Loss Boudg heorem Gve ay operatg pot suffcetly ear a flatvoltage-profle such that the load flow Jacoba matr, J ( ), s o-sgular, the, for ay feasble vector of geeralzed load flow jectos, z, the correspodg system satsfes, (6) β z (9) he proof of ths theorem s detaled Apped B. E. Corollary to the System Loss Boudg heorem Uder the codtos of the system boudg theorem, for ay feasble z ot collear wth z ( ), > β z () whle for ay feasble z collear wth z ( ), = β z () he proof of ths corollary follows drectly from Result 9 Apped A. he above lear equato () defes what s called a supportg hyper-plae of the ukow o-lear fucto. Sce power systems ormally operate relatvely ear the F, the promty codto requrg that the epaso pot be suffcetly ear the F s geerally ot a severe lmtato. Epermetal evdece preseted the results secto backs up ths statemet. F. Lear System Loss Model Based o Supportg Hyper- laes Sce (9) s a ecessary codto, ay umber of supportg hyper-plaes based o dfferet epaso pots ca be merged to a more restrctve costrat set. Let SH be a set of β ' s of whch each elemet correspods to a supportg hyper-plae of the system fucto. he, the power balace equato wth es ca be appromated by the followg lear model, = = β z; β SH () I addto, operatoal costrats such as m ma retag ts lear ature. ca easly be added to the model whle stll I. UMERICAL RESULS All quattes ths secto are per ut ecept for the bus voltage phase agles whch are degrees. A. Eamples of Supportg ad o-supportg Hyperplaes Cosder a 5-bus etwork wth the followg le data, ABLE I: ES EWORK LIE DAA; MA, K BASES. From Bus o Bus Seres R Seres X Shut B he F codto polar coordates s gve by δ = [,,,] ad = [,,,,], where the referece agle s δ (5) =. I terms of geeralzed load flow jectos, the F s descrbed by ' = [ ] = [ ],,, ad,,,,. ables II, III ad I below descrbe the supportg hyperplae correspodg to a specfed operatg pot epressed polar coordates. he symbols β ad β these tables deote respectvely the compoets of β correspodg to the geeralzed jectos, ' ad.
4 4 ABLE II: EXAMLE A: SUORIG HYER-LAE. 5-BUS SYSEM ode δ ( ') β ABLE III: EXAMLE A: SUORIG HYER-LAE. 5-BUS SYSEM ode ( ) β ABLE I: EXAMLE A: SUORIG HYER-LAE. 5-BUS SYSEM Egevalue umber λ (L) λ ( L H( β)) I ths case, the epaso pot,, show able II, epressed polar coordates, has phase agles that are ot all close to zero, oe of them beg equal to 43 degrees. hs suggests that the requremet that the epaso pot be suffcetly ear the F s ot very strget. Moreover, as called for by the theory, able I shows that the egevalues of L are all postve but oe. ote, addto, that the egevalues of the error matr, L H( β ), are all close to those of L, ad that the zero egevalue of L remas at zero for the error matr. able, I ad II show a eample of a operatg pot that does ot correspod to a supportg hyper-plae. Here, two of the egevalues of [ L H( β ) ] are egatve. A lear epaso about ths operatg pot would therefore ot be a lower boud for the system over all z. ote that ths o-supportg case, the epaso pot s relatvely far from the F as evdeced from the values of the phase agles. ABLE : EXAMLE B: O-SUORIG HYER-LAE. 5-BUS SYSEM ode δ ( ') β ABLE I: EXAMLE B: O-SUORIG HYER-LAE. 5-BUS SYSEM ode ( ) β ABLE II: EXAMLE B: O-SUORIG HYER-LAE. 5-BUS SYSEM Egevalue umber λ (L) λ ( L H( β)) A etwork wth data based o the IEEE-8-bus test system wth all bus voltages held at per ut s cosdered ow. Due to the large data set, oly some values are provded. able III shows the operatg pot polar coordates, the geeralzed jectos, ad the values of β at two buses. able IX shows a subset cosstg of the smallest ad largest egevalues of L ad of the error matr, all of whch are postve ecept for oe whch s zero. ABLE III: EXAMLE C: SUORIG HYER-LAE. 8-BUS SYSEM z ( ) β ode o o =.975 = = 4.43 = β β 4 4 β β 9 9 =.633 =.44 =.487 =.679 ABLE IX: EXAMLE C: SUORIG HYER-LAE. 8-BUS SYSEM λ (L) λ ( L H( β)) B. Cotuty ad Rage of Hyper-plae Support he rage of hyper-plae support (RHS) descrbes the set of epaso pots vald. Formally, the RHS s a set of for whch the covety property (9) s defed by,
5 5 / β z; β = ( ) ; ( ) RHS J L L H β = (3) = ( ); z= z( ); R Sce the load flow Jacoba,, s o-sgular whe s at or ear, the RHS s a set that cludes ad eteds cotuously outward from boudares are reached whe J( ) L H(β) becomes defte; () J ) all drectos. Its satsfes ether of: () becomes sgular. et, some epermetal tests to estmate the sze of the RHS are coducted. ( C. umercal Estmates of the Rage of Hyper-plae Support I order to quatfy the RHS, a eghborhood aroud the F codto s defed, j ma le δ δ δ ;, j Ω (4) hs codto lmts the agle dffereces across all les. able X shows the results of varous epermets o the 5-bus ma δ le etwork for dfferet values of, each based o, values of the epaso pot,, pcked uformly radomly wth the set defed by (4). As able X shows, the rage of estece of supportg hyper-plaes s qute broad, ecompassg bus agles betwee ±9 degrees ad le agle dffereces up to 76.5 degrees; the latter value beg the most sgfcat. Outsde the rego of support, a small but ozero percetage of the radomly geerated hyper-plaes does ot support the fucto. Moreover, the farther from that rage that a hyper-plae les, the hgher the chaces that t s o-supportg. ABLE X: RAGE OF HYER-LAE SUOR. 5-BUS SYSEM Rato of o-supportg to Supportg hyper-plaes ma δ le /, 76.5 º 3/, 8º /, 9 º he results for the 8-bus etwork wth radomly geerated epaso pots are show able XI. ABLE XI: RAGE OF HYER-LAE SUOR. 8-BUS SYSEM Rato of o-supportg to Supportg hyper-plaes ma δ le /, 68.4 º 7/, 7 º 73/, 8 º hese results dcate that f the absolute le agle dffereces are less tha 68.4º, the, the correspodg hyper-plae s lkely to be supportg. ote that epermets dcate that the rage of hyper-plae support remas essetally uaffected f the bus voltages are kept wth plus or mus 5% of omal.. ALICAIO: ECOOMIC DISACH WIH LOSSES he geeral theory preseted ths paper s sutable for use may dfferet practcal problems where es are a mportat ssue; eamples beg ut commtmet, ecoomcdspatch ad hydro-thermal coordato. ote however that the supportg hyper-plae approach s ot teded to be used to solve problems wth costrats o depedet varables such as le flows, load bus voltage magtudes, or reactve geerato levels. Such problems should be solved wth a optmal power flow For the sake of llustrato, the proposed geeral theory s tested solvg a y ecoomc dspatch problem; the system used s the prevously preseted 8-bus power system; wth 49 geeratg buses, 58 load buses ad 79 les. A. Solvg the ecoomc dspatch. ote that the ecoomc dspatch problem, whe modeled wth supportg hyper-plaes, cossts of the objectve fucto, the bouds for the geerato, oe equalty costrat per hyper-plae ad a addtoal equalty costrat statg that the summato of geerato must equal the summato of demad plus es; thus, o eplct accout for the etwork s eeded. he y ecoomc dspatch may be solved eactly usg hyper-plaes as follows (Algorthm ): Step. Italze a workg set of hyper-plaes wth the trval equalty costrat, that s,. Step. Solve the ecoomc-dspatch problem usg the workg set of hyper-plaes for the geerator outputs. Step. Solve a power flow problem whose put data are the geerato levels of all uts from Step, ecept for the slack. he soluto provdes the values of ad of the slack phase agles of all buses. Step 3. Wth the data obtaed Step compute vector β usg (5). hs vector defes a ew hyper-plae. Step 4. If the dfferece betwee the values obtaed for steps ad s uder a certa threshold, SO; slack otherwse, add the hyper-plae derved Step 3 to the set of workg hyper-plaes ad retur to Step. A alteratve algorthm that s better suted for o-le operato s ow preseted (Algorthm ): Step : Compute a good set of hyper-plaes off-le. Step : Whe o-le, for ay actual value of the vector of demads, solve a ecoomc dspatch usg the equalty costrats obtaed Step. B. Estmatg the error of the method It s clear that Algorthm, the umber of hyper-plaes calculated off-le, H, s the most mportat parameter determg the cost error. hus, a procedure s ow preseted for estmatg the mmum cost error terms of the umber of hyper-plaes used to solve a y ecoomc dspatch, error( H ). (Algorthm 3): he cocept of good depeds o the precso desred; the more hyperplaes cosdered Step, the more accurate the soluto of Step.
6 6 Step : Usg ay sutable algorthm, such as Algorthm above, solve for the eact soluto of a large umber,, of y ecoomc dspatch problems wth dfferet demad levels, ad calculate for each of them the optmal cost: Cost ( ) :.... Step : Compute H supportg hyper-plaes for a set of H operatg codtos dfferet from those foud Step. Step : Solve each of the problems from Step usg oly the H equalty costrats from Step usg Algorthm to obta the appromate costs Cost( ) :.... Comparg these costs wth the optmal costs that were obtaed Step, compute the correspodg 3 errors. Defe ε ) = Cost ( ) Cost ( ) as the error ( obtaed whe tryg to appromate the -th problem wth H hyper-plaes. Defe the mamum of these errors as a estmate of the mamum total error that wll be obtaed wheever H hyper-plaes are used to solve the ecoomc dspatch problem wth Algorthm. I Algorthm 3, ote that H s the umber of hyper-plaes used to create a model. Also, ote that s the umber of problems used to assess the accuracy of that model; the eact soluto for those problems s calculated Step. Fgure shows the mamum total cost error as a percetage of the eact o-lear cost terms of the umber of hyper-plaes used Algorthm. ote that there are four plots Fgure for two dfferet values of ( ad 4) ad two dfferet values of w (% ad 3%). arameter w dcates how ear to the base-case demad were the H radom demads. For stace, w=3% meas that, the radomly geerated demad values le sde a terval of ±3% aroud ther orgal value. Cost error [%],8,7,6,5,4,3,,, =4, w =3% =, w =3% =4, w =% =, w =% Fgure. Cost error as a fucto of the umber of hyper-plaes used Algorthm ote that the effect of addg or removg hyper-plaes as llustrated by Fgure of the paper s ot very sgfcat o the optmum dspatch cost, for eample, as few as hyperplaes gve a error of aroud.3%. I. COCLUSIOS We have show that the system fucto a power etwork s bouded below by ay umber of supportg hyper-plaes the space of geeralzed jectos. A supportg hyper-plae s defed by the lear aylor seres epaso of the system fucto wth respect to the geeralzed jectos aroud a gve operatg pot. hs supportg hyper-plae property s vald over a cotuous rage of epaso pots provded that these are suffcetly ear the flat-voltage profle (F). he requred degree of promty betwee a operatg pot defg a supportg hyper-plae ad the F, called the rage of hyper-plae support (RHS), was assessed epermetally. For typcal etworks, the RHS was foud ot to be very restrctve, partcularly for etworks where the bus voltages are closely cotrolled ear per ut. he supportg hyper-plae model was also tested as part of a ecoomc dspatch wth trasmsso es to demostrate that ths lear model provdes the same results as whe the es are treated as a o-lear fucto. Future research should look at tegratg the supportg hyper-plae model to problems such as ut commtmet wth trasmsso es. he system boudg theorem should also be eamed wth etworks cotag FACS devces or phase-shftg trasformers. I addto, t would be terestg to eted these results to load flow models wth dstrbuted slack geerato. II. AEDIX A For completeess, most of the relevat propertes of the load flow equatos rectagular coordates are preseted ad prove ths apped. For these ad other terestg propertes of the quadratc formulato of the load flow problem, see also [], [] ad [3]. For all the results that follow, we make the realstc assumpto that all trasmsso les have o-zero seres resstace ad zero shut coductace. Result : he system a power trasmsso etwork s postve or zero. roof: hs s self-evdet from eergy coservato ad the fact that the trasmsso etwork s a passve RLC etwork. Result : he system s zero f ad oly f all the seres brach currets are zero. roof: Let I ser j be the seres curret magtude of the le coectg buses ad j wth o-zero resstace R j. he correspodg trasmsso compoet, I ser jrj, s the zero f ad oly f the seres curret magtude s zero. Uder the assumpto that the shut braches are less, the system wll also be zero f ad oly f all the seres brach currets are zero. Result 3: All the seres brach currets are zero f ad oly f the comple bus voltages at all buses are equal. A sutable mmum umber for the IEEE 8-bus system wll be. 3 hs wll always be a postve umber, because the supportg hyperplaes are lower bouds of the actual fucto.
7 7 roof: he seres brach curret through the le coectg buses ad j wth o-zero seres mpedace, R j +jx j, s I ser j j = R + jx j j, whch s zero f ad oly f = j. Result 4: he system s zero f ad oly f the etwork operates at F. roof: hs follows drectly from the defto of F, Results ad 3, ad the zero shut coductace assumpto. Result 5: I rectagular coordates of the comple bus voltages, the system ca be epressed by the followg pure quadratc form, = L ( ) (A) where, = G L (A) G' roof: Epressg the comple ode voltages rectagular coordates, = e+ j f, the total real power cosumed by the etwork, whch defes the system es, s gve by, * * = =R Y jj = = j= * { Y } * {( e jf) ( G jb)( e jf) } =R =R If G s symmetrc, the = ege+ fgf ( ) = ege+ f' Gf ' ' (A3) (A4) However, f G s ot symmetrc, as may occur etworks wth phase-shftg trasformers, the equato (A4), G gets replaced, wthout of geeralty by the symmetrc matr ( G+ G )/. hs s ecessary sce, f G s ot symmetrc, ts egevalues are ot related to those of ( G+ G )/. he latter s always a symmetrc matr whose egevalues are real ad o-egatve. We therefore assume throughout ths apped that G s symmetrc. Recall from the omeclature that, wthout of geeralty, f s here set to zero, thus obtag the secod form equato (A4). Result 6: he symmetrc matr G s postve sem-defte. roof: From Result, for all f ad e. I partcular, by lettg f =, t follows from equato (A) Result 5 that = ege for ay e, whch defes G as postve semdefte. From matr theory, the postve sem-defte codto also mples that all the egevalues of G are real ad o-egatve. Result 7: he coductace matr G has oe ad oly oe zero egevalue wth egevector e e where e s a arbtrary costat. he remag - egevalues must therefore be real ad postve. = roof: o prove that G has oe o-zero egevalue, let the etwork operate at F, that s, = e. From Result 3, all brach currets must be ull, whch mples that the compoet of the et bus curret jectos due to brach currets must also be ull, other words, Y = ( G+ jb ) e = (A5) where ser = + j ser ser ser Y G B s the compoet of the etwork admttace matr due to seres elemets. ote that the shut curret vector compoets, Y = ( jb ) s ot zero at F but uder the assumpto of zero shut coductaces ths compoet does ot cotrbute to the system. Codto (A5) that the seres curret jectos are ull uder F mples that G ( e ) =, whch proves that G has a zero egevalue whose egevector s e. I order to prove that G has oly oe zero egevalue, suppose that there ests aother o-zero vector ê R satsfyg the zero egevalue codto, G ê =, ad ot of the F form e. Sce = e ˆ s a vector of voltages that does ot correspod to the F, accordg to Result 4, the system must be strctly postve, that s, >. A cotradcto the arses sce (A4) states that ˆ ˆ = ege =. Cosequetly, there caot est a vector e ˆ wth the above property ad the oly possble egevector correspodg to the zero egevalue s e. Result 8: he symmetrc matr G ' s postve defte. roof: From Result 6, G s postve sem-defte. hus for all f, f Gf. Moreover, from Result 7, f Gf = for ozero f oly f f s proportoal to. ow, f f =, ( f') G' f' = f Gf. However, wth f =, t s mpossble for f to be proportoal to. hus, wth f =, ( f') G' f' = f Gf > for all f ' ecept for the trval case whe f' =, whch s the codto for postve defteess of G '. Result 9: he matr L defed Result 5 s postve semdefte wth a sgle zero egevalue, whose egevector correspods to the F. roof: From the propertes of block-dagoal matrces, the set of egevalues of L s the uo of the respectve sets of egevalues of G ad G '. As prove Results 7 ad 8, the combed egevalues of these two matrces are all postve ecept for oe whch s zero whose egevector correspods to the F. Result : he vectors of real ad reactve power jectos, ad Q, as well as the vector of bus voltage magtudes,, are quadratc fuctos of the real ad magary comple voltage compoets,. roof: he vector of comple power jectos s, sh sh * { } I dag{ }{ Y} * S = + j Q = dag = (A6) whch rectagular form becomes, S = dag e+ j f G+ j B e+ j f } (A7) { }( )( ) { *
8 8 Separatg S to ts real ad magary parts, ad Q, gves the followg quadratc equatos, = dag{}( e Ge Bf ) + dag{}( f Gf + Be) (A8) Q = dag f Ge Bf dag e Gf + Be {}( ) {}( ) Smlarly, sce = e + f, the vector of voltage magtudes squared ca also be epressed as a quadratc, = dag e e+ dag f f (A9) { } { } Result : Ay scalar lear combato of elemets of Q, ad ca be epressed the form H where H s a real symmetrc matr depedet oly o the etwork parameters. roof: From the defto of the -dmesoal vector the omeclature ad from equato (A8) Result the real power jecto at bus s epressed as, = ( ) ˆ {}( ) dag{}( ) { }( ) dag{ }( ) = ( ) dag e Ge Bf + f Gf + Be = e dag Ge Bf + f Gf + Be Defg ow the symmetrc matr, H = { } G+ G { } { } B+ B { } dag{ } dag{ } dag{ } dag{ } dag dag dag dag B + B G+ G t follows from (A) that ca be epressed as, e ˆ e = H f f (A) (A) (A) ow, defg H as H ˆ wthout the -th row ad colum, = H (A3) Smlarly, for the reactve power jectos we obta, Q Q = H (A4) Q where H s foud by removg the -th row ad colum from, ˆ Q dag{ } B B dag{ } dag{ } G+ G dag{ } H = G dag{ } + dag{ } G dag{ } B B dag{ } (A5) Fally, for bus voltage magtudes squared we obta, where from, H = H (A6) s foud by removg the -th colum ad row dag{ } = dag{ } ˆ H (A7) ow, from equatos (A3), (A4) ad (A6), t follows that ay dvdual geeralzed power flow jecto, z, ca therefore be epressed as a pure quadratc of the form, z = [ H ] (A8) hus, α z = Hα ( ) (A9) = where we defed, H( α) = αh = (A) Result : he load flow problem wth specfed geeralzed jectos, z ; =,,, ca be epressed as a set of equatos of the form, z= J ( ) (A) where J ( ) s a by matr defed by, H J ( ) = M (A) H roof: hs result follows mmedately from equato (A9) Result. Result 3: he matr J() s the Jacoba matr of the load flow equatos rectagular coordates. roof: he load flow Jacoba matr s defed as the sestvty of z wth respect to, whch from equatos (A) ad (A) Result gves, H M H z H = = M = J ( ) (A3) H Result 4: he Jacoba matr J() s lear. roof: hs result follows drectly from (A3). Result 5: For ay, α J( ) = H(α). roof: From Result, Hα ( ) = α = α H = H = H = = α M α J( ) H III. AEDIX B: ROOF OF SYSEM LOSS BOUDIG HEOREM (A4) As show Result 5 of Apped A, the system fucto ca be wrtte as a pure quadratc form, = L (B) Furthermore, from Result 9 Apped A, L has oly oe zero egevalue whose correspodg egevector s the flat voltage profle. We eame ow the behavor of the load flow equatos ad the system fucto terms of devatos from a arbtrary operatg pot,, so that = +. Recallg the epresso for the load flow equatos Result ad the fact that the load flow Jacoba matr s lear from
9 9 Result 4, we have, z = J( + )( + ) (B) = J( )( ) + J( ) + J( ) Smlarly, the system equato (B) ca be rewrtte as, = ( ) L + ( ) L + L (B3) Sce the load flow Jacoba at the operatg pot,, s o-sgular, from (B), = [ J( )] z J( )( ) J( )( ) (B4) Substtutg (B4) to the lear term of (B3) gves, = ( ) L + [ ] L + ( ) L [ J( )] z J( )( ) J( )( ) (B5) [ ( )] z J( )( ) + [ ] L = ( ) L J From the defto of β (5), ad takg ts traspose, the eact system equato (B5) ca ow be epressed as, = β z J( )( ) + [ ] L (B6) = β z β J( )( ) + [ ] L ow, usg Result 5 wth α takg the partcular value of β, t follows from (B6) that, = β z+ [ ] [ L H(β) ] (B7) Comparg (B7) wth (8), we see that the appromato error, ε, s gve by the quadratc form, ε= [ ] [ L H( β) ] (B8) o complete the proof of the System Loss Boudg heorem, we ow show that f F, the real symmetrc matr les suffcetly ear a L H(β) s postve semdefte. hs wll mply that for all R, [ L H( β) ] satsfyg respectvely (B) ad (4), o prove that L H(β) or equvaletly that for all β z ad z (B9) s postve sem-defte frst recall from Result 9 that L =. hus, f, the L L =. Smlarly, gve that the load flow Jacoba s o-sgular at, the J ( ) L =. ow, from equato (A) Result, f β the H( β). β = [ J( ] ) L Sce L ad L H( β ) are arbtrarly close real symmetrc matrces, ther respectve sets of egevalues are real ad arbtrarly close. From Result 9, the egevalues of L are strctly postve ecept for λ ( L ) =. hs meas that all the egevalues of L H( β ) are guarateed to be postve ecept possbly λ ( L H( β )) whch although arbtrarly close to λ ( ) L =, may have shfted to a egatve value. o dscard ths opto, recall from equato (5) ad from Result 5 Apped A that, he, [ J( )] β H( β) L = = [ ] (B) L H( β) = (B) hs mples that L H( β ) stll has a zero egevalue wth egevector ad therefore that λ ( L H( β)) =. hs proves that L H(β) s postve sem-defte ad cocludes the proof of the system boudg theorem. IX. REFERECES [] A. Brooke, D. Kedrck, A. Meeraus ad R. Rama. GAMS/CLEX 7. User otes. GAMS Developmet Corporato. Washgto.. [] E. F. Hll ad W. D. Steveso. A ew Method of Determg Loss Coeffcets. IEEE ras. ower App. ad Syst. ol. AS-87, o. 7, pp , July 968. [3] S. Iwamoto ad Y. amura. A Fast Load Flow Method Retag olearty. IEEE ras. ower App. ad Syst. ol. AS-97, o. 3, pp , March 978. [4] F. D. Galaa ad M. H. Baakar. Appromato Formulae for Depedet Load Flow arables. IEEE ras. ower App. ad Syst. ol. AS-, o. 3, pp. 8-36, March 98. [5] W. Xu et al. Seres Load Flow: A ovel o-teratve Load Flow Method. IEE roc. Ge. ras. Dstrb. ol. 45, o. 3, pp. 5-56, Ja [6] A. Jag ad S. Erter. olyomal Loss Models for Ecoomc Dspatch ad Error Estmato. IEEE ras. ower Syst. ol., o. 3, pp , Aug [7] R. Baldck ad F. F. Wu. Appromato Formulas for the Dstrbuto System: the Loss Fucto ad oltage Depedece. IEEE ras. ower Delvery. ol. 6, o., pp. 5-59, Ja. 99. [8] A. L. Motto, F. D. Galaa, A. J. Coejo ad J. M. Arroyo. etwork- Costraed Mult-erod Aucto for a ool-based Electrcty Market. IEEE ras. ower Syst. ol. 7, o. 3, pp , Aug.. [9] R. A. Jabr. Modellg etwork Losses Usg Quadratc Coes. IEEE ras. ower Syst. ol., o., pp , Feb. 5. [] F. D. Galaa ad M. H. Baakar. Realzablty Iequaltes for Securty Costraed Load Flow arables. IEEE ras. Crcuts ad Syst. ol. CAS-9, o., pp , ov. 98. [] Y.. Makarov, D. J. Hll, ad I. A. Hskes. ropertes of quadratc equatos ad ther applcato to power system aalyss. Iteratoal Joural of Electrcal ower ad Eergy Systems, ol., pp ,. []. Balabaa ad. A. Bckart. Electrcal etwork heory. Joh Wley, 969. [3] J. Jarjs ad F. D. Galaa. Quattatve Aalyss of steady state stablty power etworks. IEEE ras. o ower Apparatus ad Systems, ol. AS-, o., Jauary 98l, pp X. BIOGRAHIES Sebastá de la orre (M 4) was bor Málaga, Spa. He receved the Igeero Idustral degree from the Uversdad de Málaga, Spa, 999 ad the hd degree from the Uversdad de Castlla-La Macha, Spa 3. Hs research terests clude operatos ad ecoomcs of electrcal eergy systems, restructurg of electrc systems ad developmet of ew electrcty markets.
10 Fracsco D. Galaa (F 9) receved hs B.Eg. (Hos.) from McGll Uversty, Motreal, QC, Caada followed by the S.M. ad the h.d. degrees from M.I.. He spet several years at the Brow Bover Research Ceter ad at the Uversty of Mchga. He s presetly the rofessor of Electrcal Egeerg at McGll Uversty. Hs research terest have cluded system securty, optmal power flow, epert systems, ad recetly, ope access ad competto power etworks.
Newton s Power Flow algorithm
Power Egeerg - Egll Beedt Hresso ewto s Power Flow algorthm Power Egeerg - Egll Beedt Hresso The ewto s Method of Power Flow 2 Calculatos. For the referece bus #, we set : V = p.u. ad δ = 0 For all other
More information1 Lyapunov Stability Theory
Lyapuov Stablty heory I ths secto we cosder proofs of stablty of equlbra of autoomous systems. hs s stadard theory for olear systems, ad oe of the most mportat tools the aalyss of olear systems. It may
More informationd dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin
Learzato of the Swg Equato We wll cover sectos.5.-.6 ad begg of Secto 3.3 these otes. 1. Sgle mache-fte bus case Cosder a sgle mache coected to a fte bus, as show Fg. 1 below. E y1 V=1./_ Fg. 1 The admttace
More informationThe Mathematical Appendix
The Mathematcal Appedx Defto A: If ( Λ, Ω, where ( λ λ λ whch the probablty dstrbutos,,..., Defto A. uppose that ( Λ,,..., s a expermet type, the σ-algebra o λ λ λ are defed s deoted by ( (,,...,, σ Ω.
More informationPower Flow S + Buses with either or both Generator Load S G1 S G2 S G3 S D3 S D1 S D4 S D5. S Dk. Injection S G1
ower Flow uses wth ether or both Geerator Load G G G D D 4 5 D4 D5 ecto G Net Comple ower ecto - D D ecto s egatve sg at load bus = _ G D mlarl Curret ecto = G _ D At each bus coservato of comple power
More informationSummary of the lecture in Biostatistics
Summary of the lecture Bostatstcs Probablty Desty Fucto For a cotuos radom varable, a probablty desty fucto s a fucto such that: 0 dx a b) b a dx A probablty desty fucto provdes a smple descrpto of the
More informationTESTS BASED ON MAXIMUM LIKELIHOOD
ESE 5 Toy E. Smth. The Basc Example. TESTS BASED ON MAXIMUM LIKELIHOOD To llustrate the propertes of maxmum lkelhood estmates ad tests, we cosder the smplest possble case of estmatg the mea of the ormal
More information0/1 INTEGER PROGRAMMING AND SEMIDEFINTE PROGRAMMING
CONVEX OPIMIZAION AND INERIOR POIN MEHODS FINAL PROJEC / INEGER PROGRAMMING AND SEMIDEFINE PROGRAMMING b Luca Buch ad Natala Vktorova CONENS:.Itroducto.Formulato.Applcato to Kapsack Problem 4.Cuttg Plaes
More informationECE 421/599 Electric Energy Systems 7 Optimal Dispatch of Generation. Instructor: Kai Sun Fall 2014
ECE 4/599 Electrc Eergy Systems 7 Optmal Dspatch of Geerato Istructor: Ka Su Fall 04 Backgroud I a practcal power system, the costs of geeratg ad delverg electrcty from power plats are dfferet (due to
More informationLecture 12 APPROXIMATION OF FIRST ORDER DERIVATIVES
FDM: Appromato of Frst Order Dervatves Lecture APPROXIMATION OF FIRST ORDER DERIVATIVES. INTRODUCTION Covectve term coservato equatos volve frst order dervatves. The smplest possble approach for dscretzato
More information( ) = ( ) ( ) Chapter 13 Asymptotic Theory and Stochastic Regressors. Stochastic regressors model
Chapter 3 Asmptotc Theor ad Stochastc Regressors The ature of eplaator varable s assumed to be o-stochastc or fed repeated samples a regresso aalss Such a assumpto s approprate for those epermets whch
More informationLecture 07: Poles and Zeros
Lecture 07: Poles ad Zeros Defto of poles ad zeros The trasfer fucto provdes a bass for determg mportat system respose characterstcs wthout solvg the complete dfferetal equato. As defed, the trasfer fucto
More informationFunctions of Random Variables
Fuctos of Radom Varables Chapter Fve Fuctos of Radom Varables 5. Itroducto A geeral egeerg aalyss model s show Fg. 5.. The model output (respose) cotas the performaces of a system or product, such as weght,
More informationA conic cutting surface method for linear-quadraticsemidefinite
A coc cuttg surface method for lear-quadratcsemdefte programmg Mohammad R. Osoorouch Calfora State Uversty Sa Marcos Sa Marcos, CA Jot wor wth Joh E. Mtchell RPI July 3, 2008 Outle: Secod-order coe: defto
More informationSolving Constrained Flow-Shop Scheduling. Problems with Three Machines
It J Cotemp Math Sceces, Vol 5, 2010, o 19, 921-929 Solvg Costraed Flow-Shop Schedulg Problems wth Three Maches P Pada ad P Rajedra Departmet of Mathematcs, School of Advaced Sceces, VIT Uversty, Vellore-632
More information18.413: Error Correcting Codes Lab March 2, Lecture 8
18.413: Error Correctg Codes Lab March 2, 2004 Lecturer: Dael A. Spelma Lecture 8 8.1 Vector Spaces A set C {0, 1} s a vector space f for x all C ad y C, x + y C, where we take addto to be compoet wse
More informationANALYSIS ON THE NATURE OF THE BASIC EQUATIONS IN SYNERGETIC INTER-REPRESENTATION NETWORK
Far East Joural of Appled Mathematcs Volume, Number, 2008, Pages Ths paper s avalable ole at http://www.pphm.com 2008 Pushpa Publshg House ANALYSIS ON THE NATURE OF THE ASI EQUATIONS IN SYNERGETI INTER-REPRESENTATION
More informationChapter 9 Jordan Block Matrices
Chapter 9 Jorda Block atrces I ths chapter we wll solve the followg problem. Gve a lear operator T fd a bass R of F such that the matrx R (T) s as smple as possble. f course smple s a matter of taste.
More informationSimple Linear Regression
Statstcal Methods I (EST 75) Page 139 Smple Lear Regresso Smple regresso applcatos are used to ft a model descrbg a lear relatoshp betwee two varables. The aspects of least squares regresso ad correlato
More information7.0 Equality Contraints: Lagrange Multipliers
Systes Optzato 7.0 Equalty Cotrats: Lagrage Multplers Cosder the zato of a o-lear fucto subject to equalty costrats: g f() R ( ) 0 ( ) (7.) where the g ( ) are possbly also olear fuctos, ad < otherwse
More informationG S Power Flow Solution
G S Power Flow Soluto P Q I y y * 0 1, Y y Y 0 y Y Y 1, P Q ( k) ( k) * ( k 1) 1, Y Y PQ buses * 1 P Q Y ( k1) *( k) ( k) Q Im[ Y ] 1 P buses & Slack bus ( k 1) *( k) ( k) Y 1 P Re[ ] Slack bus 17 Calculato
More informationAN UPPER BOUND FOR THE PERMANENT VERSUS DETERMINANT PROBLEM BRUNO GRENET
AN UPPER BOUND FOR THE PERMANENT VERSUS DETERMINANT PROBLEM BRUNO GRENET Abstract. The Permaet versus Determat problem s the followg: Gve a matrx X of determates over a feld of characterstc dfferet from
More informationFeature Selection: Part 2. 1 Greedy Algorithms (continued from the last lecture)
CSE 546: Mache Learg Lecture 6 Feature Selecto: Part 2 Istructor: Sham Kakade Greedy Algorthms (cotued from the last lecture) There are varety of greedy algorthms ad umerous amg covetos for these algorthms.
More informationA tighter lower bound on the circuit size of the hardest Boolean functions
Electroc Colloquum o Computatoal Complexty, Report No. 86 2011) A tghter lower boud o the crcut sze of the hardest Boolea fuctos Masak Yamamoto Abstract I [IPL2005], Fradse ad Mlterse mproved bouds o the
More informationChapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II
CEE49b Chapter - Free Vbrato of Mult-Degree-of-Freedom Systems - II We ca obta a approxmate soluto to the fudametal atural frequecy through a approxmate formula developed usg eergy prcples by Lord Raylegh
More informationEconometric Methods. Review of Estimation
Ecoometrc Methods Revew of Estmato Estmatg the populato mea Radom samplg Pot ad terval estmators Lear estmators Ubased estmators Lear Ubased Estmators (LUEs) Effcecy (mmum varace) ad Best Lear Ubased Estmators
More information{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:
Chapter 4 Exercses Samplg Theory Exercse (Smple radom samplg: Let there be two correlated radom varables X ad A sample of sze s draw from a populato by smple radom samplg wthout replacemet The observed
More informationBeam Warming Second-Order Upwind Method
Beam Warmg Secod-Order Upwd Method Petr Valeta Jauary 6, 015 Ths documet s a part of the assessmet work for the subject 1DRP Dfferetal Equatos o Computer lectured o FNSPE CTU Prague. Abstract Ths documet
More information2006 Jamie Trahan, Autar Kaw, Kevin Martin University of South Florida United States of America
SOLUTION OF SYSTEMS OF SIMULTANEOUS LINEAR EQUATIONS Gauss-Sedel Method 006 Jame Traha, Autar Kaw, Kev Mart Uversty of South Florda Uted States of Amerca kaw@eg.usf.edu Itroducto Ths worksheet demostrates
More informationLecture 3 Probability review (cont d)
STATS 00: Itroducto to Statstcal Iferece Autum 06 Lecture 3 Probablty revew (cot d) 3. Jot dstrbutos If radom varables X,..., X k are depedet, the ther dstrbuto may be specfed by specfyg the dvdual dstrbuto
More informationEstimation of Stress- Strength Reliability model using finite mixture of exponential distributions
Iteratoal Joural of Computatoal Egeerg Research Vol, 0 Issue, Estmato of Stress- Stregth Relablty model usg fte mxture of expoetal dstrbutos K.Sadhya, T.S.Umamaheswar Departmet of Mathematcs, Lal Bhadur
More informationL5 Polynomial / Spline Curves
L5 Polyomal / Sple Curves Cotets Coc sectos Polyomal Curves Hermte Curves Bezer Curves B-Sples No-Uform Ratoal B-Sples (NURBS) Mapulato ad Represetato of Curves Types of Curve Equatos Implct: Descrbe a
More informationPart 4b Asymptotic Results for MRR2 using PRESS. Recall that the PRESS statistic is a special type of cross validation procedure (see Allen (1971))
art 4b Asymptotc Results for MRR usg RESS Recall that the RESS statstc s a specal type of cross valdato procedure (see Alle (97)) partcular to the regresso problem ad volves fdg Y $,, the estmate at the
More informationCHAPTER 4 RADICAL EXPRESSIONS
6 CHAPTER RADICAL EXPRESSIONS. The th Root of a Real Number A real umber a s called the th root of a real umber b f Thus, for example: s a square root of sce. s also a square root of sce ( ). s a cube
More informationLecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model
Lecture 7. Cofdece Itervals ad Hypothess Tests the Smple CLR Model I lecture 6 we troduced the Classcal Lear Regresso (CLR) model that s the radom expermet of whch the data Y,,, K, are the outcomes. The
More informationStrong Convergence of Weighted Averaged Approximants of Asymptotically Nonexpansive Mappings in Banach Spaces without Uniform Convexity
BULLETIN of the MALAYSIAN MATHEMATICAL SCIENCES SOCIETY Bull. Malays. Math. Sc. Soc. () 7 (004), 5 35 Strog Covergece of Weghted Averaged Appromats of Asymptotcally Noepasve Mappgs Baach Spaces wthout
More informationCubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem
Joural of Amerca Scece ;6( Cubc Nopolyomal Sple Approach to the Soluto of a Secod Order Two-Pot Boudary Value Problem W.K. Zahra, F.A. Abd El-Salam, A.A. El-Sabbagh ad Z.A. ZAk * Departmet of Egeerg athematcs
More informationBounds on the expected entropy and KL-divergence of sampled multinomial distributions. Brandon C. Roy
Bouds o the expected etropy ad KL-dvergece of sampled multomal dstrbutos Brado C. Roy bcroy@meda.mt.edu Orgal: May 18, 2011 Revsed: Jue 6, 2011 Abstract Iformato theoretc quattes calculated from a sampled
More information2/20/2013. Topics. Power Flow Part 1 Text: Power Transmission. Power Transmission. Power Transmission. Power Transmission
/0/0 Topcs Power Flow Part Text: 0-0. Power Trassso Revsted Power Flow Equatos Power Flow Proble Stateet ECEGR 45 Power Systes Power Trassso Power Trassso Recall that for a short trassso le, the power
More informationLecture Notes 2. The ability to manipulate matrices is critical in economics.
Lecture Notes. Revew of Matrces he ablt to mapulate matrces s crtcal ecoomcs.. Matr a rectagular arra of umbers, parameters, or varables placed rows ad colums. Matrces are assocated wth lear equatos. lemets
More informationbest estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best
Error Aalyss Preamble Wheever a measuremet s made, the result followg from that measuremet s always subject to ucertaty The ucertaty ca be reduced by makg several measuremets of the same quatty or by mprovg
More informationPGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation
PGE 30: Formulato ad Soluto Geosystems Egeerg Dr. Balhoff Iterpolato Numercal Methods wth MATLAB, Recktewald, Chapter 0 ad Numercal Methods for Egeers, Chapra ad Caale, 5 th Ed., Part Fve, Chapter 8 ad
More information13. Parametric and Non-Parametric Uncertainties, Radial Basis Functions and Neural Network Approximations
Lecture 7 3. Parametrc ad No-Parametrc Ucertates, Radal Bass Fuctos ad Neural Network Approxmatos he parameter estmato algorthms descrbed prevous sectos were based o the assumpto that the system ucertates
More informationObjectives of Multiple Regression
Obectves of Multple Regresso Establsh the lear equato that best predcts values of a depedet varable Y usg more tha oe eplaator varable from a large set of potetal predctors {,,... k }. Fd that subset of
More informationDerivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations
Dervato of -Pot Block Method Formula for Solvg Frst Order Stff Ordary Dfferetal Equatos Kharul Hamd Kharul Auar, Kharl Iskadar Othma, Zara Bb Ibrahm Abstract Dervato of pot block method formula wth costat
More informationX X X E[ ] E X E X. is the ()m n where the ( i,)th. j element is the mean of the ( i,)th., then
Secto 5 Vectors of Radom Varables Whe workg wth several radom varables,,..., to arrage them vector form x, t s ofte coveet We ca the make use of matrx algebra to help us orgaze ad mapulate large umbers
More informationAhmed Elgamal. MDOF Systems & Modal Analysis
DOF Systems & odal Aalyss odal Aalyss (hese otes cover sectos from Ch. 0, Dyamcs of Structures, Al Chopra, Pretce Hall, 995). Refereces Dyamcs of Structures, Al K. Chopra, Pretce Hall, New Jersey, ISBN
More informationBAYESIAN INFERENCES FOR TWO PARAMETER WEIBULL DISTRIBUTION
Iteratoal Joural of Mathematcs ad Statstcs Studes Vol.4, No.3, pp.5-39, Jue 06 Publshed by Europea Cetre for Research Trag ad Developmet UK (www.eajourals.org BAYESIAN INFERENCES FOR TWO PARAMETER WEIBULL
More informationh-analogue of Fibonacci Numbers
h-aalogue of Fboacc Numbers arxv:090.0038v [math-ph 30 Sep 009 H.B. Beaoum Prce Mohammad Uversty, Al-Khobar 395, Saud Araba Abstract I ths paper, we troduce the h-aalogue of Fboacc umbers for o-commutatve
More informationNon-uniform Turán-type problems
Joural of Combatoral Theory, Seres A 111 2005 106 110 wwwelsevercomlocatecta No-uform Turá-type problems DhruvMubay 1, Y Zhao 2 Departmet of Mathematcs, Statstcs, ad Computer Scece, Uversty of Illos at
More informationAssignment 5/MATH 247/Winter Due: Friday, February 19 in class (!) (answers will be posted right after class)
Assgmet 5/MATH 7/Wter 00 Due: Frday, February 9 class (!) (aswers wll be posted rght after class) As usual, there are peces of text, before the questos [], [], themselves. Recall: For the quadratc form
More informationCOMPROMISE HYPERSPHERE FOR STOCHASTIC DOMINANCE MODEL
Sebasta Starz COMPROMISE HYPERSPHERE FOR STOCHASTIC DOMINANCE MODEL Abstract The am of the work s to preset a method of rakg a fte set of dscrete radom varables. The proposed method s based o two approaches:
More informationDescriptive Statistics
Page Techcal Math II Descrptve Statstcs Descrptve Statstcs Descrptve statstcs s the body of methods used to represet ad summarze sets of data. A descrpto of how a set of measuremets (for eample, people
More informationMATH 247/Winter Notes on the adjoint and on normal operators.
MATH 47/Wter 00 Notes o the adjot ad o ormal operators I these otes, V s a fte dmesoal er product space over, wth gve er * product uv, T, S, T, are lear operators o V U, W are subspaces of V Whe we say
More informationOrdinary Least Squares Regression. Simple Regression. Algebra and Assumptions.
Ordary Least Squares egresso. Smple egresso. Algebra ad Assumptos. I ths part of the course we are gog to study a techque for aalysg the lear relatoshp betwee two varables Y ad X. We have pars of observatos
More informationChapter 3. Linear Equations and Matrices
Vector Spaces Physcs 8/6/05 hapter Lear Equatos ad Matrces wde varety of physcal problems volve solvg systems of smultaeous lear equatos These systems of lear equatos ca be ecoomcally descrbed ad effcetly
More informationQ-analogue of a Linear Transformation Preserving Log-concavity
Iteratoal Joural of Algebra, Vol. 1, 2007, o. 2, 87-94 Q-aalogue of a Lear Trasformato Preservg Log-cocavty Daozhog Luo Departmet of Mathematcs, Huaqao Uversty Quazhou, Fua 362021, P. R. Cha ldzblue@163.com
More information1. A real number x is represented approximately by , and we are told that the relative error is 0.1 %. What is x? Note: There are two answers.
PROBLEMS A real umber s represeted appromately by 63, ad we are told that the relatve error s % What s? Note: There are two aswers Ht : Recall that % relatve error s What s the relatve error volved roudg
More information4 Inner Product Spaces
11.MH1 LINEAR ALGEBRA Summary Notes 4 Ier Product Spaces Ier product s the abstracto to geeral vector spaces of the famlar dea of the scalar product of two vectors or 3. I what follows, keep these key
More informationDKA method for single variable holomorphic functions
DKA method for sgle varable holomorphc fuctos TOSHIAKI ITOH Itegrated Arts ad Natural Sceces The Uversty of Toushma -, Mamhosama, Toushma, 770-8502 JAPAN Abstract: - Durad-Kerer-Aberth (DKA method for
More informationA COMPARATIVE STUDY OF THE METHODS OF SOLVING NON-LINEAR PROGRAMMING PROBLEM
DAODIL INTERNATIONAL UNIVERSITY JOURNAL O SCIENCE AND TECHNOLOGY, VOLUME, ISSUE, JANUARY 9 A COMPARATIVE STUDY O THE METHODS O SOLVING NON-LINEAR PROGRAMMING PROBLEM Bmal Chadra Das Departmet of Tetle
More informationTransforms that are commonly used are separable
Trasforms s Trasforms that are commoly used are separable Eamples: Two-dmesoal DFT DCT DST adamard We ca the use -D trasforms computg the D separable trasforms: Take -D trasform of the rows > rows ( )
More informationENGI 4421 Joint Probability Distributions Page Joint Probability Distributions [Navidi sections 2.5 and 2.6; Devore sections
ENGI 441 Jot Probablty Dstrbutos Page 7-01 Jot Probablty Dstrbutos [Navd sectos.5 ad.6; Devore sectos 5.1-5.] The jot probablty mass fucto of two dscrete radom quattes, s, P ad p x y x y The margal probablty
More informationThird handout: On the Gini Index
Thrd hadout: O the dex Corrado, a tala statstca, proposed (, 9, 96) to measure absolute equalt va the mea dfferece whch s defed as ( / ) where refers to the total umber of dvduals socet. Assume that. The
More informationCHAPTER VI Statistical Analysis of Experimental Data
Chapter VI Statstcal Aalyss of Expermetal Data CHAPTER VI Statstcal Aalyss of Expermetal Data Measuremets do ot lead to a uque value. Ths s a result of the multtude of errors (maly radom errors) that ca
More informationMEASURES OF DISPERSION
MEASURES OF DISPERSION Measure of Cetral Tedecy: Measures of Cetral Tedecy ad Dsperso ) Mathematcal Average: a) Arthmetc mea (A.M.) b) Geometrc mea (G.M.) c) Harmoc mea (H.M.) ) Averages of Posto: a) Meda
More informationLecture Notes Types of economic variables
Lecture Notes 3 1. Types of ecoomc varables () Cotuous varable takes o a cotuum the sample space, such as all pots o a le or all real umbers Example: GDP, Polluto cocetrato, etc. () Dscrete varables fte
More informationBERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS. Aysegul Akyuz Dascioglu and Nese Isler
Mathematcal ad Computatoal Applcatos, Vol. 8, No. 3, pp. 293-300, 203 BERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS Aysegul Ayuz Dascoglu ad Nese Isler Departmet of Mathematcs,
More informationLikewise, properties of the optimal policy for equipment replacement & maintenance problems can be used to reduce the computation.
Whe solvg a vetory repleshmet problem usg a MDP model, kowg that the optmal polcy s of the form (s,s) ca reduce the computatoal burde. That s, f t s optmal to replesh the vetory whe the vetory level s,
More informationChapter 4 (Part 1): Non-Parametric Classification (Sections ) Pattern Classification 4.3) Announcements
Aoucemets No-Parametrc Desty Estmato Techques HW assged Most of ths lecture was o the blacboard. These sldes cover the same materal as preseted DHS Bometrcs CSE 90-a Lecture 7 CSE90a Fall 06 CSE90a Fall
More informationA Study of the Reproducibility of Measurements with HUR Leg Extension/Curl Research Line
HUR Techcal Report 000--9 verso.05 / Frak Borg (borgbros@ett.f) A Study of the Reproducblty of Measuremets wth HUR Leg Eteso/Curl Research Le A mportat property of measuremets s that the results should
More informationMu Sequences/Series Solutions National Convention 2014
Mu Sequeces/Seres Solutos Natoal Coveto 04 C 6 E A 6C A 6 B B 7 A D 7 D C 7 A B 8 A B 8 A C 8 E 4 B 9 B 4 E 9 B 4 C 9 E C 0 A A 0 D B 0 C C Usg basc propertes of arthmetc sequeces, we fd a ad bm m We eed
More informationAnalysis of Variance with Weibull Data
Aalyss of Varace wth Webull Data Lahaa Watthaacheewaul Abstract I statstcal data aalyss by aalyss of varace, the usual basc assumptos are that the model s addtve ad the errors are radomly, depedetly, ad
More informationx y exp λ'. x exp λ 2. x exp 1.
egecosmcd Egevalue-egevector of the secod dervatve operator d /d hs leads to Fourer seres (se, cose, Legedre, Bessel, Chebyshev, etc hs s a eample of a systematc way of geeratg a set of mutually orthogoal
More informationn -dimensional vectors follow naturally from the one
B. Vectors ad sets B. Vectors Ecoomsts study ecoomc pheomea by buldg hghly stylzed models. Uderstadg ad makg use of almost all such models requres a hgh comfort level wth some key mathematcal sklls. I
More informationChapter 8. Inferences about More Than Two Population Central Values
Chapter 8. Ifereces about More Tha Two Populato Cetral Values Case tudy: Effect of Tmg of the Treatmet of Port-We tas wth Lasers ) To vestgate whether treatmet at a youg age would yeld better results tha
More informationChapter 5 Properties of a Random Sample
Lecture 6 o BST 63: Statstcal Theory I Ku Zhag, /0/008 Revew for the prevous lecture Cocepts: t-dstrbuto, F-dstrbuto Theorems: Dstrbutos of sample mea ad sample varace, relatoshp betwee sample mea ad sample
More informationA unified matrix representation for degree reduction of Bézier curves
Computer Aded Geometrc Desg 21 2004 151 164 wwwelsevercom/locate/cagd A ufed matrx represetato for degree reducto of Bézer curves Hask Suwoo a,,1, Namyog Lee b a Departmet of Mathematcs, Kokuk Uversty,
More informationChapter 4 Multiple Random Variables
Revew for the prevous lecture: Theorems ad Examples: How to obta the pmf (pdf) of U = g (, Y) ad V = g (, Y) Chapter 4 Multple Radom Varables Chapter 44 Herarchcal Models ad Mxture Dstrbutos Examples:
More informationLecture 5: Interpolation. Polynomial interpolation Rational approximation
Lecture 5: Iterpolato olyomal terpolato Ratoal appromato Coeffcets of the polyomal Iterpolato: Sometme we kow the values of a fucto f for a fte set of pots. Yet we wat to evaluate f for other values perhaps
More informationChapter 14 Logistic Regression Models
Chapter 4 Logstc Regresso Models I the lear regresso model X β + ε, there are two types of varables explaatory varables X, X,, X k ad study varable y These varables ca be measured o a cotuous scale as
More informationEvaluation of uncertainty in measurements
Evaluato of ucertaty measuremets Laboratory of Physcs I Faculty of Physcs Warsaw Uversty of Techology Warszawa, 05 Itroducto The am of the measuremet s to determe the measured value. Thus, the measuremet
More informationPoint Estimation: definition of estimators
Pot Estmato: defto of estmators Pot estmator: ay fucto W (X,..., X ) of a data sample. The exercse of pot estmato s to use partcular fuctos of the data order to estmate certa ukow populato parameters.
More informationLecture 2 - What are component and system reliability and how it can be improved?
Lecture 2 - What are compoet ad system relablty ad how t ca be mproved? Relablty s a measure of the qualty of the product over the log ru. The cocept of relablty s a exteded tme perod over whch the expected
More informationBlock-Based Compact Thermal Modeling of Semiconductor Integrated Circuits
Block-Based Compact hermal Modelg of Semcoductor Itegrated Crcuts Master s hess Defese Caddate: Jg Ba Commttee Members: Dr. Mg-Cheg Cheg Dr. Daqg Hou Dr. Robert Schllg July 27, 2009 Outle Itroducto Backgroud
More informationChapter 3 Sampling For Proportions and Percentages
Chapter 3 Samplg For Proportos ad Percetages I may stuatos, the characterstc uder study o whch the observatos are collected are qualtatve ature For example, the resposes of customers may marketg surveys
More informationMOLECULAR VIBRATIONS
MOLECULAR VIBRATIONS Here we wsh to vestgate molecular vbratos ad draw a smlarty betwee the theory of molecular vbratos ad Hückel theory. 1. Smple Harmoc Oscllator Recall that the eergy of a oe-dmesoal
More informationUNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS
UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS Postpoed exam: ECON430 Statstcs Date of exam: Jauary 0, 0 Tme for exam: 09:00 a.m. :00 oo The problem set covers 5 pages Resources allowed: All wrtte ad prted
More informationEVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM
EVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM Jose Javer Garca Moreta Ph. D research studet at the UPV/EHU (Uversty of Basque coutry) Departmet of Theoretcal
More informationA New Method for Solving Fuzzy Linear. Programming by Solving Linear Programming
ppled Matheatcal Sceces Vol 008 o 50 7-80 New Method for Solvg Fuzzy Lear Prograg by Solvg Lear Prograg S H Nasser a Departet of Matheatcs Faculty of Basc Sceces Mazadara Uversty Babolsar Ira b The Research
More informationMatrix Algebra Tutorial With Examples in Matlab
Matr Algebra Tutoral Wth Eamples Matlab by Klaus Moelter Departmet of Agrcultural ad Appled Ecoomcs Vrga Tech emal: moelter@vt.edu web: http://faculty.ageco.vt.edu/moelter/ Specfcally desged as a / day
More informationSimple Linear Regression
Correlato ad Smple Lear Regresso Berl Che Departmet of Computer Scece & Iformato Egeerg Natoal Tawa Normal Uversty Referece:. W. Navd. Statstcs for Egeerg ad Scetsts. Chapter 7 (7.-7.3) & Teachg Materal
More informationA New Family of Transformations for Lifetime Data
Proceedgs of the World Cogress o Egeerg 4 Vol I, WCE 4, July - 4, 4, Lodo, U.K. A New Famly of Trasformatos for Lfetme Data Lakhaa Watthaacheewakul Abstract A famly of trasformatos s the oe of several
More informationThe internal structure of natural numbers, one method for the definition of large prime numbers, and a factorization test
Fal verso The teral structure of atural umbers oe method for the defto of large prme umbers ad a factorzato test Emmaul Maousos APM Isttute for the Advacemet of Physcs ad Mathematcs 3 Poulou str. 53 Athes
More informationSolution of General Dual Fuzzy Linear Systems. Using ABS Algorithm
Appled Mathematcal Sceces, Vol 6, 0, o 4, 63-7 Soluto of Geeral Dual Fuzzy Lear Systems Usg ABS Algorthm M A Farborz Aragh * ad M M ossezadeh Departmet of Mathematcs, Islamc Azad Uversty Cetral ehra Brach,
More informationUnsupervised Learning and Other Neural Networks
CSE 53 Soft Computg NOT PART OF THE FINAL Usupervsed Learg ad Other Neural Networs Itroducto Mture Destes ad Idetfablty ML Estmates Applcato to Normal Mtures Other Neural Networs Itroducto Prevously, all
More informationQR Factorization and Singular Value Decomposition COS 323
QR Factorzato ad Sgular Value Decomposto COS 33 Why Yet Aother Method? How do we solve least-squares wthout currg codto-squarg effect of ormal equatos (A T A A T b) whe A s sgular, fat, or otherwse poorly-specfed?
More informationBayesian Inferences for Two Parameter Weibull Distribution Kipkoech W. Cheruiyot 1, Abel Ouko 2, Emily Kirimi 3
IOSR Joural of Mathematcs IOSR-JM e-issn: 78-578, p-issn: 9-765X. Volume, Issue Ver. II Ja - Feb. 05, PP 4- www.osrjourals.org Bayesa Ifereces for Two Parameter Webull Dstrbuto Kpkoech W. Cheruyot, Abel
More informationLecture Note to Rice Chapter 8
ECON 430 HG revsed Nov 06 Lecture Note to Rce Chapter 8 Radom matrces Let Y, =,,, m, =,,, be radom varables (r.v. s). The matrx Y Y Y Y Y Y Y Y Y Y = m m m s called a radom matrx ( wth a ot m-dmesoal dstrbuto,
More information1. Introduction. Keywords: Dynamic programming, Economic power dispatch, Optimization, Prohibited operating zones, Ramp-rate constraints.
A Novel TANAN s Algorthm to solve Ecoomc ower Dspatch wth Geerator Costrats ad Trasmsso Losses Subramaa R 1, Thaushkod K ad Neelakata N 3 1 Assocate rofessor/eee, Akshaya College of Egeerg ad Techology,
More information