针对频率 - 波数场的Matlab逆快速傅立叶变换,我需要进行共轭和翻转吗?

时间:2014-03-26 14:04:20

标签: matlab fft waveform

首先我描述物理学,它是在轴对称空间中,一个声源放置在原始点,一个传感器放置在源下的轴上。给出源波形,我试着得到传感器的波形。所有物质参数都是已知的,例如声速,密度。

我写了Matlab脚本来计算它,通过求解我能得到的声音传播方程 一个函数,比如A(w,k),w是频率,k是波数,这就是所谓的频率 - 波数场。我的matlab代码是这样的, 离散的w和k,得到一个A数组。首先使用FFT到k,获取空间和频率信息 然后,FFT到w,得到空间和时间信息,即不同点的波形。

虚假代码

for i_w=...
w=...
 for i_k=...
 k=...
 M=A(w,k)
 end 
wave_space_freq=ifft(M)  
end  % here can specify the only point of the sensor
wave_space_freq=ifft(wave_space_freq)

我的问题是当我使用IFFT时需要进行共轭和翻转,例如ifft(M,0,fliplr(conj(M)))。因为我看到其他人使用它们,但我不明白为什么?

1 个答案:

答案 0 :(得分:0)

如果你想要一个严格的实值结果波形(不是很复杂的虚构件),那么IFFT的输入必须是共轭对称的,例如: IFFT(dc_term中,M,0,fliplr(缀(M)))。