Reconfigurable UKF for In-Flight Magnetometer Calibration and Attitude Parameter Estimation

Size: px
Start display at page:

Download "Reconfigurable UKF for In-Flight Magnetometer Calibration and Attitude Parameter Estimation"

Transcription

1 Prprnts of th 8th IFAC World Congrss Mlano (Italy) August 8 - Sptmbr, Rconfgurabl UKF for In-Flght Magntomtr Calbraton and Atttud Paramtr Estmaton Hall E. Son*. Chngz Hajyv.** * h Graduat Unvrsty for Advancd Studs (SOKENDAI), Dpartmnt of Spac and Astronautcal Scnc 3--, Yoshnoda, Sagamhara, Kanagawa, JAPAN (l: ; -mal:rsn_son@ac.jaa.jp). ** Istanbul chncal Unvrsty Aronautcs and Astronautcs Faculty, Masla, Istanbul, URKEY (-mal: cngz@tu.du.tr). Abstract: In ths study a rconfgurabl unscntd Kalman fltr (UKF) basd algorthm for th stmaton of magntomtr bass and scal factors s proposd as a part of th atttud stmaton schm of a pco satllt. Algorthm s composd of two stags; n frst stag UKF stmats magntomtr bass and scal factors as wll as s atttud paramtrs of th satllt. Dffrntly from th stng algorthms, scal factors ar not tratd togthr wth th othr paramtrs as a part of th stat vctor; thr scal factors ar stmatd va a nwly proposd tnson for th UKF. Aftr a convrgnc rul for th bass scond stag starts and th UKF rconfgurs tslf for th stmaton of only atttud paramtrs. At ths stag fltr rgards th bass and scal factors stmatd at th ntal stag. Proposd algorthm s smulatd for atttud stmaton of a pco satllt whch has thr magntomtrs and thr rat gyros as masurmnt snsors. Kywords: Atttud algorthms, Calbraton, Kalman fltrs, Estmaton paramtrs, Satllt applcatons. INRODUCION Bcaus of advantags such as provdng contnuously avalabl two-as atttud masurmnts; rlatv low cost and almost nsgnfcant powr dmand, opratng magntomtrs as prmary snsor n small satllt mssons s a common mthod for achvng atttud nformaton. Howvr, ths snsors ar not rror fr bcaus of th bass, scalng rrors and msalgnmnts (nonorthagonalty). hs trms nhbt th fltr ffcncy and so atttud data accuracy and vn thy may brng about th fltr dvrgnc n long trms. h atttud accuracy rqurmnts dmand compnsaton for th magntomtr rrors such as msalgnmnts and bass (Wrtz, 988). Estmatng magntomtr bass and scal factors as wll as th atttud of th satllt s a proposd tchnqu to solv such problms and ncras on-board accuracy. In ltratur thr ar svral mthods for stmatng th magntomtr bas whn th atttud nowldg s not avalabl (.g. at th njcton phas whr th spaccraft s spnnng rapdly). As a rlatvly nw mthod, twostp algorthm, whch was proposd n ordr to ovrrun quartc cost functon problm mt durng atttud fr bas stmaton procss, may b gvn as an ampl (Alonso and Shustr, a). In (Alonso and Shustr, 3) t s also shown that twostp algorthm can b also mplmntd n cas of ncomplt obsrvablty of th magntomtr bas vctor. Morovr (Alonso and Shustr, b) uss algorthm to stmat scal factors and nonorthogonalty corrctons as wll as magntomtr bass. On th othr hand, Km and Bang ntgrats twostp algorthm wth th gntc algorthm, from whch th ntal stmats of th magntomtr bas stmaton ar provdd (Km and Bang, 7). wostp algorthm l mthods can b usd for gussng an ntal stmat for Kalman fltr typ stochastc atttud dtrmnaton mthods. hos a pror stmats for th stmator may b thn corrctd on an pandd stat vctor ncludng atttud paramtrs and bass. As anothr ampl for atttud-fr bas stmaton procdur, n (Crassds t al., 5; Huang and Jng, 8) t s provd that a full magntomtr calbraton can b prformd on-orbt durng typcal spaccraft msson-mod opratons by th us of ral-tm algorthms basd on both th Etndd Kalman Fltr (EKF) and Unscntd Kalman Fltr (UKF). Although ths studs can stmat th magntomtr charactrstcs such as bass, scal factors and nonorthogonalty corrctons, thy all dsrgard th atttud dynamcs of th satllt and f ths nformaton s possbl at any nstant, thr accuracs can b cdd. Hnc thy may b consdrd as a part of spaccraft mods whr atttud data s absnt. In (Styn, 995) both th magntomtr bass and scal factors ar stmatd by rcursv last squars mthod. h vctor of unnown paramtrs s formd of magntomtr bass and scal factors. Mthod s not atttud fr snc th atttud nowldg s ncssary to form th atttud matr. Howvr as a drawbac, usd atttud paramtrs ar not stmatd by th mthod tslf and th spaccraft dynamcs ar not tan nto consdraton n ths mthod. On th othr hand, a fw studs l (DaForno t al., 4; Ma, 5), whch handl magntomtr calbraton procss togthr wth th stmaton of atttud paramtrs, also do st. In (DaForno t al., 4) EKF s usd as a part of th stmaton schm. Va an UKF basd stmaton algorthm, whch has advantags ovr EKF such as absnc of Jacoban calculatons tc., ts accuracy may b surmountd. Ma and Jang solv magntomtr calbraton and atttud stmaton problm va two UKFs worng synchronously. As a Copyrght by th Intrnatonal Fdraton of Automatc Control (IFAC) 74

2 Prprnts of th 8th IFAC World Congrss Mlano (Italy) August 8 - Sptmbr, dsadvantag, ths approach rqurs a hgh computatonal ffort bcaus of two dstnct UKFs and t may b not sutabl for pco satllts whr th procssng capacty of th atttud computr s lmtd. In ths study a rconfgurabl unscntd Kalman fltr (UKF) basd algorthm s proposd for magntomtr calbraton and atttud paramtr stmaton. Algorthm s composd of two stags; n frst stag UKF stmats magntomtr bass and scal factors as wll as s atttud paramtrs of th satllt. Dffrntly from th stng algorthms, scal factors ar not tratd togthr wth th othr paramtrs as a part of th stat vctor; thr scal factors ar stmatd va a nwly proposd tnson for th UKF. At th scond stag whch starts aftr a convrgnc rul for th bas stmatons UKF rconfgurs tslf and rgardng th stmatd bass and scal factors, t runs for stmaton of only atttud paramtrs. Proposd algorthm s smulatd for atttud stmaton of a pco satllt whch has thr magntomtrs and thr rat gyros as masurmnt snsors.. MAHEMAICAL MODEL OF HE PICO SAELLIE If th nmatcs of th pco satllt s drvd n th bas of Eulr angls, thn th mathmatcal modl can b prssd wth a 9 dmnsonal systm vctor whch s mad of atttud Eulr angls (ϕ s th roll angl about as; θ s th ptch angl about y as; ψ s th yaw angl about z as) vctor, th body angular rat vctor wth rspct to th nrtal as fram, and th vctor formd of bas trms of magntomtrs. Hnc, For th frst stag whr bas trms of thr magntomtr ar stmatd; = ϕ θ ψ ω ω ω, y z bm b m b y m () z For th scond stag only atttud paramtrs ar stmatd; = ϕ θ ψ ω ω ω, () y z Hr subscrpt m for bas trms rprsnts magntomtr. Also for consstncy wth th furthr planatons, th body angular rat vctor wth rspct to th nrtal as fram should b statd sparatly as; ω = ω ω ω, BI y z whr ω BI s th angular vlocty vctor of th body fram wth rspct to th nrtal fram. Bsds, dynamc quatons of th satllt can b drvd by th us of th angular momntum consrvaton law; ω BI J = Ngg ωbi ( JωBI ), (3) dt whr J s th nrta matr conssts of prncpal momnts of nrta as J = dag ( J, J y, J z) and N gg s th vctor of gravty gradnt torqu affctng th satllt as (Shavat t al., 7); N ( J ) y Jz A3 A gg 33 μ Ngg 3 3 ( ) y = Jz J A3 A33. (4) r N gg ( J Jy) A3 A3 z Hr μ s th gravtatonal constant, r s th dstanc btwn th cntr of mass of th satllt and th Earth and A j rprsnts th corrspondng lmnt of th drcton cosn matr of; c( θ) c( ψ) c( θ) s( ψ) s( θ) ( ϕ) ( ψ) ( ϕ) ( θ) ( ψ) ( ϕ) ( ψ) ( ϕ) ( θ) ( ψ) ( ϕ) ( θ). (5) ( ϕ) ( ψ) + ( ϕ) ( θ) ( ψ) ( ϕ) ( ψ) + ( ϕ) ( θ) ( ψ) ( ϕ) ( θ) s ar th cosns and snus functons A= c s + s s c c c + s s s s c s s c s c s c c s s c c In matr A, ( ) c and ( ) succssvly. Knmatc quatons of moton of th pco satllt wth th Eulr angls can b gvn as ϕ s( ϕ) t( θ) c( ϕ) t( θ) p θ ( ϕ) ( ϕ = c s ) q. (6) ψ ( ϕ) / ( θ) ( ϕ) / ( θ) s c c c r Hr t ( ) stands for tangnt functon and p, q and r ar th componnts of ω BR vctor whch ndcats th angular vlocty of th body fram wth rspct to th rfrnc fram. ω BI and ω BR can b rlatd va, [ ] BR = BI + A. (7) ω ω ω whr ω dnots th angular vlocty of th orbt wth 3 rspct to th nrtal fram, found as = ( r ) / ω μ/. 3. MEASUREMEN SENSOR MODEL 3. h Magntomtr Modl As th satllt navgats along ts orbt, magntc fld vctor dffrs n a rlvant way wth th orbtal paramtrs. If thos paramtrs ar nown, thn, magntc fld tnsor vctor that affcts satllt can b shown as a functon of tm analytcally (Shavat t al., 7). Not that, ths trms ar obtand n th orbt rfrnc fram. M H() t = { cos( ωt) 3 cos( ε) sn( ) sn( ε) cos( ) cos( ωt) r sn( ωt) sn( ε) sn ( ωt )} (8) M H () t = 3 cos( ε) cos( ) + sn ( ε) sn ( ) cos ( ωt), (9) r M H3() t = { sn( ωt) 3 cos( ε) sn ( ) sn ( ε) cos( ) cos( ωt) r sn( ωt) sn( ε) sn( ωt )} () Hr M s th magntc dpol momnt of th Earth as 5 M = Wbm., μ s th Earth Gravtatonal 4 3 constant as μ = m / s, s th orbt nclnaton as = 97, ω s th spn rat of th Earth as 5 ω 7.9 rad / s =, ε s th magntc dpol tlt as ε =.7 and r s th dstanc btwn th cntr of mass of th satllt and th Earth as r = 6,98,4 m. hr onboard magntomtrs of pco satllt masurs th componnts of th magntc fld vctor n th body fram. hrfor for th masurmnt modl, whch charactrzs th masurmnts n th body fram, gand magntc fld 74

3 Prprnts of th 8th IFAC World Congrss Mlano (Italy) August 8 - Sptmbr, trms must b transformd by th us of drcton cosn matr, A. Ovrall masurmnt modl may b gvn as; H ( ϕθψ,,, ) ( ) t H t S Hy ( ϕ, θψ,, t) = A H() t + bm + η, () ( ϕθψ,,, ) 3 () Hz t H t whr, H () t, H ( t) and H3 ( t ) rprsnt th Earth magntc fld vctor componnts n th orbt fram as a ϕθψ,,, H ϕθψ,,, t and functon of tm, and H ( t ), y ( ) ( ϕθψ,,, ) H z t show th masurd Earth magntc fld vctor componnts n body fram as a functon of tm and varyng Eulr angls. Furthrmor, S s th dagonal scal matr as, S = dag( s s, s ), (), 33 b m s th magntomtr bas vctor as b = m bm bm y b m z and η s th zro man Gaussan wht nos wth th charactrstc of E ηη j = I 33σmδj. (3) I33s th dntty matr wth th dmnson of 3 3, σ m s th standard dvaton of ach magntomtr rror and δ j s th Kroncr symbol. 3. h Rat Gyro Modl hr rat gyros ar algnd through thr as, orthogonally to ach othr and thy supply drctly th angular rats of th body fram wth rspct to th nrtal fram. Hnc th modl for rat gyros can b gvn as; ωbi, mas = ωbi + η. (4) whr, ω BImas, s th masurd angular rats of th satllt, and η s th zro man Gaussan wht nos wth th charactrstc of E η η j = I 33σδ g j, (5) Hr, σ g s th standard dvaton of ach rat gyro random rror. 4. UKF BASED BIAS AND SCALE FACOR ESIMAION ALGORIHM 4. Unscntd Kalman Fltr In ordr to utlz Kalman fltr for nonlnar systms wthout any lnarzaton stp, th unscntd transform and so Unscntd Kalman Fltr s on of th tchnqus. UKF uss th unscntd transform, a dtrmnstc samplng tchnqu, to dtrmn a mnmal st of sampl ponts (or sgma ponts) from th a pror man and covaranc of th stat. hn, ths sgma ponts go through nonlnar transformaton. h postror man and th covaranc ar obtand from ths transformd sgma ponts (Julr t al., 995). As t s statd, UKF procdur bgns wth th dtrmnaton of n + sgma ponts wth a man of ˆ ( ) and a covaranc of P( ). For an n dmnsonal stat vctor, ths sgma ponts ar obtand by ( ) = ( ) (6) ˆ ( ) ( κ ) ( ) = ( ) + ( + ) ( ) + ( ) ˆ n κ P Q (7) ( ) ˆ + = ( ) ( + ) ( ) + ( ) n n P Q, (8) whr, ( ), ( ) and + n( ) ar sgma ponts, Q( ) s th procss nos covaranc matr, n s th stat numbr and κ s th scalng paramtr whch s usd for fn tunng and th hurstc s to chos that paramtr as n + κ = 3 (Julr t al., 995). Also, s gvn as = n. Nt stp of th UKF procss s transformng ach sgma pont by th us of systm dynamcs, ( + ) = ( ), f. (9) hn ths transformd valus ar utlzd for ganng th prdctd man and th covaranc (Crassds and Marly, 3; Son and Hajyv 9). ˆ ( ) + = κ ( + ) + n ( + ), () n + κ = { P( + ) = κ ( + ) ˆ( + ) ( + ) ˆ( + ) + κ n + n ( + ) ˆ( + ) ( + ) ˆ( + ). () = Hr, ˆ ( + ) s th prdctd man and P( + ) s th prdctd covaranc. Nonthlss, prdctd obsrvaton vctor s, ˆ ( ) + = κ ( + ) + n y y y ( + ), () n + κ = whr, y ( + ) = ( + ), ( ), h v. (3) Aftr that, obsrvaton covaranc matr s dtrmnd as, { P ( + ) = κ ( + ) ˆ( + ) ( + ) ˆ( + ) yy + κ y y y y n + n y ( ) ˆ( ) ( ) ˆ( ) y y y, (4) = whr nnovaton covaranc s P + = P + + R + (5) vv ( ) yy ( ) ( ) Hr v ( ) s th wht Gaussan masurmnt nos and t s non-corrlatd wth th procss nos. Bsds, R( +) s th masurmnt nos covaranc matr. h cross corrlaton matr can b obtand as, { P ( + ) = κ ( + ) ˆ( + ) ( + ) ˆ( + ) y + κ y y n + n ( ) ˆ( ) ( ) ˆ( ) y y. (6) = 743

4 Prprnts of th 8th IFAC World Congrss Mlano (Italy) August 8 - Sptmbr, Followng part s th updat phas of UKF algorthm. At that phas, frst by usng masurmnts, y( +), rsdual trm (or nnovaton squnc) s found as ( + ) = y ( + ) y ˆ ( + ), (7) and thn Kalman gan s computd va quaton of, ( ) y ( ) vv ( ) K + = P + P +. (8) At last, updatd stats and covaranc matr ar dtrmnd by, ˆ + + = ˆ + + K+ +, (9) ( ) ( ) ( ) ( ) ( + + ) = ( + ) ( + ) ( + ) vv ( + ) Hr, ˆ ( + + ) P( + + ) s th stmatd covaranc matr. P P K P K. (3) s th stmatd stat vctor and 4. Scal Factor Estmaton Whn th condton of th masurmnt systm dos not corrspond to th modl usd n th synthss of th fltr that mans th ncomng masurmnts ar scald and thy nd to b r-scald aftr stmaton of scal factors. If th systm oprats normally, th ral and th thortcal nnovaton covaranc matr valus match as n (3). ( + ) ( + ) ( + ) + ( + ), μ = Pyy R (3) + w j μw hr, μ w s th wdth of th movng wndow. Howvr, whn th ral masurmnts ar scald as n () th ral rror wll cd th thortcal on. In ths cas th matr bult of scal factors, S( ), s addd nto th algorthm as, ( + ) ( + ) = ( + ) + ( ) ( + ) ( ). μ = Pyy SR S (3) + w j μw hn, snc t s now that S( ) s dagonal (thr n no msalgnmnt/nonorthagonalty for magntomtrs) and also R( +) s dagonal by ts natur, t s possbl to rwrt (3), ( + ) ( + ) = ( + ) + ( ) ( ). μ + = Pyy S R (33) + w j μw Fnally, scal factor matr can b found as gvn blow, S( ) = ( + ) ( + ) Pyy ( + ) R ( + ). (34) μ w j= μw+ 4.3 Stoppng Rul for Bas Estmaton h followng stoppng rul may b ntroducd for th bas stmaton problm (Hajyv, 994; ): r = (bˆ -b ˆ ) D (bˆ - b ˆ ) = ε, (35) - - Δb - whr D Δb s th covaranc matr of th dscrpancy btwn two succssv bas stmats b ˆ and ˆb -, and ε s a prdtrmnd small numbr. Snc Kalman fltr stmats th paramtrs from a squnc of obsrvatons wth Gaussan tolrancs of th masurmnt rrors and systm nos, th fltr ylds an stmat wth an pctd valu qual to th stmatd quantty and a Gaussan dstrbuton functon. h dscrpancy bˆ -b ˆ - thn has a normal dstrbuton as wll, snc t s a lnar combnaton of two Gaussan random varabls. Wth ths consdratons n mnd t s nown that th statstc r has a χ dstrbuton wth n dgrs of frdom (n s th numbr of dmnsons of th bas vctor b ), and th thrshold valus of r can b found by dtrmnng th tabulatd valus of th χ dstrbuton for a gvn lvl of sgnfcanc. It s vdnt from rlaton (35) that th smallr th valu of r, th gratr wll b th consstncy of th stmats. Usually n th tstng of consstncy n such cass th lowr lmt of th confdnc ntrval must b qual to zro, and th uppr lmt s dtrmnd by th lvl of sgnfcanc α. o tst th consstncy of th stmats, th lvl of sgnfcanc α s adoptd, whch corrsponds to th confdnc coffcnt β = α. h thrshold χ β s spcfd n trms of ths probablty, usng th dstrbuton of th nvstgatd statstc r : { β } P χ < χ = β,< β <. (36) h stmaton procss s stoppd whn r < χ β, snc furthr obsrvatons yld nsgnfcant mprovmnt of th dntfd modl. If th quadratc form r s largr than or qual to th spcfd thrshold χ β, stmaton should b contnud. Nonthlss computaton of th covaranc matr of th dscrpancy btwn two succssv bas stmats b ˆ and ˆb - s nvstgatd n (Hajyv, ). h followng prsson for th covaranc matr D Δb s obtand n ths wor; D = P P. (37) Δb Dscrbd stoppng rul can b usd to ma a tmly dcson to stop th bas stmaton procss and t dos not rqur larg computatonal pndturs. 5. SIMULAIONS Smulatons ar ralzd for 5 sconds wth a samplng tm of Δ t =.sc. As an prmntal platform a cubsat modl s usd and th dagonal trms of th nrta matr s tan as J =. 3, J =. 3 y, J =.9 3 z, whr all non-dagonal trms ar zro snc th rotaton s about th 744

5 Prprnts of th 8th IFAC World Congrss Mlano (Italy) August 8 - Sptmbr, prncpal as of th satllt. Nonthlss th orbt of th satllt s a crcular orbt wth an alttud of ra = 55m. Othr orbt paramtrs ar sam as t s prsntd n th scton for th Magntomtr Modl (Scton 3.). Smulatons can b catgorzd n two: Estmaton of scal factors, bass and atttud paramtrs wthout stoppng rul for th fltr (wthout UKF rconfguraton) and stmaton ralzd by an UKF wth stoppng rul (wth UKF rconfguraton). Rsults gand va ths two dstnct smulaton scnaros ar compard n ordr to clarfy th ffctvnss of th proposd rconfgurabl UKF basd mthod. Frstly, va a sngl UKF algorthm t s possbl to calbrat th magntomtrs. In ths algorthm, stat vctor s formd of 9 stats; Eulr angls, angular rats and magntomtr bass. As wll scal factors ar stmatd va (34). Bas stmaton ampl s gvn n Fg.. As t s apparnt, bas trms ar stmatd accuratly. bz & bz a (tsla) rror(tsla) varanc (tsla ) Magntomtr bas bz Estmaton Kalman Estmaton Actual Valu tm(sc) Fg.. Estmaton of th bas for magntomtr algnd through z as. Bsds, n smulatons, scal factors ar mplmntd to th magntomtr masurmnts n ordr to smulat stmaton procss ( S ). In rturn, scal factors gvn at th rght-hand sd ( Ŝ ) ar stmatd va th proposd mthod (gvn rsults ar th man of scal factor stmaton sampls gand tll th bas stmaton stoppng nstant):..4 S =.3 Sˆ.356 =.5.55 As sn, stmaton rror for th scal factors s not hgh. Dffrnc btwn th actual valu and th Kalman stmaton s crtanly causd by th random procss wht nos of th masurmnts. On th othr hand, t s possbl to stmat also all othr paramtrs accuratly va ths algorthm. As sn n Fg. atttud paramtr stmaton rrors ar n accptabl lmts for a pco satllt. Not that smlar rsults hav bn obtand for all othr paramtrs. thta(dg) rror(dg) Ptch Angl Estmaton Kalman Estmaton Actual Valu varanc (dg ) tm(sc) Fg.. Ptch angl stmaton va UKF wthout rconfguraton. It s obvous that wth an UKF algorthm t s possbl to stmat all 9 paramtrs of th stat vctor () and also th scal factors prcsly. Howvr, a 9 dmnsonal stat vctor and th addtonal part to UKF for scal factor stmaton dmand a hgh computatonal ffort. In gnral, Kalman fltr s computatonal burdn ncrass n a drct rlaton wth th cub of th numbr of stats. Hnc, consttutng a stat vctor wth 9 stats mans mor computatonal burdn than a stat vctor wth 6 stats. hat s an mportant rdundancy for consdrably lmtd atttud computr of a pco satllt. As a rsult t s bttr to rconfgur fltr rgardng th scal factors and bass of magntomtrs and to procd th stmaton procdur by only stmatng Eulr angls and angular rats. As gvn n scton 4.3, a stoppng rul s proposd n ordr to rconfgur UKF aftr convrgnc of th bas stmatons to th actual valus. In that scond stag, prvously stmatd magntomtr bas trms and scal factors ar tan nto account by UKF algorthm and numbr of stats to b stmatd rducs to 6 as shown wth (). Ptch angl stmaton rsult s gvn n Fg.3. thta(dg) rror(dg) Ptch Angl Estmaton Kalman Estmaton Actual Valu varanc (dg ) -5 Bas Estmaton Stoppng Instant tm(sc) Fg. 3. Ptch angl stmaton va UKF wth rconfguraton. 745

6 Prprnts of th 8th IFAC World Congrss Mlano (Italy) August 8 - Sptmbr, Hr n Fg. 3 rconfguraton occurs at about 8 th scond and t s mard on fgur wth dashd ln. As t s apparnt from th stmaton varanc, rconfguraton brngs about an nhancmnt to fltr accuracy. So as to undrstand th ffct of rconfguraton mor clarly, absolut valus of rror ar tabulatd for two fltrs; UKF wth and wthout rconfguraton (abl ). abl. Comparson of Absolut Valus of Error for UKFs wth and wthout rconfguraton. Abs. Valus of Err. for UKF wth Rconfguraton Abs. Valus of Err. for UKF wthout Rconfguraton Par. s. 4s. s. 4s. ϕ(dg) θ(dg) ψ (dg) ω (dg/s) ωy (dg/s) ωz (dg/s) As sn from tabl rducng th sz of th stat vctor wthout any chang at masurmnt vctor brngs about an ncrmnt n th accuracy of stmaton of th stats. Nonthlss, a smallr stat vctor mans lowr computatonal burdn and rs for th dvrgnc of th fltr. 6. CONCLUSIONS In ths study a rconfgurabl unscntd Kalman fltr (UKF) basd algorthm s proposd for magntomtr calbraton and atttud paramtr stmaton. Smulaton rsults show that, t s possbl to stmat both, magntomtr bass and scal factors as wll as th atttud paramtrs va th proposd algorthm. Nvrthlss, proposd rconfgurabl UKF algorthm s advantagous bcaus smulaton rsults provs that rducng th sz of th stat vctor wthout any chang at masurmnt vctor brngs about an ncrmnt n th accuracy of stmaton of th stats. Also a smallr stat vctor mans lowr computatonal burdn and rs for th dvrgnc of th fltr. Snc magntomtr utlzaton has sgnfcanc spcally for pco satllt mssons, proposd algorthm may affct th msson prformanc and rlablty n a consdrabl dgr. ACKNOWLEDGMEN hs wor was supportd n part by UBIAK (h Scntfc and chnologcal Rsarch Councl of ury) undr Grant 8M53. REFERENCES Alonso, R. and M.D. Shustr (a). wo Stp: A Fast Robust Algorthm for Atttud- Indpndnt Magntomtr-Bas Dtrmnaton. h Journal of th Astronautcal Scncs, Vol.5, No.4, Alonso, R. and M.D. Shustr (b). Complt lnar atttud-ndpndnt magntomtr calbraton. h Journal of th Astronautcal Scncs, Vol.5, No.4, Alonso, R. and M.D. Shustr (3). Cntrng and Obsrvablty n Atttud-Indpndnt Magntomtr- Bas Dtrmnaton. h Journal of th Astronautcal Scncs, Vol.5, No., Crassds, J.L. and F.L. Marly (3). Unscntd Fltrng for Spaccraft Atttud Estmaton. Journal of Gudanc, Control and Dynamcs, Vol.6, No.4, Crassds, J.L., K. La, and R.R. Harman (5). Ral-tm atttud-ndpndnt thr-as magntomtr calbraton. Journal of Gudanc, Control, and Dynamcs, Vol.8, No., 5-. DaForno R. and t al. (4). Autonomous navgaton of mgsat : Atttud, snsor bas and scal factor stmaton by EKF and magntomtr-only masurmnt. Proc. nd AIAA Intrnatonal Communcatons Satllt Systms Confrnc and Ehbt, Calforna, USA. Hajyv, Ch.M.(994). Som Problms of Dynamc Systms Paramtrcal Idntfcaton. In Dgst of Frst Asan Control Confrnc, July 7-3 (994), Japan,. Vol.: Hajyv, Ch.M.(). Stoppng ruls formaton and faults dtcton n paramtrc dntfcaton problms. IMchE Procdngs, Part I, J.Systms & Control Engnrng, Vol.5 (I4): Huang, L. and W. Jng (8). Atttud-ndpndnt gomagntc navgaton usng onboard complt thras magntomtr calbraton. Proc. 8 IEEE Arospac Confrnc, Montana, USA, -7. Julr, S.J., J.K. Uhlmann and H.F. Durrant-Whyt (995). A Nw Approach for Fltrng Nonlnar Systms. Procdngs of Amrcan Control Confrnc, Vol.3, Km, E. and H. Bang (7). Bas stmaton of magntomtr usng gntc algorthm. Proc. Intrnatonal Confrnc of Control Automaton and Systm 7, Soul, Kora, Ma, G.F. and X.Y. Jang (5). Unscntd Kalman Fltr for Spaccraft Atttud Estmaton and Calbraton Usng Magntomtr Masurmnts, Proc. 4 th Intrnatonal Confrnc on Machn Larnng and Cybrntcs, Guangzhou, Chna, Shavat P., Q. Gong and I.M. Ross (7). NPSA I Paramtr Estmaton Usng Unscntd Kalman Fltr. Procdngs of 7 Amrcan Control Confrnc, , IEEE, Nw Yor, USA. Styn, W.H. (995). A Mult-mod Atttud Dtarmnaton and Control Systm for Small Satllts. Dssrtaton Prsntd for h Dgr of Doctor of Phlosophy, Unvrsty of Stllnbosch, South Afrca. Son, H.E. and Ch. Hajyv (9). UKF for th Idntfcaton of th Pco Satllt Atttud Dynamcs Paramtrs and th Etrnal orqus on IMU and Magntomtr Masurmnts. Proc. 4th Intrnatonal Confrnc on Rcnt Advancs n Spac chnologs, Istanbul, ury, Wrtz, J.R. (988). Spaccraft Atttud Dtrmnaton and Control, Kluwr Acadmc Publshrs, Dordrcht, Holland. 746

COMPLEX NUMBER PAIRWISE COMPARISON AND COMPLEX NUMBER AHP

COMPLEX NUMBER PAIRWISE COMPARISON AND COMPLEX NUMBER AHP ISAHP 00, Bal, Indonsa, August -9, 00 COMPLEX NUMBER PAIRWISE COMPARISON AND COMPLEX NUMBER AHP Chkako MIYAKE, Kkch OHSAWA, Masahro KITO, and Masaak SHINOHARA Dpartmnt of Mathmatcal Informaton Engnrng

More information

Review - Probabilistic Classification

Review - Probabilistic Classification Mmoral Unvrsty of wfoundland Pattrn Rcognton Lctur 8 May 5, 6 http://www.ngr.mun.ca/~charlsr Offc Hours: Tusdays Thursdays 8:3-9:3 PM E- (untl furthr notc) Gvn lablld sampls { ɛc,,,..., } {. Estmat Rvw

More information

Outlier-tolerant parameter estimation

Outlier-tolerant parameter estimation Outlr-tolrant paramtr stmaton Baysan thods n physcs statstcs machn larnng and sgnal procssng (SS 003 Frdrch Fraundorfr fraunfr@cg.tu-graz.ac.at Computr Graphcs and Vson Graz Unvrsty of Tchnology Outln

More information

8-node quadrilateral element. Numerical integration

8-node quadrilateral element. Numerical integration Fnt Elmnt Mthod lctur nots _nod quadrlatral lmnt Pag of 0 -nod quadrlatral lmnt. Numrcal ntgraton h tchnqu usd for th formulaton of th lnar trangl can b formall tndd to construct quadrlatral lmnts as wll

More information

The Hyperelastic material is examined in this section.

The Hyperelastic material is examined in this section. 4. Hyprlastcty h Hyprlastc matral s xad n ths scton. 4..1 Consttutv Equatons h rat of chang of ntrnal nrgy W pr unt rfrnc volum s gvn by th strss powr, whch can b xprssd n a numbr of dffrnt ways (s 3.7.6):

More information

CHAPTER 7d. DIFFERENTIATION AND INTEGRATION

CHAPTER 7d. DIFFERENTIATION AND INTEGRATION CHAPTER 7d. DIFFERENTIATION AND INTEGRATION A. J. Clark School o Engnrng Dpartmnt o Cvl and Envronmntal Engnrng by Dr. Ibrahm A. Assakka Sprng ENCE - Computaton Mthods n Cvl Engnrng II Dpartmnt o Cvl and

More information

ST 524 NCSU - Fall 2008 One way Analysis of variance Variances not homogeneous

ST 524 NCSU - Fall 2008 One way Analysis of variance Variances not homogeneous ST 54 NCSU - Fall 008 On way Analyss of varanc Varancs not homognous On way Analyss of varanc Exampl (Yandll, 997) A plant scntst masurd th concntraton of a partcular vrus n plant sap usng ELISA (nzym-lnkd

More information

A Note on Estimability in Linear Models

A Note on Estimability in Linear Models Intrnatonal Journal of Statstcs and Applcatons 2014, 4(4): 212-216 DOI: 10.5923/j.statstcs.20140404.06 A Not on Estmablty n Lnar Modls S. O. Adymo 1,*, F. N. Nwob 2 1 Dpartmnt of Mathmatcs and Statstcs,

More information

Analyzing Frequencies

Analyzing Frequencies Frquncy (# ndvduals) Frquncy (# ndvduals) /3/16 H o : No dffrnc n obsrvd sz frquncs and that prdctd by growth modl How would you analyz ths data? 15 Obsrvd Numbr 15 Expctd Numbr from growth modl 1 1 5

More information

EEC 686/785 Modeling & Performance Evaluation of Computer Systems. Lecture 12

EEC 686/785 Modeling & Performance Evaluation of Computer Systems. Lecture 12 EEC 686/785 Modlng & Prformanc Evaluaton of Computr Systms Lctur Dpartmnt of Elctrcal and Computr Engnrng Clvland Stat Unvrsty wnbng@.org (basd on Dr. Ra Jan s lctur nots) Outln Rvw of lctur k r Factoral

More information

ATTITUDE ESTIMATION PROCESS OF THE SENSING REMOTE SATELLITE CBERS-2 WITH UNSCENTED KALMAN FILTER AND QUATERNION INCREMENTAL, USING REAL DATA

ATTITUDE ESTIMATION PROCESS OF THE SENSING REMOTE SATELLITE CBERS-2 WITH UNSCENTED KALMAN FILTER AND QUATERNION INCREMENTAL, USING REAL DATA AIUDE ESIMAION PROCESS OF HE SENSING REMOE SAELLIE CBERS- WIH UNSCENED KALMAN FILER AND QUAERNION INCREMENAL, USING REAL DAA R. V. Garca (, H. K. Kua (, M. C. Zanard (3, and N. F. O. Matos (4 ( Unvrsdad

More information

2. Grundlegende Verfahren zur Übertragung digitaler Signale (Zusammenfassung) Informationstechnik Universität Ulm

2. Grundlegende Verfahren zur Übertragung digitaler Signale (Zusammenfassung) Informationstechnik Universität Ulm . Grundlgnd Vrfahrn zur Übrtragung dgtalr Sgnal (Zusammnfassung) wt Dc. 5 Transmsson of Dgtal Sourc Sgnals Sourc COD SC COD MOD MOD CC dg RF s rado transmsson mdum Snk DC SC DC CC DM dg DM RF g physcal

More information

The Fourier Transform

The Fourier Transform /9/ Th ourr Transform Jan Baptst Josph ourr 768-83 Effcnt Data Rprsntaton Data can b rprsntd n many ways. Advantag usng an approprat rprsntaton. Eampls: osy ponts along a ln Color spac rd/grn/blu v.s.

More information

Lucas Test is based on Euler s theorem which states that if n is any integer and a is coprime to n, then a φ(n) 1modn.

Lucas Test is based on Euler s theorem which states that if n is any integer and a is coprime to n, then a φ(n) 1modn. Modul 10 Addtonal Topcs 10.1 Lctur 1 Prambl: Dtrmnng whthr a gvn ntgr s prm or compost s known as prmalty tstng. Thr ar prmalty tsts whch mrly tll us whthr a gvn ntgr s prm or not, wthout gvng us th factors

More information

Maneuvering Target Tracking Using Current Statistical Model Based Adaptive UKF for Wireless Sensor Network

Maneuvering Target Tracking Using Current Statistical Model Based Adaptive UKF for Wireless Sensor Network Journal of Communcatons Vol. 0, No. 8, August 05 Manuvrng argt racng Usng Currnt Statstcal Modl Basd Adaptv UKF for Wrlss Snsor Ntwor Xaojun Png,, Kuntao Yang, and Chang Lu School of Optcal and Elctronc

More information

Maneuvering Target Tracking Using Current Statistical Model Based Adaptive UKF for Wireless Sensor Network

Maneuvering Target Tracking Using Current Statistical Model Based Adaptive UKF for Wireless Sensor Network Journal of Communcatons Vol. 0, No. 8, August 05 Manuvrng argt racng Usng Currnt Statstcal Modl Basd Adaptv UKF for Wrlss Snsor Ntwor Xaojun Png,, Kuntao Yang, and Chang Lu School of Optcal and Elctronc

More information

Soft k-means Clustering. Comp 135 Machine Learning Computer Science Tufts University. Mixture Models. Mixture of Normals in 1D

Soft k-means Clustering. Comp 135 Machine Learning Computer Science Tufts University. Mixture Models. Mixture of Normals in 1D Comp 35 Machn Larnng Computr Scnc Tufts Unvrsty Fall 207 Ron Khardon Th EM Algorthm Mxtur Modls Sm-Suprvsd Larnng Soft k-mans Clustrng ck k clustr cntrs : Assocat xampls wth cntrs p,j ~~ smlarty b/w cntr

More information

A NEW GENERALISATION OF SAM-SOLAI S MULTIVARIATE ADDITIVE GAMMA DISTRIBUTION*

A NEW GENERALISATION OF SAM-SOLAI S MULTIVARIATE ADDITIVE GAMMA DISTRIBUTION* A NEW GENERALISATION OF SAM-SOLAI S MULTIVARIATE ADDITIVE GAMMA DISTRIBUTION* Dr. G.S. Davd Sam Jayakumar, Assstant Profssor, Jamal Insttut of Managmnt, Jamal Mohamd Collg, Truchraall 620 020, South Inda,

More information

Jones vector & matrices

Jones vector & matrices Jons vctor & matrcs PY3 Colást na hollscol Corcagh, Ér Unvrst Collg Cork, Irland Dpartmnt of Phscs Matr tratmnt of polarzaton Consdr a lght ra wth an nstantanous -vctor as shown k, t ˆ k, t ˆ k t, o o

More information

Journal of Theoretical and Applied Information Technology 10 th January Vol. 47 No JATIT & LLS. All rights reserved.

Journal of Theoretical and Applied Information Technology 10 th January Vol. 47 No JATIT & LLS. All rights reserved. Journal o Thortcal and Appld Inormaton Tchnology th January 3. Vol. 47 No. 5-3 JATIT & LLS. All rghts rsrvd. ISSN: 99-8645 www.att.org E-ISSN: 87-395 RESEARCH ON PROPERTIES OF E-PARTIAL DERIVATIVE OF LOGIC

More information

10/7/14. Mixture Models. Comp 135 Introduction to Machine Learning and Data Mining. Maximum likelihood estimation. Mixture of Normals in 1D

10/7/14. Mixture Models. Comp 135 Introduction to Machine Learning and Data Mining. Maximum likelihood estimation. Mixture of Normals in 1D Comp 35 Introducton to Machn Larnng and Data Mnng Fall 204 rofssor: Ron Khardon Mxtur Modls Motvatd by soft k-mans w dvlopd a gnratv modl for clustrng. Assum thr ar k clustrs Clustrs ar not rqurd to hav

More information

Grand Canonical Ensemble

Grand Canonical Ensemble Th nsmbl of systms mmrsd n a partcl-hat rsrvor at constant tmpratur T, prssur P, and chmcal potntal. Consdr an nsmbl of M dntcal systms (M =,, 3,...M).. Thy ar mutually sharng th total numbr of partcls

More information

SCITECH Volume 5, Issue 1 RESEARCH ORGANISATION November 17, 2015

SCITECH Volume 5, Issue 1 RESEARCH ORGANISATION November 17, 2015 Journal of Informaton Scncs and Computng Tchnologs(JISCT) ISSN: 394-966 SCITECH Volum 5, Issu RESEARCH ORGANISATION Novmbr 7, 5 Journal of Informaton Scncs and Computng Tchnologs www.sctcrsarch.com/journals

More information

Physics of Very High Frequency (VHF) Capacitively Coupled Plasma Discharges

Physics of Very High Frequency (VHF) Capacitively Coupled Plasma Discharges Physcs of Vry Hgh Frquncy (VHF) Capactvly Coupld Plasma Dschargs Shahd Rauf, Kallol Bra, Stv Shannon, and Kn Collns Appld Matrals, Inc., Sunnyval, CA AVS 54 th Intrnatonal Symposum Sattl, WA Octobr 15-19,

More information

Chapter 6 Student Lecture Notes 6-1

Chapter 6 Student Lecture Notes 6-1 Chaptr 6 Studnt Lctur Nots 6-1 Chaptr Goals QM353: Busnss Statstcs Chaptr 6 Goodnss-of-Ft Tsts and Contngncy Analyss Aftr compltng ths chaptr, you should b abl to: Us th ch-squar goodnss-of-ft tst to dtrmn

More information

??? Dynamic Causal Modelling for M/EEG. Electroencephalography (EEG) Dynamic Causal Modelling. M/EEG analysis at sensor level. time.

??? Dynamic Causal Modelling for M/EEG. Electroencephalography (EEG) Dynamic Causal Modelling. M/EEG analysis at sensor level. time. Elctroncphalography EEG Dynamc Causal Modllng for M/EEG ampltud μv tm ms tral typ 1 tm channls channls tral typ 2 C. Phllps, Cntr d Rchrchs du Cyclotron, ULg, Blgum Basd on slds from: S. Kbl M/EEG analyss

More information

Group Codes Define Over Dihedral Groups of Small Order

Group Codes Define Over Dihedral Groups of Small Order Malaysan Journal of Mathmatcal Scncs 7(S): 0- (0) Spcal Issu: Th rd Intrnatonal Confrnc on Cryptology & Computr Scurty 0 (CRYPTOLOGY0) MALAYSIA JOURAL OF MATHEMATICAL SCIECES Journal hompag: http://nspm.upm.du.my/ournal

More information

THREE DIMENSIONAL GEOMETRY MAINTENANCE FOR FORMATION FLYING ON ELLIPTIC ORBITS

THREE DIMENSIONAL GEOMETRY MAINTENANCE FOR FORMATION FLYING ON ELLIPTIC ORBITS HREE DIMENSIONAL GEOMERY MAINENANCE FOR FORMAION FLYING ON ELLIPIC ORBIS akanao SAIKI ), Koch NASUME ) and Jun chro KAWAGUCHI ) ABSRAC ) Mtsubsh Havy Industrs, Ltd. ) Mtsubsh Elctrc Co. ) Japan Arospac

More information

A Probabilistic Characterization of Simulation Model Uncertainties

A Probabilistic Characterization of Simulation Model Uncertainties A Proalstc Charactrzaton of Sulaton Modl Uncrtants Vctor Ontvros Mohaad Modarrs Cntr for Rsk and Rlalty Unvrsty of Maryland 1 Introducton Thr s uncrtanty n odl prdctons as wll as uncrtanty n xprnts Th

More information

COMPARISON OF L1 C/A L2C COMBINED ACQUISITION TECHNIQUES

COMPARISON OF L1 C/A L2C COMBINED ACQUISITION TECHNIQUES COMPARION OF C/A C COMBINE ACQUIITION TECHNIQUE Cyrll Grnot, Kyl O Kf and Gérard achapll Poston, ocaton and Navgaton PAN Rsarch Group partmnt of Gomatcs Engnrng, Unvrsty of Calgary chulch chool of Engnrng

More information

Decentralized Adaptive Control and the Possibility of Utilization of Networked Control System

Decentralized Adaptive Control and the Possibility of Utilization of Networked Control System Dcntralzd Adaptv Control and th Possblty of Utlzaton of Ntworkd Control Systm MARIÁN ÁRNÍK, JÁN MURGAŠ Slovak Unvrsty of chnology n Bratslava Faculty of Elctrcal Engnrng and Informaton chnology Insttut

More information

Naresuan University Journal: Science and Technology 2018; (26)1

Naresuan University Journal: Science and Technology 2018; (26)1 Narsuan Unvrsty Journal: Scnc and Tchnology 018; (6)1 Th Dvlopmnt o a Corrcton Mthod or Ensurng a Contnuty Valu o Th Ch-squar Tst wth a Small Expctd Cll Frquncy Kajta Matchma 1 *, Jumlong Vongprasrt and

More information

September 27, Introduction to Ordinary Differential Equations. ME 501A Seminar in Engineering Analysis Page 1. Outline

September 27, Introduction to Ordinary Differential Equations. ME 501A Seminar in Engineering Analysis Page 1. Outline Introucton to Ornar Dffrntal Equatons Sptmbr 7, 7 Introucton to Ornar Dffrntal Equatons Larr artto Mchancal Engnrng AB Smnar n Engnrng Analss Sptmbr 7, 7 Outln Rvw numrcal solutons Bascs of ffrntal quatons

More information

Study of Dynamic Aperture for PETRA III Ring K. Balewski, W. Brefeld, W. Decking, Y. Li DESY

Study of Dynamic Aperture for PETRA III Ring K. Balewski, W. Brefeld, W. Decking, Y. Li DESY Stud of Dnamc Aprtur for PETRA III Rng K. Balws, W. Brfld, W. Dcng, Y. L DESY FLS6 Hamburg PETRA III Yong-Jun L t al. Ovrvw Introducton Dnamcs of dampng wgglrs hoc of machn tuns, and optmzaton of stupol

More information

Three-Node Euler-Bernoulli Beam Element Based on Positional FEM

Three-Node Euler-Bernoulli Beam Element Based on Positional FEM Avalabl onln at www.scncdrct.com Procda Engnrng 9 () 373 377 Intrnatonal Workshop on Informaton and Elctroncs Engnrng (IWIEE) Thr-Nod Eulr-Brnoull Bam Elmnt Basd on Postonal FEM Lu Jan a *,b, Zhou Shnj

More information

Α complete processing methodology for 3D monitoring using GNSS receivers

Α complete processing methodology for 3D monitoring using GNSS receivers 7-5-5 NATIONA TECHNICA UNIVERSITY OF ATHENS SCHOO OF RURA AND SURVEYING ENGINEERING DEPARTMENT OF TOPOGRAPHY AORATORY OF GENERA GEODESY Α complt procssng mthodology for D montorng usng GNSS rcvrs Gorg

More information

Observer Bias and Reliability By Xunchi Pu

Observer Bias and Reliability By Xunchi Pu Obsrvr Bias and Rliability By Xunchi Pu Introduction Clarly all masurmnts or obsrvations nd to b mad as accuratly as possibl and invstigators nd to pay carful attntion to chcking th rliability of thir

More information

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches. Subjct Chmistry Papr No and Titl Modul No and Titl Modul Tag 8/ Physical Spctroscopy / Brakdown of th Born-Oppnhimr approximation. Slction ruls for rotational-vibrational transitions. P, R branchs. CHE_P8_M

More information

Relate p and T at equilibrium between two phases. An open system where a new phase may form or a new component can be added

Relate p and T at equilibrium between two phases. An open system where a new phase may form or a new component can be added 4.3, 4.4 Phas Equlbrum Dtrmn th slops of th f lns Rlat p and at qulbrum btwn two phass ts consdr th Gbbs functon dg η + V Appls to a homognous systm An opn systm whr a nw phas may form or a nw componnt

More information

Folding of Regular CW-Complexes

Folding of Regular CW-Complexes Ald Mathmatcal Scncs, Vol. 6,, no. 83, 437-446 Foldng of Rgular CW-Comlxs E. M. El-Kholy and S N. Daoud,3. Dartmnt of Mathmatcs, Faculty of Scnc Tanta Unvrsty,Tanta,Egyt. Dartmnt of Mathmatcs, Faculty

More information

Ερωτήσεις και ασκησεις Κεφ. 10 (για μόρια) ΠΑΡΑΔΟΣΗ 29/11/2016. (d)

Ερωτήσεις και ασκησεις Κεφ. 10 (για μόρια) ΠΑΡΑΔΟΣΗ 29/11/2016. (d) Ερωτήσεις και ασκησεις Κεφ 0 (για μόρια ΠΑΡΑΔΟΣΗ 9//06 Th coffcnt A of th van r Waals ntracton s: (a A r r / ( r r ( (c a a a a A r r / ( r r ( a a a a A r r / ( r r a a a a A r r / ( r r 4 a a a a 0 Th

More information

ON THE COMPLEXITY OF K-STEP AND K-HOP DOMINATING SETS IN GRAPHS

ON THE COMPLEXITY OF K-STEP AND K-HOP DOMINATING SETS IN GRAPHS MATEMATICA MONTISNIRI Vol XL (2017) MATEMATICS ON TE COMPLEXITY OF K-STEP AN K-OP OMINATIN SETS IN RAPS M FARAI JALALVAN AN N JAFARI RA partmnt of Mathmatcs Shahrood Unrsty of Tchnology Shahrood Iran Emals:

More information

Transient Multiexponential Data Analysis Using A Combination of ARMA and ECD Methods

Transient Multiexponential Data Analysis Using A Combination of ARMA and ECD Methods 4 th Intrnatonal Confrnc on Computr Engnrng and Tchnology (ICCET IPCIT vol.4 ( ( IACIT Prss, ngapor Transnt ultxponntal Data Analyss Usng A Combnaton of ARA and ECD thods Abdussamad U. Jba, omoh-jmoh E.

More information

A Self-adaptive open loop architecture for weak GNSS signal tracking

A Self-adaptive open loop architecture for weak GNSS signal tracking NTERNATONAL JOURNAL OF CRCUTS, SYSTEMS AND SGNAL PROCESSNG Volum 8, 014 A Slf-adaptv opn loop archtctur for wa GNSS sgnal tracng Ao Png, Gang Ou, Janghong Sh Abstract An FFT-basd opn loop carrr tracng

More information

From Structural Analysis to FEM. Dhiman Basu

From Structural Analysis to FEM. Dhiman Basu From Structural Analyss to FEM Dhman Basu Acknowldgmnt Followng txt books wr consultd whl prparng ths lctur nots: Znkwcz, OC O.C. andtaylor Taylor, R.L. (000). Th FntElmnt Mthod, Vol. : Th Bass, Ffth dton,

More information

Heisenberg Model. Sayed Mohammad Mahdi Sadrnezhaad. Supervisor: Prof. Abdollah Langari

Heisenberg Model. Sayed Mohammad Mahdi Sadrnezhaad. Supervisor: Prof. Abdollah Langari snbrg Modl Sad Mohammad Mahd Sadrnhaad Survsor: Prof. bdollah Langar bstract: n ths rsarch w tr to calculat analtcall gnvalus and gnvctors of fnt chan wth ½-sn artcls snbrg modl. W drov gnfuctons for closd

More information

Addition of angular momentum

Addition of angular momentum Addition of angular momntum April, 0 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat th

More information

Optimal Ordering Policy in a Two-Level Supply Chain with Budget Constraint

Optimal Ordering Policy in a Two-Level Supply Chain with Budget Constraint Optmal Ordrng Polcy n a Two-Lvl Supply Chan wth Budgt Constrant Rasoul aj Alrza aj Babak aj ABSTRACT Ths papr consdrs a two- lvl supply chan whch consst of a vndor and svral rtalrs. Unsatsfd dmands n rtalrs

More information

Addition of angular momentum

Addition of angular momentum Addition of angular momntum April, 07 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat

More information

Consider a system of 2 simultaneous first order linear equations

Consider a system of 2 simultaneous first order linear equations Soluon of sysms of frs ordr lnar quaons onsdr a sysm of smulanous frs ordr lnar quaons a b c d I has h alrna mar-vcor rprsnaon a b c d Or, n shorhand A, f A s alrady known from con W know ha h abov sysm

More information

SPECTRUM ESTIMATION (2)

SPECTRUM ESTIMATION (2) SPECTRUM ESTIMATION () PARAMETRIC METHODS FOR POWER SPECTRUM ESTIMATION Gnral consdraton of aramtrc modl sctrum stmaton: Autorgrssv sctrum stmaton: A. Th autocorrlaton mthod B. Th covaranc mthod C. Modfd

More information

horizontal force output data block Hankel matrix transfer function complex frequency response function impedance matrix

horizontal force output data block Hankel matrix transfer function complex frequency response function impedance matrix AMBENT VBRATON Nomnclatur 1 NOMENCLATURE a a acclraton, coffcnt dmnsonlss corrcton factor A, B, C, D dscrt-tm stat sac modl b coffcnt c damng, stffnss C constant valu, damng matrx d damtr d dyn d stat

More information

Search sequence databases 3 10/25/2016

Search sequence databases 3 10/25/2016 Sarch squnc databass 3 10/25/2016 Etrm valu distribution Ø Suppos X is a random variabl with probability dnsity function p(, w sampl a larg numbr S of indpndnt valus of X from this distribution for an

More information

ACOUSTIC WAVE EQUATION. Contents INTRODUCTION BULK MODULUS AND LAMÉ S PARAMETERS

ACOUSTIC WAVE EQUATION. Contents INTRODUCTION BULK MODULUS AND LAMÉ S PARAMETERS ACOUSTIC WAE EQUATION Contnts INTRODUCTION BULK MODULUS AND LAMÉ S PARAMETERS INTRODUCTION As w try to vsualz th arth ssmcally w mak crtan physcal smplfcatons that mak t asr to mak and xplan our obsrvatons.

More information

Economics 600: August, 2007 Dynamic Part: Problem Set 5. Problems on Differential Equations and Continuous Time Optimization

Economics 600: August, 2007 Dynamic Part: Problem Set 5. Problems on Differential Equations and Continuous Time Optimization THE UNIVERSITY OF MARYLAND COLLEGE PARK, MARYLAND Economcs 600: August, 007 Dynamc Part: Problm St 5 Problms on Dffrntal Equatons and Contnuous Tm Optmzaton Quston Solv th followng two dffrntal quatons.

More information

Phys 774: Nonlinear Spectroscopy: SHG and Raman Scattering

Phys 774: Nonlinear Spectroscopy: SHG and Raman Scattering Last Lcturs: Polaraton of Elctromagntc Wavs Phys 774: Nonlnar Spctroscopy: SHG and Scattrng Gnral consdraton of polaraton Jons Formalsm How Polarrs work Mullr matrcs Stoks paramtrs Poncar sphr Fall 7 Polaraton

More information

167 T componnt oftforc on atom B can b drvd as: F B =, E =,K (, ) (.2) wr w av usd 2 = ( ) =2 (.3) T scond drvatv: 2 E = K (, ) = K (1, ) + 3 (.4).2.2

167 T componnt oftforc on atom B can b drvd as: F B =, E =,K (, ) (.2) wr w av usd 2 = ( ) =2 (.3) T scond drvatv: 2 E = K (, ) = K (1, ) + 3 (.4).2.2 166 ppnd Valnc Forc Flds.1 Introducton Valnc forc lds ar usd to dscrb ntra-molcular ntractons n trms of 2-body, 3-body, and 4-body (and gr) ntractons. W mplmntd many popular functonal forms n our program..2

More information

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012 Th van dr Waals intraction D. E. Sopr 2 Univrsity of Orgon 20 pril 202 Th van dr Waals intraction is discussd in Chaptr 5 of J. J. Sakurai, Modrn Quantum Mchanics. Hr I tak a look at it in a littl mor

More information

VISUALIZATION OF DIFFERENTIAL GEOMETRY UDC 514.7(045) : : Eberhard Malkowsky 1, Vesna Veličković 2

VISUALIZATION OF DIFFERENTIAL GEOMETRY UDC 514.7(045) : : Eberhard Malkowsky 1, Vesna Veličković 2 FACTA UNIVERSITATIS Srs: Mchancs, Automatc Control Robotcs Vol.3, N o, 00, pp. 7-33 VISUALIZATION OF DIFFERENTIAL GEOMETRY UDC 54.7(045)54.75.6:59.688:59.673 Ebrhard Malkowsky, Vsna Vlčkovć Dpartmnt of

More information

Lecture 23 APPLICATIONS OF FINITE ELEMENT METHOD TO SCALAR TRANSPORT PROBLEMS

Lecture 23 APPLICATIONS OF FINITE ELEMENT METHOD TO SCALAR TRANSPORT PROBLEMS COMPUTTION FUID DYNMICS: FVM: pplcatons to Scalar Transport Prolms ctur 3 PPICTIONS OF FINITE EEMENT METHOD TO SCR TRNSPORT PROBEMS 3. PPICTION OF FEM TO -D DIFFUSION PROBEM Consdr th stady stat dffuson

More information

External Equivalent. EE 521 Analysis of Power Systems. Chen-Ching Liu, Boeing Distinguished Professor Washington State University

External Equivalent. EE 521 Analysis of Power Systems. Chen-Ching Liu, Boeing Distinguished Professor Washington State University xtrnal quvalnt 5 Analyss of Powr Systms Chn-Chng Lu, ong Dstngushd Profssor Washngton Stat Unvrsty XTRNAL UALNT ach powr systm (ara) s part of an ntrconnctd systm. Montorng dvcs ar nstalld and data ar

More information

Stress-Based Finite Element Methods for Dynamics Analysis of Euler-Bernoulli Beams with Various Boundary Conditions

Stress-Based Finite Element Methods for Dynamics Analysis of Euler-Bernoulli Beams with Various Boundary Conditions 9 Strss-Basd Fnt Elmnt Mthods for Dynamcs Analyss of Eulr-Brnoull Bams wth Varous Boundary Condtons Abstract In ths rsarch, two strss-basd fnt lmnt mthods ncludng th curvatur-basd fnt lmnt mthod (CFE)

More information

Basic Electrical Engineering for Welding [ ] --- Introduction ---

Basic Electrical Engineering for Welding [ ] --- Introduction --- Basc Elctrcal Engnrng for Wldng [] --- Introducton --- akayosh OHJI Profssor Ertus, Osaka Unrsty Dr. of Engnrng VIUAL WELD CO.,LD t-ohj@alc.co.jp OK 15 Ex. Basc A.C. crcut h fgurs n A-group show thr typcal

More information

Supplementary Materials

Supplementary Materials 6 Supplmntary Matrials APPENDIX A PHYSICAL INTERPRETATION OF FUEL-RATE-SPEED FUNCTION A truck running on a road with grad/slop θ positiv if moving up and ngativ if moving down facs thr rsistancs: arodynamic

More information

Estimation of apparent fraction defective: A mathematical approach

Estimation of apparent fraction defective: A mathematical approach Availabl onlin at www.plagiarsarchlibrary.com Plagia Rsarch Library Advancs in Applid Scinc Rsarch, 011, (): 84-89 ISSN: 0976-8610 CODEN (USA): AASRFC Estimation of apparnt fraction dfctiv: A mathmatical

More information

te Finance (4th Edition), July 2017.

te Finance (4th Edition), July 2017. Appndx Chaptr. Tchncal Background Gnral Mathmatcal and Statstcal Background Fndng a bas: 3 2 = 9 3 = 9 1 /2 x a = b x = b 1/a A powr of 1 / 2 s also quvalnt to th squar root opraton. Fndng an xponnt: 3

More information

Construction of asymmetric orthogonal arrays of strength three via a replacement method

Construction of asymmetric orthogonal arrays of strength three via a replacement method isid/ms/26/2 Fbruary, 26 http://www.isid.ac.in/ statmath/indx.php?modul=prprint Construction of asymmtric orthogonal arrays of strngth thr via a rplacmnt mthod Tian-fang Zhang, Qiaoling Dng and Alok Dy

More information

ph People Grade Level: basic Duration: minutes Setting: classroom or field site

ph People Grade Level: basic Duration: minutes Setting: classroom or field site ph Popl Adaptd from: Whr Ar th Frogs? in Projct WET: Curriculum & Activity Guid. Bozman: Th Watrcours and th Council for Environmntal Education, 1995. ph Grad Lvl: basic Duration: 10 15 minuts Stting:

More information

Lecture 14. Relic neutrinos Temperature at neutrino decoupling and today Effective degeneracy factor Neutrino mass limits Saha equation

Lecture 14. Relic neutrinos Temperature at neutrino decoupling and today Effective degeneracy factor Neutrino mass limits Saha equation Lctur Rlc nutrnos mpratur at nutrno dcoupln and today Effctv dnracy factor Nutrno mass lmts Saha quaton Physcal Cosmoloy Lnt 005 Rlc Nutrnos Nutrnos ar wakly ntractn partcls (lptons),,,,,,, typcal ractons

More information

Background: We have discussed the PIB, HO, and the energy of the RR model. In this chapter, the H-atom, and atomic orbitals.

Background: We have discussed the PIB, HO, and the energy of the RR model. In this chapter, the H-atom, and atomic orbitals. Chaptr 7 Th Hydrogn Atom Background: W hav discussd th PIB HO and th nrgy of th RR modl. In this chaptr th H-atom and atomic orbitals. * A singl particl moving undr a cntral forc adoptd from Scott Kirby

More information

An Overview of Markov Random Field and Application to Texture Segmentation

An Overview of Markov Random Field and Application to Texture Segmentation An Ovrvw o Markov Random Fld and Applcaton to Txtur Sgmntaton Song-Wook Joo Octobr 003. What s MRF? MRF s an xtnson o Markov Procss MP (D squnc o r.v. s unlatral (causal: p(x t x,

More information

:2;$-$(01*%<*=,-./-*=0;"%/;"-*

:2;$-$(01*%<*=,-./-*=0;%/;-* !"#$%'()%"*#%*+,-./-*+01.2(.*3+456789*!"#$%"'()'*+,-."/0.%+1'23"45'46'7.89:89'/' ;8-,"$4351415,8:+#9' Dr. Ptr T. Gallaghr Astrphyscs Rsarch Grup Trnty Cllg Dubln :2;$-$(01*%

More information

Radial Cataphoresis in Hg-Ar Fluorescent Lamp Discharges at High Power Density

Radial Cataphoresis in Hg-Ar Fluorescent Lamp Discharges at High Power Density [NWP.19] Radal Cataphorss n Hg-Ar Fluorscnt Lamp schargs at Hgh Powr nsty Y. Aura, G. A. Bonvallt, J. E. Lawlr Unv. of Wsconsn-Madson, Physcs pt. ABSTRACT Radal cataphorss s a procss n whch th lowr onzaton

More information

Spectral stochastic finite element analysis of structures with random field parameters under bounded-but-uncertain forces

Spectral stochastic finite element analysis of structures with random field parameters under bounded-but-uncertain forces Southrn Cross Unvrsty Publcatons@SCU 23rd Australasan Confrnc on th Mchancs of Structurs and Matrals 24 Spctral stochastc fnt lmnt analyss of structurs wth random fld paramtrs undr boundd-but-uncrtan forcs

More information

Lecture 3: Phasor notation, Transfer Functions. Context

Lecture 3: Phasor notation, Transfer Functions. Context EECS 5 Fall 4, ctur 3 ctur 3: Phasor notaton, Transfr Functons EECS 5 Fall 3, ctur 3 Contxt In th last lctur, w dscussd: how to convrt a lnar crcut nto a st of dffrntal quatons, How to convrt th st of

More information

VLSI Implementation and Performance Evaluation of Low Pass Cascade & Linear Phase FIR Filter

VLSI Implementation and Performance Evaluation of Low Pass Cascade & Linear Phase FIR Filter Intrnatonal Journal of Engnrng and Tchncal Rsarch IJETR ISS: 3-869, Volum-3, Issu-6, Jun 5 VLSI Implmntaton and Prformanc Evaluaton of Low Pass Cascad & Lnar Phas Fltr Jaya Gupta, Arpan Shah, Ramsh Bhart

More information

Decision-making with Distance-based Operators in Fuzzy Logic Control

Decision-making with Distance-based Operators in Fuzzy Logic Control Dcson-makng wth Dstanc-basd Oprators n Fuzzy Logc Control Márta Takács Polytchncal Engnrng Collg, Subotca 24000 Subotca, Marka Orškovća 16., Yugoslava marta@vts.su.ac.yu Abstract: Th norms and conorms

More information

Higher order derivatives

Higher order derivatives Robrto s Nots on Diffrntial Calculus Chaptr 4: Basic diffrntiation ruls Sction 7 Highr ordr drivativs What you nd to know alrady: Basic diffrntiation ruls. What you can larn hr: How to rpat th procss of

More information

COMPLIANCE ANALYSIS, OPTIMISATION AND COMPARISON OF A NEW 3PUS-PU MECHANISM. B. Wei

COMPLIANCE ANALYSIS, OPTIMISATION AND COMPARISON OF A NEW 3PUS-PU MECHANISM. B. Wei Intrnatonal Journal of Automotv and Mchancal Engnrng (IJAME) ISSN: 9-869 (Prnt); ISSN: 8-66 (Onln); Volum 7, pp. 9-99, Januar-Jun Unvrst Malasa Pahang DOI: http://d.do.org/.58/jam.7...9-99 COMPLIANCE ANALYSIS,

More information

UNIT 8 TWO-WAY ANOVA WITH m OBSERVATIONS PER CELL

UNIT 8 TWO-WAY ANOVA WITH m OBSERVATIONS PER CELL UNIT 8 TWO-WAY ANOVA WITH OBSERVATIONS PER CELL Two-Way Anova wth Obsrvatons Pr Cll Structur 81 Introducton Obctvs 8 ANOVA Modl for Two-way Classfd Data wth Obsrvatons r Cll 83 Basc Assutons 84 Estaton

More information

EE 570: Location and Navigation: Theory & Practice

EE 570: Location and Navigation: Theory & Practice EE 570: Locaton and Navgaton: Thor & Practc Navgaton Snsors and INS Mchanaton Thursda 8 F 013 NMT EE 570: Locaton and Navgaton: Thor & Practc Sld 1 of 10 Navgaton Snsors and INS Mchanaton Navgaton Equatons

More information

A New Fast Acquisition Algorithm for GPS Receivers

A New Fast Acquisition Algorithm for GPS Receivers A Nw Fast Acuston Algorthm for GS cvrs Hung Sok So *, Chansk ark **, and Sang Jong L *** * pt. of Elctroncs Engnrng, Chungnam Natonal Unvrsty, ajon 35-764 Kora (l : 8-4-85-399; Fax : 8-4-83-4494 ; E-mal:

More information

The graph of y = x (or y = ) consists of two branches, As x 0, y + ; as x 0, y +. x = 0 is the

The graph of y = x (or y = ) consists of two branches, As x 0, y + ; as x 0, y +. x = 0 is the Copyright itutcom 005 Fr download & print from wwwitutcom Do not rproduc by othr mans Functions and graphs Powr functions Th graph of n y, for n Q (st of rational numbrs) y is a straight lin through th

More information

HORIZONTAL IMPEDANCE FUNCTION OF SINGLE PILE IN SOIL LAYER WITH VARIABLE PROPERTIES

HORIZONTAL IMPEDANCE FUNCTION OF SINGLE PILE IN SOIL LAYER WITH VARIABLE PROPERTIES 13 th World Confrnc on Earthquak Engnrng Vancouvr, B.C., Canada August 1-6, 4 Papr No. 485 ORIZONTAL IMPEDANCE FUNCTION OF SINGLE PILE IN SOIL LAYER WIT VARIABLE PROPERTIES Mngln Lou 1 and Wnan Wang Abstract:

More information

GPC From PeakSimple Data Acquisition

GPC From PeakSimple Data Acquisition GPC From PakSmpl Data Acquston Introducton Th follong s an outln of ho PakSmpl data acquston softar/hardar can b usd to acqur and analyz (n conjuncton th an approprat spradsht) gl prmaton chromatography

More information

Lesson 7. Chapter 8. Frequency estimation. Bengt Mandersson LTH. October Nonparametric methods: lesson 6. Parametric methods:

Lesson 7. Chapter 8. Frequency estimation. Bengt Mandersson LTH. October Nonparametric methods: lesson 6. Parametric methods: Otmal Sgnal Procssng Lsson 7 Otmal Sgnal Procssng Chatr 8, Sctrum stmaton onaramtrc mthods: lsson 6 Chatr 8. Frquncy stmaton Th rodogram Th modfd Prodogram (ndong Aragng rodogram Bartltt Wlch Th nmum aranc

More information

Einstein Equations for Tetrad Fields

Einstein Equations for Tetrad Fields Apiron, Vol 13, No, Octobr 006 6 Einstin Equations for Ttrad Filds Ali Rıza ŞAHİN, R T L Istanbul (Turky) Evry mtric tnsor can b xprssd by th innr product of ttrad filds W prov that Einstin quations for

More information

Emotion Recognition from Speech Using IG-Based Feature Compensation

Emotion Recognition from Speech Using IG-Based Feature Compensation Computatonal Lngustcs and Chns Languag Procssng Vol. 12, No. 1, March 2007, pp. 65-78 65 Th Assocaton for Computatonal Lngustcs and Chns Languag Procssng Emoton Rcognton from Spch Usng IG-Basd Fatur Compnsaton

More information

Fakultät III Univ.-Prof. Dr. Jan Franke-Viebach

Fakultät III Univ.-Prof. Dr. Jan Franke-Viebach Unv.Prof. r. J. FrankVbach WS 067: Intrnatonal Economcs ( st xam prod) Unvrstät Sgn Fakultät III Unv.Prof. r. Jan FrankVbach Exam Intrnatonal Economcs Wntr Smstr 067 ( st Exam Prod) Avalabl tm: 60 mnuts

More information

Chapter 13 GMM for Linear Factor Models in Discount Factor form. GMM on the pricing errors gives a crosssectional

Chapter 13 GMM for Linear Factor Models in Discount Factor form. GMM on the pricing errors gives a crosssectional Chaptr 13 GMM for Linar Factor Modls in Discount Factor form GMM on th pricing rrors givs a crosssctional rgrssion h cas of xcss rturns Hors rac sting for charactristic sting for pricd factors: lambdas

More information

EDGE PEDESTAL STRUCTURE AND TRANSPORT INTERPRETATION (In the absence of or in between ELMs)

EDGE PEDESTAL STRUCTURE AND TRANSPORT INTERPRETATION (In the absence of or in between ELMs) I. EDGE PEDESTAL STRUCTURE AND TRANSPORT INTERPRETATION (In th absnc of or n btwn ELMs) Abstract W. M. Stacy (Gorga Tch) and R. J. Grobnr (Gnral Atomcs) A constrant on th on prssur gradnt s mposd by momntum

More information

CHAPTER 33: PARTICLE PHYSICS

CHAPTER 33: PARTICLE PHYSICS Collg Physcs Studnt s Manual Chaptr 33 CHAPTER 33: PARTICLE PHYSICS 33. THE FOUR BASIC FORCES 4. (a) Fnd th rato of th strngths of th wak and lctromagntc forcs undr ordnary crcumstancs. (b) What dos that

More information

Elements of Statistical Thermodynamics

Elements of Statistical Thermodynamics 24 Elmnts of Statistical Thrmodynamics Statistical thrmodynamics is a branch of knowldg that has its own postulats and tchniqus. W do not attmpt to giv hr vn an introduction to th fild. In this chaptr,

More information

Fine Structure Calculation of Energy Levels, Oscillator Strengths and Lifetimes in Se-Like Ions

Fine Structure Calculation of Energy Levels, Oscillator Strengths and Lifetimes in Se-Like Ions ISSN: 77-75 ISO 9:8 Crtfd Intrnatonal Journal of Engnrng and Innovatv Tchnology (IJEIT) Volum, Issu 6, Dcmbr Fn Structur Calculaton of Enrgy Lvls, Oscllator Strngths and Lftms n S-Lk Ions M. Mossad, Abr

More information

Discrete Shells Simulation

Discrete Shells Simulation Dscrt Shlls Smulaton Xaofng M hs proct s an mplmntaton of Grnspun s dscrt shlls, th modl of whch s govrnd by nonlnar mmbran and flxural nrgs. hs nrgs masur dffrncs btwns th undformd confguraton and th

More information

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002 3.4 Forc Analysis of Linkas An undrstandin of forc analysis of linkas is rquird to: Dtrmin th raction forcs on pins, tc. as a consqunc of a spcifid motion (don t undrstimat th sinificanc of dynamic or

More information

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula 7. Intgration by Parts Each drivativ formula givs ris to a corrsponding intgral formula, as w v sn many tims. Th drivativ product rul yilds a vry usful intgration tchniqu calld intgration by parts. Starting

More information

4037 ADDITIONAL MATHEMATICS

4037 ADDITIONAL MATHEMATICS CAMBRIDGE INTERNATIONAL EXAMINATIONS GCE Ordinary Lvl MARK SCHEME for th Octobr/Novmbr 0 sris 40 ADDITIONAL MATHEMATICS 40/ Papr, maimum raw mark 80 This mark schm is publishd as an aid to tachrs and candidats,

More information

ANALYTICITY THEOREM FOR FRACTIONAL LAPLACE TRANSFORM

ANALYTICITY THEOREM FOR FRACTIONAL LAPLACE TRANSFORM Sc. Rs. hm. ommn.: (3, 0, 77-8 ISSN 77-669 ANALYTIITY THEOREM FOR FRATIONAL LAPLAE TRANSFORM P. R. DESHMUH * and A. S. GUDADHE a Prof. Ram Mgh Insttt of Tchnology & Rsarch, Badnra, AMRAVATI (M.S. INDIA

More information

JEE-2017 : Advanced Paper 2 Answers and Explanations

JEE-2017 : Advanced Paper 2 Answers and Explanations DE 9 JEE-07 : Advancd Papr Answrs and Explanatons Physcs hmstry Mathmatcs 0 A, B, 9 A 8 B, 7 B 6 B, D B 0 D 9, D 8 D 7 A, B, D A 0 A,, D 9 8 * A A, B A B, D 0 B 9 A, D 5 D A, B A,B,,D A 50 A, 6 5 A D B

More information