clear all; close all;hold on;
om=0:pi/128:2*pi;
z=exp(i*om);
%z=exp(i*om).^(4);
Qz=3/8 * z .^(2) - 9/4 * z .^(1) +19/4 - 9/4 * z .^(-1) + 3/8 * z .^(-2);
Qn=[3/8 , -9/4 , 19/4 , -9/4 ,  3/8];
%Qz=z.^10 .* (3/8 * z .^(2) - 9/4 * z .^(1) +19/4 - 9/4 * z .^(-1) + 3/8 * z .^(-2))
plot(om,abs(Qz),'cyan');
plot(om,real(Qz),'blue');
plot(om,imag(Qz),'red');
plot(om,unwrap(angle(Qz)),'black');

PnDb3=conv (conv([ 1 3 3 1],[1 3 3 1]),Qn)
% nasledne staci spravit P(z), resp. P(OMEGA) z Pn
% ...
PnDb2=conv (conv([ 1 2 1],[1 2 1]), [-1 4 -1])
% nasledne staci spravit P(z), resp. P(OMEGA) z Pn
% ...