FILTRY

N=55; %dlugosc odp. impulsowej = rzad filtru(liczba nieparzysta)
fp = 4000;
fn = fp/2;% czestotliwosc nyquista
fs = 60;% skladowe sygnalu
A = ones(size(fs)); % amplitudy skladowych

NP = 400; % granice pasma
t = linspace(0,N/fp,NP);% os czasu[s]
f = linspace(0,2,NP); % os czestotliwosci znormalizowana wzgledem fn
% generowanie sygnalu testowego
f1=100;f2=320;
s1=sin(2*pi*f1*t);
s2=sin(2*pi*f2*t);
syg=s1+s2;
syg=syg+randn(1,length(syg));%sygnał zaszumiony

a=1;
%górno przepustowo
%dolno przepustowy typ_filtru = 'low'
%pasmowo przepustowy Wa = .4; Wb = .5; Wn3 = [Wa Wb]; % granice pasma

Wn = 0.4;%czestotliwosc odciecia
typ_filtru = 'high';
Wn= 0.45; %czestotliwosc odciecia
b = fir1(N-1,Wn,typ_filtru); % projekt filtru
y = filter(b,a,syg);
%y2=freqz(syg,b,a,fs);

subplot(3,2,1)
plot(t,syg)
title 'Syganał Wejściowy'

subplot(3,2,2)
plot(f,abs(fft(syg)))
title 'Modul widma sygnalu wejsciowego'

subplot(3,2,3)
plot(f,abs(fft(y)))
title 'Modul widma sygnalu wyjsciowego'

subplot(3,2,4)
plot(t,y)
title 'Sygnal wyjsciowy'

%subplot(3,2,5)
%plot(f,20*log10(abs(fft(b,NP))))
%title 'Skala Logrytmiczna w [dB]'