Can anyone verify this code for spectral magnitude error?

%%Input image
origimage = imread('greekchurch','jpg');
// Smoothing via gaussian lowpass
window = fspecial('gaussian', 3, 0.5);
distimage= origimage*window
// Fastfourier transform of image
fftA = fft2(double(origimage));
fftB = fft2(double(distimage));
[r c] = size(origimage);
for i=1:r
for j=i:j
result=abs(fftA)-abs(fftB);
result = result*result/rc;
end
end

[Walter Roberson]

That is not valid MATLAB code. MATLAB does not support // comments. It might be Octave code.

[aliya]

No it is a Matlab code that's just a typo mistake

[Walter Roberson]

If it is MATLAB code then I can verify immediately that it is incorrect. MATLAB code never uses // in any context other than strings or already in a comment.

If it is only approximately the code, then we can only approximately evaluate it; in that case flip a coin, Yes or No, it is valid or it isn't: either would be approximately correct.

[aliya]

No, i just have tried this by following concept of SME. By coding perspective i am not sure about it.

[Image Analyst]

It seems like you're only reading comments, not any answers below. Any reason why not? Do I need to move my suggestion up here?

Answer 1

[Image Analyst]

What do you want to achieve? Are you sure you want to subtract then square, and not square then subtract?