matlab音频信号噪音去除

Question

下面的带阻滤波器部分看不太懂，不知道是不是设置问题，无法正常运行
clear;
close all;
clc;
[x,FS,]=audioread('musicnoise.wav'); % 将 WAV 文件转换成变量
N=length(x); % 计算数据点数

% 原信号波形频谱分析
tx=(0:N-1)/FS; % 计算原信号数据点时刻
subplot(3,2,1);plot(tx,x); % 绘制原信号波形
title(' 原信号波形图 '); % 加标题
xf=fft(x); % 求原信号频谱
fx=(0:N/2)*FS/N; % 确定频谱图频率刻度
subplot(3,2,2);plot(fx,abs(xf(1:N/2+1))); % 绘制原信号频谱
title(' 原信号频谱图 '); % 加标题
disp('播放原声')
sound(x,FS)
pause(11)

fpd=3500;fsd=3550;fsu=4000;fpu=4050;Rp=1;AS=100;
fed=(fpd+fsd)/2;
fcu=(fpu+fsu)/2;
df=min((fsd-fpd),(fpu-fsu));
M=3396;L=(M+mod(M,2))/2;
wsd=fsd/FS;
wsu=fsu/FS;
F=[0:1/L:1];
A=[ones(1,fix(wsdM)),zeros(1,fix(wsuM)-fix(wsdM)),ones(1,L+1-fix(wsuM))];
B=fir2(M,F,A);
y_fil=fftfilt(B,x);
Y_fil=abs(fft(y_fil));Y_fil=Y_fil(1:length(Y_fil)/2);
figure;
subplot(3,2,1);
plot(tx,y_fil);
grid on;
xlabel('时间');
title('过滤后波形');
% subplot(3,2,2);
% plot(fx,Y_fil);
% grid on;
% xlabel('频率');
% title('过滤后频谱');
sound(Y_fil,FS);
麻烦改一下吧

Answer 1

你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答


本次提问扣除的有问必答次数，已经为您补发到账户，我们后续会持续优化，扩大我们的服务范围，为您带来更好地服务。