Inverse Fourier transform

  • Thread starter Telemachus
  • Start date
  • #1
833
30
Hi there. I have some trouble with this. I have to find the inverse fourier transform for: [tex]\frac{e^{i 6\omega}}{\omega}[/tex]
So I'm using a table, then:
[tex]F^{-1}\left ( \frac{e^{i 6\omega}}{\omega}\right )=F^{-1}\left ( e^{i 6\omega}\right ) * F^{-1}\left ( \frac{1}{\omega}\right )=2\pi\left[ \delta (t+6)*\frac{i}{2}sg(t)\right][/tex] where (*) represents the convolution.

Well, everything fine till there, but when I tried to corroborate my result with mathematica I get:
attachment.php?attachmentid=36378&stc=1&d=1307828102.gif


I don't know whats wrong.
 

Attachments

  • ft.GIF
    ft.GIF
    2 KB · Views: 635

Answers and Replies

  • #3
833
30
I don't know which integral you are referring. I'm using that [tex]F(sg(t))=\frac{2}{i \omega}[/tex] and [tex]F(\delta(t-t_0))=e^{-i \omega t_0}[/tex]

Then [tex]F^{-1}(e^{i6\omega})=\delta(t+6)[/tex] and [tex]F^{-1}\left (\frac{1}{\omega}\right )=F^{-1}\left ( \frac{i}{2}\frac{2}{i \omega}\right )=\frac{i}{2}sg(t)[/tex]
 
  • #4
hunt_mat
Homework Helper
1,745
26
The convolution is given by:
[tex]
f\star g(x)=\int_{-\infty}^{\infty}f(y)(g(x-y)dy
[/tex]
Then the Fourier transform is just the product of the Fourier transforms of f and g. So in this case you have [itex]\hat{f}(\omega )=e^{6\omega i}[/itex] and [itex]\hat{g}(\omega )=\omega^{-1}[/itex] and so the founction you're tring to find is:
[tex]
i\sqrt{\frac{\pi}{2}}\int_{\infty}^{\infty}\delta (y+6)\textrm{sign}(y-x)dy
[/tex]
 
  • #5
833
30
I follow your reasoning till here:
[tex]
f\star g(x)=\int_{-\infty}^{\infty}f(y)(g(x-y)dy
[/tex]

But I don't know how you get this from there:

[tex]
i\sqrt{\frac{\pi}{2}}\int_{\infty}^{\infty}\delta (y+6)\textrm{sign}(y-x)dy
[/tex]
Particularly this: [tex]i\sqrt{\frac{\pi}{2}}[/tex]
I don't know where you get that from.

I get from here: [tex]2\pi\left[ \delta (t+6)*\frac{i}{2}sg(t)\right]=i \pi \int_{\infty}^{\infty}\delta (t+6)\textrm{sign}(x-t)dy[/tex] so I make [tex]x_0=-6[/tex], then:
[tex]i \pi \int_{\infty}^{\infty}\delta (x+6)\textrm{sign}(t-x)dx=i \pi \textrm{Sign}(t+6)[/tex]
 
Last edited:
  • #6
764
71
Edit: incorrect
 
Last edited:
  • #7
833
30
I think it's fine. Look, I have the table right here.

[tex]F^{-1}(e^{-i\omega t_0})=\delta(t-t_0)[/tex]
Then I have [tex]t_0=-6[/tex]

Thats what I've used.
 
  • #8
764
71
Yes, I think I made a mistake. I don't know how Mathematica arrives at that result.
 
  • #9
833
30
Yes, I don't know neither. And actually the problem is I don't know where I have the mistake, I think that my solution is fine, maybe I'm not doing it well in mathematica, I'm thinking that there are different ways for defining a Fourier transform. You can use the frequencies or the angular frequency. Maybe I'm just not doing it well with mathematica, perhaps I have to define something else, thats a possibility.
 
  • #10
Ray Vickson
Science Advisor
Homework Helper
Dearly Missed
10,706
1,722
Hi there. I have some trouble with this. I have to find the inverse fourier transform for: [tex]\frac{e^{i 6\omega}}{\omega}[/tex]
So I'm using a table, then:
[tex]F^{-1}\left ( \frac{e^{i 6\omega}}{\omega}\right )=F^{-1}\left ( e^{i 6\omega}\right ) * F^{-1}\left ( \frac{1}{\omega}\right )=2\pi\left[ \delta (t+6)*\frac{i}{2}sg(t)\right][/tex] where (*) represents the convolution.

Well, everything fine till there, but when I tried to corroborate my result with mathematica I get:
attachment.php?attachmentid=36378&stc=1&d=1307828102.gif


I don't know whats wrong.

Here is what I get using Maple:


fw:=exp(I*6*w)/w;
exp(6 I w)
fw := ----------
w


> ft:=invfourier(fw,w,t);

ft := 1/2 I (2 Heaviside(6 + t) - 1)

> convert(%,piecewise);

{ -1/2 I t < -6
{
{ undefined I t = -6
{
{ 1/2 I -6 < t

RGV
 
  • #11
833
30
Thats even weirder.
 
  • #12
vela
Staff Emeritus
Science Advisor
Homework Helper
Education Advisor
15,095
1,675
If you're just wondering about the constants out front, you have to remember there are various conventions of where to put the [itex]2\pi[/itex]. Mathematica has a parameter called, appropriately enough, FourierParameters that you can set depending on which convention you want to use.
 
  • #13
833
30
Thank you vela :) do you have any idea about how mathematica gets to that? I think I should read more carefully the help content on this topic, but I don't have to much time for this, so maybe if you know the underlying formalism for this I'll appreciate if you let me know :P

Does it have any sense the solution with the heaviside function that Ray found using Maple? thats not a constant difference, but I'm aware of the similarities between the unit step function and the sign function, they differ only for the negative real values (I have no idea about what happens on the complex field).
 
Last edited:
  • #14
vela
Staff Emeritus
Science Advisor
Homework Helper
Education Advisor
15,095
1,675
By default, this is what Mathematica uses for the inverse Fourier transform:

[tex]f(t) = \frac{1}{\sqrt{2\pi}}\int_{-\infty}^\infty F(\omega)e^{-i\omega t}\,d\omega[/tex]
So Mathematica would give you, for example,

[tex]F^{-1}[e^{i6\omega}] = \frac{1}{\sqrt{2\pi}}\int_{-\infty}^\infty e^{i6\omega}e^{-i\omega t}\,d\omega = \sqrt{2\pi}\left[\frac{1}{2\pi}\int_{-\infty}^\infty e^{i(6-t)\omega}\,d\omega\right][/tex]
where the quantity in the square brackets is the Dirac delta function [itex]\delta(t-6)[/itex].

From looking at the results you found in the table, I think you're using the convention

[tex]f(t) = \int_{-\infty}^\infty F(\omega)e^{i\omega t}\,d\omega[/tex]
If you look in the Mathematica documentation on integral transforms, there's a list of what you want to set FourierParameters to to follow the convention you're using in your class.

As to the unit step and sign functions, they're related by a vertical shift and vertical scaling. You should be able to figure that out pretty easily.
 
  • #15
hunt_mat
Homework Helper
1,745
26
The definition of the integral transform I am using is:
[tex]
\hat{f}(\omega )=\int_{-\infty}^{\infty}f(x)e^{-i\omega x}dx
[/tex]
All I did was take the information from the wikipedia page and use that. The coefficient [itex]i\sqrt{\pi /2}[/itex] came from examining the fourier transform of the sign function.
 
  • #16
833
30
Thank you vela.
hunt_mat, I think wikipedia uses the formalism that vela just mentioned. You used the Fourier transform for the sign that you found in wikipedia or the one I've posted?
 
  • #17
hunt_mat
Homework Helper
1,745
26
having said that now, I am a little unsure of the sign I used.
 
  • #18
Ray Vickson
Science Advisor
Homework Helper
Dearly Missed
10,706
1,722
That is also the definition that Maple uses. Of course, that means that the inverse FT has the form (1/(2*pi))* \int \hat{f}(w) exp(i*w*x) dw.

RGV
 
  • #19
833
30
Thank you all.
 

Related Threads on Inverse Fourier transform

  • Last Post
Replies
2
Views
2K
  • Last Post
Replies
4
Views
2K
Replies
2
Views
20K
  • Last Post
Replies
9
Views
980
  • Last Post
Replies
2
Views
1K
  • Last Post
Replies
1
Views
992
  • Last Post
Replies
3
Views
3K
  • Last Post
Replies
3
Views
3K
Replies
1
Views
875
Replies
2
Views
1K
Top