Fourier pseudo spectral method 这是要用到的方法

请问有人知道运用Fourier pseudospectral method 怎么写梯度算子的程序吗，Matlab
以下是我写的程序；但不对，希望能有盆友给我指出，非常感谢（哭！）
1：Nx=128;
2：Lx=1;
3：x = 0 + Lx/Nx*(0:Nx-1);
4：y = 0 + Lx/Nx*(0:Nx-1);
5：[xx,yy] = meshgrid(x,y);
6：v=[0:Nx/2-1,0,-Nx/2+1:-1];
7：v1=diag(v);
8：k_x=v1*(2*pi/Lx);
9：u=cos(xx).*cos(yy);
10：u1=fft2(u);
11：ku=k_x.*u1;
12：u_x=real(ifft2(ku)); %u对x求导