I want to implement an approach for calculating Phase spectrum of a set of coordinates. In the paper, it is described as 3 equations. However, I have no idea about the function variables (x,y) in the equation. I want to implement it via Matlab. Perhaps, it is easier for someone who knows pure signal processing well.
In (1), the 2D signal series (the coordinates) is represented as a 2D dirac delta function. Varialbe, theta_i, is the direction of each coordinate. The fourier transform of (1) is in the form of (2). Then (3) is the phase spectrum. So, maybe it can be implemented via a few lines in Matlab. But I just don't know it. Could anyone here give any help? Thanks.
This looks to me is essentially
theta is the function exp(j\theta_i) defined at each (xi,yi) position
Of course, it depends on how xi, yi, u, and v are defined.