%% ---------------------- FUNKCIJA: LOG_ANVELOPA_EN ---------------------- %%

function [ym_en]=log_anvelopa_en(M,ym);


% Funkcija koja sluzi za ljepsi prikaz izlaznih signala.
% Radi se estimacija energije (E(x^2)) svakog kanala, i konvolucija niskopropusnim
% filtrom (firgauss), koji ima impulsni odziv duzine nekoliko puta vece
% nego sto je najsira latica impulsnih odziva ulaznih filtara. Time
% dobivamo anvelopu signala. Na kraju se jos signal logaritmira.

b=firgauss(4,256);  %K=4,N=256, lenght(b)=K*(N-1)+1
b=b/sum(b);

for i=1:M+1,
    ym_en(:,i)=conv(b,ym(:,i).^2);
end;

% Nakon konvolucije signal koji dobijemo je duzine ulaznog signala + reda
% firgauss filtra, pa ga smanjimo na pocetku i nakraju za pola reda filtra.

ym_en=ym_en((length(b)-1)/2+1:length(ym_en)-(length(b)-1)/2,:);

% Logaritmiranje i mnozenje s 10 (jer je x^2).

ym_en=10*log10(ym_en);