Tb=1;
f=0.01:0.01:3;
x=pi*f*Tb;
%sp_hdb3_per=(32*cos(2*x)-40)./(465*(1025-64*cos(10*x))).*(75.5-1797*cos(2*x)+1073*cos(4*x)+320*cos(6*x)-288*cos(8*x)+(565125.5-580957*cos(2*x)+15215*cos(4*x))./(16*cos(6*x)+24*cos(4*x)+44*cos(2*x)-85));
sinc=sin(x/2)./(x/2);
sp_5B6B_per=1./192*(205-44*cos(2*x)-23*cos(4*x)-12*cos(6*x)-2*cos(8*x)-4*cos(10*x)-(195+300*cos(2*x)+195*cos(4*x)+180*cos(6*x)+150*cos(8*x)+60*cos(10*x))./(17-8*cos(12*x)));
sinc_patrat=sinc.*sinc;
sp_5B6B=sp_5B6B_per.*sinc_patrat;
figure(1);plot(f,sp_5B6B,'g');grid;title('DSP pt 5B6B');xlabel('frecventa');ylabel('DSP');