This page provides resources for the visualizations and the algorithms in the research manuscript:

# "Representation of the inverse of a frame multiplier"Peter Balazs and Diana T. Stoeva

## Abstract

Certain mathematical objects appear in a lot of scientific disciplines, like physics, signal processing and, naturally, mathematics. In a general setting they can be described as frame multipliers, consisting of analysis, multiplication by a fixed sequence (called the symbol), and synthesis. They are not only interesting mathematical objects, but also important for applications, for example for the realization of time-varying filters. In this paper we show a surprising result about the inverse of such operators, if existing, as well as new results about a core concept of frame theory, dual frames. We show that for semi-normalized symbols, the inverse of any invertible frame multiplier can always be represented as a frame multiplier with dual frames and reciprocal symbol. Furthermore, one of those dual frames is uniquely defined and the other one can be arbitrarily chosen. We investigate sufficient conditions for the special case, when both dual frames can be chosen to be the canonical duals. In connection to the above, we show that the set of dual frames determines a frame uniquely. Furthermore, for a given frame, the union of all coefficients of its dual frames is dense in l2. We investigate invertible Gabor multipliers; we show that the inverse of every invertible lattice-invariant operator (in particular, every invertible Gabor frame multiplier with a constant symbol (1)) can be represented as a Gabor frame multiplier with a constant symbol (1). Finally we give a numerical example for the invertibility of multipliers in the Gabor case.

## Figures:

### Figure 1 in the manuscript (an illustrative example to visualize actions of a multiplier on a musical signal)

In Figure 1 we show a visualization of a multiplier Mm,Φ,Ψ in the time-frequency plane. We consider a music signal f and the action of a multiplier Mm,Φ,Ψ on f. For f we use a 2 seconds long excerpt of the "Jump" from Van Halen (click here to listen the signal). For a time-frequency representation of the musical signal f (TOP LEFT) we use a 'painless' Gabor frame Ψ (a 80 ms Hanning window with 12,5% overlap). By manual estimation, we determine the symbol m that should describe the time-frequency region of the singer's voice. This region is then multiplied by 0.01, the rest by 1 (TOP RIGHT) (see the symbol here). Finally, we show a time-frequency representation of the modified signal Mm,Φ,Ψf (BOTTOM). To listen the modified signal,

 (TOP LEFT) The time-frequency representation of the music signal f (TOP RIGHT) The symbol m, found by a (manual) estimation of the time-frequency region of the singer's voice. (BOTTOM) Time-Frequency representation of Mm,Φ,Ψf.

### Figure 2 in the manuscript (numerical visualization of inversion of multipliers)

Here we use the same signal f and the same multiplier Mm,Φ,Ψ as in Figure 1. Note that all the elements of the symbol m fulfill mn,k∊{1,10-2}. Since m is semi-normalized, the multiplier M is analytically invertible [1]. However, the operator is badly conditioned, the condition number is around 99. The signal f is approximately 2 seconds long, using a sampling rate of 44100. Thus, the signal is a 128148-dimensional vector.

Starting from

we compare two approaches numerically:

• a 'naive' inversion

<>(corresponding to the approach raised in Question [Q2] in the manuscript)

• the 'iterative' inversion

Clearly, the naive approach has strong artifacts. The error is especially big at the boundaries of the constant region of the symbols. The chosen atoms are well localized in time-frequency, so that within the interior of the constant regions, this inversion works well. This could be expected as we have shown in the manuscript that constant symbols allow this kind of inversion for equivalent frames.

The iterative inversion worked well with an error of 3%. This could, naturally, be decreased by investing more calculation time. But also in the chosen setting for the iterative inversion (100 iterations in iframemul [2]) no difference can be seen in the time-frequency representation, as well as no audible difference can be detected.

 (TOP LEFT) The time-frequency representation of the result of the 'naive' inversion . (TOP RIGHT) The time-frequency representation of the error of the 'naive' inversion, i.e. . (BOTTOM LEFT) The time-frequency representation of the iterative inversion . (BOTTOM RIGHT) The time-frequency representation of  the error of the iterative inversion .

## Algorithms:

The above visualizations are done using algorithms in the Matlab/Octave toolbox Linear Time-Frequency Analysis (LTFAT) [2] (version 1.4.0 and above). In order to run the script, provided below, first you need to install the the Matlab/Octave toolbox LTFAT, freely available at Sourceforge.

The Matlab script for producing Figures 1 and 2 is available for download here. To run the scrip, you need to have the following two files in the same folder:

Running the script, the output is the Figures 1 and 2.

## References:

• [1] D. T. Stoeva and P. Balazs, "Invertibility of multipliers", Applied and Computational Harmonic Analysis, vol. 33 (2012), 292-299.
• [2] P. Søndergaard, B. Torrésani, and P. Balazs,"The linear time frequency analysis toolbox", International Journal of Wavelets, Multiresolution and Information Processing, vol. 10, no. 4, pp. 1250032, 2012.

## Upcoming Events

#### International Day of Mathematics 2020

14th March 2020

10:00 am – 6:00 pm

Prechtl-Saal, TU Wien
Karlsplatz 13, 1040 Wien

#### ARI-Guest-Talk: Imitation of novel conspecific and human speech sounds in the killer whale (Orcinus orca)

23. March 2020

Imitation of novel conspecific and human speech sounds in the killer whale (Orcinus orca) - José Francisco Zamorano Abramson

14.00 o'clock,

Seminar Room, Wohllebengasse 12-14 / Ground Floor

#### ARI-Guest-Talk: Causal Inference in Neuroimaging

31. March 2020

Causal Inference in Neuroimaging - Moritz Grosse-Wentrup

14.00 o'clock,

Seminar Room, Wohllebengasse 12-14 / Ground Floor

#### ARI-Guest-Talk: Inference and efficient coding in natural auditory scenes

15. April 2020

Inference and efficient coding in natural auditory scenes - Wiktor Młynarski

11.00 o'clock,

Seminar Room, Wohllebengasse 12-14 / Ground Floor

#### 23. Internationaler Tag gegen Lärm 2020

Aktionstag des ÖAW-Institut für Schallforschung

am 29. April 2020, 9:30 - 17:30 Uhr

www.oeaw.ac.at/tgl2020

#### Applied Harmonic Analysis and Friends

Strobl20

June 1st - 5th 2020

Strobl, AUSTRIA

#### 23rd International Conference on Digital Audio EffectsDAFx2020

DAFx20

Sep 8th - 12th, 2020

Vienna, AUSTRIA

https://dafx2020.mdw.ac.at

## News

• #### Deutsch in Österreich - Teil 2

28. February 2020

Das an der Universität Salzburg und dem ISF beheimatete Teilprojekt 02 des Sonderforschungsbereichs Deutsch in Österreich läuft jetzt schon 4 Jahre und wurde erst unlängst um 4 weitere Jahre...

• #### Interview mit Peter Balazs

31. January 2020

Im folgenden Interview von Peter Balazs gibt es Interessantes zur aktuellen Schallforschung und dem diesjährigen Jahr des Schall zu erfahren. Das Auftaktevent des Internationalen Jahr des Schalls...

• #### Mein Englisch klingt besser als Deins!

11. February 2020

Heute spricht Eva Reinisch vom Institut für Schallforschung der ÖAW über die kürzlich veröffentlichte Studie, die sich u.a. mit der Selbstüberschätzung von Englisch Lernenden beschäftigt hat. In der...

04. February 2020

The Austrian Academy of Sciences (ÖAW), Austria’s central non-university research and science institution is offering a position as a Academy Scientist (f*m) (40 hours per week) at the Acoustics...

• #### Öffentlicher Vortrag: Mathematik und Akustik

09. January 2020

Am Donnerstag 9. Jänner 2020 hält Peter Balazs seinen öffentlichen Vortrag über Mathematik und Akustik an der TU Wien im Rahmen des Forum Mathematik....

• #### Youtube Video: Mathematics and Acoustics

15. January 2020

Dr. Peter Balazs lecture on mathematics and acoustics held at the Technical University of Vienna on 9th January 2020.

• #### SynesthesiaSynästhesie

27. November 2019

Last Thursday it was time. The opening of the exhibition of Hermann Nitsch's orgies Myterien Theater Synesthesia took place there. The spectrogram by Peter Balazs and Anton Noll with the institute's...