Scribd
Cargar
33 vistas6 páginas

Examen Parcial 2

Este documento presenta un examen parcial de Introducción a Matlab con 5 problemas que deben resolverse usando código de Matlab. El examen cuenta con un total de 100 puntos y debe entregarse el 28 de febrero de 2023. Los problemas incluyen resolver ecuaciones diferenciales, encontrar raíces de una ecuación no lineal, resolver un sistema de ecuaciones lineales, hallar transformadas de Laplace y analizar la estabilidad de un sistema de control.

Cargado por

joski ortiz
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
33 vistas6 páginas

Examen Parcial 2

Este documento presenta un examen parcial de Introducción a Matlab con 5 problemas que deben resolverse usando código de Matlab. El examen cuenta con un total de 100 puntos y debe entregarse el 28 de febrero de 2023. Los problemas incluyen resolver ecuaciones diferenciales, encontrar raíces de una ecuación no lineal, resolver un sistema de ecuaciones lineales, hallar transformadas de Laplace y analizar la estabilidad de un sistema de control.

Cargado por

joski ortiz
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
Está en la página 1/ 6

Introduccion al Matlab

Segundo Examen Parcial


Universidad Fidelitas
Facultad de Ingeniería

Nombre Puntos
Obtenidos

Porcentaje de Nota
la Nota:

Total de Puntos: 100


Fecha entrega: 28 de febrero de 2023
Porcentaje de la Nota: 30%
Tiempo máximo: 3 horas
Nota: El examen es individual. Incluir graficas y códigos necesarios para la
resolución del problema. Solo se aceptan consultas respecto al formato del
examen o algún código de error. Se puede subir el documento en formato
Word o PDF

II Parte. Desarrollo
Resolución en la casa (Fecha de Entrega 11 de Marzo 2023)
Total de puntos 90%
Realice en Matlab una secuencia de código para resolver lo solicitado. Se
debe copiar el código generado en el espacio en blanco. En caso de
graficación incluir la imagen.
Problema 2. 1. Total de puntos 18 puntos (6% cada acierto). Proceda a
resolver las ecuaciones diferenciales con el comando dsolve en Matlab
'' '
a ) y =cos ( x ), y (0 )=−1 , y (0)=3
b ) y ' −5 y =exp( x )
'' ' ''
c ) y −5 y +9 y =ln(t )

%Ejercicio 2.1
disp('Ejercicio 2.1')
a=dsolve('D2y=cos(x)','Dy(0)=3','y(0)=-1')
b=dsolve('Dy-5*y=exp(x)')
c=dsolve('D3y-5*D2y+9*y=log(t)')

Problema 2. 2. Total de puntos 18


Dada la ecuación no lineal, resuelva con fzero y con el valor inicial de x0=1

x 2−4 tan( x )+exp(6 x )−7 x =0


%Ejercicio 2.2
disp('Ejercicio 2.2')
fa=@(x)x^2-4*tan(x)+exp(6*x)-7*x
xo=1;
x=fzero(fa,xo)

Problema 2.3 Total de puntos 18 Dado el sistema de ecuaciones:

4 a+6 b−6 c +7 d +9 e=−10


10 a−9 b−13 c +5 d−8 e=9
−8 a+10 b−15 c +9 d−2 e=4
5 a+3 b+4 c +8 d−4 e=20
5 a+10 b+20 c +5 d−6 e=6
%Ejercicio 2.3
Matriz_A=[4 6 -6 7 9;10 -9 -13 5 -8;-8 10 -15 9 -2;5 3 4 8 -4;5 10 20 5 -
6];
Matriz_B=[-10;9;4;20;6];
Matriz_Final=linsolve(Matriz_A,Matriz_B);
a=Matriz_Final(1,1)
b=Matriz_Final(2,1)
c=Matriz_Final(3,1)
d=Matriz_Final(4,1)

Página 2 de 6
=Matriz_Final(5,1)

e Proceda a resolver empleando el comando linsolve()


Problema 2. 4. Total de puntos 18. (6% cada acierto).
Halle la transformada de Laplace de las siguientes funciones en el dominio
del tiempo

a )t 3 sen(6 t )
b ) e−3 t *cos ( t )
c )t 3 −4 t 2 +5 t−6

%Ejercicio 2.4
syms t positive
a4=laplace(t^3*sin(6*t))
b4=laplace(exp(-3*t)*cos(t))
c4=laplace(t^3-4*t^2+5*t-6)

Problema 2. 5. Total de puntos 18


Sea un sistema de control con una planta P, controlador C y una función de
lazo cerrado H dadas por las funciones de transferencia:
5
P=
s3 −9 s2 −9 s−10
C=s−10
H=1
a) Aplique y obtenga la grafica de la respuesta al escalón y al impulso en Matlab
únicamente para la planta P. Comente sobre la estabilidad.
%Ejercicio 2.5. Parte a
P=tf([5],[1 -9 -9 -10])
C=tf([1 -10],[0 1]);
H=1;
figure(1) , subplot(2,1,1) , step(P) , grid on , title('Escalón')
subplot(2,1,2) , impulse(P) , grid on , grid on , title('Impulso')
pole(P)

Página 3 de 6
El sistema de la planta muestra un comportamiento inestable, al usar el
comando pole(), notamos que uno de sus polos está en la parte positiva, esto
nos confirma la inestabilidad de la planta.

b) Aplique y obtenga la grafica de la respuesta al escalón y al impulso en Matlab


para el sistema en lazo cerrado. El controlador C esta en serie con la planta P y el
bloque serie está en retroalimentación con la función H. Comente sobre el efecto
de la estabilidad en el nuevo sistema
%Ejercicio 2.5. Parte b
P=tf([5],[1 -9 -9 -10])
C=tf([1 -10],[0 1]);
H=1;

FT1=series(C,P);
FT2=feedback(FT1,1,-1)

figure(2) , subplot(2,1,1) , step(FT2) , title('Escalon 2') , grid on


subplot(2,1,2) , impulse(FT2) , title('Impulso 2') , grid on

Página 4 de 6
Este sistema con retroalimentación, al igual que la planta presenta
inestabilidad, podemos comprobarlo con el comando pole(), al ver que posee
un polo en la parte positiva

%Ejercicio 2.5. Parte b


P=tf([5],[1 -9 -9 -10])
C=tf([1 -10],[0 1]);
H=1;

FT1=series(C,P);
FT2=feedback(FT1,1,-1)

figure(2) , subplot(2,1,1) , step(FT2) , title('Escalon 2') , grid on


subplot(2,1,2) , impulse(FT2) , title('Impulso 2') , grid on

pole(FT2)

Página 5 de 6
Suerte y muchos éxitos

Página 6 de 6

También podría gustarte