clear all;
om=linspace(0,2*pi,1000);
H=(1-2*exp(-i*om))./((1+3*i*exp(-i*om)).*(1-3*i*exp(-i*om)));
figure;
plot(abs(H));