clf;clear all
rzad_filtru=10;
f_gr=60;
fs=600;
N=1200;
t=(1:N)/fs;
f1=40; f2=100; %czestotliwosci sinusoid
A1=1; A2=1; %amplitudy sinusoid
s=A1*sin(2*pi*f1*t)+A2*sin(2*pi*f2*t)+randn(1,N); %sygnal testowy
[b,a]=butter(rzad_filtru,2*f_gr/fs) %wspolczynnik b i a ]
%
[H,wf]=freqz(b,a,N) %transmitancja filtru
s_filtr=filter(b,a,s); %sygnal po filtracji

%wykkresy----------------------------------
subplot(2,1,1)
fn=(wf/(pi))*fs; % czestotliwosc unormowana  (0-1)
plot(fn/2, abs(H), 'r') ; grid on
xlabel('Czestotliwosc[Hz');
ylabel('char.ampl.filtru|H(f)|');
HM=abs(H); %modul funkcji transmitancji (char.amplitudowa)
Hmax=max(HM); %wartosc max modulu funkcji transmitancji
subplot(2,1,2)
plot(fn/2,20*log10