Filtracja sygnalu

clc
clear
xdel(winsid())
[signal,Fs] = wavread("D:\Lukasz\Dokumenty\Studia\APS\Lab13\jagiello.wav");//bez szumu
[noise,Fs] = wavread("D:\Lukasz\Dokumenty\Studia\APS\Lab13\jagiello_15000.wav");//z szumem
//playsnd(signal,Fs); //odtworzenie dźwięku
//sound(signal,Fs);
signal=signal(1,:);//wyciągamy jeden kanał
noise=noise(1,:); //----//---------
//określamy zakres częstotliwości:
Fmin=20;//Hz
Fmax=20000;//Hz

scf()
analyze(signal, Fmin, Fmax, Fs, length(signal))
//można podejrzeć interesujące nas pasmo:
//Fm=4000;
//analyze(signal, Fmin, Fm, Fs, length(signal))
scf()
analyze(noise, Fmin, Fmax, Fs, length(noise))

h=wfir();//obliczamy współczynniki filtra
filt= convol(h,noise)//przefiltrowany sygnał (splot)

scf()
analyze(filt,Fmin, Fmax, Fs, length(filt))
playsnd(filt,Fs);