[求助]使用傅里叶变换重建波前数据

[chenzheyi]

f是一个方形实数矩阵，大小为N*N，实际为一波面数据f=f（x，y）。
现要求实现以下公式：
F[∂2f/∂x2 +∂2f/∂y2 ]=-(u2+v2)F(u,v)，此即为频率域的拉普拉斯算子。
其中F(u,v)是f(x,y)是二维傅里叶变换。
等式左边是对f在x方向上与y方向上的二次偏导之和作傅里叶变换。
希望高手能实现，即给定f的数据矩阵，利用此公式求得F(u,v),并对此作傅里叶逆变换求得重建的f’，与原先的f作比较。