Diffraction by dielectric wedges: high frequency and time domain solutions

Size: px
Start display at page:

Download "Diffraction by dielectric wedges: high frequency and time domain solutions"

Transcription

1 UNIVERSITÀ DEGLI STUDI DI SALERNO Facolà d Ingegnera Dparmeno d Ingegnera dell Informazone ed Elerca e Maemaca Applcaa Doorao d Rcerca n Ingegnera dell Informazone XIII Cclo Nuova Sere TESI DI DOTTORATO Dffracon by delecrc wedges: hgh frequency and me doman soluons CANDIDATO: MARCELLO FRONGILLO COORDINATORE: PROF. MAURIZIO LONGO TUTOR: PROF. GIOVANNI RICCIO Anno Accademco 14 15

2

3 A u professor che ho nconrao nella ma va. Una pare d ognuno d loro è presene n queso lavoro. Soprauo, a mglor professor che ho avuo: mo padre e ma madre.

4

5 Conens Chaper 1 Inroducon 1 Chaper Hgh frequency echnques 5.1 Geomercal Theory of Dffracon 5. Geomercal Opcs 8..1 General expresson for he feld 8.. Reflecon from surfaces 11.3 Dffraced Feld Dffracon by perfecly conducng surfaces Slope dffracon Dffracon by fne conducvy surfaces: UTD heursc soluons.4 Physcal Opcs 1.5 Unform Asympoc Physcal Opcs approach: 3 Chaper 3 Hgh frequency dffracon by an arbrary-angled delecrc wedge Dffracon by delecrc wedges: sae of he ar Dffracon by an arbrary-angled delecrc wedge GO feld model for < ϕ < π/ Equvalen problems Dffraced feld: UAPO soluons for < ϕ < π/ Formulaon of he problem for π/ < ϕ < π-α Numercal resuls Formulaon of he problem for H-polarzed plane waves Parcular cases: rgh- and obuse-angled delecrc wedges Comparsons wh a rgh-angled delecrc wedge Comparsons wh an obuse-angled delecrc wedge 93 Chaper 4 Tme Doman soluons for dffracon by delecrc wedges 98

6 4.1 Tme Doman dffracon by an acue-angled delecrc wedge Formulaon of he problem for < ϕ < π/ Tme Doman UAPO dffracon coeffcens Numercal examples 18 Chaper 5 Conclusons and fuure works Summary Fuure sudes 119 Appendx A Acue angled wedge: equvalen surface currens 1 Bblography 137

7

8

9 Chaper 1 Inroducon The knowledge of he propagaon characerscs of he elecromagnec felds s fundamenal n he analyss and plannng of anennas and opcal sysems. The propagaon predcon models based on ray-racng are he mos employed echnques n modern rado communcaon sysems. The ray-based models allow one o calculae magnude and phase of he receved elecromagnec feld and also he delay of each ray due o he propagaon mechansms. These las are ofen complex and can be generally arbued o he phenomena of reflecon, ransmsson, dffracon and scaerng. As well known, when he dmensons of he sysems are large n erms of he elecromagnec wavelengh, dffracon conrbuons due o maeral dsconnues can be neglgble and mus be accuraely calculaed. By usng he numercal dscrezaon mehods such as he Fne-Dfference Tme-Doman (FDTD) mehod, he Fne Elemen Mehod (FEM), and he Mehod of Momens s possble o oban relable soluons for he scaerng problems. Unforunaely, such echnques have he drawback of he ncreasng wh frequency of he number of unknowns and so he compuaons become neffcen or nracable a he hgh frequences. Furhermore, he physcal undersandng of he propagaon mechansms s dffcul wh hs knd of approach. The asympoc mehods resul o be an effcen alernave for solvng hgh frequency scaerng problems. The Geomercal Theory of Dffracon (GTD) [1] has receved grea aenon n he las decades because of s smplcy and srengh. I s based on he rayorened heory of he Geomercal Opcs (GO) and descrbes he dffracon a he edge of a srucure lke a local phenomenon a he hgh frequences basng on a generalzaon of he Ferma prncple. The major GTD dsadvanage s he presence of feld sngulares a he GO shadow boundares. Therefore, a unform verson of GTD was

10 Chaper 1 Inroducon elaboraed, he Unform Theory of Dffracon (UTD) []. I uses he same GTD rays bu furnshes unform asympoc soluons for he dffraced feld ha are vald also n he ranson regons. As consequence of he local properes of hgh frequency elecromagnec felds, he dffracon phenomenon depends only on he geomerc and elecromagnec feaures of he obsacle n proxmy of he dffracon pon. Accordngly, for many canoncal srucures (half-plane, wedge, ec.) he correspondng dffracon coeffcens, assocang he dffraced feld o he ncden feld, are avalable n leraure. I s mporan o hghlgh ha he rgorous, closed form dffracon coeffcens have been derved only for parcular canoncal confguraons and maeral. In [1], [] he auhors produced he exac dffracons coeffcens for srucures wh perfecly conducng surfaces. In he nex research acves he maeral dsconnues have been modelled by mpedance or ransmssve boundary condons, and he soluons are generally expressed n erms of he Maluzhnes funcon ([6], [14], [19]-[1]). Ths las s almos complcaed and can be easly evaluaed only for parcular geomerc confguraons. Therefore, her applcaon o he new praccal engneerng problems s lmed. Ths work has he purpose o produce approxmae bu relable, closed form Unform Asympoc Physcal Opcs (UAPO) soluons for several canoncal problems of dffracon concernng wedges made of delecrc maerals. UAPO dffracon coeffcens are que accurae and easy o handle snce expressed n erms of he GO response of he srucure and he ranson funcon of he UTD. Furhermore, her knowledge allows one o evaluae he correspondng dffracon coeffcens n he me doman (TD) for he same geomeres. The hess s srucured as follows. In Chaper are presened some of he mos common hgh frequency echnques for analysng wave propagaon n presence of obsacles. In parcular GO, GTD, Physcal Opcs (PO) mehods are dscussed and he UAPO echnque s nroduced. In Chaper 3 he sae of research abou dffracon by delecrc wedges s examned and a hgh-frequency soluon for he problem of plane wave dffracon by an arbrary-angled delecrc wedge s shown. Alhough such a problem has been already ackled n leraure [6]-[3], he soluon proposed here s fully formulaed n he UTD

11 3 framework and hen, s a user-frendly worhwhle alernave for he consdered dffracon problem. I s obaned by a decomposon of he consdered scaerng problem no wo sub-problems, respecvely exernal and nernal o he wedge regon. For each of hem, s consdered a PO approxmaon of he elecrc and magnec equvalen surface currens n he radaon negral and a unform asympoc evaluaon of hs las s performed. The found UAPO soluon has been mplemened n a compuer code, and he smulaons show ha compensaes he GO feld dsconnues a he shadow boundares. Moreover, s accuracy s weghed by he good agreemen acheved wh numercal resuls obaned va he commercal ool Comsol Mulphyscs based on FEM and va a FDTD compuer code. Chaper 4 s devoed o he sudy of he dffracon of plane waves by an arbrary-angled delecrc wedge n he TD framework. Dfferenly from he Frequency Doman (FD) framework, n he TD one no closed form soluons are avalable n leraure for hs knd of problem, excep hose concernng rgh- and obuse-angled lossless wedges [58], [59]. The nverse Laplace ransform s appled o he FD- UAPO soluons for he dffracon coeffcens found n Chaper 3 by akng advanage of her UTD-lke naure. Then, he ransen dffraced feld s evaluaed va a convoluon negral. The resuls provded by smulaons for hese soluons have verfed her goodness. The TD-UAPO formulaon has he same advanages and lmaons of he FD-UAPO ones, arsng from he use of a PO approxmaon. Anyway, a he wrng me of hs hess he TD- UAPO formulaon s he only one analycal soluons o TD scaerng problems nvolvng penerable wedges. Alhough unavodably approxmae, he FD-UAPO soluons are smple, effcen and que accurae for solvng several dffracon problems nvolvng penerable obsacles. Furhermore, he correspondng TD-UAPO soluons are relable and he only avalable n closed form for solvng he same problems n he TD framework. Accordngly, hey can be surely useful from an engneerng vewpon when no rgorous and effcen soluons for he dffraced feld are avalable.

12

13 Chaper Hgh frequency echnques The nen of hs chaper s o dscuss some of he mos common hgh frequency echnques for analysng he propagaon of elecromagnec waves n presence of obsacles. As well-known, when an objec s nsered n he pah of a propagang elecromagnec wave modfes he feld dsrbuon n he surroundng space. Then, he oal feld a a gven observaon pon s gven by he superposon of he ncden feld and he scaered feld. The evaluaon of he scaered feld can be performed va dfferen analycal echnques. In parcular, Geomercal Opcs (GO), Geomercal Theory of Dffracon (GTD) and Physcal Opcs (PO) are llusraed n he followng, because hey consue he bases for he UAPO mehod shown a las..1 Geomercal Theory of Dffracon When he dmensons of a radang elemen or a scaerng objec are large n erms of wavelengh, s possble o use he hgh frequency asympoc echnques o solve n a que smple way many problems ha would oherwse be mahemacally nracable. In he pas years GTD, defned by Keller [1] and laer developed by Kouyoumjan and Pahak [], has receved grea aenon. I s an exenson of GO whch s an approxmae and well-known hgh frequency echnque based on he heory of propagaon of elecromagnec waves along rays. Accordng o GO, only drec, refleced and refraced rays exs; as a consequence he feld exhbs sharp ransons n correspondence of he shadow boundares and s null n he space regons no drecly llumnaed by he ncden feld. In he framework of GTD, a new ype of rays, namely he dffraced rays, s nroduced o overcome he GO naccuraces.

14 6 Chaper Hgh frequency echnques The GTD s based on he prncple of localy and on a generalzaon of he Ferma s prncple. The frs one saes ha, a he hgh frequences, reflecon, refracon and dffracon are local phenomena dependng only on he geomerc and elecromagnec properes of an objec n proxmy of he reflecon, refracon and dffracon pon, respecvely. Therefore, many dffracon coeffcens relang he dffraced feld o he ncden feld have been deermned for several canoncal problems. For example, he smples confguraons are he wedges of nfne dmenson, made of perfecly conducng maeral and llumnaed by unform plane waves. Is dffracon coeffcens have been deermned by rgorous soluons of Maxwell s equaons. As consequence, he basc mehod of GTD o ackle and solve a scaerng problem s o decompose he orgnal geomery n smpler geomerc confguraons, whose dffracon coeffcens are known. In fac, s possble o oban he fnal soluon as he superposon of he conrbuons relevan o each canoncal problem. The generalzaon of he Ferma s prncple allows one he ray racng and he deermnaon of he dffracon pons. As wellknown, saes ha n a homogeneous medum a ray moves along he shores pah from he source o he observaon pon. Such a prncple can be also generalzed o dffraced rays. The applcaon of GTD o solve scaerng problems s smple n prncple bu some problems may arse. A frs, he number of rays o consder grows rapdly wh he geomercal complexy of he problem. Moreover, even dffracon coeffcens relaed o smple srucures are sll unknown. Fnally, s dffcul o deermnae he lower frequency lm beyond whch he resuls obaned by GTD can be consdered accurae. Accordng o GTD, when an elecromagnec wave a hgh frequency s ncden on a runcaed surface, generaes a refleced wave, an edge dffraced wave and a surface wave propagang along he srucure. Ths las, can also be exced n he shadow regon of a curved surface. Such a phenomenon s represened n Fg..1, where a plane perpendcular o he edge a he dffracon pon D s shown. A ray mpngng on he edge a he pon D produces edge dffraced rays (ed) and surface rays (sr). Wh reference o convex geomeres, he surface ray produce rays dffraced from he surface a any pon Q

15 .1 Geomercal Theory of Dffracon 7 along s pah. As a resul, ES s he borderlne beween he rays dffraced from he edge and he rays dffraced from he surface and s angen o he surface n D. Moreover, SB bounds he zone llumnaed by ncden feld and RB ha llumnaed by refleced rays. The srucure n Fg..1 s non penerable oherwse also refraced rays and dffraced rays n he nner regon mus be consdered. When boh he surfaces are llumnaed by he ncden feld, here are no shadow boundares bu wo reflecon boundares arse. Snce he behavour of he opc rays s dfferen n he wo regons delmed by a boundary (ES, SB, RB), a ranson regon wheren he feld exhbs an abrup varaon s found near each boundary. Fgure.1 Incden, refleced and dffraced rays. Shadow boundares. I s assumed n he followng ha he feld sources and he observaon pon are far enough from boh he surface and he boundary ES n order o neglec he conrbuons due o he surface rays. Consequenly, he oal elecrc feld can be expressed as: wheren space whch s null n he shadow regon, r d E E E E (.1) E denoes he ncden elecrc feld, namely he feld free d E s he dffraced feld. r E s he refleced feld and

16 8 Chaper Hgh frequency echnques. Geomercal Opcs..1 General expresson for he feld In lne wh he GO, rays can be nerpreed as flux lnes of he power densy and he varaon of he feld nensy along a ray can be deermned by applyng he prncple of conservaon of he power flux n a ube of rays. Le s suppose ha a pon source emanaes soropcally waves and denoe wh L and L L wavefrons a he me and, respecvely (see Fg..). Fgure. Prmary and secondary wavefrons of a radaed wave. I s possble o consder he ube of rays beween he cross-seconal areas d a P and d a P hrough whch he power flux s consan. Therefore, he followng deny holds: wheren Sd Sd (.)

17 . Geomercal Opcs 9 E S (.3) s he power densy and he nrnsc mpedance of he medum. Eq. (.) can be rewren as: or E d E d (.4) E d E (.5) d I can be useful o express d and d n erms of he curvaure rad of he wavefrons. To hs end, le us consder he general case of an asgmac ube of rays shown n Fg..3. Fgure.3 Asgmac ube of rays. The surfaces d and d have prncpal rad of curvaure 1, and 1 s, s, respecvely, where s denoes he dsance along

18 1 Chaper Hgh frequency echnques he ray pah from P o P. Snce boh d and d subend he same sold angle, resuls: d d s s 1 1 (.6) or equvalenly E E 1 s s 1 (.7) I s evden n Fg..3 ha rays focus (cross) a dfferen pons. As a maer of fac, he ube of rays degeneraes no a lne and he feld nensy approaches nfny when s 1 or s. The locus of pons where hs occurs s called causc. I mus be sressed ha felds have always a fne value and, as a consequence, neher GO or GTD can be used o predc he feld srengh a causc pons. As he phase varaon along a ray s gven by ks, he elecrc feld along a ray can be wren as follows: E E A s e (.8) () jks wheren E denoes he elecrc feld a a reference pon P ( s ), k s he propagaon consan and As () s he spaal aenuaon facor, also known as spreadng facor. Ths las reduces o he followng expressons n he specal case of plane wave ( 1, ), cylndrcal wave ( 1, ) and sphercal wave ( 1, ): 1 plane-wave As ( ) / s cylndrcal-wave / s sphercal-wave (.9)

19 . Geomercal Opcs 11 I s worh o noe ha when a causc s crossed, s changes he sgn and hs corresponds o a phase varaon of π/. Eq. (.8) allows evaluang he feld a a gven pon P n erms of s value a a reference pon P. I has been derved by usng he prncple of conservaon of energy n a ube of rays. Alhough s a vald hgh frequency approxmaon for lgh waves, could be no accurae for elecromagnec waves of lower frequences... Reflecon from surfaces GO laws such as Snell s law can be used o deermne he feld refleced from surfaces. Le us consder a unform plane wave mpngng on a planar surface S assumed perfecly conducng. Ths case s llusraed n Fg..4, where s ˆ and s ˆr denoe he un vecors n he ncdence and reflecon drecons, respecvely. The oal feld a he reflecon pon Q on S mus mee he boundary condon: ˆn beng he un vecor normal o S. r r n ˆ E E (.1) Fgure.4 Reflecon from a planar surface. I s opporune o express he ncden and refleced felds E and n erms of her componens parallel and perpendcular o he ncdence plane (ha s he plane defned by s ˆ and ˆn ): r E

20 1 Chaper Hgh frequency echnques E Ee ˆ Eeˆ (.11) r r r r r E Ee ˆ Ee ˆ (.1) Equaon (.1) can be rewren n a compac form: r E E R r E E (.13) 1 where R s he dyadc reflecon coeffcen for a perfecly 1 conducng surface. Wh reference o Fg..4, ê s a un vecor normal o he ncdence plane, and eˆ, eˆr are un vecors parallel o he ncdence plane defned by: eˆ eˆ sˆ (.14) r r eˆ eˆ sˆ (.15) In agreemen wh he prncple of localy, he feld refleced from a surface does no change when hs las s no more planar, beng always possble o consder he local plane angen a Q r. Obvously, he spreadng facor changes o accoun for he varaon of he refleced wavefron. Therefore, s possble o express he refleced feld as follows: wheren r jks r E () s R E ( Q ) A() s e (.16) r r 1 As () (.17) r r 1 s s

21 . Geomercal Opcs 13 and s s he dsance beween he reflecon pon Q r and he r r observaon pon P. In addon, 1 and are he prncpal curvaure rad of he refleced wavefron n Q. They are relaed o he prncpal curvaure rad of he ncden wavefron 1. and o hose relevan o he surface S ( a 1 and a ). I can be shown ha 1 r and r are gven by []: r r 1 f r 1 1 f1 (.18) (.19) wheren he parameers f 1 and f n he parcular bu more neresng case of sphercal wave ncdence smplfy o: 1 1 sn sn f1, cos a1 a 1 sn sn 4 cos a a a a 1 1 (.) In eq. (.), s he ncdence angle, 1 and are he angles formed by s ˆ and he drecons assocaed o he prncpal rad of curvaure of S. In he far feld approxmaon, eq. (.16) reduces o: jks r r r e E () s R E ( Qr ) 1 (.1) s I s mmedae o verfy ha n he case of plane wave ncdence resuls:

22 14 Chaper Hgh frequency echnques r r 1 aa 1 (.) Accordngly, GO fals predcng he refleced feld when one or boh he curvaure rad of S are nfne..3 Dffraced Feld In agreemen wh eq. (.8), he dffraced feld can be expressed as: ' E () s E ( O') e d d jks s ' s (.3) d E ( O ') beng he dffraced feld a a reference pon O '. I s convenen o choose he dffracon pon Q d as reference pon. Le us consder he case of dffracon from an edge whch s, obvously, a causc for he dffraced rays. Accordngly, he lm d lm E ( O' ) ' ' (.4) exss and s no zero. Moreover, he dffraced feld mus be proporonal o he ncden feld n Q d and hus s possble o wre: d lm E ( O' ) ' DE ( Qd ) ' (.5) wheren D s usually referred o as dyadc dffracon coeffcen. As a consequence, he feld dffraced from he pon Q d on he edge s gven by: E () s DE ( Qd ) e s d jks (.6)

23 .3 Dffraced Feld 15 beng he dsance beween he causc a he edge and he second causc for he dffraced rays. In [] has been demonsraed ha: 1 1 nˆe sˆ sˆ (.7) e a sn ' where e s he curvaure radus of he ncden wavefron n Q d, aken n he plane conanng he ncden ray and he un vecor ê angen o he edge n Q d (see Fg..5). Moreover, a s he curvaure radus of he edge and n ˆe s he un vecor normal o he edge dreced away from he cener of curvaure. When he dsance s posve, he dffraced ray encouners no causc along s pah. Indeed, such a dsance s negave only when he second causc s locaed on he same sde of he observaon pon. I s opporune o remark ha he dffraced feld expressed by eq. (.6) s undefned a he causc pons and crossng such pons mples a phase varaon of π/. Fgure.5 Dffracon from a curved edge. In he case of verex dffracon, dffraced rays are orgnaed from a pon causc. Accordngly, s possble o wre: d lm E ( O' ) ' DE ( Qv ) (.8) '

24 16 Chaper Hgh frequency echnques and he feld dffraced from he verex can be evaluaed as: jks d e E () s DE ( Qv ) (.9) s.3.1 Dffracon by perfecly conducng surfaces The canoncal problem of plane wave dffracon by a perfecly conducng edge wh planar surfaces has been reaed by many auhors bu he mos popular resuls have been presened by Keller [1]. The relaed dffracon coeffcens are non-unform, namely no vald a he GO shadow boundares. Ths lmaon has been overcome hanks o he soluon provded by Kouyoumjan [], [3]. Le us consder a radang source a a pon S n presence of a perfecly conducng wedge and observe he feld a a pon P. In accordance wh he Ferma s prncple, he dffracon pons can be deermned by mnmzng he dsance SQd P. Ths leads o he law of dffracon: sˆ eˆsˆ eˆ (.3) where ê s he un vecor dreced along he edge, s ˆ and ŝ are he un vecors n he ncdence and dffracon drecons, respecvely. As consequence of eq. (.3), when a ray mpnges on he edge a oblque ncdence formng an angle ' wh respec o ê, he dffraced rays le on a concal surface whose sem-aperure angle s sll '. The ncden ray drecon s usually descrbed by a couple of angles ', ' and he poson of he dffraced ray on he cone s fxed by he angle (Fg..6). As suggesed by Kouyoumjan and Pahak, s convenen o defne he ncdence plane (ha conanng he ncden ray and he un vecor ê ) and he dffracon plane (ha conanng he dffraced ray and he un vecor ê ). A reference sysem n he ncdence plane s nroduced and fxed by he un vecors ˆ', ˆ ', s ˆ. In a smlar way, a reference sysem fxed by ˆ, ˆ, ŝ s nroduced for he

25 .3 Dffraced Feld 17 dffraced ray. The un vecors ˆ' and ˆ are perpendcular o he ncdence and dffracon planes, respecvely. The un vecors ˆ ' and ˆ are parallel o he ncdence and dffracon planes and defned by: ˆ ' sˆ ˆ ' (.31) ŝ ˆ (.3) Fgure.6 Reflecon and dffracon by a wedge. When adopng hese ray-fxed reference sysems for descrbng he ncden and dffraced felds, resuls: d d E Ds E ' jks e d D d E h E ' s( s) (.33) wheren, sh D are he scalar dffracon coeffcens for he Drchle (sof) and Neumann (hard) boundary condons, respecvely.

26 18 Chaper Hgh frequency echnques The unform dffracon coeffcens for a perfecly conducng wedge [] are repored n he followng for reader s convenence: where: j/4 e Dsh, (, ', ') n ksn ' ' co FkLa ' n ' co FkLa ' n ' co FkLa ' n ' co FkLa ' (.34) n jx j F( x) j xe e d (.35) x s he UTD ranson funcon, a and a are defned by: a nn x cos (.36) N beng he neares negers sasfyng he followng relaon: n N x (.37) The dsance parameer L n (.34) can be deermned by mposng he connuy of he oal feld a he shadow boundares [4], hus obanng:

27 .3 Dffraced Feld 19 where 1, e s 1 sn 1s s s ' L e (.38) are he prncpal rad of curvaure of he ncden wavefron and e s he curvaure radus of he ncden wavefron aken n he ncdence plane. Dffracon coeffcens for he half-plane can be easly found by seng n = n eq. (.34). I can be easly verfed ha n correspondence of a shadow boundary one of he coangen funcons becomes sngular whereas s produc wh he correspondng ranson funcon s fne. Grazng ncdence mus be consdered separaely. In hs case, Ds and he expresson of D h n eq. (.34) mus be mulpled by he facor 1/. When he wedge surfaces are curved (see Fg..1), s possble o consder n accordance wh he prncple of localy a wedge wh planar surfaces angen o he faces of he curved wedge. Moreover, he dffracon coeffcens (.34) can be sll appled bu he parameer L mus be properly modfed o accoun for he curvaure of he refleced wavefron []..3. Slope dffracon In addon o he edge dffracon conrbuon (.33), s ofen necessary o consder a hgher order erm whch s proporonal o he normal dervave of he ncden feld a Q d. Ths conrbuon may be no neglgble when he amplude of he ncden feld a Q d s small. Such ype of dffracon s usually referred o as slope dffracon. By smlar argumens, s possble o consder hgh order slope dffracon erms. When only he frs slope dffracon erm s consdered, he oal dffraced feld reads as: 1 D U Q Ud UQd A() s e jk ' n sh, d jks (.39)

28 Chaper Hgh frequency echnques where U d denoes a componen of he dffraced feld, Dsh, / ' are he slope dffracon coeffcens and U / n s he normal dervave of he ncden feld n Q d. Compac expressons for he slope dffracon coeffcens relevan o a perfecly conducng wedge can be found n [3]..3.3 Dffracon by fne conducvy surfaces: UTD heursc soluons In he pas years, he UTD has been successfully appled o solve a large varey of elecromagnec wave neracon problems wh perfecly conducng surfaces, such as he analyss of radaon characerscs of smple and complex anenna sysems [5]. However, s exenson o fne conducvy srucures s n general a non-rval problem and sll subjec of connung research. Ths neres s jusfed snce n many applcaon felds, such as rado propagaon, Radar Cross Secon (RCS) predcon, analyss of novel anenna sysems, dffracon from no perfecly conducng objecs s very frequen and he use of dffracon coeffcens (.34) could lead o naccurae resuls. As a consequence, s essenal o have relable soluons for he feld dffraced by no perfecly conducng srucures n order o make accurae feld predcons. I s opporune o noe ha rgorous and exac dffracon coeffcens have been developed for dfferen surface mpedance srucures wh he Maluzhnes heory [6]. These are raher cumbersome and, due o her complexy, ofen no easy o use n roune applcaons such as propagaon predcon ools. Thus, he dffculy arsng n usng such soluons forces smplfcaons o be made and o use approxmae echnques. Heursc soluons have been proposed n he leraure wh reference o varous scaerng problems. They are no based on a rgorous soluon of Maxwell s equaons bu on a suable modfcaon o he dffracon coeffcens (.34) wh he am o compensae he geomercal opcs feld dsconnues a he shadow boundares. Moreover, hey are easy o handle and compuaonally smple. The problem of plane wave scaerng by a hn lossless delecrc slab has been reaed n boh he normal and oblque ncdence case [7]. The derved dffracon coeffcens are an

29 .4 Physcal Opcs 1 exenson of dffracon coeffcens for he perfecly conducng halfplane. Followng he same dea n [7], wo-dmensonal dffracon coeffcens have been proposed by Luebbers for he fne conducvy wedge [8] and are currenly used n many ray racng ools. Based on hs soluon, UTD slope dffracon coeffcens have been laer developed by he same auhor [9]. As dscussed by Luebbers, he accurae use of hese dffracon coeffcens s resrced o applcaons nvolvng wedges wh large neror angles, and o he observaon n proxmy of shadow boundares. Alernave heursc UTD dffracon coeffcens for no perfecly conducng wedges have been proposed n he recen years o mprove he Luebbers soluon accuracy whn a wder angular regon [1]-[1]..4 Physcal Opcs The feld scaered n he far zone by a perfecly conducng objec llumnaed by an elecromagnec wave s gven by: s S ˆˆ S, E jk I RR J G r r' ds (.4) where J s denoes he curren dsrbuon nduced on he surface S of he objec, k and are he wavenumber and he mpedance of free jk space, r r ' Grr (, ') e 4 r r' s he Green funcon and r, r ' denoe he poson vecors a he observaon and negraon pons respecvely, ˆR s he un vecor from he radang elemen a r ' o he observaon pon and I s he (3 3) deny marx. Physcal Opcs approxmaes he curren dsrbuon J s by usng geomercal opcs and provdes accurae resuls only when he objec s far enough from he sources so ha he ncden feld can be descrbed n erms of wavefrons and rays. Accordng o hs approxmaon, he surface curren s null n he shadow regon. Whereas, a any pon P on he objec surface n he llumnaed regon can be deermned by assumng ha he ncden elecrc feld s

30 Chaper Hgh frequency echnques refleced n he same way as would be from he nfne plane angen o he surface a P [13]. When he surface S s perfecly conducng, he surface curren dsrbuon can be evaluaed accordng o: ˆ ˆ r ˆ J nh n H H n H (.41) S where ˆn s he un vecor normal o he surface a he ncdence pon and H s he ncden magnec feld. As a maer of fac, when a plane wave mpnges on a perfecly conducng surface, resuls for boh polarzaons (see Fgs..7(a) and.7(b)): r nˆh nˆ H (.4) (a) (b) Fgure.7 Reflecon from a perfecly conducng surface. Parallel polarzaon (a). Perpendcular polarzaon (b). PO can be also employed o evaluae he feld scaered from fne conducvy objecs. In such a case, boh elecrc J s and magnec J ms surface currens mus be aken no accoun n he radaon negral and so eq. (.4) becomes: ˆˆ ˆ s ms, s E jk I RR J J R G r r' ds S (.43) The accuracy of he resuls aanable wh he PO mehod depends on he degree of approxmaon for he surface currens and on he observaon drecon. As a maer of fac, when he conrbuons due

31 .5 Unform Asympoc Physcal Opcs approach 3 o he pars of he objec no drecly llumnaed by he ncden feld are no neglgble, some lmaons n he feld predcon are expeced..5 Unform Asympoc Physcal Opcs approach In he recen years, Unform Asympoc Physcal Opcs soluons have been developed o solve varous dffracon problems [14]-[18]. The basc dea o oban such soluons s he use of a PO approxmaon of he equvalen surface currens nduced by an ncden feld on a srucure. As well-known, he scaered feld s expressed by means of he radaon negral (eqns. (.4) and (.43)) and ncludes boh he geomercal opcs and he dffracon conrbuons. The applcaon of he seepes descen mehod and a unform asympoc evaluaon of he radaon negral allow dervng he dffracon coeffcens n closed form. These las are expressed n erms of he UTD ranson funcon and perfecly compensae he GO feld dsconnues. I s opporune o pon ou ha hese soluons are nevably approxmae snce he surface currens n he radaon negral are no exac bu based on a PO approxmaon. In spe of hs hey are que accurae when compared wh more rgorous soluons and are smple, easy o handle and o mplemen. For hese reason, hey are poenally useful for praccal applcaons when no exac analycal soluons are avalable or when hese las canno be evaluaed n an effcen way. Alhough heursc and effcen soluons have been proposed n he leraure (see Subsec..3.3) o solve dffracon problems havng no exac soluons, hey are no based on a rgorous analycal procedure and, herefore, should be used wh consderable aenon. The key pons showng how o consruc a UAPO soluon are revewed n he followng wh reference o he problem of plane wave dffracon by a penerable half-plane surrounded by free-space [18]. Le us consder he oblque ncdence of an arbrarly polarzed plane wave over a penerable half-plane surrounded by free-space (see Fg..8). The ncdence drecon s fxed by he angles ' and '. In parcular, ' s a measure of he ncdence drecon skewness wh

32 4 Chaper Hgh frequency echnques respec o he edge ( ' / corresponds o he normal ncdence). The observaon drecon s specfed by he angles and. I s convenen o nroduce he ray-fxed coordnae sysems sˆ, ˆ ', ˆ ' and s ˆ, ˆ, ˆ for he source and he observaon pons, respecvely (see Fg..8) and o assume he observaon pon P on he Keller s dffracon cone,.e. '. b^ s^ ^ f z y x b' ^ s ' ^ ^ f' b' Fgure.8 Dffracon by he half-plane edge. A any observaon pon P, he oal elecrc feld E s gven by he s superposon of he ncden feld E and he scaered feld E. By applyng he equvalence heorem, he surface currens nduced by an ncden plane wave on he surface S of he half-plane can be s nerpreed as sources of E and herefore, n he far-feld approxmaon, he scaered feld can be expressed by means of he followng radaon negral (.44) s ˆˆ PO PO E jk ˆ ( I RR) J s J ms RG( r, r ') ds S

33 .5 Unform Asympoc Physcal Opcs approach 5 n whch he observaon and source pons are denoed by r xxˆ y yˆzzˆ zzˆ and r' x' xˆ z' zˆ ' z' zˆ. By usng a PO approxmaon for he equvalen elecrc and magnec surface currens on S and expressng he felds n erms of her parallel and perpendcular componens, resuls: PO Js * jk ( x'sn 'cos ' z 'cos ') Jse Jms PO * jk ( 'sn 'cos ' z 'cos ') Jmse (.45) (.46) and he hree-dmensonal Green funcon can be wren as: jk rr ' e Grr (, ') 4 rr' jk ' zz ' e 4 ' ' z z (.47) To evaluae he edge dffracon confned o he Keller cone for whch ', s possble o approxmae ˆR by he un vecor ŝ n he dffracon drecon [19] Rˆ sˆsn 'cos xˆsn 'sn yˆ cos ' zˆ (.48) Accordngly, resuls: e s jk ˆˆ * * E ( I ss) Js Jmssˆ 4 jk ' zz ' jkx'sn 'cos ' z 'cos ' e ' zz' * * ( I ss ˆˆ) J s J mssˆ I s dz' dx' (.49) The expressons of he PO surface currens n erms of he ncden elecrc feld are here obaned by assumng such currens as

34 6 Chaper Hgh frequency echnques equvalen sources orgnaed by he dsconnues of he angenal GO feld componens across he layer,.e., ˆ ˆ S S r jk x'sn 'cos ' z 'cos ' s PO r J y H H y H H H yˆ H H H e (.5) PO ˆ r ms ˆ S S r jk x'sn 'cos ' z 'cos ' E E E yˆ e J E E y E E E y (.51) As well-known, s convenen o work n he sandard plane of ncdence and o consder he GO feld componens parallel ( ) and perpendcular ( ) o. Therefore, resuls: r cos ˆ r J ˆ s E E E e E E E ms r cos ˆ r ˆ (.5) J EE E e E E E (.53) wheren E and E are he ncden elecrc feld componens (a he orgn) parallel and perpendcular o he ordnary plane of ncdence, s he ncdence angle n such a plane, eˆ sˆ' yˆ sˆ' yˆ s he un vecor normal o he ordnary plane of ncdence, s ˆ' beng he un vecor of he ncdence drecon, and ˆ yˆ eˆ. The refleced and ransmed feld componens can be expressed n erms of he ncden feld componens by means of he reflecon marx R and ransmsson marx T, and her expressons can be found n [18]. As well-known, n he hgh-frequency approxmaon, he PO negral (.44) exended o S can be reduced asympocally o a sum of ray feld conrbuons from (solaed) neror saonary phase

35 .5 Unform Asympoc Physcal Opcs approach 7 pons on S and an edge dffraced feld conrbuon. To hs end, s necessary o perform he evaluaon of he followng negral: jk I s 4 jk ' zz ' jkx'sn 'cos ' jkz 'cos ' e e e dz' dx' ' zz' (.54) As repored n [19], by makng he subsuon z' z ' snh and usng one of he negral represenaons of he zeroh order () Hankel funcon of he second knd H resuls: e jkz ' cos ' e jk ' zz ' ' zz' dz' jk z cos ' () j e H k 'sn ' (.55) The nvolved Hankel funcon can be now wren wh he useful negral represenaon () 1 H k 'sn ' e d C jk 'sn 'cos (.56) where C s he negraon pah n he complex -plane (see Fg..9). The angle s beween he llumnaed face and he vecor ', and he sgn (+) apples f y ( y ). If, accordng o

36 8 Chaper Hgh frequency echnques he geomery shown n Fg..1, 'sn sn and 'cos cos x ', so obanng 1 jk sn ' cos 'sn ' () jkx'sn ' cos C H k e e d (.57) Im[a] C p Re[a] Fgure.9 Inegraon pah n he complex α-plane. and hen: C C I s k e 4 jkz cos ' jksn ' cos jkx'sn ' coscos ' e e d dx' 4 k e jkz cos ' jksn 'cos jkx'sn ' coscos ' e e dx' d (.58)

37 .5 Unform Asympoc Physcal Opcs approach 9 By applyng he Sommerfeld-Maluzhnes nverson formula, resuls: jk x 1 e dx' (.59) jk sn ' cos cos ' 'sn ' cos cos ' y P r r - r' f q x z r' Fgure.1 Inegraon pah n he complex -plane. so ha Is e 1 e sn ' j cos cos ' C jkz cos ' jksn ' cos d (.6) where he sgn (+) apples n he range ( ). Such an negral can be evaluaed by usng he Seepes Descen Mehod [19]. To hs end, he negraon pah C s closed wh he Seepes Descen Pah (SDP) passng hrough he pernen saddle pon s as shown n Fg..11. Accordng o he Cauchy resdue heorem, he conrbuon relaed o he negraon along C (dsored for he presence of sngulares n he negrand) s equvalen o he sum of he negral along he SDP and he resdue conrbuons Re s p assocaed wh all hose poles ha are nsde he closed pah C+SDP,.e.,

38 3 Chaper Hgh frequency echnques 1 f I s g e d j C 1 f Re s p g e d j SDP Re s p I (.61) n whch 1 f I g e d j SDP f e s ffs g e d j (.6) SDP g jk e z cos ' 1 A sn ' cos cos ' cos cos ' (.63) k f jsn ' cos. Noe ha s ypcally large, p ' and s ( s ) f ( ). Moreover, by pung ' j" and mposng ha Im f Im f s and Re f Re f s, he consdered SDP s descrbed by: and 1 1 s s gd ' sgn( '')cos ( '') (.64) cosh '' where gd " s he Gudermann funcon. By usng now he change of varable f f s, eq. (.6) can be wren as

39 .5 Unform Asympoc Physcal Opcs approach 31 I G e d (.65) wheren 1 f s d G g e (.66) j d Im[a] SDP C ap Re[a] as p When Fgure.11 SDP n he complex -plane. p s approachng s, he funcon G canno be expanded n a Taylor seres. To overcome hs drawback s convenen o regularse he negrand n (.65) usng he Mulplcave Mehod. I requres nroducng he regularsed funcon wh G p G p (.67) sn ' 1 cos ' p f s f p j

40 3 Chaper Hgh frequency echnques jsn 'cos ' j (.68) and s a measure of he dsance beween p and s. Accordngly, Snce p e I Gp d p (.69) G s analyc near, can be expanded n a Taylor seres. By reanng only he frs erm (.e., he snce 1, resuls: n whch G I F j 1 -order erm) p p (.7) p Gp 1 f s d G g s e j d p jsn ' 1 Ae j 4 e j cos cos ' sn ' (.71) and j F j j e e d (.7) s he UTD ranson funcon []. By subsung (.68) and (.71) n (.7), he explc form of he asympoc evaluaon of I s:

41 .5 Unform Asympoc Physcal Opcs approach 33 I j 4 jksn ' z cos ' e e k sn ' sn ' cos cos ' ' F k sn 'cos j 4 jks e e 1 k s sn ' cos cos ' ' F kssn 'cos (.73) where he denes s sn ' and z scos ' are used on he dffracon cone. The above analyc resul conrbues o he UAPO dffraced feld o be added o he GO feld and s referred o as a unform asympoc soluon because I s well-behaved when p s. In he GTD framework, he marx formulaon for he scaered feld can be wren as s ' E E s E M I s E E ' s (.74) As a consequence he UAPO edge dffracon conrbuon d E s: d jk ' ' s E E E d e E M I D d E E ' E ' s (.75) where he UAPO soluon for he dffracon marx D D D ' ' D D ' '

42 34 Chaper Hgh frequency echnques M11 M /4 1 j 1 e M 1 M k sn ' (cos cos ' ) ' F ks sn ' cos (.76) where he sgn ( ) apples f ( ). The explc expresson of he coeffcens M j (, j = 1, ) can be found n [18]. Accordngly, he UAPO soluons have he same ease of handlng of hose derved n he UTD framework and, n addon, hey have he nheren advanage of provdng he dffraced feld from he knowledge of he GO response of he srucure. In oher words, s suffcen o make explc he reflecon and ransmsson coeffcens relaed o he consdered srucure for obanng he UAPO dffracon coeffcens.

43

44 Chaper 3 Hgh frequency dffracon by an arbrary-angled delecrc wedge The problem of dffracon by a delecrc wedge has grea relevance for praccal applcaons. The accurae pah loss predcon n rado wave propagaon envronmens requres a correc characerzaon of dffracon conrbuons arsng from he presence of edges and corners. However, he research acvy focused manly on mpenerable srucures a hgh frequency and rgorous soluons have been repored n he leraure for perfecly conducng wedges or mpedance wedges. Sarng from hese soluons, here have been some aemps o exend her valdy o he penerable srucures. Ths mehods resul o be complcaed and manageable jus for smple confguraons, so has been necessary o buld up new approaches for hs knd of problems. The lack of an exac and, a he same me, compuaonally effcen soluon whch can be employed n ray racng ools forces smplfcaons o be made and o apply approxmae soluons. The am of hs chaper s o provde a UAPO soluon for he feld dffraced by he edge of a lossless arbrary-angled delecrc wedge n he case of normal ncdence. The soluon s derved by means of he decomposon of he consdered scaerng problem no wo subproblems relevan o exernal and nernal regons of he wedge. For each of hem, proper equvalen currens, whch can be nerpreed as sources for he scaered felds, are deermned by akng no accoun he penerable naure of he srucure. Then, unform asympoc evaluaons of he radaon negrals allow dervng closed form expressons for he dffraced feld. As demonsraed by numercal examples, he here developed UAPO soluon compensaes he GO feld dsconnues n he exeror and neror regons. Furhermore,

45 3.1 Dffracon by delecrc wedges: sae of he ar 37 s accuracy s confrmed by he good agreemen aaned wh resuls provded by numercal mehods. 3.1 Dffracon by delecrc wedges: sae of he ar The problem of plane wave dffracon by wedges has receved grea aenon due o he mporance of s soluons n rado propagaon plannng, analyss and desgn of radang srucures and wavegude heory. The frs sudes have concerned wedges wh perfecly conducng surfaces [], [6] or mpedances faces [14], [], [1]. Heursc soluons have also been proposed n [8], [9] for solvng dffracon problems n he case of delecrc srucures. Alhough effcen, hey are no based on a rgorous analycal procedure and, herefore, should be used wh consderable aenon. Rawlns [] presened an approxmae soluon for he feld produced when an elecromagnec plane wave s dffraced by an arbrary-angled delecrc wedge, whose refracve ndex s near uny. I was obaned from an applcaon of he Konorovch-Lebedev ransform and a formal Neumann-ype expanson. The resuls were n agreemen wh hose derved by he same auhor wh reference o a rgh-angled wedge [3]. The dffraced feld of an E-polarzed plane wave by a delecrc wedge was formulaed n erms of negral equaons by Bernsen n [4]. These were ransformed no Fredholm negral equaons and solved by erave mehods for lmed values of he delecrc consan. Joo e al. [5] proposed an asympoc soluon for he feld dffraced by a rgh-angled delecrc wedge based on a correcon of PO approxmaon o he edge dffracon. The correcon n he far-feld zone was calculaed by solvng a dual seres equaon agreeable o smple numercal evaluaon. The exenson of he approach o arbrary-angled wedges was addressed n wo companon papers [6] and [7]. Burge e al. [8], sarng from a PO verson of he GTD, provded an edge coeffcen for he exernal and nernal dffracon by an arbrary-angled delecrc wedges. Rouvere e al. n [9] mproved he Luebbers heursc soluon [8] n he UTD framework [] by addng wo new erms ha compensae he dsconnuy creaed by he feld ransmed hrough he srucure. Anoher heursc soluon for he dffracon coeffcen

46 38 Chaper 3 Hgh frequency dffracon by an arbrary-angled delecrc wedge of a penerable wedge has been proposed n [3] by Bernard e al.. I s obaned by a modfcaon of he exac UTD dffracon coeffcen relaed o he meallc wedge. The FDTD mehod was appled n [31] and [3], provdng numercal resuls for he dffracon coeffcen of rgh-angled wedges made of perfec elecrcal conducor (PEC), lossless delecrc and lossy delecrc maeral. Anoher neresng aemp o solve he problem of he dffracon by a rgh-angled delecrc wedge was proposed by Radlow [33]. He used muldmensonal Wener-Hopf equaons o model he problem, bu he facorzaon of hese equaons needs funcon-heorec echnques employng wo complex varables ha are cumbersome o handle. Baes [34] nroduced erave formulas o oban a se of basc wave-funcons useful o represen he sources suaed on he plane of symmery of an arbrary-angled delecrc wedge. Seo and Ra [35] proposed a remarkable mehodology o evaluae he scaerng from a lossy delecrc wedge: refleced and refraced GO felds are obaned by nhomogeneous plane wave racng n he lossy medum, whle PO approxmaon n he edge dffraced felds are accuraely correced by addng he mulpole lne sources a he edge of he wedge o sasfy he exncon heorem. The unknown expanson coeffcens used a he p of he wedge o make dsappearng he exncon negral are evaluaed numercally. Salem e al. [36] exended Rawlns approach [] o a more general form of excaon (lne source) of a delecrc wedge. They expressed he Konorovch- Lebedev ransform funcons as a Neumann seres and represened he scaered feld as a Bessel funcon seres, so exendng soluons o he case of real valued wavenumbers and arbrarly posoned source and observer. Danele and Lombard [37] analyzed he problem of dffracon of an ncden plane wave by an soropc penerable wedge n he general case of skew ncdence. They used generalzed Wener- Hopf equaons, and he soluon s obaned usng analycal and numercal-analycal approaches ha reduce he Wener-Hopf facorzaon o Fredholm negral equaons of second knd. Vaslev e al. [38], [39] developed a numercal approach based on he mehod of negral equaons o solve he problem. They represened he unknown surface currens on he delecrc wedge as a sum of unform and non-unform componens (usng he Ufmsev s mehod), makng

47 3. Dffracon by an arbrary-angled delecrc wedge 39 possble he dervaon of he numercal soluon. Budaev [4] developed a hybrd mehodology o solve he problem of dffracon by a delecrc wedge n an exac sense, whch combnes analycal and numercal echnques. UAPO soluons were proposed n explc closed form n [41] and [4] for evaluang he feld dffraced by rgh- and obuse-angled lossless delecrc wedges n he nner and ouer regons. The UAPObased approach has been also appled o acue-angled wedges n [43] consderng a specfc range of ncdence drecons for boh cases of E- and H-polarzed plane wave. Then, n [44] he analyss has been exended o all possble cases of ncdence drecon provdng generalzed UAPO soluons whch are vald for he dffracon by a wedge wh arbrary aperure. 3. Dffracon by an arbrary-angled delecrc wedge In hs secon he UAPO soluons for he dffraced feld orgnaed by an arbrary-angled lossless delecrc wedge are deermned for an ncden plane wave. The geomery used as reference n he followng s relevan o a wedge havng an acue nernal apex angle. Ths choce s jusfed by he complexy of he propagaon mechansms (mulple nernal and exernal rays) whch nclude hose concernng rgh- and obuse-angled wedges [41], [4], as wll be demonsraed a he end of hs chaper. Le us consder he problem of he dffracon of an E-polarzed plane-wave by he edge of an acue-angled delecrc wedge (see Fg. 3.1). Is nernal apex angle s and s maeral s lossless, nonmagnec ( r 1), wh relave permvy r and propagaon consan kd kd r. The wedge surfaces are denoed by S and S n. A Caresan reference sysem x, yz, s nroduced wh he y-axs normal o he face S and he z-axs dreced along he edge. The ncdence drecon s assumed perpendcular o he edge (see Fg. 3.1) and defned by he angle ' ( ' corresponds o he grazng ncdence wh respec o S ).

48 4 Chaper 3 Hgh frequency dffracon by an arbrary-angled delecrc wedge Fgure 3.1 Geomery of he dffracon problem. The observaon pon s denoed by P(, ). The srucure dvdes he space n wo regons: he exeror regon and he neror regon GO feld model for < ϕ < π/ Closed form expressons for he GO response when ' are derved n hs secon. Only he face S s drecly llumnaed by he ncden feld. A ransmed ray exss n he nner regon and, n absence of oal nernal reflecon on S n, s ransmed hrough hs las. The refleced ray ravels oward S and gves rse o followng reflecon and ransmsson mechansms nvolvng S and S n. When he value of he nernal angle of ncdence makes o occur he nverson he rays begn o go far from he apex. Accordngly, wo seres of nernal and exernal rays exs: he frs conans rays (black arrows n Fg. 3.) ravellng owards he apex and he second comprses hose (red arrows n Fg. 3.) gong far from he apex, as llusraed n Fg. 3.. The exernal rays exs unl he nernal oal reflecon doesn occur. The number of pre- and pos-nverson rays s deermned by he values of ', and r. In order o beer undersand he ray propagaon nsde and ousde he wedge, he reader can refer o Fg If ' s he exernal ncdence angle, he wave

49 3. Dffracon by an arbrary-angled delecrc wedge 41 peneraes no he wedge hrough S wh he ransmsson angle L sn sn and propagaes oward S n. 1 r ' S Fgure 3. Inernal reflecon and exernal ransmsson. The ray mpnges on S n wh he ncdence angle 1 and, for he r Snell s laws, s refleced ( 1 1 ) and ransmed hrough S n wh he angle 1, as shown n Fg Accordng o he consdered values of ' and r, he oal nernal reflecon can occur and no ransmed wave hrough S n exss. In such a case, an evanescen wave propagang along S n and aenuang n he drecon normal o arses. The condon ensurng he absence of oal reflecon a P 1 s: c 1 arcsn 1 r (3.1) c where s he correspondng crcal angle. By denong wh ˆ s cos ', sn ', he un vecor n he ncdence drecon, resuls: jksˆ r jk cos ˆ ' E Ee z Ee zˆ (3.) S n

50 4 Chaper 3 Hgh frequency dffracon by an arbrary-angled delecrc wedge wheren E s he amplude of he ncden feld a he orgn. The feld E r refleced by S can be deermned as: r jk dr E R E P e (3.3) where R s he reflecon coeffcen a P and dr s he dsance ' S Fgure 3.3 Rays ransmssons/reflecons hrough he wedge from P o P along he refleced ray: dr cos ' x cos' (3.4) n whch x r s he coordnae of P. As a resul, eq. (3.3) reads as: As regards he feld y x sˆ1 1 E 1 sˆr E r 1 1 r r jk cos' E REe z (3.5) ˆ E L ransmed hrough S, resuls: P L jk d r E x r P sˆr1 E L P 1 1 ˆ L s S n E T E P e (3.6) ˆ s E

51 3. Dffracon by an arbrary-angled delecrc wedge 43 T beng he ransmsson coeffcen a P and d s he dsance from P o P along he ransmed ray: As a consequence, resuls: L sn L sn L d x r (3.7) cos sn sn L L jk sn d L jkx r ' r jkd ˆ E T E e e z (3.8) TEe zˆ The feld E r 1 ransmed hrough S and hen refleced a P 1 ( 1, 1 ) on S n wh 1 s gven by: r1 1 L jkddr E R E P1 e 1 (3.9) wheren R 1 s he reflecon coeffcen a P 1 and dr 1 s he dsance from P 1 o P along he ray refleced from S n : Now, s smple o verfy ha 1 1 sn 1 sn 1 dr (3.1) 1 L (3.11) and so eq. (3.9) reads as r1 jk sn d 1 1 ˆ E RT E e z (3.1) If 1 c, he feld E 1 ransmed hrough S n can be evaluaed as:

52 44 Chaper 3 Hgh frequency dffracon by an arbrary-angled delecrc wedge 1 1 jk 1 d E TE P e 1 (3.13) wheren T 1 s he ransmsson coeffcen a P 1 and d 1 s he dsance from P 1 o P along he ransmed ray: d1 1sn1 sn 1 (3.14) 1 where 1 sn r sn1 s he ransmsson angle a P 1. Therefore, eq. (3.13) reads as: The feld r 1 jksn 1 1 ˆ E TT E e z (3.15) E refleced a, P on S wh s gven by: r r jk dr E RE P e 1 d (3.16) wheren R s he reflecon coeffcen a P and dr s he dsance from P o P along he refleced ray: sn sn dr (3.17) where. I can be easly verfed ha: L r jk sn d 1 ˆ E R RT E e z (3.18) If c, he feld E ransmed hrough S can be evaluaed as:

53 3. Dffracon by an arbrary-angled delecrc wedge 45 r jk d E T E P e (3.19) 1 wheren T s he ransmsson coeffcen a P and d s he dsance from P o P along he ransmed ray: 1 Snce sn r sn sn sn d (3.) eq. (3.19) reads as: jksn 1 ˆ E TRTEe z (3.1) The ray ransmed hrough S ravels oward he apex undergong N In L reflecons/ransmssons ( In denoes he neger par of he argumen) before movng away from. Accordngly, r n n L n for n 1,,..., N and he oal nernal reflecon c doesn happen f n. Wh reference o he frs se of rays propagang owards he apex, he feld conrbuons for n 3,..., N can be evaluaed by reerang he prevous approach. The expressons of he correspondng nernal ( E ) and exernal ( E ex L ) GO felds are: L n N n L r E n L E E n1 sn N n sn jkd L jkd n ET e Rp e n1 p1 neven

54 46 Chaper 3 Hgh frequency dffracon by an arbrary-angled delecrc wedge and N n jk sn d n Rp e zˆ (3.) n1 p1 nodd N 1 ex r n L n1 jkcos( ') zw ˆ (, ') jkcos( ') E E E E Ee REe zw ˆ (, ') N1 n1 jksn p n ET R p Te n n1 p1 nodd 3 W pn, Uc n zˆ N 1 n1 jksnqn ET Rp Te n W, qnuc nzˆ n1 p1 neven (3.3) where N N N 1 f N odd f N even (3.4) he funcon W () s defned as:

55 3. Dffracon by an arbrary-angled delecrc wedge 47 W, 1 U () s he un sep funcon: 1 f 1 (3.5) elsewhere U c c n n 1 elsewhere (3.6) and he parameers p and q have been defned as: 1 f nn p (3.7) 1 elsewhere 1 f nn1 q (3.8) 1 elsewhere The expresson (3.4) allows one o keep n accoun boh he possble cases of nverson so ha (3.) and (3.3) are generally vald, he un sep funcon W () s useful o keep no accoun he shadow boundares and he parameers p and q have been nroduced o have smlar expressons for he wo ses of neracon pons, as wll be shown n he followng. A crucal neracon pon s he N 1 -h one whch corresponds o he nernal ncdence angle N 1 N1 L. The nverson occurs a hs pon and he ray ravels far from he apex, as llusraed n Fg The use of he parameers p (eq. 3.7) and q (eq. 3.8) allows one o consder he frs conrbuon of hs se as arsng from S. Then, s opporune o descrbe hs seres of neracons sarng from hs ray. The angle R N r 1 N 1 plays now he role of L for evaluang he feld conrbuons. The number of neracon pons of he ray propagang far from he apex s

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

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

II. Light is a Ray (Geometrical Optics)

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

More information

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

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

( ) () 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

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

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

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

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

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

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

2.1 Constitutive Theory

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

More information

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

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

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

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

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

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours NATONAL UNVERSTY OF SNGAPORE PC5 ADVANCED STATSTCAL MECHANCS (Semeser : AY 1-13) Tme Allowed: Hours NSTRUCTONS TO CANDDATES 1. Ths examnaon paper conans 5 quesons and comprses 4 prned pages.. Answer all

More information

Chapter 2 Linear dynamic analysis of a structural system

Chapter 2 Linear dynamic analysis of a structural system Chaper Lnear dynamc analyss of a srucural sysem. Dynamc equlbrum he dynamc equlbrum analyss of a srucure s he mos general case ha can be suded as akes no accoun all he forces acng on. When he exernal loads

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

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

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

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

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

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

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

Part II CONTINUOUS TIME STOCHASTIC PROCESSES

Part II CONTINUOUS TIME STOCHASTIC PROCESSES Par II CONTINUOUS TIME STOCHASTIC PROCESSES 4 Chaper 4 For an advanced analyss of he properes of he Wener process, see: Revus D and Yor M: Connuous marngales and Brownan Moon Karazas I and Shreve S E:

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

[ ] 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

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

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

More information

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 Finite Element Method for the Analysis of Non-Linear and Dynamic Systems

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Swss Federal Insue of Page 1 The Fne Elemen Mehod for he Analyss of Non-Lnear and Dynamc Sysems Prof. Dr. Mchael Havbro Faber Dr. Nebojsa Mojslovc Swss Federal Insue of ETH Zurch, Swzerland Mehod of Fne

More information

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria IOSR Journal of Mahemacs (IOSR-JM e-issn: 78-578, p-issn: 9-765X. Volume 0, Issue 4 Ver. IV (Jul-Aug. 04, PP 40-44 Mulple SolonSoluons for a (+-dmensonalhroa-sasuma shallow waer wave equaon UsngPanlevé-Bӓclund

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

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair TECHNI Inernaonal Journal of Compung Scence Communcaon Technologes VOL.5 NO. July 22 (ISSN 974-3375 erformance nalyss for a Nework havng Sby edundan Un wh ang n epar Jendra Sngh 2 abns orwal 2 Deparmen

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

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

. 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

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

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

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

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

More information

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

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

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

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

More information

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

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

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

. 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

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

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

A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION

A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION S19 A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION by Xaojun YANG a,b, Yugu YANG a*, Carlo CATTANI c, and Mngzheng ZHU b a Sae Key Laboraory for Geomechancs and Deep Underground Engneerng, Chna Unversy

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

3. OVERVIEW OF NUMERICAL METHODS

3. OVERVIEW OF NUMERICAL METHODS 3 OVERVIEW OF NUMERICAL METHODS 3 Inroducory remarks Ths chaper summarzes hose numercal echnques whose knowledge s ndspensable for he undersandng of he dfferen dscree elemen mehods: he Newon-Raphson-mehod,

More information

Anisotropy and oblique total transmission at a planar negative-index interface

Anisotropy and oblique total transmission at a planar negative-index interface Ansoropy and oblque oal ransmsson a a planar negave-ndex nerface Le Zhou, C.T. Chan and P. Sheng Deparmen of Physcs, The Hong Kong Unversy of Scence and Technology Clear Waer Bay, Kowloon, Hong Kong, Chna

More information

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

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

More information

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

Robust and Accurate Cancer Classification with Gene Expression Profiling

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

More information

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

Sampling Procedure of the Sum of two Binary Markov Process Realizations

Sampling Procedure of the Sum of two Binary Markov Process Realizations Samplng Procedure of he Sum of wo Bnary Markov Process Realzaons YURY GORITSKIY Dep. of Mahemacal Modelng of Moscow Power Insue (Techncal Unversy), Moscow, RUSSIA, E-mal: gorsky@yandex.ru VLADIMIR KAZAKOV

More information

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

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

More information

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

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

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

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

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

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

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

More information

H = d d q 1 d d q N d d p 1 d d p N exp

H = d d q 1 d d q N d d p 1 d d p N exp 8333: Sacal Mechanc I roblem Se # 7 Soluon Fall 3 Canoncal Enemble Non-harmonc Ga: The Hamlonan for a ga of N non neracng parcle n a d dmenonal box ha he form H A p a The paron funcon gven by ZN T d d

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

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

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

, t 1. Transitions - this one was easy, but in general the hardest part is choosing the which variables are state and control variables

, t 1. Transitions - this one was easy, but in general the hardest part is choosing the which variables are state and control variables Opmal Conrol Why Use I - verss calcls of varaons, opmal conrol More generaly More convenen wh consrans (e.g., can p consrans on he dervaves More nsghs no problem (a leas more apparen han hrogh calcls of

More information

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

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

More information

APPROXIMATE ANALYTIC SOLUTIONS OF A NONLINEAR ELASTIC WAVE EQUATIONS WITH THE ANHARMONIC CORRECTION

APPROXIMATE ANALYTIC SOLUTIONS OF A NONLINEAR ELASTIC WAVE EQUATIONS WITH THE ANHARMONIC CORRECTION THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY Seres A OF THE ROMANIAN ACADEMY Volume 6 Number /5 pp 8 86 APPROXIMATE ANALYTIC SOLUTIONS OF A NONLINEAR ELASTIC WAVE EQUATIONS WITH THE ANHARMONIC

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

Lecture 6: Learning for Control (Generalised Linear Regression)

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

More information

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

Appendix H: Rarefaction and extrapolation of Hill numbers for incidence data

Appendix H: Rarefaction and extrapolation of Hill numbers for incidence data Anne Chao Ncholas J Goell C seh lzabeh L ander K Ma Rober K Colwell and Aaron M llson 03 Rarefacon and erapolaon wh ll numbers: a framewor for samplng and esmaon n speces dversy sudes cology Monographs

More information

Polymerization Technology Laboratory Course

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

More information

ESTIMATIONS OF RESIDUAL LIFETIME OF ALTERNATING PROCESS. COMMON APPROACH TO ESTIMATIONS OF RESIDUAL LIFETIME

ESTIMATIONS OF RESIDUAL LIFETIME OF ALTERNATING PROCESS. COMMON APPROACH TO ESTIMATIONS OF RESIDUAL LIFETIME Srucural relably. The heory and pracce Chumakov I.A., Chepurko V.A., Anonov A.V. ESTIMATIONS OF RESIDUAL LIFETIME OF ALTERNATING PROCESS. COMMON APPROACH TO ESTIMATIONS OF RESIDUAL LIFETIME The paper descrbes

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

Stochastic Maxwell Equations in Photonic Crystal Modeling and Simulations

Stochastic Maxwell Equations in Photonic Crystal Modeling and Simulations Sochasc Maxwell Equaons n Phoonc Crsal Modelng and Smulaons Hao-Mn Zhou School of Mah Georga Insue of Technolog Jon work wh: Al Adb ECE Majd Bade ECE Shu-Nee Chow Mah IPAM UCLA Aprl 14-18 2008 Parall suppored

More information

Anisotropic Behaviors and Its Application on Sheet Metal Stamping Processes

Anisotropic Behaviors and Its Application on Sheet Metal Stamping Processes Ansoropc Behavors and Is Applcaon on Shee Meal Sampng Processes Welong Hu ETA-Engneerng Technology Assocaes, Inc. 33 E. Maple oad, Sue 00 Troy, MI 48083 USA 48-79-300 whu@ea.com Jeanne He ETA-Engneerng

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

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

Transient Response in Electric Circuits

Transient Response in Electric Circuits Transen esponse n Elecrc rcus The elemen equaon for he branch of he fgure when he source s gven by a generc funcon of me, s v () r d r ds = r Mrs d d r (')d' () V The crcu s descrbed by he opology equaons

More information

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

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

More information

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

MEEN Handout 4a ELEMENTS OF ANALYTICAL MECHANICS

MEEN Handout 4a ELEMENTS OF ANALYTICAL MECHANICS MEEN 67 - Handou 4a ELEMENTS OF ANALYTICAL MECHANICS Newon's laws (Euler's fundamenal prncples of moon) are formulaed for a sngle parcle and easly exended o sysems of parcles and rgd bodes. In descrbng

More information

The Elastic Wave Equation. The elastic wave equation

The Elastic Wave Equation. The elastic wave equation The Elasc Wave Eqaon Elasc waves n nfne homogeneos soropc meda Nmercal smlaons for smple sorces Plane wave propagaon n nfne meda Freqency, wavenmber, wavelengh Condons a maeral dsconnes nell s Law Reflecon

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

Dual Approximate Dynamic Programming for Large Scale Hydro Valleys

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

More information

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

ELASTIC MODULUS ESTIMATION OF CHOPPED CARBON FIBER TAPE REINFORCED THERMOPLASTICS USING THE MONTE CARLO SIMULATION

ELASTIC MODULUS ESTIMATION OF CHOPPED CARBON FIBER TAPE REINFORCED THERMOPLASTICS USING THE MONTE CARLO SIMULATION THE 19 TH INTERNATIONAL ONFERENE ON OMPOSITE MATERIALS ELASTI MODULUS ESTIMATION OF HOPPED ARBON FIBER TAPE REINFORED THERMOPLASTIS USING THE MONTE ARLO SIMULATION Y. Sao 1*, J. Takahash 1, T. Masuo 1,

More information

Econ107 Applied Econometrics Topic 5: Specification: Choosing Independent Variables (Studenmund, Chapter 6)

Econ107 Applied Econometrics Topic 5: Specification: Choosing Independent Variables (Studenmund, Chapter 6) Econ7 Appled Economercs Topc 5: Specfcaon: Choosng Independen Varables (Sudenmund, Chaper 6 Specfcaon errors ha we wll deal wh: wrong ndependen varable; wrong funconal form. Ths lecure deals wh wrong ndependen

More information

ON THE WEAK LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS

ON THE WEAK LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS ON THE WEA LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS FENGBO HANG Absrac. We denfy all he weak sequenal lms of smooh maps n W (M N). In parcular, hs mples a necessary su cen opologcal

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

THE method of moment (MOM) is widely used to extract. Fast Green Function Evaluation for Method of Moment. arxiv: v1 [cs.

THE method of moment (MOM) is widely used to extract. Fast Green Function Evaluation for Method of Moment. arxiv: v1 [cs. 1 Fas Green Funcon Evaluaon for Mehod of Momen Shunchuan Yang, Member, IEEE, Dongln Su, Member, IEEE arxv:1901.04162v1 [cs.ce] 14 Jan 2019 Absrac In hs leer, an approach o accelerae he marx fllng n mehod

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