任选一副数字图像（图像格式自己定），用MATLAB编程对其进行图像变换（处理）。急急急

%傅里叶变换
clear
clc
A=imread('b.jpg');
imshow(A);
title 原图
figure;
[M,N]=size(A);
A=double(A);
%傅里叶正变换
B=zeros(M,N);
C=zeros(M,N);
for m=1:M
for n=1:N
A(m,n)=A(m,n)*(-1)^(m+n);
end
end
for v=1:N
for y=1:N
B(:,v)=B(:,v)+A(:,y)*exp(-i*2*pi*v*y/N);
end
end
B=B/N;
for u=1:M
for x=1:M
C(u,:)=C(u,:)+B(x,:)*exp(-i*2*pi*u*x/M);
end
end
C=C/M;
imshow(log(abs(C)+1));
title 频域图像