The type 3 nonuniform FFT and its applications
|
|
- Clyde Julian McDaniel
- 5 years ago
- Views:
Transcription
1 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, Seodaemoongu, Seoul, , Korea b Courant Intitute of Mathematical Science, New York Univerity, New York, NY 10012, USA Received 16 Augut 2004 received in revied form 11 October 2004 accepted 10 December 2004 Available online 22 January Abtract The nonequipaced or nonuniform fat Fourier tranform (NUFFT) arie in a variety of application area, including imaging proceing and the numerical olution of partial differential equation. In it mot general form, it take a input an irregular ampling of a function and eek to compute it Fourier tranform at a nonuniform ampling of frequency location. Thi i ometime referred to a the NUFFT of type 3. Like the fat Fourier tranform, the amount of work required i of the order O(N log N), where N denote the number of ampling point in both the phyical and pectral domain. In thi hort note, we preent the eential idea underlying the algorithm in imple term. We alo illutrate it utility with application to problem in magnetic reonance imagin and heat flow. Ó 2005 Elevier Inc. All right reerved. Keyword: Nonuniform fat Fourier tranform Fourier integral Heat equation Magnetic reonance imaging 1. Introduction The nonuniform or unequipaced fat Fourier tranform (NUFFT) arie in a variety of application area, including imaging proceing and the numerical olution of partial differential equation. In it implet form, one i given an irregular ampling of a function at N data point and eek to compute the coefficient of the correponding Fourier erie. Following Dutt and Rokhlin [2], we refer to thi a a type 1 tranform. The adjoint of thi procedure i that of evaluating a given Fourier erie at a et of nonuniform target point, which we refer to a a type 2 tranform. Like the fat Fourier tranform (FFT) for equipaced data, the NUFFT reduce the cot of the computation from O(N 2 ) operation to O(N log N) operation. Thi hort note i not intended a a review article, and we refer the reader to a ampling of the relevant * Correponding author. Tel.: fax: addree: jylee@math.ewha.ac.kr, jylee@cim.nyu.edu (J.-Y. Lee) /$ - ee front matter Ó 2005 Elevier Inc. All right reerved. doi: /j.jcp
2 2 J.-Y. Lee, L. Greengard / Journal of Computational Phyic 206 (2005) 1 5 literature in the paper [1,3,4,7,8]. Additional citation can be found in our earlier paper [5], where we decribe imple and efficient implementation of thee NUFFT. Here, we conider the mot general uch tranform, where both the phyical and Fourier pace ampling are nonuniform and dicu ome of it application. More preciely, in d dimenion, we conider the computation of F k ¼ XN 1 f j e i kx j ð1þ at N location k. Thi i ometime referred to a a type 3 verion of the NUFFT. The algorithm i not new. Decription can be found, for example, in [1,2,7,8]. Neverthele, it i worth decribing the eential feature of the method in imple term. Moreover, the type 3 NUFFT ha intereting application which do not appear to be widely appreciated. 2. The nonuniform FFT of type 3 We can think of (1) a a dicretization of the continuou Fourier tranform, F ðþ ¼ 1 Z 1 Z 1 f ðxþe ix dx ð2þ ðþ d 1 1 uing nonuniformly ampled dicretization point and evaluated at nonuniformly ampled frequencie. All exiting nonuniform FFT are baed, in eence, on combining a local interpolation cheme with the tandard FFT. The type 3 NUFFT i no exception. For the ake of implicity, however, we limit the detailed dicuion to the one dimenional cae. Note firt that Eq. (1) can be interpreted a the continuou Fourier tranform of the function p f ðxþ ¼ ffiffiffiffiffi X N 1 f j dðx x j Þ ð3þ evaluated at the point = k. We obviouly cannot ample the function f(x) on a uniform meh, o we begin by convolving f(x) with the (one-dimenional) Gauian g ðxþ ¼e x2 =4. Thu, we define f (x) by f ðxþ ¼f g ðxþ ¼ 1 Z 1 ffiffiffi f ðyþg ðx yþdy: ð4þ 1 Since f i now a mooth, infinitely differentiable function, it can be well-reolved by a uniform meh in x whoe pacing i determined by the parameter. For thi, we aume that jx j j 6 X. We define the dicrete equipaced ample of f by f ðnd x Þ¼ XN 1 f j g ðnd x x j Þ: Let u now carry out a le intuitive anti-diffuion tep, whoe purpoe will become clear hortly. For thi, we define pthe (continuou) Fourier tranform of g (x) byg (). A traightforward calculation how that G ðþ ¼ ffiffiffiffiffi 2 e 2. We let ðxþ ¼ f ðxþ G r ðxþ ¼ p 1 ffiffiffiffiffi e rx2 f ðxþ: ð6þ 2r The Fourier tranform of f r, namely ðþ ¼p 1 ffiffiffiffiffi Z 1 1 ðxþe ix dx ð5þ
3 J.-Y. Lee, L. Greengard / Journal of Computational Phyic 206 (2005) can be computed with high accuracy on a uniform grid in with pacing D uing the tandard FFT on a ufficiently fine grid. That i, ðmd Þp Dx ffiffiffiffiffi X n ðnd x Þe imndxd : ð7þ Remark 1. Here, we aume that j k j 6 S. Becaue of the rapid decay of the Gauian in (5), we can ignore contribution to f (nd x ) from point x j more than a certain ditance away with an exponentially mall error. We will denote by m p the number of grid point to which we extend the influence of a point ource in each direction. Given m p, we chooe D x p 1 S R D p 1 X þm p D x R, and M r ¼ D x D, where the overampling parameter R > 1. The actual value of R and m p can be optimized once the accuracy requirement are known. The next tep i to recover the value F ( k ) by convolving ðþ with g r (): i.e., F ð k Þ¼ g r ðk Þ¼ p 1 ffiffiffiffiffi Z 1 1 ðþg r ð k Þd D ffiffiffi X m ðmd Þ g r ð k md Þ: ð8þ Thi follow from the convolution theorem and the fact that we already carried out the deconvolution tep in (6). Once the value F ( k ) are known, it remain only to correct for the initial moothing to obtain the deired value F( k ). An elementary calculation how that F ð k Þ¼p 1 ffiffiffiffiffi e 2 k F ð k Þ: ð9þ 2 Thi follow again from the convolution theorem. The actual implementation require a complete pecification of all detail. We do not repeat the analyi of [2] here, but ummarize the relevant reult a follow: if we let 1 p ffiffiffi 2 p 1 D p ffiffi X þ m p D x 2 D x p ffi D x r ffi D ffi S 2 ð 2 1Þ 2 ð 2 1Þ ð10þ then carrying out the convolution in (5) and (8) with m p = 9 yield about ix digit of accuracy. Carrying out the convolution with m p = 18 yield about twelve digit of accuracy. Efficient implementation of thee tep can be carried out uing the fat Gauian gridding algorithm of [5]. The higher dimenional verion involve more notation but are obviou extenion of the one-dimenional cheme. Appropriate value of M r, D x and D can be choen for each dimenion eparately. m p m p M r ¼ D x D Remark 2. Greater efficiency can be achieved by uing convolution kernel other than Gauian [1,3,7], but at an increaed torage cot compared to the cheme of [5]. 3. Numerical example The NUFFT of type 1, 2 and 3 have been implemented (in Fortran) with fat gridding in one, two and three dimenion. Detailed experiment were preented in [5] for the firt two type in one and two dimenion. Here, we illutrate the performance of the type 3 tranform in the context of two application.
4 4 J.-Y. Lee, L. Greengard / Journal of Computational Phyic 206 (2005) 1 5 (a) 5 Source Point (b) 40 Frequency Sampling Fig. 1. The left-hand figure how a et of ampling point for a function f(x) and the right-hand figure how a et of dicretization point in Fourier pace. Example 1. An important application of the type 3 NUFFT i to heat flow in exterior domain. We have hown [6] that the free-pace heat kernel i well-repreented in the Fourier domain uing a highly nonuniform ampling one that cluter exponentially to the origin (Fig. 1(b)). In order to project a function onto thee Fourier mode, one mut evaluate um of the form F ðk j Þ XN w n f ðx n Þe ik jx n where the point x n are the ampling point for f(x). (The tenor product dicretization hown i eay to contruct but not entirely optimal. More ophiticated quadrature following the general approach of [6] would achieve a modet reduction in the total number of node.) Suppoe for example that f(x) i a ingular heat ource concentrated on the curve depicted in Fig. 1(a). With N = 22,500 point and a pectral dicretization uing node, the type 3 NUFFT with ix digit of accuracy require 0.4 while the direct calculation require about 110. If the data were uniformly paced, the tandard FFT would require only about Example 2. A econd important application of the nonuniform FFT i to magnetic reonance imaging (MRI) [9]. Under the aumption of a perfectly homogeneou magnetic field, the MRI hardware i able to acquire the Fourier tranform of a particular tiue property at elected point in the frequency domain. The ignal produced during the readout phae at time t i given by Z ðtþ ¼ qðxþe ikðtþx dx where x =(x 1, x 2 ) i a point in the two-dimenional image plane. In other word, (t) i preciely the value of the Fourier tranform ^q at the location k(t) =(k 1 (t), k 2 (t)). In the preence of a (known or unknown) field inhomogeneity given by /(x), however, we have Z ðtþ ¼ qðxþe ikðtþx e i/ðxþt dx: ð11þ Suppoe now that one eek to model the ignal (t) from known function of pace q(x) and /(x). Thi require the computation of
5 J.-Y. Lee, L. Greengard / Journal of Computational Phyic 206 (2005) ðt j Þ XN w n qðx n Þe ikðtjþxn e i/ðxnþtj ¼ XN w n qðx n Þe ikjxn where K j =(k 1 (t j ), k 2 (t j ), t j ) and X n ¼ðx 1 n x2 n /ðx nþþ. Thu, by embedding the data point in a higher dimenional pace, one can carry out the tranformation uing a type 3 NUFFT. We cite one timing reult, for a cae where the (x 1, x 2 ) variable are dicretized uing a regular meh on the unit box, and k-pace i travered by a piral trajectory imilar to that hown in Fig. 1(a) up to a maximum frequency of 60. The time interval of the readout phae wa caled i to [0,1], and /(x) wa allowed to vary in the interval [ 5p,5p]. Under thee condition, direct evaluation required about 140, while the NUFFT required 3.8 to obtain ix digit of accuracy. 4. Concluion We have preented a imple verion of the type 3 nonuniform FFT. It can be ued to approximate the continuou Fourier tranform when neither the patial nor the Fourier domain pacing i regular. Further, it allow for the evaluation of more general integral operator uch a (11) by embedding them in a higher dimenional pace. Acknowledgement Thi work wa upported by the Applied Mathematical Science Program of the US Department of Energy under Contract DEFGO200ER Reference [1] G. Beylkin, On the fat Fourier tranform of function with ingularitie, Appl. Comput. Harmonic Anal. 2 (1995) [2] A. Dutt, V. Rokhlin, Fat Fourier tramform for nonequipaced data, SIAM J. Sci. Comput. 14 (1993) [3] J.A. Feler, B.P. Sutton, Nonuniform Fat Fourier tranform uing min max interpolation, IEEE Tran. Signal Proc. 51 (2003) [4] K. Fourmont, Non-equipaced Fat Fourier tranform with application to tomography, J. Fourier Anal. Appl. 9 (2003) [5] L. Greengard, J.-Y. Lee, Accelerating the nonuniform fat fourier tranform, SIAM Rev. 46 (2004) [6] L. Greengard, P. Lin, Spectral approximation of the free-pace heat kernel, Appl. Comput. Harmonic Anal. 9 (2000) [7] A. Nielony, G. Steidl, Approximate factorization of Fourier matrice with nonequipaced knot, Lin. Alg. Appl. 366 (2003) [8] D. Pott, G. Steidl, M. Tache, Fat Fourier Tranform for Nonequipaced Data: A Tutorial, Modern Sampling Theory: Mathematic and Application, in: J.J. Benedetto, P. Ferreira (Ed.), Birkhauer Boton, Applied and Numerical Harmonic Analyi Serie, 2001, Chapter 12, pp [9] B.P. Sutton, D.C. Noll, J.A. Feler, Fat iterative recontruction for MRI in the preence of field inhomogeneitie, IEEE Tran. Med. Imag 22 (2003)
Determination of the local contrast of interference fringe patterns using continuous wavelet transform
Determination of the local contrat of interference fringe pattern uing continuou wavelet tranform Jong Kwang Hyok, Kim Chol Su Intitute of Optic, Department of Phyic, Kim Il Sung Univerity, Pyongyang,
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 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 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 informationBUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL
J. Autral. Math. Soc. Ser. B 4(999), 332 349 BUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL J. LEE and J.-M. VANDEN-BROECK 2 (Received 22 April 995; revied 23 April 996)
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 informationCHAPTER 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 informationThe continuous time random walk (CTRW) was introduced by Montroll and Weiss 1.
1 I. CONTINUOUS TIME RANDOM WALK The continuou time random walk (CTRW) wa introduced by Montroll and Wei 1. Unlike dicrete time random walk treated o far, in the CTRW the number of jump n made by the walker
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 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 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 informationComputers and Mathematics with Applications. Sharp algebraic periodicity conditions for linear higher order
Computer and Mathematic with Application 64 (2012) 2262 2274 Content lit available at SciVere ScienceDirect Computer and Mathematic with Application journal homepage: wwweleviercom/locate/camwa Sharp algebraic
More informationEE 477 Digital Signal Processing. 4 Sampling; Discrete-Time
EE 477 Digital Signal Proceing 4 Sampling; Dicrete-Time Sampling a Continuou Signal Obtain a equence of ignal ample uing a periodic intantaneou ampler: x [ n] = x( nt ) Often plot dicrete ignal a dot or
More informationChapter 4: Applications of Fourier Representations. Chih-Wei Liu
Chapter 4: Application of Fourier Repreentation Chih-Wei Liu Outline Introduction Fourier ranform of Periodic Signal Convolution/Multiplication with Non-Periodic Signal Fourier ranform of Dicrete-ime Signal
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 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 informationA SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES. Sanghyun Cho
A SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES Sanghyun Cho Abtract. We prove a implified verion of the Nah-Moer implicit function theorem in weighted Banach pace. We relax the
More informationWavelet Analysis in EPR Spectroscopy
Vol. 108 (2005) ACTA PHYSICA POLONICA A No. 1 Proceeding of the XXI International Meeting on Radio and Microwave Spectrocopy RAMIS 2005, Poznań-Bȩdlewo, Poland, April 24 28, 2005 Wavelet Analyi in EPR
More informationAdvanced Digital Signal Processing. Stationary/nonstationary signals. Time-Frequency Analysis... Some nonstationary signals. Time-Frequency Analysis
Advanced Digital ignal Proceing Prof. Nizamettin AYDIN naydin@yildiz.edu.tr Time-Frequency Analyi http://www.yildiz.edu.tr/~naydin 2 tationary/nontationary ignal Time-Frequency Analyi Fourier Tranform
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 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 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 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 informationFactor Analysis with Poisson Output
Factor Analyi with Poion Output Gopal Santhanam Byron Yu Krihna V. Shenoy, Department of Electrical Engineering, Neurocience Program Stanford Univerity Stanford, CA 94305, USA {gopal,byronyu,henoy}@tanford.edu
More informationChapter 1 Basic Description of Laser Diode Dynamics by Spatially Averaged Rate Equations: Conditions of Validity
Chapter 1 Baic Decription of Laer Diode Dynamic by Spatially Averaged Rate Equation: Condition of Validity A laer diode i a device in which an electric current input i converted to an output of photon.
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 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 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 informationTuning of High-Power Antenna Resonances by Appropriately Reactive Sources
Senor and Simulation Note Note 50 Augut 005 Tuning of High-Power Antenna Reonance by Appropriately Reactive Source Carl E. Baum Univerity of New Mexico Department of Electrical and Computer Engineering
More informationChapter 2: Problem Solutions
Chapter 2: Solution Dicrete Time Proceing of Continuou Time Signal Sampling à 2.. : Conider a inuoidal ignal and let u ample it at a frequency F 2kHz. xt 3co000t 0. a) Determine and expreion for the ampled
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 informationSECTION x2 x > 0, t > 0, (8.19a)
SECTION 8.5 433 8.5 Application of aplace Tranform to Partial Differential Equation In Section 8.2 and 8.3 we illutrated the effective ue of aplace tranform in olving ordinary differential equation. The
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 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 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 informationLecture 2: The z-transform
5-59- Control Sytem II FS 28 Lecture 2: The -Tranform From the Laplace Tranform to the tranform The Laplace tranform i an integral tranform which take a function of a real variable t to a function of a
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 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 informationUnbounded solutions of second order discrete BVPs on infinite intervals
Available online at www.tjna.com J. Nonlinear Sci. Appl. 9 206), 357 369 Reearch Article Unbounded olution of econd order dicrete BVP on infinite interval Hairong Lian a,, Jingwu Li a, Ravi P Agarwal b
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 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 informationCoordinate independence of quantum-mechanical q, qq. path integrals. H. Kleinert ), A. Chervyakov Introduction
vvv Phyic Letter A 10045 000 xxx www.elevier.nlrlocaterpla Coordinate independence of quantum-mechanical q, qq path integral. Kleinert ), A. Chervyakov 1 Freie UniÕeritat Berlin, Intitut fur Theoretiche
More informationSOLUTIONS TO ALGEBRAIC GEOMETRY AND ARITHMETIC CURVES BY QING LIU. I will collect my solutions to some of the exercises in this book in this document.
SOLUTIONS TO ALGEBRAIC GEOMETRY AND ARITHMETIC CURVES BY QING LIU CİHAN BAHRAN I will collect my olution to ome of the exercie in thi book in thi document. Section 2.1 1. Let A = k[[t ]] be the ring of
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 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 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 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 informationOnline supplementary information
Electronic Supplementary Material (ESI) for Soft Matter. Thi journal i The Royal Society of Chemitry 15 Online upplementary information Governing Equation For the vicou flow, we aume that the liquid thickne
More informationarxiv: v1 [math.ca] 23 Sep 2017
arxiv:709.08048v [math.ca] 3 Sep 07 On the unit ditance problem A. Ioevich Abtract. The Erdő unit ditance conjecture in the plane ay that the number of pair of point from a point et of ize n eparated by
More informationA Short Note on Hysteresis and Odd Harmonics
1 A Short Note on yterei and Odd armonic Peter R aoput Senior Reearch Scientit ubocope Pipeline Engineering outon, X 7751 USA pmaoput@varco.com June 19, Abtract hi hort note deal with the exitence of odd
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 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 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 informationSTOCHASTIC GENERALIZED TRANSPORTATION PROBLEM WITH DISCRETE DISTRIBUTION OF DEMAND
OPERATIONS RESEARCH AND DECISIONS No. 4 203 DOI: 0.5277/ord30402 Marcin ANHOLCER STOCHASTIC GENERALIZED TRANSPORTATION PROBLEM WITH DISCRETE DISTRIBUTION OF DEMAND The generalized tranportation problem
More informationTheoretical Computer Science. Optimal algorithms for online scheduling with bounded rearrangement at the end
Theoretical Computer Science 4 (0) 669 678 Content lit available at SciVere ScienceDirect Theoretical Computer Science journal homepage: www.elevier.com/locate/tc Optimal algorithm for online cheduling
More informationFAST METHOD FOR HIGH-FREQUENCY ACOUSTIC SCATTERING FROM RANDOM SCATTERERS
International Journal for Uncertainty Quantification, 1 (2): 99 117 (2011) FAST METHOD FOR HIGH-FREQUENCY ACOUSTIC SCATTERING FROM RANDOM SCATTERERS Paul Tuji, 1, Dongbin Xiu, 2 & Lexing Ying 3 1 ICES,
More information0 of the same magnitude. If we don t use an OA and ignore any damping, the CTF is
1 4. Image Simulation Influence of C Spherical aberration break the ymmetry that would otherwie exit between overfocu and underfocu. One reult i that the fringe in the FT of the CTF are generally farther
More informationarxiv: v1 [math.mg] 25 Aug 2011
ABSORBING ANGLES, STEINER MINIMAL TREES, AND ANTIPODALITY HORST MARTINI, KONRAD J. SWANEPOEL, AND P. OLOFF DE WET arxiv:08.5046v [math.mg] 25 Aug 20 Abtract. We give a new proof that a tar {op i : i =,...,
More informationModeling the scalar wave equation with Nyström methods
GEOPHYSICS, VOL. 71, NO. 5 SEPTEMBER-OCTOBER 200 ; P. T151 T158, 7 FIGS. 10.1190/1.2335505 Modeling the calar wave equation with Nytröm method Jing-Bo Chen 1 ABSTRACT High-accuracy numerical cheme for
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 informationSimple Observer Based Synchronization of Lorenz System with Parametric Uncertainty
IOSR Journal of Electrical and Electronic Engineering (IOSR-JEEE) ISSN: 78-676Volume, Iue 6 (Nov. - Dec. 0), PP 4-0 Simple Oberver Baed Synchronization of Lorenz Sytem with Parametric Uncertainty Manih
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 informationMolecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions
Original Paper orma, 5, 9 7, Molecular Dynamic Simulation of Nonequilibrium Effect ociated with Thermally ctivated Exothermic Reaction Jerzy GORECKI and Joanna Natalia GORECK Intitute of Phyical Chemitry,
More informationNCAAPMT Calculus Challenge Challenge #3 Due: October 26, 2011
NCAAPMT Calculu Challenge 011 01 Challenge #3 Due: October 6, 011 A Model of Traffic Flow Everyone ha at ome time been on a multi-lane highway and encountered road contruction that required the traffic
More informationPDF hosted at the Radboud Repository of the Radboud University Nijmegen
PDF hoted at the Radboud Repoitory of the Radboud Univerity Nijmegen The following full text i an author' verion which may differ from the publiher' verion. For additional information about thi publication
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 informationTRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL
GLASNIK MATEMATIČKI Vol. 38583, 73 84 TRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL p-laplacian Haihen Lü, Donal O Regan and Ravi P. Agarwal Academy of Mathematic and Sytem Science, Beijing, China, National
More informationNew bounds for Morse clusters
New bound for More cluter Tamá Vinkó Advanced Concept Team, European Space Agency, ESTEC Keplerlaan 1, 2201 AZ Noordwijk, The Netherland Tama.Vinko@ea.int and Arnold Neumaier Fakultät für Mathematik, Univerität
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 informationLTV System Modelling
Helinki Univerit of Technolog S-72.333 Potgraduate Coure in Radiocommunication Fall 2000 LTV Stem Modelling Heikki Lorentz Sonera Entrum O heikki.lorentz@onera.fi Januar 23 rd 200 Content. Introduction
More informationCompact finite-difference approximations for anisotropic image smoothing and painting
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
More informationSUPPLEMENTARY INFORMATION
DOI: 10.1038/NPHOTON.014.108 Supplementary Information "Spin angular momentum and tunable polarization in high harmonic generation" Avner Fleicher, Ofer Kfir, Tzvi Dikin, Pavel Sidorenko, and Oren Cohen
More informationLecture 23 Date:
Lecture 3 Date: 4.4.16 Plane Wave in Free Space and Good Conductor Power and Poynting Vector Wave Propagation in Loy Dielectric Wave propagating in z-direction and having only x-component i given by: E
More informationCOPRIME ARRAYS AND SAMPLERS FOR SPACE-TIME ADAPTIVE PROCESSING
COPRIE ARRAYS AND SAPLERS FOR SPACE-IE ADAPIVE PROCESSING Chun-Lin Liu 1 and P. P. Vaidyanathan 2 Dept. of Electrical Engineering, 136-93 California Intitute of echnology, Paadena, CA 91125, USA E-mail:
More informationA note on the bounds of the error of Gauss Turán-type quadratures
Journal of Computational and Applied Mathematic 2 27 276 282 www.elevier.com/locate/cam A note on the bound of the error of Gau Turán-type quadrature Gradimir V. Milovanović a, Miodrag M. Spalević b, a
More informationFourier-Conjugate Models in the Corpuscular-Wave Dualism Concept
International Journal of Adanced Reearch in Phyical Science (IJARPS) Volume, Iue 0, October 05, PP 6-30 ISSN 349-7874 (Print) & ISSN 349-788 (Online) www.arcjournal.org Fourier-Conjugate Model in the Corpucular-Wae
More informationEvolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis
Proceeding of 01 4th International Conference on Machine Learning and Computing IPCSIT vol. 5 (01) (01) IACSIT Pre, Singapore Evolutionary Algorithm Baed Fixed Order Robut Controller Deign and Robutne
More informationOn the Isomorphism of Fractional Factorial Designs 1
journal of complexity 17, 8697 (2001) doi:10.1006jcom.2000.0569, available online at http:www.idealibrary.com on On the Iomorphim of Fractional Factorial Deign 1 Chang-Xing Ma Department of Statitic, Nankai
More informationLecture 8: Period Finding: Simon s Problem over Z N
Quantum Computation (CMU 8-859BB, Fall 205) Lecture 8: Period Finding: Simon Problem over Z October 5, 205 Lecturer: John Wright Scribe: icola Rech Problem A mentioned previouly, period finding i a rephraing
More informationObserving Condensations in Atomic Fermi Gases
Oberving Condenation in Atomic Fermi Gae (Term Eay for 498ESM, Spring 2004) Ruqing Xu Department of Phyic, UIUC (May 6, 2004) Abtract Oberving condenation in a ga of fermion ha been another intereting
More informationDynamic response of a double Euler Bernoulli beam due to a moving constant load
Journal of Sound and Vibration 297 (26) 477 49 JOURNAL OF SOUND AND VIBRATION www.elevier.com/locate/jvi Dynamic repone of a double Euler Bernoulli beam due to a moving contant load M. Abu-Hilal Department
More informationin a circular cylindrical cavity K. Kakazu Department of Physics, University of the Ryukyus, Okinawa , Japan Y. S. Kim
Quantization of electromagnetic eld in a circular cylindrical cavity K. Kakazu Department of Phyic, Univerity of the Ryukyu, Okinawa 903-0, Japan Y. S. Kim Department of Phyic, Univerity of Maryland, College
More informationCalculation of the temperature of boundary layer beside wall with time-dependent heat transfer coefficient
Ŕ periodica polytechnica Mechanical Engineering 54/1 21 15 2 doi: 1.3311/pp.me.21-1.3 web: http:// www.pp.bme.hu/ me c Periodica Polytechnica 21 RESERCH RTICLE Calculation of the temperature of boundary
More informationKERNEL-RESOLVENT RELATIONS FOR AN INTEGRAL EQUATION
Ø Ñ Å Ø Ñ Ø Ð ÈÙ Ð Ø ÓÒ DOI: 1.478/v117-11-3-7 Tatra Mt. Math. Publ. 48 (11), 5 4 KERNEL-RESOLVENT RELATIONS FOR AN INTEGRAL EQUATION Theodore A. Burton ABSTRACT. Weconideracalarintegralequationz(t) =
More informationPacific Journal of Mathematics
Pacific Journal of Mathematic OSCILLAION AND NONOSCILLAION OF FORCED SECOND ORDER DYNAMIC EQUAIONS MARIN BOHNER AND CHRISOPHER C. ISDELL Volume 230 No. March 2007 PACIFIC JOURNAL OF MAHEMAICS Vol. 230,
More informationTrajectory Planning and Feedforward Design for High Performance Motion Systems
Trajectory Planning and Feedforward Deign for High Performance Motion Sytem Paul Lambrecht, Matthij Boerlage, Maarten Steinbuch Faculty of Mechanical Engineering, Control Sytem Technology Group Eindhoven
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 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 informationSupplementary Figures
Supplementary Figure Supplementary Figure S1: Extraction of the SOF. The tandard deviation of meaured V xy at aturated tate (between 2.4 ka/m and 12 ka/m), V 2 d Vxy( H, j, hm ) Vxy( H, j, hm ) 2. The
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 information4.6 Principal trajectories in terms of amplitude and phase function
4.6 Principal trajectorie in term of amplitude and phae function We denote with C() and S() the coinelike and inelike trajectorie relative to the tart point = : C( ) = S( ) = C( ) = S( ) = Both can be
More informationRandom Sparse Linear Systems Observed Via Arbitrary Channels: A Decoupling Principle
Random Spare Linear Sytem Oberved Via Arbitrary Channel: A Decoupling Principle Dongning Guo Department of Electrical Engineering & Computer Science Northwetern Univerity Evanton, IL 6008, USA. Chih-Chun
More informationSAMPLING. Sampling is the acquisition of a continuous signal at discrete time intervals and is a fundamental concept in real-time signal processing.
SAMPLING Sampling i the acquiition of a continuou ignal at dicrete time interval and i a fundamental concept in real-time ignal proceing. he actual ampling operation can alo be defined by the figure belo
More informationRandom vs. Deterministic Deployment of Sensors in the Presence of Failures and Placement Errors
Random v. Determinitic Deployment of Senor in the Preence of Failure and Placement Error Paul Baliter Univerity of Memphi pbalitr@memphi.edu Santoh Kumar Univerity of Memphi antoh.kumar@memphi.edu Abtract
More informationModule 4: Time Response of discrete time systems Lecture Note 1
Digital Control Module 4 Lecture Module 4: ime Repone of dicrete time ytem Lecture Note ime Repone of dicrete time ytem Abolute tability i a baic requirement of all control ytem. Apart from that, good
More informationCHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS
Copyright 22 IFAC 5th Triennial World Congre, Barcelona, Spain CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS Tritan Pérez Graham C. Goodwin Maria M. Serón Department of Electrical
More informationApproximate Analytical Solution for Quadratic Riccati Differential Equation
Iranian J. of Numerical Analyi and Optimization Vol 3, No. 2, 2013), pp 21-31 Approximate Analytical Solution for Quadratic Riccati Differential Equation H. Aminikhah Abtract In thi paper, we introduce
More informationDimensional Analysis A Tool for Guiding Mathematical Calculations
Dimenional Analyi A Tool for Guiding Mathematical Calculation Dougla A. Kerr Iue 1 February 6, 2010 ABSTRACT AND INTRODUCTION In converting quantitie from one unit to another, we may know the applicable
More informationFinite Element Analysis of a Fiber Bragg Grating Accelerometer for Performance Optimization
Finite Element Analyi of a Fiber Bragg Grating Accelerometer for Performance Optimization N. Baumallick*, P. Biwa, K. Dagupta and S. Bandyopadhyay Fiber Optic Laboratory, Central Gla and Ceramic Reearch
More informationThe Influence of the Load Condition upon the Radial Distribution of Electromagnetic Vibration and Noise in a Three-Phase Squirrel-Cage Induction Motor
The Influence of the Load Condition upon the Radial Ditribution of Electromagnetic Vibration and Noie in a Three-Phae Squirrel-Cage Induction Motor Yuta Sato 1, Iao Hirotuka 1, Kazuo Tuboi 1, Maanori Nakamura
More information