Retomo este tema para practicar un poco, resulta que encontré esto por casualidad y me puse a recordar los sistemas de control:
http://real2electronics.blogspot.com.es/2011/07/maqueta-de-control-pid-con-arduino.htmlhe echo ese mismo montaje con un envase de aceitunas, y con seis resistencias de 10homios, dos grupos de tres en serie puestas en paralelo, total 15homios. Y para medir la temperatura NTC de 10K al positivo, con una resistencia de 3k3 a masa, midiendo en la de fija de 3k3.
Con eso obtengo los siguientes datos:
tiempo=[0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210];
%son voltios de un divisor, ntc 10k al positivo, 3k3 a masa.
%medido en la de 3k3 a masa.
temperatura=[4.77 4.69 4.61 4.51 4.40 4.29 4.14 4.04 3.94 3.86 3.78 3.70 3.64 3.60 3.55 3.51 3.49 3.46 3.40 3.36 3.34 3.34];
plot(tiempo,temperatura)
%en eje x tiempo. Eje y voltios o temperatura.
Me voy a curvefitting y obtengo lo siguiente:
syms t
a = 3.406;
b = -0.1962;
c = 0.3603;
d = 0.5446;
yt = a*exp(b*t) + c*exp(d*t)
YS = laplace (yt)
%yt =
%(1703*exp(-(981*t)/5000))/500 + (3603*exp((2723*t)/5000))/10000
%YS =
%1703/(500*(s + 981/5000)) + 3603/(10000*(s - 2723/5000))
introduzco la función de tranasferencia multiplicado por 's':
g1=tf([1703 0],[500 500* 981/5000])
g2=tf([3603 0],[10000 -(10000*( 2723/5000))])
g=g1+g2
%rltool(g)
step(g)
y resulta un sistema inestable, tiene un polo en positivo.
Para obtener la función de la planta he puesto tensión en las resistencias, he esperado que se calienten, y luego he puesto en marcha el ventilador y cada 10 segundos he tomado un valor de voltaje en la resistencia de 3k3, es correcto ¿no? ¿alguna idea?
Saludos.