clear ; clc ;
close all;
% Sygnal informacyjny
f0 = 1e5 ; % czestotliwosc podstawowa
T = 1/f0 ; % okres sygnalu informacyjnego
To = 3*T ; % czas obserwacji sygnalu
N = 2048 ; % liczba probek
t = [0:N-1]/N*To ; % os czasu (0.50 pkt)
m = square(t*2*pi*f0,50) ; % sygnal informacyjny (1.00 pkt)
plot(t,m)  % wykres sygnalu informacyjnego
pause;

% Filtr dolnoprzepustowy

dt = To/N ; % krok czasowy
fs = 0.25 ; % czestotliwosc probkowania (0.50 pkt)
fp = 0.10 ;% czestotliwosc przepuszczania (0.50 pkt)
fz = 0.3 ; % czestotliwosc zaporowa (0.50 pkt)
M = ceil(1.8/(2*(fz-fp))) + 3
fg = 0.5*(fp+fz) ;
b =  firlpf(fg, M); % licznik transmitancji filtru (1.00 pkt)
a = 1 ; % mianownik transmitancji filtru (1.00 pkt)
[Hejw, fnorm] = freqchar(b)  ; % transmitancja filtru (0.50 pkt)
plot(fnorm*fs,10*log10(abs(Hejw)),[fp,fz]*fs,[-3,-15],'o') ; % porownanie transmitancji z gabarytem
pause ;

s = filter(b,a,m);
plot(t,s)