cps lab6 zad2

clear variables; close all;

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

figure; spectrogram(y, 4096, 4096-512, [0:5:2000], Fs);

r = [0 9 5 8 7];

butt = load('butter.mat');
z = butt.z; p = butt.p; k = butt.k;
[zd,pd,kd] = bilinear(z,p,k,Fs);
[bd, ad] = zp2tf(zd, pd, kd);
yf = filter(bd, ad, y);

figure; spectrogram(yf, 4096, 4096-512, [0:5:2000], Fs);

t = (1:1:length(y))./Fs;

figure; plot(t,y, t, yf);