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

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);