Jason (jcreed) wrote,

When I awoke from troubled dreams this morning I found myself thinking about the discrete fourier transform. If you take a signal, and phase-shift each frequency component by some phase proportional to the square of its frequency, shouldn't that give you a kind of Laplace operator? The reasoning being that this is like saying that energy is proportional to the square of momentum, so the the time derivative of the result should look like the square of the momentum operator, hence a second spatial derivative. But if I do this with discrete time and space, it's not really a derivative, just some funny analogue of it.

In mathematica I can do
Delta[x_, y_] = If[x == y, 1, 0]; M = 100; K = 100;
H[s_] = s^2/(2 M K);
f[x_] = (1/(2 M)) Sum[
    Exp[2 I Pi (x s - H[s])/(2 M)], {s, -M, M - 1}];
rt = Table[{x, K*K*(Re[f[x]] - Delta[x, 0])}, {x, -20, 20}];
st = Table[{x, Re[f[x]] }, {x, -20, 20}];
it = Table[{x, K *Im[f[x]]}, {x, -20, 20}]; ListPlot[{st , it, rt }, 
 Joined -> True, Filling -> Axis, PlotRange -> (1 {-0.6, 1})]

and it yields this graph

which seems to be convergent as I increase M, the number of points in the world, and K, a scaling factor that when bigger makes every momentum count as less energetic. Maybe it is effectively acting like hbar (or its reciprocal) or something? I dunno.

Doing some numerical fiddling, it seems the spikes in the imaginary part (the red graph) just one step away from the origin on either side have the value of exactly 1/π --- not sure what causes that.
Tags: math, mathematica, physics

  • (no subject)

    Something that's bugged me for a long time is this: How many paths, starting at the origin, taking N steps either up, down, left or right, end up at…

  • (no subject)

    Still sad that SAC seems to end up being as complicated as it is. Surely there's some deeper duality between…

  • (no subject)

    I had already been meaning to dig into JaneSt's "Incremental" library, which bills itself as a practical implementation (in ocaml) of the ideas in…

  • Post a new comment


    Anonymous comments are disabled in this journal

    default userpic

    Your reply will be screened

    Your IP address will be recorded