语音滤波去电流噪音.zip

% 去噪 clear all close all clc dbstop if error warning off disp(['Started at ',datestr(now,'yyyy-mm-dd HH:MM:SS')]) %% 读取 [y,Fs]=audioread('WhyIHitYou.mp3'); y1=y(:,1); y2=y(:,2); %% 切比雪夫带通 % Wp = [60 200]/500; % Ws = [50 250]/500; % Rp = 3; % Rs = 40; % [n,Wp]=cheb1ord(Wp,Ws,Rp,Rs); %% 切比雪夫低通 Wp = 60/500; Ws = 300/500; Rp = 3; Rs = 60; [n,Wp]=cheb1ord(Wp,Ws,Rp,Rs); [b,a]=cheby1(n,Rp,Wp); %% 滤波 z1=filter(b,a,y1); z2=filter(b,a,y2); z=[z1,z2]; %% 存储 audiowrite('new.wav',z,Fs) disp(['Finished at ',datestr(now,'yyyy-mm-dd HH:MM:SS')])