%zistovanie casovo-frekvencneho okna funkcii
x=-5:0.02:5
phi=(1-x.*x).*exp(-x.*x/2); % morletov wavelet
plot(x,phi);

syms t
%sym_phi=(1-t*t)*exp(-t*t/2)
sym_phi=1/(2+t*t)
t0=int(t*sym_phi*sym_phi,-inf,inf)
sigma_tmp=int((t-t0)^2 * sym_phi^2,-inf,inf)

energia_t=int(sym_phi^2,-inf,inf)

vysledok_sigma_t=sigma_tmp/energia_t