how to plot the frequency response of the following transfer function.?

i have written the code as below

a0=4964-1102+273-39;
%a0=0.3038094552593955-0.068156473531124598+0.016118874664139887-0.0036378272766718939;
%a1=-0.068156473531124598+0.016118874664139887-0.0036378272766718939;
%a2=0.016118874664139887-0.0036378272766718939;
%a3=-0.0036378272766718939;
a1=-1102+273-39;
a2=273-39;
a3=-39;
for w=1:45:180;
s(w)=(1+exp(-1j*w));
x(w)=(a0*exp(-1j*6*w)+a1*(exp(-1j*4*w)-exp(-1j*8*w))+a2*(exp(-1j*2*w)-exp(-1j*10*w))+a3*(1-exp(-1j*12*w)));
y(w)=(a0*exp(-1j*6*w)-a1*(exp(-1j*4*w)-exp(-1j*8*w))-a2*(exp(-1j*2*w)-exp(-1j*10*w))-a3*(1-exp(-1j*12*w)))*exp(-1j*w);
*z(w)=s(w)*(x(w)+y(w));*
end

[Ryan G]

I don't see a transfer function. It looks like the for loop is attempting to achieve frequency domain response manually. If you have a transfer function you can use bode, margin, pzmap, rlocus and nyquist to observe the frequency response.

[mahesh sasanur]

H(z)=(1+z^-1){[(a0*z^-6+a1*((z^-4)-(z^-8))+a2*((z^-2)-(z^-10))+a3*(1-(z^-12))]+ (z^-1)*[(a0*z^-6-a1*((z^-4)-(z^-8))-a2*((z^-2)-(z^-10))-a3*(1-(z^-12))]}

[mahesh sasanur]

the transfer functin is given below H(z)=(1+z^-1){[(a0*z^-6+a1*((z^-4)-(z^-8))+a2*((z^-2)-(z^-10))+a3*(1-(z^-12))]+ (z^-1)*[(a0*z^-6-a1*((z^-4)-(z^-8))-a2*((z^-2)-(z^-10))-a3*(1-(z^-12))]}

how to odtian frequency response. a0=h1+h3+h5+h7 a1=h3+h5+h7 a2=h5+h7 a3=h7 h1, h3,h5, h7 are filter tap co-efficients