匀速直线运动模糊

changccaiyang

小弟模拟匀速直线运动模糊始终无法得到正确结果，还望高人指点，在此谢过。
    一维退化函数如下：H(u,v)=Tsinc(pi*u*a)*exp(-j*pi*u*a),式中T为曝光时间，当前图像在x方向以速度x0=at/T做匀速直线运动.若y分量也以速度y0=bt/T做匀速直线运动,退化函数变为：H(u,v)=Tsinc(pi*(u*a+v*b))*exp(-j*pi*(u*a+v*b)).式中sinc(x)=sin(x)/x.
        具体过程先对图像f(x,y)傅里叶变换为F(U,V)，将退化函数H(U,V)离散化后与F(U,V)相乘得G(U,V)，然后对G(U,V)进行傅里叶反变换得到退化图像。但是在实现过程中一维退化函数能得到正确的结果；二维始终得不到意想的结果。不知道是不是二维退化函数离散化的方法有问题?
以下是二维离散化过程,图像大小为m*n
H=zeros(m,n);
for u=0:m-1
    for v=0:n-1
         H(u+1,v+1)=T*sinc(u*a+v*b)*exp(-i*pi*(u*a+v*b));      
%          H1(u+1,v+1)=T*sinc(v*a)*exp(-i*pi*v*a);
    end
end