求MATLAB 代码 注解

a=zeros(300,300); % 建一个300X300的零矩阵
a(50:250,100:150)=255; % 50-250行，100-150列，赋值 255
subplot(2,2,1) %画图1
imshow(real(a)); % 显示a的实部（其实就是a)
fftI=fft2(a); % 傅立叶变换 a
b=angle(fftI); % 角度值 取出来
subplot(2,2,2); % 画图2
imshow(real(b)); % 显示b的实部
D=ifft2(b); % 对b傅立叶逆变换
subplot(2,2,3);
imshow(real(D)) % 再画一次