quinta-feira, 23 de dezembro de 2010

PRAZER, EU SOU A PROGRAMAÇÃO E VIM PRA FERRAR A TUA VIDA

f = logspace(1,5,10);
w = 2*pi*f;
R = 10e3;
C = 10e-9;
w0 = 1/(R*C);
H = (j*w)./(10e3+j*w);

f_medidas = [50 80 100 200 500 800 1000 2000 5000 8000 10000 20000];
magH_medidas = [0.036 0.051 0.068 0.14 0.335 0.048 0.58 0.82 1 1 1 1]; %valores medidos em bancada
faseH_medidas = [90 86. 85.5 78.46 66.31 62.6 58.37 38.4 19.72 15.48 10.71 7.29]; %valores medidos em bancada

figure(1)
traco = semilogx(f,abs(H),'k')
set(traco,'LineWidth',1.5)
xlabel('Frequência (Hz)','FontSize',16);
ylabel('Magnitude','FonteSize',16);
title('Magnitude da resposta em frequência','Fontsize',16)
grid
hold
marcador = plot(f_medidas, magH_medidas,'kx');
set(marcador,'MarkSize',10);
set(gca,'FontSize',14);
axis([10 1e5 0 1.2])
legend('Curva Teórica','Medidas Práticas')

Nenhum comentário: