




tf=1
t=0.01:0.01:1
Dteta2=-6*pi./(tf.^3).*t.^2+6*pi/(tf.^2).*t

bin_cal=[ -9.0745160 -9.0997273 -9.1313981 -9.1691866 -9.2126621 -9.2613030 -9.3144929 -9.3715138 -9.4315397 -9.4936272 -9.5567070 -9.6195736 -9.6808761 -9.7391092 -9.7926066 -9.8395362 -9.8779000 -9.9055385 -9.9201426 -9.9192730 -9.9003896 -9.8608912 -9.7981675 -9.7096619 -9.5929456 -9.4458004 -9.2663071 -9.0529358 -8.8046310 -8.5208867 -8.2018041 -7.8481253 -7.4612385 -7.0431508 -6.5964290 -6.1241099 -5.6295869 -5.1164796 -4.5884974 -4.0493075 -3.5024160 -2.9510722 -2.3982000 -1.8463595 -1.2977377 -7.5416495e-001 -2.1715267e-001  3.1205649e-001  8.3243044e-001  1.3430892  1.8432506  2.3321846  2.8091775  3.2735085  3.7244358  4.1611944  4.5830009  4.9890654  5.3786067  5.7508710  6.1051508  6.4408035  6.7572682  7.0540806  7.3308851  7.5874439  7.8236434  8.0394975  8.2351483  8.4108647  8.5670380  8.7041762  8.8228972  8.9239193  9.0080520  9.0761853  9.1292788  9.1683504  9.1944646  9.2087216  9.2122459  9.2061758  9.1916531  9.1698137  9.1417783  9.1086445  9.0714795  9.0313127  8.9891301  8.9458688  8.9024121  8.8595852  8.8181511  8.7788071  8.7421814  8.7088294  8.6792300  8.6537820  8.6327998  8.6165080]


for i=1:100
   pot(i)=Dteta2(i).*abs(bin_cal(i))
end

subplot(2,1,1)
plot(t,bin_cal)
title('Binário')
grid on

subplot(2,1,2)
plot(t,abs(pot))
title('Potência')
grid on

