关于fourier的疑问

bigkevin86

简单介绍哈算法，主要还是编程的问题，望高人解答。
w(x,y)是初始相位函数，Δw(x,y)是相位差函数，算法的目的是想通过相位差函数，恢复相位函数。
波前相位重建的fourier方法，Δw(x,y)=w(x+s/2,y)-w(x-s/2,y)
由fourier变化可得到：         F[w(x,y)](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的函数，不知道哪里有问题，望达人解答。