首先我描述物理学,它是在轴对称空间中,一个声源放置在原始点,一个传感器放置在源下的轴上。给出源波形,我试着得到传感器的波形。所有物质参数都是已知的,例如声速,密度。
我写了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)))。因为我看到其他人使用它们,但我不明白为什么?
答案 0 :(得分:0)
如果你想要一个严格的实值结果波形(不是很复杂的虚构件),那么IFFT的输入必须是共轭对称的,例如: IFFT(dc_term中,M,0,fliplr(缀(M)))。