voicefilter.m

来自「对声音中的杂音进行过滤」· M 代码 · 共 20 行

function out=voicefilter(v)
Fs=22050;
[b,a] = ellip(8,0.05,80,[150 4000]*2/Fs);
subplot(3,1,1);
[H,w] = freqz(b,a,512);
plot(w*Fs/(2*pi),abs(H));
xlabel('Frequency (Hz)');
ylabel('Mag. of frequency response');
subplot(3,1,2);
plot(v);
xlabel('Time');
ylabel('Mag.');
subplot(3,1,3);

sf = filter(b,a,v);
plot(sf);
xlabel('Time');
ylabel('Mag.');
out=sf;