UNIVERSIDAD NACIONAL DE
SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PROCESOS
ESCUELA PROFESIONAL DE INGENIERÍA
METALÚRGICA
LABORATORIO Nº 10:
IMPLEMENTACION DE ECUACIONES DE SEGUNDO
ORDEN EN SIMULINK
PRESENTADO POR:
CCAHUANA LAYME RICHARD
Docente: ING. PEDRO NUÑEZ RAMIREZ
AREQUIPA - PERU
2021
Control de Procesos
PRÁCTICA DE LABORATORIO N° 10
IMPLEMENTACION DE ECUACIONES DE SEGUNDO ORDEN EN SIMULINK
1. OBJETIVO
La práctica tiene por objetivo resolver ecuaciones de segundo orden en el programa Simulink,
empleando mayor cantidad de recursos del programa.
2. EQUIPO Y MATERIALES.
1. Laptop o Computadora personal.
2. Servicio de Internet
3. Software MathLab – Simulink
4. Apuntes de la sesión de clase teórica
3. PROCEDIMIENTO.
IMPLEMENTAR LA SIGUIENTE ECUACION. Despejamos la Segunda derivada
previamente.
a. Ejercicio 1
Implementar un nuevo proyecto en Simulink e incorporar los elementos para la solución
En los parámetros de simulación del modelo indicamos 5 seg, el método de solución de la
ecuación diferencial SOLVER lo dejamos en automático.
2
Control de Procesos
3
Control de Procesos
b. Ejercicio 2. Implementar la ecuación siguiente ecuación en SIMULINK
4
Control de Procesos
CONCLUYENDO
Con el Simulik obtenemos la solución de ecuaciones de segundo orden que pueden
representar a diferentes modelos.
5
Control de Procesos
Elaborar 3 ejercicios similares de implementación de ecuaciones diferenciales asociadas a
sistemas de segundo orden con la aplicación del programa SIMULINK.
Proyecto 1
ECUACIONES DE SEGUNDO ORDEN
y ' ' =4 y' +2 y−2
6
Control de Procesos
Proyecto 2
y ' ' =8 y ' −16 y +1
7
Control de Procesos
Proyecto 3
y ' ' =−3 y ' −2 y+ 30 e2
8
Control de Procesos
Proyecto 4
1
x '' = ( x+5−2∗x ' )
3
4. RESULTADOS.
Se logro mostrar con resultados gráficos según los datos a diferentes soluciones de
ecuaciones diferenciales y comprobar su ejecución correcta dentro del programa Simulink.
9
Control de Procesos
5. CONCLUSIONES.
Se logro comprobar y cumplir los objetivos de la práctica en el programa Simulink para
control de procesos, en donde en el proyecto 4 se aprecia de que las curvas varían según
el orden de diferencial, es decir mientras se aumenta el integrador la curva tiende a
linealizarse (estabilidad).
6. BIBLIOGRAFIA.
[1] Simulación de Ecuación diferencial en Simulink - Matlab
[Link]
[2] Solución de ecuaciones diferenciales en Matlab y Simulink
[Link]
[3] Ecuaciones Diferenciales en Simulink de Matlab
[Link]
[4] Introducción a la simulación de sistemas dinámicos en “MatlabSimulink” a partir de
la ecuación diferencial que rige su comportamiento:
[Link]
%20a%20la%20simulaci%C3%B3n%20de%20sistemas%20din%C3%A1micos%20en
%20Matlab-Simulink%20a%20partir%20de%20la%20ecu....pdf?sequence=1&isAllowed=y
COMPLEMENTO
SOLUCIÓN ANALÍTICA DE ECUACIÓN DIFERENCIAL DE SEGUNDO ORDEN MEDIANTE MATLAB
EDO DE SEGUNDO ORDEN EN MATLAB
EMPLEO DE LIVE SCRIPT DE MATH LAB
Resultados de practica [Link] de Mathlab
CODIGO EN MATHLAB
%% RESPUESTA DE UN SISTEMA DE SEGUNDO ORDEN
%% a(d2/dt2).x+b.d/dt)x+c.x=u
%clear all
%close all
syms x(t) %x(t)=2*t^2
%diff(x,t,2) %syms
x(t) a=1; b=1; c=1;
10
Control de Procesos
u=1;
eqn=a*diff(x,t,2)+b*diff(x,t)+c*x==u
%condiciones inicales de un istema de 2do orden
% una condicon para x y otra para la derivada de x
%voy definir una variable de x que sea igual a la derivada de x
%luego declara el vector de condiciones iniciales de x en cero que seria
% %la derivada de x en su valor inicial cero
dx=diff(x,t)
ic=[dx(0);x(0)]==[0;0]
% utilizo el comando dsolve que recibe a la ecuación diferencial y a
las condiciones iniciales x(t)=dsolve(eqn,ic)
% luego variamos las condiciones iniciales [10;1]
%reemplazar en la la ecuación de condiciones iniciales ic
%ic=[dx(0);x(0)]==[10;1]
%sin condiciones iniciales x(t)=dsolve(eqn)
%x suponemos que x representa posicion y x punto representa velocidad Dx=diff(x)
%para graficar empleamos el comando fplot para la x y la derivada de x
%fplot(x)
%fplot(Dx)
%si las queremos colocar en una sola grafica empleamos hold
on fplot(x,[0,20]) hold on fplot(Dx,[0,20]) hold off
legend('x','Derivada de x')
D2x=diff(Dx,t) % SEGUNDA DERIVADA DE X
% pero si cambiamos un parámetro que ya no sea constante y sea en función del
tiempo ejm sen t
%a=1; b=1; c=1; u=1; en b remplazamos por sin(t) ya no se encuentra una solución
cerrada
%a=1; b=sin(t); c=1; u=1;
11
Control de Procesos
SALIDAS DEL LIVE EDITOR DE MATHLAB
12
Control de Procesos
Con la leyenda añadida
Si cambiamos un parámetro por una función del tiempo ya no encontramos una solución
X(t)=[empty sym]
13
Control de Procesos
Dx(t) =[empty sym]
OTRA APLICACION
EDO de segundo orden + Matlab
Solución de una ecuación diferencial de segundo orden por el método de Euler y su implementación
en Matlab
Euler EDO Orden 2 + Matlab
1. [Link]
2. [Link]
14
