如何用matlab求解快速正弦离散函数

在u是一个已知函数，M=N且已知的情况下，如何用MATLAB编写出来呢，请各位帮帮忙！！
maltab中的命令只有 idst ,但是idst只有一个求和符号 ！
请各位赐教！

可以使用双重循环实现，下面是例子：

M = 100;
N = 100;
u = ones(M,N);
k = 1;
l = 1;
ukl = 0;
for ii = 1:M-1
    for jj = 1:N-1
        ukl = ukl + (2/M)*(2/N)*u(ii,jj)*sin(ii*k*pi/M)*sin(jj*l*pi/N);
    end
end
ukl