Computations of Discrete Gabor Coefficients. Sigang Qiu
|
|
- Pamela Mills
- 6 years ago
- Views:
Transcription
1 !fie<?r ~oz ~~ G\.~t')M4'ffJ Characterizations of Gabor-Gram Matrices and Efficient Computations of Discrete Gabor Coefficients Sigang Qiu Department of Mathematics, University of Connecticut V-9, 196 Auditorium Road, Storrs, CT Tel: (203) Fax: (203) ABSTRACT The motivation of this paper is to give a direct way for determining the crucial Gabor coefficients. Based on the characterized Gabor-Gram matrix structures, we are able to propose fast and direct computational algorithms. We derive also a new way for dual Gabor atoms. 1 Introduction and Preliminaries The theory of Gabor transform is concerned with the main problem of representing a given signal into a series of building blocks (the Gabor elementary functions) and obtained from a basic one (Gabor atom or Gabor wavelet) by appropriate time-frequency shifts along a TF-Iattice. The importance of the Gabor transform is that Gabor coefficients reveal the frequency localization of a signal (or image) very well. Since the Gabor elementary functions are not orthogonal to each other, it is usually very complicated and difficult to determine the Gabor coefficients. By using a so-called bi-orthogonal function,,1 many practical algorithms have been presented in the recent studies on the subject3,7-11,18,19 mostly via the Zak transform and matrix decomposition methods. However, this, plays only an auxiliary role. Generally, we are only required to determine the Gabor coefficients, and do not need to know what,looks like. Our main motivation in this paper is to give a directly iterative method for the computation of Gabor coefficients. In the following, we fix some notations and introduce the concept of Gabor-Gram matrices. In Section 2, we present some equivalent formulations of the bi-orthogonality. Section 3 is the main part of this paper. We present structures of Gabor-Gram matrices and conjugate gradient (CG)-algorithms for calculations of Gabor coefficients and dual gabor atoms.
2 1.1 Notations In order to summarize the results, we have to fix some notations. Essentially, we will keep using the natations in our paper.14 For the self-containedness, we recall some of them. We consider signals as N-periodic row-vectors in ~N. The inner product of two signals is given as (m, y) = L~=-;/ x(k)y(k). We identify linear mapping from the right: m 1-+ m * A. "*,, is used for matrix multiplication. For a matrix A, we use A'to denote the conjugated and transposed matrix of A. We use F and F-1 to denote the discrete Fourier transform (DFT) and the inverse discrete Fourier transform (IDFT), respectively. They are performed by the efficient fast Fourier transform (FFT) and the inverse fast Fourier transform (IFFT). For a signal x E ~N, F(m) = m * FN; and F-1(m) = m * FN, where FN =" (e~) is the Fourier matrix of order N. We also use a: to denote F(x). NxN For a given Gabor atom 9 and lattice constants (a, b), we call (ii, b) := (~, If) the dual lattice constants of (a, b). A discrete Gabor representation of a signal x E ~N with respect to a Gabor triple (g, a, b) is of the form X 0.-1 b-1 = L: L: cnmgnm, n=om=o (1) where gnm = MmbTnag. {Cn,m}n,m are called Gabor coefficients. 1.2 Gabor-Gram Matrix Given a Gabor triple (g,a,b), we use GAB(g,a,b) to denote the Gabor basic matrix with the following block form. 9 Mbg GAB(g,a,b) = M(b-1)bg Tag MbTag M(b_1)bTag (2) MmbTnag where M(b_l)bTOi-1)ag MmbTnag = Mmb (g(o - na), g(l- na),..., g(n - 1- na»
3 for m = 0.1,..., b - 1; n = 0,1,..., a - 1, is an N-elements row vector and Ml is the modulation operator. Remark: 1. It is clear that the row vectors of GAB(g.,b) form the Gabor family {gnm}n,m' 2. It is easy to see that GAB(g,,b) and the Gabor basic matrix GAB(g, a, b) introduced in the previous papers12 14 are exactly the same, except for the different orders of the row vectors. Both the row vectors of GAB(g..b) and GAB(g, a, b) form the Gabor family {gnm}' 3. It is not difficult to check that [GAB(g,,b)]' * [GAB(g,.b)] = [GAB(g, a, b)]' * [GAB(g, a, b)]. (3) Therefore, the Gabor matrix corresponding to (g, a, b) can be formulated by Eq. (3) either with GAB(g, a, b) or with GAB(g.,b)' That is why we also call GAB(g..b) a Gabor basic matrix. We use the Gabor basic matrix GAB(g..b) to define the Gabor-Gram matrix. DEFINITION 1.1 (GABOR-GRAM MATRIX). For a Gabor triple (g,a,b), we call the N x N (N = ~:) matrix GM(g..b) = [GAB(g.,b)] * [GAB(g.,bl' the Gabor-Gram matrix associated to (g, a, b). We generally call GM(9i,92,.b) = [GAB(g!,,b)] * [GAB(92.,bl' the Gabor-Gram matrix associated to Gabor triples (g1, a, b) and (g2,a,b). In addition, we use the following rotation operator 0". DEFINITION 1.2. For a given nu.mber r which divides N, we define 0"(', r) as the following rotation operator on (CN. where y = (y(j))f:~1 is given by (j : (CN -+ (CN; X = (X(j));=~1 -+ y, Y ( J') = { x(j + -1) r - 1) if 1=0.. N for J = kr + 1 wzth k = 0, 1,..., and 1 = 0,1,..., r #0 r Definition 1.2 depends only on the coordinates of the given vector x. 0"(" r) is well-defined for both row and column vectors. It is not difficult to check that the inverse of o{, r) is defined as follows. where y = (y(j))f=~l is given by (j-1 : (f}n -+ (CN; X = (x(j));=~1 -+ y, Y (')=={ J x(j-r+l) + 1) if l=r-1 1 # r - 1 for j = kr + 1 with k = 0,1,..., 1'f - 1 and l = 0, 1,..., r - 1.
4 For a matrix A E q;nxq column vector of A. and a number r which divides N, we define o-(a, r) by applying 0-(., r) to each Similarly, we have no difficulty to see that 0--1(A, r) is obtained by applying 0--1 to each column vector of A. 2 Matrix-Formulation of the Bi-orthogonality Using the introduced Gabor basic matrix, we are able to state the Gabor bi-orthogonality18 and the equivalent conditions in matrix forms. THEOREM 1. Given a Gabor triple (g, a, b) and another signal h, the following statements are equivalent: where lr denotes the r x r identity matrix, for rein. ab GM(g,h,b,ii) = GM(h,g,b,ii) = N * lab, [ ]' ab h * GAB(g,b,ii) = N * eab, where er E q;r denotes the 'Unit vector with the first entry being ]' ab 9 * [GAB(h,b,ii) = N * eab GM(h.,U,a, H b) = GM(. 9, h.,a, - b) = N*lab. h.* [GAB(g,Q.,b) ]' = ab * eab. 9 * [GAB(h.,ii,b)]' = ab * eab Moreover, (g, a, b) generates a Gabor frame if and only if there exists a signal h such that one of the above conditions is satisfied. For the higher dimensional cases, if Gabor functions are separable, we have no difficulty to state the results similar to Theorem 1. The last three conditions are based on the fact that (g, a, b) generates a frame if and only if (9, b, a) generates a frame.12
5 3 Gabor-Gram Matrix Characterization We present observations on the Gabor-Gram matrix structures in the first part. This is the basis for algorithms proposed in the second part of this section. 3.1 Gabor-Gram Matrix Structure THEOREM 2. Let us write the Gabor-Gram-matrix GM as the following block-matrix form GM = [B1,B2,, BN], (4) where Bk for k = 1,2,..., IV is the k-th column vector of GM. Set B = [Bl,Bb+1,..., B(ii-l)bH], (5) then the IV x IV matrix GM is completely determined by the IV x a small block matrix B via cyclic rotations operator 0'(" b). p'recisely, if we write then B(s) = [Bs,Bb+S'..., B(ii-l)b+S], for s = 2,...,b, B(s) = (7s-1(B, b). 3.2 Algorithms Gabor-Gram Matrix-Vector (GGMV) Multiplication In this section, we are presenting a useful Gabor-Gram matrix-vector (GGMV') multiplication: For any vector Xc E q:n, we show that Yc = Xc * GM can be obtained by the Gabor-Gram block matrix B only. ALGORITHM 1. Under the assumptions of Theorem 2, let Xc be an IV-entries vector. Then Yc = Xc * GM can be computed by Xc and the Gabor-Gram block matrix B only. Yc = (YcU»f:o is given by Yc(j) = O'-k(xc,b) * B(s), where j = sb + k for k = 0, 1,..., b - 1, s = 0,1,...,a-I, and B(s) for k = 0,1,...,a-I is the s + 1-th column vector of B. Proof. Using the notations of Theorem 2, for j = sb + k with k = 0,1,..., b - 1 and s = 0,1,..., a, Theorem 2 implies that B 81+1<:+1 = O'k (B 8;;+1)' Thus, We finished the proof. 0 Yc(j) = Yc(sb + k) = Xc * B sbh = Xc * O'k(B sbh) = Xc * (7k (B(s»'
6 3.2.2 Gabor Coefficients For an arbitrary Gabor triple (g, a, b), the following CG-algorithm gives a general way for computing the Gabor coefficients. The algorithm holds for oversampling (ab < N), critical sampling (ab = N) and undersampling (ab > N) cases. We do not assume that (g, a, b) generates a frame. In this case, we determine the best approximations by linear combinations of Gabor familif1s, lfi:j""-, ':l~,,{_} ALGORITHM 2 (GABOR COEFFICIENTS). For a signal x E (CN, the Gabor coeficients C = {cnm} with respect to (g, a, b) are determined via the following procedu.res: 1st step: Caculate Xc = X * [GAB(g,a,b)]" 2nd step: Compu.te the Gabor-Gram block matrix B. 3rd step: Obtain the Gabor coefficients C = {cn,m} by solving the linear system Xc = C * GM via the CG-method.6 Remark: Algorithm 2 is based on the fact that if we can get the solution for C from the linear system C * GM = Xc, i.e., C[GAB(g,a,b)J * [GAB(g,a,b)l' = x * GAB(g,a,b)' then it is easy to see that C * [GAB(g,a,b)] * S = x * S. where S = [GAB(g,a,b)l' * [GAB(g,a,b)] is the associated Gabor matrix. If (g, a, b) generates a Gabor frame, since S is invertible, we deduce that x = C * [GAB(g,a,b)J. This is the same as Gabor representation Eq. (1). For technical details, we refer to the subsequent papers. The 1st step and 2nd step are performed efficiently by the short time Fourier transform (STFT). Although GM is generally not invertible (e.g., for the oversampling cases), we show that the standard CG-method6 can be applied to the 3rd step for the undersampled cases. to determine the Gabor coefficients. Algorithm 2 works fast The matrix-vector multiplications involved in the 3rd step are performed by the efficient GGMV-multiplication with the precalculated Gabor-Gram block B Dual Gabor Atom If a given Gabor triple (g,a,b) generates a frame, by Theorem 1, there exists a Gabor signal g (we call it a dual Gabor atom such that 9 * {GAB(g,b,ii)l' = ~ eab. This implies that --1 g = {g * GAB( g, -b,a-) * GM } * [GAB( 9, b-,a-)]
7 ab --1 = {Neab * GM } * [GAB(g,b,a)1 where GM = [GAB(g,b,a)]' * [GAB(g,b,ei)] is the Gabor-Gram matrix corresponding to (g, b, ii). Using the CG-algorithm, we solve the linear equation tl * GM = ~e. Then we obtain h = tl * GAB(g,a,b) by STFT. If (g, a, b) does not generate a frame, we will use 9 = {~eab *pinv(gm)} * [GAB(g,b,ei)] to get the generalized dual Gabor atom (GDGA),13 where pinv(gm) is the pseudo-inverse of GM. 4 Further Results In this section we report some of the further results obtained after our submission of the proposal. For more detailed information, we refer to the subsequent papers. THEOREM 3 (16). Under the assumptions of Theorem 2, if we write GM = (Uk,I)NxN as the following block form, D2 1 D2 2. D2 a GM=.'.' ( Dei,l Dl,l Dei,2 Dl,2... Dei,a D1,ei 1 (6) where Dp,q = (U(P-l)b+l:,(q-l)b+l) b,b then Dp,q are circulant matrices of order b. for p, q = 1,2,...,a. THEOREM 4 (16). Gabor-Gram matrices are simultaneously unitarily similar to banded matrices. Precisely, under the assumptions of Theorem 3, GM = u' * G M * U is a banded sparse matrix with all the possible nonzero entries being in k-th (for k = 0, ±(b - 1),..., ±(b - l)(a - 1») diagonals, where U = Ia 0 Fb, Ia is the a x a identity matrix and Fb is the Fourier matrix of order b. We call GM the unitary equivalent Gabor-Gram matrix (UEGM) associated with Gabor triple (g, a, b). Based on Theorem 3 and Theorem 4, a new fast direct way of determining the Gabor coefficients has been proposed.16 5 Numerical Illustrations We illustrate some of the numerical results done with MATLAB at a SUN4 workstation. Figure 1 illustrates a Gabor atom of signal length N = 560 and the associated dual Gabor atoms calculated by the stated in Section and the under-oversampling relationship established in.is Case (0) is a oversampling
8 a where we take lattice constants pair (a, b) = (28,16). Case (U) is undersampling with lattice constants (a, b) = (40,35). Figure 2 illustrates a Chirp signal and the Gabor coefficients corresponding to Case (0) obtained by Algorithm 2. The reconstruction (relative) error is in order of Figure 3 is for the Case (U). It shows a signal from the Gabor space and the Gabor coefficients. reconstruction (relative) error is also in order of The Figure 4 shows the explicit sparsity of the unitary equivalent Gabor-Gram matrix corresponding to three cases. The case (01) and (02) are oversampling with lattice constants pair (a,b) = (6,6), (8,4). Case (u) is undersampling with (a, b) = (8,8). The signal length N is Conclusion We have introduced concepts of basic Gabor matrix and Gabor-Gram matrices. Simply equivalent conditions for the bi-orthogonality have been formulated. We have developed the structural properties of the Gabor-Gram matrices. A direct algorithm based on the CG-procedure has been presented for the determination of the Gabor coefficients. The algorithm is applied for undersampling as well as critical and oversampling cases. It works even if a given Gabor triple (g, a, b) does not generate a frame. In this case, we are able to use the calculated Gabor coefficients to obtain the best approximation by linear combinations of Gabor families. A new algorithm is also proposed for the dual Gabor atom. In addition, we have summarized some of the latest new results. This enables us to give much efficient computational algorithms. 7 REFERENCES [1] M. J. Bastiaans, Gabor's signal expansion of a signal into Gaussian elementary signals, Proc. IEEE, ,1980. [2] M. J. Bastiaans, Gabor's Signal Expansion and Its Relation to Sampling of the Sliding- Window Spectrum, Advanced Topics in Shannon Sampling and Interpolation Theory (edited by R. J. Marks II), 1-35, [3] J. Benedetto, Gabor representations and wavelets, Commutative Harmonic Analysis (edited by D. Colella), Contemp. Math. 19, 9-27, [4] E. C. Boman, fast Algorithms for Toeplitz Equations, Ph.D thesis, University of Connecticut, [5] C. E. Heil and D. F. Walnut, Continuous and Discrete Wavelet Transforms, SIAM Review 31 (4), , [6] G. Golub and C. van Loan, Matrix Computations, The Johns Hopkins University Press, Baltimore and London, 1993.
9 [7] 1. Daubechies, The wavelet transform, time-frequency localization, and signal an31ysis, IEEE Trans. Inform. Theory, 36, , [8] A. J. E. M. Janssen, Gabor representation of generalized functions, J. Math. Anal. Appl., 83, , [9] A. J. E. M. Janssen, Duality and biorthogonality for Weyl-Heisenberg frames, Journal of Fourier Analysis and Applications, 1(3), [10] R. S. Orr, Tbe Order of Computation for Finite Discrete Gabor Transform, IEEE Trans. on Signal Processing, 41(2), , [11] S. Qian and D. Chen, Discrete Gabor Transform, IEEE Trans. on Signal Processing, 41(7), , [12] S. Qiu and H. G. Feichtinger, Discrete Gabor Structure and Optimal Representation, to appear in IEEE Trans. on Signal Processing, October, [13] S. Qiu, Generalized Dual Gabor Atoms and Best Approximations by Gabor Family, submitted for publication. [14] S. Qiu, Superfast Computations of Dual and Tight Gabor Atoms, SPIE-Proceeding: San Diego'95. [15] S. Qiu, On the Best Approximation by Undersampled Gabor Families, submitted to IEEE Trans. on Signal Processing. [16] S. Qiu, Discrete Gabor Transforms: the Gabor-Gram Matrix Approacb, submitted for publication. [17] M. Unser, Fast Gabor-Like Windowed Fourier and Continuous Wavelet Transforms, IEEE Signal Processing Letters, 1(5), [18] J. Wexler and S. Raz, Discrete Gabor Expansion, Signal Processing, 21(1990), [19] M. Zibulski and Y. Y. Zeevi, Oversampling in tbe Gabor Scbeme, IEEE Trans. on Signal Processing, 41(8), , August, 1993.
10 I 1 Gabor 12 atom dual Gabor atom (0) C) 1 -o.sl O.st (0): (U): oversampling undersampling 2 0.5; generalized dual Gabor atom (U) Figure 1: Gabor atom and the a.gsociated dual Gabor atoms corresponding to a oversampling ca.geand a undersampling case, where lattice constants (a, b) are (28,16) and (40,35). The signal length N is 560. abor coefficients I (0) 20 Of- -11-, ,., v 400, , V V V 10 t ''> } ([l..'n~ <~ "-,.'" T.... Chirp signal --=--.:;::: \ ~ ~ A Figure 2: A Chirp signal and the Gabor coefficients for Case (0). The reconstruction relative error is in order of
11 1 signal from Gabor space o Gabor coefficients (U) Figure 3: A signal from the Gabor space and the Gabor coefficients. It is the undersampled case (U). The reconstruction relative error is in order of I 0.5 (01): (a, b) = (6, 6) o (02): (a, b) =.(8, 4) -0.5 (u): (a, b) = (8, 8) UEGM matrix (02) o ~ '....,, " ', 20 ~ ~ ' ~~ '~"'",., ~ ~ ~....., 40~.,. 60 ~~:::::~:::::::::::::::::::::::::::~::::::::::::::::::: o nz = ~_1 UEGM matrix (01) o nz = UEGM matrix (u) 01_. ~ ~ ~ ~ ~ _. _. _. _. a_. e a t -Cl. -_. - -_. e - -. a_. a_. _ _.. 20~ _.. e. a _. _. _ 30 te -. o nz = 216 Figure 4: The unitarily equivalent matrices structures of the Gabor-Gram matrices corresponding to three cases. (01) and (02) are oversampling and (u) is undersampling.
Approximately 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 informationHyperbolic Secants Yield Gabor Frames
Applied and Computational Harmonic Analysis 1, 59 67 ( doi:1.16/acha.1.376, available online at http://www.idealibrary.com on Hyperbolic Secants Yield Gabor Frames A. J. E. M. Janssen Philips Research
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 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 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 informationGeneralized shift-invariant systems and frames for subspaces
The Journal of Fourier Analysis and Applications Generalized shift-invariant systems and frames for subspaces Ole Christensen and Yonina C. Eldar ABSTRACT. Let T k denote translation by k Z d. Given countable
More informationHyperbolic secants yield Gabor frames
arxiv:math/334v [math.fa] 3 Jan 3 Hyperbolic secants yield Gabor frames A.J.E.M. Janssen and Thomas Strohmer Abstract We show that (g,a,b) is a Gabor frame when a >,b >,ab < and g (t) = ( πγ) (coshπγt)
More informationGABOR FRAMES AND OPERATOR ALGEBRAS
GABOR FRAMES AND OPERATOR ALGEBRAS J-P Gabardo a, Deguang Han a, David R Larson b a Dept of Math & Statistics, McMaster University, Hamilton, Canada b Dept of Mathematics, Texas A&M University, College
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 informationGabor Frames. Karlheinz Gröchenig. Faculty of Mathematics, University of Vienna.
Gabor Frames Karlheinz Gröchenig Faculty of Mathematics, University of Vienna http://homepage.univie.ac.at/karlheinz.groechenig/ HIM Bonn, January 2016 Karlheinz Gröchenig (Vienna) Gabor Frames and their
More informationIntroduction to Hilbert Space Frames
to Hilbert Space Frames May 15, 2009 to Hilbert Space Frames What is a frame? Motivation Coefficient Representations The Frame Condition Bases A linearly dependent frame An infinite dimensional frame Reconstructing
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 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 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 informationSOME TOPICS ON WAVELETS
2 SOME TOPICS ON WAVELETS RYUICHI ASHINO 1. Introduction Let us consider music. If the music is recorded from a live broadcast onto tape, we have a signal, that is, a function f(x). The time-frequency
More informationBANACH FRAMES GENERATED BY COMPACT OPERATORS ASSOCIATED WITH A BOUNDARY VALUE PROBLEM
TWMS J. Pure Appl. Math., V.6, N.2, 205, pp.254-258 BRIEF PAPER BANACH FRAMES GENERATED BY COMPACT OPERATORS ASSOCIATED WITH A BOUNDARY VALUE PROBLEM L.K. VASHISHT Abstract. In this paper we give a type
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 informationPhase Retrieval from Local Measurements: Deterministic Measurement Constructions and Efficient Recovery Algorithms
Phase Retrieval from Local Measurements: Deterministic Measurement Constructions and Efficient Recovery Algorithms Aditya Viswanathan Department of Mathematics and Statistics adityavv@umich.edu www-personal.umich.edu/
More informationSo reconstruction requires inverting the frame operator which is often difficult or impossible in practice. It follows that for all ϕ H we have
CONSTRUCTING INFINITE TIGHT FRAMES PETER G. CASAZZA, MATT FICKUS, MANUEL LEON AND JANET C. TREMAIN Abstract. For finite and infinite dimensional Hilbert spaces H we classify the sequences of positive real
More informationA Novel Fast Computing Method for Framelet Coefficients
American Journal of Applied Sciences 5 (11): 15-157, 008 ISSN 1546-939 008 Science Publications A Novel Fast Computing Method for Framelet Coefficients Hadeel N. Al-Taai Department of Electrical and Electronic
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 informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationDecomposition of Riesz frames and wavelets into a finite union of linearly independent sets
Decomposition of Riesz frames and wavelets into a finite union of linearly independent sets Ole Christensen, Alexander M. Lindner Abstract We characterize Riesz frames and prove that every Riesz frame
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 1: Course Overview & Matrix-Vector Multiplication Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 20 Outline 1 Course
More informationarxiv: v1 [math.na] 1 Sep 2018
On the perturbation of an L -orthogonal projection Xuefeng Xu arxiv:18090000v1 [mathna] 1 Sep 018 September 5 018 Abstract The L -orthogonal projection is an important mathematical tool in scientific computing
More informationPOINT VALUES AND NORMALIZATION OF TWO-DIRECTION MULTIWAVELETS AND THEIR DERIVATIVES
November 1, 1 POINT VALUES AND NORMALIZATION OF TWO-DIRECTION MULTIWAVELETS AND THEIR DERIVATIVES FRITZ KEINERT AND SOON-GEOL KWON,1 Abstract Two-direction multiscaling functions φ and two-direction multiwavelets
More informationDORIN ERVIN DUTKAY AND PALLE JORGENSEN. (Communicated by )
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9939(XX)0000-0 OVERSAMPLING GENERATES SUPER-WAVELETS arxiv:math/0511399v1 [math.fa] 16 Nov 2005 DORIN ERVIN DUTKAY
More informationA fast randomized algorithm for overdetermined linear least-squares regression
A fast randomized algorithm for overdetermined linear least-squares regression Vladimir Rokhlin and Mark Tygert Technical Report YALEU/DCS/TR-1403 April 28, 2008 Abstract We introduce a randomized algorithm
More informationA Method for Constructing Diagonally Dominant Preconditioners based on Jacobi Rotations
A Method for Constructing Diagonally Dominant Preconditioners based on Jacobi Rotations Jin Yun Yuan Plamen Y. Yalamov Abstract A method is presented to make a given matrix strictly diagonally dominant
More informationON SAMPLING RELATED PROPERTIES OF B-SPLINE RIESZ SEQUENCES
ON SAMPLING RELATED PROPERTIES OF B-SPLINE RIESZ SEQUENCES SHIDONG LI, ZHENGQING TONG AND DUNYAN YAN Abstract. For B-splineRiesz sequencesubspacesx span{β k ( n) : n Z}, thereis an exact sampling formula
More informationFrame Diagonalization of Matrices
Frame Diagonalization of Matrices Fumiko Futamura Mathematics and Computer Science Department Southwestern University 00 E University Ave Georgetown, Texas 78626 U.S.A. Phone: + (52) 863-98 Fax: + (52)
More informationFRAME DUALITY PROPERTIES FOR PROJECTIVE UNITARY REPRESENTATIONS
FRAME DUALITY PROPERTIES FOR PROJECTIVE UNITARY REPRESENTATIONS DEGUANG HAN AND DAVID LARSON Abstract. Let π be a projective unitary representation of a countable group G on a separable Hilbert space H.
More informationFinite Frame Quantization
Finite Frame Quantization Liam Fowl University of Maryland August 21, 2018 1 / 38 Overview 1 Motivation 2 Background 3 PCM 4 First order Σ quantization 5 Higher order Σ quantization 6 Alternative Dual
More informationThe Lanczos and conjugate gradient algorithms
The Lanczos and conjugate gradient algorithms Gérard MEURANT October, 2008 1 The Lanczos algorithm 2 The Lanczos algorithm in finite precision 3 The nonsymmetric Lanczos algorithm 4 The Golub Kahan bidiagonalization
More informationContents. Acknowledgments
Table of Preface Acknowledgments Notation page xii xx xxi 1 Signals and systems 1 1.1 Continuous and discrete signals 1 1.2 Unit step and nascent delta functions 4 1.3 Relationship between complex exponentials
More information5. Orthogonal matrices
L Vandenberghe EE133A (Spring 2017) 5 Orthogonal matrices matrices with orthonormal columns orthogonal matrices tall matrices with orthonormal columns complex matrices with orthonormal columns 5-1 Orthonormal
More informationTransforms and Orthogonal Bases
Orthogonal Bases Transforms and Orthogonal Bases We now turn back to linear algebra to understand transforms, which map signals between different domains Recall that signals can be interpreted as vectors
More informationFrames and operator representations of frames
Frames and operator representations of frames Ole Christensen Joint work with Marzieh Hasannasab HATA DTU DTU Compute, Technical University of Denmark HATA: Harmonic Analysis - Theory and Applications
More informationClarkson Inequalities With Several Operators
isid/ms/2003/23 August 14, 2003 http://www.isid.ac.in/ statmath/eprints Clarkson Inequalities With Several Operators Rajendra Bhatia Fuad Kittaneh Indian Statistical Institute, Delhi Centre 7, SJSS Marg,
More informationAPPROXIMATING THE INVERSE FRAME OPERATOR FROM LOCALIZED FRAMES
APPROXIMATING THE INVERSE FRAME OPERATOR FROM LOCALIZED FRAMES GUOHUI SONG AND ANNE GELB Abstract. This investigation seeks to establish the practicality of numerical frame approximations. Specifically,
More informationNotes on basis changes and matrix diagonalization
Notes on basis changes and matrix diagonalization Howard E Haber Santa Cruz Institute for Particle Physics, University of California, Santa Cruz, CA 95064 April 17, 2017 1 Coordinates of vectors and matrix
More informationLinear Independence of Finite Gabor Systems
Linear Independence of Finite Gabor Systems 1 Linear Independence of Finite Gabor Systems School of Mathematics Korea Institute for Advanced Study Linear Independence of Finite Gabor Systems 2 Short trip
More informationThe Theory of Wavelets with Composite Dilations
The Theory of Wavelets with Composite Dilations Kanghui Guo 1, Demetrio Labate 2, Wang Q Lim 3, Guido Weiss 4, and Edward Wilson 5 1 Department of Mathematics, Southwest Missouri State University, Springfield,
More informationApplied and Computational Harmonic Analysis
Appl. Comput. Harmon. Anal. 32 (2012) 139 144 Contents lists available at ScienceDirect Applied and Computational Harmonic Analysis www.elsevier.com/locate/acha Letter to the Editor Frames for operators
More information1 Digital Signal Processing Template
This is the file: http://www.univie.ac.at/nuhag/feicours/ss/digsigbook.pdf Digital Signal Processing Template Summary. Keywords. test. Subsection.2 Subsection 2.3 MATLAB CODE VERBATIM PART.4 Exercises.
More informationThe Density Theorem and the Homogeneous Approximation Property for Gabor Frames
The Density Theorem and the Homogeneous Approximation Property for Gabor Frames Christopher Heil School of Mathematics, Georgia Institute of Technology, Atlanta, GA 30332 USA heil@math.gatech.edu Summary.
More informationOperators with Closed Range, Pseudo-Inverses, and Perturbation of Frames for a Subspace
Canad. Math. Bull. Vol. 42 (1), 1999 pp. 37 45 Operators with Closed Range, Pseudo-Inverses, and Perturbation of Frames for a Subspace Ole Christensen Abstract. Recent work of Ding and Huang shows that
More informationInverses. Stephen Boyd. EE103 Stanford University. October 28, 2017
Inverses Stephen Boyd EE103 Stanford University October 28, 2017 Outline Left and right inverses Inverse Solving linear equations Examples Pseudo-inverse Left and right inverses 2 Left inverses a number
More informationMS&E 318 (CME 338) Large-Scale Numerical Optimization
Stanford University, Management Science & Engineering (and ICME MS&E 38 (CME 338 Large-Scale Numerical Optimization Course description Instructor: Michael Saunders Spring 28 Notes : Review The course teaches
More informationMath 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination
Math 0, Winter 07 Final Exam Review Chapter. Matrices and Gaussian Elimination { x + x =,. Different forms of a system of linear equations. Example: The x + 4x = 4. [ ] [ ] [ ] vector form (or the column
More informationWavelets: Theory and Applications. Somdatt Sharma
Wavelets: Theory and Applications Somdatt Sharma Department of Mathematics, Central University of Jammu, Jammu and Kashmir, India Email:somdattjammu@gmail.com Contents I 1 Representation of Functions 2
More informationA NEW IDENTITY FOR PARSEVAL FRAMES
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9939(XX)0000-0 A NEW IDENTITY FOR PARSEVAL FRAMES RADU BALAN, PETER G. CASAZZA, DAN EDIDIN, AND GITTA KUTYNIOK
More informationValidity of WH-Frame Bound Conditions Depends on Lattice Parameters
Applied and Computational Harmonic Analysis 8, 104 112 (2000) doi:10.1006/acha.1999.0281, available online at http://www.idealibrary.com on Validity of WH-Frame Bound Conditions Depends on Lattice Parameters
More informationON WEIGHTED PARTIAL ORDERINGS ON THE SET OF RECTANGULAR COMPLEX MATRICES
ON WEIGHTED PARTIAL ORDERINGS ON THE SET OF RECTANGULAR COMPLEX MATRICES HANYU LI, HU YANG College of Mathematics and Physics Chongqing University Chongqing, 400030, P.R. China EMail: lihy.hy@gmail.com,
More informationDUALITY PRINCIPLE IN g-frames
Palestine Journal of Mathematics Vol. 6(2)(2017), 403 411 Palestine Polytechnic University-PPU 2017 DUAITY PRINCIPE IN g-frames Amir Khosravi and Farkhondeh Takhteh Communicated by Akram Aldroubi MSC 2010
More informationOn Riesz-Fischer sequences and lower frame bounds
On Riesz-Fischer sequences and lower frame bounds P. Casazza, O. Christensen, S. Li, A. Lindner Abstract We investigate the consequences of the lower frame condition and the lower Riesz basis condition
More information1 Singular Value Decomposition
1 Singular Value Decomposition Factorisation of rectangular matrix (generalisation of eigenvalue concept / spectral theorem): For every matrix A C m n there exists a factorisation A = UΣV U C m m, V C
More informationShort-time Fourier transform for quaternionic signals
Short-time Fourier transform for quaternionic signals Joint work with Y. Fu and U. Kähler P. Cerejeiras Departamento de Matemática Universidade de Aveiro pceres@ua.pt New Trends and Directions in Harmonic
More informationA short introduction to frames, Gabor systems, and wavelet systems
Downloaded from orbit.dtu.dk on: Mar 04, 2018 A short introduction to frames, Gabor systems, and wavelet systems Christensen, Ole Published in: Azerbaijan Journal of Mathematics Publication date: 2014
More informationSymmetric Wavelet Tight Frames with Two Generators
Symmetric Wavelet Tight Frames with Two Generators Ivan W. Selesnick Electrical and Computer Engineering Polytechnic University 6 Metrotech Center, Brooklyn, NY 11201, USA tel: 718 260-3416, fax: 718 260-3906
More informationbuer overlfows at intermediate nodes in the network. So to most users, the behavior of a packet network is not characterized by random loss, but by un
Uniform tight frames for signal processing and communication Peter G. Casazza Department of Mathematics University of Missouri-Columbia Columbia, MO 65211 pete@math.missouri.edu Jelena Kovacevic Bell Labs
More informationM.A. Botchev. September 5, 2014
Rome-Moscow school of Matrix Methods and Applied Linear Algebra 2014 A short introduction to Krylov subspaces for linear systems, matrix functions and inexact Newton methods. Plan and exercises. M.A. Botchev
More informationAn Introduction to HILBERT-HUANG TRANSFORM and EMPIRICAL MODE DECOMPOSITION (HHT-EMD) Advanced Structural Dynamics (CE 20162)
An Introduction to HILBERT-HUANG TRANSFORM and EMPIRICAL MODE DECOMPOSITION (HHT-EMD) Advanced Structural Dynamics (CE 20162) M. Ahmadizadeh, PhD, PE O. Hemmati 1 Contents Scope and Goals Review on transformations
More informationFrom Bernstein approximation to Zauner s conjecture
From Bernstein approximation to Zauner s conjecture Shayne Waldron Mathematics Department, University of Auckland December 5, 2017 Shayne Waldron (University of Auckland) Workshop on Spline Approximation
More informationDouble Preconditioning for Gabor Frames
IEEE TRANSACTIONS ON SIGNAL PROCESSING 1 Double Preconditioning for Gabor Frames Peter Balazs*, Member, IEEE, Hans G. Feichtinger, Mario Hampejs and Günther Kracher Abstract We present an application of
More informationON THE HÖLDER CONTINUITY OF MATRIX FUNCTIONS FOR NORMAL MATRICES
Volume 10 (2009), Issue 4, Article 91, 5 pp. ON THE HÖLDER CONTINUITY O MATRIX UNCTIONS OR NORMAL MATRICES THOMAS P. WIHLER MATHEMATICS INSTITUTE UNIVERSITY O BERN SIDLERSTRASSE 5, CH-3012 BERN SWITZERLAND.
More informationFast Angular Synchronization for Phase Retrieval via Incomplete Information
Fast Angular Synchronization for Phase Retrieval via Incomplete Information Aditya Viswanathan a and Mark Iwen b a Department of Mathematics, Michigan State University; b Department of Mathematics & Department
More informationMain matrix factorizations
Main matrix factorizations A P L U P permutation matrix, L lower triangular, U upper triangular Key use: Solve square linear system Ax b. A Q R Q unitary, R upper triangular Key use: Solve square or overdetrmined
More informationA fast algorithm of two-level banded Toeplitz systems of linear equations with application to image restoration
NTMSCI 5, No. 2, 277-283 (2017) 277 New Trends in Mathematical Sciences http://dx.doi.org/ A fast algorithm of two-level banded Toeplitz systems of linear equations with application to image restoration
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences)
AMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences) Lecture 1: Course Overview; Matrix Multiplication Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical
More informationMATH 425-Spring 2010 HOMEWORK ASSIGNMENTS
MATH 425-Spring 2010 HOMEWORK ASSIGNMENTS Instructor: Shmuel Friedland Department of Mathematics, Statistics and Computer Science email: friedlan@uic.edu Last update April 18, 2010 1 HOMEWORK ASSIGNMENT
More informationON WEIGHTED PARTIAL ORDERINGS ON THE SET OF RECTANGULAR COMPLEX MATRICES
olume 10 2009, Issue 2, Article 41, 10 pp. ON WEIGHTED PARTIAL ORDERINGS ON THE SET OF RECTANGULAR COMPLEX MATRICES HANYU LI, HU YANG, AND HUA SHAO COLLEGE OF MATHEMATICS AND PHYSICS CHONGQING UNIERSITY
More informationThe inverse of a tridiagonal matrix
Linear Algebra and its Applications 325 (2001) 109 139 www.elsevier.com/locate/laa The inverse of a tridiagonal matrix Ranjan K. Mallik Department of Electrical Engineering, Indian Institute of Technology,
More informationA Non-Linear Lower Bound for Constant Depth Arithmetical Circuits via the Discrete Uncertainty Principle
A Non-Linear Lower Bound for Constant Depth Arithmetical Circuits via the Discrete Uncertainty Principle Maurice J. Jansen Kenneth W.Regan November 10, 2006 Abstract We prove a non-linear lower bound on
More informationON ORTHOGONAL REDUCTION TO HESSENBERG FORM WITH SMALL BANDWIDTH
ON ORTHOGONAL REDUCTION TO HESSENBERG FORM WITH SMALL BANDWIDTH V. FABER, J. LIESEN, AND P. TICHÝ Abstract. Numerous algorithms in numerical linear algebra are based on the reduction of a given matrix
More informationORTHONORMAL SAMPLING FUNCTIONS
ORTHONORMAL SAMPLING FUNCTIONS N. KAIBLINGER AND W. R. MADYCH Abstract. We investigate functions φ(x) whose translates {φ(x k)}, where k runs through the integer lattice Z, provide a system of orthonormal
More informationFinite Frames and Graph Theoretical Uncertainty Principles
Finite Frames and Graph Theoretical Uncertainty Principles (pkoprows@math.umd.edu) University of Maryland - College Park April 13, 2015 Outline 1 Motivation 2 Definitions 3 Results Outline 1 Motivation
More informationOn Frame Wavelet Sets and Some Related Topics
On Frame Wavelet Sets and Some Related Topics Xingde Dai and Yuanan Diao Abstract. A special type of frame wavelets in L 2 (R) or L 2 (R d ) consists of those whose Fourier transforms are defined by set
More informationFilter Bank Frame Expansions With Erasures
IEEE TRANSACTIONS ON INFORMATION THEORY, VOL 48, NO 6, JUNE 2002 1439 Filter Bank Frame Expansions With Erasures Jelena Kovačević, Fellow, IEEE, Pier Luigi Dragotti, Student Member, IEEE, and Vivek K Goyal,
More informationRepresentation of doubly infinite matrices as non-commutative Laurent series
Spec. Matrices 217; 5:25 257 Research Article Open Access María Ivonne Arenas-Herrera and Luis Verde-Star* Representation of doubly infinite matrices as non-commutative Laurent series https://doi.org/1.1515/spma-217-18
More informationG-frames in Hilbert Modules Over Pro-C*-algebras
Available online at http://ijim.srbiau.ac.ir/ Int. J. Industrial Mathematics (ISSN 2008-5621) Vol. 9, No. 4, 2017 Article ID IJIM-00744, 9 pages Research Article G-frames in Hilbert Modules Over Pro-C*-algebras
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 informationSTABILITY OF INVARIANT SUBSPACES OF COMMUTING MATRICES We obtain some further results for pairs of commuting matrices. We show that a pair of commutin
On the stability of invariant subspaces of commuting matrices Tomaz Kosir and Bor Plestenjak September 18, 001 Abstract We study the stability of (joint) invariant subspaces of a nite set of commuting
More informationDUALITY FOR FRAMES ZHITAO FAN, ANDREAS HEINECKE, AND ZUOWEI SHEN
DUALITY FOR FRAMES ZHITAO FAN, ANDREAS HEINECKE, AND ZUOWEI SHEN Abstract. The subject of this article is the duality principle, which, well beyond its stand at the heart of Gabor analysis, is a universal
More informationMath Linear Algebra II. 1. Inner Products and Norms
Math 342 - Linear Algebra II Notes 1. Inner Products and Norms One knows from a basic introduction to vectors in R n Math 254 at OSU) that the length of a vector x = x 1 x 2... x n ) T R n, denoted x,
More informationTHE KADISON-SINGER PROBLEM AND THE UNCERTAINTY PRINCIPLE
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9939(XX)0000-0 THE KADISON-SINGER PROBLEM AND THE UNCERTAINTY PRINCIPLE PETER G. CASAZZA AND ERIC WEBER Abstract.
More informationSobolev Duals of Random Frames
Sobolev Duals of Random Frames C. Sinan Güntürk, Mark Lammers, 2 Alex Powell, 3 Rayan Saab, 4 Özgür Yılmaz 4 Courant Institute of Mathematical Sciences, New York University, NY, USA. 2 University of North
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 informationRadu Balan. Siemens Corporate Research. 755 College Road East. Princeton, NJ ABSTRACT
Multiplexing of Signals using Superframes Radu Balan Siemens Corporate Research 755 College Road East Princeton, NJ 854 ABSTRACT In this paper I intend to present the concept of superframes and its use,
More informationSPARSE SIGNAL RESTORATION. 1. Introduction
SPARSE SIGNAL RESTORATION IVAN W. SELESNICK 1. Introduction These notes describe an approach for the restoration of degraded signals using sparsity. This approach, which has become quite popular, is useful
More informationEFFICIENT ALGORITHMS FOR THE DISCRETE GABOR TRANSFORM WITH A LONG FIR WINDOW
EFFICIENT ALGORITHMS FOR THE DISCRETE GABOR TRANSFORM WITH A LONG FIR WINDOW PETER L. SØNDERGAARD Abstract. The Discrete Gabor Transform (DGT) is the most commonly used signal transform for signal analysis
More informationarxiv:math/ v1 [math.fa] 5 Aug 2005
arxiv:math/0508104v1 [math.fa] 5 Aug 2005 G-frames and G-Riesz Bases Wenchang Sun Department of Mathematics and LPMC, Nankai University, Tianjin 300071, China Email: sunwch@nankai.edu.cn June 28, 2005
More informationEfficient and Accurate Rectangular Window Subspace Tracking
Efficient and Accurate Rectangular Window Subspace Tracking Timothy M. Toolan and Donald W. Tufts Dept. of Electrical Engineering, University of Rhode Island, Kingston, RI 88 USA toolan@ele.uri.edu, tufts@ele.uri.edu
More information5742 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 55, NO. 12, DECEMBER /$ IEEE
5742 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 55, NO. 12, DECEMBER 2009 Uncertainty Relations for Shift-Invariant Analog Signals Yonina C. Eldar, Senior Member, IEEE Abstract The past several years
More information2 PETER G. CASAZZA, MANUEL T. LEON In this paper we generalize these results to the case where I is replaced by any positive selfadjoint invertible op
FRAMES WITH A GIVEN FRAME OPERATOR PETER G. CASAZZA, MANUEL T. LEON Abstract. Let S be a positive self-adjoint invertible operator on an N-dimensional Hilbert space H N and let M N. We give necessary and
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2018 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationAN INVERSE EIGENVALUE PROBLEM AND AN ASSOCIATED APPROXIMATION PROBLEM FOR GENERALIZED K-CENTROHERMITIAN MATRICES
AN INVERSE EIGENVALUE PROBLEM AND AN ASSOCIATED APPROXIMATION PROBLEM FOR GENERALIZED K-CENTROHERMITIAN MATRICES ZHONGYUN LIU AND HEIKE FAßBENDER Abstract: A partially described inverse eigenvalue problem
More informationCoding the Matrix Index - Version 0
0 vector, [definition]; (2.4.1): 68 2D geometry, transformations in, [lab]; (4.15.0): 196-200 A T (matrix A transpose); (4.5.4): 157 absolute value, complex number; (1.4.1): 43 abstract/abstracting, over
More informationTensor-Tensor Product Toolbox
Tensor-Tensor Product Toolbox 1 version 10 Canyi Lu canyilu@gmailcom Carnegie Mellon University https://githubcom/canyilu/tproduct June, 018 1 INTRODUCTION Tensors are higher-order extensions of matrices
More informationA Guided Tour from Linear Algebra to the Foundations of Gabor Analysis
A Guided Tour from Linear Algebra to the Foundations of Gabor Analysis Hans G. Feichtinger, Franz Luef, and Tobias Werther NuHAG, Faculty of Mathematics, University of Vienna Nordbergstrasse 15, A-1090
More information