A Method for Recovering Arbitrary Graphs

Size: px
Start display at page:

Download "A Method for Recovering Arbitrary Graphs"

Transcription

1 A Mehod for Recovering Arbirary Graphs Nicolas Addingon January 4, 005 Absrac We develop a mehod for recovering arbirary graphs, a generalizaion of he sar k mehod ha avoids spurious parameers. They key is o recognize ha he inverse problem amouns o undoing he Schur complemen and o analyze he residue erm BC 1 B. Conens 1 Inroducion Graphs wih One Inerior Node.1 Basic Recovery Technique Example: The Kie Quadrilaerals: A Picorial Tool The Square Roo Trick The Sar K Mehod Wih Marices Wih Quadrilaerals Example: The Marshmallow A General Recovery Mehod Shorcomings of he Sar K Mehod Looking a Muliple Layers Recovery wih Larger Submarices Tesing for Singulariy Recovering an Arbirary Graph Example: Circles, 4 Rays Connecions and Deerminans Direcions for Fuure Research 11 6 Acknowledgemens 13 1

2 1 Inroducion We assume ha he reader is familiar wih he basics of Curis and Morrow s inverse problem [1]: graphs, Kirchhoff marices, response marices, and he Schur complemen. The forward problem is o ake a Kirchhoff marix K and find is response marix Λ. This is easily done using he Schur complemen: if we wrie K in he block form ( ) A B B C where C represens he inerior-oinerior connecions, Λ = A BC 1 B = K/C, he Schur complemen of C in K. A imes i will be convenien, albei awkward-sounding, o say K Schur C raher han he Schur complemen of C in K. I will also be convenien o refer o BC 1 B as he residue erm. The inverse problem is o ake a graph and is response marix and find he Kirchhoff marix, ha is, knowing K Schur C and knowing which enries of K are zero, o recover he remaining enries of K. The inverse problem, hen, amouns o undoing he Schur complemen. This is difficul because he Schur complemen desroys a lo of informaion: if a graph has eigh boundary nodes and eigh inerior nodes, we are rying o recover a marix from an 8 8 wih he help of a few zeros. The inverse problem has been solved for circular planar graphs [1] and annular graphs have been sudied some [3]; in his paper we are ineresed in he recoverabiliy of arbirary graphs. Our plan is o se up a sequence of inermediae marices beween he Kirchhoff marix and he response marix and pull he informaion in he response marix back hrough hese inermediae marices o he Kirchhoff marix. In we sudy he simples case, graphs wih a single inerior node. In 3 we describe he sar k mehod, which amouns o ieraing he one-node case. In 4 we develop a general recovery mehod, which is our main resul. In 5 we sugges direcions for fuure research. Graphs wih One Inerior Node.1 Basic Recovery Technique Consider a graph wih n boundary nodes and one inerior node. Fix an ordering of he nodes wih he boundary nodes firs. Le γ 1,...,γ n be he conducances of he edges joining he inerior node o he boundary nodes 1,...,n, le γ = (γ 1 γ n ) be he column vecor of hese conducances, and le σ = γ γ n be heir sum. Then when we wrie K in block form, A γ K = γ σ so Λ = A γγ σ. We wish o view his in he following manner: Λ is jus he upper lef corner of K superimposed wih a residue R = γγ σ from he eliminaion of he inerior node. Le us inroduce he noaion K for he upper lef corner of K (he size of which will always be clear from conex), so Λ = K +R. If we imagine laying hese marices fla, we can wrie R + K Λ. (.1.1) Recall ha K and Λ are symmeric, and observe ha R mus be as well. We are given all of Λ and he posiion of he zeros in K, so wherever here is a zero in K, we can read he corresponding enry of R direcly from Λ. Now γ 1 γ 1 γ γ 1 γ 3... γ 1 γ n R = γγ σ = 1 γ 1 γ γ γ γ 3... γ γ n γ 1 γ 3 γ γ 3 γ 3... γ γ n σ γ 1 γ n γ γ n γ 3 γ n... γ n

3 so any submarix of R has deerminan 0. Thus if we know hree enries in some submarix of R, we can recover he fourh. If here are sufficienly many well-placed zeros in K, we may be able o recover all of R his way. Once we have R, we can ge all bu he las row and column of K since K = Λ R. Since K is a Kirchhoff marix, is rows and columns sum o zero, so we can recover he las row and column as well.. Example: The Kie Consider he kie graph (Figure 1). We are given Λ = (λ ij ) and we wish o recover K. We know ha k 13, k 3, k 4, and k 34 are all zero since here are no edges beween he corresponding pairs of nodes, so we know he following enries of R:?? λ 13? R =?? λ 3 λ 4 λ 13 λ 3? λ 34.? λ 4 λ 34? Since r 13 r 4 r 14 r 3 = 0, we recover r 14 = λ13λ4 λ 3. We can recover all of R similarly: R = λ 13 λ 4 λ 13λ 4 λ λ 3λ 34 λ 34 λ 13λ 4 13 λ 3 λ 13λ 4 λ 34 λ 3λ 4 λ 34 λ 3 λ 4 λ 13 λ 3 λ 3λ 34 λ 4 λ 34 λ 13λ 4 λ 3 λ 4 λ 34 λ 4λ 34 λ 3 Observe ha here are several ways o recover a given enry, for example r 11 = r13r14 r 34 = r1r13 r 3 = r14 r 44, bu hese all give he same resul. (Observe also ha he indices cancel in each expression.) I appears ha in general, if we begin wih a se of enries in R which is independen, ha is, if no enry in our se can be goen from he ohers, hen no maer how we recover he remaining enries we ge he same values. Now ha we have R, we can recover he upper lef corner of K: λ 11 λ13 λ 4 λ 3λ 34 λ 1 λ13λ4 λ 34 0 λ 14 λ13λ4 λ 3? λ 1 λ13λ4 λ 34 λ λ3λ4 λ ? K = 0 0 λ 33 λ3λ34 λ 4 0? λ 14 λ13λ4 λ λ 14 λ4λ34 λ 3?????? and since he rows and columns of K sum o 0, we can recover all of K. Observe ha k 1 = λ 1 λ 13λ 4 λ 1 λ 13 λ 4 λ 34 = λ 34 and ha he deerminan in he las expression corresponds o he -connecion from nodes 1 and 4 o and 3. Since here is an edge beween nodes 1 and, k 1 is negaive, so his deerminan is posiive, and he connecion exiss. Removing his edge would break his connecion, bu i would also make k 1 = 0 and hence make his deerminan 0. Thus we see he link beween -connecions and deerminans very clearly by sudying graphs wih one inerior node..3 Quadrilaerals: A Picorial Tool We can idenify a nework (i.e. a graph wih conducances) wih is Kirchhoff marix, and a response marix is he Kirchhoff marix of a complee graph, so we can undersand aking he Schur complemen of he Kirchhoff marix as replacing he inerior nodes and he edges inciden o hem wih a response-equivalen complee graph (Figure ). λ 34. 3

4 For graphs wih one inerior node, we can idenify he superimposed complee graph wih he off-diagonal enries of he residue R, since R = Λ excep on he diagonal if here are no boundary-o-boundary connecions. Thus our saemen Λ = K + R can be inerpreed picorially as in Figure 3. In his inerpreaion, a submarix of R ( rik ) r il r jk r jl corresponds o a quadrilaeral in he superimposed complee graph as in Figure 4, so if we know hree sides of his quadrilaeral we can recover a fourh. Quadrilaerals provide a quick way o check he recoverabiliy of a graph. To show ha he kie is recoverable, we begin by drawing he picure in Figure 5. Where here are doubled edges, knowing he response marix gives us heir sum, so we need o separae ha number ino he par ha came from he boundary-o-boundary connecion and he par ha came from he superimposed complee graph. For he edge 1, we know he oher hree sides of he quadrilaeral 1 3 4, and he quadrilaeral gives us he edge 1 4. Now we know he boundary-o-boundary connecions and he superimposed complee graph, i.e. he off-diagonal enries of R, so we know we can recover he diagonal R, so he kie is recoverable. If insead we were rying o recover he bowie, we would draw he picure in Figure 6. Any quadrilaeral wih 1 4 as a side also has 3 as a side, so we canno ge hree sides of any quadrilaeral o recover 1 4 or 3, and he bowie is no recoverable..4 The Square Roo Trick In.1 we said ha since any submarix of R is singular, we can recover an enry in a submarix given he oher hree. Since R is symmeric, however, more is rue. If our submarix is principal and we know only he wo diagonal enries, we can recover he off-diagonal enries: 0 = r ii r ij = r iir jj rij, r ij r jj so r ij = r ii r jj ; he enries of R are necessarily negaive. Observe again ha he indices cancel. The Sar of David graph (Figure 7) can be recovered using his square roo rick bu no wihou i. 3 The Sar K Mehod 3.1 Wih Marices In his secion we shall assume ha our graph has 1 boundary nodes and 4 inerior nodes o avoid drowning he reader in a sea of index variables; wha is mean in general should be clear. Le K 16 = K be he Kirchhoff marix of he original nework. Le K 15 equal K 16 Schur is lower righ 1 1 corner. Le R 16 be he residue γγ σ from his Schur complemen, i.e. he residue from he eliminaion of node 16, so K 15 = K 16 + R 16. Similarly, le K 14 equal K 15 Schur is lower righ 1 1 corner and R 15 be he residue, so K 14 = K 15 + R 15. Coninue in his fashion down o K 1. Proposiion 1. K 1 equals K 16 Schur is lower righ 4 4 corner, ha is, K 1 = Λ. Proof. In [], Crabree and Haynsworh show ha if C is an inverible square submarix of M and D is an inverible square submarix of D hen M/D = (M/C)/(C/D). I is easy o check ha he lower righ 1 1 corner of K 15 equals he lower righ corner of K 16 Schur is lower righ 1 1 corner, so K 14 equals K 16 Schur is lower righ corner. The desired resul is obained by inducion. We are given he graph, or equivalenly, he zeros of K = K 16. From his we can deermine he zeros of R 16 and hence he zeros of K 15, and so on down o K 13. Given Λ = K 1 and he zeros of K 13 we may be able, using our one-inerior-node mehod, o recover R 14 and hence K 14. From his we may be able o recover K 15, and so on back up o K 16. If his process succeeds, we know ha he graph is recoverable. If a some poin during his process we ge suck, ha is, we know a few enries of some R n bu canno recover any more using submarices, we can paramerize one of he unknown enries and coninue. The 4

5 number of parameers gives a measure of he irrecoverabiliy of our graph. I may also happen ha a one sage we do no have enough informaion o coninue, so we inroduce a parameer, bu a a laer sage we have more informaion han we need and are able o use submarices o eliminae he parameer. An example will be given in Wih Quadrilaerals In.3, we replaced he inerior node and he edges inciden o i, or he sar around he inerior node, wih he response-equivalen complee graph, or k. Now ha we are dealing wih several inerior nodes, we make his sar k replacemen several imes and obain a sequence of inermediae graphs as in Figure 8. The K n of he previous secion are he Kirchhoff marices of hese inermediae graphs; specifically, K n is he inermediae graph wih n nodes. If each inermediae graph is recoverable, he original graph is recoverable. Here again, quadrilaerals provide a quick check of recoverabiliy. Consider he hexagons on a erahedron in Figure 9. Wih he quadrilaeral we can separae he wo edges 1, and by symmery all he doubled edges, so we can recover all four inerior nodes. This example is ineresing because i is a flower: i has neiher a boundary-o-boundary connecion nor a boundary spike (a boundary node conneced o a single inerior node). Curis and Morrow have shown ha no circular planar flower is recoverable; his flower is, of course, no circular planar. Jeff Russell and Tracy Lovejoy s oy drum (Figure 10) is anoher example of a recoverable flower. 3.3 Example: The Marshmallow Consider he marshmallow graph (Figure 11), so called because he presenaion in which he auhor inroduced i involved a marshmallow model. For his graph,? ? 0 0? 0 0 0??? ? 0 0? 0 0?? 0????? λ 11 λ 1 λ 13 λ 14 λ 15 K 7 = 0 0 0? 0?? K 6 = 0????? λ 1 λ λ 3 λ 4 λ ??? 0????? K 5 = λ 13 λ 3 λ 33 λ 34 λ 35?????? 0 0????? λ 14 λ 4 λ 34 λ 44 λ 45 λ?????? 15 λ 5 λ 35 λ 45 λ 55 0???? 0? so? λ 1 λ 13 λ 14 λ 15 λ 1???? R 6 = λ 13???? λ 14????. λ 15???? No furher enries of R 6 can be recovered, bu if we paramerize he 1,1 enry by < 0, we can ge everyhing else: λ 1 λ 13 λ 14 λ 15 λ λ 1 λ 1λ 13 λ 1λ 14 λ 1λ 15 1 λ R 6 = λ 1λ 13 λ 13 λ 13λ 14 λ 13λ λ 13λ 14 λ 14 λ 14λ 15 λ λ 1λ λ λ 1λ λ 13λ 15 λ 14λ 15 From here we can ge K 6 and complee he recovery. Along he way, however, an ineresing hing happens. In R 5, rows and 3 by columns 4 and 5 are ( ) λ4 λ1λ14 λ 5 λ1λ15 = 1 λ 1 λ 14 λ 4 λ 1 λ 15 λ 5 λ 34 λ13λ14 λ 35 λ13λ15 λ 13 λ 14 λ 34 λ 13 λ 15 λ 35 λ 15 5

6 bu since his is a submarix of an R n, i is singular, so by he six-erm ideniy, 0 = 1 λ 1 λ 14 λ 4 λ 1 λ 15 λ 5 λ 13 λ 14 λ 34 λ = 1 λ 14 λ 15 λ 1 λ 4 λ 5 13 λ 13 λ 34 λ 35. λ 15 λ 35 If ( λ 4 λ 5 ) λ 34 λ 35 is inverible, his allows us o solve for, ha is, was a spurious parameer. Observe ha ( λ 4 λ 5 ) λ 34 λ 35 corresponds o he -connecion from nodes and 3 o 4 and 5. This connecion exiss, bu he corresponding deerminan may be zero because we can also connec and 3 o 5 and 4 (he opposie permuaion). For mos graphs in which spurious parameers arise, i is possible o decide wheher he necessary deerminan is inverible or no by examining connecions; he marshmallow is somehing of a pahalogical case in his respec. If ( λ 4 λ 5 ) ( λ 34 λ 35 is inverible, he marshmallow is recoverable. By symmery, if λ3 λ 5 ) ( λ 34 λ 45 or λ3 λ 4 ) λ 35 λ 45 is inverible, i is also recoverable. If all hree of hese are singular, i is possible o exhibi a one-parameer family of Kirchhoff marices wih he same response marix. The marshmallow, hen, is almos always recoverable. Tracy s J (Figure 1) has he same propery. I is ineresing o compare his sor of graph wih he -o-1 graphs sudied by Ernie Esser and Tracy Lovejoy [5], which are almos never recoverable he former are recoverable unless cerain submarices of Λ are singular, whereas he laer are recoverable only if cerain submarices of Λ are singular. 4 A General Recovery Mehod 4.1 Shorcomings of he Sar K Mehod The sar k mehod has wo main shorcomings. Firs, i only works unequivocally when each of he inermediae graphs is recoverable. This is ofen rue for small graphs, such as he opha (Figure 8), and graphs wih few connecions, such as hexagons on a erahedron, bu i fails for many denser graphs. There is no order of inerior nodes for which he sar k mehod can recover Ernie Esser s circles, 4 rays (Figure 13a) or even he well-conneced graph on 5 nodes shown in Figure 14a, which is circular planar. Boh hese graphs are known o be recoverable by oher means. Second, he sar k mehod is very sensiive o he order in which he inerior nodes are eliminaed. Figure 15a is recoverable using he sar k mehod, bu Figure 15b is no. In he previous secion, he recovery of he marshmallow required one parameer, bu if we swich he order of nodes 6 and 7, i requires four parameers (hree of which end up being spurious). Parameers which come and go provide a way around hese shorcomings, and i may be ha all recoverable graphs can be recovered using he sar k mehod and spurious parameers. This mehod is inelegan, however, and does no reflec wha is really going on. If a spurious parameer arises in he sar k recovery of a graph, i is because one of he inermediae graphs was no recoverable bu he original graph was; some informaion was presen in he chain of inermediae graphs ha could no be seen by looking only from one graph o he nex. In he previous secion we saw how when we look across several seps, 3 3 deerminans arise from deerminans of deerminans by way of he six-erm ideniy. In his secion we will use 3 3 and larger deerminans direcly, wihou he help of spurious parameers. 4. Looking a Muliple Layers In.1 we wroe Λ = K + R, or when we laid K and R fla, R + K Λ We were given cerain enries of he marices in his sum: he zeros of K and R and all he enries of Λ. We used verical informaion knowing k ij and λ ij we could find r ij, and knowing r ij and λ ij we could find. 6

7 k ij and horizonal informaion he singulariy of cerain submarices of R and he fac ha he rows and columns of K summed o zero o recover all he enries of everyhing, and paricular o recover K. In 3.1 we had or if we lay he marices fla, K 1 = K 13 + R 13 = K 14 + R 14 + R 13 = K 15 + R 15 + R 14 + R 13 = K 16 + R 16 + R 15 + R 14 + R 13 R 13 R 14 R 15 R 16 + K 16 K 1. (4..1) Here we did he same hing: we began wih he known zeros of K 16 (and hence of he R n ) and he known enries of Λ = K 1 and used verical and horizonal informaion o peel off he R n one a a ime and recover K 16. The following sor of hing happens ofen: he i,j enries of K 16, R 16, and R 15 are zero bu he i,j enries of R 13 and R 14 are no, so we can recover he sum of hese wo enries, bu we do no know wha porion of he number comes from R 13 and wha porion from R 14. To make use of his informaion, we inroduce a marix R14 13 = R 14 + R 13 ; we ypically know more abou his marix han we do abou R 13 and R 14 separaely. This is our plan: in addiion o considering he single layers R 13,...,R 16, we will also consider all sums of coniguous blocks of hem, for example R15 13 = R 15 + R 14 + R 13 and R16 15 = R 16 + R 15 and even K 14 = K 16 + R 16 + R 15. These muli-layer gadges add up nicely: R16 13 = R R 13 = R R14 13 = R 16 + R Thus we undersand he verical informaion for he R m n. To recover a graph, we wish o use horizonal informaion abou he R m n o recover heir enries from a few given enries. For single layers, we used he fac ha submarices were singular. For muliple layers, we will use larger submarices. 4.3 Recovery wih Larger Submarices The following Lemma was used acily in 3.3 and will be used exensively in wha follows. Lemma. If an n n marix M is singular, we know all bu one enry m ij, and he cofacor M ij is inverible, we can recover he unknown enry. Proof. By cofacor expansion along he i h row, 0 = dem = ( 1) i+1 m i1 de M i1 + + ( 1) i+j m ij de M ij + + ( 1) i+n m in dem in, so since de M ij 0, we can solve for m ij. In.4, we used he symmery of R o do slighly beer for marices. I is unclear wheher i is worh doing his for larger marices. Suppose ha a square marix M is singular and symmeric, wih he block form a B x M = B C D x D e 7

8 where he corners a, x, and e are 1 1. We wish o solve for x. By he six-erm ideniy, 0 = M C = a B B C C = a B B C C = a B B C C D D D D e B x C D B C x D D e B x C D D ( e B 0 C D + ( 1)n x C ) so if C is inverible, x = ( 1)n+1 C ( B 0 C D ± a B B C C D D ) e. If he plus-or-minus erm is nonzero, we have wo possible values for x. If one is posiive and one negaive, in our applicaion we will know which one we wan, bu here is no obvious way of elling wheher he wo soluions are of he same or opposie sign. On he oher hand, if he plus-or-minus erm is zero, B C D x = 0, so we can find x from ha. I is no known wheher any graphs can be recovered using his generalized square roo rick bu no wihou i. We guess ha here are none, bu suspec ha his rick is a work in -o-1 graphs. 4.4 Tesing for Singulariy To make use of Lemma, we need o be able o es submarices of he Rn m for singulariy. We coninue o assume ha our graph has 16 nodes. We defined R 15 (for example) as he residue erm γγ σ from K 15 Schur is lower righ 1 1 corner. We wish o show ha R15 14 (for example), which we defined as R 15 +R 14, is he residue from K 15 Schur is lower righ corner. Le C15 14 be ha corner: rows 14 o 15 by columns 14 o 15 of K 15. Le K 15 be rows 1 o 13 by columns 1 o 13 and B be rows 1 o 13 by columns 14 o 15, so K 15 B K 15 =. B C15 14 Now K 13 = K 15 + R 15 + R 14 = K 15 + R15, 14 bu by Proposiion 1, K 13 = K 15 B (C15) 14 1 B as well, so R15 14 = B (C15) 14 1 B, i.e. he residue from he eliminaion of nodes 15 and 14. K 13 is he Schur complemen of C15 14 in K 15. R15 14 is he Schur complemen of C15 14 in a closely relaed marix: wih B as above, define 0 B Z15 14 =, B C15 14 ha is, K 15 wih he upper lef corner replaced wih zeros. Now R = 0 B (C 14 15) 1 B = Z 14 15/C Proposiion 3. Le M be he submarix of R15 14 consising of rows r 1,...,r n by columns c 1,...,c n. Le N be he submarix of Z15 14 consising of rows r 1,...,r n,14,15 by columns c 1,...,c n,14,15. Then M is singular if and only if N is. Proof. M = N/C15, 14 so de M = (den)/(de C15). 14 C15 14 is no singular since i is a principal proper submarix of a Kirchhoff marix. 8

9 In oher words, o decide wheher a square submarix M of R15 14 is singular, we can ake he same submarix of Z15 14 (which is jus K 15 wih he upper lef corner suppressed), ack on he rows and columns ha would be chopped off in aking he Schur complemen of C15, 14 and decide if ha is singular. We know which enries of K 15 are zero, which are posiive (he diagonal enries), and which are negaive (he off-diagonal enries ha are no zero). We can ry o decide wheher a submarix N of Z15 14 is singular as follows. Expand he deerminan as a polynomial in he enries of N (hence of K 15 ). If he erms are all zero, N is singular. If he erms are all posiive or all negaive, N is inverible. If some of he erms are posiive and some negaive, he signs of he enries alone are no enough o decide. 4.5 Recovering an Arbirary Graph Now we are in a posiion o ouline a mehod for recovering he Kirchoff marix of an arbirary graph from is response marix. Suppose our graph has M boundary nodes and N nodes alogeher. 1. Wrie down he signs of all he enries of K N, which we know from he graph. From hese, deermine he signs of all he enries of K M,...,K N 1. We will use hese o es submarices for singulariy using Proposiion 3.. Make empy marices K M,...K N, R M+1,...,R N, and R m n, M < m < n N of he appropriae sizes: K n is n n, and R n and R m n are n 1 n 1. Fill in he zeros of all hese, which can be derived from he zeros of K N. Fill in he enries of K M, he response marix. 3. Whenever we know wo of hree enries from somehing of he form R16 13 = R R 13 K 13 = K 15 + R15, 14 recover he hird. 4. Whenever we know all bu one enry of a submarix of any marix, if he submarix is singular and he cofacor of he unknown enry is inverible, recover he unknown enry using Lemma. Also use he basic square roo rick from Whenever we know all bu one enry in a row of a K n, recover i using he fac ha he rows of Kirchhoff marices sum o zero. 6. If a any poin no more enries can be recovered bu some are sill missing, paramerize an unknown enry. The firs single layer R n wih unknown enries (firs in he sense ha n is leas) seems o be he bes place o parameerize. Observe ha if we resric our recovery o he single layers R n and omi he muli-layers R m n, his is exacly he sar k mehod. This mehod, while powerful (here is no known recoverable graph ha i fails o recover), is impracical o work by hand for large graphs, as we shall see in he nex secion, and unnecessary for small graphs, where he sar k mehod usually works. I is ideally suied o a compuer, however. The auhor has implemened he mehod as a C++ program, which is available a hp:// reu/. 4.6 Example: Circles, 4 Rays Consider circles, 4 rays (Figure 13a), which has 8 boundary nodes and 8 inerior nodes. An ad hoc proof of his graph s recoverabiliy is given by Ernie Esser in [3]. Here we recover i using our general mehod. These are he signs of he enries in K 16, he Kirchhoff marix: K 16 = or

10 We are given K 8, he response marix. Firs we wish o recover R 9, he residue from he eliminaion of node 9. This is an 8 8 marix, all of whose enries are negaive. Now R = K 9 = K 16 + R = so we know he following enries of R 9 = K 8 K 9 :???? λ 15??????? λ 5??????? λ 35??? R 9 = K 8 K 9 =???? λ 45??? λ 15 λ 5 λ 35 λ 45? λ 56 λ 57 λ 58.???? λ 56??????? λ 57??????? λ 58??? To recover R 9, i suffices o recover he 5,5 enry, for hen we can he res using submarices as we did in.. We will recover he 5,5 enry by way of R16. 9 Consider he submarix M of R16 9 consising of rows 1,, 4, 5 by columns 3, 5, 7, 8. Since R16 9 = K 8 K 16 and he submarix 1,,4,5 3,5,7,8 of K 16 is zero excep a 5,5, we have λ 13 λ 15 λ 17 λ 18 M = λ 3 λ 5 λ 7 λ 8 λ 34 λ 45 λ 47 λ 48 λ 35? λ 57 λ 58 I suffices o show ha 1,,4,5 3,5,7,8 is singular bu 1,,4 3,7,8 is inverible, for hen we can recover he 5,5 enry by Lemma ; hen since R 9 = R16 9 R16, 10 and he 5,5 enry of R16 10 is 0, we can recover he 5,5 enry of K 9. From 4.4 we know ha he submarix 1,,4 3,7,8 of R16 9 is inverible if and only if he submarix 1,,4,9 16 3,7,8,9 16 of K 16 wih he upper lef 3 3 corner changed o zeros is inverible: By cofacor expansion along he op row, his submarix is inverible if and only if he submarix,4,9 16 3,7,8,9 1,14 16 is inverible; ha is, since here is a single nonzero enry in he op row, we can cross ou row 1 by column 13 (he reader is encouraged o ge a pencil and do his). Similarly, we can cross ou 14, 4 16, 15 3 (going down he lef-hand column now), 11 7, 1 8, 13 9, and This leaves us wih he submarix 9, 14, 16 10, 1, 15: which is necessarily inverible, so he submarix 1,,4 3,7,8 of R 9 16 is inverible. For 1,,4,5 3,5,7,8, we play he same game, bu afer crossing ou a few hings we find a row of all zeros, so ha marix is singular. Thus R 9 is recoverable, so K 9 is recoverable. The recovery of K 10,...,K 16 is similar. Exercise. Show ha he well-conneced graph on 5 nodes (Figure 14a) is recoverable. (Hin: Recover he 1,1 enry of R 6 9 using he submarix 1,,3 1,4,5.) 10

11 4.7 Connecions and Deerminans In he previous example, we saw ha a cerain submarix of R 9 16 was inverible because a relaed submarix of Z 9 16 was. Because of he placemen of zeros in K 16, he submarix of R 9 16 was equal o a submarix of K 8 and he submarix of Z 9 16 was equal o a submarix of K 16. Bu we know ha he subdeerminan 1,,4 3,7,8 of K 8, he response marix, corresponds o a 3-connecion from nodes 1,, and 4 o 3, 7, and 8. The calculaion we jus made shows us how. The reader is encouraged o mark Figure 13a wih a pencil. Firs we crossed ou row 1 by column 13. A 3-connecion from 1,, 4 o 3, 7, 8, mus include he edge from node 1 o node 13; from node 1, we canno go anywhere else. Similarly, we mus go from o 14, 4 o 16, 15 o 3, 11 o 7, and 1 o 8. Now we could go from 13 o 14, 16, or 9, bu if we wen o 14 or 16 we would collide wih he connecions coming from and 4, so we mus in fac go o 9. Similarly, we mus ge o 11 from 10. Now here are only wo ways o complee he connecion: eiher 9 o 10, 14 o 15, and 16 o 1, or 9 o 1, 14 o 10, and 16 o 15; hese correspond o he wo non-zero erms in he deerminan of How subdeerminans of Rn m ranslae ino connecions like his has no been sudied carefully. We can say a few hings. The lower index n indicaes he inermediae graph hrough which he connecion is going; subdeerminans of R14, 9 for example, correspond o connecions no hrough he original graph bu hrough he inermediae graph wih 14 nodes (where nodes 15 and 16 have been sar k d ). The upper index m indicaes how many nodes are considered boundary nodes; in R16, 9 nodes 1 8 were boundary nodes, bu for R16, 11 nodes 1 10 would ac as boundary nodes. The fac ha any submarix of a single layer, say R 11, is singular reflecs he fac ha here are no -connecions hrough a sar. Similarly, for an Rn m comprising k = n m + 1 layers, any k + 1 k + 1 submarix is singular since here are no k + 1 connecions hrough a graph wih k inerior nodes. There are many quesions, however. Wha does i mean o replace he upper righ corner of a Kirchhoff marix wih zeros? If he zeros in K 16 had no been so foruiously placed, how would a subdeerminan of R16 9 differ from a subdeerminan of K 8? Wha does i mean when a subdeerminan inersecs he diagonal does i make sense o speak of he connecion from nodes 1,, and 3 o 1, 4, and 5, as we seem o have been doing in 3.3? A circular planar graph is recoverable if and only if i is criical, ha is, if deleing or conracing any edge (replacing is conduciviy wih 0 or ) would break some connecion. Wha is he analoguous resul for non-planar graphs? To recover a piece of informaion wih our mehod, we need some deerminan o be zero, bu become nonzero when we crossed off he unknown enry s row and column (i.e. consider he unknown enry s cofacor); his mus be closely relaed o criicaliy. 5 Direcions for Fuure Research Our sudy on undoing he Schur complemen raises many quesions. 1. If our mehod succeeds, he graph is recoverable. Is he converse rue? Are here any recoverable graphs for which his mehod fails?. If we have o inroduce a parameer, does i maer where we inroduce i? Suppose we can recover some bu no all enries of he K n, R n, and R m n. We paramerize an unknown enry α, which allows us o recover one or more addiional enries β. We would hope ha paramerizing any β would have allowed us o recover α; ha is, any enry in he se we gain is as good as any oher. This does no appear o be he case. If we paramerize enries of single layers a he op of he pile (4..1), say R 14 if R 13 is compleely recovered, his seems o give us more han if we paramerize near he boom, say in K 16. Things propogae down beer han hey propogae up. Thus we suspec ha we have overlooked some informaion available o us, jus as he sar k mehod overlooked he muli-layers R m n. 11

12 One soluion may be o look a hings like his. Using he noaion of.1, ) ( γiγ j σ γ iγ k σ is singular. Half of his marix comes from R and half from K; i is as hough siing in R in (.1.1), we have peered over he ledge, down ino K. Now γ j is he sum of he j h row of R, so more generally, we can consider sums of enries in R: he marix ( ) rik + r il r im γ j γ k r jk + r jl r jm is singular, for example, and as we have seen, we ofen know more abou he sum han abou he summands, and invesigaing his can be fruiful. Some hings can be said abou he row sums of he R m n as well. We do no wan o use hese facs piecemeal, however. Wha is he general phenomenon ha is a work here? 3. Is any of he informaion we are using redundan? Can our recovery mehod be simplified? In paricular, is i necessary o consider submarices of he inermediae K n or do he he R m n suffice? Are here any graphs ha can be recovered by looking a he K n bu no by looking a he R m n alone? 4. Is our recovery mehod independen of he order of inerior nodes? Tha is, if i succeeds for one order, does i succeed for all possible orders? This appears o be rue. Some orders are nicer han ohers, however. I is much easier o recover circles, 4 rays wih he order shown in Figure 13b (due o Jeff Russell) han wih he usual order (Figure 13a). Using he sar k mehod, he former requires fewer spurious parameers. Figure 14b is similar. Why is his? Wha makes one order nicer han anoher? How should one choose an ordering of he inerior nodes o make he recovery as fas as possible? 5. During he recovery process, mos of he work is done a he surface, in he single-layer R n s, bu we usually need a leas one deerminan buried in he hickes R m n s. In 4.6, his was he submarix 1,,4,5 3,5,7,8 of R 9 16, which is a residue marix eigh layers hick. Currenly, he only way o find hese is for a human o have some inuiion abou he graph or for he compuer o search exhausively. Is here some way o make his search smarer? How can we know in advance which deeply buried submarices will be helpful? Where should we look for hem? 6. If we have o inroduce parameers, i appears ha he graph is no recoverable. However, for -o-1 graphs, which are almos recoverable, he parameer is almos spurious presumably some deerminan somewhere forces i o ake only finiely many values. Where exacly does his -o-1 ness come from in he R m n? Wha role does he generalized square roo rick of 4.3 play? Are here any 3-o-1 graphs? 7. The quesions abou deerminans, connecions, and criical graphs from he end of he previous secion. 8. Aside from he square roo rick, we did no use he symmery of our marices. How easily can our mehod be applied o direced graphs? How can we undersand ordinary (undireced) graphs as direced graphs wih some addiional consrains, jus like layered neworks? Wha abou verex conduciviy neworks? 9. Michael Goff makes an exensive sudy of neworks wih negaive conduciviies in [4]. The Schur complemen sill plays a key role for such neworks, bu hey can be quie pahalogical principal proper submarices of Kirchhoff marices need no be inverible, for example. Our mehod relied heavily on he signs of he enries of Kirchhoff marices being nice. Can i be modified o work wih negaive conduciviies? 10. Wha does our mehod allow us o say abou special classes of graphs, such as hree-dimensional laices? 1

13 6 Acknowledgemens The work on he sar k mehod was done in collaboraion wih Jeff Russell, whose exposiion can be found in [6]. The idea of inroducing parameers is due o Jeff Giansiracusa. The auhor would also like o hank Sam Coskey and Jim Morrow for many helpful discussions and suggesions. References [1] E. B. Curis and J. A. Morrow, Inverse problems for elecrical neworks, World Scienific, (000). [] D. Crabree and E. Haynsworh, An ideniy for he Schur complemen of a marix, Proc. Amer. Mah Soc. (1969), [3] Ernie Esser, On solving he inverse conduciviy problem for annular neworks, (000). [4] Michael Goff, Recovering neworks wih signed conduciviies, (003). [5] Tracy Lovejoy, Applicaions of he sar k ool, (003). [6] Jeff Russell, Sar and K solve he inverse problem, (003). 13

14

15

16

17

18

19

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

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

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

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n Lecure 3 - Kövari-Sós-Turán Theorem Jacques Versraëe jacques@ucsd.edu We jus finished he Erdős-Sone Theorem, and ex(n, F ) ( /(χ(f ) )) ( n 2). So we have asympoics when χ(f ) 3 bu no when χ(f ) = 2 i.e.

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

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

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

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

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

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

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

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

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

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

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

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

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

THE MATRIX-TREE THEOREM

THE MATRIX-TREE THEOREM THE MATRIX-TREE THEOREM 1 The Marix-Tree Theorem. The Marix-Tree Theorem is a formula for he number of spanning rees of a graph in erms of he deerminan of a cerain marix. We begin wih he necessary graph-heoreical

More information

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu ON EQUATIONS WITH SETS AS UNKNOWNS BY PAUL ERDŐS AND S. ULAM DEPARTMENT OF MATHEMATICS, UNIVERSITY OF COLORADO, BOULDER Communicaed May 27, 1968 We shall presen here a number of resuls in se heory concerning

More information

Math 333 Problem Set #2 Solution 14 February 2003

Math 333 Problem Set #2 Solution 14 February 2003 Mah 333 Problem Se #2 Soluion 14 February 2003 A1. Solve he iniial value problem dy dx = x2 + e 3x ; 2y 4 y(0) = 1. Soluion: This is separable; we wrie 2y 4 dy = x 2 + e x dx and inegrae o ge The iniial

More information

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

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

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

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

More information

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively: XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of

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

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

Homework sheet Exercises done during the lecture of March 12, 2014

Homework sheet Exercises done during the lecture of March 12, 2014 EXERCISE SESSION 2A FOR THE COURSE GÉOMÉTRIE EUCLIDIENNE, NON EUCLIDIENNE ET PROJECTIVE MATTEO TOMMASINI Homework shee 3-4 - Exercises done during he lecure of March 2, 204 Exercise 2 Is i rue ha he parameerized

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

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

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

Chapter 6. Systems of First Order Linear Differential Equations

Chapter 6. Systems of First Order Linear Differential Equations Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

Some Ramsey results for the n-cube

Some Ramsey results for the n-cube Some Ramsey resuls for he n-cube Ron Graham Universiy of California, San Diego Jozsef Solymosi Universiy of Briish Columbia, Vancouver, Canada Absrac In his noe we esablish a Ramsey-ype resul for cerain

More information

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,

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

5.1 - Logarithms and Their Properties

5.1 - Logarithms and Their Properties Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We

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

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

= ( ) ) or a system of differential equations with continuous parametrization (T = R

= ( ) ) or a system of differential equations with continuous parametrization (T = R XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of

More information

Math 115 Final Exam December 14, 2017

Math 115 Final Exam December 14, 2017 On my honor, as a suden, I have neiher given nor received unauhorized aid on his academic work. Your Iniials Only: Iniials: Do no wrie in his area Mah 5 Final Exam December, 07 Your U-M ID # (no uniqname):

More information

6.01: Introduction to EECS I Lecture 8 March 29, 2011

6.01: Introduction to EECS I Lecture 8 March 29, 2011 6.01: Inroducion o EES I Lecure 8 March 29, 2011 6.01: Inroducion o EES I Op-Amps Las Time: The ircui Absracion ircuis represen sysems as connecions of elemens hrough which currens (hrough variables) flow

More information

Notes for Lecture 17-18

Notes for Lecture 17-18 U.C. Berkeley CS278: Compuaional Complexiy Handou N7-8 Professor Luca Trevisan April 3-8, 2008 Noes for Lecure 7-8 In hese wo lecures we prove he firs half of he PCP Theorem, he Amplificaion Lemma, up

More information

Math 116 Practice for Exam 2

Math 116 Practice for Exam 2 Mah 6 Pracice for Exam Generaed Ocober 3, 7 Name: SOLUTIONS Insrucor: Secion Number:. This exam has 5 quesions. Noe ha he problems are no of equal difficuly, so you may wan o skip over and reurn o a problem

More information

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

Explore 2 Proving the Vertical Angles Theorem

Explore 2 Proving the Vertical Angles Theorem Explore 2 Proving he Verical Angles Theorem The conjecure from he Explore abou verical angles can be proven so i can be saed as a heorem. The Verical Angles Theorem If wo angles are verical angles, hen

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

) were both constant and we brought them from under the integral.

) were both constant and we brought them from under the integral. YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates Biol. 356 Lab 8. Moraliy, Recruimen, and Migraion Raes (modified from Cox, 00, General Ecology Lab Manual, McGraw Hill) Las week we esimaed populaion size hrough several mehods. One assumpion of all hese

More information

Explaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015

Explaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015 Explaining Toal Facor Produciviy Ulrich Kohli Universiy of Geneva December 2015 Needed: A Theory of Toal Facor Produciviy Edward C. Presco (1998) 2 1. Inroducion Toal Facor Produciviy (TFP) has become

More information

Sections 2.2 & 2.3 Limit of a Function and Limit Laws

Sections 2.2 & 2.3 Limit of a Function and Limit Laws Mah 80 www.imeodare.com Secions. &. Limi of a Funcion and Limi Laws In secion. we saw how is arise when we wan o find he angen o a curve or he velociy of an objec. Now we urn our aenion o is in general

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

1 Review of Zero-Sum Games

1 Review of Zero-Sum Games COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any

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

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

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

More information

Math From Scratch Lesson 34: Isolating Variables

Math From Scratch Lesson 34: Isolating Variables Mah From Scrach Lesson 34: Isolaing Variables W. Blaine Dowler July 25, 2013 Conens 1 Order of Operaions 1 1.1 Muliplicaion and Addiion..................... 1 1.2 Division and Subracion.......................

More information

Optimality Conditions for Unconstrained Problems

Optimality Conditions for Unconstrained Problems 62 CHAPTER 6 Opimaliy Condiions for Unconsrained Problems 1 Unconsrained Opimizaion 11 Exisence Consider he problem of minimizing he funcion f : R n R where f is coninuous on all of R n : P min f(x) x

More information

Monochromatic Infinite Sumsets

Monochromatic Infinite Sumsets Monochromaic Infinie Sumses Imre Leader Paul A. Russell July 25, 2017 Absrac WeshowhahereisaraionalvecorspaceV suchha,whenever V is finiely coloured, here is an infinie se X whose sumse X+X is monochromaic.

More information

Math 315: Linear Algebra Solutions to Assignment 6

Math 315: Linear Algebra Solutions to Assignment 6 Mah 35: Linear Algebra s o Assignmen 6 # Which of he following ses of vecors are bases for R 2? {2,, 3, }, {4,, 7, 8}, {,,, 3}, {3, 9, 4, 2}. Explain your answer. To generae he whole R 2, wo linearly independen

More information

A NOTE ON THE STRUCTURE OF BILATTICES. A. Avron. School of Mathematical Sciences. Sackler Faculty of Exact Sciences. Tel Aviv University

A NOTE ON THE STRUCTURE OF BILATTICES. A. Avron. School of Mathematical Sciences. Sackler Faculty of Exact Sciences. Tel Aviv University A NOTE ON THE STRUCTURE OF BILATTICES A. Avron School of Mahemaical Sciences Sacler Faculy of Exac Sciences Tel Aviv Universiy Tel Aviv 69978, Israel The noion of a bilaice was rs inroduced by Ginsburg

More information

Brock University Physics 1P21/1P91 Fall 2013 Dr. D Agostino. Solutions for Tutorial 3: Chapter 2, Motion in One Dimension

Brock University Physics 1P21/1P91 Fall 2013 Dr. D Agostino. Solutions for Tutorial 3: Chapter 2, Motion in One Dimension Brock Uniersiy Physics 1P21/1P91 Fall 2013 Dr. D Agosino Soluions for Tuorial 3: Chaper 2, Moion in One Dimension The goals of his uorial are: undersand posiion-ime graphs, elociy-ime graphs, and heir

More information

Chapter Floating Point Representation

Chapter Floating Point Representation Chaper 01.05 Floaing Poin Represenaion Afer reading his chaper, you should be able o: 1. conver a base- number o a binary floaing poin represenaion,. conver a binary floaing poin number o is equivalen

More information

Intermediate Macro In-Class Problems

Intermediate Macro In-Class Problems Inermediae Macro In-Class Problems Exploring Romer Model June 14, 016 Today we will explore he mechanisms of he simply Romer model by exploring how economies described by his model would reac o exogenous

More information

18 Biological models with discrete time

18 Biological models with discrete time 8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so

More information

Then. 1 The eigenvalues of A are inside R = n i=1 R i. 2 Union of any k circles not intersecting the other (n k)

Then. 1 The eigenvalues of A are inside R = n i=1 R i. 2 Union of any k circles not intersecting the other (n k) Ger sgorin Circle Chaper 9 Approimaing Eigenvalues Per-Olof Persson persson@berkeley.edu Deparmen of Mahemaics Universiy of California, Berkeley Mah 128B Numerical Analysis (Ger sgorin Circle) Le A be

More information

Expert Advice for Amateurs

Expert Advice for Amateurs Exper Advice for Amaeurs Ernes K. Lai Online Appendix - Exisence of Equilibria The analysis in his secion is performed under more general payoff funcions. Wihou aking an explici form, he payoffs of he

More information

Displacement ( x) x x x

Displacement ( x) x x x Kinemaics Kinemaics is he branch of mechanics ha describes he moion of objecs wihou necessarily discussing wha causes he moion. 1-Dimensional Kinemaics (or 1- Dimensional moion) refers o moion in a sraigh

More information

Seminar 4: Hotelling 2

Seminar 4: Hotelling 2 Seminar 4: Hoelling 2 November 3, 211 1 Exercise Par 1 Iso-elasic demand A non renewable resource of a known sock S can be exraced a zero cos. Demand for he resource is of he form: D(p ) = p ε ε > A a

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

Fishing limits and the Logistic Equation. 1

Fishing limits and the Logistic Equation. 1 Fishing limis and he Logisic Equaion. 1 1. The Logisic Equaion. The logisic equaion is an equaion governing populaion growh for populaions in an environmen wih a limied amoun of resources (for insance,

More information

THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX

THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX J Korean Mah Soc 45 008, No, pp 479 49 THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX Gwang-yeon Lee and Seong-Hoon Cho Reprined from he Journal of he

More information

4.1 - Logarithms and Their Properties

4.1 - Logarithms and Their Properties Chaper 4 Logarihmic Funcions 4.1 - Logarihms and Their Properies Wha is a Logarihm? We define he common logarihm funcion, simply he log funcion, wrien log 10 x log x, as follows: If x is a posiive number,

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

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

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page Assignmen 1 MATH 2270 SOLUTION Please wrie ou complee soluions for each of he following 6 problems (one more will sill be added). You may, of course, consul wih your classmaes, he exbook or oher resources,

More information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

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

Phys1112: DC and RC circuits

Phys1112: DC and RC circuits Name: Group Members: Dae: TA s Name: Phys1112: DC and RC circuis Objecives: 1. To undersand curren and volage characerisics of a DC RC discharging circui. 2. To undersand he effec of he RC ime consan.

More information

3.1 More on model selection

3.1 More on model selection 3. More on Model selecion 3. Comparing models AIC, BIC, Adjused R squared. 3. Over Fiing problem. 3.3 Sample spliing. 3. More on model selecion crieria Ofen afer model fiing you are lef wih a handful of

More information

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

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

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

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

Physics 101 Fall 2006: Exam #1- PROBLEM #1

Physics 101 Fall 2006: Exam #1- PROBLEM #1 Physics 101 Fall 2006: Exam #1- PROBLEM #1 1. Problem 1. (+20 ps) (a) (+10 ps) i. +5 ps graph for x of he rain vs. ime. The graph needs o be parabolic and concave upward. ii. +3 ps graph for x of he person

More information

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,

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

Electrical and current self-induction

Electrical and current self-induction Elecrical and curren self-inducion F. F. Mende hp://fmnauka.narod.ru/works.hml mende_fedor@mail.ru Absrac The aricle considers he self-inducance of reacive elemens. Elecrical self-inducion To he laws of

More information

Random Walk with Anti-Correlated Steps

Random Walk with Anti-Correlated Steps Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and

More information

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8) I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression

More information

Physics for Scientists & Engineers 2

Physics for Scientists & Engineers 2 Direc Curren Physics for Scieniss & Engineers 2 Spring Semeser 2005 Lecure 16 This week we will sudy charges in moion Elecric charge moving from one region o anoher is called elecric curren Curren is all

More information

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion

More information

Section 5: Chain Rule

Section 5: Chain Rule Chaper The Derivaive Applie Calculus 11 Secion 5: Chain Rule There is one more ype of complicae funcion ha we will wan o know how o iffereniae: composiion. The Chain Rule will le us fin he erivaive of

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

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n

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

Section 1.2 Angles and Angle Measure

Section 1.2 Angles and Angle Measure Sec.. ngles and ngle Measure LSSIFITION OF NGLES Secion. ngles and ngle Measure. Righ angles are angles which. Sraigh angles are angles which measure measure 90. 80. Every line forms a sraigh angle. 90

More information

Introduction to Probability and Statistics Slides 4 Chapter 4

Introduction to Probability and Statistics Slides 4 Chapter 4 Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random

More information

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits EEE25 ircui Analysis I Se 4: apaciors, Inducors, and Firs-Order inear ircuis Shahriar Mirabbasi Deparmen of Elecrical and ompuer Engineering Universiy of Briish olumbia shahriar@ece.ubc.ca Overview Passive

More information