f=100000;
T=1/f;
N=256;
t=[0:N-1]/N*T;
%t = linspace (0,1,1000);
t3=[t t+T t+2*T t+3*T];
z=square(2*pi*f*t,50);
z3=[z z z z];
plot(t3,z3);
axis([0 4*T -1.5 1.5])
XLABEL('3 Okresy');
YLABEL('Prostokatny');
pause;
X=abs(fft(z));
plot(X);
XLABEL('Nieznormalizowana czestotliwosc');
YLABEL('Harmoniczne');
pause;
funkcja=pltdbmag(X);
pause;
fp=0.016;
fz=0.032;

M = ceil(6.6/(2*(fz-fp))) + 8 		% hamming
b = firlpf((fp+fz)/2,M);
b = b .* hamming(M)';
Hejw = freqchar(b);
clgf;
HdB = pltdbmag(Hejw);
title('charakterystyka amplitudowa w dB - fp=0.2, fz=0.25, hamming');
pause;
Np = round(fp*max(size(HdB))) + 1;
Rpmax = -min(HdB(1:Np))
Nz = round(fz*max(size(HdB)));
Half = round(max(size(HdB))/2);
Azmin = -max(HdB(Nz:Half))
pause;
%mojeeee
%Az = 30;
%df = (fz - fp)
%M = 215%
%M = ceil((Az-7.95)/(14.36*df)) + 1
%pause;% kaiser
%beta = 0.5842*(Az-21)^0.4+0.07886*(Az-21);
%b = firlpf((fp+fz)/2,M);
%b = b .* kaiser(M,beta)';
%Hejw = freqchar(b);
%clgf;
%HdB = pltdbmag(Hejw);
%title('charakterystyka amplitudowa w dB - fp=0.2, fz=0.25, kaiser');
%pause;
%Np = round(fp*max(size(HdB))) + 1;
%Rpmax = -min(HdB(1:Np))
%Nz = round(fz*max(size(HdB)));
%Half = round(max(size(HdB))/2);
%Azmin = -max(HdB(Nz:Half))


filtered_signal=filter(b,1,z);
X2=abs(fft(filtered_signal));
plot(X2)