W. Kozek and H.G. Feichtinger. NUHAG, Dept. of Mathematics, University of Vienna. J. Scharinger. Altenberger Str. 69, A-4040 Linz, Austria ABSTRACT
|
|
- Josephine Hensley
- 5 years ago
- Views:
Transcription
1 GABOR ANALYSIS AND LINEAR SYSTEM IDENTIFICATION 3 W. Kozek and H.G. Feichtinger NUHAG, Dept. of Mathematics, University of Vienna Strudlhofg.4, A{1090 Vienna, Austria ( kozek@tyche.mat.univie.ac.at) J. Scharinger Institute of Systems Science, Johannes Kepler University Linz Altenberger Str. 69, A-4040 Linz, Austria ( js@cast.uni-linz.ac.at) ABSTRACT Many practically important linear systems are slowly time{varying. In this paper we consider the frequency{domain identication of slowly time{varying systems with limited impulse response (underspread systems). We review time{frequency representations of LTV systems and its connection to the Gabor expansion of signals. It is shown how the matched multiwindow Gabor expansion can be used for the optimum frequency{domain identication and canonical realization of underspread LTV systems. 1. INTRODUCTION Nonparametric identication of linear time{varying (LTV) systems is a problem of particular relevance for communication engineering. For instance, LTV wave propagation channels are characterized by known partial dierential equations but highly variable and unstructured side constraints. Moreover, in control engineering a nonparametric approach may be appropriate in case of modelling uncertainty as a rst step towards model selection. Linear system theory is dominated by the treatment of the time{invariant case. For LTI systems the frequency domain characterization (transfer function) is matched in the sense that the "basis" functions of the Fourier transform are eigenfunctions of the system. The eigenfunctions of a general LTV system are a priori unknown, but in case of slow time{variation one may assume that windowed sinusoids are approximate eigenfunctions of the system. According to this assumption one characterizes the LTV system by a time{varying transfer function, i.e., a joint function of time and frequency. System identication can then 3 Funding by the FWF grant S7001-MAT and OENB grant 4913.
2 be performed by the observation of windowed versions of the input and output signal leading to point estimates of the transfer function. When the input signal is approximately white noise over some frequency band of interest, then one can obtain a global estimate of the transfer function for this band. Obviously, such a procedure cannot work for a general LTV systems with arbitrarily fast time{ variations. However, detailed theoretical results in the sense of a bias/variance analysis seem to exist only for the stationary case. In this paper we report about recent results on the frequency{domain identication of LTV systems. A brief review of nonparametric LTV system theory is followed by a discussion of the Gabor expansion and its application to the realization of underspread systems. We present explicit bias and variance results for the transfer function estimation and show that optimum estimators can be formulated in terms of Gabor expansions of the input/output signal. 2. LTV SYSTEM REPRESENTATION We consider LTV systems which act as integral operators on L 2 {signals. The input{output relation is given by 1 : (Hx) (t) := h(t; )x(t 0 )d = k H (t; s)x(s)ds (1) where k H (t; s) is the kernel (impulse response) and h(t; ) := k H (t; t 0 ) is the delay{spread function [10]. For engineering applications, the delay{spread function is often advantageous since it depends separately on absolute time t and time delay. I.e., for an LTI system h(t; ) is invariant w.r.t. t. Apart from limit case considerations we assume square{integrable kernels that correspond to Hilbert{Schmidt (HS) operators [7] with an inner product hg; Hi := k G (t; s)k H (t; s)dt ds (2) and associated HS norm khk 2 := hh; Hi (overline denotes complex conjugation). Spreading Function. The spreading function of an LTV system is dened as the Fourier transform of h(t; ) w.r.t. absolute time [10]: S H (; ) := h(t; )e 0j2t dt; (3) it is generally complex{valued even in case of a real{valued impulse response. Slowly time{varying systems are characterized by a limited {support of S H (; ). In a more symmetric way, we can interpret S H (; ) as time{frequency{ shifting representation of an LTV system in the sense that the output signal is 1 All sums and integrals go from 01 to 1 unless otherwise stated.
3 given by a superposition of time{shifted (delayed) and frequency{shifted (modulated) versions of the input signal: (Hx) (t) = S H (; )x(t 0 )e j2t d d; (4) For later use we introduce a time{shift operator (T x) (t) := x(t 0 ); and a frequency shift operator: (M x) (t) := x(t)e j2t : This allows to reformulate (3) in a signal free formulation H = S H (; )M T d d; (5) and we can formally dene the spreading function via an inner product S H (; ) = hh ; M T i : (6) The spreading function is a unique system representation since the map h(t; ) 7! S H (; ) is unitary. Time{Varying Transfer Function. The spreading function characterizes the energy preserving eects, time delay and frequency shift, in a pointwise clear manner. However, for the study of ltering eects of an LTV system it is more natural to represent the system by a time{varying transfer function. Such a transfer function can be dened by a Fourier transform w.r.t. the lag variable [11]: H (t; f) := h(t; )e 0j2f d; (7) which reduces to the usual time{independent transfer function in case of an LTI system. The time{varying transfer function is the symplectic Fourier transform of the spreading function: H (t; f) = S H (; )e 0j2(f 0t) d d: (8) In order to formulate the system identication problem in a compact manner we introduce time{frequency shifting of operators by P(t; f) := M f T t PT (0t) M (0f ) ; (9) i.e., P(t; f) is unitarily equivalent to P and its eigensignals are time{frequency{ shifted versions of that of P. With other words, when P selects signals temporally and spectrally localized about zero, then P(t; f) selects signals temporally and spectrally localized at (t; f). The transfer function of P(; ) gives a clear picture: P (;) (t; f) = P (t 0 ; f 0 ) (10) The kernel of P(; ) is given by: (P(; )) (t; s) = k P (t0; s0)e j2(t0s) ; when k P (t; s) is the kernel of P.
4 With this notational convention we can introduce a time{frequency{parametrized system representation dual to (5) (continuous Weyl{Heisenberg expansion of an LTV system) [4]: H = H (t; f)p 0 (t; f)dt df (11) where P 0 is a \prototype" operator dened by H (t; f) = hh; P 0 (t; f)i () P0 (t; f) = (t)(f): (12) The label time{varying transfer function suggests that H (t; f) reects time{ frequency{selective ltering in a manner just as the time{independent LTI system's transfer function characterizes frequency{selective ltering. However, this interpretation would require that the prototype operator P 0 (t; f) performs a perfect time{frequency localization in obvious contradiction to Heisenberg's uncertainty principle. Hence, one may already expect that for a general LTV system H (t; f) does not have much practical relevance. 3. UNDERSPREAD SYSTEMS Slowly time{varying systems with limited time delay can be characterized by a compactly supported spreading function ( [a;b] (:) denotes the indicator function of the interval [a; b] on the real line): S H (; ) = S H (; ) [00 ; 0 ]() [00 ; 0 ]() with H = ; (13) where 0 characterizes the maximum frequency shift (proportional to the maximum velocity of time{variations) and 0 characterizes the maximum time shift. Systems with H 1 will be called underspread [4]. An important practical example is the mobile radio propagation channel [8]. Transfer Function Calculus. LTI systems can be viewed as a limit case of underspread systems with 0! 0. Moreover, one can show that all of the well{ known properties of an LTI system's transfer function hold approximately for the time{varying transfer function of an underspread system [6]. Here, we mention only one important property of the transfer function for two jointly underspread (sub)systems G, H: the transfer function of the composite system is essentially given by the product of the subsystems' transfer functions 2 k GH 0 G H k 2 < c 3 H kgk2 khk 2 ; (14) where c is a system independent constant. Two LTI systems are jointly underspread with H = 0 such that (14) implies the well known identity: GH (f) = G (f) H (f): We emphasize that for overspread systems ( H > 1) the error bound in (14) gets dominant, and, more generally, the time{varying transfer 2 Jointly underspread means that both systems satisfy the same spreading constraint.
5 function loses all of the important properties of the LTI system's transfer function. Discrete Weyl{Heisenberg Expansion. According to (8) the transfer function of an underspread system is a 2D lowpass function. Hence, using the sampling theorem of 2D functions [1] and (10) one can formulate a discretized version of (11) based on a dierent prototype operator: H = 1 X H k X l H (kt; lf )P (kt; lf ) (15) with the sampling grid and the prototype operator matched to the underspread support as follows T F = 0 0 and T F = 1 H ; (16) S P (; ) = [00 ; 0 ]() [00 ; 0 ]() () H (t; f) = hh; P (t; f)i (17) The validity of a discrete expansion makes the transfer function calculus to a feasible framework for signal processing or control engineering problems where one needs to compute composite systems or inverse systems in a numerically ecient way. 4. GABOR EXPANSION The time{varying transfer function is a joint function of time and frequency. This suggests that time{frequency signal representations will play an important role in the realization of transfer function estimators that are based on input{output observations. In this section we give a brief review of linear time{frequency signal representations and their application to LTV system realization. Short{Time Fourier Transform. The short{time Fourier transform (STFT) of a signal x(t) is dened as [9, 3] (ST F T x) (t; f) := x(s)(s 0 t)e 0j2f s ds = hx; M f T t i ; where (t) is the analysis window. window and is highly redundant. The STFT depends signicantly on the Gabor Expansion. The Gabor expansion of a signal x(t) is dened as [3, 12] X X x(t) = (G x) [k; l] (M lf T kt g) (t) (18) k l where the coecients can be dened as a sampled STFT: (G x) [k; l] := (ST F T x) (kt; lf ); (19)
6 and g(t) is called synthesis window. Note that the validity of (18) imposes severe constraints on the choice of the analysis/synthesis windows [3, 12]. System Representation via Multiplicative Modication. Using our notation for time{frequency{shifted operators (see (9)) it is straightforward to write (18) in a signal{free formulation: I X X = P g; (kt; lf ); k l where I denotes the identity operator and (P g; ) (t; s) := g(t)(s). For a synthesis of LTV systems it is natural to consider multiplicative perturbations of this resolution of the identity [2]: H X X G = M G [k; l]p g; (kt; lf ); (20) k l where M G [k; l] is a multiplicator function that characterizes the time{frequency ltering behaviour of H G in a manner similar to a sampled transfer function. Eq.(20) is just a discrete version of multiplicative modication of the STFT [9] where the system is given by H ST F T = M(t; f)p g; (t; f)dt df: (21) Based upon appropriately selected analysis and synthesis windows one can show that any underspread system can be realized by such an STFT lter [6]. This does not hold true for the Gabor lters as dened by (20). One way to overcome the basic restrictions of Gabor lters is the introduction of multiwindow methods [5, 12]. Consider the singular value decomposition of the prototype operator in (15) 1X P = s m P um;vm (22) m=1 here s m are the singular values, u m (t); v m (t) are the singular signals, and (P um;vm ) (t; s) := u m (t)v m (s). Now, when we insert (22) in (15) comparison with (20) shows that we can realize any underspread system by multiplicative modication of a multiwindow Gabor expansion. The windows are the singular functions of the prototype operator thus determined by the spreading constraint (13), and the multiplicator function is given by the samples of H (t; f). A precise realization would require an innite number of windows. In practice, one has to truncate (22) in the sense of a nite{rank approximation leading to an approximation of (15) in the following way: NX X X fh = s m H (kt; lf )P um;vm (kt; lf ): (23) m=1 k l For a given N, this approximation is optimal (in a HS sense) when the indices of (22) correspond to nonincreasing order of the singular values.
7 5. TRANSFER FUNCTION ESTIMATION We now show that the above introduced prototype operator P also establishes an optimum estimator of the time{varying transfer function H (t; f) under convenient statistical assumptions. System Identication Setup. We assume an observation of (i) a circular complex, stationary white Gaussian input process and (ii) the corresponding output process contaminated by statistically independent, circular complex, stationary additive white Gaussian noise. We consider point estimation of the time{varying transfer function dened by its prototype operator P 0 in the form: H (t; f) = hh; P 0 (t; f)i For notational simplicity we assume an input process with normalized variance Efx(t)x(s)g = (t 0 s): According to the above assumptions, the output signal (Hx) (t) is a realization of a zero{mean, circular complex Gaussian nonstationary process. We assume a noiseless observation of the input signal x(t) and a noisy observation of the output signal: y(t) = (Hx) (t) + n(t); where n(t) is circular complex, zero{mean, stationary white Gaussian noise with covariance: Efn(t)n(s)g = n 2 (t 0 s): Estimator. The estimator is dened as a sesquilinear form based on a dierent prototype operator b P: b H (t; f) = D y; b P(t; f)xe Based on the statistical independence of x(t) and n(t) it is straightforward to show that the expectation of the estimate results in The bias is then given by: E n b H (t; f) o = D H; b P(t; f) E : B(t; f) := E n b H (t; f) o 0 H (t; f) = D H; e P(t; f) E ; where we have introduced a bias operator dened as: ep := b P 0 P 0 : There exist two essentially dierent strategies for unbiased estimation: (I) Unbiased estimators with P b = P 0. (II) Unbiased estimators with P b 6= P but nonoverlapping spreading functions of the bias operator and the system: S H (; )Se P (; ) 0: (24)
8 The following variance analysis will show that the latter is the practically relevant case. Variance. The variance of the estimator is dened as V 2 (t; f) := E b H (t; f) 2 0 E n bh (t; f) o 2 ; in general it depends on time and frequency. One can show [4] that the variance of the estimator is given by ( 3 denotes the adjoint operator): V 2 (t; f) = D HH 3 ; bp b P 3 (t; f) E + 2 nk b Pk 2 (25) Note that the constant term is proportional to the HS norm of the prototype operator. This shows that we can already exclude nite variance unbiased estimation of the time{varying transfer function of a general LTV system (without a priori knowledge) since the prototype operator P 0 (see (12)) is not HS. The variance is depends on time and frequency. However, in practice we would like to use one and the same point estimator for any point of the time{frequency plane. In order to optimize an unbiased time{frequency invariant estimator we have to consider global measures for the variance. First, we have a tight bound for the maximum variance: V 2 (t; f) k b Pk 2 khk2 + 2 n Second, we compute an integrated variance by integrating over the time{ frequency dependent term of (25): V 2 0 := t f V 2 (t; f) 0 2 nk b Pk 2 dt df = k b Pk 2 khk 2 Both the integrated and the maximum variance are proportional to the HS norm of b P. Hence, an optimum time{frequency independent unbiased estimator has to minimize the HS norm of its prototype operator without introducing bias. However, for underspread LTV systems with known spreading constraint this optimization is seen to be trivial when formulated in terms of Sb P (; ): The spreading function must be equal to one within the support of S P (; ) and arbitrary elsewhere (cf. (24)). Since the HS norm of b P is equal to the L 2 {norm of Sb P (; ), we have to minimize its support which just leads back to the prototype operator introduced in (15) bp opt = P : Note that our bias/variance analysis holds for arbitrary support of the system's spreading function. Estimator Realization via Gabor Analysis. For a practical (digital) implementation we have to reformulate the estimator in terms of the STFT of the :
9 input and observed output. To this end, note that our specic sesquilinear form can be written as a product of the STFTs of the output and input signal: hy; P g; (t; f)xi = hy; M f T t i hx; M f T t gi = (ST F T y) (t; f)(st F T g x) (t; f): Then, as discussed in the previous section we introduce a truncated singular value decomposition of the prototype operator such that a canonical estimator of the critically sampled transfer function is given by: b H (kt; lf ) = NX m=1 s m (G vm y) [k; l](g um x) [k; l]; i.e., up to a conjugation we have a product of the multiwindow Gabor coecients of the input signal and the output signal. In a digital implementation, this formulation is structurally simple and permits the use of FFTs. Observe that the estimator realization is structurally equivalent to the system realization via multiplicative modication of the Gabor expansion (see (23)). Consequently, one can use the Gabor coecients not only for the system identi- cation but also for a subsequent processing of the input or output signal. This makes our approach a candidate for adaptive signal processing matched to nonstationary environments. 6. CONCLUSIONS We have studied the frequency domain identication of linear time{varying (LTV) systems. The discussion has been restricted to underspread systems whose frequency{domain characterization achieves the essential properties of the frequency{domain characterization of LTI systems. We have pointed out the realization of such underspread systems via multiplicative modication of a multiwindow Gabor expansion. A canonical class of point estimators in terms of the input and the output signal has been introduced. Under the assumption of a white{noise input signal and observation of the output signal subject to additive white noise, one can derive closed form expressions for the bias and variance of the point estimate. We have then specialized to unbiased time{frequency{ invariant estimators and evaluated their performance in the sense of a maximum and integrated variance. The theoretically optimum estimator has been determined by its spreading function: It is an indicator function with the support given by the spreading constraint of the system at hand. These theoretical estimators can be approximately realized as the product of a matched multiwindow Gabor expansion of the known input signal and the observed output signal.
10 7. REFERENCES [1] D.E. Dudgeon and R.M. Merserau. Multidimensional Signal Processing. Prentice Hall, Englewood Clis (NJ), [2] S. Farkash and S. Raz. Time Variant Filtering via the Gabor Representation. Signal Processing V, pages 509{512, [3] H.G. Feichtinger and K. Grochenig. Gabor wavelets and the Heisenberg group: Gabor Expansions and Short Time Fourier transform from the group theoretical point of view. In C.K. Chui, editor, Wavelets { A Tutorial in Theory and Applications, pages 359{397. Academic Press, Boston, [4] W. Kozek. On underspread WSSUS channels: Calculus of transfer functions identication and signal design. Submitted to IEEE Trans. Signal Processing, Special Issue on Advanced Communications. [5] W. Kozek. On the generalized Weyl correspondence and its application to time{frequency analysis of linear time{varying systems. In Proc. IEEE Int.Symp.Time{Frequency/Time{Scale Analysis, pages 167{170, Victoria (Canada), [6] W. Kozek. Matched Weyl{Heisenberg expansions of nonstationary environments. PhD thesis, Vienna University of Technology, in preparation. [7] A.W. Naylor and G.R. Sell. Linear Operator Theory in Engineering and Science. Springer, New York, [8] J.D. Parsons. The Mobile Radio Propagation Channel. Pentech Press, London, [9] M.R. Portno. Time-Frequency Representation of Digital Signals and Systems Based on Short-Time Fourier Analysis. IEEE Trans. Sign. Proc., 28:55{69, [10] K.A. Sostrand. Mathematics of the time{varying channel. Proceedings NATO Advanced Study Institute on Signal Procesing with Emphasis on Underwater Acoustics, 2:25{1{25{20, [11] L. adeh. Frequency analysis of variable networks. Proc. of IRE, 67:291{299, March [12] M. ibulski and Y.Y. eevi. Multiwindow Gabor-type Schemes by Frame Methods. Technical Report CC PUB 101, Technion{Isreal Institute of Technology, Center for Communication, Haifa, Isreal, 1995.
IEEE Trans. Information Theory, vol. 52, no. 3, Mar. 2006, pp , Copyright IEEE 2006
IEEE Trans. Information Theory, vol. 52, no. 3, Mar. 2006, pp. 1067 1086, Copyright IEEE 2006 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 52, NO. 3, MARCH 2006 1067 Nonstationary Spectral Analysis Based
More informationTIME-FREQUENCY ANALYSIS: TUTORIAL. Werner Kozek & Götz Pfander
TIME-FREQUENCY ANALYSIS: TUTORIAL Werner Kozek & Götz Pfander Overview TF-Analysis: Spectral Visualization of nonstationary signals (speech, audio,...) Spectrogram (time-varying spectrum estimation) TF-methods
More information1 Introduction Spectral analysis of (wide-sense) stationary random processes by means of the power spectral density (PSD) is a useful concept. owever,
IEEE Trans. Signal Proc., vol. 45, no. 6, pp. 150{1534, June 1997. Generalized Evolutionary Spectral Analysis and the Weyl Spectrum of Nonstationary andom Processes SP-865, nal version Gerald Matz y, Franz
More informationElec4621 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis
Elec461 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis Dr. D. S. Taubman May 3, 011 In this last chapter of your notes, we are interested in the problem of nding the instantaneous
More informationCovariant Time-Frequency Representations. Through Unitary Equivalence. Richard G. Baraniuk. Member, IEEE. Rice University
Covariant Time-Frequency Representations Through Unitary Equivalence Richard G. Baraniuk Member, IEEE Department of Electrical and Computer Engineering Rice University P.O. Box 1892, Houston, TX 77251{1892,
More informationDigital Filter Design Using Non-Uniform Sampling
Proceedings of the 6th WSEAS International Conference on Applications of Electrical Engineering, Istanbul, Turkey, May 27-29, 2007 120 Digital Filter Design Using Non-Uniform Sampling K. BUSAWON, WAI-PANG
More informationSYSTEM RECONSTRUCTION FROM SELECTED HOS REGIONS. Haralambos Pozidis and Athina P. Petropulu. Drexel University, Philadelphia, PA 19104
SYSTEM RECOSTRUCTIO FROM SELECTED HOS REGIOS Haralambos Pozidis and Athina P. Petropulu Electrical and Computer Engineering Department Drexel University, Philadelphia, PA 94 Tel. (25) 895-2358 Fax. (25)
More informationExamples. 2-input, 1-output discrete-time systems: 1-input, 1-output discrete-time systems:
Discrete-Time s - I Time-Domain Representation CHAPTER 4 These lecture slides are based on "Digital Signal Processing: A Computer-Based Approach, 4th ed." textbook by S.K. Mitra and its instructor materials.
More informationA Fast Algorithm for. Nonstationary Delay Estimation. H. C. So. Department of Electronic Engineering, City University of Hong Kong
A Fast Algorithm for Nonstationary Delay Estimation H. C. So Department of Electronic Engineering, City University of Hong Kong Tat Chee Avenue, Kowloon, Hong Kong Email : hcso@ee.cityu.edu.hk June 19,
More informationApproximately dual frames in Hilbert spaces and applications to Gabor frames
Approximately dual frames in Hilbert spaces and applications to Gabor frames Ole Christensen and Richard S. Laugesen October 22, 200 Abstract Approximately dual frames are studied in the Hilbert space
More informationOPTIMUM ARRAY PROCESSING AND REPRESENTATION OF NONSTATIONARY RANDOM SCATTERING (Invited Paper)
OPTIMUM ARRAY PROCESSING AND REPRESENTATION OF NONSTATIONARY RANDOM SCATTERING (Invited Paper) LEON H. SIBUL, MAEX L, FOWLER, AND GAVIN J. HARBISON Applied Research Laboratory The Pennsylvania State University
More informationA New Subspace Identification Method for Open and Closed Loop Data
A New Subspace Identification Method for Open and Closed Loop Data Magnus Jansson July 2005 IR S3 SB 0524 IFAC World Congress 2005 ROYAL INSTITUTE OF TECHNOLOGY Department of Signals, Sensors & Systems
More informationGaussian Basics Random Processes Filtering of Random Processes Signal Space Concepts
White Gaussian Noise I Definition: A (real-valued) random process X t is called white Gaussian Noise if I X t is Gaussian for each time instance t I Mean: m X (t) =0 for all t I Autocorrelation function:
More informationFilter Banks for Image Coding. Ilangko Balasingham and Tor A. Ramstad
Nonuniform Nonunitary Perfect Reconstruction Filter Banks for Image Coding Ilangko Balasingham Tor A Ramstad Department of Telecommunications, Norwegian Institute of Technology, 7 Trondheim, Norway Email:
More informationComparison of DDE and ETDGE for. Time-Varying Delay Estimation. H. C. So. Department of Electronic Engineering, City University of Hong Kong
Comparison of DDE and ETDGE for Time-Varying Delay Estimation H. C. So Department of Electronic Engineering, City University of Hong Kong Tat Chee Avenue, Kowloon, Hong Kong Email : hcso@ee.cityu.edu.hk
More informationA Multi-window Fractional Evolutionary Spectral Analysis
A Multi-window Fractional Evolutionary Spectral Analysis YALÇIN ÇEKİÇ, AYDIN AKAN, and MAHMUT ÖZTÜRK University of Bahcesehir, Department of Electrical and Electronics Engineering Bahcesehir, 49, Istanbul,
More informationFourier Analysis Linear transformations and lters. 3. Fourier Analysis. Alex Sheremet. April 11, 2007
Stochastic processes review 3. Data Analysis Techniques in Oceanography OCP668 April, 27 Stochastic processes review Denition Fixed ζ = ζ : Function X (t) = X (t, ζ). Fixed t = t: Random Variable X (ζ)
More informationDigital Image Processing Lectures 15 & 16
Lectures 15 & 16, Professor Department of Electrical and Computer Engineering Colorado State University CWT and Multi-Resolution Signal Analysis Wavelet transform offers multi-resolution by allowing for
More informationRelative Irradiance. Wavelength (nm)
Characterization of Scanner Sensitivity Gaurav Sharma H. J. Trussell Electrical & Computer Engineering Dept. North Carolina State University, Raleigh, NC 7695-79 Abstract Color scanners are becoming quite
More information2 DEMETRIO LABATE the mask of the lter, since it eights selectively the dierent time{frequency components of the signal. In our approach, the pseudodi
Preprint (revised) 07/18/00 TIME{FREQUENCY ANALYSIS OF PSEUDODIFFERENTIAL OPERATORS Demetrio Labate Abstract. In this paper e apply a time{frequency approach to the study of pseudodierential operators.
More informationIntroduction to Gabor Analysis
Theoretical and Computational Aspects Numerical Harmonic Group under the supervision of Prof. Dr. Hans Georg Feichtinger 30 Oct 2012 Outline 1 2 3 4 5 6 7 DFT/ idft Discrete Given an input signal f of
More information7 The Waveform Channel
7 The Waveform Channel The waveform transmitted by the digital demodulator will be corrupted by the channel before it reaches the digital demodulator in the receiver. One important part of the channel
More informationApproximately dual frame pairs in Hilbert spaces and applications to Gabor frames
arxiv:0811.3588v1 [math.ca] 21 Nov 2008 Approximately dual frame pairs in Hilbert spaces and applications to Gabor frames Ole Christensen and Richard S. Laugesen November 21, 2008 Abstract We discuss the
More informationOle Christensen 3. October 20, Abstract. We point out some connections between the existing theories for
Frames and pseudo-inverses. Ole Christensen 3 October 20, 1994 Abstract We point out some connections between the existing theories for frames and pseudo-inverses. In particular, using the pseudo-inverse
More information3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE
3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE 3.0 INTRODUCTION The purpose of this chapter is to introduce estimators shortly. More elaborated courses on System Identification, which are given
More informationNoise Reduction in Oversampled Filter Banks Using Predictive Quantization
IEEE TRANSACTIONS ON INFORMATION THEORY, VOL 47, NO 1, JANUARY 2001 155 Noise Reduction in Oversampled Filter Banks Using Predictive Quantization Helmut Bölcskei, Member, IEEE, Franz Hlawatsch, Member,
More informationDigital Image Processing
Digital Image Processing, 2nd ed. Digital Image Processing Chapter 7 Wavelets and Multiresolution Processing Dr. Kai Shuang Department of Electronic Engineering China University of Petroleum shuangkai@cup.edu.cn
More informationREGLERTEKNIK AUTOMATIC CONTROL LINKÖPING
Time-domain Identication of Dynamic Errors-in-variables Systems Using Periodic Excitation Signals Urban Forssell, Fredrik Gustafsson, Tomas McKelvey Department of Electrical Engineering Linkping University,
More informationTopic 3: Fourier Series (FS)
ELEC264: Signals And Systems Topic 3: Fourier Series (FS) o o o o Introduction to frequency analysis of signals CT FS Fourier series of CT periodic signals Signal Symmetry and CT Fourier Series Properties
More informationGaussian process for nonstationary time series prediction
Computational Statistics & Data Analysis 47 (2004) 705 712 www.elsevier.com/locate/csda Gaussian process for nonstationary time series prediction Soane Brahim-Belhouari, Amine Bermak EEE Department, Hong
More informationFilter Banks with Variable System Delay. Georgia Institute of Technology. Atlanta, GA Abstract
A General Formulation for Modulated Perfect Reconstruction Filter Banks with Variable System Delay Gerald Schuller and Mark J T Smith Digital Signal Processing Laboratory School of Electrical Engineering
More informationEQUIVALENCE OF DFT FILTER BANKS AND GABOR EXPANSIONS. Helmut Bolcskei and Franz Hlawatsch
SPIE Proc. Vol. 2569 \Wavelet Applications in Signal and Image Processing III" San Diego, CA, July 995. EQUIVALENCE OF DFT FILTER BANKS AND GABOR EPANSIONS Helmut Bolcskei and Franz Hlawatsch INTHFT, Technische
More informationApplications and fundamental results on random Vandermon
Applications and fundamental results on random Vandermonde matrices May 2008 Some important concepts from classical probability Random variables are functions (i.e. they commute w.r.t. multiplication)
More informationButterworth Filter Properties
OpenStax-CNX module: m693 Butterworth Filter Properties C. Sidney Burrus This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 3. This section develops the properties
More informationChapter 7 Interconnected Systems and Feedback: Well-Posedness, Stability, and Performance 7. Introduction Feedback control is a powerful approach to o
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 7 Interconnected
More informationOctober 7, :8 WSPC/WS-IJWMIP paper. Polynomial functions are renable
International Journal of Wavelets, Multiresolution and Information Processing c World Scientic Publishing Company Polynomial functions are renable Henning Thielemann Institut für Informatik Martin-Luther-Universität
More informationThe Bias-Variance dilemma of the Monte Carlo. method. Technion - Israel Institute of Technology, Technion City, Haifa 32000, Israel
The Bias-Variance dilemma of the Monte Carlo method Zlochin Mark 1 and Yoram Baram 1 Technion - Israel Institute of Technology, Technion City, Haifa 32000, Israel fzmark,baramg@cs.technion.ac.il Abstract.
More informationMultidimensional partitions of unity and Gaussian terrains
and Gaussian terrains Richard A. Bale, Jeff P. Grossman, Gary F. Margrave, and Michael P. Lamoureu ABSTRACT Partitions of unity play an important rôle as amplitude-preserving windows for nonstationary
More informationNOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION. M. Schwab, P. Noll, and T. Sikora. Technical University Berlin, Germany Communication System Group
NOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION M. Schwab, P. Noll, and T. Sikora Technical University Berlin, Germany Communication System Group Einsteinufer 17, 1557 Berlin (Germany) {schwab noll
More informationFRAMES AND TIME-FREQUENCY ANALYSIS
FRAMES AND TIME-FREQUENCY ANALYSIS LECTURE 5: MODULATION SPACES AND APPLICATIONS Christopher Heil Georgia Tech heil@math.gatech.edu http://www.math.gatech.edu/ heil READING For background on Banach spaces,
More informationEE538 Final Exam Fall :20 pm -5:20 pm PHYS 223 Dec. 17, Cover Sheet
EE538 Final Exam Fall 005 3:0 pm -5:0 pm PHYS 3 Dec. 17, 005 Cover Sheet Test Duration: 10 minutes. Open Book but Closed Notes. Calculators ARE allowed!! This test contains five problems. Each of the five
More informationCommunications and Signal Processing Spring 2017 MSE Exam
Communications and Signal Processing Spring 2017 MSE Exam Please obtain your Test ID from the following table. You must write your Test ID and name on each of the pages of this exam. A page with missing
More informationIntroduction to Wavelet. Based on A. Mukherjee s lecture notes
Introduction to Wavelet Based on A. Mukherjee s lecture notes Contents History of Wavelet Problems of Fourier Transform Uncertainty Principle The Short-time Fourier Transform Continuous Wavelet Transform
More informationDISCRETE-TIME SIGNAL PROCESSING
THIRD EDITION DISCRETE-TIME SIGNAL PROCESSING ALAN V. OPPENHEIM MASSACHUSETTS INSTITUTE OF TECHNOLOGY RONALD W. SCHÄFER HEWLETT-PACKARD LABORATORIES Upper Saddle River Boston Columbus San Francisco New
More informationA Radon Transform Based Method to Determine the Numerical Spread as a Measure of Non-Stationarity
A Radon Transform Based Method to Determine the Numerical Spread as a Measure of Non-Stationarity ROBERT A. HEDGES and BRUCE W. SUTER Air Force Research Laboratory IFGC 525 Brooks Rd. Rome, NY 13441-455
More informationPART 1. Review of DSP. f (t)e iωt dt. F(ω) = f (t) = 1 2π. F(ω)e iωt dω. f (t) F (ω) The Fourier Transform. Fourier Transform.
PART 1 Review of DSP Mauricio Sacchi University of Alberta, Edmonton, AB, Canada The Fourier Transform F() = f (t) = 1 2π f (t)e it dt F()e it d Fourier Transform Inverse Transform f (t) F () Part 1 Review
More informationThe DFT as Convolution or Filtering
Connexions module: m16328 1 The DFT as Convolution or Filtering C. Sidney Burrus This work is produced by The Connexions Project and licensed under the Creative Commons Attribution License A major application
More informationA WAVELET FILTER BANK WHICH MINIMIZES A NOVEL TRANSLATION INVARIANT DISCRETE UNCERTAINTY MEASURE. P. Tay, J.P. Havlicek, and V.
A WAVELET FILTER BAK WHICH MIIMIZES A OVEL TRASLATIO IVARIAT DISCRETE UCERTAITY MEASURE P. Tay, J.P. Havlicek, and V. DeBrunner School of Electrical and Computer Engineering University of Oklahoma, orman,
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationCompressive Identification of Linear Operators
Compressive Identification of Linear Operators Reinhard Heckel and Helmut Bölcskei Department of Information Technology and Electrical Engineering, ETH Zurich E-mail: heckel,boelcskei}@nari.ee.ethz.ch
More informationBasics about Fourier analysis
Jérôme Gilles UCLA PART ONE Fourier analysis On the menu... Introduction - some history... Notations. Fourier series. Continuous Fourier transform. Discrete Fourier transform. Properties. 2D extension.
More informationLecture 1 Some Time-Frequency Transformations
Lecture 1 Some Time-Frequency Transformations David Walnut Department of Mathematical Sciences George Mason University Fairfax, VA USA Chapman Lectures, Chapman University, Orange, CA 6-10 November 2017
More informationMultiscale Image Transforms
Multiscale Image Transforms Goal: Develop filter-based representations to decompose images into component parts, to extract features/structures of interest, and to attenuate noise. Motivation: extract
More informationDetection and Quantification of the Influence of Time-Variation in Frequency Response Function Measurements Using Arbitrary Excitations
Detection and Quantification of the Influence of Time-Variation in Frequency Response Function Measurements Using Arbitrary Excitations R. Pintelon, E. Louarroudi, and J. Lataire, Vrije Universiteit Brussel,
More informationTHE ZAK TRANSFORM (ZT), which is also known
Copyright IEEE 1997 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 45, NO. 4, APRIL 1997 851 Discrete Zak Transms, Polyphase Transms, and Applications Helmut Bölcskei, Student Member, IEEE, and Franz Hlawatsch,
More informationCOMPLEX WAVELET TRANSFORM IN SIGNAL AND IMAGE ANALYSIS
COMPLEX WAVELET TRANSFORM IN SIGNAL AND IMAGE ANALYSIS MUSOKO VICTOR, PROCHÁZKA ALEŠ Institute of Chemical Technology, Department of Computing and Control Engineering Technická 905, 66 8 Prague 6, Cech
More informationOn the Hilbert Transform of Wavelets
On the Hilbert Transform of Wavelets Kunal Narayan Chaudhury and Michael Unser Abstract A wavelet is a localized function having a prescribed number of vanishing moments. In this correspondence, we provide
More informationExperimental evidence showing that stochastic subspace identication methods may fail 1
Systems & Control Letters 34 (1998) 303 312 Experimental evidence showing that stochastic subspace identication methods may fail 1 Anders Dahlen, Anders Lindquist, Jorge Mari Division of Optimization and
More informationwhich the theoretical results were obtained. Image Signals The geometry of visual scenes under perspective projection generally yields complex image s
NonFourier Motion in The Fourier Spectrum Steven S. Beauchemin GRASP Laboratory Department of Computer and Information Science University of Pennsylvania Philadelphia PA 191048 Abstract Typically, image
More informationIntroduction to Spectral and Time-Spectral Analysis with some Applications
Introduction to Spectral and Time-Spectral Analysis with some Applications A.INTRODUCTION Consider the following process X t = U cos[(=3)t] + V sin[(=3)t] Where *E[U] = E[V ] = 0 *E[UV ] = 0 *V ar(u) =
More informationA general theory of discrete ltering. for LES in complex geometry. By Oleg V. Vasilyev AND Thomas S. Lund
Center for Turbulence Research Annual Research Briefs 997 67 A general theory of discrete ltering for ES in complex geometry By Oleg V. Vasilyev AND Thomas S. und. Motivation and objectives In large eddy
More informationDigital Image Processing Lectures 13 & 14
Lectures 13 & 14, Professor Department of Electrical and Computer Engineering Colorado State University Spring 2013 Properties of KL Transform The KL transform has many desirable properties which makes
More informationTowards a Mathematical Theory of Super-resolution
Towards a Mathematical Theory of Super-resolution Carlos Fernandez-Granda www.stanford.edu/~cfgranda/ Information Theory Forum, Information Systems Laboratory, Stanford 10/18/2013 Acknowledgements This
More informationMATTHIAS P ATZOLD, MEMBER, IEEE, ULRICH KILLAT, MEMBER, IEEE, FRANK LAUE, YINGCHUN LI. Technical University of Hamburg-Harburg.
On the Statistical Properties of Deterministic Simulation Models for Mobile Fading Channels MATTHIAS P ATZOLD, MEMBER, IEEE, ULRICH KILLAT, MEMBER, IEEE, FRANK LAUE, YINGCHUN LI Technical University of
More informationIntroduction to time-frequency analysis. From linear to energy-based representations
Introduction to time-frequency analysis. From linear to energy-based representations Rosario Ceravolo Politecnico di Torino Dep. Structural Engineering UNIVERSITA DI TRENTO Course on «Identification and
More informationContents. 0.1 Notation... 3
Contents 0.1 Notation........................................ 3 1 A Short Course on Frame Theory 4 1.1 Examples of Signal Expansions............................ 4 1.2 Signal Expansions in Finite-Dimensional
More informationMathematical Institute, University of Utrecht. The problem of estimating the mean of an observed Gaussian innite-dimensional vector
On Minimax Filtering over Ellipsoids Eduard N. Belitser and Boris Y. Levit Mathematical Institute, University of Utrecht Budapestlaan 6, 3584 CD Utrecht, The Netherlands The problem of estimating the mean
More informationThese outputs can be written in a more convenient form: with y(i) = Hc m (i) n(i) y(i) = (y(i); ; y K (i)) T ; c m (i) = (c m (i); ; c m K(i)) T and n
Binary Codes for synchronous DS-CDMA Stefan Bruck, Ulrich Sorger Institute for Network- and Signal Theory Darmstadt University of Technology Merckstr. 25, 6428 Darmstadt, Germany Tel.: 49 65 629, Fax:
More informationPerformance Comparison of Two Implementations of the Leaky. LMS Adaptive Filter. Scott C. Douglas. University of Utah. Salt Lake City, Utah 84112
Performance Comparison of Two Implementations of the Leaky LMS Adaptive Filter Scott C. Douglas Department of Electrical Engineering University of Utah Salt Lake City, Utah 8411 Abstract{ The leaky LMS
More informationExpressions for the covariance matrix of covariance data
Expressions for the covariance matrix of covariance data Torsten Söderström Division of Systems and Control, Department of Information Technology, Uppsala University, P O Box 337, SE-7505 Uppsala, Sweden
More informationOn the Space-Varying Filtering
On the Space-Varying Filtering LJubiša Stanković, Srdjan Stanković, Igor Djurović 2 Abstract Filtering of two-dimensional noisy signals is considered in the paper. Concept of nonstationary space-varying
More information= w 2. w 1. B j. A j. C + j1j2
Local Minima and Plateaus in Multilayer Neural Networks Kenji Fukumizu and Shun-ichi Amari Brain Science Institute, RIKEN Hirosawa 2-, Wako, Saitama 35-098, Japan E-mail: ffuku, amarig@brain.riken.go.jp
More informationComparative study of different techniques for Time-Frequency Analysis
Comparative study of different techniques for Time-Frequency Analysis A.Vishwadhar M.Tech Student Malla Reddy Institute Of Technology And Science,Maisammaguda, Dulapally, Secunderabad. Abstract-The paper
More information2.1 Basic Concepts Basic operations on signals Classication of signals
Haberle³me Sistemlerine Giri³ (ELE 361) 9 Eylül 2017 TOBB Ekonomi ve Teknoloji Üniversitesi, Güz 2017-18 Dr. A. Melda Yüksel Turgut & Tolga Girici Lecture Notes Chapter 2 Signals and Linear Systems 2.1
More informationHexagonal QMF Banks and Wavelets. A chapter within Time-Frequency and Wavelet Transforms in Biomedical Engineering,
Hexagonal QMF Banks and Wavelets A chapter within Time-Frequency and Wavelet Transforms in Biomedical Engineering, M. Akay (Editor), New York, NY: IEEE Press, 99. Sergio Schuler and Andrew Laine Computer
More informationAbstract. Jacobi curves are far going generalizations of the spaces of \Jacobi
Principal Invariants of Jacobi Curves Andrei Agrachev 1 and Igor Zelenko 2 1 S.I.S.S.A., Via Beirut 2-4, 34013 Trieste, Italy and Steklov Mathematical Institute, ul. Gubkina 8, 117966 Moscow, Russia; email:
More informationLecture Notes 5: Multiresolution Analysis
Optimization-based data analysis Fall 2017 Lecture Notes 5: Multiresolution Analysis 1 Frames A frame is a generalization of an orthonormal basis. The inner products between the vectors in a frame and
More informationBand-limited Signal Concentration in Time-frequency
Band-limited Signal Concentration in Time-frequency Liying Wei, Rodney A. Kennedy and Tharaka A. Lamahewa College of Engineering and Computer Science, The Australian National University, Canberra, Australia.
More informationA Strict Stability Limit for Adaptive Gradient Type Algorithms
c 009 IEEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional A Strict Stability Limit for Adaptive Gradient Type Algorithms
More informationAtomic decompositions of square-integrable functions
Atomic decompositions of square-integrable functions Jordy van Velthoven Abstract This report serves as a survey for the discrete expansion of square-integrable functions of one real variable on an interval
More informationDenoising Gabor Transforms
1 Denoising Gabor Transforms James S. Walker Abstract We describe denoising one-dimensional signals by thresholding Blackman windowed Gabor transforms. This method is compared with Gauss-windowed Gabor
More informationCover page. : On-line damage identication using model based orthonormal. functions. Author : Raymond A. de Callafon
Cover page Title : On-line damage identication using model based orthonormal functions Author : Raymond A. de Callafon ABSTRACT In this paper, a new on-line damage identication method is proposed for monitoring
More information6.435, System Identification
System Identification 6.435 SET 3 Nonparametric Identification Munther A. Dahleh 1 Nonparametric Methods for System ID Time domain methods Impulse response Step response Correlation analysis / time Frequency
More informationEstimation of the Optimum Rotational Parameter for the Fractional Fourier Transform Using Domain Decomposition
Estimation of the Optimum Rotational Parameter for the Fractional Fourier Transform Using Domain Decomposition Seema Sud 1 1 The Aerospace Corporation, 4851 Stonecroft Blvd. Chantilly, VA 20151 Abstract
More informationGabor Deconvolution. Gary Margrave and Michael Lamoureux
Gabor Deconvolution Gary Margrave and Michael Lamoureux = Outline The Gabor idea The continuous Gabor transform The discrete Gabor transform A nonstationary trace model Gabor deconvolution Examples = Gabor
More informationA Banach Gelfand Triple Framework for Regularization and App
A Banach Gelfand Triple Framework for Regularization and Hans G. Feichtinger 1 hans.feichtinger@univie.ac.at December 5, 2008 1 Work partially supported by EUCETIFA and MOHAWI Hans G. Feichtinger hans.feichtinger@univie.ac.at
More informationFrom Fractional Brownian Motion to Multifractional Brownian Motion
From Fractional Brownian Motion to Multifractional Brownian Motion Antoine Ayache USTL (Lille) Antoine.Ayache@math.univ-lille1.fr Cassino December 2010 A.Ayache (USTL) From FBM to MBM Cassino December
More informationShift Covariant Time-Frequency Distributions of Discrete. Signals. Jerey C. O'Neill. of the requirements for the degree of. Doctor of Philosophy
Shift Covariant Time-Frequency Distributions of Discrete Signals by Jerey C. O'Neill A dissertation submitted in partial fulllment of the requirements for the degree of Doctor of Philosophy (Electrical
More informationLecture 3 Kernel properties and design in Cohen s class time-frequency distributions
Lecture 3 Kernel properties and design in Cohen s class time-frequency distributions Time-frequency analysis, adaptive filtering and source separation José Biurrun Manresa 22.02.2011 Time-Frequency representations
More informationIEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 56, NO. 1, JANUARY
IEEE TRANSACTIONS ON INFORMATION THEORY, VOL 56, NO 1, JANUARY 2010 367 Noncoherent Capacity of Underspread Fading Channels Giuseppe Durisi, Member, IEEE, Ulrich G Schuster, Member, IEEE, Helmut Bölcskei,
More informationSIGNAL AND IMAGE RESTORATION: SOLVING
1 / 55 SIGNAL AND IMAGE RESTORATION: SOLVING ILL-POSED INVERSE PROBLEMS - ESTIMATING PARAMETERS Rosemary Renaut http://math.asu.edu/ rosie CORNELL MAY 10, 2013 2 / 55 Outline Background Parameter Estimation
More informationGabor multipliers for pseudodifferential operators and wavefield propagators
Gabor multipliers for pseudodifferential operators and wavefield propagators Michael P. Lamoureux, Gary F. Margrave and Peter C. Gibson ABSTRACT Gabor multipliers We summarize some applications of Gabor
More informationImproved system blind identification based on second-order cyclostationary statistics: A group delay approach
SaÅdhanaÅ, Vol. 25, Part 2, April 2000, pp. 85±96. # Printed in India Improved system blind identification based on second-order cyclostationary statistics: A group delay approach P V S GIRIDHAR 1 and
More informationEFFECTS OF ILL-CONDITIONED DATA ON LEAST SQUARES ADAPTIVE FILTERS. Gary A. Ybarra and S.T. Alexander
EFFECTS OF ILL-CONDITIONED DATA ON LEAST SQUARES ADAPTIVE FILTERS Gary A. Ybarra and S.T. Alexander Center for Communications and Signal Processing Electrical and Computer Engineering Department North
More informationA REDUCED MULTIPLE GABOR FRAME FOR LOCAL TIME ADAPTATION OF THE SPECTROGRAM
A REDUCED MULTIPLE GABOR FRAME FOR LOCAL TIME ADAPTATION OF THE SPECTROGRAM Marco Liuni, Università di Firenze, Dip. di Matematica U. Dini, Viale Morgagni, 67/a - 5034 Florence - ITALY IRCAM - CNRS STMS,
More informationOperator identification and Feichtinger s algebra
SAMPLING THEORY IN SIGNAL AND IMAGE PROCESSING c 23 SAMPLING PUBLISHING Vol. 1, No. 1, Jan. 22, pp. -5 ISSN: 153-6429 Operator identification and Feichtinger s algebra Götz E. Pfander School of Engineering
More informationWavelet Transform. Figure 1: Non stationary signal f(t) = sin(100 t 2 ).
Wavelet Transform Andreas Wichert Department of Informatics INESC-ID / IST - University of Lisboa Portugal andreas.wichert@tecnico.ulisboa.pt September 3, 0 Short Term Fourier Transform Signals whose frequency
More informationChapter 8 The Discrete Fourier Transform
Chapter 8 The Discrete Fourier Transform Introduction Representation of periodic sequences: the discrete Fourier series Properties of the DFS The Fourier transform of periodic signals Sampling the Fourier
More informationsize within a local area. Using the image as a matched lter, a nite number of scales are searched within a local window to nd the wavelet coecients th
In Proceeding of the 3rd International Workshop on Digital Mammography Chicago, U.S.A. 1996, pp. 447-450 Segmentation of Masses Using Continuous Scale Representations Andrew Laine a, Walter Huda b, Dongwei
More informationN.G.Bean, D.A.Green and P.G.Taylor. University of Adelaide. Adelaide. Abstract. process of an MMPP/M/1 queue is not a MAP unless the queue is a
WHEN IS A MAP POISSON N.G.Bean, D.A.Green and P.G.Taylor Department of Applied Mathematics University of Adelaide Adelaide 55 Abstract In a recent paper, Olivier and Walrand (994) claimed that the departure
More information