5.2

clear all
close all
w3dB = 2*pi*100;
N = [2 , 4, 6, 8];

    for k = 1:N(1)
    p1(k) = w3dB * exp(1i*((pi/2) + (1/2)*(pi/N(1))+(k-1)*(pi/N(1))));
    end

    for k = 1:N(2)
    p2(k) = w3dB * exp(1i*((pi/2) + (1/2)*(pi/N(2))+(k-1)*(pi/N(2))));
    end

    for k = 1:N(3)
    p3(k) = w3dB * exp(1i*((pi/2) + (1/2)*(pi/N(3))+(k-1)*(pi/N(3))));
    end

    for k = 1:N(4)
    p4(k) = w3dB * exp(1i*((pi/2) + (1/2)*(pi/N(4))+(k-1)*(pi/N(4))));
    end

    w = 0:1000;


    a = poly(p2);
    b = prod(-p2);

    H = polyval(b,w)./polyval(a,w);


    %step(H)
    impulse(H)

%     figure(1)
%     plot(real(p2),imag(p2),'or');grid;