filtr52

N1=10; N2=20; N3=40; N4=80; h0=1;
fg = 500;
fs= 8000;
v = 2*fg/fs;
m1=1:floor(N1/2);
m2=1:floor(N2/2);
m3=1:floor(N3/2);
m4=1:floor(N4/2);
h1=(sin(pi*v*m1))./(pi*v*m1);
h2=(sin(pi*v*m2))./(pi*v*m2);
h3=(sin(pi*v*m3))./(pi*v*m3);
h4=(sin(pi*v*m4))./(pi*v*m4);
h1=[fliplr(h1) h0 h1];
h2=[fliplr(h2) h0 h2];
h3=[fliplr(h3) h0 h3];
h4=[fliplr(h4) h0 h4];
b1 = h1/sum(h1);
b2 = h2/sum(h2);
b3 = h3/sum(h3);
b4 = h4/sum(h4);
[h10,w10] = freqz(b1,1,2^16,fs);
[h20,w20] = freqz(b2,1,2^16,fs);
[h40,w40] = freqz(b3,1,2^16,fs);
[h80,w80] = freqz(b4,1,2^16,fs);
figure();
plot(w10,20*log10(abs(h10)) );
hold on;
plot(w20,20*log10(abs(h20)));
plot(w40,20*log10(abs(h40)));
plot(w80,20*log10(abs(h80)));
xlabel('Częstotliwość[Hz]');
ylabel('Magnituda[dB]');
title('Charakterystyka amplitudowa')
legend('Rzad 10','Rzad20','Rzad 40','Rzad 80')
axis([0 4000 -60 5]);
hold off;

% b)
 N1 = 10;
N2 = 20;
N3 = 40;
N4 = 80;
h0=1;
m1=1:floor(N1/2);
m2=1:floor(N2/2);
m3=1:floor(N3/2);
m4=1:floor(N4/2);
h1=(sin(pi*v*m1))./(pi*v*m1);
h2=(sin(pi*v*m2))./(pi*v*m2);
h3=(sin(pi*v*m3))./(pi*v*m3);
h4=(sin(pi*v*m4))./(pi*v*m4);
h1=[fliplr(h1) h0 h1];
h2=[fliplr(h2) h0 h2];
h3=[fliplr(h3) h0 h3];
h4=[fliplr(h4) h0 h4];
b1 = h1/sum(h1);
b2 = h2/sum(h2);
b3 = h3/sum(h3);
b4 = h4/sum(h4);

     [h10,w10] = freqz(b1,1,2^16,fs);
     [h20,w20] = freqz(b2,1,2^16,fs);
     [h40,w40] = freqz(b3,1,2^16,fs);
     [h80,w80] = freqz(b4,1,2^16,fs);
     f1=angle(h10);
     f2=angle(h20);
     f3=angle(h40);
     f4=angle(h80);
     g1 = unwrap(f1);
     g2 = unwrap(f2);
     g3 = unwrap(f3);
     g4 = unwrap(f4);
      subplot(2,2,1);
      plot(w10,g1);
xlabel('Częstotliwość[Hz]');
ylabel('Faza [x pi rad]');
 title('Charakterystyka fazowa N = 10');
 subplot(2,2,2);
 plot(w10,g2);
xlabel('Częstotliwość[Hz]');
 Ylabel('Faza [x pi rad]');
title('Charakterystyka fazowa N = 20');
subplot(2,2,3);
plot(w10,g3);
xlabel('Częstotliwość[Hz]');
Ylabel('Faza [x pi rad]');
title('Charakterystyka fazowa N = 40');
subplot(2,2,4);
plot(w10,g4);
xlabel('Częstotliwość[Hz]');
Ylabel('Faza [x pi rad]');
title('Charakterystykafazowa N = 80');