From: Roland Kruse on
The problem can be efficiently solved in the spatiall frequency domain as demonstrated here:
http://matlabdb.mathematik.uni-stuttgart.de/download.jsp?MP_ID=31