MPC for Coupled Station Keeping, Attitude Control, and Momentum Management of Low-Thrust Geostationary Satellites

Size: px
Start display at page:

Download "MPC for Coupled Station Keeping, Attitude Control, and Momentum Management of Low-Thrust Geostationary Satellites"

Transcription

1 MITSUBISHI ELECTRIC RESEARCH LABORATORIES h:// MPC for Couled Saon Keeng, Aude Conrol, and Momenum Managemen of Low-Thrus Geosaonary Saelles Walsh, A.; D Carano, S.; Wess, A. TR26-47 July 26 Absrac Ths aer develos a model redcve conrol (MPC) olcy for smulaneous saon keeng, aude conrol, and momenum managemen of a nadr-onng geosaonary saelle equed wh hree reacon wheels and four gmbaled elecrc hrusers ha are locaed on he an-nadr face of he saelle. The MPC olcy works n combnaon wh an nnerloo SO3q-based aude conroller ha ensures he saelle manans a nadr-onng aude. The MPC olcy s able o manan he saelle s oson whn a rescrbed laude and longude wndow, whle mnmzng he v requred by he hrusers. The MPC olcy also enforces hruser onng consrans and manages he saelle s sored momenum. Wh reference o smulaon resuls, we exlan how he MPC s uned for saon keeng, he need for an nner-loo aude conroller, and how hese searae sysems work ogeher o acheve all he conroller s objecves. 26 Amercan Conrol Conference (ACC) Ths work may no be coed or reroduced n whole or n ar for any commercal urose. Permsson o coy n whole or n ar whou aymen of fee s graned for nonrof educaonal and research uroses rovded ha all such whole or aral coes nclude he followng: a noce ha such coyng s by ermsson of Msubsh Elecrc Research Laboraores, Inc.; an acknowledgmen of he auhors and ndvdual conrbuons o he work; and all alcable orons of he coyrgh noce. Coyng, reroducon, or reublshng for any oher urose shall requre a lcense wh aymen of fee o Msubsh Elecrc Research Laboraores, Inc. All rghs reserved. Coyrgh c Msubsh Elecrc Research Laboraores, Inc., 26 2 Broadway, Cambrdge, Massachuses 239

2

3 MPC for Couled Saon Keeng, Aude Conrol, and Momenum Managemen of Low-Thrus Geosaonary Saelles Alex Walsh, Sefano D Carano : and Avsha Wess ; Absrac Ths aer develos a model redcve conrol (MPC) olcy for smulaneous saon keeng, aude conrol, and momenum managemen of a nadr-onng geosaonary saelle equed wh hree reacon wheels and four gmbaled elecrc hrusers ha are locaed on he an-nadr face of he saelle. The MPC olcy works n combnaon wh an nnerloo SO3q-based aude conroller ha ensures he saelle manans a nadr-onng aude. The MPC olcy s able o manan he saelle s oson whn a rescrbed laude and longude wndow, whle mnmzng he v requred by he hrusers. The MPC olcy also enforces hruser onng consrans and manages he saelle s sored momenum. Wh reference o smulaon resuls, we exlan how he MPC s uned for saon keeng, he need for an nner-loo aude conroller, and how hese searae sysems work ogeher o acheve all he conroller s objecves. I. INTRODUCTION Recenly, here has been sgnfcan neres n he ulzaon of low-hrus, hgh-secfc-mulse elecrc engnes for saelle conrol boh for roellan-effcen nerlaneary rajecory desgn [], [2], and for roune orbal saon keeng [3] [8]. Elecrc hrusers have consderably hgher secfc mulse han ha of chemcal hrusers, allowng for a reducon n roellan mass, hereby lowerng launch coss and enablng larger ayloads [9]. Convenonal chemcal hrusers roduce hundreds of newons for orb conrol [], and n geosaonary Earh orb (GEO), may be fred oenloo once every wo weeks o comensae for saelle drf []. In conras, he maxmum hrus currenly generaed by sae-of-he-ar elecrc roulson sysems s n he hundreds of mllnewons [2], whch for GEO, s on he order of he erurbaon forces [3]. Thus, he ulzaon of elecrc hrusers oses new conrol challenges, requrng near connuous oeraon n order o rovde he necessary mared mulse, and sugges he use of closed-loo feedback conrol [3], [6], [7]. The roblem of mananng a low-hrus saelle n GEO s reaed n [6], where he auhors roose a model redcve conrol (MPC) olcy for smulaneous saon keeng and momenum managemen. The auhors develo her conrol for a saelle equed wh hree reacon wheels and sx dual-axs hrusers mouned on each face of he saelle, enablng ure hruss and orques aled n any drecon and Ph.D. Canddae, Dearmen of Aerosace Engneerng, Unversy of Mchgan, Ann Arbor, MI 489, USA. Emal: aexwalsh@umch.edu : Senor Prncal Member Research Saff, Msubsh Elecrc Research Laboraores, Cambrdge, MA 239, USA. Emal: dcarano@merl.com ; Member Research Saff, Msubsh Elecrc Research Laboraores, Cambrdge, MA 239, USA. Emal: wess@merl.com abou any axs. In hs aer, we exand on [6] o consder a realsc hruser confguraon wh gmbaled hrusers locaed on he saelle s an-nadr face [3], [4], [5]. Ths arrangemen leaves room on he saelle s oher faces for anennas and solar anels, whou he rsk of hruser lume mngemen. However, he confguraon also oses new challenges for conrol desgn, severely resrcng he conrol acon and makng more dffcul o rejec orbal erurbaons. Furhermore, o boh effcenly counerac he erurbaon forces and allow for momenum unloadng, he hrusers are ermed o gmbal. In dong so, hey may also generae undesrable orques on he saelle, sgnfcanly coulng he saelle orbal and aude dynamcs. Thus, he need o rea concurrenly boh he fas aude dynamcs and he slow orbal dynamcs resen comuaonal challenges n resource consraned hardware. Wh fas dynamcs, a small dscrezaon me-se s requred. However, a long MPC redcon horzon s necessary for saon keeng. A small dscrezaon me se accomaned by a long redcon horzon resuls n a large and comuaonally challengng omzaon roblem. The man conrbuon of hs aer s he desgn of an nner-ouer loo conrol archecure ha overcomes he aforemenoned challenges, where an ouer-loo MPC olcy oeraes n conjuncon wh an nner-loo SO3q-based aude conroller for oson, aude, and momenum conrol. The aer begns wh he descron of he saelle model n Secon II, followed by saon-keeng MPC, model lnearzaon, and consrans n Secon III. The nnerloo aude conroller s nroduced n Secon IV. Closng remarks follow he man conrbuon n Secon V. The followng noaon wll be used n all dervaons. A frame of reference F a s defned by a se of hree orhonormal dexral bass vecors, and can be wren as a T vecrx a F, where F r a a 2 a a 3 s. A hyscal vecor descrbng he oson of a on y 2 relave o on y s y wren as v 2y. A hyscal vecor can be resolved n any reference frame, such as, v a v a ` a 2 v a2 ` a 3 v a3 T rv a v a2 v a3 sf a F v a a, where he subscr a on he column marx of comonens v a ndcaes ha he hyscal T vecor v s resolved n F a. Snce v F v T a a F v b b, alyng a lef do-roduc wh b F yelds v b b F F T v a a C ba v a. The marx C ba P SO3q s a drecon cosne marx (DCM) defned by C T ba b F F, where a SO3q C P R 3ˆ3 C T C, de C u, and s he deny marx. Prncal roaons abou he a bass vecor

4 by an angle α are denoed by C ba C αq. The cross oeraor mas a hree-dmensonal vecor o a skew symmerc marx, ha s qˆ : R 3 Ñ so3q, where so3q U P R 3ˆ3 U U T u. The uncross ma does he oose where q v : so3q Ñ R 3, where vˆq v P R 3. The skew-symmerc rojecon oeraor s P a Uq 2 U UT q, for all U P R nˆn. II. SATELLITE MODEL Consder he saelle deced n Fg. ha has hree axsymmerc reacon wheels aached o a rgd bus n an orhogonal and mass balanced confguraon. Four elecrc hrusers are locaed on he an-nadr face of he saelle. The hrusers rovde hrus n a lne concden wh her oson and he cener of mass of he saelle, bu are able o gmbal away from hs nomnal drecon. The saelle s cener of mass s denoed by c and has a bus-fxed frame F. The frame F g s he Earh Cenred Ineral (ECI) frame, w s an unforced arcle, and s assumed ha w s collocaed wh he cener of he Earh. Each hruser has a frame F, locaed a on for hrusers,..., 4. w f f 2 3 F g ϕ hrus 25 P r c 2 F c m θ hrus 4 b cenrod.25 m 5 m 7.5 m 2» r cb 2.25 f 2.5 fl m 3.75 Fg. : Saelle B wh laform P showng elecrc hrusers. A force f s shown when he hrus s erfecly n lne wh he cener of mass of he saelle. F 3 The equaons of moon of he saelle are gven by J Bc 9ωg g ` µ rcw g r cw 3 a g ` f hrus m B :r cw ` ω gˆ g g, `JBc ω g ` J s 9γ ` J s η τ ` τ hrus, 9C g ω gˆ C g, :γ η, () where r cw g s he oson of he saelle and γ s a column marx conanng he angle of roaon of each reacon wheel. The vecor ω g s he angular velocy of F relave o F g resolved n F. The marx J Bc s he momen of nera of he saelle B relave o s cener of mass, resolved n F. Each hruser ales a force f a each on, and s resolved n a frame dedcaed o each hruser F. The reacon wheel array has a momen of nera J s and he wheels are conrolled wh an acceleraon η. The erm a g reresens he exernal erurbaons on he saelle due o Earh s non shercal gravaonal feld, solar and lunar gravaonal aracon, solar radaon ressure (SRP), and are defned n [6, Eq. (2)]. The erm τ reresens he SRP erurbaon orque, whch assumes oal absoron, and s gven by [6,. 229]. The hrusers are laced on he rear sde of he saelle, and each hruser has a oson r b. The nomnal hrus drecon s along he vecor r c 2. The frame F has s un vecor defned such ha 2 r c { r c. The DCM C s defned usng he geomery of he saelle. Each force f s aled a on. The oal force due o he hrusers aled o he saelle resolved n F g s gven by f hrus g f g C T g C T f. (2) The orque roduced by he elecrc hrusers relave o he saelle s cener of mass, resolved n he bus fxed frame F s gven by τ hrus r cˆ C T f. (3) The gmbal angle of he hruser s a quany ha can be evaluaed once he conroller has found he conrol nu f. The gmbal angles comrse wo angles, α and α 3, whch are deced n Fg. 2. The angles are defned such ha f can be descrbed by wo rncal roaons, yeldng f C T 3 α 3 qc T α q 2 f. (4) The negave sgn s due o he convenon ha a zero gmbal angle, he hrus s n he 2 drecon. The conroller comues he hrus comonens resolved n F, raher han a magnude and wo angles, because he laer roduces a nonlnear roblem and makes he omzaon much more dffcul. F 2 3 α3 Fg. 2: A saelle hruser wh gmbal angle and lnear onng consran. In hs fgure, he force f does no sasfy he onng consran as s ousde he blue shaded yramd. α f III. MPC FOR STATION KEEPING, MOMENTUM MANAGEMENT, AND ATTITUDE CONTROL The man conrbuon of hs aer s he desgn of an archecure ha allows for smulaneous saon keeng, aude conrol, and momenum managemen for he saelle shown n Fg.. The conroller s objecves are o

5 ) manan he saelle n he saon keeng wndow, 2) mnmze fuel consumon of he saelle, 3) ensure hrus remans whn rescrbed consrans, 4) manan a nadr-onng aude, and 5) manage he reacon wheel seeds. Followng [6], an MPC olcy s deal for hs ask as s naurally able o handle hese objecves. The MPC conroller solves a recedng-horzon fne-me omal conrol roblem ha s based on a sysem model subjec o onwse-nme sae and conrol consrans and a user-defned cos funcon. By usng lnearzed equaons of moon, lnear equaly and nequaly consrans, and quadrac coss on he saes and conrol, he MPC can be formulaed as a quadrac rogrammng (QP) roblem. A QP can be solved quckly and effcenly, whch makes sued for saelles where on-board comung ower may be lmed. A. Lnearzaon for MPC The dynamcs () are resolved n Hll s frame F h and lnearzed abou a nomnal crcular orb wh mean moon n, angular velocy ω, n a nadr-onng confguraon defned by F d wh zero reacon wheel seed. The aude error s aramerzed by a 3 2 Euler angle sequence gven by C d C g C T dg C φqc 2 θqc 3 ψq, where he angles ψ, θ and φ corresond o he yaw, ch and roll of he saelle. The Euler angles are condensed o a column marx δθ rδφ δθ δψs T. Defnng Ω dag 3n 2,, n 2 q, he lnearzed equaons of moon are gven by δ:r Ωδr 2ωˆ δ9r ` a h ` C T dh m B δθ 9 ωˆ δθ ` δω (5b) δ 9ω J Bc ωˆ JBc J Bc ω qˆqδω J Bc ωˆ J sδ 9γ J Bc J s η ` C T f (5a) r cˆ C T f. (5c) The dscree-me lnear model wh me se s x k` Ax k ` B w,d w k ` Bu k, (6) where w k models he dsurbances a h. Equaon (6) s used along a redcon horzon k,..., N. Durng he redcon horzon, we can aroxmae he dsurbances w k, k,..., N as he dsurbances assocaed wh he saelle a s desred oson, snce he saelle wll always be ke very close o he desred oson durng correc oeraon. The conrol nu u k s u k rf T... f T 4 4 η T s T. B. Consrans The lnearzed dynamcs are resolved n F h, wh δr rr r 2 r 3 s T. As such, he saon keeng consrans are gven by δr 2 ď r anλ,max q, δr 3 ď r anλ 2,max q, where λ,max and λ 2,max are he maxmum olerable longude and laude errors, resecvely, and where r s he nomnal orbal radus. The hrus onng consran s enforced by consranng he hrus o reman n he neror of four lanes, and s descrbed by he elemen-wse consrans D f ě, (7) where each of he four rows of D conans a normal vecor descrbng a lane. The four lanes n (7) are vsualzed n Fg. 2. The hrus magnude consran s a norm consran f ď f max, whch s nonlnear. Ths can be aroxmaed as lnear by mosng he consran ˇˇf ˇ ď f max, where f max r s T f max, and f max s he maxmum allowable hrus. Noe ha whle hs aroxmaon works well n our case, several oher aroxmaons of he norm are ossble. The aude error of he saelle s also subjec o he elemen-wse consran δθ ď θ max, where θ max s he maxmum allowable aude error. C. MPC for Saon Keeng Wh Aude Conrol The MPC olcy solves he fne-horzon omal conrol roblem a me subjec o N mn x T N Px ÿ N ` x T k Qx k ` u T k Ru k, (8) U k x k` Ax k ` B w,d w k ` Bu k, k,..., N, ˇ x xq, ˇδr2k ˇˇ ˇ ˇ ď r anλ,max q, ˇδθk ˇˇ ď θ max, ˇδr3k ˇˇ ď r anλ 2,max q, ˇ D f k ě, ˇf k ˇˇˇ ď fmax,,..., 4, where N s he redcon horzon, U u,..., u N u, and U u,..., u N u s he omzer of (8). The marces Q ě and R ą are sae and conrol weghs, and P ą s he ermnal cos deermned from he soluon of he Dscree Algebrac Rcca Equaon (DARE) for he nfne horzon roblem. The conrol s seleced as uq u. Snce he ermnal cos P s a soluon of he DARE, local sably of he equlbrum s guaraneed. Near equlbra, f sae and conrol consrans are nacve, he soluon of (8) s equvalen o ha of an LQR. Feasbly of (8) s guaraneed by relaxng he sae consrans o sof consrans wh he nroducon of slack varables. D. Smulaon Resuls The urose of he frs smulaon s o roduce resuls smlar o [6], demonsrang ha MPC aled o he nonlnear dynamcs () s able sasfy conrol objecves -5. The saelle s mass s m B 4 kg. The mean surface area s S 2 m 2, surface reflecance c refl.6, and he solar facng area s S facng 37.5 m 2. Each reacon wheel has a mass m Wk 2 kg, a radus r Wk.75 m, and a hegh of l Wk.2 m. For smlcy, he hree reacon wheels are assumed o concde wh he saelle s cener of mass. Earh s gravaonal arameer s µ E 3986 km 3 /s 2, he Sun s s µ S km 3 /s 2, he Moon s s µ M km 3 /s 2. Earh s J 2 arameer s J ,

6 s radus s R E 6378 km, and he solar radaon ressure consan s C sr N/m 2. The smulaon me s se o 9 mnues wh 5 mnues. The saelle mus reman n a saon keeng wndow of. n laude and longude. The maxmum hrus s f max. N, wh a onng angle α α The maxmum allowable error on he Euler angles s.2. The weghng marces for he saes n he smulaon are Q r Q 9r 3, Q θ 4, Q ω 2, and Q 9γ 5. The nal condon s se on January, 2 wh a semmajor axs of km and all oher orbal elemens zero. The aude nal condons are C g C 3 8 q, ω g ω, and 9γq r5 5 5s T rad/s, whch corresonds o C d and ω d. The resuls are shown n Fg. 3. As n [6], he MPC olcy s able o drve he reacon wheel seed o zero (Fg. 3a). However, he omal conrol u k causes he onng consran o be volaed, and hs s due o several reasons. Frs, he dynamcs of he nonlnear sysem and he lnear sysem dverge wh nonzero reacon wheel seeds. For examle, when he reacon wheel seeds ncrease from o 5 rad/s, a lghly damed mode aears a.5 rad/s, and a hgher damed mode aears a 95.3 rad/s. Snce here s a large dfference beween he lnearzed and nonlnear sysems, he MPC olcy redcs ha he aude consrans are sasfed when he lnear dynamcs are roagaed forward n me, as shown by he do-dashed lne n Fg. 3b. However, when he conrol s aled o he nonlnear sysem, he consrans are n fac volaed. Second, n [6], he hrusers are able o creae ure orques ha are decouled from body forces. Wh he hruser confguraon n hs aer, a hruser orque smulaneously creaes a ne body force. Ths effec coules he orbal and aude dynamcs, whch comlcaes conrol. Ths smulaon denfes addonal key roblems wh an MPC conrol olcy ha handles boh aude conrol and saon keeng. For aude conrol n [6], mn, and here 5 mn are used o handle he relavely fas aude dynamcs. Saon keeng dynamcs are much slower han he aude dynamcs, and hus he dscrezaon of he orbal dynamcs do no requre small me ses. To ake advanage of he saon keeng dynamcs, a horzon of a leas half an orb s needed. A small me se ( 5 mn) couled wh a large redcon horzon (5 h) dramacally ncreases he sze of he omzaon roblem (N 8). Moreover, hs omzaon would need o be solved every 5 mnues nsead of every hour, whch would requre a consderable ncreased use of onboard resources. IV. MPC WITH INNER-LOOP ATTITUDE CONTROL To overcome he roblems dscussed n he revous secon, an nner-loo aude conroller s roosed, wh conrol archecure descrbed n Fg. 4. The nner-loo conrols he aude va he reacon wheels, whle he MPC olcy conrols he closed-loo saelle/aude conroller sysem. Therefore, he lnearzaon of he oen-loo sysem (5) s no longer adequae, and he MPC requres he closed-loo saelle/aude conroller o be lnearzed. Seed (rad/s) δφ (deg) δθ (deg) δψ (deg) Tme (mn) (a) Reacon wheel seeds decreasng from 5 rad/s γ γ 2 γ 3 Aude Consrans Predced Aude Tme (mn) (b) Aude error (n sold blue), aude consrans (n dashed red), and redced aude from lnearzaon (n do-dashed yellow). Fg. 3: Smulaon for MPC conroller whou nner-loo conrol and 5 mn. A non-adave form of he conroller roosed n [7] s used for he nner-loo aude conroller as rovdes almos global asymoc sably for aude rackng and s caable of harmonc dsurbance rejecon. For an exlanaon of almos global, see [7]. Wh resec o he aude conroller, he SRP orque and he orques due o he elecrc hrusers are consdered o be aude dsurbances. I s assumed ha he orque dsurbance s he ouu of an LTI sysem A d, C d. Defnng S P a C d q v, we oban S 9 P a ω dˆ dsurbance τ s esmaed by C d q v, and ω d ω g C d ω dg d. The 9ˆd A dˆd ` B d ω d ` K Sq, (9a) ˆτ ds C dˆd, (9b) where K K T ą s a gan and where B d s desgned such ha (9) s osve real. Nex defne ν ω gˆ J Bc ω g ` J s 9γq J Bc K S 9 ` ω dˆ ω g q, ν 2 ˆτ ds, () ν 3 K v ω d ` K Sq K S, where K v K T v ą and K K T ą. The aude conrol law s gven by η J s ν ` ν 2 ` ν 3 q. ()

7 Subsung he lnearzaon of () no (5c) yelds δ 9ω τ hrus ` K ωˆ ` r K ` ωˆ JBc ωˆ ωˆ ` J Bc K v ωˆ K vk K q K v sδω J Bc C dˆd ı δθ, (2) where he SRP orque dsurbance has been gnored. The erm C dˆd s he ouu of (), and so he lnearzaon of (9) s also requred, and s gven by 9ˆd A dˆd ` B d δω ` B d K ωˆ qδθ, ˆτ ds C dˆd. (3a) (3b) The lnearzed redcon model for MPC s based on (5a), (5b), (2), (3), wh saes x rδr T δ9r T δθ T δω T δ 9γ T ˆd T s T. The change o he MPC olcy (8) s ha whle he regulaon of reacon wheel seeds 9γ s sll an objecve, her acceleraon η s no longer a conrol nu deermned as ar of he omzaon, bu s now rescrbed by he nner-loo conrol (). Ths ensures ha he MPC oeraes by conrollng he hrusers n a way ha causes he aude conroller () o reduce he reacon wheel seeds, and s ossble because he MPC olcy uses a closed-loo aude model ha caures he acon of he aude conroller. MPC Lnearzaon for MPC abou nomnal saes f Nonlnear Dynamcs Dsurbances Orbal Dynamcs Aude Dynamcs Aude Conrol Fg. 4: MPC wh nner-loo conroller. A. Smulaon Resuls The smulaon wh he nner-loo conroller has almos all he same arameers as n Secon III-D. The weghng marces of hese smulaons are Q r dag,, q, Q 9r and R 4 o ensure ha conrol s heavly enalzed. The sae wegh s laced on longude and laude devaons. I s ossble o use mnmal weghng on he saes because he saon keeng consrans are bul no he MPC olcy. The weghng for he aude relaed saes are Q θ Q ω 3, Q 9γ 2, and Qˆd. The marx B d s deermned by choosng Q d.5, solvng A T d P d `P d A d Q d for P d, and hen seng B d C TP d d. The horzon for MPC s chosen o be 5 hours o exlo he erodcy of he orbal dsurbances. In addon, he dynamcs are dscrezed usng h, as n [6], whch s now ossble because he fas aude dynamcs are managed by he nner-loo conroller, and hus N 5. The nal condons are he same as Secon III-D, exce wh 9γ. Cd! d 9 r 9r! 9 ˆd The smulaon s run for 425 orbs, bu only he las year (365 orbs) s used for analyss o solae he seady sae of he sysem (see Fg. 5). Fg. 5a shows ha he saelle remans whn he saon keeng wndow. Fg. 5c shows he aude error, whch remans whn.2. Fg. 5d shows ha he reacon wheel seed s conrolled and does no ncrease whou bound, as s he case n smulaons where momenum managemen s no erformed. Fg. 5b shows he cumulave v roduced by each hruser, and he v due o he ne force of he hrusers, resolved n F h. The v 3 s slghly less han 6 m/s, whch s execed due o he ou-of-lane erurbaons. The oal v.97 m/s s 22.3 m/s more han requred for saon keeng only. The exra hrus s needed o counerac he solar radaon ressure orque and o manage he momenum of he reacon wheels. Fg. 6 shows he las fve orbs are loed for enhanced deal. The reacon wheel seeds are cyclc over an orb as execed. The hrusers are held consan over a one hour erod, and he MPC soluon resuls n only wo hrusers beng acve a a me, as shown n Fg. 6d. Eher he wo Norh hrusers (forces a and 2 ), or he wo Souh hrusers (forces a 3 and 4 ) are acve smulaneously. Noe ha due o he naure of he MPC olcy, he hrus s held consan for one hour, and hus he gmbal angles are also consan. Gmbal angles do no change from `36 o 36 from one hour o he nex, bu do so n ncremens. The gmbal angles are ofen a her lms, whch ndcaes ha he hrusers are almos always creang some sor of orque, and mlyng ha he reacon wheels are always n use. V. CLOSING REMARKS Auonomous closed-loo feedback conrol for saon keeng oeraons wll ncrease he safey, robusness, and relably of saelles. In hs aer, we resened an MPC olcy wh an nner-loo SO3q-based aude conroller ha s able o sasfy saon keeng consrans, whle also rovdng momenum managemen. The SO3q-based aude conroller s also able o kee he saelle s aude error whn he onng consrans. All hs can be acheved for a saelle confguraon wh he hrusers locaed on he an-nadr face of he saelle, whch s convenen for saelle desgn bu challengng for conrol. REFERENCES [] J. A. Sarek and I. V. Kolmanovsky, Nonlnear model redcve conrol sraegy for low hrus sacecraf mssons, Omal Conrol Alcaons and Mehods, vol. 35, no.,. 2, 24. [2] O. Abdelkhalk and E. Taher, Aroxmae on-off low-hrus sace rajecores usng fourer seres, Journal of Sacecraf and Rockes, vol. 49, no. 5, , 22. [3] D. Losa, M. Lovera, J.-P. Marmora, T. Dargen, and J. Amalrc, Saon keeng of geosaonary saelles wh on-off elecrc hrusers, n Proc. 26 IEEE Inernaonal Conference on Conrol Alcaons, , 26. [4] A. Sukhanov and A. Prado, On one aroach o he omzaon of low-hrus saon keeng manoeuvres, Advances n Sace Research, vol. 5, no., , 22. [5] Y. Ulybyshev, Long-erm saon keeng of sace saon n lunar halo orbs, Journal of Gudance, Conrol, and Dynamcs, vol. 38, no. 6,. 63 7, 24.

8 . 3 v n frames F, oal v =.973m/s.5 v (m/s) τ w (N m) Seed (rad/s) Aude error (deg) Aude error (deg) - -2 Laude (deg) Longude (deg) (a) Laude and longude error. δφ δθ δψ Tme (orb) (c) Euler angle aude error. Seed (rad/s) v (m/s) v v2 v3 v n Hll s frame Fh Tme (orb) (b) v from each hruser and ne v resolved n F h. 4 2 γ γ 2 γ Tme (orb) (d) Reacon wheel seeds. Fg. 5: One year smulaon for MPC conroller wh nner-loo conrol Tme (orb) (a) Reacon wheel seeds and nu. γ γ 2 γ 3 J sη J sη 2 J sη 3-3 δφ δθ δψ Tme (orb) (c) Aude error. Maxmum error s.2. α (deg) α3 (deg) f (mn) 5 f f 2 f 3 f Tme (orb) (b) Gmbal angles for each hruser wh 36 consran Tme (orb) (d) Thrus magnude. Maxmum hrus s mn. Fg. 6: Close-u of 5 orbs dealng reacon wheels and hrusers [6] A. Wess, U. Kalabc, and S. D Carano, Model redcve conrol for smulaneous saon keeng and momenum managemen of low-hrus saelles, n Proc. 25 Amercan Conrol Conference, , 25. [7] U. Kalabc, A. Wess, S. D Carano, and I. Kolmanovsky, Saonkeeng and momenum-managemen on halo orbs around L2: Lnear-quadrac feedback and model redcve conrol aroaches, n 25h AAS/AIAA Sace Flgh Mechancs Meeng, AAS 5-37, 25. [8] A. Garull, A. Gannraan, M. Leomann, and F. Scorecc, Auonomous low-earh-orb saon-keeng wh elecrc roulson, Journal of Gudance, Conrol, and Dynamcs, vol. 34, no. 6, , 2. [9] M. Marnez-Sanchez and J. E. Pollard, Sacecraf elecrc roulson an overvew, Journal of Proulson and Power, vol. 4, no. 5, , 998. [] J. R. Werz and W. J. Larson, eds., Sacecraf Msson Analyss and Desgn. Mcrocosm Press, 3rd ed., 999. [] E. M. Soo, Handbook of Geosaonary Orbs. Srnger, 994. [2] D. Herman, G. C. Soulas, J. L. Van Noord, and M. J. Paerson, NASA s evoluonary xenon hruser (NEXT) long-duraon es resuls, Journal of Proulson and Power, vol. 28, , 25/9/ [3] R. R. Bae, D. D. Mueller, and J. E. Whe, Fundamenals of Asrodynamcs. Dover Publcaons Inc., 97. [4] J. Clmer, Boeng: World s frs all-elecrc roulson saelle begns oeraons, Inerne: h://boeng.medaroom.com/ndex.h?s=2295&em=2956, Se. 25, [22 Se. 25]. [5] B. M. Anzel, Mehod and aaraus for a saelle saon keeng. U.S. Paen , Augus 22, 995. [6] A. H. J. de Ruer, C. J. Damaren, and J. R. Forbes, Sacecraf Dynamcs and Conrol: An Inroducon. Chcheser, Wes Sussex, U.K.: John Wley & Sons, Ld, 23. [7] A. Wess, I. Kolmanovsky, D. S. Bernsen, and A. Sanyal, Inerafree sacecraf aude conrol usng reacon wheels, Journal of Gudance, Conrol, and Dynamcs, vol. 36, no. 5, , 23.

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

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

More information

Motion in Two Dimensions

Motion in Two Dimensions Phys 1 Chaper 4 Moon n Two Dmensons adzyubenko@csub.edu hp://www.csub.edu/~adzyubenko 005, 014 A. Dzyubenko 004 Brooks/Cole 1 Dsplacemen as a Vecor The poson of an objec s descrbed by s poson ecor, r The

More information

Graduate Macroeconomics 2 Problem set 5. - Solutions

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

More information

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

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

More information

A Cell Decomposition Approach to Online Evasive Path Planning and the Video Game Ms. Pac-Man

A Cell Decomposition Approach to Online Evasive Path Planning and the Video Game Ms. Pac-Man Cell Decomoson roach o Onlne Evasve Pah Plannng and he Vdeo ame Ms. Pac-Man reg Foderaro Vram Raju Slva Ferrar Laboraory for Inellgen Sysems and Conrols LISC Dearmen of Mechancal Engneerng and Maerals

More information

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

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

More information

FI 3103 Quantum Physics

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

More information

Solution in semi infinite diffusion couples (error function analysis)

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

More information

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

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

More information

Variants of Pegasos. December 11, 2009

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

More information

Mechanics Physics 151

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

More information

On One Analytic Method of. Constructing Program Controls

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

More information

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

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

More information

Laser Interferometer Space Antenna (LISA)

Laser Interferometer Space Antenna (LISA) aser nerferomeer Sace Anenna SA Tme-elay nerferomery wh Movng Sacecraf Arrays Massmo Tno Je Proulson aboraory, Calforna nsue of Technology GSFC JP 8 h GWAW, ec 7-0, 00, Mlwaukee, Wsconsn WM Folkner e al,

More information

Mechanics Physics 151

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

More information

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

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

More information

Advanced time-series analysis (University of Lund, Economic History Department)

Advanced time-series analysis (University of Lund, Economic History Department) Advanced me-seres analss (Unvers of Lund, Economc Hsor Dearmen) 3 Jan-3 Februar and 6-3 March Lecure 4 Economerc echnues for saonar seres : Unvarae sochasc models wh Box- Jenns mehodolog, smle forecasng

More information

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation

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

More information

Linear Response Theory: The connection between QFT and experiments

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

More information

PHYS 705: Classical Mechanics. Canonical Transformation

PHYS 705: Classical Mechanics. Canonical Transformation PHYS 705: Classcal Mechancs Canoncal Transformaon Canoncal Varables and Hamlonan Formalsm As we have seen, n he Hamlonan Formulaon of Mechancs,, are ndeenden varables n hase sace on eual foong The Hamlon

More information

Density Matrix Description of NMR BCMB/CHEM 8190

Density Matrix Description of NMR BCMB/CHEM 8190 Densy Marx Descrpon of NMR BCMBCHEM 89 Operaors n Marx Noaon Alernae approach o second order specra: ask abou x magnezaon nsead of energes and ranson probables. If we say wh one bass se, properes vary

More information

Density Matrix Description of NMR BCMB/CHEM 8190

Density Matrix Description of NMR BCMB/CHEM 8190 Densy Marx Descrpon of NMR BCMBCHEM 89 Operaors n Marx Noaon If we say wh one bass se, properes vary only because of changes n he coeffcens weghng each bass se funcon x = h< Ix > - hs s how we calculae

More information

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

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

More information

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

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

More information

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

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

More information

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

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

More information

A New Generalized Gronwall-Bellman Type Inequality

A New Generalized Gronwall-Bellman Type Inequality 22 Inernaonal Conference on Image, Vson and Comung (ICIVC 22) IPCSIT vol. 5 (22) (22) IACSIT Press, Sngaore DOI:.7763/IPCSIT.22.V5.46 A New Generalzed Gronwall-Bellman Tye Ineualy Qnghua Feng School of

More information

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

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

More information

Response of MDOF systems

Response of MDOF systems Response of MDOF syses Degree of freedo DOF: he nu nuber of ndependen coordnaes requred o deerne copleely he posons of all pars of a syse a any nsan of e. wo DOF syses hree DOF syses he noral ode analyss

More information

Imperfect Information

Imperfect Information Imerfec Informaon Comlee Informaon - all layers know: Se of layers Se of sraeges for each layer Oucomes as a funcon of he sraeges Payoffs for each oucome (.e. uly funcon for each layer Incomlee Informaon

More information

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

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

More information

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

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

More information

PHYS 1443 Section 001 Lecture #4

PHYS 1443 Section 001 Lecture #4 PHYS 1443 Secon 001 Lecure #4 Monda, June 5, 006 Moon n Two Dmensons Moon under consan acceleraon Projecle Moon Mamum ranges and heghs Reerence Frames and relae moon Newon s Laws o Moon Force Newon s Law

More information

FTCS Solution to the Heat Equation

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

More information

Chapter 6: AC Circuits

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

More information

Mechanics Physics 151

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

More information

Chapter Lagrangian Interpolation

Chapter Lagrangian Interpolation Chaper 5.4 agrangan Inerpolaon Afer readng hs chaper you should be able o:. dere agrangan mehod of nerpolaon. sole problems usng agrangan mehod of nerpolaon and. use agrangan nerpolans o fnd deraes and

More information

Robustness Experiments with Two Variance Components

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

More information

CS286.2 Lecture 14: Quantum de Finetti Theorems II

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

More information

NPTEL Project. Econometric Modelling. Module23: Granger Causality Test. Lecture35: Granger Causality Test. Vinod Gupta School of Management

NPTEL Project. Econometric Modelling. Module23: Granger Causality Test. Lecture35: Granger Causality Test. Vinod Gupta School of Management P age NPTEL Proec Economerc Modellng Vnod Gua School of Managemen Module23: Granger Causaly Tes Lecure35: Granger Causaly Tes Rudra P. Pradhan Vnod Gua School of Managemen Indan Insue of Technology Kharagur,

More information

Foundations of State Estimation Part II

Foundations of State Estimation Part II Foundaons of Sae Esmaon Par II Tocs: Hdden Markov Models Parcle Flers Addonal readng: L.R. Rabner, A uoral on hdden Markov models," Proceedngs of he IEEE, vol. 77,. 57-86, 989. Sequenal Mone Carlo Mehods

More information

Chapters 2 Kinematics. Position, Distance, Displacement

Chapters 2 Kinematics. Position, Distance, Displacement Chapers Knemacs Poson, Dsance, Dsplacemen Mechancs: Knemacs and Dynamcs. Knemacs deals wh moon, bu s no concerned wh he cause o moon. Dynamcs deals wh he relaonshp beween orce and moon. The word dsplacemen

More information

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

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

More information

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

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

More information

An ant colony optimization solution to the integrated generation and transmission maintenance scheduling problem

An ant colony optimization solution to the integrated generation and transmission maintenance scheduling problem JOURNAL OF OTOELECTRONICS AND ADVANCED MATERIALS Vol. 0, No. 5, May 008,. 46-50 An an colony omzaon soluon o he negraed generaon and ransmsson manenance schedulng roblem. S. GEORGILAKIS *,. G. VERNADOS

More information

Decentralised Sliding Mode Load Frequency Control for an Interconnected Power System with Uncertainties and Nonlinearities

Decentralised Sliding Mode Load Frequency Control for an Interconnected Power System with Uncertainties and Nonlinearities Inernaonal Research Journal of Engneerng and echnology IRJE e-iss: 2395-0056 Volume: 03 Issue: 12 Dec -2016 www.re.ne p-iss: 2395-0072 Decenralsed Sldng Mode Load Frequency Conrol for an Inerconneced Power

More information

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study)

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study) Inernaonal Mahemacal Forum, Vol. 8, 3, no., 7 - HIKARI Ld, www.m-hkar.com hp://dx.do.org/.988/mf.3.3488 New M-Esmaor Objecve Funcon n Smulaneous Equaons Model (A Comparave Sudy) Ahmed H. Youssef Professor

More information

2/20/2013. EE 101 Midterm 2 Review

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

More information

THEORETICAL AUTOCORRELATIONS. ) if often denoted by γ. Note that

THEORETICAL AUTOCORRELATIONS. ) if often denoted by γ. Note that THEORETICAL AUTOCORRELATIONS Cov( y, y ) E( y E( y))( y E( y)) ρ = = Var( y) E( y E( y)) =,, L ρ = and Cov( y, y ) s ofen denoed by whle Var( y ) f ofen denoed by γ. Noe ha γ = γ and ρ = ρ and because

More information

Set point control in the state space setting

Set point control in the state space setting Downloaded from orb.du.dk on: Apr 28, 28 Se pon conrol n he sae space seng Poulsen, Nels Kjølsad Publcaon dae: 29 Documen Verson Publsher's PDF, also known as Verson of record Lnk back o DU Orb Caon (APA):

More information

CONSISTENT ESTIMATION OF THE NUMBER OF DYNAMIC FACTORS IN A LARGE N AND T PANEL. Detailed Appendix

CONSISTENT ESTIMATION OF THE NUMBER OF DYNAMIC FACTORS IN A LARGE N AND T PANEL. Detailed Appendix COSISE ESIMAIO OF HE UMBER OF DYAMIC FACORS I A LARGE AD PAEL Dealed Aendx July 005 hs verson: May 9, 006 Dane Amengual Dearmen of Economcs, Prnceon Unversy and Mar W Wason* Woodrow Wlson School and Dearmen

More information

P R = P 0. The system is shown on the next figure:

P R = P 0. The system is shown on the next figure: TPG460 Reservor Smulaon 08 page of INTRODUCTION TO RESERVOIR SIMULATION Analycal and numercal soluons of smple one-dmensonal, one-phase flow equaons As an nroducon o reservor smulaon, we wll revew he smples

More information

Thruster Modulation for Unsymmetric Flexible Spacecraft with Consideration of Torque Arm Perturbation

Thruster Modulation for Unsymmetric Flexible Spacecraft with Consideration of Torque Arm Perturbation hruer Modulaon for Unymmerc Flexble Sacecraf wh onderaon of orue rm Perurbaon a Shgemune anwak Shnchro chkawa a Yohak hkam b a Naonal Sace evelomen gency of Jaan 2-- Sengen ukuba-h barak b eo Unvery 3--

More information

THERMODYNAMICS 1. The First Law and Other Basic Concepts (part 2)

THERMODYNAMICS 1. The First Law and Other Basic Concepts (part 2) Company LOGO THERMODYNAMICS The Frs Law and Oher Basc Conceps (par ) Deparmen of Chemcal Engneerng, Semarang Sae Unversy Dhon Harano S.T., M.T., M.Sc. Have you ever cooked? Equlbrum Equlbrum (con.) Equlbrum

More information

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

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

More information

AQM Algorithm Based on Kelly s Scheme Using Sliding Mode Control

AQM Algorithm Based on Kelly s Scheme Using Sliding Mode Control 009 Amercan Conrol Conference Hya Regency Rverfron, S. Lous, MO, USA June 0-, 009 WeC06.6 AQM Algorhm Based on Kelly s Scheme Usng Sldng Mode Conrol Nannan Zhang, Georg M. Dmrovsk, Yuanwe Jng, and Syng

More information

Lecture 11 SVM cont

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

More information

Scattering at an Interface: Oblique Incidence

Scattering at an Interface: Oblique Incidence Course Insrucor Dr. Raymond C. Rumpf Offce: A 337 Phone: (915) 747 6958 E Mal: rcrumpf@uep.edu EE 4347 Appled Elecromagnecs Topc 3g Scaerng a an Inerface: Oblque Incdence Scaerng These Oblque noes may

More information

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

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

More information

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

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

More information

Comparison of Differences between Power Means 1

Comparison of Differences between Power Means 1 In. Journal of Mah. Analyss, Vol. 7, 203, no., 5-55 Comparson of Dfferences beween Power Means Chang-An Tan, Guanghua Sh and Fe Zuo College of Mahemacs and Informaon Scence Henan Normal Unversy, 453007,

More information

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

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

More information

DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL

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

More information

. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue.

. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue. Mah E-b Lecure #0 Noes We connue wh he dscusson of egenvalues, egenvecors, and dagonalzably of marces We wan o know, n parcular wha condons wll assure ha a marx can be dagonalzed and wha he obsrucons are

More information

MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES. Institute for Mathematical Research, Universiti Putra Malaysia, UPM Serdang, Selangor, Malaysia

MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES. Institute for Mathematical Research, Universiti Putra Malaysia, UPM Serdang, Selangor, Malaysia Malaysan Journal of Mahemacal Scences 9(2): 277-300 (2015) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homeage: h://ensemumedumy/journal A Mehod for Deermnng -Adc Orders of Facorals 1* Rafka Zulkal,

More information

Department of Economics University of Toronto

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

More information

Comb Filters. Comb Filters

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

More information

Tight results for Next Fit and Worst Fit with resource augmentation

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

More information

Volatility Interpolation

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

More information

Beyond Balanced Growth : Some Further Results

Beyond Balanced Growth : Some Further Results eyond alanced Growh : Some Furher Resuls by Dens Sec and Helmu Wagner Dscusson Paer o. 49 ay 27 Dskussonsberäge der Fakulä für Wrschafswssenschaf der FernUnversä n Hagen Herausgegeben vom Dekan der Fakulä

More information

Inverse Joint Moments of Multivariate. Random Variables

Inverse Joint Moments of Multivariate. Random Variables In J Conem Mah Scences Vol 7 0 no 46 45-5 Inverse Jon Momens of Mulvarae Rom Varables M A Hussan Dearmen of Mahemacal Sascs Insue of Sascal Sudes Research ISSR Caro Unversy Egy Curren address: Kng Saud

More information

First-order piecewise-linear dynamic circuits

First-order piecewise-linear dynamic circuits Frs-order pecewse-lnear dynamc crcus. Fndng he soluon We wll sudy rs-order dynamc crcus composed o a nonlnear resse one-por, ermnaed eher by a lnear capacor or a lnear nducor (see Fg.. Nonlnear resse one-por

More information

ISSN MIT Publications

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

More information

Should Exact Index Numbers have Standard Errors? Theory and Application to Asian Growth

Should Exact Index Numbers have Standard Errors? Theory and Application to Asian Growth Should Exac Index umbers have Sandard Errors? Theory and Applcaon o Asan Growh Rober C. Feensra Marshall B. Rensdorf ovember 003 Proof of Proposon APPEDIX () Frs, we wll derve he convenonal Sao-Vara prce

More information

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

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

More information

How about the more general "linear" scalar functions of scalars (i.e., a 1st degree polynomial of the following form with a constant term )?

How about the more general linear scalar functions of scalars (i.e., a 1st degree polynomial of the following form with a constant term )? lmcd Lnear ransformaon of a vecor he deas presened here are que general hey go beyond he radonal mar-vecor ype seen n lnear algebra Furhermore, hey do no deal wh bass and are equally vald for any se of

More information

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

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

More information

Motion of Wavepackets in Non-Hermitian. Quantum Mechanics

Motion of Wavepackets in Non-Hermitian. Quantum Mechanics Moon of Wavepaces n Non-Herman Quanum Mechancs Nmrod Moseyev Deparmen of Chemsry and Mnerva Cener for Non-lnear Physcs of Complex Sysems, Technon-Israel Insue of Technology www.echnon echnon.ac..ac.l\~nmrod

More information

CHAPTER 10: LINEAR DISCRIMINATION

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

More information

Energy Storage Devices

Energy Storage Devices Energy Sorage Deces Objece of Lecure Descrbe he consrucon of a capacor and how charge s sored. Inroduce seeral ypes of capacors Dscuss he elecrcal properes of a capacor The relaonshp beween charge, olage,

More information

Cubic Bezier Homotopy Function for Solving Exponential Equations

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

More information

An introduction to Support Vector Machine

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

More information

Dynamic Poverty Measures

Dynamic Poverty Measures heorecal Economcs Leers 63-69 do:436/el34 Publshed Onlne November (h://wwwscrporg/journal/el) Dynamc Povery Measures Absrac Eugene Kouass Perre Mendy Dara Seck Kern O Kymn 3 Resource Economcs Wes Vrgna

More information

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION

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

More information

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy Arcle Inernaonal Journal of Modern Mahemacal Scences, 4, (): - Inernaonal Journal of Modern Mahemacal Scences Journal homepage: www.modernscenfcpress.com/journals/jmms.aspx ISSN: 66-86X Florda, USA Approxmae

More information

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

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

More information

Lecture 2 M/G/1 queues. M/G/1-queue

Lecture 2 M/G/1 queues. M/G/1-queue Lecure M/G/ queues M/G/-queue Posson arrval process Arbrary servce me dsrbuon Sngle server To deermne he sae of he sysem a me, we mus now The number of cusomers n he sysems N() Tme ha he cusomer currenly

More information

Genetic Algorithm in Parameter Estimation of Nonlinear Dynamic Systems

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

More information

Outline. Energy-Efficient Target Coverage in Wireless Sensor Networks. Sensor Node. Introduction. Characteristics of WSN

Outline. Energy-Efficient Target Coverage in Wireless Sensor Networks. Sensor Node. Introduction. Characteristics of WSN Ener-Effcen Tare Coverae n Wreless Sensor Newors Presened b M Trà Tá -4-4 Inroducon Bacround Relaed Wor Our Proosal Oulne Maxmum Se Covers (MSC) Problem MSC Problem s NP-Comlee MSC Heursc Concluson Sensor

More information

Normal Random Variable and its discriminant functions

Normal Random Variable and its discriminant functions Noral Rando Varable and s dscrnan funcons Oulne Noral Rando Varable Properes Dscrnan funcons Why Noral Rando Varables? Analycally racable Works well when observaon coes for a corruped snle prooype 3 The

More information

E c. i f. (c) The built-in potential between the collector and emitter is. 18 ae bicb

E c. i f. (c) The built-in potential between the collector and emitter is. 18 ae bicb haer 8 nergy and Dagram of T 8 a & b or he gven dong concenraons, one comues f - = -05 ev, 049 ev, and 099 ev n he emer, base and collecor, resecvely Also wh a >> d, he - deleon wdh wll le almos eclusvely

More information

Homework 8: Rigid Body Dynamics Due Friday April 21, 2017

Homework 8: Rigid Body Dynamics Due Friday April 21, 2017 EN40: Dynacs and Vbraons Hoework 8: gd Body Dynacs Due Frday Aprl 1, 017 School of Engneerng Brown Unversy 1. The earh s roaon rae has been esaed o decrease so as o ncrease he lengh of a day a a rae of

More information

On computing differential transform of nonlinear non-autonomous functions and its applications

On computing differential transform of nonlinear non-autonomous functions and its applications On compung dfferenal ransform of nonlnear non-auonomous funcons and s applcaons Essam. R. El-Zahar, and Abdelhalm Ebad Deparmen of Mahemacs, Faculy of Scences and Humanes, Prnce Saam Bn Abdulazz Unversy,

More information

OP = OO' + Ut + Vn + Wb. Material We Will Cover Today. Computer Vision Lecture 3. Multi-view Geometry I. Amnon Shashua

OP = OO' + Ut + Vn + Wb. Material We Will Cover Today. Computer Vision Lecture 3. Multi-view Geometry I. Amnon Shashua Comuer Vson 27 Lecure 3 Mul-vew Geomer I Amnon Shashua Maeral We Wll Cover oa he srucure of 3D->2D rojecon mar omograh Marces A rmer on rojecve geomer of he lane Eolar Geomer an Funamenal Mar ebrew Unvers

More information

Endogeneity. Is the term given to the situation when one or more of the regressors in the model are correlated with the error term such that

Endogeneity. Is the term given to the situation when one or more of the regressors in the model are correlated with the error term such that s row Endogeney Is he erm gven o he suaon when one or more of he regressors n he model are correlaed wh he error erm such ha E( u 0 The 3 man causes of endogeney are: Measuremen error n he rgh hand sde

More information

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

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

More information

Relative controllability of nonlinear systems with delays in control

Relative controllability of nonlinear systems with delays in control Relave conrollably o nonlnear sysems wh delays n conrol Jerzy Klamka Insue o Conrol Engneerng, Slesan Techncal Unversy, 44- Glwce, Poland. phone/ax : 48 32 37227, {jklamka}@a.polsl.glwce.pl Keywor: Conrollably.

More information

Pavel Azizurovich Rahman Ufa State Petroleum Technological University, Kosmonavtov St., 1, Ufa, Russian Federation

Pavel Azizurovich Rahman Ufa State Petroleum Technological University, Kosmonavtov St., 1, Ufa, Russian Federation VOL., NO. 5, MARCH 8 ISSN 89-668 ARN Journal of Engneerng and Aled Scences 6-8 Asan Research ublshng Nework ARN. All rghs reserved. www.arnjournals.com A CALCULATION METHOD FOR ESTIMATION OF THE MEAN TIME

More information

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

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

More information

Keywords: Hedonic regressions; hedonic indexes; consumer price indexes; superlative indexes.

Keywords: Hedonic regressions; hedonic indexes; consumer price indexes; superlative indexes. Hedonc Imuaon versus Tme Dummy Hedonc Indexes Erwn Dewer, Saeed Herav and Mck Slver December 5, 27 (wh a commenary by Jan de Haan) Dscusson Paer 7-7, Dearmen of Economcs, Unversy of Brsh Columba, 997-873

More information