miércoles, 21 de noviembre de 2012

EJERCICIO 27
Creación de SUBSISTEMAS en SIMULINK

Un sistema puede ser creado de dos formas:

1. Encerrando los bloques involucrados y creando un nuevo bloque SUBSISTEMA

2. Intercalando un bloque SUBSISTEMA vacio y posteriormente colocar los bloques correspondientes al mismo.



Se ha convertido en un subsistema el conjunto de bloques encerrados en el cuadro negro negro.




Aqui se ve el contenido del subsystema creado, es importante colocarle los elementos de entrada y salida (in, out)


Respuesta obtenida en el osciloscopio



miércoles, 7 de noviembre de 2012

EJERCICIO 26


Movimiento Armónico Forzado

Aqui se desarrolla la simulación numérica de la ecuación diferencial y''+2y'+10y=F(t)
F(t)=5*Cos(4*t + 45)





martes, 6 de noviembre de 2012

EJERCICIO 25

Sistema armónico amortiguado


Una masa de 8 Ib de peso estira 2 ft un resorte. Si una fuerza de amortiguamiento numéricamente
igual a 2 veces la velocidad instantánea actúa sobre el contrapeso, deduzca la ecuación del movimiento si la masa se suelta de la posición de equilibrio con una velocidad hacia arriba de 3 ft/s.


 x''+8x'+l6x=0
La ecuación auxiliar  es m^2 + 8m + 16 = (m + 4)2 = 0, de forma que rn1 = m2 = -4.
Luego el sistema es críticamente amortiguado y
x(t) = cl*e^-4*t+ c2*t*e^-4*t
Al aplicar las condiciones iniciales x(O) = 0 y x’(O) = -3 vemos, a su vez, que cl = 0 y c2 = -3.
Así, la ecuación del movimiento es
x(t) = -3*t*e^-4*t




Sistema armónico subamortiguado

 x''+2x'+10x=0

x(t) = (-2cos(3t)-(2/3)sen(3t) )*e^-t 





lunes, 5 de noviembre de 2012

EJERCICIO 24

Simulación de un circuito resistivo con SimPowerSystem de Simulink.

R= 6 Ohm
Fuente de poder de amplitud 24V, fase 75, frecuencia 60 Herz

En este caso la corriente esta en fase con el voltaje.

SimPowerSystem de Simulink permite simular el circuito eléctrico tal como lo tendríamos en el laboratorio. Esto permite obtener los valores de corriente y voltaje sin realizar la simulación numérica y colocar todos las ecuaciones propias del circuito.

En este caso en particular se tiene un circuito sencillo con una resistencia, utilizada de RCL Branch. Un medidor de corriente y voltaje de measurents y una fuente de poder de power.

La colocación de los elementos que realizan la medición es como se hace en el laboratorio, el medidor de corriente se coloca en serie con el circuito y el medidor de voltaje en paralelo.



viernes, 2 de noviembre de 2012

EJERCICIO 23

CIRCUITO R-C con una señal de entrada cuadrada






EJERCICIO 22

CIRCUITO R-C con fuente alterna senoidal con amplitud de 0.5 y frase 45º



EJERCICIO 21

Simulación de un circuito R-C. Procedimiento de carga del condensador.

rc.gif (1854 bytes)
Qo=45.10^-6 C
R=58Kohm
C=1.5x10^-6 F
V=30 volt

La ecuación diferencial a resolver es de la foma Q'=(V/R)- (1/RC)Q

Siendo la carga máxima en un tiempo muy grande (matemáticamente en t igual al infinito) a C*V


La intensidad de corriente es maxima en i=V/R









EJERCICIO 20

Simulación de un circuito R-C. Procedimiento de descarga.

Qo=45.10^-6 C
R=58Kohm
C=1.5x10^-6 F

Se sabe que Q'= - Q(1/ RC))



La salida SCOPE representa la carga en función en función del tiempo. Validamos nuestro modelo considerando t=0 obteniendo que la carga resulta ser la condición Qo. De igual forma para la corriente en función del tiempo, en la salida SCOPE1
EJERCICIO 19

Usar la interface GUI de Matlab en un programa.

Inicialmente se llama desde el workspace la aplicación guide. Posteriormente se tiene una ventana para diseñar nuestro cuadro de dialogo.

Las propiedades de cada elemento se modifican en la opción Property Inspector


Guide genera dos archivos, uno para la figura y otro .m que guardara el código fuente de la aplicación


jueves, 1 de noviembre de 2012

EJERCICIO 18

Simulación de una Ecuación Diferencial Ordinaria de 1er orden.

y'+ky=0    
y' = -ky

Condición inicial y(0) = 1


EJERCICIO 14

Simulación de la Sucesion de Fibonacci. (Mal llamada serie de Fibonacci)


En matemática, la sucesión de Fibonacci es la siguiente sucesión infinita de números naturales:
0,1,1,2,3,5,8,13,21,34,55,89,144,233,377
\ldots \,

La sucesión inicia con 0, y a partir de ahí cada elemento, es la suma de los dos anteriores(0,1,1,2,3,5,8...)



2f


En la era digital la Sucesión de Fibonacci se escribe así:
1 0 1 1 0 1 0 1 1 0 1 1 0 1 0 1 1 0 1 …






EJERCICIO 17

Modele la rata de flujo del agua en un tanque.






Tenemos que la ecuación que simula nuestro ejercicio es (A/K)Q' + Q=Qi, Tx'+x = ku(t) para resolverla se utiliza el operador Transformada de Laplace 's' y a partir de allí diseñamos nuestra simulación.


Qi es la tasa de entrada de agua medida en m^3/seg
Q es la tasa de salida medida en m^3/seg
A/K en donde A es el área del tanque y K es una constante proporcionalidad
k representa el comportamiento de estado del fluido, para un estado estacionario k=1







Podemos variar los parametros de la funcion Step y los parametros del operador de Laplace. Se observa que con un numerador mayor la salida de agua es mayor y viceversa.
EJERCICIO 16

Considere una señal de entrada senoidal y otra cosenoidal. Obtenga la respuesta en ventanas de salida de WORSKPACE






Recordar utilizar los bloques de To go workspace, clock, nombrar las variables y llamarlas en workspace en forma vectorial


EJERCICIO 15

Considere una señal de entrada senoidal y otra cosenoidal. Obtenga la respuesta en donde ventanas diferentes






El block Sine Wave 1, tiene una amplitud de 3 y una fase de Pi/2. Es decir, que tenemos una función Coseno con amplitud 3



miércoles, 31 de octubre de 2012

EJERCICIO 13

BLOQUE SIMOUT en SIMULINK Matlab en simulación de Temperatura en Farenheit







Este bloque permite visualizar en el WORKSPACE los valores guardados en la variable YOUT, TOUT
 u otra variable definida.
Los cambios en los parámetros de simulación permitirán obtener mayor o menor cantidad de valores.

En este ejemplo simularemos la ecuación de la Temperatura en grados Farenheit con valores de entrada en grados Celcius.






En Workspace llamamos a la variable yout y obtenemos los valores de la temperatura entre 0 y 15 grados celcius. Se puede aumentar el rango aumentando el tiempo de simulación.







Utilizando en el workspace de Matlab plot(tout,yout) se obtiene una gráfica que se puede guardar en un archivo jpg, modificarle los colores del fondo, colocarle etiquetas, etc.



martes, 30 de octubre de 2012

EJERCICIO 12

Modelo de ecuaciones diferenciales sencillas

d(y)/dt = 1

Solución de la ecuación diferencial: y(t) = t +A







d^2(y)/dt^2= 1



Solución de la ecuación diferencial: y(t) = 0.5*t^2+A*t+B Y y'(t)= A*t




EJERCICIO 11

PÉNDULO SIMPLE

Modelar el comportamiento de un péndulo simple.

El péndulo simple se rige por una ecuación diferencial, que se obtuvo luego de realizar las consideraciones físicas. Se considero sen(theta)=theta para theta menor a 1

Con una solución de la forma:


Considerando g/l=5, condiciones iniciales theta(t=0) = Pi/2 y para la derivada d(theta(t=0)) = 2*Pi.
Theta es el angulo inicial y su derivada la velocidad angular, luego de tener la simulación se puede probar con otros parámetros de borde y verificar el comportamiento del sistema.















lunes, 29 de octubre de 2012

EJERCICIO 10

BLOQUE RAMP

Simula expresiones de la forma y(t)= m*t+b


Slope Pendiente de la rampa.
Start time Tiempo en el cual se inicia la rampa.
Initial output Valor inicial de la rampa.

Para este caso en particular la salida es una funcion y(t)= t, ya q m=1 y b=0

EJERCICIO 9

Simulacion utilizando bloque CONSTANT como un vector escalar 1-D




EJERCICIO 8
Calculos de temperatura

function t_f=temperatura_celsius(x);
t_f=(5/9)*x+32
% esta funcion determina la temperatura en farenheit dada en celsius



>> x=[10 20 30 40]



x =

    10    20    30    40



>> temperatura_celsius(x)

t_f =

   37.5556   43.1111   48.6667   54.2222


ans =

   37.5556   43.1111   48.6667   54.2222

EJERCICIO 7

Archivos tipo m.

Con expresiones functions.


function p=f(x);
p=x.^2
% esta funcion determina el cuadrado de x

Esta función creada por usuario se ha guardado en un archivo tipo .m, llamado f.m.
Es importante tomar el nombre que genera Matlab por defecto con el fin de no crear conflictos si se tienen gran cantidad de funciones.

>> x=[1 2 3 4 5 6 7 8 9 10]

x =

     1     2     3     4     5     6     7     8     9    10

>> f(x)

p =

     1     4     9    16    25    36    49    64    81   100


ans =

     1     4     9    16    25    36    49    64    81   100


>> x=[1 2; 1 2]

x =

     1     2
     1     2

>> f(x)

p =

     1     4
     1     4


ans =

     1     4
     1     4

>> x=[1 2 3; 4 5 6; 7 8 9]

x =

     1     2     3
     4     5     6
     7     8     9

>> f(x)

p =

     1     4     9
    16    25    36
    49    64    81


ans =

     1     4     9
    16    25    36
    49    64    81


EJERCICIO 6

Realize un modelo para calcular la temperatura en Farenheit según entradas con temperaturas en grados celsius.




 Sabemos que la expresión general para calcular la temperatura viene dada por F=5/9C+32

 Utilizamos los bloques:

MUX: Permite ingresar valores de entrada, estos pueden ser escalares, matrices, funciones.
DEMUX: Permite obtener los valores de salida
DISPLAY: Muestra en pantalla los valores obtenidos







Se puede crear en Matlab un archivo .m creando una function.




EJERCICIO 5
Generar un modelo con una señal de onda cosenoidal.

Para general una señal de onda cosenoidal se utiliza el bloque para una señal senoidal y posteriormente de realizar un desfase de la onda en Pi/2

Podemos utilizar una sola hola de trabajo.



Se tiene como gráficas de salida: