function [msk_vystup] = MSK_demod(msk_signal,msk_signal_sum,fb,w,N)

M = 100;                          
krok = 1/(fb*M); 

% rozdelenie na I a Q zlozky

t = krok*(-M:M*N-1);
In_phase = msk_signal_sum.*cos(w*t);
Quadrature = msk_signal_sum.*sin(w*t);

% filtracia I(t) a Q(t) zloziek

[n,Wn] = buttord(0.1, 0.2, 1, 10); 
[b,a] = butter(n,Wn);             

In = filter(b,a,In_phase); 
Qu = filter(b,a,Quadrature);

% scitanie I(t) a Q(t) a detekcia

vy = In + Qu;

tmp2(1:N) = (M:M:M*N);
for i=1:N
    if  (vy(tmp2(i)) > 0)
        vystup(i) = 1;
    else
        vystup(i) = 0;
    end
end
msk_vystup = vystup;