关于fourier的疑问
简单介绍哈算法，主要还是编程的问题，望高人解答。w(x,y)是初始相位函数，Δw(x,y)是相位差函数，算法的目的是想通过相位差函数，恢复相位函数。
波前相位重建的fourier方法，Δw(x,y)=w(x+s/2,y)-w(x-s/2,y)
由fourier变化可得到： F(v)=T(v)*F[Δw(x,y)](v) v是空间域x对应的频域变量。
T(v)=1/2*i*sin(pi*s*v) T(v)是剪切转换函数
以下是程序：
clear all;
syms x y z v;
s=0.4;%剪切量
a=(x^2+y^2-2/3)*y+4/3;%原始相位
b=(((x+s/2)^2+y^2-2/3)*y+4/3)-(((x-s/2)^2+y^2-2/3)*y+4/3);%相位差函数
b
B=fourier(b,x,v);%对相位差函数进行fourier变换
C=1/(2*i*sin(pi*s*v));%剪切转换函数
D=B*C;
d=ifourier(D,v,x);%反变换
有算法知道，d应该和a是一个函数，如果通过取值来验证的话，相差也不应该太大，我仿真出来， d是一个只关于y的函数，不知道哪里有问题，望达人解答。
页:
[1]