2015-05-04 22 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