Entrada # 3

Programa # 1 Estabilidad

Para esta ocasión tenemos que demostrar la estabilidad de nuestras funciones de transferencia, antes de comenzar les muestro nuestra función:

A continuación tenemos el comportamiento de los polos y ceros en forma gráfica de la función de transferencia:

Como vemos en la imagen vemos en la imagen, la función tiende a ser estable hasta un cierto punto en el cuál su comportamiento se vuelve inestable.

Después tenemos el código utilizado como vemos en los comentarios, los valores de la Resistencia y el Calor son Constantes, estos son inicializados a 1 para que no apliquen un gran cambio en el resultado:

	#La función es: (c*m(To)/s)/ (s*Q(s)*I*R)
	function estabilidad
	s = tf('s'); %Variable s
	c = 1; %Calor especifico del cautín, como es constante va 1
	m = 3; %Masa del cautin, como no tengo el peso va uno x
	To = 30; %Temperatura inicial del cautin, uno de prueba
	I = 100;%32;%12;%2; %Intensidad de la corriente, va uno de prueba
	R = 1; %Resistencia del material de la punta del cautin, como es constante va 1
	Q = 4;%Potencial de calor
	sys = transferencia(s, c, m, To, I, R, Q)
	step(sys, To); %Usando solamente step
	tiempo = 0:.001:30;
	ruido = sin(tiempo);
	lsim(sys, ruido, tiempo);
	rlocus(sys);
	endfunction
	%Como para sacar la función de tranferencia use varios pasos, todo esto lo paso a una función
	function paso3 = transferencia(s, c, m, To, I, R, Q)
	paso3 = 1;
	paso1 = tf(c*m*To);
	paso1 = tf(paso1/s);
	paso2 = tf(s*Q*s*I*R);
	paso3 = tf(paso1/paso2);
	return
	endfunction

view raw estabilidad.m hosted with ❤ by GitHub

El resultado del código lo mostramos a continuación, el resultado no muestra el comportamiento con distintas Intensidades de entrada y con sus respectiva gráfica cuando se le aplica un ruido (este ruido fue realizado mediante la función seno y un tiempo de 0 y 20 con .001 entre cada uno):

Con 2

                                                                        Entradas

                                                                       Con Ruido

Con 12

                                                                                Entradas

                                                                               Con Ruido

Con 32

                                                                               Entradas

                                                                             Con Ruido

Con 100

                                                                                Entradas

                                                                               Con Ruido

Como vemos en todas las imágenes anteriores, las gráficas nunca se detienen en un punto en especifico, por lo tanto se mantiene la propuesta de que nuestra función tiende a ser inestable.