2015-05-04 21 views
-1

由于优化,我改变了矩阵在频域中的大小。我不确定如何补偿比例因子。目前我这样做,但我不确定它是正确的:fft2/ifft2的比例因子改变矩阵之间的大小

B1 = fft2(A1)/numel(A1); 
B2 = B1(idx1, idx2); %Cut out some parts of B1 
A2 = fft2(B2); %I want the sign change 

这是否正确缩放?如果是这样,请显示。

编辑:用numel改变乘法除法

回答

2

如果你的应用fftifft顺序,你可能不担心缩放,他们会互相抵消的意思ifft(fft(A))正好等于A。但是,在你的情况下,当尺寸发生变化,只是做

B1 = fft2(A1); 
B2 = B1(idx1, idx2); 
A2 = numel(B2)/numel(B1)*ifft2(B2); 

注意,在代码中,你已经如上图所示,你没有应用傅立叶逆变换第二次,这可能会导致错误。

+0

你能解释缩放吗?我很困惑fft2和ifft2的缩放比例。 我使用fft2的原因是因为它对应于我正在解决的标志应该翻转的实际问题(成像方程)。 – Forss