单音干扰频域matlab代码

当处理音频信号时，单音干扰是一种常见的问题。在频域中，可以使用傅里叶变换来分析和处理音频信号。下面是一个示例的MATLAB代码，用于处理单音干扰的频域滤波：

% 读取音频文件
[y, fs] = audioread('audio.wav');

% 进行傅里叶变换
Y = fft(y);

% 计算频率轴
f = (0:length(Y)-1)*fs/length(Y);

% 找到干扰频率的位置
interference_freq = 1000; % 设置干扰频率
[~, index] = min(abs(f - interference_freq));

% 将干扰频率处的幅度置零
Y(index) = 0;

% 进行逆傅里叶变换
y_filtered = ifft(Y);

% 播放滤波后的音频
sound(y_filtered, fs);

上述代码中，首先读取了音频文件，并进行了傅里叶变换得到频域表示。然后，通过找到干扰频率在频率轴上的位置，将该位置处的幅度置零。最后，通过逆傅里叶变换将滤波后的信号转换回时域，并播放滤波后的音频。