function [H, W] = spektrum(h, N)
%   SPEKTRUM pocita FFT a zobrazi dvoustranne spektrum
%   [H, W] = spektrum(h, N)
%        h : casova posloupnost delky L
%        N : delka FFT [-pi,pi); (N >= L)
%        H : komplexni hodnoty spektra
%        w : frekvence

N          = fix(N);
L          = length(h);  
h          = h(:);                    
W          = (2*pi/N) * [ 0:(N-1) ]';
stred      = ceil(N/2) + 1;
W(stred:N) = W(stred:N) - 2*pi;  
W          = fftshift(W);
H          = fftshift( fft( h, N ) );