Help with Discrete Sine Transform

  • #1
21
0
Hi,
I am a neophyte in Discrete Fourier Transform and I am procticing with discrete Sine-transform.
Specifically I want to calculate $$ \mathcal{F}_s \lbrace x \cdot e^{-x^2} \rbrace = \int\limits_{0}^{\infty } x \cdot e^{-x^2} \cdot \sin\left(x\right)\, \mathrm{d}x
= \frac{1}{2} \pi^2 \omega e^{- \frac{1}{4} \pi^2 \omega^2} $$
I am using the FFTw library and with a 64 points I get a very poor result. Can anybody please help to improve it?
This is my C++ code, plot and numerical data are attached.
Thank you very much in advance
 
Last edited:

Answers and Replies

  • #2
scottdave
Science Advisor
Homework Helper
Insights Author
1,801
754
I don't see any attachments.
Your itex "command is not working as you intend to produce laTex. I think you should use double number-signs, before and after. I think this is what you want.
$$ \mathcal{F}_s \lbrace x \cdot e^{-x^2} \rbrace = \int\limits_{0}^{\infty } x \cdot e^{-x^2} \cdot \sin\left(x\right)\, \mathrm{d}x
= \frac{1}{2} \pi^2 \omega e^{- \frac{1}{4} \pi^2 \omega^2} $$

You can always click the preview button before posting a message. There is a latex help guide available. https://www.physicsforums.com/help/latexhelp/
 
  • #4
mathman
Science Advisor
7,890
460
Your original equation has a typo (sin(?)). It would be helpful to give more information about the graphs - which is which?
 
  • #5
21
0
Hi mathman,
the original equation is $$ f(x) = x \cdot e^{-x^2}$$ . Its analytical sine-transform is given by
$$ \mathcal{F}_s \lbrace f(x) \rbrace (\omega) = \int\limits_{0}^{\infty } x \cdot e^{-x^2} \cdot \sin\left(x\right)\, \mathrm{d}x
= \frac{1}{2} \pi^2 \omega e^{- \frac{1}{4} \pi^2 \omega^2}$$
I uploaded the plot again with a legend
 

Attachments

  • #6
mathman
Science Advisor
7,890
460
I presume the argument for the sin is ωx. I suspect that you need a finer mesh for the Fourier transform for small ω.
 
  • #7
scottdave
Science Advisor
Homework Helper
Insights Author
1,801
754
I presume the argument for the sin is ωx. I suspect that you need a finer mesh for the Fourier transform for small ω.
It has been awhile since I have done these, but I was kind of wondering how the ω just popped in when it is not in the integral.
 
  • #9
I like Serena
Homework Helper
6,577
176
Hi kaniello,

There seems to be a bit of a mixup with factors pi.
W|A (link) reports the sine transform to be:
$$\mathscr F_s\{xe^{-x^2}\}(\omega) = \frac 1{2\sqrt 2}\omega e^{-\frac 14 \omega^2}$$
Of course we can have a different normalization constant than the one given, but not a different argument to the exponential function.

Anyway, if I do a discrete fourier transform on the data, I don't get the dip that you're seeing between 2 and 7.
So it appears there is something wrong with that FFTW routine that you're using.
Unfortunately I'm not familiar with it.

I do notice that the dip might be explained if the cosine transform somehow crept in - that one does have a dip between 2 and 7.
And I guess it might also be caused by a filter (window function) that is applied before applying the FFT.
Have you tried using an FFT routine from a different library?
 
Last edited:
  • #10
21
0
Dear I like Serena,
thanks a lot for your reply.

Which library did you use to perform the sine transform of ## xe^{-x^2} ## ? It would be really interesting to compare the results.

I do not apply any filter to my Input data but I will try to understand from the FFTw Website if it is somehow executed automatically
 
  • #11
I like Serena
Homework Helper
6,577
176
You're welcome.
I didn't use any library - for a quick check I just used Excel's built in Fourier transform (one of the Data Analysis Tools).
 
  • #12
scottdave
Science Advisor
Homework Helper
Insights Author
1,801
754
Hello @kaniello . I just thought I would add some information for you, since you did state you are new to Discrete Fourier Transforms.
Are you aware of how sampling frequency affects the ability to regain the original signal? Every time varying signal is composed of various sinusoidal frequencies.
For the highest frequency in the time varying signal, you must sample at a rate at least twice of that frequency.

You may find this Wikipedia article helpful, to start. https://en.wikipedia.org/wiki/Nyquist_frequency

In real-world digital signal processing, we want to make sure that we are sampling at a high enough frequency. That is why a low pass filter is used on the signal before sampling - to filter out the higher frequencies that can lead to "aliasing".
 

Related Threads on Help with Discrete Sine Transform

  • Last Post
Replies
15
Views
2K
  • Last Post
Replies
2
Views
3K
Replies
5
Views
7K
Replies
1
Views
2K
Replies
1
Views
1K
Replies
8
Views
4K
Replies
10
Views
2K
Replies
2
Views
1K
Replies
2
Views
4K
Replies
8
Views
3K
Top