Spectrum of Fractal Interpolation Functions
|
|
- Jade McDowell
- 6 years ago
- Views:
Transcription
1 Spectrum of Fractal Interpolation Functions ikolaos Vasiloglou, Member, IEEE and Petros Maragos, Fellow, IEEE June 8, Abstract In this paper we compute the Fourier spectrum of the Fractal Interpolation Functions FIFs as introduced by Michael Barnsley We show that there is an analytical way to compute them In this paper we attempt to solve the inverse problem of FIF by using the spectrum results in a sequence of a sets that converges in the attractor of the operator A, which satisfies the condition A = W(A) = n= w n(a) Any system that uses the Hutchinson operator in order to generate iteratively the attractor A is called Iterated Function System (IFS) Iterated Function Systems Original Translation & Scaling The affine transform performs translation stretching and rotation on a given set In the special case of two dimensions the affine transform on a set S in -D space is described by the equation: [ ][ ] [ ] a b x e w(x, y) = + c d y f where (x, y) S The effects of an affine transform on a set are depicted in fig The union of affine transformations is called the Hutchinson operator: W = n= w n For a specified metric the distance h(a, B) between two sets A, B can be defined Under certain conditions [] the Hutchinson operator is contractive, h(w(a), W(B)) sh(a, B), s < Successive iterations with Hutchinson operator on a random set Submitted: Sep This research work was supported by the Greek Secretariat for Research and Technology and by the European Union under the program EΠET-98 with Grant # 98ΓT6, when both authors where with with the Department of Electrical & Computer Engineering, ational Technical University of Athens ikolaos Vasiloglou is now with Analytics 5 LLC Atlanta Georgia USA Petros Maragos is with the Department of Electrical & Computer Engineering, ational Technical University of Athens, Zografou, 577 Athens, Greece nvasil@ieeeorg,maragos@csntuagr Translation & Rotation Scaling & Reflection Figure : An affine transform can translate rotate and flip a -dimensional shape Fractal Interpolation Functions Fractal Interpolation Functions (FIF) is a special case of the -dimensional IFS and maintain all their characteristics FIF attractors are continuous functions that can be used to model continuous signals FIF interpolate a given set of + points (x n, y n ), n =, x < x < x < x < < x
2 The FIF that interpolates the above set is comprised of affine maps: [ ] [ ][ ] [ ] x an x en w n = +, n =,, y d n y f n c n satisfies the conditions of (), has a unique fixed point, the function F C, (TF)(x) = F(x), x [x, x ] The necessary condition is that the Interpolation Function passes from the + initial points, [ ] [ ] x xn w n = y y n [ ] [ ] x xn and w n = y y n n =,, () We call order of the FIF, the number of the affine maps The conditions provide 4 equations for 5 parameters, so d n, the vertical scaling factor is chosen to be the free parameter If we solve the above equations for a n, c n, e n, f n in terms of d n, we find :, The following restrictions guarantee the contractivity of T operator: a n < d n < Barnsley s operator is very useful because the attractor of an FIF can be generated with the iterative application of T on an initial signal F() : F m+ = TF m, F = lim m T m (F ), F C a n = x n x n x x, () e n = x x n x x n x x, () c n = y n y n x x d n(y y ) x x, (4) f n = x y n x y n x x d n x y x y x x (5) Let the real numbers a n, c n, e n, f n be defined by (-5) Barnsley [] introduced the operator T for the class C of continuous functions, T : C C by F (x) F (x) Attractor and (TF)(x) = c n l n (x) + d n F(l n (x)) + f n x [x n, x n ], n =,,, l n : [x, x ] [x n, x n ] the invertible transformation The above operator: l n (x) = a n x + e n is a contraction mapping according to Hausdorff metrics, Figure : Formation of the FIF attractor after successive iterations Discrete FIF One of the most interesting properties of the FIF is that its attractor is independent of the initial signal (initiator) If the initiator is a continuous signal, for example the linear interpolation between the given interpolation points, all the instances throughout all the steps of the iterations will be continuous signals, fig On the other hand if the initiator is a single point the attractor formed after infinite iterations
3 will be a continuous signal, but signals instances of the FIF throughout the iterations will be discrete signals, fig Although during all the iterations the instances are discrete signals or strictly mathematically speaking finite countable sets, the attractor is a continuous signal, an infinite and uncountable set Instances of a the formation of an FIF from a discrete initiator are shown in fig Figure : Generation of a discrete FIF It is essential to show that a good choice of the initiator is the + given interpolation points (x n, y n ), n =, After the first iteration in each of the subintervals between the points, new points are generated fig Let these points be (x s, y s ), s =,,( ) otice that all these ++( ) = + points belong to the attractor of FIF This wouldn t be true if the initiator was a different set apart from the given interpolation points By repeating this procedure after m iterations we get an M(m)-point discrete sequence that is a sampling of the FIF s attractor, with sampling period T s = M If the initiator included any other irrelevant point, that would not be mapped to an attractor s point after a finite number of iterations It can be proved that if the initiator is not the set of the initial interpolation points, then the error of the formed sequence after the mth iteration from the attractor decreases and goes to zero as m The number of the attractor s samples is : ( ) = ( ) + ( ) = + M = M(m) = m + (6) The discrete signal formed after the mth iteration over the discrete initiator, is called discrete FIF : f[n] = F(x n ), n =,,,M (7) Before expanding Barnsley operator T introduced above for the discrete FIF it is necessary to make clear that according to the strict definition of Fractals the discrete FIFs are not fractal sets because they are finite Discrete FIFs must be considered as approximations of the continuous Since the Barnsley operator assumes infinite resolution, it doesn t apply in discrete signals So for discrete FIF the following modified operator is used: T f[n] = (d k f[ n e k ] +c k ( n e k a a ) + f n) (8) k= (u(n e k ) u(n e k )) The Symbol f[n] defines that within two successive samples of f[n], zeros have been interpolated, fig 4 The term f n denotes the parameter of FIF as defined in (5) and should not be confused with the discrete signal f[n] Similarly to the continuous case, a discrete FIF of m iterations can be constructed with the following procedure: f m+ [n] = T f m [n] (9) It is obvious that when m the discrete FIF becomes continuous otice that T depends on the number of iterations m More specifically the parameters of FIF as defined in (-5) depend on the x n The x n change value because of the upsampling
4 F = F = and the points are evenly spaced Then(-5) becomes, a n =, () e n = n, () c n = F n F n, () f n = F n () 5 The above equations show that the FIF parameters are decoupled from each other This is very impor- because they can be estimated independently tant Moreover it is clear that if a n parameter is estimated Figure 4: Barnsley operator for discrete signals The then all parameters can be found directly, except for initiator is upsampled and interpolated according to d n (9 Any FIF can be transformed to an equivalent FIF that satisfies the first two conditions without loosing its fundamental properties More specifically fig 6 The generation of an FIF can be represented in terms shows how a given FIF can be transformed so as to of a linear system as shown in fig 5 satisfy the conditions for the first and the last point It is convenient to use an auxiliary affine map w aux that will rotate, scale and translate the given one ĝ(ω) ˆf(ω) The w aux transform is invertible and does not affect the intrinsic parameters a n, d n ˆq(ω) 5 5 Figure 5: Block diagram of the FIF Genaration 5 5 Computation of FIF s Fourier Spectrum In order to simplify (-5), it is very convenient to adopt the following assumptions x =, x =, Figure 6: An example of FIF out of range [, ] Applying an affine transform we can tie it at points (, ) and (, ) 4
5 Let A = {(x, y) : F(x) = y}, be the original attractor and w n(a) = A (4) n= where w n, [ ] [ ] [ ] [ ] x an x en w n = +, n =,, G(x) = y d n y f n c n and u(x) = { x < x G(x) is the piecewise linear function between the interpolation points : n= c n ( x e n a +f n )(u(x e n ) u(x e n )) The new transformed attractor is A = {(x, y ), F (x) = y } The new points are connected with the initial [ x y ] = [ ] [ aaux x d aux y c aux ] [ ] eaux + f aux That means A = w aux A These new points belong to a new FIF Setting in (4) A = waux A and applying the map w aux, with w auxw n w n= aux (A ) = A, (5) [ waux = /a aux caux a auxd aux /d aux ] By setting f aux = e aux =, d aux =, the affine maps of the new FIF are : w n [ ] [ x y = ] a n (c n d n c aux + c aux a n )/a aux d n [ ] [ ] x a y + aux e n c aux e n + f n otice that the new FIF has the same order and the same d n parameters Spectrum of Continuous FIF The application of the above simplifications to the Barnsley operator T results in the following equation for the FIF: FOURIER SPECTRUM FREQUECY (Hz) Figure 7: An example of Q(Ω) Applying continuous fourier transform [5]: F(Ω) = + F(x) F(Ω) F(x)e iωx dx F(Ω) = G(Ω) + af(aω) d n e iωen, (7) n= G(Ω) = n n= ( e iω e iω n (iω) ( cn a + iω(f n cnen a )) + cn a n n e iω n n e iω iω ) We define the function : (8) F(x) = n= (d nf( x en a )+ +c n ( x en a ) + f n )(u(x e n ) u(x e n )) (6) Q(Ω) = a n= (n ) iω d n e = a d n+ e iωn (9) n= 5
6 otice that the above function is the discrete time fourier transform of the discrete sequence {d, d, d }, so we deduce that Q(Ω kπ) = Q((k + )π Ω), k =,, (fig 7) FOURIER SPECTRUM (db) FOURIER SPECTRUM (db) FOURIER SPECTRUM (db) FREQUECY (Hz) FREQUECY (Hz) FREQUECY (Hz) Figure 8: Left column: Fractal interpolation between points (, ), (5, ), (5, 4), (75, 5), (, ) after,,5 iterations right column: Corresponding spectrums The Fourier spectrum satisfies the following equation: F(Ω) = Q(Ω)F(aΩ) + G(Ω) () Through Barsley Operator in frequency domain the Fourier spectrum can be computed iteratively, fig 8 F m+ (Ω) = Q(Ω)F m (aω) + G(Ω) () F (Ω) = After infinite iterations: i F(Ω) = G(Ωa i ) Q(Ωa j ) () i= j= Spectrum of Discrete FIF From the () it is obvious that the FIF signals are not band-limited As a result the spectrum of a discrete FIF is aliased Assume T s = M has been chosen as the sampling period, then the Discrete Time Fourier Transform (DTFT) [5] is : ˆf(ω) = M f[n]e iωn, ω = ΩT s n= As in continuous case for the computation of the DTFT the T operator is used in frequency domain: ˆf m+ (ω) = ˆf m (aω) d k ae iωep + ĝ m (ω) () p= Posing the analogy between the continuous and the discrete case we define the q function : ˆq(ω) = ˆq m (ω) = a d p e iωep p= d p e iωep, p= and e p = p (M ), so ˆq(ω) = p= d p+e iωp (M ) p =,,, = p= d p+e ( iωp)(m m i) i=, p =,,, (4) ˆf m+ (ω) = ˆq m (ω) ˆf m (aω) + ĝ m (ω) (5) After m iterations m i ˆf m (ω) = ĝ m i (ωa i ) ˆq i j+ (ωa j ) (6) i= j= The Discrete Fourier Transform can be computed after keeping the frequency values between and π and sampling the spectrum at ω =, π M, π M,, (M )π M, fig 9 In the continuous case Q(Ω) has period π In the discrete case where ω [, π] the ˆq(ω) function has period π M 6
7 FOURIER SPECTRUM (db) FREQUECY (Hz) Figure 9: 5-point FIF attractor with d n = 74, 8, 77, 85, 88 and the spectrum after the 5th iteration of linear equations ˆq() ˆq[ω ] e iω e ( )iω ˆq[ω ] = e iω e ( )iω ˆq[ω ] 5 5 e iω R e ( )iω R d d d d (7) FIF parameter estimation using spectrum FIF modelling of signals has been proposed by Mazel [], using information of signal in time domain In this section the spectrum of signal is used to estimate its parameters, provided the FIF order is known Given the signal f[n], n =,,,M, let be FIF s estimated order We expect its spectrum to satisfy the following equation: (b) (a) (c) ˆq[ω k ] = ˆf m+ [ω k ] ĝ[ω k ], ˆf m [aω k ] (d) (e) aω k =, π M, π π,, (M ) M M All values that zero denominator are excluded otice that f m+ [n] and f m [n] signals must have the same length Considering also that f m+ [] = f m [] and f m+ [M(m + )] = f m [M(m)], where M(m), M(m + ) their lengths In order to equalize their lengths it is necessary to interpolate with zeros f m [n] Because of the term ˆf m [aω k ] the computation of FFT must be on kπ M, k =,,M This is done by padding the signal with M M zeros The d n are determined by the solution of the system Figure : (a) FIF without noise In (b),(c) FIF contaminated with noise SR = and SR = In (d),(e) reconstructed FIF 4 Results The above algorithm was tested in FIF signals contaminated with white noise and the parameters extracted were very close to the real, fig The main advantage of the method is the use of FFT It is well known that FFT is quite simple and easily implemented The first disadvantage of the method is that 7
8 it cannot find the FIF order In the above experiments we tried for =,, until the estimated ˆq function satisfied the periodicity and symmetry conditions mentioned earlier The second and most important problem is that it is very sensitive in the window effect of the FFT It is known that although a part of FIF signal is self affine it is not an FIF Trying to model it as an FIF results in wrong estimations In the example of fig, it is evident that the period of the ˆq function has been expanded Having chosen the order, the period of ˆq is known π M In the right plot of fig it is evident that the period of ˆq is much higher than the expected But although the spectrum is not the best method for solving the inverse problem it can be used as powerful analysis tool As shown above it can reveal the FIF nature of a signal and it can also help in the prediction of a missing part of a time series, given that it belongs to class of FIF [4] H DI Abarbanel, Analysis of Observed Chaotic Data, Springer Verlag, ew York 996 [5] A VOppenheim, R W Schafer, Deiscrete-Time Signal Processing, Prentice Hall, FOURIER TRASFORM 5 5 FOURIER TRASFORM ORMALIZED FREQUECY (ω/π) x ORMALIZED FREQUECY (ω/π) x 4 Figure : Left : ˆq(ω) Right : ˆq(ω) of the corrupted FIF References [] B B Mandelbrot, The Fractal Geometry of ature WH Freeman and Company, ew York, 98 [] M F Barnsley, Fractals Everywhere, Academic Press, 99 [] D S Mazel, M H Hayes, Using Iterated Function Systems to Model Discrete Sequences, IEEE Transactions on Signal Processing, Vol 4 (7), pp 74-74, July 99 8
9
10 FOURIER SPECTRUM (db) FREQUECY (Hz)
11 TIME
Spectrum of Fractal Interpolation Functions
Spectrum of Frctl Interpoltion Functions ikolos Vsiloglou, Member, IEEE nd Petros Mrgos, Fellow, IEEE October, 8 Abstrct In this pper we compute the Fourier spectrum of the Frctl Interpoltion Functions
More informationIterated function systems on functions of bounded variation
Iterated function systems on functions of bounded variation D. La Torre, F. Mendivil and E.R. Vrscay December 18, 2015 Abstract We show that under certain hypotheses, an Iterated Function System on Mappings
More informationThe Fourier Transform Method
The Fourier Transform Method R. C. Trinity University Partial Differential Equations April 22, 2014 Recall The Fourier transform The Fourier transform of a piecewise smooth f L 1 (R) is ˆf(ω) = F(f)(ω)
More informationApproximation of the attractor of a countable iterated function system 1
General Mathematics Vol. 17, No. 3 (2009), 221 231 Approximation of the attractor of a countable iterated function system 1 Nicolae-Adrian Secelean Abstract In this paper we will describe a construction
More information1 From Fourier Series to Fourier transform
Differential Equations 2 Fall 206 The Fourier Transform From Fourier Series to Fourier transform Recall: If fx is defined and piecewise smooth in the interval [, ] then we can write where fx = a n = b
More informationTHE CONLEY ATTRACTORS OF AN ITERATED FUNCTION SYSTEM
Bull. Aust. Math. Soc. 88 (2013), 267 279 doi:10.1017/s0004972713000348 THE CONLEY ATTRACTORS OF AN ITERATED FUNCTION SYSTEM MICHAEL F. BARNSLEY and ANDREW VINCE (Received 15 August 2012; accepted 21 February
More informationTHE SHIFT SPACE FOR AN INFINITE ITERATED FUNCTION SYSTEM
THE SHIFT SPACE FOR AN INFINITE ITERATED FUNCTION SYSTEM ALEXANDRU MIHAIL and RADU MICULESCU The aim of the paper is to define the shift space for an infinite iterated function systems (IIFS) and to describe
More informationLecture 4: Fourier Transforms.
1 Definition. Lecture 4: Fourier Transforms. We now come to Fourier transforms, which we give in the form of a definition. First we define the spaces L 1 () and L 2 (). Definition 1.1 The space L 1 ()
More informationConstruction of Smooth Fractal Surfaces Using Hermite Fractal Interpolation Functions. P. Bouboulis, L. Dalla and M. Kostaki-Kosta
BULLETIN OF THE GREEK MATHEMATICAL SOCIETY Volume 54 27 (79 95) Construction of Smooth Fractal Surfaces Using Hermite Fractal Interpolation Functions P. Bouboulis L. Dalla and M. Kostaki-Kosta Received
More informationSelf-Referentiality, Fractels, and Applications
Self-Referentiality, Fractels, and Applications Peter Massopust Center of Mathematics Research Unit M15 Technical University of Munich Germany 2nd IM-Workshop on Applied Approximation, Signals, and Images,
More informationGENERALIZED CUBIC SPLINE FRACTAL INTERPOLATION FUNCTIONS
SIAM J NUMER ANAL Vol 44, No, pp 55 7 c Society for Industrial and Applied Mathematics GENERALIZED CUBIC SPLINE FRACTAL INTERPOLATION FUNCTIONS A K B CHAND AND G P KAPOOR Abstract We construct a generalized
More informationNOTES ON BARNSLEY FERN
NOTES ON BARNSLEY FERN ERIC MARTIN 1. Affine transformations An affine transformation on the plane is a mapping T that preserves collinearity and ratios of distances: given two points A and B, if C is
More informationFourier transforms. R. C. Daileda. Partial Differential Equations April 17, Trinity University
The Fourier Transform R. C. Trinity University Partial Differential Equations April 17, 214 The Fourier series representation For periodic functions Recall: If f is a 2p-periodic (piecewise smooth) function,
More informationFourier Series. Suppose f : [0, 2π] C and: f(x) = n=0. a n cos nx. How could we find the a n if we know f? Having a look at cos :
Fourier Series Suppose f : [0, 2π] C and: f(x) = n=0 a n cos nx How could we find the a n if we know f? Having a look at cos : cos(0x) cos(1x) cos(2x) Average 1 Average 0 Average 0 2π 0 f(x) dx = n=0 2π
More informationAnalytic Continuation of Analytic (Fractal) Functions
Analytic Continuation of Analytic (Fractal) Functions Michael F. Barnsley Andrew Vince (UFL) Australian National University 10 December 2012 Analytic continuations of fractals generalises analytic continuation
More informationHIDDEN VARIABLE BIVARIATE FRACTAL INTERPOLATION SURFACES
Fractals, Vol. 11, No. 3 (2003 277 288 c World Scientific Publishing Company HIDDEN VARIABLE BIVARIATE FRACTAL INTERPOLATION SURFACES A. K. B. CHAND and G. P. KAPOOR Department of Mathematics, Indian Institute
More informationFrom Fourier to Wavelets in 60 Slides
From Fourier to Wavelets in 60 Slides Bernhard G. Bodmann Math Department, UH September 20, 2008 B. G. Bodmann (UH Math) From Fourier to Wavelets in 60 Slides September 20, 2008 1 / 62 Outline 1 From Fourier
More informationToday s lecture. The Fourier transform. Sampling, aliasing, interpolation The Fast Fourier Transform (FFT) algorithm
Today s lecture The Fourier transform What is it? What is it useful for? What are its properties? Sampling, aliasing, interpolation The Fast Fourier Transform (FFT) algorithm Jean Baptiste Joseph Fourier
More informationarxiv: v2 [math.ca] 4 Jun 2017
EXCURSIONS ON CANTOR-LIKE SETS ROBERT DIMARTINO AND WILFREDO O. URBINA arxiv:4.70v [math.ca] 4 Jun 07 ABSTRACT. The ternary Cantor set C, constructed by George Cantor in 883, is probably the best known
More informationThree hours THE UNIVERSITY OF MANCHESTER. 24th January
Three hours MATH41011 THE UNIVERSITY OF MANCHESTER FOURIER ANALYSIS AND LEBESGUE INTEGRATION 24th January 2013 9.45 12.45 Answer ALL SIX questions in Section A (25 marks in total). Answer THREE of the
More informationFractals. R. J. Renka 11/14/2016. Department of Computer Science & Engineering University of North Texas. R. J. Renka Fractals
Fractals R. J. Renka Department of Computer Science & Engineering University of North Texas 11/14/2016 Introduction In graphics, fractals are used to produce natural scenes with irregular shapes, such
More informationFractals and Dimension
Chapter 7 Fractals and Dimension Dimension We say that a smooth curve has dimension 1, a plane has dimension 2 and so on, but it is not so obvious at first what dimension we should ascribe to the Sierpinski
More informationA Class of Möbius Iterated Function Systems
1 A Class of Möbius Iterated Function Systems arxiv:180.0143v1 [math.ds] 30 Jan 018 Gökçe ÇAKMAK, Ali DENİZ, Şahin KOÇAK Abstract We give a procedure to produce Möbius iterated function systems (MIFS)
More informationON THE DIAMETER OF THE ATTRACTOR OF AN IFS Serge Dubuc Raouf Hamzaoui Abstract We investigate methods for the evaluation of the diameter of the attrac
ON THE DIAMETER OF THE ATTRACTOR OF AN IFS Serge Dubuc Raouf Hamzaoui Abstract We investigate methods for the evaluation of the diameter of the attractor of an IFS. We propose an upper bound for the diameter
More information1 The Complex Fourier Series
The Complex Fourier Series Adding Complexity to life The complex Fourier series is in some ways a superior product, at least for those people who are not terrified by complex numbers. Suppose fx) is a
More informationTHE LIMIT PROCESS (AN INTUITIVE INTRODUCTION)
The Limit Process THE LIMIT PROCESS (AN INTUITIVE INTRODUCTION) We could begin by saying that limits are important in calculus, but that would be a major understatement. Without limits, calculus would
More informationLecture 34. Fourier Transforms
Lecture 34 Fourier Transforms In this section, we introduce the Fourier transform, a method of analyzing the frequency content of functions that are no longer τ-periodic, but which are defined over the
More informationFractals: A Mathematical Framework
Fractals: A Mathematical Framework John E Hutchinson Department of Mathematics School of Mathematical Sciences Australian National University (e-mail: JohnHutchinson@anueduau) Abstract We survey some of
More informationEXAMINATION: MATHEMATICAL TECHNIQUES FOR IMAGE ANALYSIS
EXAMINATION: MATHEMATICAL TECHNIQUES FOR IMAGE ANALYSIS Course code: 8D Date: Thursday April 8 th, Time: 4h 7h Place: AUD 3 Read this first! Write your name and student identification number on each paper
More informationITERATED FUNCTION SYSTEMS WITH A GIVEN CONTINUOUS STATIONARY DISTRIBUTION
Fractals, Vol., No. 3 (1) 1 6 c World Scientific Publishing Company DOI: 1.114/S18348X1517X Fractals Downloaded from www.worldscientific.com by UPPSALA UNIVERSITY ANGSTROM on 9/17/1. For personal use only.
More informationHABILITATION THESIS. New results in the theory of Countable Iterated Function Systems
HABILITATION THESIS New results in the theory of Countable Iterated Function Systems - abstract - Nicolae-Adrian Secelean Specialization: Mathematics Lucian Blaga University of Sibiu 2014 Abstract The
More informationRECURRENT ITERATED FUNCTION SYSTEMS
RECURRENT ITERATED FUNCTION SYSTEMS ALEXANDRU MIHAIL The theory of fractal sets is an old one, but it also is a modern domain of research. One of the main source of the development of the theory of fractal
More informationJim Lambers ENERGY 281 Spring Quarter Lecture 5 Notes
Jim ambers ENERGY 28 Spring Quarter 27-8 ecture 5 Notes These notes are based on Rosalind Archer s PE28 lecture notes, with some revisions by Jim ambers. Fourier Series Recall that in ecture 2, when we
More informationDiscrete Fourier Transform
Discrete Fourier Transform DD2423 Image Analysis and Computer Vision Mårten Björkman Computational Vision and Active Perception School of Computer Science and Communication November 13, 2013 Mårten Björkman
More informationMULTIPLIERS ON SPACES OF FUNCTIONS ON A LOCALLY COMPACT ABELIAN GROUP WITH VALUES IN A HILBERT SPACE. Violeta Petkova
Serdica Math. J. 32 (2006), 215 226 MULTIPLIERS ON SPACES OF FUNCTIONS ON A LOCALLY COMPACT ABELIAN ROUP WITH VALUES IN A HILBERT SPACE Violeta Petkova Communicated by S. L. Troyansky We prove a representation
More informationIEICE Transactions on Fundamentals of Electronics, Communications and Computer Sciences, Vol. E82-A, No. 8, pp , August, 1999.
5 PAPER Special Section on Digital Signal Processing An Encoding Algorithm for IFS Coding of Homogeneous Fractal Images Using Univariate Polynomial Manipulation Toshimizu ABIKO y, Student Member and Masayuki
More informationSelf-Simi lar Structure in Hilbert's Space-Fi l 1 ing Curve
40 MATHEMATICS MAGAZINE Self-Simi lar Structure in Hilbert's Space-Fi l 1 ing Curve MARK MCCLURE University of North Carolina at Asheville Asheville, North Carolina 28801 mcmcclur@bulldog.unca.edu Hilbert's
More informationGeneration of Fractals via Self-Similar Group of Kannan Iterated Function System
Appl. Math. Inf. Sci. 9, No. 6, 3245-3250 (2015) 3245 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.12785/amis/090652 Generation of Fractals via Self-Similar
More informationVALUATION THEORY, GENERALIZED IFS ATTRACTORS AND FRACTALS
VALUATION THEORY, GENERALIZED IFS ATTRACTORS AND FRACTALS JAN DOBROWOLSKI AND FRANZ-VIKTOR KUHLMANN Abstract. Using valuation rings and valued fields as examples, we discuss in which ways the notions of
More informationITERATED FUNCTION SYSTEMS WITH CONTINUOUS PLACE DEPENDENT PROBABILITIES
UNIVERSITATIS IAGELLONICAE ACTA MATHEMATICA, FASCICULUS XL 2002 ITERATED FUNCTION SYSTEMS WITH CONTINUOUS PLACE DEPENDENT PROBABILITIES by Joanna Jaroszewska Abstract. We study the asymptotic behaviour
More informationA sufficient condition for the existence of the Fourier transform of f : R C is. f(t) dt <. f(t) = 0 otherwise. dt =
Fourier transform Definition.. Let f : R C. F [ft)] = ˆf : R C defined by The Fourier transform of f is the function F [ft)]ω) = ˆfω) := ft)e iωt dt. The inverse Fourier transform of f is the function
More informationLanguage-Restricted Iterated Function Systems, Koch Constructions, and L-systems
Language-Restricted Iterated Function Systems, Koch Constructions, and L-systems Przemyslaw Prusinkiewicz and Mark Hammel Department of Computer Science University of Calgary Calgary, Alberta, Canada T2N
More informationDISCRETE FOURIER TRANSFORM
DD2423 Image Processing and Computer Vision DISCRETE FOURIER TRANSFORM Mårten Björkman Computer Vision and Active Perception School of Computer Science and Communication November 1, 2012 1 Terminology:
More informationsine wave fit algorithm
TECHNICAL REPORT IR-S3-SB-9 1 Properties of the IEEE-STD-57 four parameter sine wave fit algorithm Peter Händel, Senior Member, IEEE Abstract The IEEE Standard 57 (IEEE-STD-57) provides algorithms for
More informationCompensation. June 29, Arizona State University. Edge Detection from Nonuniform Data via. Convolutional Gridding with Density.
Arizona State University June 29, 2012 Overview gridding with density compensation is a common method for function reconstruction nonuniform Fourier data. It is computationally inexpensive compared with
More informationMath 6120 Fall 2012 Assignment #1
Math 6120 Fall 2012 Assignment #1 Choose 10 of the problems below to submit by Weds., Sep. 5. Exercise 1. [Mun, 21, #10]. Show that the following are closed subsets of R 2 : (a) A = { (x, y) xy = 1 },
More informationFourier-like Transforms
L 2 (R) Solutions of Dilation Equations and Fourier-like Transforms David Malone December 6, 2000 Abstract We state a novel construction of the Fourier transform on L 2 (R) based on translation and dilation
More informationTime-Dependent Statistical Mechanics A1. The Fourier transform
Time-Dependent Statistical Mechanics A1. The Fourier transform c Hans C. Andersen November 5, 2009 1 Definition of the Fourier transform and its inverse. Suppose F (t) is some function of time. Then its
More informationExamination paper for TMA4122/TMA4123/TMA4125/TMA4130 Matematikk 4M/N
Department of Mathematical Sciences Examination paper for TMA4122/TMA4123/TMA4125/TMA4130 Matematikk 4M/N Academic contact during examination: Markus Grasmair Phone: 97 58 04 35 Code C): Basic calculator.
More informationCh.11 The Discrete-Time Fourier Transform (DTFT)
EE2S11 Signals and Systems, part 2 Ch.11 The Discrete-Time Fourier Transform (DTFT Contents definition of the DTFT relation to the -transform, region of convergence, stability frequency plots convolution
More informationx[n] = x a (nt ) x a (t)e jωt dt while the discrete time signal x[n] has the discrete-time Fourier transform x[n]e jωn
Sampling Let x a (t) be a continuous time signal. The signal is sampled by taking the signal value at intervals of time T to get The signal x(t) has a Fourier transform x[n] = x a (nt ) X a (Ω) = x a (t)e
More informationUSING THE RANDOM ITERATION ALGORITHM TO CREATE FRACTALS
USING THE RANDOM ITERATION ALGORITHM TO CREATE FRACTALS UNIVERSITY OF MARYLAND DIRECTED READING PROGRAM FALL 205 BY ADAM ANDERSON THE SIERPINSKI GASKET 2 Stage 0: A 0 = 2 22 A 0 = Stage : A = 2 = 4 A
More informationWavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing
Wavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing Qingtang Jiang Abstract This paper is about the construction of univariate wavelet bi-frames with each framelet being symmetric.
More informationMTH 202 : Probability and Statistics
MTH 202 : Probability and Statistics Lecture 5-8 : 15, 20, 21, 23 January, 2013 Random Variables and their Probability Distributions 3.1 : Random Variables Often while we need to deal with probability
More informationFINAL PROJECT TOPICS MATH 399, SPRING αx 0 x < α(1 x)
FINAL PROJECT TOPICS MATH 399, SPRING 2011 MARIUS IONESCU If you pick any of the following topics feel free to discuss with me if you need any further background that we did not discuss in class. 1. Iterations
More informationThe Gaussian transform. ALECU, Teodor, VOLOSHYNOVSKYY, Svyatoslav, PUN, Thierry. Abstract
Proceedings Chapter The aussian transform ALECU, Teodor, VOLOSHYNOVSKYY, Svyatoslav, PUN, Thierry Abstract This paper introduces the general purpose aussian Transform, which aims at representing a generic
More informationZhenjiang, Jiangsu, , P.R. China (Received 7 June 2010, accepted xx, will be set by the editor)
ISSN 1749-3889 print), 1749-3897 online) International Journal of Nonlinear Science Vol.132012) No.3,pp.380-384 Fractal Interpolation Functions on the Stability of Vertical Scale Factor Jiao Xu 1, Zhigang
More informationLEAST Mean Squares Algorithm (LMS), introduced by
1 Hard Threshold Least Mean Squares Algorithm Lampros Flokas and Petros Maragos arxiv:1608.0118v1 [cs.sy] 3 Aug 016 Abstract This work presents a new variation of the commonly used Least Mean Squares Algorithm
More informationFourier Series. Fourier Transform
Math Methods I Lia Vas Fourier Series. Fourier ransform Fourier Series. Recall that a function differentiable any number of times at x = a can be represented as a power series n= a n (x a) n where the
More information1.1. MEASURES AND INTEGRALS
CHAPTER 1: MEASURE THEORY In this chapter we define the notion of measure µ on a space, construct integrals on this space, and establish their basic properties under limits. The measure µ(e) will be defined
More informationEconomics 204 Summer/Fall 2011 Lecture 5 Friday July 29, 2011
Economics 204 Summer/Fall 2011 Lecture 5 Friday July 29, 2011 Section 2.6 (cont.) Properties of Real Functions Here we first study properties of functions from R to R, making use of the additional structure
More informationNotes on Fourier Series and Integrals Fourier Series
Notes on Fourier Series and Integrals Fourier Series et f(x) be a piecewise linear function on [, ] (This means that f(x) may possess a finite number of finite discontinuities on the interval). Then f(x)
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationSC434L_DVCC-Tutorial 6 Subband Video Coding
SC434L_DVCC-Tutorial 6 Subband Video Coding Dr H.R. Wu Associate Professor Audiovisual Information Processing and Digital Communications Monash University http://www.csse.monash.edu.au/~hrw Email: hrw@csse.monash.edu.au
More informationHow many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation?
How many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation? (A) 0 (B) 1 (C) 2 (D) more than 2 (E) it depends or don t know How many of
More informationCHAOS/FRACTAL
CHAOS/FRACTAL 8 6.8.6 4.4. 4 5..4.6.8.4 8...8 6.8.6 - - 4.6.4.4 -.. - -6-5 -4 - - -..95.955.96.965.97.975.98..4.6.8 Presentation in Dynamical System. By Fred Khoury Introduction This paper is divided into
More informationMINIMAL DEGREE UNIVARIATE PIECEWISE POLYNOMIALS WITH PRESCRIBED SOBOLEV REGULARITY
MINIMAL DEGREE UNIVARIATE PIECEWISE POLYNOMIALS WITH PRESCRIBED SOBOLEV REGULARITY Amal Al-Rashdan & Michael J. Johnson* Department of Mathematics Kuwait University P.O. Box: 5969 Safat 136 Kuwait yohnson1963@hotmail.com*
More informationSignals, Instruments, and Systems W5. Introduction to Signal Processing Sampling, Reconstruction, and Filters
Signals, Instruments, and Systems W5 Introduction to Signal Processing Sampling, Reconstruction, and Filters Acknowledgments Recapitulation of Key Concepts from the Last Lecture Dirac delta function (
More informationThe Central Limit Theorem
The Central Limit Theorem (A rounding-corners overiew of the proof for a.s. convergence assuming i.i.d.r.v. with 2 moments in L 1, provided swaps of lim-ops are legitimate) If {X k } n k=1 are i.i.d.,
More informationP-adic Functions - Part 1
P-adic Functions - Part 1 Nicolae Ciocan 22.11.2011 1 Locally constant functions Motivation: Another big difference between p-adic analysis and real analysis is the existence of nontrivial locally constant
More informationREAL AND COMPLEX ANALYSIS
REAL AND COMPLE ANALYSIS Third Edition Walter Rudin Professor of Mathematics University of Wisconsin, Madison Version 1.1 No rights reserved. Any part of this work can be reproduced or transmitted in any
More informationThe Hausdorff Measure of the Attractor of an Iterated Function System with Parameter
ISSN 1479-3889 (print), 1479-3897 (online) International Journal of Nonlinear Science Vol. 3 (2007) No. 2, pp. 150-154 The Hausdorff Measure of the Attractor of an Iterated Function System with Parameter
More information6. Duals of L p spaces
6 Duals of L p spaces This section deals with the problem if identifying the duals of L p spaces, p [1, ) There are essentially two cases of this problem: (i) p = 1; (ii) 1 < p < The major difference between
More informationHomework 1 Real Analysis
Homework 1 Real Analysis Joshua Ruiter March 23, 2018 Note on notation: When I use the symbol, it does not imply that the subset is proper. In writing A X, I mean only that a A = a X, leaving open the
More informationV(t) = Total Power = Calculating the Power Spectral Density (PSD) in IDL. Thomas Ferree, Ph.D. August 23, 1999
Calculating the Power Spectral Density (PSD) in IDL Thomas Ferree, Ph.D. August 23, 1999 This note outlines the calculation of power spectra via the fast Fourier transform (FFT) algorithm. There are several
More informationFractals. Krzysztof Gdawiec.
Fractals Krzysztof Gdawiec kgdawiec@ux2.math.us.edu.pl 1 Introduction In the 1970 s Benoit Mandelbrot introduced to the world new field of mathematics. He named this field fractal geometry (fractus from
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science : Discrete-Time Signal Processing
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.34: Discrete-Time Signal Processing OpenCourseWare 006 ecture 8 Periodogram Reading: Sections 0.6 and 0.7
More informationBEE604 Digital Signal Processing
BEE64 Digital Signal Processing Copiled by, Mrs.S.Sherine Assistant Professor Departent of EEE BIHER. COTETS Sapling Discrete Tie Fourier Transfor Properties of DTFT Discrete Fourier Transfor Inverse Discrete
More informationIntroduction to Proofs in Analysis. updated December 5, By Edoh Y. Amiran Following the outline of notes by Donald Chalice INTRODUCTION
Introduction to Proofs in Analysis updated December 5, 2016 By Edoh Y. Amiran Following the outline of notes by Donald Chalice INTRODUCTION Purpose. These notes intend to introduce four main notions from
More informationSOLUTIONS TO THE FINAL EXAM
SOLUTIONS TO THE FINAL EXAM Short questions 1 point each) Give a brief definition for each of the following six concepts: 1) normal for topological spaces) 2) path connected 3) homeomorphism 4) covering
More informationOptimal Sojourn Time Control within an Interval 1
Optimal Sojourn Time Control within an Interval Jianghai Hu and Shankar Sastry Department of Electrical Engineering and Computer Sciences University of California at Berkeley Berkeley, CA 97-77 {jianghai,sastry}@eecs.berkeley.edu
More informationThe Order in Chaos. Jeremy Batterson. May 13, 2013
2 The Order in Chaos Jeremy Batterson May 13, 2013 Fractal geometry will make you see everything differently. There is a danger in reading further. You risk the loss of your childhood vision of clouds,
More informationMATH FALL 2014 HOMEWORK 10 SOLUTIONS
Problem 1. MATH 241-2 FA 214 HOMEWORK 1 SOUTIONS Note that u E (x) 1 ( x) is an equilibrium distribution for the homogeneous pde that satisfies the given boundary conditions. We therefore want to find
More informationOn Convergence of Sequences of Measurable Functions
On Convergence of Sequences of Measurable Functions Christos Papachristodoulos, Nikolaos Papanastassiou Abstract In order to study the three basic kinds of convergence (in measure, almost every where,
More informationNotions such as convergent sequence and Cauchy sequence make sense for any metric space. Convergent Sequences are Cauchy
Banach Spaces These notes provide an introduction to Banach spaces, which are complete normed vector spaces. For the purposes of these notes, all vector spaces are assumed to be over the real numbers.
More informationi. v = 0 if and only if v 0. iii. v + w v + w. (This is the Triangle Inequality.)
Definition 5.5.1. A (real) normed vector space is a real vector space V, equipped with a function called a norm, denoted by, provided that for all v and w in V and for all α R the real number v 0, and
More informationELLIPTIC CURVES BJORN POONEN
ELLIPTIC CURVES BJORN POONEN 1. Introduction The theme of this lecture is to show how geometry can be used to understand the rational number solutions to a polynomial equation. We will illustrate this
More informationSINGULAR INTEGRALS ON SIERPINSKI GASKETS
SINGULAR INTEGRALS ON SIERPINSKI GASKETS VASILIS CHOUSIONIS Abstract. We construct a class of singular integral operators associated with homogeneous Calderón-Zygmund standard kernels on d-dimensional,
More informationRecovery of Compactly Supported Functions from Spectrogram Measurements via Lifting
Recovery of Compactly Supported Functions from Spectrogram Measurements via Lifting Mark Iwen markiwen@math.msu.edu 2017 Friday, July 7 th, 2017 Joint work with... Sami Merhi (Michigan State University)
More informationREAL VARIABLES: PROBLEM SET 1. = x limsup E k
REAL VARIABLES: PROBLEM SET 1 BEN ELDER 1. Problem 1.1a First let s prove that limsup E k consists of those points which belong to infinitely many E k. From equation 1.1: limsup E k = E k For limsup E
More informationFractal Geometry Time Escape Algorithms and Fractal Dimension
NAVY Research Group Department of Computer Science Faculty of Electrical Engineering and Computer Science VŠB- TUO 17. listopadu 15 708 33 Ostrava- Poruba Czech Republic Basics of Modern Computer Science
More informationMathematical Methods for Physics and Engineering
Mathematical Methods for Physics and Engineering Lecture notes for PDEs Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 The integration theory
More informationSome Background Material
Chapter 1 Some Background Material In the first chapter, we present a quick review of elementary - but important - material as a way of dipping our toes in the water. This chapter also introduces important
More informationMATH 118, LECTURES 27 & 28: TAYLOR SERIES
MATH 8, LECTURES 7 & 8: TAYLOR SERIES Taylor Series Suppose we know that the power series a n (x c) n converges on some interval c R < x < c + R to the function f(x). That is to say, we have f(x) = a 0
More informationDiscontinuous order preserving circle maps versus circle homeomorphisms
Discontinuous order preserving circle maps versus circle homeomorphisms V. S. Kozyakin Institute for Information Transmission Problems Russian Academy of Sciences Bolshoj Karetny lane, 19, 101447 Moscow,
More informationLecture 9 February 2, 2016
MATH 262/CME 372: Applied Fourier Analysis and Winter 26 Elements of Modern Signal Processing Lecture 9 February 2, 26 Prof. Emmanuel Candes Scribe: Carlos A. Sing-Long, Edited by E. Bates Outline Agenda:
More informationSinusoidal Modeling. Yannis Stylianou SPCC University of Crete, Computer Science Dept., Greece,
Sinusoidal Modeling Yannis Stylianou University of Crete, Computer Science Dept., Greece, yannis@csd.uoc.gr SPCC 2016 1 Speech Production 2 Modulators 3 Sinusoidal Modeling Sinusoidal Models Voiced Speech
More informationOne-to-one functions and onto functions
MA 3362 Lecture 7 - One-to-one and Onto Wednesday, October 22, 2008. Objectives: Formalize definitions of one-to-one and onto One-to-one functions and onto functions At the level of set theory, there are
More informationExamination paper for TMA4130 Matematikk 4N
Department of Mathematical Sciences Examination paper for TMA4130 Matematikk 4N Academic contact during examination: Morten Nome Phone: 90 84 97 83 Examination date: 13 December 2017 Examination time (from
More informationTHE HUTCHINSON BARNSLEY THEORY FOR INFINITE ITERATED FUNCTION SYSTEMS
Bull. Austral. Math. Soc. Vol. 72 2005) [441 454] 39b12, 47a35, 47h09, 54h25 THE HUTCHINSON BARNSLEY THEORY FOR INFINITE ITERATED FUNCTION SYSTEMS Gertruda Gwóźdź- Lukawska and Jacek Jachymski We show
More information