音频信号的傅立叶变换滤波

Question

我有一个数组，其中包含输入音频信号的傅里叶变换（对应于各种频率的幅度）。我希望选择特定的信号范围不使用内置功能。我已经执行了以下简单的操作：

[audio_in,audio_freq_sampl]=audioread('F:\Signals and Systems\Take Me Home Country Roads (John Denver Cover).wav');
Length_audio=length(audio_in);
df=audio_freq_sampl/Length_audio;
frequency_audio=-audio_freq_sampl/2:df:audio_freq_sampl/2-df;
figure
FFT_audio_in=fft(audio_in);

n = length(FFT_audio_in);
init = 30000;
fin = 40000;  
conji= mod((n-init+2),n) ;
conjf= mod((n-fin+2),n) ;
fs_1(1:n) = 0.0 ;
fs_1(init:fin) = FFT_audio_in(init:fin);
fs_1(conji:conjf) = FFT_audio_in(conji:conjf);

plot(frequency_audio,abs(fs_1));

正如我们在这里看到的，只有一个峰值。另一个应该在该范围内绘图的另一端可见。

这首歌可以在这里找到 - https://www.youtube.com/watch?v=WF046Z5tPJE

歌曲在阅读前必须转换为.wav文件。

上面的代码应该给我一个包含两个与频率范围相对应的小峰值的图 - (init , fin) 和 (conji , conjf)。但是，我得到的峰值仅对应于第一个范围。 这两个范围都在数组的大小范围内 - FFT_audio_in。

Answer 1

错误在于以下代码行：

n = length(FFT_audio_in);
init = 30000;
fin = 40000;  
conji= mod((n-init+2),n) ;
conjf= mod((n-fin+2),n) ;
fs_1(1:n) = 0.0 ;
fs_1(init:fin) = FFT_audio_in(init:fin);
fs_1(conji:conjf) = FFT_audio_in(conji:conjf);

原来在上面的代码行中，conji

n = length(FFT_audio_in);
init = n/4;
fin = n/4 + 10000;  
conji= mod((n-init+2),n) ;
conjf= mod((n-fin+2),n) ;
init_2 = min(conji , conjf);
fin_2 = max(conji , conjf);
fs_1(1:n) = 0.0 ;
fs_1(init:fin) = FFT_audio_in(init:fin);
fs_1(init_2:fin_2) = FFT_audio_in(init_2:fin_2);