I want to find 2D-Fourier transformation of the function given below

f = HeavisideTheta[y1]*HeavisideTheta[y2 - y1]

For the purpose, I use built-in function in two ways as below,

FourierTransform[f, {y2, y1}, {q2, q1}] 
FourierTransform[f, {y1, y2}, {q1, q2}] 

The outputs are different depending on the order of variables. Please, See below.

-(1/(2 π q1 q2 + 2 π q2^2)) + (I DiracDelta[q2])/(2 q1) + 
(I DiracDelta[q1 + q2])/(2 q2) + 1/2 π DiracDelta[q2] DiracDelta[q1 + q2]

-(1/(2 π q1 q2 + 2 π q2^2)) + (I DiracDelta[q2])/(2 q1) - (I DiracDelta[q1 + q2])/(2 q1)

I thought they should be same. I don't know why they are different.

Please help me to figure out which one is correct.

  • $\begingroup$ I think this is a bug. If you do a one dimensional transformation from the time domain to the frequency domain of a standard function with a HeavisideTheta you get a DiracDelta in addition to the correct result. Could an expert confirm that you cannot have a DiracDelta in the frequency domain where functions are analytic except at poles and branch points. A work around could be to use the Laplace transform and convert the answer. $\endgroup$ – Hugh Mar 5 '14 at 12:34
  • $\begingroup$ Try defining g[y1_, y2_] = HeavisideTheta[y1]*HeavisideTheta[y2 - y1] and see if you get what you expect. $\endgroup$ – b.gates.you.know.what Mar 5 '14 at 12:52
  • $\begingroup$ I guess it might be related to this and this. $\endgroup$ – xzczd Mar 6 '14 at 8:24

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.