EIGENVALUE PROBLEM IN PRECISE MAGNET ALIGNMENT

Size: px
Start display at page:

Download "EIGENVALUE PROBLEM IN PRECISE MAGNET ALIGNMENT"

Transcription

1 EIGENVALUE PROBLEM IN PRECISE MAGNET ALIGNMENT K. Endo and K. Mshma KEK and SKD (Graduate Unversty for Advanced Studes), Tsukuba, Japan Abstract For the algnment of the accelerator components at the desgned coordnates, a matrx connectng observables and coordnates can be treated as an egenvalue problem and has ts own characterstcs dependng on the structure of the survey network. Once ts egenvectors are determned, the msalgnment of the magnets s expressed by ther lnear combnaton. By decomposng the magnet msalgnments nto the egenmodes and evaluatng each contrbuton of the egenmode to the closed orbt dstorton (COD), t s possble to fnd whch mode affects sgnfcantly to the beam n the synchrotron. Ths approach s appled to the synchrotron rngs wth the typcal symmetrc and asymmetrc magnet arrangements to whch several survey networks are adopted and gves the useful nformaton for the magnet precse algnment, to say, the msalgnment modes whch are less sgnfcant to the beam can be safely omtted. In ths work, three knds of survey networks, whch were appled for the real synchrotron, are treated n detals. 1. INTRODUCTION An algnment matrx equaton s obtaned from the trangulaton of the survey networks for the precse algnment of the synchrotron magnets. The horzontal msalgnments, the spatal devatons of magnet from the planned postons n the horzontal plane, are expressed by (x, y) or (r, θ) coordnates wth an orgn at the center of the synchrotron rng. From the survey of the physcal dstances such as the dstance of the adjacent quadrupoles and the perpendcular from the quadrupole to the straght lne stretched between the quadrupoles at both sdes, the devatons of these dstances from the desgn values are expressed approxmately as a lnear functon of the msalgnments of the related quadrupoles. These relatons for all quadrupoles are merged nto a global matrx equaton termed here as an algnment equaton. The matrx depends on both the magnet lattce desgn and the survey network and connects the magnet msalgnment wth the devatons of the surveyed dstances, or the survey varables. As the magnet lattce s fxed, the varous knds of the survey network formed on the targets of quadrupoles or the representatve ponts close to each quadrupole ntroduce the varance n the matrx of the algnment equaton. The survey network or mesh should be determned accordng to the magnet lattce structure and the easness of the survey measurement n the narrow synchrotron tunnel. If the target s placed on each quadrupole magnet, an accurate drlled hole where the target s nserted wth a mnmal clearance s requred. The msalgnment s corrected by movng the

2 quadrupole so that the target sts at the desgned coordnate. If the representatve ponts are used, each quadrupole s algned by referrng these ponts. Large colldng beam synchrotron wth the expermental nsertons s composed of the arced and straght sectons. The arced secton has a regular magnet lattce and quadrupoles are placed equdstant, but the expermental straght secton has dfferent quadrupole arrangements because t accommodates the RF cavtes and mn-beta quadrupole magnet system to squeeze the beam at the collson pont. The synchrotron rng n whch the quadrupoles are not placed on the common crcle s classfed here as an asymmetrc rng (nhomogeneous rng). To the contrary the hypothetcal rng where all quadrupoles st on the common crcle s called here as a symmetrc rng (homogeneous rng).. DERIVATION OF THE ALIGNMENT EQUATION Assumng that ether symmetrc or asymmetrc synchrotron rng has N quadrupoles and ther N survey targets ( = 1,,, N ) are on the desgned orbt, at least two survey varables S and P are necessary to each quadrupole. These varables are expressed as a functon of R and Θ as S = F(R,Θ ) and P = G(R,Θ ), (1) where the subscrpt of R and Θ s consdered for all possble range. For an example, f S s a dstance between the adjacent quadrupole magnets, S = F(R, R +1,Θ ). Dfferentatng the above relatons to the frst order, F F F F s S = R + Θ = r + (θ +1 θ ) () R Θ R Θ and G G G G p P = R + Θ = r + (θ +1 θ ), (3) R Θ R Θ where s and p are the devatons from the desgn value whch are determned from the precse survey, and r and θ the quadrupole dsplacements to the radal and azmuthal drecton, respectvely. Dervng relatons () and (3) for all quadrupoles, p1 a11 a1n c11 c1n r1 : : : : : : p N an1 ann cn1 c NN r N = (4) s1 b11 b1n d11 d1n θ1 : : : : : : s N bn1 bnn dn1 d NN θ N where the N x N matrx s obtaned for the above relatons [1][]. If 3 varables S, P and Q are consdered, the matrx becomes 3N x N and the left column vector has 3N components. In any case the matrx equaton (4) s expressed smply as p = (A)r (5) and solved by the least squares method as

3 h (A T )p = ( A T A)r (6) T T where ( A ) s the transposed matrx of ( A ) and ( AA) s the N x N matrx. The soluton s gven by r = (A T A) 1 h, (7) f the matrx s not sngular. In general, however, the matrx ( A) s sngular. Therefore one must consder to reduce the rank of the matrx by fxng several coordnates, for an example, r 1 =, r = and θ 1 = wthout loosng generalty. The number of coordnates to be fxed depends on the adopted survey network. Another method to solve (6) s to use the generalzed nverse matrx. Even f the matrx s sngular, ts nverse matrx s determned unquely followng the method of the sngular value decomposton. Replacng ( A T A) wth (H) and assumng (H) + the generalzed nverse matrx, the soluton of (6) s smply gven by r = (H) + h. (8) In ths manuscrpt the sngular matrx of (4) s reduced to a lower rank to obtan the non-sngular matrx so as to analyze ts propertes. From the relaton (6) the precse magnet algnment problem can be treated as an egenvalue problem, (H)v = λ v, (9) where λ s an egenvalue belongng to the egenvector v. If all components of egenvector are multpled by the coeffcents and are combned lnearly, the dsplacements of the quadrupole magnets ( r and θ ) are obtaned as a lnear combnaton of the egenvectors, r = c k v k, (1) k where c k s the coeffcent. From (6) and (1), (H)c k v k = h. (11) k The lnear coeffcent of (1) s c k = (v T k h)/ λ k, (1) T where v k s the transpose of the k-th egenvector v k. The j-th component of the dsplacement vector s (v r j = kt h) (v v kj or r j = kt A T p) v kj, (13) k λ k where v kj s the j-th component of the egenvector v k [3]. 3 SURVEY NETWORKS k λ k In many cases of the synchrotron rng, all magnets are algned n the narrow tunnel and the survey network for the measurements of dstances between magnets s restrcted. Nowadays very precse nstruments such as a laser tracker and a mekometer are avalable, but the maxmum length to be measured s lmted to mantan accuracy. Thus the survey network must be establshed consderng the dstrbuton of the survey ponts whch can attan the requred accuracy. From (13), the dsplacement error r due to the survey error p s estmated by

4 a pv vv m km m, kj km rj = vkj = am k λk m k λk Then the standard devaton of m kj km = k, m λk r j s expressed as follows, p. (14) a v v σ r σ j p, (15) where σ p s the expected error for each p, the statstcal error for repeated dstance measurements between the same ponts. Between the standard devatons of the coordnates and dstance measurement errors for the survey network there s a relaton medated by the network structure [4]. The smaller coordnate error s expected f omttng the contrbutons by egenmodes wth small egenvalues. For large synchrotrons there are not so many varatons for the survey networks, because the number of the collmaton axes extended from each pont to another s lmted. The typcal survey networks consdered here for the synchrotron are (Case #1) Short chord and perpendcular, (Case #) Long chord and perpendcular, and (Case #3) Short chord and perpendculars, as shown n Fg.1, where the survey varables are gven by the sold lnes. Dfferences of the numbers of egenmodes and vector elements are compared n Table R-1 θ -1 P R S θ θ +1 Θ -1 Θ +1 R+1-1 R-1 θ -1 S +1-1 P,1 P, P +1 - L -1 + R R R-1 R+1 R+1 R- θ θ +1 R+ θ θ -1 θ +1 θ + Θ -1 Θ θ - Θ - Θ -1 Θ Θ +1 Case#1 Case# Case#3 Fg.1 Survey varables (sold lnes) for the dfferent survey networks. The matrx equatons for 3 cases are derved for both symmetrc and asymmetrc rng, assumng both rngs have the same dameter (48 m) and the same number of the quadrupole magnets (39 quads). The asymmetrc rng s an electron-postron colldng synchrotron called TRISTAN of KEK whch has the long expermental nsertons at both sdes of 4 nteracton ponts. To the contrary, the symmetrc rng s the hypothetcal synchrotron of whch all quadrupole magnets locate on the same crcle at equdstant all around the rng.

5 The egenvectors for largest and smallest egenvalues are shown for the symmetrc 3 cases n Fg. and for the asymmetrc 3 cases n Fg.3. In these fgures the radal and azmuthal dsplacements are gven n seperate both n the order of quadrupole magnet confguraton. Table 1 Dfferences of the numbers of egenmodes and vector elements for 3 survey networks. N s the number of the quadrupole magnets. Case #1 Case # Case #3 Number of egenmodes N-3 N-4 N-3 Number of the egenvector elements of each egenmode N-3 N-4 N # Symmetrc Case# # # Symmetrc Case# # # Symmetrc Case#1.1 # # Symmetrc Case# # Egenvector #779 Egenvector #78.1 # Symmetrc Case# # # Symmetrc Case# # Egenvector #781 Egenvector #78 Egenvector #78 Egenvector # # Symmetrc Case# # # Symmetrc Case# # Egenvector #78 Egenvector #781.1 # Symmetrc Case# # # Symmetrc Case# # Egenvector #78.1 # Symmetrc Case# # Egenvector #781 Egenvector # # Symmetrc Case# # Egenvector #781 Fg. Egenvectors wth largest and smallest egenvalues for the symmetrc cases. has the largest egenvalue and egenvector #781 (#78 for case#) the smallest one. These egenvectors gve the magnet msalgnment patterns to the radal and azmuthal drectons n the order of magnet confguraton. Egenvectors for (a) Upper: Symmetrc case#1, (b) Mddle: Symmetrc case#, and (c) Lower: Symmetrc case#3. In every case the same egenvector s obtaned for the smallest egenvalue. It gves the quadrupole dsplacements formng a snusodal undulaton along the synchrotron rng n both radal and azmuthal drectons. Generally for larger egenvalue the number of the snusodal dsplacement

6 undulatons along the rng ncreases accordng to the magntude of the egenvalue. They gve the characterstc envelope patterns of the snusodal dsplacements for large egenvalues. However, these msalgnment undulatons are obtaned wthout consderng the betatron moton of the partcle. As each egenvector has ts own pattern of the magnet msalgnment, f the betatron oscllaton number becomes smlar to the quadrupole magnet algnment pattern gves the relatve large COD to the beam. Asymmetrc Case#1. # # Asymmetrc Case#. # # Asymmetrc Case# # # Asymmetrc Case#1..6 # # Asymmetrc Case#..6 # # Asymmetrc Case# # # Egenvector #779 Egenvector #78 Egenvector #78.1 # Asymmetrc Case# # Asymmetrc Case# # # # Asymmetrc Case# # Egenvector #781 Egenvector #781 Egenvector #78 Egenvector #78 Egenvector #779 Egenvector # # Asymmetrc Case# # Asymmetrc Case# # # # Asymmetrc Case# # Egenvector #78 Egenvector #781 Egenvector #781 Fg.3 Egenvectors wth largest and smallest egenvalues for the asymmetrc cases. has the largest egenvalue and egenvector #781 (#78 for the case#) the smallest one. These egenvectors gve the magnet msalgnment patterns to the radal and azmuthal drectons n the order of magnet confguraton. Egenvectors for (a) Upper: Asymmetrc case#1, (b) Mddle: Asymmetrc case#, and (c) Lower: Asymmetrc case#3. 4. RMS CLOSED ORBIT DISTORSION (COD) To estmate the rms COD the magntude of the msalgnment s normalzed to N ( x ) =.1mm, (16) N =1

7 after gvng the radal and azmuthal msalgnments to every quadrupole magnet. If N =1 ( x ) / N =.1mm s assumed only for the radal dsplacements, an extraordnary azmuthal dsplacement s ntroduced even for the small radal msalgnment. To avod ths problem, the normalzaton over both radal and azmuthal dsplacements s adopted. For the calculaton of COD the beam optcs code s used [5]. An essental feature of the COD calculaton depends on the transfer matrces of all machne elements ncludng quadrupoles, dpoles and drft spaces that the crculatng beam encounters when travelng the synchrotron orbt. The transfer matrx through one perod can be rewrtten by usng the twss parameters and the betatron phase advances. Once the twss parameters are obtaned at a pont, those at other ponts are easly calculated. A matrx to calculate COD at the beam poston montors due to the msalgnments of the quadrupole magnets can be derved by usng the twss parameters. If each egenvector s appled to ths matrx, ts contrbuton to the beam s estmated by usng the rms value of COD. Every quadrupole s shfted to the radal and longtudnal drectons accordng to the assumed msalgnment n the horzontal plane and the drft space lengths are adjusted to make up for the longtudnal shft of quadrupoles before the COD calculaton. By nvestgatng the relatons between the msalgnment modes and the beam orbt deformatons, the underlyng rules for the accelerator algnment ssue wll be found. In the synchrotron the charged partcles perform the betatron oscllaton whch s determned by the magnetc feld gradent and the confguraton of the quadrupole magnets. The COD s defned as the orbt of the gravtatonal center of the crculatng partcles. If the undulaton number of the msalgnment pattern concdes wth that of the betatron oscllaton, partcles suffer the Lorenz force due to the systematcally varyng magnetc feld of the dsplaced quadrupole magnet and wll result n the large orbt dstorton The rms COD s for each egenmode The rms COD s defned here as N k 1 < r >= ( k ) N = 1 r, (17) where r k s the radal COD at the -th beam poston montor and < r k > the correspondng rms COD orgnated from the k-th egenvector. The rms COD may be used to estmate the effect to the beam orbt by each egenmode Symmetrc cases In general, the rms COD responses are dfferent dependng on the survey networks. However, there s no structural dfference between cases #1 and # for egenvalues as shown n Fg.4. The undulaton number of the egenvector (egenmode #687) shown n Fg.5 s very close to the betatron oscllaton frequency. The egenmode #687 of the symmetrc case# has the same undulaton pattern as Fg.5. The largest peak correspondng to the egenmode #178 has the egenvector shown n Fg.6. Ths egenmode has about 4 tmes of the undulaton number of the egenmode #687 whch s close to the horzontal betatron frequency. As almost all dots shown n Fg.7 ft on the phase of the

8 msalgnment pattern of the egenmode #178 whch has the largest rms COD response, the egenmode #178 s expected to have a large contrbuton to COD. The same relaton concernng to the largest rms COD peak holds on other cases. The egenmode #176 of the symmetrc case# have the same undulaton pattern except for the azmuthal part whch has 71 undulatons. rmscod (arbtrary) Egenmode #178 Symmetrc Case#1 Tune hor/ver = 37./39.1 rmscod Egenmode number Egenmode #687 & # rms COD (arbtrary) Egenmode #176 Symmetrc Case# Tune hor/ver = 37./39.1 rms COD Egenmode number Egenmode #687 & # (a) Symmetrc case#1 (b) Symmetrc case# Fg.4 The rms COD responses of the symmetrc cases #1 and #. Egenvector # oscllatons Symmetrc Case#1 37 oscllatons # Azmuthal poston Radal & s Fg.5 Undulaton pattern of the egenvector #687 ( nd peak) of the symmetrc case#1. The survey network of the symmetrc case #3 has the 3 rd peak as shown n Fg.8 whch does not exst n cases #1 and #. The nd peak corresponds to the egenmode #65 havng the undulaton number close to the betatron frequency (Fg.9). But the 1 st and 3 rd peaks have larger undulaton numbers as shown n Fg.1 and Fg.11, respectvely.

9 oscllatons Symmetrc Case#1 149 oscllatons # Radal poston Radal & s Fg.6 Undulaton pattern of the egenvector #178 (1 st peak) of the symmetrc case#1.1 Egenmode #687 (red and blue dots) s superposed on the egenvector #178 (sold lne) Symmetrc case# A part of radal poston Fg.7 Egenmode #178 (sold lne) has 4 tmes of the undulaton number of the egenmode #687 (red and blue dots) whch s close to the horzontal betatron frequency n the symmetrc case#1. Almost all dots ft on the phase of the msalgnment pattern of the egenmode #178 whch has the largest rms COD response. Only one fourth of the crcumference s shown for the explanaton. #178 rms COD (arbtrary) Egenmode #18 rms COD Egenmode #468 Symmetrc Case#3 Tune hor/ver = 37./39.1 Egenmode # Egenmode number correspondng to egenvalue Fg.8 The rms COD response of the symmetrc case #3.

10 Egenmode #65.1 Egenmode of the nd peak of the rms COD.5 Symmetrc Case# oscllatons #65 36 oscllatons Radal & s Fg.9 Undulaton pattern of the egenvector #65 ( nd peak) of the symmetrc case# Egenmode #18 for the largest peak of the rms COD Symmetrc Case#3 # oscllatons Radal & s Fg.1 Undulaton pattern of the egenvector #18 (1 st peak) of the symmetrc case#3. Egenmode # Egenmode #468 for the 3rd peak of the rms COD #468 Symmetrc Case#3 89 oscllatons 66 oscllatons Radal poston Radal & s Fg.11 Undulaton pattern of the egenvector #468 (3 rd peak) of the symmetrc case#3.

11 4.1. Asymmetrc cases In the asymmetrc cases there appear many peaks n the rms COD response spectra. It seems that the rregularty n dstance between quadrupole magnets plays a role. As shown n Fg.1 the cases #1 and # have almost the same rms COD responses and the peaks around the egenmode #677 correspond to the egenvector close to the betatron frequency as shown n Fg.13. The egenmode for the largest peak (Fg.14) has a smlar egenvector as Fg.6. 5 Egenmode # Egenmode #174 5 rmscod (arbtrary) 4 3 Asymmetrc Case#1 Tune hor/ver = 37./39.1 Egenmode #59 Egenmode # rmscod (arbtrary) 4 3 Asymmetrc Case# Tune hor/ver = 37./39.1 Egenmode # rmscod 1 1 rmscod Egenmode number Egenmode number (a) Asymmetrc case#1 (b) Asymmetrc case# Egenvector #677 Fg.1 The rms COD responses of the asymmetrc cases #1 and #..1 4 oscllatons Asymmetrc Case# oscllatons -.5 # Radal & s Fg.13 Undulaton pattern of the egenvector #677 of the asymmetrc case#1.

12 oscllatons.5 Asymmetrc Case#1 148 oscllatons -.5 # Azmutrhal poston Radal & s Fg.14 Undulaton pattern of the egenvector #176 of the asymmetrc case#1. rmscod (arbtrary) Egenmode #19 #149 rmscod Asymmetrc Case#3 Tune hor/ver = 37./39.1 #454 #63 #67 # Egenmode number Fg.15 The rms COD response of the asymmetrc case #3. Egenvector # oscllatons Asymmetrc Case#3 4 oscllatons # Radal & s Fg.16 Undulaton pattern of the egenvector #69 of the asymmetrc case#3.

13 For the asymmetrc case#3 has the rms COD response spectrum dfferent from Fg.1 as shown n Fg.15. The peaks of the rms COD responses for the egenmodes correspondng to the betatron frequency shft to larger egenvalues. The undulaton number of the egenmode #69 (Fg.16) s less than the betatron frequency whch corresponds to the egenmode #63 (Fg.17). Comparng the cases #1, # and #3, there are bg dfferences as shown n Fg.1 and Fg.15. However, t seems lttle effect to COD by omttng the correctons of egenmodes (more than #7) wth small egenvalues. In ths way the beam senstve egenmodes, whch must be corrected, are selected through the nvestgaton of the egenmode problem. Egenvector # oscllatons Asymmetrc Case#3 4 oscllatons # Radal & s Fg.17 Undulaton pattern of the egenvector #63 of the asymmetrc case#3. 4. Solvng the egenmode based algnment problem If several egenmodes are selected as the quadrupole algnment errors, the correspondng egenvectors can be combned lnearly assumng arbtrary coeffcents and transformed to the survey varables such as the perpendcular and short chord lengths n the case#1. Accordng to (13) the dfference of the survey data from the desgn values can be processed to obtan the quadrupole msalgnments. An example for the asymmetrc case#1 s gven n Fg.18, where the 677-th egenmode s assumed as the msalgnment mode and the dfference between the msalgnment and the correspondng egenvector s shown by the blue lne. Most sgnfcant contrbuton to the msalgnment s the 677-th egenmode as shown n Fg.18 by the peak of the green lne. Other egenmodes at the rght sde of the absolute coeffcent (green lne) also contrbute the msalgnment but these egenmodes arse by solvng the sngular matrx problem by the rank reducton method. The survey varables, perpendcular and short chord length errors, are obtaned from the 677-th egenvector as gven n Fg.19. Even f two egenmodes are lnearly combned, the same procedure can be appled to obtan the msalgnment. Combnng the 176-th and the 677-th egenmodes wth the equal weght but the dfferent sgns for the asymmetrc case#1 (Fg.), both egenmodes have sgnfcant contrbutons to the msalgnment and egenmodes at the rght sde of the absolute coeffcent (green lne) also have relatvely large contrbutons whch are also ntroduced through the rank reducton of the sngular matrx. The dfference between the calculated msalgnment and the combned egenvector

14 arses from the exctaton of the lower undulaton frequency modes wth small egenvalues. The egenmodes wth the undulaton frequency much lower than the betatron frequency and wth the small rms COD response do not contrbute to the beam moton and can be safely dsregarded from the msalgnment correcton. If the random msalgnment s assumed to obtan the perpendcular and short chord errors, most egenmodes contrbute equally to the msalgnment as shown n Fg.1 for the asymmetrc case#1. Soluton & Dfference (=Sol-677th egenvector) Asymmetrc Case#1, 677-th Egenvector Sol Sol-677thEV Abs.coeffcent Radal and azmuthal postons of quadrupole magnets or Egenmode number for Absolute coeffcents Fg.18 Msalgnment due to the 677-th egenmode of the asymmetrc case#1. The red lne s the quadrupole msalgnment and the blue lne s the dfference between the calculated msalgnment and the 677-th egenvector. The green lne s the absolute coeffcents combnng all egenvectors to make up the msalgnment of the red lne Abs.coeffcent 677-th Egenvector thEV Asymmetrc Case#1, 677-th Egenvector Perpendcular&ShortChord Perpendcular error Short chord error Radal and azmuthal postons of quadrupole magnets Perpendcular & Short Chord Errors Fg.19 Comparson of the 677-th egenvector and the correspondng survey varables obtaned from ths egenvector for the asymmetrc case#1. The perpendcular and short chord errors (blue lne) gve the soluton (red lne) n Fg.18.

15 Soluton & Dfference(=Sol-176th-677thEV) Asymmetrc Case#1, 176th + 677th Egenvector Soluton Sol-176thEV-677thEV Abs.coeffcent Radal and azmuthal postons of quadrupole magnets or Egenmode number for Absolute coeffcents Abs.coeffcent Fg. Msalgnment due to the combnaton of the 176-th and the 677-th egenmodes. Gven random msalgnment (m) Perpendcular & Shortchord eroor Contrbuton of each egenmode Gven random msalgnment.1 mm (rms) Asymmetrc case#1 Random msalgnment Egenmode Asymmetrc case#1 Soluton - Gven random msalgnment -.6 Random msalgnment Quadrupole magnet.1..5 Perpendcular & Shortchord eroor Asymmetrc case#1 Ramdom msalgnment.15 Perpendcular error Short chord error Soluton Radal msalgnment Azmuthal msalgnment Quadrupole magnet Sol-Gven random msalgnment (m) Soluton (a) Assumed random msalgnment (d) Dfference between soluton and assumed random msalgnment: (d) = (c) - (a) (b) Perpendcular and Shott chord error obtaned from the assumed random msalgnment error (c) Solved by egenvalue method usng errors of (b) (e) Contrbuton of egenmode to the soluton of (c) Random dsplacement error for Asymmetrc case#1 Fg.1 Random msalgnment of.1 mm (rms) s assumed for the asymmetrc case#1. Most egenmodes contrbute equally to the msalgnment and t also gves the smlar dfference between the assumed and calculated msalgnment.

16 4.3 Effect of the betatron tunes to rms COD The peaks of the rms COD response shown n Fg.4, Fg.8, Fg.1 and Fg.15 shft, f the horzontal betatron frequency changes. Examples for the asymmetrc and symmetrc case#1 are gven n Fg.. The sharp rms COD response peaks are obtaned for the symmetrc cases whch are better to get the relaton between the betatron tune and egenmode. If the horzontal tune ncreases, the maxmum rms COD response peaks shft towards larger egenmode number whch has smaller egenvalue and the rms COD response peaks n the symmetrc case#1 close to the horzontal betatron tune shft to the smaller egenmode number whch has larger undulaton frequency. 8 Symmetrc case#1 8 Asymmetrc case#1 7 7 Egenmode to rmscod response Peak rmscod response Egenmode close to QH frequency Egenmode to rmscod response Egenmode of max. rmscod response Egenmode close to QH freq Horzontal tune (QH) (a) Symmetrc case# Horzontal tune (QH) (b) Asymmetrc case#1 Fg. Dependence of the rms COD response peaks on the horzontal betatron tune for the symmetrc and asymmetrc case#1. 5. CONCLUSION As the egenvectors of the matrx whch s derved for the survey network correspond to the components of the msalgnment, t s possble to estmate ther effects on COD by treatng each egenvector as a msalgnment mode n the beam optcs code. In ths study the followng results were obtaned. 1) There are some peaks n the rms COD response spectrum correspondng to the egenmodes of whch snusodal undulaton numbers are close to or multple of the horzontal betatron oscllaton frequency. These peaks shft dependng on the horzontal tune.

17 ) For the nhomogeneous quadrupole confguraton, there are many rms COD response peaks, whereas for the homogeneous one only few peaks appear. 3) In practce, the msalgnment modes wth relatvely large rms COD response should be corrected. Consderng the possble extent of the betatron tune varaton, the msalgnment correcton should be made for all egenmodes that the large rms COD response s expected. 4) Lower modes, havng small egenvalue whch have small number of the snusodal undulatons, are nevtably ntroduced through the reducton of the rank when obtanng non-sngular matrx, but these modes can be left uncorrected because ther effects on the rms COD response are very small. 5) The largest peak of the rms COD response corresponds to the undulaton number close to 4 tmes of the horzontal betatron frequency. Ths msalgnment pattern wll ft well to the msalgnment mode of whch undulaton number s close to the horzontal betatron frequency. 6. REFERENCES [1] K. Endo and K. Mshma, Survey Meshes for Synchrotron Magnet Algnment, Proc. EPAC, Venna, pp [] K. Mshma and K. Endo, COD Response to Msalgnment Mode of Synchrotron Quadrupoles, APAC1, Bejng, pp [3] J.H. Wlknson, The Algebrac Problem, Oxford Unv. Press, G. Goertzel and N. Trall, Some Mathematcal Methods of Physcs, McGraw-Hll Book Co. Inc., 196. [4] K. Mshma, Study on the Precse Algnment of Accelerator (n Japanese), thess for dssertaton,. [5] A.S. Kng, M.J. Lee and W.W. Lee, MAGIC, A Computer Code for Desgn Studes of Insertons and Storage Rng, SLAC-183, It was modfed greatly to apply for the COD calculaton due to the quad dsplacements.

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

Formulas for the Determinant

Formulas for the Determinant page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use

More information

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system. Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and

More information

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is. Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion Assessment of Ste Amplfcaton Effect from Input Energy Spectra of Strong Ground Moton M.S. Gong & L.L Xe Key Laboratory of Earthquake Engneerng and Engneerng Vbraton,Insttute of Engneerng Mechancs, CEA,

More information

Statistics Chapter 4

Statistics Chapter 4 Statstcs Chapter 4 "There are three knds of les: les, damned les, and statstcs." Benjamn Dsrael, 1895 (Brtsh statesman) Gaussan Dstrbuton, 4-1 If a measurement s repeated many tmes a statstcal treatment

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

COMPLEX NUMBERS AND QUADRATIC EQUATIONS COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

Mathematical Preparations

Mathematical Preparations 1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the

More information

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2 Physcs 607 Exam 1 Please be well-organzed, and show all sgnfcant steps clearly n all problems. You are graded on your wor, so please do not just wrte down answers wth no explanaton! Do all your wor on

More information

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.

More information

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department

More information

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown

More information

Negative Binomial Regression

Negative Binomial Regression STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent

More information

Norms, Condition Numbers, Eigenvalues and Eigenvectors

Norms, Condition Numbers, Eigenvalues and Eigenvectors Norms, Condton Numbers, Egenvalues and Egenvectors 1 Norms A norm s a measure of the sze of a matrx or a vector For vectors the common norms are: N a 2 = ( x 2 1/2 the Eucldean Norm (1a b 1 = =1 N x (1b

More information

x = , so that calculated

x = , so that calculated Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to

More information

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

This model contains two bonds per unit cell (one along the x-direction and the other along y). So we can rewrite the Hamiltonian as:

This model contains two bonds per unit cell (one along the x-direction and the other along y). So we can rewrite the Hamiltonian as: 1 Problem set #1 1.1. A one-band model on a square lattce Fg. 1 Consder a square lattce wth only nearest-neghbor hoppngs (as shown n the fgure above): H t, j a a j (1.1) where,j stands for nearest neghbors

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

MEASUREMENT OF MOMENT OF INERTIA

MEASUREMENT OF MOMENT OF INERTIA 1. measurement MESUREMENT OF MOMENT OF INERTI The am of ths measurement s to determne the moment of nerta of the rotor of an electrc motor. 1. General relatons Rotatng moton and moment of nerta Let us

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

Kinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017

Kinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017 17/0/017 Lecture 16 (Refer the text boo CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlnes) Knematcs of Fluds Last class, we started dscussng about the nematcs of fluds. Recall the Lagrangan and Euleran

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Analyss of Varance and Desgn of Exerments-I MODULE III LECTURE - 2 EXPERIMENTAL DESIGN MODELS Dr. Shalabh Deartment of Mathematcs and Statstcs Indan Insttute of Technology Kanur 2 We consder the models

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Chapter 13: Multiple Regression

Chapter 13: Multiple Regression Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to

More information

CHALMERS, GÖTEBORGS UNIVERSITET. SOLUTIONS to RE-EXAM for ARTIFICIAL NEURAL NETWORKS. COURSE CODES: FFR 135, FIM 720 GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET. SOLUTIONS to RE-EXAM for ARTIFICIAL NEURAL NETWORKS. COURSE CODES: FFR 135, FIM 720 GU, PhD CHALMERS, GÖTEBORGS UNIVERSITET SOLUTIONS to RE-EXAM for ARTIFICIAL NEURAL NETWORKS COURSE CODES: FFR 35, FIM 72 GU, PhD Tme: Place: Teachers: Allowed materal: Not allowed: January 2, 28, at 8 3 2 3 SB

More information

A Robust Method for Calculating the Correlation Coefficient

A Robust Method for Calculating the Correlation Coefficient A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal

More information

Problem Points Score Total 100

Problem Points Score Total 100 Physcs 450 Solutons of Sample Exam I Problem Ponts Score 1 8 15 3 17 4 0 5 0 Total 100 All wor must be shown n order to receve full credt. Wor must be legble and comprehensble wth answers clearly ndcated.

More information

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

11. Dynamics in Rotating Frames of Reference

11. Dynamics in Rotating Frames of Reference Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons

More information

A how to guide to second quantization method.

A how to guide to second quantization method. Phys. 67 (Graduate Quantum Mechancs Sprng 2009 Prof. Pu K. Lam. Verson 3 (4/3/2009 A how to gude to second quantzaton method. -> Second quantzaton s a mathematcal notaton desgned to handle dentcal partcle

More information

Chapter 3 Describing Data Using Numerical Measures

Chapter 3 Describing Data Using Numerical Measures Chapter 3 Student Lecture Notes 3-1 Chapter 3 Descrbng Data Usng Numercal Measures Fall 2006 Fundamentals of Busness Statstcs 1 Chapter Goals To establsh the usefulness of summary measures of data. The

More information

Report on Image warping

Report on Image warping Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

8.592J: Solutions for Assignment 7 Spring 2005

8.592J: Solutions for Assignment 7 Spring 2005 8.59J: Solutons for Assgnment 7 Sprng 5 Problem 1 (a) A flament of length l can be created by addton of a monomer to one of length l 1 (at rate a) or removal of a monomer from a flament of length l + 1

More information

Feb 14: Spatial analysis of data fields

Feb 14: Spatial analysis of data fields Feb 4: Spatal analyss of data felds Mappng rregularly sampled data onto a regular grd Many analyss technques for geophyscal data requre the data be located at regular ntervals n space and/or tme. hs s

More information

π e ax2 dx = x 2 e ax2 dx or x 3 e ax2 dx = 1 x 4 e ax2 dx = 3 π 8a 5/2 (a) We are considering the Maxwell velocity distribution function: 2πτ/m

π e ax2 dx = x 2 e ax2 dx or x 3 e ax2 dx = 1 x 4 e ax2 dx = 3 π 8a 5/2 (a) We are considering the Maxwell velocity distribution function: 2πτ/m Homework Solutons Problem In solvng ths problem, we wll need to calculate some moments of the Gaussan dstrbuton. The brute-force method s to ntegrate by parts but there s a nce trck. The followng ntegrals

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

Georgia Tech PHYS 6124 Mathematical Methods of Physics I Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends

More information

Army Ants Tunneling for Classical Simulations

Army Ants Tunneling for Classical Simulations Electronc Supplementary Materal (ESI) for Chemcal Scence. Ths journal s The Royal Socety of Chemstry 2014 electronc supplementary nformaton (ESI) for Chemcal Scence Army Ants Tunnelng for Classcal Smulatons

More information

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011 Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected

More information

Singular Value Decomposition: Theory and Applications

Singular Value Decomposition: Theory and Applications Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real

More information

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t 8.5: Many-body phenomena n condensed matter and atomc physcs Last moded: September, 003 Lecture. Squeezed States In ths lecture we shall contnue the dscusson of coherent states, focusng on ther propertes

More information

Uncertainty as the Overlap of Alternate Conditional Distributions

Uncertainty as the Overlap of Alternate Conditional Distributions Uncertanty as the Overlap of Alternate Condtonal Dstrbutons Olena Babak and Clayton V. Deutsch Centre for Computatonal Geostatstcs Department of Cvl & Envronmental Engneerng Unversty of Alberta An mportant

More information

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

Chapter 9: Statistical Inference and the Relationship between Two Variables

Chapter 9: Statistical Inference and the Relationship between Two Variables Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,

More information

Title: Radiative transitions and spectral broadening

Title: Radiative transitions and spectral broadening Lecture 6 Ttle: Radatve transtons and spectral broadenng Objectves The spectral lnes emtted by atomc vapors at moderate temperature and pressure show the wavelength spread around the central frequency.

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Comparison of Regression Lines

Comparison of Regression Lines STATGRAPHICS Rev. 9/13/2013 Comparson of Regresson Lnes Summary... 1 Data Input... 3 Analyss Summary... 4 Plot of Ftted Model... 6 Condtonal Sums of Squares... 6 Analyss Optons... 7 Forecasts... 8 Confdence

More information

Number of cases Number of factors Number of covariates Number of levels of factor i. Value of the dependent variable for case k

Number of cases Number of factors Number of covariates Number of levels of factor i. Value of the dependent variable for case k ANOVA Model and Matrx Computatons Notaton The followng notaton s used throughout ths chapter unless otherwse stated: N F CN Y Z j w W Number of cases Number of factors Number of covarates Number of levels

More information

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016 U.C. Berkeley CS94: Spectral Methods and Expanders Handout 8 Luca Trevsan February 7, 06 Lecture 8: Spectral Algorthms Wrap-up In whch we talk about even more generalzatons of Cheeger s nequaltes, and

More information

Linear Correlation. Many research issues are pursued with nonexperimental studies that seek to establish relationships among 2 or more variables

Linear Correlation. Many research issues are pursued with nonexperimental studies that seek to establish relationships among 2 or more variables Lnear Correlaton Many research ssues are pursued wth nonexpermental studes that seek to establsh relatonshps among or more varables E.g., correlates of ntellgence; relaton between SAT and GPA; relaton

More information

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA 4 Analyss of Varance (ANOVA) 5 ANOVA 51 Introducton ANOVA ANOVA s a way to estmate and test the means of multple populatons We wll start wth one-way ANOVA If the populatons ncluded n the study are selected

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

More information

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11)

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11) Gravtatonal Acceleraton: A case of constant acceleraton (approx. hr.) (6/7/11) Introducton The gravtatonal force s one of the fundamental forces of nature. Under the nfluence of ths force all objects havng

More information

LECTURE 9 CANONICAL CORRELATION ANALYSIS

LECTURE 9 CANONICAL CORRELATION ANALYSIS LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of

More information

THE CURRENT BALANCE Physics 258/259

THE CURRENT BALANCE Physics 258/259 DSH 1988, 005 THE CURRENT BALANCE Physcs 58/59 The tme average force between two parallel conductors carryng an alternatng current s measured by balancng ths force aganst the gravtatonal force on a set

More information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

More information

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS) Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998

More information

Module 3: Element Properties Lecture 1: Natural Coordinates

Module 3: Element Properties Lecture 1: Natural Coordinates Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers

More information

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George

More information

Structure from Motion. Forsyth&Ponce: Chap. 12 and 13 Szeliski: Chap. 7

Structure from Motion. Forsyth&Ponce: Chap. 12 and 13 Szeliski: Chap. 7 Structure from Moton Forsyth&once: Chap. 2 and 3 Szelsk: Chap. 7 Introducton to Structure from Moton Forsyth&once: Chap. 2 Szelsk: Chap. 7 Structure from Moton Intro he Reconstructon roblem p 3?? p p 2

More information

The Feynman path integral

The Feynman path integral The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space

More information

In this section is given an overview of the common elasticity models.

In this section is given an overview of the common elasticity models. Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp

More information

ρ some λ THE INVERSE POWER METHOD (or INVERSE ITERATION) , for , or (more usually) to

ρ some λ THE INVERSE POWER METHOD (or INVERSE ITERATION) , for , or (more usually) to THE INVERSE POWER METHOD (or INVERSE ITERATION) -- applcaton of the Power method to A some fxed constant ρ (whch s called a shft), x λ ρ If the egenpars of A are { ( λ, x ) } ( ), or (more usually) to,

More information

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on

More information

Chapter 8. Potential Energy and Conservation of Energy

Chapter 8. Potential Energy and Conservation of Energy Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal

More information

Comparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy

Comparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy Comparatve Studes of Law of Conservaton of Energy and Law Clusters of Conservaton of Generalzed Energy No.3 of Comparatve Physcs Seres Papers Fu Yuhua (CNOOC Research Insttute, E-mal:fuyh1945@sna.com)

More information

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 3, Issue 1, July 2013

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 3, Issue 1, July 2013 ISSN: 2277-375 Constructon of Trend Free Run Orders for Orthogonal rrays Usng Codes bstract: Sometmes when the expermental runs are carred out n a tme order sequence, the response can depend on the run

More information

THEOREMS OF QUANTUM MECHANICS

THEOREMS OF QUANTUM MECHANICS THEOREMS OF QUANTUM MECHANICS In order to develop methods to treat many-electron systems (atoms & molecules), many of the theorems of quantum mechancs are useful. Useful Notaton The matrx element A mn

More information

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Adiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram

Adiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram Adabatc Sorpton of Ammona-Water System and Depctng n p-t-x Dagram J. POSPISIL, Z. SKALA Faculty of Mechancal Engneerng Brno Unversty of Technology Techncka 2, Brno 61669 CZECH REPUBLIC Abstract: - Absorpton

More information

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force. Secton 1. Dynamcs (Newton s Laws of Moton) Two approaches: 1) Gven all the forces actng on a body, predct the subsequent (changes n) moton. 2) Gven the (changes n) moton of a body, nfer what forces act

More information