Master's Thesis - 5

In question 15 I want to make a picture how looks like for . I want to take at least two-dimensional vectors, as the inner product is too trivial otherwise. But then is a function on , which is rather difficult to plot. So I want to fix and plot , which is then a function on with values in and therefore easy to plot. As the inner product is , it is in the 2-dimensional case a simple plane-equation ; of course, the inner product is linear and goes through the origin! It cuts the zero-plane in an angle which is dependent on the orientation of the fixed : The zero-line evolves where the vectors are orthogonal to the fixed , i.e., .

This zero-line will now change its angle on the zero-plane while varying the fixed vector . This gives a picture how it behaves as a whole. In Octave, I set

[xx yy]=meshgrid(linspace(-2,2,50));

and then made a plot corresponding to by

mesh(xx+yy)

The resulting 3D-plot shows a boring inclined plane. To better see the location and orientation of the zero-line, I plot the absolute values:

mesh(abs(xx+yy))

To see the orientation even better, I chose

imagesc(abs(xx+yy))
imagesc(abs(xx+2*yy))
imagesc(abs(-xx+4*yy))
corresponding to the further vectors and . This is what I get:

The turning around of that plane now corresponds to the behavior of the mentioned inner product function on . (Beware that the absolute values used in the plots hide the asymmetry of the inner product function!) The values are always real and go into the exponent of the exponential function. then “rotates fastest” as runs along to a fixed , and is constant = 1 if it runs orthogonally; Beware that there’s a difference whether it is constant = 1 or if it reaches that value again and again. So, in the integral of the FT, the exponential function simply covers all frequencies via exhausting all values of , and the FT is then a function in , which takes out a single orientation of those frequencies; As already mentioned some time ago, the term “frequency” isn’t useful anymore in higher dimensions. For two-dimensional signals, i.e. images, a “frequency” can be interpreted as a pattern of lines: The narrower the pattern, the higher the frequency is. The difference to one-dimensional signals is, that those frequencies may now have an angle. It is not enough to just look at a single frequency of that line-pattern, but also to all their orientations! Only with using all orientations of patterns, an image can be analyzed according to the contained frequencies. The 2D-FT of an image therefore tells you what patterns occur in what orientation. To be more precise, you select the frequency by the length and the orientation by its angle , because the inner product as a function in grows faster with growing length of . For me, this is now a milestone in understanding the 2D-FT.

When I started on my way to my Masters degree seriously in January 2006, several questions began to arise in the course of time, and they continue to do so. While there were some very fundamental questions at the beginning (because my last contact to FA and basic analysis had been quite some time ago), they matured to some more advanced questions today. Here, I give a list of them, in their order of appearance. The answers, however, and as far as they have already been provided anyway, did not come up in the same order. Some questions are so general that they are to be understood as a challenge to look into the books and get familiar with it.

This list will continue to grow in separate entries. Also, for the sake of documentation, currently missing answers will only be given separately.

Here I give my personal reading approach to get familiar with the topics which the NuHAG is working on. This is just one way to get started, and I will mention some pros and cons I noticed. There is always some redundancy, but in my eyes that was of great worth.

A good way to get some first impressions is to read A Short Introduction to Gabor Analysis by Thomas Strohmer, which is the introductory chapter of [Feichtinger/Strohmer] and is also available online. You could also read the introductory chapters of [Blatter], which will make clear what the central ideas are in general. After having read that, you also want to understand it all.

I’m starting with the last chapter in Christensen and plan to have finished it on Nov. 5th. I still have some pages left in ch. 11 of Gröchenig as well as the whole ch. 12+13; I plan to have finished them by the mid of November. Then there’s the whole time of the day available for selected chapters of HGFei/Strohmer and—primarily—the GA4IP-stuff.

I also got myself into the MATLAB/Octave toolboxes and went through the diary (i.e. the MATLAB logfile) of HGFei and the examples of P. Søndergaard. (BTW, Søndergaard will be in Vienna at Nov. 8th for a talk at NuHAG.) Tomorrow I’ll visit the IT-staff at the institute to get a mathematics/MATLAB computing account. This term, I’ll visit HGFei’s lecture on Selected Topics in Modern Analysis, which takes place Monday and Thursday 10-11 a.m. Although this hour is “in the middle” of my working time at SOL4-IT, I’ll have to be there regularly, at any cost.

HGFei did some talks on Wiener Amalgam spaces and Banach Gelfand triples as well as introductions into their MATLAB toolboxes in the course of the so-called September ‘06 introduction, where some guests and participants of the EUCETIFA and Initiativkolleg were taking part. In addition, he sent me some first material about Gabor Analysis for image processing. Maybe this week or next I’ll visit Martin P. and ask him to get me a local computing account.

I’m almost done with the most relevant chapters in Christensen, but still not through the important chapters 11+12 of Gröchenig about modulation spaces. Meanwhile, I’ve bought G. Strang’s Linear Algebra in a German translation, as other literature I own about Linear Algebra—and even the course I took part at the beginning of my study—never mentioned topics as the pseudo-inverse or singular values and the singular value decomposition of matrices, the discrete version of frame theory.

Yesterday, I had a 1.5 hours talk with HGFei about further proceedings. I told him about my status and what I did the months before. He said that he could partly integrate me into the NuHAG from September on, where I’ll be attendant one afternoon and Friday morning every week. As working title for my master’s thesis we fixed “Gabor Analysis for Image Processing”: Some results have already been delivered in that area, and it should be my task to collect these together and give a survey. He will look for corresponding literature. During my stay at NuHAG, I could also be involved in the arising Computational Science project.

I’ll have to look if I should buy the Matlab student version. Also, the TF-toolbox of Peter Søndergaard seems to have been updated - it’s partly compatible to Octave.

I didn’t proceed much since the end of my holidays, because of managing photos and writing a travelogue, and my mum was celebrating an anniversary. I also had to hold a Linux course surprisingly, because one of the heads of my employer was saying goodbye.

The time plans will prolong for about one month, but now I get into it again.

On Thursday, I got another reply from HGFei to my “quality demand” mail. He told me that if I were interested to apply for the Initiativkolleg—a promoted doctoral program in association of NuHAG and DIANA—, one of the criteria I’d have to fulfill is to have the master thesis finished by January 2007. I replied that this still is my target date, but I don’t feel so sure about whether I’m qualified for pure mathematical research. He suggested to file the master thesis in May/June and asked me to pre-register as an applicant for the Initiativkolleg, as this doesn’t force me to anything, and all ways are kept open.

At the beginning of June, some colleagues of him will come to Vienna, and he will give some introductorial talks about the research at NuHAG, which I should listen to. As I’m on holiday until June 12th, I hope I won’t be missing much.