hehe1

clc;
load 'dtmf53'
fd=941;
fg=1209;

fp=f_probkowanie;

deltaW=2*pi*(fg-fd)/fp;
wc=(((fd+fg)/2)/fp)*2*pi;
A=40;
M=511; %6.2*pi/deltaW
n=1:M;
d=(M-1)/2;

hd=sin(wc*(n-(d*0.99999)))./(wc*(n-(d*0.99999))+0.0000000000001);
okno=hamming(M);
h=hd.*okno';
h=h./sum(h);
fwp=zeros(1,M);
fwp(d)=1;
fgp=fwp-h;
widmo=20*log10(abs(fft(fgp,fp)));
n1=0:fp-1;
sig=filter(fgp,1,signal);
t=0:1/fp:(numel(signal)-1)/fp;

%przed filtracja figure(1); subplot(2,1,1); plot(n1,widmo); grid; subplot(2,1,2); plot(20*log10(abs(fft(fgp,fp))));grid;
figure(2); plot(t,sig); %postac czasowa
figure(3); subplot(2,1,1);plot(fp,20*log10(abs(fft(sig))));grid;subplot(2,1,2);plot(fp,20*log10(abs(fft(signal))));grid; %po filtracji