在 impulse response 的使用中遇到了小小的问题,具体如图,第一行是原始信号,第二行是响应的信号,在通过将原始信号和相应信号分别转换到频域求得 transfer function 后,ifft 得到 impulse response。第三行是原始信号和 ir 进行卷积得到的结果,可以看到在信号最开始的阶段,有一个意料外的响应。:
这一段的代码如下:
% src: 原始信号
% resp: 响应的信号
Ssrc = fft(src);
Sresp = fft(rest);
Stf = Ssrc./Sresp;
ir = ifft(Stf);
z = conv(src,ir);
先说短的结论,把最后一条命令 conv
换为 cconv
即可。
和 Daniel 探讨后,问题出在 得到的 ir 上。
对于这个ir,是不能直接进行卷积的,而应该先做一个 circular shift ( ir = circsh