Compact finite-difference approximations for anisotropic image smoothing and painting
|
|
- Jeffry Cannon
- 5 years ago
- Views:
Transcription
1 CWP-593 Compact finite-difference approximation for aniotropic image moothing and painting Dave Hale Center for Wave Phenomena, Colorado School of Mine, Golden CO 80401, USA ABSTRACT Finite-difference approximation are uccinctly repreented by their tencil, a et of weight that when applied to adjacent ample of a function approximate ome differential operator. In image proceing the ample are pixel or voxel, and the differential operator mut be inverted for moothing or painting application. For efficiency in uch application requiring invere, the finitedifference tencil hould be compact, uing only a mall 3 3 et of nine pixel or a et of twenty-even voxel. From 2 2 and approximation to gradient operator I obtain 3 3 (9-point) and (27-point) tencil that approximate Laplacian operator. The latter may include tenor coefficient that make them aniotropic. By deriving finite-difference tencil for aniotropic Laplacian in thi way, that i, from approximation to gradient, we guarantee that our approximation to aniotropic Laplacian are ymmetric and poitive-emidefinite. And by chooing the gradient approximation carefully, dicretization error can be made iotropic to leading order. Key word: finite-difference numerical method 1 INTRODUCTION A Laplacian i a ymmetric differential operator that roughen. When applied to a function, it remove any contant bia or linear trend while enhancing higher frequencie. Thi behavior i apparent in the Fourier tranform of the Laplacian operator. For if F (k) denote the Fourier tranform of a function f(x) then f(x) F (k), (1) f(x) F (k). (2) The minu ign in make thi ymmetric differential operator poitive-emidefinite (SPS); it Fourier tranform = k T k i non-negative and implie amplification of higher wavenumber. An iotropic Laplacian roughen a function equally in all direction; it Fourier tranform depend on only the magnitude k of the wavenumber vector k. To roughen a function in one direction more than other, we may include a ymmetric poitiveemidefinite (SPS) tenor D in our Laplacian operator: D f(x) k T Dk F (k). (3) For thi aniotropic Laplacian operator, the direction of maximum roughening i aligned with that eigenvector of D having larget eigenvalue. In practical application involving image, the function f(x) i ampled and partial derivative mut be approximated with finite difference. If we arrange all image ample into a ingle column vector f, then an aniotropic Laplacian may be approximated by D f G T DGf, (4) where G i ome pare matrix that repreent a finitedifference approximation to a gradient. In particular, for a 2-D image with N = N 1 N 2 ample, G i a 2N N matrix G1 G =, (5) G 2 with N N pare matrix component G 1 and G 2 that
2 76 D. Hale approximate partial derivative along the firt and econd image dimenion, repectively. In the finite-difference approximation of equation 4, the tenor field D i repreented by a block-diagonal matrix D11 D 12 D =, (6) D 12 D 22 compoed of N N diagonal matrice D 11, D 12 and D 22. The 3N poibly non-zero element of thoe matrice determine the direction and extent of roughening, which may vary from ample to ample. 1.1 Image moothing If an aniotropic Laplacian roughen in certain direction, then the invere of an aniotropic Laplacian mooth in thoe ame direction, and one application of uch an invere operator i tructure-oriented image moothing (e.g., Fehmer and Höcker, 2003). Given an input image f, an aniotropically moothed output image g may be obtained a the olution to: (I + G T DG)g = f. (7) The addition of the identity matrix I to the aniotropic Laplacian G T DG ha two conequence. Firt, it make the ymmetric matrix I + G T DG poitive-definite (SPD) o that a olution g can be found efficiently uing conjugate-gradient iteration. Second, it implie that no moothing occur in eigenvector direction for which correponding eigenvalue of the tenor in D are zero. The direction and extent of moothing via equation 7 are determined entirely by the element of the block-diagonal matrix D. Thoe 3N value may be computed from the tructure of the input image f (van Vliet et al., 1995), which may vary from ample to ample. Indeed, our deire to perform patially-varying moothing i the primary reaon that we cannot imply ue fat Fourier tranform to olve equation Image painting Another application of aniotropic Laplacian i image painting, in which miing ample of a painting p are computed to be conitent with the tenor in D and known ample of p. Following Claerbout (1992), let M denote a diagonal making matrix that decribe the location of the miing ample: ( 1 ; p[i] miing M[i, i] = (8) 0 ; p[i] known, and let K be the diagonal complement of M uch that K + M = I. Then to compute the miing part Mp of p we may olve (K + MG T DGM)p = (K MG T DGK)p. (9) Figure 1. A eimic image of channel (provided by Joe Stefani of Chevron). To facilitate the olution of equation 9, I contructed the matrix on the left-hand-ide to be ymmetric and poitive-definite (SPD). Any olution p to equation 9 i alo a olution to an alternative ymmetric poitiveemidefinite (SPS) ytem analogou to that propoed by Claerbout (1992, p. 178): MG T DGMp = MG T DGKp. (10) Uing K + M = I, thi ytem i equivalent to MG T DGp = 0. (11) which imply tate that the miing part of p hould vary linearly in patially-varying direction determined by eigen-decompoition of the tenor in D. In other word, olution of equation 9 yield a tructure-oriented multi-dimenional linear interpolation of the known ample Kp of the painting p. In practice, thoe known ample could be pecified interactively or determined from other information. For eimic image of the earth uburface, the known ample could correpond to location where well log are available. Figure 1 and 2 how an example of image painting via equation 9. In thi example, I painted one pixel red near the center of the image p hown in Figure 1. The olution to equation 9 i hown in Figure 2. Although thi painting i an extreme example in which only a ingle pixel i known, it demontrate that paint flow to other ample according to the image tructure repreented by the tenor in D. The painted image in Figure 2 may vary depending on the finite-difference approximation G 1 and G 2 in equation 9. In fact, my firt attempt at image painting yielded reult more like the one hown in Figure 3, in which artifact roughly orthogonal to image feature are
3 Finite-difference approximation 77 Figure 2. The image of Figure 1 painted with a finitedifference approximation to an aniotropic Laplacian. A ingle pixel inide the white circle wa contrained to be red. All other pixel were painted uing tructure etimated from the image. Figure 4. The central portion of the image of Figure 3. Painting artifact have a checkerboard pattern correponding to high wavenumber near the patial Nyquit frequencie. aniotropic Laplacian. The purpoe of thi paper i to explain thee artifact and to derive the improved approximation that were ued to obtain the painting in Figure 2. I begin by deigning a family of 2-D finite-difference approximation to iotropic Laplacian. The deign method ued facilitate the derivation of approximation to aniotropic Laplacian. I then ue thi ame method to obtain new 3-D finite-difference approximation for both iotropic and aniotropic Laplacian. Both 2-D and 3-D approximation are compact in that they approximate derivative for any image ample uing only nearet neighbor ample: 9 ample for 2-D and 27 ample for 3-D Laplacian. Both 2-D and 3-D approximation maintain the SPS property of G T DG, which enure SPD ytem of equation 7 and 9. Compact SPD finite-difference approximation lead to the mot efficient method for olving uch equation. Figure 3. The image of Figure 1 painted with a poor finitedifference approximation to an aniotropic Laplacian. A ingle pixel inide the white circle wa contrained to be red. All other pixel were painted uing tructure etimated from the image. apparent. A cloeup view in Figure 4 how that thee artifact have a checkerboard pattern correponding to high wavenumber near the patial Nyquit frequencie. 2 ISOTROPIC LAPLACIANS Let u firt conider approximation to iotropic Laplacian in two dimenion: G T G = G T 1 G 1 + G T 2 G 2. (12) Finite-difference approximation G 1 and G 2 to partial derivative are mot uccinctly decribed by tencil that pecify the weight applied to adjacent image ample. To obtain compact tencil for G T G, I conider Artifact uch a thoe highlighted in Figure 4 are caued by poor finite-difference approximation to
4 78 D. Hale approximation G 1 and G 2 with the following tencil: r G 1 =, G 2 =, r r r r 0, r + = 1. (13) The condition r + = 1 enure that G T 1 G 1 and G T 2 G 2 are 2nd-order approximation, auming unit image ampling interval. Gathering image ample with weight provided by the tencil in equation 13 i equivalent to multiplication by the matrice G 1 and G 2. Scattering with the ame weight i equivalent to multiplication by their tranpoe G T 1 and G T 2. Here i a mall fragment of a computer program that both gather and catter in thi way to compute g = G T Gf: for (int i2=1; i2<n2; ++i2) { for (int i1=1; i1<n1; ++i1) { float f1r = f[i2 ][i1 ]-f[i2 ][i1-1]; float f1 = f[i2-1][i1 ]-f[i2-1][i1-1]; float f2r = f[i2 ][i1 ]-f[i2-1][i1 ]; float f2 = f[i2 ][i1-1]-f[i2-1][i1-1]; float g1 = r*f1r+*f1; float g2 = r*f2r+*f2; // gather end float g1r = g1*r; // catter begin float g1 = g1*; float g2r = g2*r; float g2 = g2*; g[i2 ][i1 ] = g1r+g2r; g[i2 ][i1-1] -= g1r-g2; g[i2-1][i1 ] += g1-g2r; g[i2-1][i1-1] -= g1+g2; Although thi program fragment ignore boundarie, it can be eaily modified to implement zero-value or zerolope boundary condition. Gathering and cattering in thi way uing the 2 2 tencil of equation 13 i equivalent to imply gathering with the following 3 3 tencil: G T G = r (14) Compact finite-difference approximation to Laplacian are often pecified in thi way, a a 3 3 tencil for G T G (e.g., Patra and Karttunen, 2005). I focu intead on the 2 2 tencil for G 1 and G 2, the component of G, becaue thoe will facilitate compact SPS finite-difference approximation G T DG to aniotropic Laplacian with patially-varying tenor in D. The condition r 0 and r + = 1 imply a family of finite-difference approximation for which 0 r 1/4. Figure 5 how three different member of thi family, correponding to r = 0, 1/12, and 1/4, with their Fourier tranform /6-2/3-1/6-2/3 10/3-2/3-1/6-2/3-1/6-1/2 0-1/ /2 0-1/2 (a) r = 0 (b) r = 1/12 (c) r = 1/4 Figure 5. Stencil and Fourier tranform for three different finite-difference approximation to Laplacian. Fourier amplitude diplayed are clipped between zero (dark blue) and four (dark red). Wavenumber and are in the range [ π, π] radian per ample. The three approximation diplayed in Figure 5 are comparable near the origin, for low wavenumber and, where contour are nearly circular with amplitude that well approximate the ideal Fourier tranform = Thi iotropy can be important in application uch a image moothing and painting, where reult hould not depend on ampling direction. For higher wavenumber, error are viible in all three approximation, but contour for r = 1/12 appear mot nearly circular. In fact, error for thi middle approximation are well-known to be iotropic to 2nd order, and among all compact 9-point tencil thi approximation to the iotropic Laplacian operator i optimal in that ene (e.g., Patra and Karttunen, 2005). Let u define dicretization error a a function E(k) of wavenumber k by the Fourier tranform pair G T G [1 E(k)]. (15)
5 Finite-difference approximation 79 The Fourier tranform of our finite-difference tencil G T G (equation 14) i the product of the ideal factor and a econd factor [1 E(k)]. The econd factor include the error E(k) that we would like to be zero or, failing that, at leat iotropic to O( ). From the erie expanion of the Fourier tranform of G T G we find that for r = 1/12 the dicretization error E(k) = k O(k4 ) (16) i indeed iotropic to O( ). For r 1/12, the firt 2nd-order term of thi error i aniotropic; it depend on the wavenumber vector k and not jut it magnitude k. The coefficient r and of G 1 and G 2 for the optimal approximation are eaily found from the condition r 0, r + = 1 and r = 1/12. They are r = 1/2 + 1/ 6, = 1/2 1/ 6. The coefficient for r = 1/4 are imply r = = 1/2. Thee coefficient imply a imple averaging of finite difference in G 1 and G 2 defined by equation 13. Oberve however that the Fourier tranform diplayed in Figure 5c implie that thi approximation attenuate high wavenumber near Nyquit, in addition to low wavenumber near the origin. Thi obervation can help u undertand the image painting artifact in Figure 4. When the finite-difference operator of Figure 5c i applied to an input contant or checkerboard image f, the output image g = G T Gf i zero. The invere of uch an operator will tend to amplify the ame contant or checkerboard feature. 3 ANISOTROPIC LAPLACIANS From finite-difference approximation G 1 and G 2 we may alo contruct approximation D G T DG = ˆG»D T 1 G T 11 D 12 G1 2 D 12 D 22 G 2 (17) to aniotropic Laplacian. We imply gather with G, multiply by D, and catter with G T. The correponding computer program might look like thi: for (int i2=1; i2<n2; ++i2) { for (int i1=1; i1<n1; ++i1) { float f1r = f[i2 ][i1 ]-f[i2 ][i1-1]; float f1 = f[i2-1][i1 ]-f[i2-1][i1-1]; float f2r = f[i2 ][i1 ]-f[i2-1][i1 ]; float f2 = f[i2 ][i1-1]-f[i2-1][i1-1]; float f1 = r*f1r+*f1; float f2 = r*f2r+*f2; float g1 = f1*d11[i2][i1]+f2*d12[i2][i1]; float g2 = f1*d12[i2][i1]+f2*d22[i2][i1]; float g1r = g1*r; float g1 = g1*; float g2r = g2*r; float g2 = g2*; g[i2 ][i1 ] = g1r+g2r; g[i2 ][i1-1] -= g1r-g2; g[i2-1][i1 ] += g1-g2r; g[i2-1][i1-1] -= g1+g2; The gather-catter ymmetry in thi program enure an SPS implementation of G T DG. A above, thi implementation i equivalent to one that gather with a 9-point tencil. However, the latter implementation would require a more complex evaluation of tencil coefficient that vary from ample to ample inide the innermot loop. 3.1 The problem Although the approximation implemented by thi computer program i SPS, it i a rather poor approximation except for the pecial cae r = 1/4. To ee the problem, let u firt conider the approximation with r = 0, for three different et of contant coefficient d 11, d 12 and d 22 of the tenor D. Figure 6a how the 9-point tencil for the cae d 11 = 1, d 12 = d 22 = 0. A expected, thi tencil approximate a econd derivative in the direction of the vertical axi x 1. The problem i apparent in the other two tencil hown in Figure 6b and 6c; thee tencil approximate econd derivative in the direction x 1 = x 2 and x 1 = x 2, repectively. Both tencil approximate the appropriate derivative; but they do o with very different dicretization error. The tencil in Figure 6c hould ideally be a rotated verion of the tencil in Figure 6b, but thi clearly i not the cae. Image moothing or painting with uch tencil will yield urpriing difference for feature oriented at angle of 45 degree and +45 degree with repect to the ampling grid. Stencil for r = 1/12 how the ame problem, a indicated in Figure 7. Although r = 1/12 i optimal for the iotropic cae d 11 = d 22 = 1, d 12 = 0 hown in Figure 5, thi value i clearly not optimal for the aniotropic cae hown in Figure 7. Stencil for r = 1/4, hown in Figure 8, do how the expected (±45-degree) ymmetry and in thi ene are better approximation to aniotropic Laplacian. However, they unfortunately alo exhibit the ame attenuation at high wavenumber near Nyquit that we have already een in Figure 5. Thi attenuation of high wavenumber in approximation to aniotropic Laplacian will caue amplification of thoe ame wavenumber in image moothing and painting application. Thee wavenumber near the patial Nyquit frequencie are apparent in Figure 4, where they appear a a checkboard pattern of painting artifact. Thee high wavenumber are enhanced when we invert approximation like thoe hown in Figure 8.
6 80 D. Hale /12-5/6-1/ /6 5/3 1/ (a) d 11 = 1, d 12 = d 22 = 0-1/12-5/6-1/12 (a) d 11 = 1, d 12 = d 22 = /2-1/6-2/3 1/ /3 7/3-2/3 1/2-1 0 (b) d 11 = d 12 = d 22 = 1/2 1/3-2/3-1/6 (b) d 11 = d 12 = d 22 = 1/ / / /2 0 0 (c) d 11 = d 12 = d 22 = 1/2-1/2 0 0 (c) d 11 = d 12 = d 22 = 1/2 Figure 6. Finite-difference approximation G T DG for r = 0. Figure 7. Finite-difference approximation G T DG for r = 1/ A olution Recall that the approximation for r = 1/12 i in one ene optimal; it i the mot iotropic approximation of the iotropic Laplacian. For aniotropic Laplacian we can improve thi approximation by conidering two additional approximation G 1 and G 2 to the gradient operator. G 1a = G 1b = r r r r, G 2a = r, G 2b = r r r,, (18) In hindight, we had no reaon above to prefer one or the other of the approximation G 1a or G 1b ; likewie for G 2a and G 2b. In equation 13 above, I choe the tencil G 1a and G 2a, but I might jut a well have choen G 1b or G 2b. With two choice for each of the finite-difference approximation G 1 and G 2, we can form a total of four different G T DG. Let u average all four to obtain D 1»D ˆGT 1a G T 11 D 12 G1a 2a 4 D 12 D 22 G 2a + ˆG»D T 1a G T 11 D 12 G1a 2b D 12 D 22 G 2b + ˆG»D T 1b G T 11 D 12 G1b 2a D 12 D 22 G 2a + ˆG»D «T 1b G T 11 D 12 G1b 2b. (19) D 12 D 22 G 2b Thi average of four SPS matrice i clearly SPS. For r = 1/12 and contant tenor coefficient, thi average yield the 9-point tencil hown in Figure 9. Thee tencil exhibit deired ymmetrie the tencil in Figure 9c i a rotated verion of that in Figure 9b and they do not attenuate wavenumber near the patial Nyquit frequencie. The finite-difference approx-
7 Finite-difference approximation 81-1/4-1/2-1/4-1/12-5/6-1/12 1/2 1 1/2 1/6 5/3 1/6-1/4-1/2-1/4 (a) d 11 = 1, d 12 = d 22 = 0-1/12-5/6-1/12 (a) d 11 = 1, d 12 = d 22 = 0-1/ /3-1/3 1/ /3 5/3-1/ /2 (b) d 11 = d 12 = d 22 = 1/2 1/6-1/3-1/3 (b) d 11 = d 12 = d 22 = 1/ /2 1/6-1/3-1/ /3 5/3-1/3-1/2 0 0 (c) d 11 = d 12 = d 22 = 1/2-1/3-1/3 1/6 (c) d 11 = d 12 = d 22 = 1/2 Figure 8. Finite-difference approximation G T DG for r = 1/4. Figure 9. Finite-difference approximation obtained by averaging four different G T DG for r = 1/12. imation in equation 19 i the one that I ued in the image painting example of Figure 2. The ame average of four approximation for r = 0 yield a lightly different finite-difference approximation that wa derived from a mixed finite-element method by Arbogat et al. (1997). I favor the approximation hown here with r = 1/12 becaue a dicued above it error are more iotropic. Implementation of the approximation in equation 19 need not require a factor of four increae in computational cot, becaue the averaging can be performed analytically. The derivation i tediou but the reult i imple, even for the cae of patially varying tenor coefficient. For that cae, define ymmetric tenor D at the four corner of an image ample with indice (i 1, i 2) a follow: D(i 1 1, 2 i2 1 ) a00 b 00 2 b 00 c 00 D(i 1 1, 2 i2 + 1 ) a10 b 10 2 b 10 c 10 D(i 1 + 1, 2 i2 1 ) a01 b 01 2 b 01 c 01 D(i 1 + 1, 2 i2 + 1 ) a11 b 11 2 b 11 c 11 Then, uing the finite-difference approximation of equation 19, the 9-point tencil centered on the image ample
8 82 D. Hale with indice (i 1, i 2) i r b 00 a 00 a 10 b 10 a 00 + a 01 + a 10 + a 11 c 00 c 01 +b 00 b 01 b 10 + b 11 c 10 c 11 +c 00 + c 01 + c 10 + c 11 b 01 a 01 a 11 b 11 a 00 + a 10 a 00 c 00 +c 00 + c 10 a 10 c 10 a 00 + a 01 a 00 a 01 a 10 a 11 a 10 + a 11 +c 00 + c 01 c 00 c 01 c 10 c 11 +c 10 + c 11 a 01 c 01 a 01 + a 11 a 11 c 11 +c 01 + c 11 For r = 0 only the firt part of thi tencil i ignificant, and thi part i the finite-difference approximation of Arbogat et al. (1997). Addition of the econd part with r = 1/12 yield a finite-difference approximation with error that are more iotropic. For contant tenor coefficient a = d 11, b = d 12 and c = d 22 thi tencil i conitent with the three example in Figure 9. For patially varying tenor coefficient, the variou um and difference in thi 9-point tencil maintain the important SPS property of our finite-difference approximation. A hown above, a impler way to maintain thi SPS property i to gather and catter conitently, a illutrated by the following program fragment: for (int i2=1; i2<n2; ++i2) { for (int i1=1; i1<n1; ++i1) { float a = 0.5f*d11[i2][i1]; float b = 0.5f*d12[i2][i1]; float c = 0.5f*d22[i2][i1]; float t = 2.0f*r*(a+c); float fpp = f[i2 ][i1 ]; float fpm = f[i2 ][i1-1]; float fmp = f[i2-1][i1 ]; float fmm = f[i2-1][i1-1]; float apppm = (a-t)*(fpp-fpm); float ampmm = (a-t)*(fmp-fmm); float bppmm = (b+t)*(fpp-fmm); float bpmmp = (b-t)*(fpm-fmp); float cppmp = (c-t)*(fpp-fmp); float cpmmm = (c-t)*(fpm-fmm); g[i2 ][i1 ] = apppm+bppmm+cppmp; g[i2 ][i1-1] -= apppm+bpmmp-cpmmm; g[i2-1][i1 ] += ampmm+bpmmp-cppmp; g[i2-1][i1-1] -= ampmm+bppmm+cpmmm; Note that thi implementation require for each image ample only one evaluation of tenor element D LAPLACIANS The methodology ued above to derive a finite-difference approximation to a 2-D aniotropic Laplacian extend naturally to three dimenion. The two important tep are (i) Find approximation G 1, G 2, and G 3 to component of the gradient uch that G T G = G T 1 G 1 + G T 2 G 2 + G T 3 G 3 ha (to econd order) iotropic dicretization error. (ii) Average all poible combination of thee component in approximation G T DG to the aniotropic Laplacian. The combination G T DG will have the form G T DG = ˆG D 11 D 12 D 13 G 1 T 1 G T 2 G T 4 3 D 12 D 22 D G 2 5. (20) D 13 D 23 D 33 G 3 Compact tencil for G 1, G 2, and G 3 are rotated verion of each other: G 1 = -t - t G 2 = -t t G 3 = - -t - - -r - -r r - -r r r t 0, rt = 2, r t = 1. (21) In thi notation, the back part of each tencil i pecified left of the front part. The condition rt = 2 follow from deired ymmetrie in G T 1 G 1, G T 2 G 2 and G T 3 G 3. The condition r t = 1 enure that G T 1 G 1, G T 2 G 2 and G T 3 G 3 are 2nd-order approximation. Thee tencil for G 1, G 2, and G 3 lead to compact point tencil for G T DG. In the iotropic cae where D = I, we have t G T G = G T 1 G 1 + G T 2 G 2 + G T 3 G 3 (22) For thi cae Patra and Karttunen (2005) give condition for which compact 3-D tencil have (to 2nd order) iotropic dicretization error, and they lit example of tencil that meet thoe condition. In three dimenion, unlike two dimenion, more than one uch tencil i poible. However, none of the iotropic tencil they cite can be expreed in the form of equation 22, in term of approximation to component of gradient. In other word, they do not correpond to approximation G 1, G 2, and G 3 that we could ue in the aniotropic form of equation 20. By carefully chooing the coefficient r, and t ubject to the condition cited above, I found a new finite- r,,,
9 Finite-difference approximation 83 difference approximation with the deired form. The coefficient are t = 5/12 1/ 6, r = (1 t) 2, = rt. Thee coefficient lead to the following 27-point tencil: G T G = c 1 c 2 c 1 c 2 c 3 c 2 c 1 c 2 c 1 c 2 c 3 c 2 c 3 c 4 c 3 c 2 c 3 c 2 c 1 c 2 c 1 c 2 c 3 c 2, c 1 c 2 c 1 c 1 = 1/48, c 2 = 1/8, c 3 = 5/12, c 4 = 25/6. (23) Thi tencil meet the condition for iotropic dicretization error pecified by Patra and Karttunen (2005), who cite a different 27-point tencil propoed by Spotz and Carey (1995): c 1 = 1/30, c 2 = 1/10, c 3 = 7/15, c 4 = 64/15. A noted above, thi different tencil doe not have correponding approximation G 1, G 2 and G 3 that we can ue in equation 20. Before inerting the tencil G 1, G 2 and G 3 of equation 21 into equation 20, we firt recognize that each of thee 3-D tencil can be written in four different way. Thoe four way are analogou to the two different way we expreed the 2-D tencil G 1 and G 2 in equation 18. Since each of the 3-D tencil G 1, G 2 and G 3 can be written four different way, we have a total of 64 = combination of the form of equation 20. By averaging all 64 combination, we obtain a compoite 3-D finite-difference approximation for an aniotropic Laplacian. Thi large number of combination to average can be reduced if we conider only pairwie combination uch a G T 1 D 11G 1 (four way) or G T 1 D 12G 2 (ixteen way). In any cae, the averaging can again be performed analytically and, although too lengthy to be included in thi paper, the reult i analogou to that obtained above by averaging 2-D approximation. 5 CONCLUSION The method decribed here for finite-difference approximation of aniotropic Laplacian i traightforward. We begin with compact finite-difference approximation to component of the gradient. We then average ymmetric poitive-emidefinite combination of thoe gradient approximation to obtain the deired approximation to aniotropic Laplacian. Finite-difference approximation obtained in thi way are guaranteed to be ymmetric and poitive-emidefinite. Thi deign method yield 9-point tencil for 2-D approximation and 27-point tencil for 3-D approximation of aniotropic Laplacian. Our 2-D approximation i a generalization of that propoed by Arbogat et al. (1997), in that for one particular gradient approximation we obtain the coefficient of their 9-point tencil. However, an alternative gradient approximation in our method lead to a 9-point tencil with iotropic dicretization error. Arbogat et al. (1997) alo propoed a 19-point tencil for a 3-D approximation, and our 27-point tencil i not a generalization of their. Indeed, I have not found a way to modify the deign propoed here to obtain uch a 19-point tencil. A 19-point tencil i attractive becaue it implie lower computational cot. In both 2-D and 3-D approximation, required ymmetry in the reulting 9-point or 27-point finitedifference tencil implie only one free parameter in the correponding gradient approximation. I chooe thi parameter to obtain tencil with iotropic dicretization error. All finite-difference method exhibit dicretization error. By chooing method for which thoe error are iotropic to leading order, we may reduce artifact aociated with ampling grid in application uch a image moothing and image painting. REFERENCES Arbogat, T., M.F. Wheeler and I. Yotov, 1997, Mixed finite element for elliptic problem with tenor coefficient a cell-centered finite-difference: SIAM Journal of Numerical Analyi, 34, Claerbout, J.F., 1992, Earth ounding analyi proceing veru inverion: Blackwell Scientific Publication. Fehmer, G.C., C.F.W. Höcker, 2003, Fat tructural interpretation with tructure-oriented filtering: Geophyic, 68, Patra, M. and M. Karttunen, 2005, Stencil with iotropic dicretization error for differential operator: Numerical Method for Partial Differential Equation, 22, Spotz, W.F. and G.F. Carey, 1996, A high-order compact formulation for the 3D Poion equation: Numerical Method for Partial Differential Equation, 12, van Vliet, L.J., and P.W. Verbeek, 1995, Etimator for orientation and aniotropy in digitized image: Proceeding of the firt annual conference of the Advanced School for Computing and Imaging ASCI 95, Heijen (The Netherland),
10 84 D. Hale
CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS
CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3
More informationinto a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get
Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }}
More informationChapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog
Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou
More informationBogoliubov Transformation in Classical Mechanics
Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How
More informationLecture 10 Filtering: Applied Concepts
Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering
More informationDesign By Emulation (Indirect Method)
Deign By Emulation (Indirect Method he baic trategy here i, that Given a continuou tranfer function, it i required to find the bet dicrete equivalent uch that the ignal produced by paing an input ignal
More informationFast Convolutional Sparse Coding (FCSC)
Fat Convolutional Spare Coding (FCSC) Bailey ong Department of Computer Science Univerity of California, Irvine bhkong@ic.uci.edu Charle C. Fowlke Department of Computer Science Univerity of California,
More information7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281
72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition
More informationDesign of Digital Filters
Deign of Digital Filter Paley-Wiener Theorem [ ] ( ) If h n i a caual energy ignal, then ln H e dω< B where B i a finite upper bound. One implication of the Paley-Wiener theorem i that a tranfer function
More informationSolutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam
BSc - Sample Examination Digital Control Sytem (5-588-) Prof. L. Guzzella Solution Exam Duration: Number of Quetion: Rating: Permitted aid: minute examination time + 5 minute reading time at the beginning
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction
More informationControl Systems Analysis and Design by the Root-Locus Method
6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If
More informationMATEMATIK Datum: Tid: eftermiddag. A.Heintz Telefonvakt: Anders Martinsson Tel.:
MATEMATIK Datum: 20-08-25 Tid: eftermiddag GU, Chalmer Hjälpmedel: inga A.Heintz Telefonvakt: Ander Martinon Tel.: 073-07926. Löningar till tenta i ODE och matematik modellering, MMG5, MVE6. Define what
More informationS_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS
S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS by Michelle Gretzinger, Daniel Zyngier and Thoma Marlin INTRODUCTION One of the challenge to the engineer learning proce control i relating theoretical
More informationLecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)
Lecture 7: Analytic Function and Integral (See Chapter 4 in Boa) Thi i a good point to take a brief detour and expand on our previou dicuion of complex variable and complex function of complex variable.
More informationOne Class of Splitting Iterative Schemes
One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi
More informationSource slideplayer.com/fundamentals of Analytical Chemistry, F.J. Holler, S.R.Crouch. Chapter 6: Random Errors in Chemical Analysis
Source lideplayer.com/fundamental of Analytical Chemitry, F.J. Holler, S.R.Crouch Chapter 6: Random Error in Chemical Analyi Random error are preent in every meaurement no matter how careful the experimenter.
More informationSocial Studies 201 Notes for November 14, 2003
1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationClustering Methods without Given Number of Clusters
Clutering Method without Given Number of Cluter Peng Xu, Fei Liu Introduction A we now, mean method i a very effective algorithm of clutering. It mot powerful feature i the calability and implicity. However,
More informationA Simple Approach to Synthesizing Naïve Quantized Control for Reference Tracking
A Simple Approach to Syntheizing Naïve Quantized Control for Reference Tracking SHIANG-HUA YU Department of Electrical Engineering National Sun Yat-Sen Univerity 70 Lien-Hai Road, Kaohiung 804 TAIAN Abtract:
More informationMath Skills. Scientific Notation. Uncertainty in Measurements. Appendix A5 SKILLS HANDBOOK
ppendix 5 Scientific Notation It i difficult to work with very large or very mall number when they are written in common decimal notation. Uually it i poible to accommodate uch number by changing the SI
More informationAvoiding Forbidden Submatrices by Row Deletions
Avoiding Forbidden Submatrice by Row Deletion Sebatian Wernicke, Jochen Alber, Jen Gramm, Jiong Guo, and Rolf Niedermeier Wilhelm-Schickard-Intitut für Informatik, niverität Tübingen, Sand 13, D-72076
More informationThe Hassenpflug Matrix Tensor Notation
The Haenpflug Matrix Tenor Notation D.N.J. El Dept of Mech Mechatron Eng Univ of Stellenboch, South Africa e-mail: dnjel@un.ac.za 2009/09/01 Abtract Thi i a ample document to illutrate the typeetting of
More informationSocial Studies 201 Notes for March 18, 2005
1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationGNSS Solutions: What is the carrier phase measurement? How is it generated in GNSS receivers? Simply put, the carrier phase
GNSS Solution: Carrier phae and it meaurement for GNSS GNSS Solution i a regular column featuring quetion and anwer about technical apect of GNSS. Reader are invited to end their quetion to the columnit,
More informationGain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays
Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,
More informationSampling and the Discrete Fourier Transform
Sampling and the Dicrete Fourier Tranform Sampling Method Sampling i mot commonly done with two device, the ample-and-hold (S/H) and the analog-to-digital-converter (ADC) The S/H acquire a CT ignal at
More informationImproving the Efficiency of a Digital Filtering Scheme for Diabatic Initialization
1976 MONTHLY WEATHER REVIEW VOLUME 15 Improving the Efficiency of a Digital Filtering Scheme for Diabatic Initialization PETER LYNCH Met Éireann, Dublin, Ireland DOMINIQUE GIARD CNRM/GMAP, Météo-France,
More informationAn Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem
An Inequality for Nonnegative Matrice and the Invere Eigenvalue Problem Robert Ream Program in Mathematical Science The Univerity of Texa at Dalla Box 83688, Richardon, Texa 7583-688 Abtract We preent
More informationSingular perturbation theory
Singular perturbation theory Marc R. Rouel June 21, 2004 1 Introduction When we apply the teady-tate approximation (SSA) in chemical kinetic, we typically argue that ome of the intermediate are highly
More informationA Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems
A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement
More informationEfficient Methods of Doppler Processing for Coexisting Land and Weather Clutter
Efficient Method of Doppler Proceing for Coexiting Land and Weather Clutter Ça gatay Candan and A Özgür Yılmaz Middle Eat Technical Univerity METU) Ankara, Turkey ccandan@metuedutr, aoyilmaz@metuedutr
More informationConvex Hulls of Curves Sam Burton
Convex Hull of Curve Sam Burton 1 Introduction Thi paper will primarily be concerned with determining the face of convex hull of curve of the form C = {(t, t a, t b ) t [ 1, 1]}, a < b N in R 3. We hall
More informationDigital Control System
Digital Control Sytem - A D D A Micro ADC DAC Proceor Correction Element Proce Clock Meaurement A: Analog D: Digital Continuou Controller and Digital Control Rt - c Plant yt Continuou Controller Digital
More informationCHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL
98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i
More informationLecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004
18.997 Topic in Combinatorial Optimization April 29th, 2004 Lecture 21 Lecturer: Michel X. Goeman Scribe: Mohammad Mahdian 1 The Lovaz plitting-off lemma Lovaz plitting-off lemma tate the following. Theorem
More informationLecture 9: Shor s Algorithm
Quantum Computation (CMU 8-859BB, Fall 05) Lecture 9: Shor Algorithm October 7, 05 Lecturer: Ryan O Donnell Scribe: Sidhanth Mohanty Overview Let u recall the period finding problem that wa et up a a function
More informationLearning Multiplicative Interactions
CSC2535 2011 Lecture 6a Learning Multiplicative Interaction Geoffrey Hinton Two different meaning of multiplicative If we take two denity model and multiply together their probability ditribution at each
More informationLecture 4 Topic 3: General linear models (GLMs), the fundamentals of the analysis of variance (ANOVA), and completely randomized designs (CRDs)
Lecture 4 Topic 3: General linear model (GLM), the fundamental of the analyi of variance (ANOVA), and completely randomized deign (CRD) The general linear model One population: An obervation i explained
More informationProblem Set 8 Solutions
Deign and Analyi of Algorithm April 29, 2015 Maachuett Intitute of Technology 6.046J/18.410J Prof. Erik Demaine, Srini Devada, and Nancy Lynch Problem Set 8 Solution Problem Set 8 Solution Thi problem
More informationIEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation
IEOR 316: Fall 213, Profeor Whitt Topic for Dicuion: Tueday, November 19 Alternating Renewal Procee and The Renewal Equation 1 Alternating Renewal Procee An alternating renewal proce alternate between
More informationQuestion 1 Equivalent Circuits
MAE 40 inear ircuit Fall 2007 Final Intruction ) Thi exam i open book You may ue whatever written material you chooe, including your cla note and textbook You may ue a hand calculator with no communication
More informationUSPAS Course on Recirculated and Energy Recovered Linear Accelerators
USPAS Coure on Recirculated and Energy Recovered Linear Accelerator G. A. Krafft and L. Merminga Jefferon Lab I. Bazarov Cornell Lecture 6 7 March 005 Lecture Outline. Invariant Ellipe Generated by a Unimodular
More informationSuggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall
Beyond Significance Teting ( nd Edition), Rex B. Kline Suggeted Anwer To Exercie Chapter. The tatitic meaure variability among core at the cae level. In a normal ditribution, about 68% of the core fall
More informationChapter 7. Root Locus Analysis
Chapter 7 Root Locu Analyi jw + KGH ( ) GH ( ) - K 0 z O 4 p 2 p 3 p Root Locu Analyi The root of the cloed-loop characteritic equation define the ytem characteritic repone. Their location in the complex
More informationWhat lies between Δx E, which represents the steam valve, and ΔP M, which is the mechanical power into the synchronous machine?
A 2.0 Introduction In the lat et of note, we developed a model of the peed governing mechanim, which i given below: xˆ K ( Pˆ ˆ) E () In thee note, we want to extend thi model o that it relate the actual
More informationCorrection for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002
Correction for Simple Sytem Example and Note on Laplace Tranform / Deviation Variable ECHE 55 Fall 22 Conider a tank draining from an initial height of h o at time t =. With no flow into the tank (F in
More informationSIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II. Solutions to Assignment 3 February 2005.
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuit II Solution to Aignment 3 February 2005. Initial Condition Source 0 V battery witch flip at t 0 find i 3 (t) Component value:
More informationIII.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SUBSTANCES
III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SBSTANCES. Work purpoe The analyi of the behaviour of a ferroelectric ubtance placed in an eternal electric field; the dependence of the electrical polariation
More informationOptimal Coordination of Samples in Business Surveys
Paper preented at the ICES-III, June 8-, 007, Montreal, Quebec, Canada Optimal Coordination of Sample in Buine Survey enka Mach, Ioana Şchiopu-Kratina, Philip T Rei, Jean-Marc Fillion Statitic Canada New
More informationDYNAMIC MODELS FOR CONTROLLER DESIGN
DYNAMIC MODELS FOR CONTROLLER DESIGN M.T. Tham (996,999) Dept. of Chemical and Proce Engineering Newcatle upon Tyne, NE 7RU, UK.. INTRODUCTION The problem of deigning a good control ytem i baically that
More informationAccuracy of Symmetric Partitioned Runge-Kutta Methods for Differential Equations on Lie-Groups
Accuracy of Symmetric Partitioned Runge-Kutta Method for Differential Equation on Lie-Group WUB / 11-23 BUW-IMACM 11-19, Michael Günther, Franceco Knechtli and Michèle Wandelt Bergiche Univerität Wuppertal,
More informationDigital Control System
Digital Control Sytem Summary # he -tranform play an important role in digital control and dicrete ignal proceing. he -tranform i defined a F () f(k) k () A. Example Conider the following equence: f(k)
More information1. The F-test for Equality of Two Variances
. The F-tet for Equality of Two Variance Previouly we've learned how to tet whether two population mean are equal, uing data from two independent ample. We can alo tet whether two population variance are
More informationSpring 2014 EE 445S Real-Time Digital Signal Processing Laboratory. Homework #0 Solutions on Review of Signals and Systems Material
Spring 4 EE 445S Real-Time Digital Signal Proceing Laboratory Prof. Evan Homework # Solution on Review of Signal and Sytem Material Problem.. Continuou-Time Sinuoidal Generation. In practice, we cannot
More informationAn estimation approach for autotuning of event-based PI control systems
Acta de la XXXIX Jornada de Automática, Badajoz, 5-7 de Septiembre de 08 An etimation approach for autotuning of event-baed PI control ytem Joé Sánchez Moreno, María Guinaldo Loada, Sebatián Dormido Departamento
More informationMarch 18, 2014 Academic Year 2013/14
POLITONG - SHANGHAI BASIC AUTOMATIC CONTROL Exam grade March 8, 4 Academic Year 3/4 NAME (Pinyin/Italian)... STUDENT ID Ue only thee page (including the back) for anwer. Do not ue additional heet. Ue of
More informationμ + = σ = D 4 σ = D 3 σ = σ = All units in parts (a) and (b) are in V. (1) x chart: Center = μ = 0.75 UCL =
Our online Tutor are available 4*7 to provide Help with Proce control ytem Homework/Aignment or a long term Graduate/Undergraduate Proce control ytem Project. Our Tutor being experienced and proficient
More informationEE Control Systems LECTURE 14
Updated: Tueday, March 3, 999 EE 434 - Control Sytem LECTURE 4 Copyright FL Lewi 999 All right reerved ROOT LOCUS DESIGN TECHNIQUE Suppoe the cloed-loop tranfer function depend on a deign parameter k We
More informationCodes Correcting Two Deletions
1 Code Correcting Two Deletion Ryan Gabry and Frederic Sala Spawar Sytem Center Univerity of California, Lo Angele ryan.gabry@navy.mil fredala@ucla.edu Abtract In thi work, we invetigate the problem of
More informationChapter 5 Consistency, Zero Stability, and the Dahlquist Equivalence Theorem
Chapter 5 Conitency, Zero Stability, and the Dahlquit Equivalence Theorem In Chapter 2 we dicued convergence of numerical method and gave an experimental method for finding the rate of convergence (aka,
More informationCHAPTER 6. Estimation
CHAPTER 6 Etimation Definition. Statitical inference i the procedure by which we reach a concluion about a population on the bai of information contained in a ample drawn from that population. Definition.
More informationSMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD
SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD S.P. Teeuwen, I. Erlich U. Bachmann Univerity of Duiburg, Germany Department of Electrical Power Sytem
More informationWhite Rose Research Online URL for this paper: Version: Accepted Version
Thi i a repoitory copy of Identification of nonlinear ytem with non-peritent excitation uing an iterative forward orthogonal leat quare regreion algorithm. White Roe Reearch Online URL for thi paper: http://eprint.whiteroe.ac.uk/107314/
More informationTHE THERMOELASTIC SQUARE
HE HERMOELASIC SQUARE A mnemonic for remembering thermodynamic identitie he tate of a material i the collection of variable uch a tre, train, temperature, entropy. A variable i a tate variable if it integral
More information"HIP Modeling Methodology Based on the Inherent Process Anisotropy
"HIP Modeling Methodology Baed on the Inherent Proce Aniotropy Victor Samarov, Vaily Golovehkin, Charle Barre, ( LNT PM, Syneretch P/M, Inc., 65 Monarch treet, Garden Grove CA, USA, 984) Abtract The net
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. VIII Decoupling Control - M. Fikar
DECOUPLING CONTROL M. Fikar Department of Proce Control, Faculty of Chemical and Food Technology, Slovak Univerity of Technology in Bratilava, Radlinkého 9, SK-812 37 Bratilava, Slovakia Keyword: Decoupling:
More informationEuler-Bernoulli Beams
Euler-Bernoulli Beam The Euler-Bernoulli beam theory wa etablihed around 750 with contribution from Leonard Euler and Daniel Bernoulli. Bernoulli provided an expreion for the train energy in beam bending,
More informationCDMA Signature Sequences with Low Peak-to-Average-Power Ratio via Alternating Projection
CDMA Signature Sequence with Low Peak-to-Average-Power Ratio via Alternating Projection Joel A Tropp Int for Comp Engr and Sci (ICES) The Univerity of Texa at Autin 1 Univerity Station C0200 Autin, TX
More informationAnnex-A: RTTOV9 Cloud validation
RTTOV-91 Science and Validation Plan Annex-A: RTTOV9 Cloud validation Author O Embury C J Merchant The Univerity of Edinburgh Intitute for Atmo. & Environ. Science Crew Building King Building Edinburgh
More informationEffects of vector attenuation on AVO of offshore reflections
GEOPHYSICS, VOL. 64, NO. 3 MAY-JUNE 1999); P. 815 819, 9 FIGS., 1 TABLE. Effect of vector attenuation on AVO of offhore reflection J. M. Carcione ABSTRACT Wave tranmitted at the ocean bottom have the characteritic
More informationProperties of Z-transform Transform 1 Linearity a
Midterm 3 (Fall 6 of EEG:. Thi midterm conit of eight ingle-ided page. The firt three page contain variou table followed by FOUR eam quetion and one etra workheet. You can tear out any page but make ure
More informationPreemptive scheduling on a small number of hierarchical machines
Available online at www.ciencedirect.com Information and Computation 06 (008) 60 619 www.elevier.com/locate/ic Preemptive cheduling on a mall number of hierarchical machine György Dóa a, Leah Eptein b,
More informationIntroduction to Laplace Transform Techniques in Circuit Analysis
Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found
More informationMacromechanical Analysis of a Lamina
3, P. Joyce Macromechanical Analyi of a Lamina Generalized Hooke Law ij Cijklε ij C ijkl i a 9 9 matri! 3, P. Joyce Hooke Law Aume linear elatic behavior mall deformation ε Uniaial loading 3, P. Joyce
More informationHomework #7 Solution. Solutions: ΔP L Δω. Fig. 1
Homework #7 Solution Aignment:. through.6 Bergen & Vittal. M Solution: Modified Equation.6 becaue gen. peed not fed back * M (.0rad / MW ec)(00mw) rad /ec peed ( ) (60) 9.55r. p. m. 3600 ( 9.55) 3590.45r.
More informationAutomatic Control Systems. Part III: Root Locus Technique
www.pdhcenter.com PDH Coure E40 www.pdhonline.org Automatic Control Sytem Part III: Root Locu Technique By Shih-Min Hu, Ph.D., P.E. Page of 30 www.pdhcenter.com PDH Coure E40 www.pdhonline.org VI. Root
More informationFinding the location of switched capacitor banks in distribution systems based on wavelet transform
UPEC00 3t Aug - 3rd Sept 00 Finding the location of witched capacitor bank in ditribution ytem baed on wavelet tranform Bahram nohad Shahid Chamran Univerity in Ahvaz bahramnohad@yahoo.com Mehrdad keramatzadeh
More informationStandard Guide for Conducting Ruggedness Tests 1
Deignation: E 69 89 (Reapproved 996) Standard Guide for Conducting Ruggedne Tet AMERICA SOCIETY FOR TESTIG AD MATERIALS 00 Barr Harbor Dr., Wet Conhohocken, PA 948 Reprinted from the Annual Book of ASTM
More informationCoupling of Three-Phase Sequence Circuits Due to Line and Load Asymmetries
Coupling of Three-Phae Sequence Circuit Due to Line and Load Aymmetrie DEGO BELLAN Department of Electronic nformation and Bioengineering Politecnico di Milano Piazza Leonardo da inci 01 Milano TALY diego.ellan@polimi.it
More informationRaneNote BESSEL FILTER CROSSOVER
RaneNote BESSEL FILTER CROSSOVER A Beel Filter Croover, and It Relation to Other Croover Beel Function Phae Shift Group Delay Beel, 3dB Down Introduction One of the way that a croover may be contructed
More informationAdvanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment
Journal of Multidiciplinary Engineering Science and Technology (JMEST) ISSN: 59- Vol. Iue, January - 5 Advanced D-Partitioning Analyi and it Comparion with the haritonov Theorem Aement amen M. Yanev Profeor,
More informationLecture 7: Testing Distributions
CSE 5: Sublinear (and Streaming) Algorithm Spring 014 Lecture 7: Teting Ditribution April 1, 014 Lecturer: Paul Beame Scribe: Paul Beame 1 Teting Uniformity of Ditribution We return today to property teting
More informationHigh-field behavior: the law of approach to saturation (Is there an equation for the magnetization at high fields?)
High-field behavior: the law of approach to aturation (I there an equation for the magnetization at high field? In the high-field region the magnetization approache aturation. The firt attempt to give
More informationarxiv: v2 [nucl-th] 3 May 2018
DAMTP-207-44 An Alpha Particle Model for Carbon-2 J. I. Rawlinon arxiv:72.05658v2 [nucl-th] 3 May 208 Department of Applied Mathematic and Theoretical Phyic, Univerity of Cambridge, Wilberforce Road, Cambridge
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Laplace Tranform Paul Dawkin Table of Content Preface... Laplace Tranform... Introduction... The Definition... 5 Laplace Tranform... 9 Invere Laplace Tranform... Step Function...4
More informationHalliday/Resnick/Walker 7e Chapter 6
HRW 7e Chapter 6 Page of Halliday/Renick/Walker 7e Chapter 6 3. We do not conider the poibility that the bureau might tip, and treat thi a a purely horizontal motion problem (with the peron puh F in the
More informationStochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions
Stochatic Optimization with Inequality Contraint Uing Simultaneou Perturbation and Penalty Function I-Jeng Wang* and Jame C. Spall** The John Hopkin Univerity Applied Phyic Laboratory 11100 John Hopkin
More informationPhysics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2014
Phyic 7 Graduate Quantum Mechanic Solution to inal Eam all 0 Each quetion i worth 5 point with point for each part marked eparately Some poibly ueful formula appear at the end of the tet In four dimenion
More informationChapter 4. The Laplace Transform Method
Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination
More informationDesign spacecraft external surfaces to ensure 95 percent probability of no mission-critical failures from particle impact.
PREFERRED RELIABILITY PAGE 1 OF 6 PRACTICES METEOROIDS & SPACE DEBRIS Practice: Deign pacecraft external urface to enure 95 percent probability of no miion-critical failure from particle impact. Benefit:
More informationON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang
Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang
More informationSERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lectures 41-48)
Chapter 5 SERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lecture 41-48) 5.1 Introduction Power ytem hould enure good quality of electric power upply, which mean voltage and current waveform hould
More informationFUNDAMENTALS OF POWER SYSTEMS
1 FUNDAMENTALS OF POWER SYSTEMS 1 Chapter FUNDAMENTALS OF POWER SYSTEMS INTRODUCTION The three baic element of electrical engineering are reitor, inductor and capacitor. The reitor conume ohmic or diipative
More informationBy Xiaoquan Wen and Matthew Stephens University of Michigan and University of Chicago
Submitted to the Annal of Applied Statitic SUPPLEMENTARY APPENDIX TO BAYESIAN METHODS FOR GENETIC ASSOCIATION ANALYSIS WITH HETEROGENEOUS SUBGROUPS: FROM META-ANALYSES TO GENE-ENVIRONMENT INTERACTIONS
More informationUnified Design Method for Flexure and Debonding in FRP Retrofitted RC Beams
Unified Deign Method for Flexure and Debonding in FRP Retrofitted RC Beam G.X. Guan, Ph.D. 1 ; and C.J. Burgoyne 2 Abtract Flexural retrofitting of reinforced concrete (RC) beam uing fibre reinforced polymer
More informationThe type 3 nonuniform FFT and its applications
Journal of Computational Phyic 206 (2005) 1 5 Short Note The type 3 nonuniform FFT and it application June-Yub Lee a, *, Lelie Greengard b a Department of Mathematic, Ewha Woman Univerity, 11-1 Daehyundong,
More informationECE382/ME482 Spring 2004 Homework 4 Solution November 14,
ECE382/ME482 Spring 2004 Homework 4 Solution November 14, 2005 1 Solution to HW4 AP4.3 Intead of a contant or tep reference input, we are given, in thi problem, a more complicated reference path, r(t)
More informationOnline Appendix for Managerial Attention and Worker Performance by Marina Halac and Andrea Prat
Online Appendix for Managerial Attention and Worker Performance by Marina Halac and Andrea Prat Thi Online Appendix contain the proof of our reult for the undicounted limit dicued in Section 2 of the paper,
More informationMath 273 Solutions to Review Problems for Exam 1
Math 7 Solution to Review Problem for Exam True or Fale? Circle ONE anwer for each Hint: For effective tudy, explain why if true and give a counterexample if fale (a) T or F : If a b and b c, then a c
More information