老式唱机播放音乐时，常会听到“噼啪”类噪声。使用MATLAB为带有这类噪声的音频文件，采用FIR数字滤波器滤除噪声。要求回放滤波前后的语音信号，并从时域和频域比较处理前后的差异。并绘制出滤波前后的波形图和频谱图，导出滤波后的音频。

好的，我会尽力回答您的问题。以下是详细的步骤：

1. 读取音频文件：使用MATLAB中的audioread函数读取带有噪声的音频文件。

[y, Fs] = audioread('noisy_audio.wav');

其中，y是音频信号，Fs是采样率。

2. 设计FIR数字滤波器：使用MATLAB中的fir1函数设计FIR数字滤波器。可以根据需要选择不同的滤波器阶数和截止频率。

order = 100;
fc = 2000/(Fs/2);
b = fir1(order, fc, 'low');

其中，order是滤波器阶数，fc是截止频率，'low'表示低通滤波器。

3. 滤波处理：使用MATLAB中的filter函数对音频信号进行滤波处理。

filtered_y = filter(b,1,y);

其中，b是滤波器系数。

4. 回放滤波前后的语音信号：使用MATLAB中的sound函数回放滤波前后的音频信号。

sound(y, Fs); % 播放滤波前的音频
sound(filtered_y, Fs); % 播放滤波后的音频

5. 比较处理前后的差异：可以通过比较滤波前后的音频信号的时域波形和频域谱图来分析处理前后的差异。

% 绘制滤波前后的波形图
t = (0:length(y)-1)/Fs;
subplot(2,1,1);
plot(t, y);
title('Original audio');
xlabel('Time (s)');
ylabel('Amplitude');
subplot(2,1,2);
plot(t, filtered_y);
title('Filtered audio');
xlabel('Time (s)');
ylabel('Amplitude');

% 绘制滤波前后的频谱图
nfft = 2^nextpow2(length(y));
Y = fft(y,nfft)/length(y);
f = Fs/2*linspace(0,1,nfft/2+1);
F_Y = 2*abs(Y(1:nfft/2+1));
F_Y = 20*log10(F_Y);
F_Y(F_Y<-60) = -60;
F_Y = max(F_Y,-60);
F_Y_filtered = 2*abs(fft(filtered_y,nfft)/length(filtered_y));
F_Y_filtered = 20*log10(F_Y_filtered);
F_Y_filtered(F_Y_filtered<-60) = -60;
F_Y_filtered = max(F_Y_filtered,-60);
figure;
subplot(2,1,1);
plot(f, F_Y);
title('Original audio spectrum');
xlabel('Frequency (Hz)');
ylabel('Magnitude (dB)');
subplot(2,1,2);
plot(f, F_Y_filtered);
title('Filtered audio spectrum');
xlabel('Frequency (Hz)');
ylabel('Magnitude (dB)');

6. 导出滤波后的音频：使用MATLAB中的audiowrite函数将滤波后的音频保存为WAV格式文件。

audiowrite('filtered_audio.wav', filtered_y, Fs);

这样，您就可以得到滤波前后的音频信号，并可以通过比较时域波形和频域谱图来分析处理前后的差异。同时，您也可以导出滤波后的音频文件。