Hi,(adsbygoogle = window.adsbygoogle || []).push({});

I want to do the following calculation:

[tex]

c_{m,n} = \int_{-\infty}^{\infty} t^m \phi(t-n)\,dt

[/tex]

I know three things:

- I know the values of [tex]c_{m,n}[/tex] for m={0,1} and the first 4 for m={2,3}
- I know that my [tex]\phi[/tex] is symmetric, i.e. [tex]\phi(t) = \phi(-t)[/tex]
- The fourier transform of [tex]\phi(t)[/tex], i.e. [tex]\Phi(\omega)[/tex]

This is how I tried solving the problem:

[tex]

c_{m,n} = \int t^m \phi(t-n)\,dt =

\int t^m \phi(-(n-t))\,dt[/tex]

because of the symmetry

[tex] = \int t^m \phi(n-t)\,dt[/tex]

...is exactly the definition of a convolution...

[tex](t^m * \phi)(n) = \mathcal{F}^{-1}\left\{ \mathcal{F}\left\{t^m\right\} \Phi(\omega)\right\} = \mathcal{F}^{-1}\left\{ j^m \sqrt{2\pi} \delta^{(n)}(\omega) \Phi(\omega) \right\} =

j^m \sqrt{2\pi} \mathcal{F}^{-1}\left\{\delta^{(n)}(\omega) \Phi(\omega) \right\}

[/tex]

Now, writing the inverse fourier transform gives:

[tex]

j^m \sqrt{2\pi} \frac{1}{\sqrt{2\pi}} \int_{-\pi}^{\pi} \delta^{(n)}(\omega) \Phi(\omega) e^{j\omega n}\,d\omega = j^m \int_{-\infty}^{\infty} \delta^{(n)}(\omega) \underbrace{\Phi(\omega) e^{j\omega n}}_{f(\omega)}\,d\omega

[/tex]

Now it is well known that

[tex]\int \delta^{(n)}(x) f(x)\,dx = (-1)^n f^{(n)}(0)[/tex]

so that

[tex]

j^m \int_{-\infty}^{\infty} \delta^{(n)}(\omega) f(\omega)\,d\omega = j^m (-1)^m \left. \frac{\partial^m}{\partial \omega^m} f(\omega) \right|_{\omega = 0}

[/tex]

So the whole procedure should reduce to calculate the derivative of [tex]f(\omega)[/tex] and set the result to zero. I fill in my known [tex]\Phi(\omega)[/tex]:

[tex]

f(\omega) = \mathrm{sinc}^2 (\omega/2) \sqrt{ 1 - \frac{2}{3} \sin^2(w/2)} \cdot e^{j\omega n}

[/tex]

Problem 1: Even the 0th derivative is only computable via a limit. Calculating the limit with [tex]\omega \rightarrow 0[/tex] gives the result:

[tex]\sqrt{\frac{2}{2+\cos(2)}}[/tex]

But this is not true. The result should be 1 for [tex]m=0[/tex].

The result for [tex]m=1[/tex] should be linear: {0,1,2,3,...}. But as you can verifiy the result with my approach is constant.

Also the result for m={2,3} should give a a quadratic and cubic function. But with my approach it stays constant; Even worse, the solution is not always real!

.... So there is anywhere a mistake. Can anybody help me where?

Thank you,

divB

**Physics Forums | Science Articles, Homework Help, Discussion**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Integrating derivate of Dirac: Where is my fault?

**Physics Forums | Science Articles, Homework Help, Discussion**