ios - 通过反卷积从录制的声音中消除已知音频以估计背景声级

Question

我有 2 个信号，一个包含在扬声器上播放的音频数据。第二个包含同时记录扬声器的麦克风数据。

到目前为止我做了什么:通过相关性在时域中对齐信号。对两个信号的重叠部分应用 fft 并将一个除以另一个以实现反卷积。
我做错了什么，因为生成的音频数据没有用。

这是我的代码:

       //put both signals in split complex vectors
       vDSP_ctoz((DSPComplex *)file, 2, &fftFileData, 1, nOver2);
       vDSP_ctoz((DSPComplex *)mic, 2, &fftMicData, 1, nOver2);


        //fft of both file and mic data
        vDSP_fft_zrip(fftSetup, &fftFileData, 1, log2n, FFT_FORWARD);
        vDSP_fft_zrip(fftSetup, &fftMicData, 1, log2n, FFT_FORWARD);


        //divide file data by mic data for deconvolution???
        vDSP_zvdiv(&fftFileData, 1, &fftMicData, 1, &fftMicData, 1, nOver2);  


        //inverse fft of mic-fft-data
        vDSP_fft_zrip(fftSetup, &fftMicData, 1, log2n, FFT_INVERSE);

        //scale back signal
        vDSP_vsmul(fftMicData.realp, 1, &scale, fftMicData.realp, 1, nOver2);
        vDSP_vsmul(fftMicData.imagp, 1, &scale, fftMicData.imagp, 1, nOver2);

        //copy back to float array
        vDSP_ztoc(&fftMicData, 1, (COMPLEX *) result, 2, nOver2);

编辑做一点澄清:感谢@Sammio2 我现在知道了，反卷积很好地描述了我的问题:

f*g=h

h 是我记录的信号，包括

f，我希望恢复我的信号

g，我知道另外记录了我的播放信号，但很可能被扬声器修改了->麦克风往返

现在我需要任何方法来恢复 f，它是除 g 之外的所有记录的声音。

重要提示:最后我不需要 f 的明确信号，只需要有关其响度或存在级别的信息。基本上是记录的往返信号 g 之外的噪声水平。

我应该如何继续收集我想要的噪音水平信息？

我希望这有助于理解我的问题。到目前为止谢谢!

Answer 1

vDSP_zvsub 的长度参数是要处理的复杂元素的数量，而不是元素数量的对数。您应该传递 nOver2 而不是 log2n。

这仅仅解决了编程方面的问题。其他答案解决信号处理问题。特别地，FFT 是线性的:给定信号 X 和 Y 以及常数 a 和 b，FFT(a•X+b•Y) = a•FFT(X)+b•FFT(Y)。逆 FFT 也是线性的。因此，两个信号的 FFT 之差的逆 FFT 不会给您带来与直接减去两个信号不同的结果，除了常见的浮点舍入误差。