mercoledì, ottobre 21, 2009

FM09010 (I)

So, let us start we the new lecture: Fourier methods and their applications to neuroscience.

Preliminary disclaimer: there are many many approaches to the topic. I will follow this one. In fact, Osgood's lecture is by far better than what I can hope to do in my life, but he has a slow pace which is not suitable for PhD students at the BCCN.

In the first lecture we explained the basic idea, which is simple: you have a complex periodic function (where complex stays for both "complicated" and "not real") and you want to represent it as weighted sum of simpler (but also not real) periodic functions. For definiteness, let us say that all functions have period 1.

We can choose to use trigonometric functions or complex exponentials. We will choose the complex exponentials for different reasons. To be precise, we want to express any periodic function of period 1 as the weighted (probably infinite) sum of the functions

e_n(t):=e^{2\pi nit}
where n ranges in the integers.
The difficult problem is: how to find the weights for the sum? It goes in the following way. You start from what you are looking for

f(t)=\sum_{n \in \mathbb Z} c_n e^{2\pi nit}

Here, the numbers c_n are the (unknown!) weights of the linear combination that you want to represent the function f. We now try to isolate a single coefficient, say the kth, to know whether is possible to get a formula for a single coefficient depending only on f and not on the other coefficients. So you get

f(t)- \sum_{n\neq k} c_n e^{2\pi nit} = c_ke^{2\pi kit}

You divide by the complex exponential obtaining

f(t)e^{-2\pi kit}- \sum_{(n-k)\neq k} c_n e^{2\pi nit} = c_k

For this you have to use all properties of the complex exponentials. Which is, by the way, only one. And in fact it is the same as for the real exponential: "sum is mapped into multiplication". This is one motivation for using complex exponentials instead of trigonometric functions.

You now have the problem that you have expressed one c_k in terms of the other ones, but you want to have c_k without its companions. Algebra gave us everything she could, so let us try with Analysis. He suggests to integrate between 0 and 1 (which is the period of all functions here) obtaining

\int_0^1 f(t)e^{-2\pi kit} dt- \int_0^1\sum_{(n-k)\neq k} c_n e^{2\pi nit}dt = c_k

The c_k, integrated, gives c_k! And now the magic. Compute the integral of the complex exponentials: it gives 0.
So, you are left with the famous expression

\int_0^1 f(t)e^{-2\pi kit} dt- \int_0^1\sum_{(n-k)\neq k} c_n e^{2\pi nit}dt = c_k=:\hat{f}(k)

where we have introduced the symbol \hat{f}(k). We will call this number by the suggestive name of the kth Fourier coefficient of f.

2 commenti:

Bluebeardburns ha detto...

Scusa, perchè non hai scritto a cosa corrispondono i simboli? Io non ho capito cosa è c: sarebbe il peso?

Lap(l)aciano ha detto...

Ciao caro,

pensavo fosse chiaro dopo la prima equazione. Comunque si, sono i pesi della combinazione lineare; ora aggiusto la spiegazione.

Grazie!
Stefano