Lecture 4 Excitation & Acquisition Lecture Notes by Assaf Tal. To Measure A Signal, Spins Must Be Excited

Size: px
Start display at page:

Download "Lecture 4 Excitation & Acquisition Lecture Notes by Assaf Tal. To Measure A Signal, Spins Must Be Excited"

Transcription

1 Lecure 4 Eciaion & Acquisiion Lecure Noes b Assaf Tal To easure A Signal, Spins us e Ecied The Saic Nuclear agneic Field Is Too Weak To e Reliabl Deeced We ve previousl calculaed he bulk magneic momen of mm 3 of waer in a 3 T magneic field and found i o be abou ~ -8 J/T Imagine his voel is inside he human and mus be deeced in a coil wrapped around he sa, cm awa from i The magneic field of he voel is dipolar r 3 rr ˆ ˆ 4 r 3 If we consider a sphere of radius r around, he field will be maimal a he poles: aimal field where is value will be (puing r= m, = -8 J/T, =4-7 N/A ): 3 r 3 ~ T= pt r This is an incredibl small magneic field The eremel weak magneic fields produced b neural currens are - orders of magniude larger, and even hose require eremel specialied hardware based on super-cooled magneic field deecors known as SQUIDs These are used in magneoencephalograph (EG) machines, bu he are no sensiive enough o deec he sor of signals we re ineresed in I is possible o consruc sensiive enough deecors, bu even if we pu aside he engineering complei and cos of hese in addiion o he RI magne iself, we are sill faced wih furher problems: How can we separae he in field creaed b he nuclear spins from he much larger sources of magneic fields creaed b oher phenomena in he, eg neural currens? And how can we, b deecing he magneic fields ouside he, deduce he disribuion of spins wihin he? Tha is, can we image (r) b measuring (r) ouside he? These sor of problems are known as inverse problems and are ofen incredibl difficul o solve properl (he same issue plagues EG) Is here a beer wa o deec he RI signal? The answer is es, and i is linked o he pession of spins around he main field Pessing Spins Induce easurable Currens In The Receiver Coils The R in RI sands for resonance As we ve remarked previousl, an spin will pess abou a consan magneic field In paricular, if we place he spins in a saic, srong magneic field sa, he 3 Tesla field of a pical RI scanner i will pess For a hdrogen nucleus, his pession frequenc would be kh T 7 H mt Using a coil and he law of iproci we can measure he ime-depdenden flu induced b he spin This dnamic measuremen will creae a significanl larger signal han a saic one We are, however, faced wih a parado: a hermal equilibrium, he bulk magneic momen is parallel o, and hence he pession is degenerae remains saic (even hough he microscopic momens, m, will pess) For us o observe rue pession, mus make some non-ero angle wih :

2 he origin) The epression for he magneic field creaed b a loop of curren a is cener is well known from basic magneism: A hermal equilibrium, and are colinear, and no pession is observed R I nˆ R I Onl b creaing some nonero angle beween and ha is, onl b eciing - can we observe pession Creaing such an angle is called eciaion Puing ha aside for a momen, le s calculae he volage induced in a coil pu around he pessing spin via Farada s law and he law of iproci Here we ake a wo circular coils of radius R in he and planes, wih a poin magneic momen placed a he origin and performing some aion as a funcion of ime: where R is he ring s radius, I he curren, and ˆn a uni vecor normal o he plane of he ring We ake uni curren (I=) and so: d d ˆ R R cos R Number Time For =4-7 Vs/(Am), =7 H for pons a 3 Tesla, R = 5 cm (head coil), and = -8 J/T (previousl calculaed magneiaion of ml of waer), we ge 3 V, around he righ order of magniude for he volages deeced in magneic resonance The componens of he pessing spin have a sinusoidal ime dependence: cos sin LH aion abou Remember ha he law of iproci ells us ha he volage induced in he coil can be calculaed via v () dm r, where is he field creaed b a uni curren he loop a he posiion of he magneic momen (a This is a small bu deecable volage level wih oda s elecronics, and his is he basis of signal epion in modern RI The smaller he radius of he coils, R, he beer: alwas build coils ha are as small as possible! Furhermore, he signal is proporional o =, and increases wih (alhough an eac analsis of he SNR will awai a laer chaper) The Spins Can e Ecied Wih A Resonan RF Field A Small RF Field Can Have a Large Effec if i is Resonan The discussion of he previous secion has shown ha, in order o induce non-ero volage, we mus il he magneiaion vecor awa from equilibrium and have i pess Indeed, he basic R eperimen can be described as follows:

3 Thermal Equilibrium: A hermal equilibrium, he spins are aligned along and do no pess Eciaion: The spins are somehow ecied, ha is, iled o some angle wih respec o This usuall happens quickl and relaaion can be negleced Pession & Deecion: Once iled, he pess and give off a ime dependen magneic field The magneic field induces a volage in a nearb RF coil via Farada s law We can also furher manipulae he spins wih magneic fields during his period o bring ou paricular conras pes We usuall have a ime ~ T before decoherence "eas up" he observable pessing magneiaion Thermaliaion: Relaaion processes kick in The ransverse magneiaion decas wih a ime consan T while he longiudinal magneiaion builds up back up due o T relaaion If we wai for a ime 5T, he magneiaion will be back a is hermal equilibrium value Each such block (ecie-acquire-wai) is called a scan I is in fac no mandaor o wai for a ime 5T for he spins o reurn o hermal equilibrium; we'll see laer on ha waiing a shorer amoun of ime has boh benefis (shorer scan imes) and disadvanages (less signal per scan) For now, however, we'll assume ha is he case, so is equal o and poins along he -ais before he beginning of each scan We ve alread remarked ha RF << How can we hope o non-negligibl ecie he spins wih such a weak RF field? The answer is ha we use a resonan field ha oscillaes a he Larmor frequenc Namel, we are going o solve he loch equaions seing G=, and RF ("on resonance irradiaion") This means we will need o solve he loch equaions wih a ime dependen magneic field Alhough a numerical soluion is possible, we will emplo a frame ransformaion rick which will enable us o solve his problem wihou an Transforming o a Frame Which Roaes A The Same Frequenc As The RF Field akes i Appear Saic: The Roaing Frame In he laboraor frame, his amouns o solving he loch equaions wih a complicaed imedependen magneic field The loch equaions are easier o solve in a frame which aes around he -ais wih a frequenc given b RF We ackle his as follows: consider a saic (laboraor) frame wih ime independen, fied uni vecors ˆ, ˆ, ˆ, and a aing frame wih uni vecors ˆ', ˆ', ˆ' ' ' ' wih ˆ sin cos ˆ RF RF RF If he aing frame is aing wih an angular veloci abou an ais given b he uni vecor ˆn, hen each of he aes of he aing frame pess abou he vecor = ˆn This means Schemaic for a simple single eperimen ("scan") (no drawn o scale) Ecie he spins This usuall akes << T, T Le spins pess while acquiring a signal unil i decas due o T Wai for magneiaion o "build back up" due o T before eciing again Time

4 each obes a pession equaion idenical (formall) o he loch equaion: dˆ ' ˆ ' ω dˆ ' ˆ ' ω dˆ ' ˆ ' ω The magneiaion vecor can be epressed in eiher frame: ˆ ˆ ˆ ˆ' ˆ' ˆ' For eample, if () is,,, ˆ sin cos ˆ RF RF RF as given above, while he aing frame aes a he angular frequenc = RF abou he -ais ( = ˆ RF ) hen he componens of in he wo frames are: cos RF, sin RF,, Noe he componens change, bu he vecor is frame-independen since i is a geomerical quani Differeniaing () wih respec o ime, we obain d d d,, d, ˆ' ˆ' ˆ' dˆ' dˆ' dˆ',,, d d,, d, ˆ' ˆ' ˆ' ˆ' ˆ' ˆ' ω d ω,,, Equaing, we obain: d d ω This is pisel he loch equaion bu wih an ecive field ω The above equaion is rue for an aing frame However, in RI, when we speak of he aing frame, we will be referring o a frame which aes a a consan angular veloci RF abou he -ais according o he lef hand rule For he aing frame, ω ˆ ˆ RF When epressed in he aing frame, he componens of he ecive field ω are:,,, (in he aing frame: RF ) If we selec RF we are on resonance: he RF irradiaes he spins a he same frequenc as heir naural frequenc, In his case: On resonance: RF In he aing frame: RF An Analog From echanics Imagine he earh going around he sun in a circle: On he oher hand, he loch equaion sas

5 Gravi This can be undersood b an observer in space he following wa: he Earh wans o go forward bu gravi pulls i inward, curving is pah ino a circle In ec, he Earh is coninuousl falling ino he sun, bu escaping doom hanks o is angenial veloci All his is all a consequence of Newon s second law, F=ma Ne, imagine how hings would look o an observer sanding on he sun and aing wih i Neglecing for he ime being he weaher on he surface, he Earh would appear saionar o such an observer: Gravi If ha observer would r o use Newon s law F=ma o undersand his world he would fail: according o F=F gravi =ma, earh should be falling owards he sun, bu i isn! The ruh is ha when ou ransform o a aing frame ou need o add a ficiious force Tha is, ou need o presuppose a force which doesn arise ou of an phsical source, called he cenripeal force, o eplain how i is possible for he earh o remain saionar: Cenripeal So, in mechanics when ou r o undersand hings in a aing frame ou need o do wo hings: Undersand how hings in he real frame would look in he aing frame (eg, he Earh would remain sill) Add ficiious forces (eg, he cenripeal force) A similar hing happens when ou go o a aing frame in magneic resonance, aing wih he same angular veloci as he RF field: Firs, he RF field appears saionar in he aing field which maches is aion frequenc (ie because RF ) Now we need o add he cor ficiious "force" - field, o be pise - given b fic To see, imaging a saic spin in he lab frame, wih no magneic field Now ransform o a frame aing wih an angular veloci abou he -ais In his frame, he spin would appear o ae wih an angular veloci fic, as if here was a ficiious field fic presen along he -ais The ulk agneiaion Pesses Around The Effecive Field In The Roaing Frame We've seen he magneiaion vecor obes he loch equaions in he aing frame, onl swapping he field for an ecive field, ω and epressing ha field in he aing frame basis (ie as i would appear o an observer aing wih he frame) This means pesses abou in he aing frame Saring from hermal equilibrium a ime =, poins along (aken o coincide wih he -ais) in boh he laboraor and he aing frames, which are also assumed o coincide for =: Gravi

6 ' ' Now we urn on he resonan RF field in he laboraor frame: ˆ sin cos ˆ RF RF RF This field aes in he -plane in he lab frame, and appears saionar in he aing frame Furhermore, if we assume our irradiaion is on resonance, RF=, he ecive field in he aing frame has no -componen: RF Lab Frame ˆ Lab Frame RF The magneiaion pesses abou he ais in he aing frame We can hus creae an angle we'd like beween i and he -ais, depending on how long we le i pess and how srong is Le's assume we have RF on for jus enough ime for he magneiaion o il o he plane - ha is, creae a 9 angle beween and Deducing he moion of in he lab frame is now merel a maer of ransforming back o he lab frame, which simpl aes a an angular veloci - relaive o he aing frame Tha is, in he lab frame performs a spiral as i descends and aes: ' ' The magneic field in he laboraor frame has a large -componen and a small, aing componen (no shown o scale) In he aing frame, assuming RF is on resonance (RF==) he ecive field is saic ' Ro Frame A ime = (hermal equilibrium), poins along he -ais in boh frames Ro Frame ' ' Seing The Radiofrequenc (RF) Pulse s (Area)=(Duraion)(Ampliude) Ses The Flip Angle We see he spins will perform a aion abou he -ais in he aing frame a a frequenc = Noe his is no he same as RF, (one is he ampliude of RF, he second is is oscillaing frequenc) Afer a ime, will have creaed an angle : Noe ha ampliude of RF duraion of RF This relaion is rue onl on resonance, when RF =, where has no -componen To ip he magneiaion ono he ais, we wai a ime 9 such ha: or Lab Frame Shown here is he rajecor of he magneiaion in he lab (lef) and aing (righ) frames The wo frames are conneced b a simple aion RF(), In he original laboraor (unaing) frame he spins eecue addiional moions, bu he imporan hing o realie is ha a spin which is in ' ' ' Ro Frame ' Ro Frame '

7 he plane in he aing frame, mus also be in he -plane in he laboraor frame (alhough where in he plane is a differen sor!) Number Time We ve remarked ha,ma ~ T for an RI scanner For pons, one would need 9 ~5 mso ecie he 3 spins ono he -plane For C, Signal 9 Is Acquired ~ ms Via Farada s Law Relaaion Can e Negleced During Eciaion Since os Pulses Are Shorer Than T, T Our calculaions in he previous secion have shown ha eciaion mosl happens on he imescale of milliseconds in RI, which is much shorer han T, T Hence, o an ecellen approimaion, relaaion ecs can be negleced for mos pulses and mos issue pes in he We will make some remarks abou he ecs of relaaion laer on bu, in general, will neglec i unless specificall saed oherwise The Phase of he Pulse Deermines The Phase of he Ecied agneiaion We have so far modeled RF in he lab frame as: cos c RF, lab sin c Since we have full conrol over he and componen we have no problem modulaing boh () and adding a ime-dependen phase () o he RF field: cos RF, sin Le's keep () and () fied Then he consan phase ()= is called he phase of he pulse, and is equal o he angle he RF field makes wih he -ais deermines where he RF pulse will poin in he ransverse plane The phase of he magneiaion is defined as he angle made b he ransverse componen of he magneiaion vecor (ie is projecion on he plane) wih he -ais ecause he magneiaion ges ipped a righ angles o he RF field following he lef hand rule, he relaion beween he pulse's and magneiaion's phase is given b: m The sandard noaion for a consan RF pulse hen assumes he form, where is is flip angle and is (consan) phase The following convenions are also used: RF : 9 : 8 : 7 : Some eamples are shown below (magneiaion is assumed o sar ou from, and is he blue vecor; he RF is he red vecor): RF, lab cos c sin c In he aing frame, his will look like his : R, where To prove his, use lab R c is a RH aion mari abou he -ais b an c angle c (he aing frame aes wih a lef handed aion and angular frequenc c; in i, i appears he RF field aes a he same angular frequenc bu in he opposie diion) There is a bi of algebra and rigonomer involved bu he proof is sraighforward

8 (9) pulse (45) pulse On he oher hand, he pulse's duraion, 9 magneiaion in plane (phase: 9) (9)- pulse 9 magneiaion in plane (phase: ) Eciaion Flip Angles < 9 Are Used To inimie The Duraion, Decrease Power Deposiion A The Cos Of SNR An eciaion pulse need no ip he spins b 9, and can creae an angle beween and he main field The disadvanage of his is is reduced signal: in our simple model we've seen ha he volage, d, R is proporional o he ime derivaive of (reoriening he coil would inroduce he ime derivaive of, and would no change our conclusions) The magniude of will be proporional o he flip angle Hence, he signal iself will also be proporional o sin() and decrease wih he flip angle : 45 magneiaion in plane (phase: 9) (45) pulse 45 magneiaion in plane (phase: 8), is proporional o he flip angle and decreases linearl (assuming we keep fied) Anoher advanage of shor pulses is ha he have reduced specific absorpion rae (SAR) Some of he RF energ is absorbed in he paien s issue and causes undesired heaing The amoun of SAR is proporional o he square of and he pulse s duraion: consan pulse, duraion T SAR We observe SAR reduces linearl wih he flip angle The amoun of SAR is limied b mos modern scanners' hardware based on our undersanding of he ec of SAR on biological issues odern RF coils deposi power on par wih modern cell-phones and are generall considered safe as long as guidelines are observed Off-Resonan Eciaion And The Concep of Selecive Eciaion The -Field Can Var as a Funcion of Posiion, Which Leads o Non-ero Offses In The Roaing Frame So far our approach has been o make disappear b moving o a aing frame a a frequenc, in which he ficiious field negaes compleel: lab signal, sin We will see laer on his decrease is acuall miigaed in mos sequences where pulses are applied rapidl (on he order of, or faser han T) and don afford he magneiaion enough ime o reurn o hermal equilibrium before he ne eciaion However, when varies as a funcion of posiion, r, i is impossible o make he -componen of he field disappear a ever poin: lab r r Some sources of variaion could include:

9 Imperfecions in he main magne field Suscepibili arifacs in he sample, in which he eernal field induces microscopic magneic momens which hemselves disor he main field (in all diions, bu predominanl in he diion on ) 3 Some paiens migh have meal implans which disor he magneic field again, in man diions, bu heir ec is mos pronounced along 4 Ofen we inenionall creae hese inhomogeneiies, as is he case wih gradien coils, in which we creae a linear dependence of he -field on posiion: G r For eample, when a gradien is urned on, lab G r G r, and, when he RF is urned on: G r If we have some form of spaial inhomogenei due o hardware imperfecions/suscepibili arifacs, we could wrie i as lab r, and in he aing frame is -componen will be lab r Number Time A gradien will creae a range of frequencies given b G over a spaial region of wih Across a mm piel, his will be G 4 H Suscepibili arifacs a 3 Tesla will creae spaial variaions across he head on he order of hundreds of H, mosl in regions where air-issue inerfaces eis such as he prefronal core, close o he oral cavi or ears, and so forh In he case he -componen is no compleel eroed ou, we mus anale and undersand he case for which The quani will be referred o as he offse of he spins Since we will be looking a a specific poin in space we can assume is jus a consan A Qualiaive Soluion: All Pulses Are Selecive Wih A Finie andwih Given ~ / I is fairl simple o divide our analsis ino wo ereme cases: in one, <<, and we can neglec i, obaining: We hus over he previous case in which we ecie he spins as usual, as if he were on resonance On he oher ereme,, and he RF eciaion will have no ec, resuling in no eciaion We can guess and erapolae beween hese wo eremes, saing ha here is a cuoff o he ec of when ~ In oher words, a range of offses ~ will be ecied This is known as he bandwih of he pulse: he range of offses (or frequencies) i will ecie W This can also be undersood graphicall, b ploing he pession cone of he spins abou he ecive field, saring ou from hermal equilibrium (ie along he -ais):

10 = poins in he -plane sa, along he -ais in his eample - and he spin pesses in a circle in he plane << sars iling up in he plane, causing he aion cone o sar folding = In his dividing case, he - and - componens of are equal The pession cone jus ouches he plane >> now becomes close o he -ais The pession cone becomes ver narrow: even if we wai for a long amoun of ime he spins will no sra far from he -ais A Consan Pulse's Eciaion Profile Le us eplore wha happens when our pulse is no hard and has a finie duraion We ll look a a 9 pulse, alhough our conclusions will appl o an pulse flip angle For a 9 flip angle, he duraion mus be 9 4 However, as previousl discussed, his onl ensures a 9 flip angle for spins when = As we increase, he magniude and diion of and is diion var, and he corresponding final posiion of he magneiaion - assumed o sar ou from hermal equilibrium along he -ais - varies In fac, even if is applied along, spins no a he cener do no even remain in he plane anmore The following diagram shows he ecive field and he pise rajecor raced b he magneiaion vecor during he pulse's duraion, 9, for he four cases oulined previousl: = >> A Hard Pulse Is One Wih High Peak Power And An Infinie andwih When we are ineresed in flipping all of he spins ono he -plane regardless of heir offse we mus creae a bandwih larger han he range of offses in our sample Since W, his means we need o have a ver high >>range of offses in our sample The duraion of an (on-resonance) 9-pulse, 9 4 We see ha hard pulses are shor and have a high peak power Such pulses are called hard pulses in R jargon The ideal hard pulse is one for which, duraion, such ha duraion equals he desired flip angle = << Insead of hese picorial diagrams, one can plo he componens of as a funcion of he offse, This is known as he pulse's response or pulse profile Such a response is ploed below for = kh ( 35 T), 9=5 ms ( applied along he -ais as in he above diagrams, ie has phase):

11 pulse we had =, = Le us see wha happens if we var () in a sinc-like manner: (kh) ime (ms) (kh) ime (ms) In he above, The magneiaion was assumed o have an arbirar magniude of uni, and he verical ais sreches from The dashed red lines signif he poins a which kh, which define he bandwih of he pulse I is quie clear ha he concep of bandwih has some arbirariness o i since he profile of and are no sharp The profile of acuall looks somewha wider han, which is a resul of he magneiaion vecor s consan magniude, and he relaion beween and (kh) For eample, if =9, hen 8 and 436 So, even if is almos unperurbed, migh sill appear o be quie siable, leading o is wider profile Also noe he eensive wiggles in ouside he slice, indicaing ha some eciaion occurs even for >> We will deal wih his in a momen b inroducing shaped pulses Shaping The Pulse Affecs The Pulse's Profile odern RF ransmiers have he capabili of shaping he RF pulse, RF i ; ha is, conrolling is - and - componens Such pulses are called shaped pulses For he consan The new pulse mainains he same peak and same area (and hence flip angle) as he angular pulse, bu is necessaril longer (since he negaive lobes of he sinc derac from he area) The frequenc response of his pulse can be calculaed b solving he loch equaions numericall, ielding: (kh) The response is shown using he same scaling and plo range as he angular pulse for a "fair" comparison The dashed red lines represen he same bandwih (= ) calculaed for he angular pulse The ensuing response is significanl beer-behaved, wih less wiggles and sharper ransiion lines We won' go ino he heor of shaped pulses in his course, bu we will remark wihou proof ha for ip angles up unil abou 9 he profile of resembles he Fourier ransform of RF ()

12 The Slice s Cener Can e Shifed Sinusoidall odulaing The Pulse s Phase The pulses discussed so far ecie a bandwih abou a cenral frequenc There is a ver simple wa o shif he cener of he ecied slice, b modulaing he pulse s phase wih a linear erm ahemaicall, his means: c shif cos e ff sin c To undersand wh his shifs he pulse, imagine shif being given in he aing frame, where i aes around he -ais according o he lef hand rule wih angular frequenc c performing a second aing frame ransformaion, ino a frame which aes wih c relaive o he original ( firs ) aing frame, we fi and add an addiional ficiious field: nd aing shif frame c becomes i i cos i cos i e i c c sin c c sin c The comple noaion makes he ec seem much simpler (ie onl mulipling he comple RF i c waveform b ) e Slice Selecion Ofen, We Are Ineresed In Eciing A Single Slice We now come o he firs form of spaial selecivi in RI: selecive eciaion, in which onl a par of he sample is ecied We will confine ourselves o he simple scenario of one-dimensional eciaion, meaning selecivel eciing a range of posiions along a fied ais: In his nd aing frame he -componen has a fied offse All of our previous argumens can be repeaed, bu now he cener of he profile would no occur a, bu raher a a frequenc defined b, or: c c Curren RI hardware enables one o conrol he phase and ampliude of he RF pulse as a funcion of ime, making shifing he profile s cener eas (we can generae an pracical frequenc c ) I also places almos no demands on he hardware 3 I s ineresing o noe ha, using comple noaion, our consan original field Wihou such selecive eciaion, an RF pulse would ecie he enire sample Alhough his can be and is someimes done, a slice-selecive approach also has is own meris Appling A Pulse In The Presence Of A Gradien Will Ecie A Slice In he join presence of a gradien and an RF irradiaion, he ecive field in he aing frame is: 3 This rick for shifing he pulse's profile and he relevan analsis oulined, also holds for non-consan pulses r, G r

13 We will assume for simplici our gradien is consan and urned on along he -diion, so G G ˆ, and so: G RF, RF, G The gradien creaes a linearl increasing offse along : Offse Pos G G Finall, jus as he profile of a pulse can be shifed in frequenc space as a funcion of he offse, hus he slice's cener can also be shifed b simpl modulaing he RF pulse's shape: The gradien assigns frequencies o posiions via G, and hence an pulse ha ecies a range of frequencies W will, in he presence of a gradien, ecie a range of posiions given b: W G G This is, in fac, he slice hickness The ecied slice will be perpendicular o he diion of he gradien vecor G For eample appling he same pulse wih a gradien in a differen diion - sa, a 45 o he -ais - will ecie a slice ha is iself iled b 45, since now our gradien will creae a linear correspondence beween frequenc and he +=cons planes: Acquisiion The Time Evoluion Of The agneiaion In The Roaing Frame In The Absence Of RF Irradiaion In he absence of RF irradiaion, he loch equaions consis onl of a -field, made up of is posiion-dependen offse due o gradiens, inhomogeneiies and so forh: r r We would like o solve he loch equaions: macro, macro, T macro, macro, T macro, macro, T Subsiuing he field above, we obain: r r T T T

14 The firs hing o noe is ha he ransverse () and longiudinal () componens of () are decoupled: he ransverse magneiaion pesses in he -plane and decas wih a ime consan T, while he longiudinal componen builds up owards hermal equilibrium wih a ime consan T, and is unaffeced b he -componen of he field We can solve for independen of he gradiens, field inhomogeneiies, ec he soluion having been alread inroduced in lecure 4: T / T / e e The equaions for and are mied ogeher, or coupled in mahemaical erms We will herefore rea he ransverse magneiaion ( and ) as one eni: i We now mulipl he second equaion b i and add up he firs wo equaions: d i d i r, i T The LHS is merel he ime derivaive of Similar simplificaions can be made for he RHS, and we obain: d This equaion is of he form wih i r, T d a air, T If a were consan, is soluion would be simple: a e Once a is ime dependen we mus break down he ime inerval ino small seps such ha a() is consan in each ime sep Then he soluion in each inerval is a e The full soluion is obained b concaenaing he shor-ime soluions: a a 3 a 4 3 a3 e e e e e e e aa aaa a3aaa This can be coninued b inducion, wih he sum urning ino an inegral as : a ' ' e Reurning back o our case and subsiuing he epression for a()= r,, we obain: ', ' / T,, r r r e Usuall, = will correspond o he poin in ime righ afer he eciaion pulse, so (r,) will be he ransverse magneiaion righ afer eciaion Noe ha we haven' reall said if we're in he aing or lab frame, and our derivaion - and final epression - are valid for boh Since lab r, r,, we can also deduce ha lab r e r, i, which makes sense, because he lab and aing frames differ b a consan aion abou he - ais b an angular frequenc The Acquired Signal In RI Is Proporional To The Transverse agneiaion We now sae a ver fundamenal relaion in RI: he acquired signal is proporional o he inegral of he magneiaion, weighed b he eiver profile:

15 r r s dv, where is he spaial field dependence of he eiver coil when we drive a uni curren hrough i In he remainder of his secion we derive his equaion, bu ou can skip his wihou loss of coninui Proof: Armed wih our epression for he ime course of he magneiaion, we urn o deriving a usable epression for he acquired signal (volage) in our eiver coil We've seen in Lecure 3 ha, for an eiver coil, d r, v r dv d d d dv Here and are boh measured in he laboraor frame Firs, he ime derivaive of he -componen of he magneiaion, which changes wih a ime consan T (on he order of H) is much smaller han he - and - componens, which pess wih a frequenc of H Therefore, we can neglec he -componen o an ecellen approimaion: d d v dv Since we have an epression for in he lab frame, we can differeniae i Since varies much faser han T and G()r, we can neglec boh erms and obain 4 : so d lab r, lab i r, v Re i r r, dv The ineresing signal is modulaed b a rapid i (~ H) phase erm, e, which is "hidden" inside I is beneficial o ge rid of i for boh convenience, as well as o lessen he burden on he analog-o-digial converer (which needs o deal wih slower varing signals) This is called demodulaion To do his, v is spli ino he idenical copies, wih one muliplied b cos( ) and he oher b sin( ), and each is hen passed hrough a low-pass filer (LPF) v sin( ) cos( ) LPF LPF ADC ADC Real channel Imaginar channel We define i, i and noe: so i i i Re Using and wriing cos sin e e e e i i i i i Therefore: lab d v Re dv 4 This assumpion needs o be modified when suding solids wih ver shor Ts (on he order of microseconds) However, such shor Ts are rarel observed in RI since he lead o signals well below he noise levels

16 lab lab i r, r, r, e i r, r, r, e we obain, righ before he LPF: i i Re im v i e e dv i i Re re v e e dv Since he LPF removes he fas changing componen, we obain, afer he LPFs: i Re im v i e dv sin dv i Re re v e dv cos dv We hen form he comple signal in he compuer: Since s v i v re im lab i lab e dv, and since he aing and lab frame magneiaion vecors are relaed via lab i r, r, e, we can simplif: r r s dv, We've omied he consan of proporionali since he acual measured signal's magniude will depend anwa on he elecronics, amplifiers and so on

Get: Nuclear (equilibrium) magnetization M 0. (Magnitude dictated by Boltzmann distribution)

Get: Nuclear (equilibrium) magnetization M 0. (Magnitude dictated by Boltzmann distribution) 9: Relaaion of nuclear magneiaion. How is he R signal deeced?. Wha is he quanum-mechanical equivalen of he roaing frame? 3. Wha is he roaing frame descripion good for? 4. How can he reurn of he magneiaion

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

NMR Spectroscopy: Principles and Applications. Nagarajan Murali 1D - Methods Lecture 5

NMR Spectroscopy: Principles and Applications. Nagarajan Murali 1D - Methods Lecture 5 NMR pecroscop: Principles and Applicaions Nagarajan Murali D - Mehods Lecure 5 D-NMR To full appreciae he workings of D NMR eperimens we need o a leas consider wo coupled spins. omeimes we need o go up

More information

Traveling Waves. Chapter Introduction

Traveling Waves. Chapter Introduction Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

15. Bicycle Wheel. Graph of height y (cm) above the axle against time t (s) over a 6-second interval. 15 bike wheel

15. Bicycle Wheel. Graph of height y (cm) above the axle against time t (s) over a 6-second interval. 15 bike wheel 15. Biccle Wheel The graph We moun a biccle wheel so ha i is free o roae in a verical plane. In fac, wha works easil is o pu an exension on one of he axles, and ge a suden o sand on one side and hold he

More information

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17 EES 16A Designing Informaion Devices and Sysems I Spring 019 Lecure Noes Noe 17 17.1 apaciive ouchscreen In he las noe, we saw ha a capacior consiss of wo pieces on conducive maerial separaed by a nonconducive

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation Hea (iffusion) Equaion erivaion of iffusion Equaion The fundamenal mass balance equaion is I P O L A ( 1 ) where: I inpus P producion O oupus L losses A accumulaion Assume ha no chemical is produced or

More information

Biomedical Imaging. Nuclear Magnetic Resonance. Patrícia Figueiredo IST,

Biomedical Imaging. Nuclear Magnetic Resonance. Patrícia Figueiredo IST, Biomedical Imaging Nuclear agneic Resonance Parícia Figueiredo IST, 213-214 The wide specrum of medical imaging echniques (F. Deconinck, Vrije Universi, Belgium). Overview Nuclear magneism Precession:

More information

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem) Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =

More information

1. Kinematics I: Position and Velocity

1. Kinematics I: Position and Velocity 1. Kinemaics I: Posiion and Velociy Inroducion The purpose of his eperimen is o undersand and describe moion. We describe he moion of an objec by specifying is posiion, velociy, and acceleraion. In his

More information

Relaxation. T1 Values. Longitudinal Relaxation. dm z dt. = " M z T 1. (1" e "t /T 1 ) M z. (t) = M 0

Relaxation. T1 Values. Longitudinal Relaxation. dm z dt. =  M z T 1. (1 e t /T 1 ) M z. (t) = M 0 Relaxaion Bioengineering 28A Principles of Biomedical Imaging Fall Quarer 21 MRI Lecure 2 An exciaion pulse roaes he magneiaion vecor away from is equilibrium sae (purely longiudinal). The resuling vecor

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

HW6: MRI Imaging Pulse Sequences (7 Problems for 100 pts)

HW6: MRI Imaging Pulse Sequences (7 Problems for 100 pts) HW6: MRI Imaging Pulse Sequences (7 Problems for 100 ps) GOAL The overall goal of HW6 is o beer undersand pulse sequences for MRI image reconsrucion. OBJECTIVES 1) Design a spin echo pulse sequence o image

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow 1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering

More information

Chapters 6 & 7: Trigonometric Functions of Angles and Real Numbers. Divide both Sides by 180

Chapters 6 & 7: Trigonometric Functions of Angles and Real Numbers. Divide both Sides by 180 Algebra Chapers & : Trigonomeric Funcions of Angles and Real Numbers Chapers & : Trigonomeric Funcions of Angles and Real Numbers - Angle Measures Radians: - a uni (rad o measure he size of an angle. rad

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

NMR Spectroscopy: Principles and Applications. Nagarajan Murali 2D NMR Heteronuclear 2D Lecture 7

NMR Spectroscopy: Principles and Applications. Nagarajan Murali 2D NMR Heteronuclear 2D Lecture 7 NMR pecroscop: Principles and Applicaions Nagarajan Murali D NMR Heeronuclear D Lecure 7 Heero Nuclear D-NMR Two dimensional NMR can be used o correlae NMR signals arising from differen nuclei such as

More information

From Complex Fourier Series to Fourier Transforms

From Complex Fourier Series to Fourier Transforms Topic From Complex Fourier Series o Fourier Transforms. Inroducion In he previous lecure you saw ha complex Fourier Series and is coeciens were dened by as f ( = n= C ne in! where C n = T T = T = f (e

More information

- If one knows that a magnetic field has a symmetry, one may calculate the magnitude of B by use of Ampere s law: The integral of scalar product

- If one knows that a magnetic field has a symmetry, one may calculate the magnitude of B by use of Ampere s law: The integral of scalar product 11.1 APPCATON OF AMPEE S AW N SYMMETC MAGNETC FEDS - f one knows ha a magneic field has a symmery, one may calculae he magniude of by use of Ampere s law: The inegral of scalar produc Closed _ pah * d

More information

This is an example to show you how SMath can calculate the movement of kinematic mechanisms.

This is an example to show you how SMath can calculate the movement of kinematic mechanisms. Dec :5:6 - Kinemaics model of Simple Arm.sm This file is provided for educaional purposes as guidance for he use of he sofware ool. I is no guaraeed o be free from errors or ommissions. The mehods and

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

IB Physics Kinematics Worksheet

IB Physics Kinematics Worksheet IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

Maxwell s Equations and Electromagnetic Waves

Maxwell s Equations and Electromagnetic Waves Phsics 36: Waves Lecure 3 /9/8 Maxwell s quaions and lecromagneic Waves Four Laws of lecromagneism. Gauss Law qenc all da ρdv Inegral From From he vecor ideni da dv Therefore, we ma wrie Gauss Law as ρ

More information

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3. Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies

More information

V The Fourier Transform

V The Fourier Transform V he Fourier ransform Lecure noes by Assaf al 1. Moivaion Imagine playing hree noes on he piano, recording hem (soring hem as a.wav or.mp3 file), and hen ploing he resuling waveform on he compuer: 100Hz

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A

More information

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4. PHY1 Elecriciy Topic 7 (Lecures 1 & 11) Elecric Circuis n his opic, we will cover: 1) Elecromoive Force (EMF) ) Series and parallel resisor combinaions 3) Kirchhoff s rules for circuis 4) Time dependence

More information

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan Ground Rules PC11 Fundamenals of Physics I Lecures 3 and 4 Moion in One Dimension A/Prof Tay Seng Chuan 1 Swich off your handphone and pager Swich off your lapop compuer and keep i No alking while lecure

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information

PHYSICS 149: Lecture 9

PHYSICS 149: Lecture 9 PHYSICS 149: Lecure 9 Chaper 3 3.2 Velociy and Acceleraion 3.3 Newon s Second Law of Moion 3.4 Applying Newon s Second Law 3.5 Relaive Velociy Lecure 9 Purdue Universiy, Physics 149 1 Velociy (m/s) The

More information

k B 2 Radiofrequency pulses and hardware

k B 2 Radiofrequency pulses and hardware 1 Exra MR Problems DC Medical Imaging course April, 214 he problems below are harder, more ime-consuming, and inended for hose wih a more mahemaical background. hey are enirely opional, bu hopefully will

More information

Basic Circuit Elements Professor J R Lucas November 2001

Basic Circuit Elements Professor J R Lucas November 2001 Basic Circui Elemens - J ucas An elecrical circui is an inerconnecion of circui elemens. These circui elemens can be caegorised ino wo ypes, namely acive and passive elemens. Some Definiions/explanaions

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

Electromagnetic Induction: The creation of an electric current by a changing magnetic field.

Electromagnetic Induction: The creation of an electric current by a changing magnetic field. Inducion 1. Inducion 1. Observaions 2. Flux 1. Inducion Elecromagneic Inducion: The creaion of an elecric curren by a changing magneic field. M. Faraday was he firs o really invesigae his phenomenon o

More information

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15. SMT Calculus Tes Soluions February 5,. Le f() = and le g() =. Compue f ()g (). Answer: 5 Soluion: We noe ha f () = and g () = 6. Then f ()g () =. Plugging in = we ge f ()g () = 6 = 3 5 = 5.. There is a

More information

CHAPTER 12 DIRECT CURRENT CIRCUITS

CHAPTER 12 DIRECT CURRENT CIRCUITS CHAPTER 12 DIRECT CURRENT CIUITS DIRECT CURRENT CIUITS 257 12.1 RESISTORS IN SERIES AND IN PARALLEL When wo resisors are conneced ogeher as shown in Figure 12.1 we said ha hey are conneced in series. As

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

Second-Order Differential Equations

Second-Order Differential Equations WWW Problems and Soluions 3.1 Chaper 3 Second-Order Differenial Equaions Secion 3.1 Springs: Linear and Nonlinear Models www m Problem 3. (NonlinearSprings). A bod of mass m is aached o a wall b means

More information

Chapter 3 Kinematics in Two Dimensions

Chapter 3 Kinematics in Two Dimensions Chaper 3 KINEMATICS IN TWO DIMENSIONS PREVIEW Two-dimensional moion includes objecs which are moing in wo direcions a he same ime, such as a projecile, which has boh horizonal and erical moion. These wo

More information

Outline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters?

Outline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters? of Topics wih MATLAB Shan He School for Compuaional Science Universi of Birmingham Module 6-3836: Compuaional Modelling wih MATLAB Wha will we learn from his lecure Aim of analsis: Aim of analsis. Some

More information

SOLUTIONS TO ECE 3084

SOLUTIONS TO ECE 3084 SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

Math Wednesday March 3, , 4.3: First order systems of Differential Equations Why you should expect existence and uniqueness for the IVP

Math Wednesday March 3, , 4.3: First order systems of Differential Equations Why you should expect existence and uniqueness for the IVP Mah 2280 Wednesda March 3, 200 4., 4.3: Firs order ssems of Differenial Equaions Wh ou should epec eisence and uniqueness for he IVP Eample: Consider he iniial value problem relaed o page 4 of his eserda

More information

Second Order Linear Differential Equations

Second Order Linear Differential Equations Second Order Linear Differenial Equaions Second order linear equaions wih consan coefficiens; Fundamenal soluions; Wronskian; Exisence and Uniqueness of soluions; he characerisic equaion; soluions of homogeneous

More information

Echocardiography Project and Finite Fourier Series

Echocardiography Project and Finite Fourier Series Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

6.003 Homework #9 Solutions

6.003 Homework #9 Solutions 6.00 Homework #9 Soluions Problems. Fourier varieies a. Deermine he Fourier series coefficiens of he following signal, which is periodic in 0. x () 0 0 a 0 5 a k sin πk 5 sin πk 5 πk for k 0 a k 0 πk j

More information

ln 2 1 ln y x c y C x

ln 2 1 ln y x c y C x Lecure 14 Appendi B: Some sample problems from Boas Here are some soluions o he sample problems assigned for Chaper 8 8: 6 Soluion: We wan o find he soluion o he following firs order equaion using separaion

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder#

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder# .#W.#Erickson# Deparmen#of#Elecrical,#Compuer,#and#Energy#Engineering# Universiy#of#Colorado,#Boulder# Chaper 2 Principles of Seady-Sae Converer Analysis 2.1. Inroducion 2.2. Inducor vol-second balance,

More information

EE243 Advanced Electromagnetic Theory Lec # 13: Waveguides and sources

EE243 Advanced Electromagnetic Theory Lec # 13: Waveguides and sources Applied M Fall 6, Neureuher Lecure #3 er /8/6 43 Advanced lecromagneic Theor Lec # 3: Waveguides and sources Source Free Region: ecor Poenials A and F Single direcion componen of A and F Give TM and T

More information

Non-uniform circular motion *

Non-uniform circular motion * OpenSax-CNX module: m14020 1 Non-uniform circular moion * Sunil Kumar Singh This work is produced by OpenSax-CNX and licensed under he Creaive Commons Aribuion License 2.0 Wha do we mean by non-uniform

More information

Chapter 7: Solving Trig Equations

Chapter 7: Solving Trig Equations Haberman MTH Secion I: The Trigonomeric Funcions Chaper 7: Solving Trig Equaions Le s sar by solving a couple of equaions ha involve he sine funcion EXAMPLE a: Solve he equaion sin( ) The inverse funcions

More information

( ) = b n ( t) n " (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2.

( ) = b n ( t) n  (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2. Andrei Tokmakoff, MIT Deparmen of Chemisry, 3/14/007-6.4 PERTURBATION THEORY Given a Hamilonian H = H 0 + V where we know he eigenkes for H 0 : H 0 n = E n n, we can calculae he evoluion of he wavefuncion

More information

EE100 Lab 3 Experiment Guide: RC Circuits

EE100 Lab 3 Experiment Guide: RC Circuits I. Inroducion EE100 Lab 3 Experimen Guide: A. apaciors A capacior is a passive elecronic componen ha sores energy in he form of an elecrosaic field. The uni of capaciance is he farad (coulomb/vol). Pracical

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion

More information

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

Today: Falling. v, a

Today: Falling. v, a Today: Falling. v, a Did you ge my es email? If no, make sure i s no in your junk box, and add sbs0016@mix.wvu.edu o your address book! Also please email me o le me know. I will be emailing ou pracice

More information

Chapter 10 INDUCTANCE Recommended Problems:

Chapter 10 INDUCTANCE Recommended Problems: Chaper 0 NDUCTANCE Recommended Problems: 3,5,7,9,5,6,7,8,9,,,3,6,7,9,3,35,47,48,5,5,69, 7,7. Self nducance Consider he circui shown in he Figure. When he swich is closed, he curren, and so he magneic field,

More information

2.5. The equation of the spring

2.5. The equation of the spring 2.5. The equaion of he spring I moun a spring wih a weigh on i. I ie he op of he spring o a sick projecing ou he op end of a cupboard door, and fasen a ruler down he edge of he door, so ha as he spring

More information

Applications of the Basic Equations Chapter 3. Paul A. Ullrich

Applications of the Basic Equations Chapter 3. Paul A. Ullrich Applicaions of he Basic Equaions Chaper 3 Paul A. Ullrich paullrich@ucdavis.edu Par 1: Naural Coordinaes Naural Coordinaes Quesion: Why do we need anoher coordinae sysem? Our goal is o simplify he equaions

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion

More information

Section 7.4 Modeling Changing Amplitude and Midline

Section 7.4 Modeling Changing Amplitude and Midline 488 Chaper 7 Secion 7.4 Modeling Changing Ampliude and Midline While sinusoidal funcions can model a variey of behaviors, i is ofen necessary o combine sinusoidal funcions wih linear and exponenial curves

More information

6.003 Homework #9 Solutions

6.003 Homework #9 Solutions 6.003 Homework #9 Soluions Problems. Fourier varieies a. Deermine he Fourier series coefficiens of he following signal, which is periodic in 0. x () 0 3 0 a 0 5 a k a k 0 πk j3 e 0 e j πk 0 jπk πk e 0

More information

Unit Root Time Series. Univariate random walk

Unit Root Time Series. Univariate random walk Uni Roo ime Series Univariae random walk Consider he regression y y where ~ iid N 0, he leas squares esimae of is: ˆ yy y y yy Now wha if = If y y hen le y 0 =0 so ha y j j If ~ iid N 0, hen y ~ N 0, he

More information

And the solution to the PDE problem must be of the form Π 1

And the solution to the PDE problem must be of the form Π 1 5. Self-Similar Soluions b Dimensional Analsis Consider he diffusion problem from las secion, wih poinwise release (Ref: Bluman & Cole, 2.3): c = D 2 c x + Q 0δ(x)δ() 2 c(x,0) = 0, c(±,) = 0 Iniial release

More information

Notes 04 largely plagiarized by %khc

Notes 04 largely plagiarized by %khc Noes 04 largely plagiarized by %khc Convoluion Recap Some ricks: x() () =x() x() (, 0 )=x(, 0 ) R ț x() u() = x( )d x() () =ẋ() This hen ells us ha an inegraor has impulse response h() =u(), and ha a differeniaor

More information

and v y . The changes occur, respectively, because of the acceleration components a x and a y

and v y . The changes occur, respectively, because of the acceleration components a x and a y Week 3 Reciaion: Chaper3 : Problems: 1, 16, 9, 37, 41, 71. 1. A spacecraf is raveling wih a veloci of v0 = 5480 m/s along he + direcion. Two engines are urned on for a ime of 84 s. One engine gives he

More information

Lecture #8 Redfield theory of NMR relaxation

Lecture #8 Redfield theory of NMR relaxation Lecure #8 Redfield heory of NMR relaxaion Topics The ineracion frame of reference Perurbaion heory The Maser Equaion Handous and Reading assignmens van de Ven, Chapers 6.2. Kowalewski, Chaper 4. Abragam

More information

This document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC

This document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC his documen was generaed a 1:4 PM, 9/1/13 Copyrigh 213 Richard. Woodward 4. End poins and ransversaliy condiions AGEC 637-213 F z d Recall from Lecure 3 ha a ypical opimal conrol problem is o maimize (,,

More information

Introduction to AC Power, RMS RMS. ECE 2210 AC Power p1. Use RMS in power calculations. AC Power P =? DC Power P =. V I = R =. I 2 R. V p.

Introduction to AC Power, RMS RMS. ECE 2210 AC Power p1. Use RMS in power calculations. AC Power P =? DC Power P =. V I = R =. I 2 R. V p. ECE MS I DC Power P I = Inroducion o AC Power, MS I AC Power P =? A Solp //9, // // correced p4 '4 v( ) = p cos( ω ) v( ) p( ) Couldn' we define an "effecive" volage ha would allow us o use he same relaionships

More information

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

More information

The Contradiction within Equations of Motion with Constant Acceleration

The Contradiction within Equations of Motion with Constant Acceleration The Conradicion wihin Equaions of Moion wih Consan Acceleraion Louai Hassan Elzein Basheir (Daed: July 7, 0 This paper is prepared o demonsrae he violaion of rules of mahemaics in he algebraic derivaion

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler EECE 3 Signals & Sysems Prof. Mark Fowler Noe Se # Wha are Coninuous-Time Signals??? /6 Coninuous-Time Signal Coninuous Time (C-T) Signal: A C-T signal is defined on he coninuum of ime values. Tha is:

More information

10.6 Parametric Equations

10.6 Parametric Equations 0_006.qd /8/05 9:05 AM Page 77 Secion 0.6 77 Parameric Equaions 0.6 Parameric Equaions Wha ou should learn Evaluae ses of parameric equaions for given values of he parameer. Skech curves ha are represened

More information

Linear Time-invariant systems, Convolution, and Cross-correlation

Linear Time-invariant systems, Convolution, and Cross-correlation Linear Time-invarian sysems, Convoluion, and Cross-correlaion (1) Linear Time-invarian (LTI) sysem A sysem akes in an inpu funcion and reurns an oupu funcion. x() T y() Inpu Sysem Oupu y() = T[x()] An

More information

Velocity is a relative quantity

Velocity is a relative quantity Veloci is a relaie quani Disenangling Coordinaes PHY2053, Fall 2013, Lecure 6 Newon s Laws 2 PHY2053, Fall 2013, Lecure 6 Newon s Laws 3 R. Field 9/6/2012 Uniersi of Florida PHY 2053 Page 8 Reference Frames

More information

2. The following diagram shows a circular loop of wire in a uniform magnetic field that points out of the page.

2. The following diagram shows a circular loop of wire in a uniform magnetic field that points out of the page. 1. Two elecromagneic waves ravel hrough emp space. Wave A as a wavelengh of 700 nm (red ligh), while Wave B has a wavelengh of 400 nm (blue ligh). Which saemen is rue? A) Wave A ravels faser because i

More information

Lab #2: Kinematics in 1-Dimension

Lab #2: Kinematics in 1-Dimension Reading Assignmen: Chaper 2, Secions 2-1 hrough 2-8 Lab #2: Kinemaics in 1-Dimension Inroducion: The sudy of moion is broken ino wo main areas of sudy kinemaics and dynamics. Kinemaics is he descripion

More information

1. VELOCITY AND ACCELERATION

1. VELOCITY AND ACCELERATION 1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under

More information

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k Challenge Problems DIS 03 and 0 March 6, 05 Choose one of he following problems, and work on i in your group. Your goal is o convince me ha your answer is correc. Even if your answer isn compleely correc,

More information

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit V() R L C 513 Elecrical Circuis Tools Used in Lab 13 Series Circuis Damped Vibraions: Energy Van der Pol Circui A series circui wih an inducor, resisor, and capacior can be represened by Lq + Rq + 1, a

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Hall effect. Formulae :- 1) Hall coefficient RH = cm / Coulumb. 2) Magnetic induction BY 2

Hall effect. Formulae :- 1) Hall coefficient RH = cm / Coulumb. 2) Magnetic induction BY 2 Page of 6 all effec Aim :- ) To deermine he all coefficien (R ) ) To measure he unknown magneic field (B ) and o compare i wih ha measured by he Gaussmeer (B ). Apparaus :- ) Gauss meer wih probe ) Elecromagne

More information

ES.1803 Topic 22 Notes Jeremy Orloff

ES.1803 Topic 22 Notes Jeremy Orloff ES.83 Topic Noes Jeremy Orloff Fourier series inroducion: coninued. Goals. Be able o compue he Fourier coefficiens of even or odd periodic funcion using he simplified formulas.. Be able o wrie and graph

More information