Deber 2 Sys

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 44

1. Para las siguientes ecuaciones diferenciales determine la respuesta.

a) y ' ( t ) + y ( t )=2 x ' ( t ) + x ( t ) , x ( t )=4 e−2 t u ( t ) , y ( 0 )=2

 Hallar la entrada del sistema


2 (−8 e−2 t u ( t ) + 4 e−2 t δ ( t ) ) + 4 e−2 t u (t )
−16 e−2 t u ( t ) +8 e−2 t δ ( t ) + 4 e−2 t u (t )
−12 e−2 t u ( t ) +8 e−2 t δ ( t )

 Hallar la respuesta al impulso


'
y ( t ) + y ( t )=0 , y ( 0 ) =1
s+1=0
−t
s=−1 y= A 1 e
 Aplicamos condición inicial
 y ( 0 )= A 1 e−t A1=1
 h ( t )=e−t u ( t )
 Aplicamos Convolución

y ( t ) =∫ ( −12e−2 τ u ( τ )+ 8 e−2 τ δ ( τ ) ) ( e−( t−τ ) u ( t−τ ) ) dτ
−∞

y 1 (t )= ∫ −12 e−2 τ u ( τ ) ∙e−t +τ u ( t −τ ) dτ
−∞
t
t
y 1 (t )=−12e
−t
∫ e− τ dτ=12 e−t|e−τ|0
0

y 1 (t )=12 e ( e −1 ) =12e−2 t −12 e−t


−t −t


y 2 ( t )= ∫ 8 e
−2 τ −t+ τ
δ (τ ) e u ( t−τ ) dτ
−∞
t
y 2 ( t )=8 e
−τ
∫ e−τ δ ( τ ) dτ
−∞

 Aplicamos la propiedad

∫ f ( t ) δ ( t−a ) dt=f ( a )
−∞

y 2 ( t )=8 e ( e )=8 e
−t −0 −t

y ( t ) = y 1+ y 2
−2 t −t −t
y ( t ) =12e −12 e +8 e
−2 t −t
y 0 ( t ) =12e −4 e

 Encontrar la respuesta a entrada cero


y T ( t )=Y P ( t ) +Y C ( t )
'
y ( t ) + y ( t )=0 y ( 0 )=2
s=−1
−t
y h ( t ) =A 1 e

 Aplicamos condiciones iniciales


y ( 0 )= A 1 e 0=2 A1=2
−t
y h ( t ) =2 e
 La respuesta es
−2 t −t −t
y ( t ) =12e −4 e +2 e
 Simplificamos
−2 t −t
y ( t ) =12e −2 e
 Graficamos

b) y ' ' ( t )+ 3 y ' ( t ) +2 y ( t )=2 x ' ( t ) + x ( t ) , x ( t )=4 u ( t ) , y ( 0 )=2 y ' ( 0 )=1

 Hallar la entrada del sistema


2 ( 0 u ( t ) + 4 δ ( t ) +4 u ( t ) +8 δ ( t ) +4 u ( t ) )
 Hallar la respuesta al impulso
y ' ' ( t )+ 3 y ' ( t ) +2 y ( t )=0 , y ( 0 ) =1 y ' ( 0 )=1
2
s +3 s +2=0
−3 ± √ 3 −4 ( 1 )( 2 )
2
s=
2 (1 )
−3+1
s1= →−2
2
−3−1
s2= →−1
2

−2 t −t
h ( t )= A1 e + A2 e
 Aplicamos condición inicial
h ( 0 ) =A 1 + A 2=0
'
h ( 0 )=−2 A1 −A 2=1
− A1=1→ A 1=−1
−1+ A2=0→ A2=1
h ( t )=−e−2 t + e−t u ( t )
 Aplicamos Convolución

y ( t ) =∫ ( (−e ) u ( t−τ )) ∙ ( 8 δ ( τ ) +4 u ( τ ) ) dτ
−2 t +2 τ −t +τ
+e
−∞
∞ t
y 1= ∫ −e ∫ e 2 t δ ( τ ) dτ
−2 t +2 τ −2 t
∙8 δ ( τ ) u (t−τ ) dτ=−8 e
−∞ 0

 Aplicamos la propiedad

∫ f ( t ) δ ( t−a ) dt=f ( a )
−∞

y 1=−8 e
−2 t
( e0 ) ; y 1=−8 e−2t


y 2= ∫ −4 e
−2 t+ 2 F
∗u ( t−T ) u ( T ) dT
−∞

t t
y 2=∫ −4 e ∗e dT ; y 2=−4 ∫ e
−2 t 2T −2 t 2T
∗e dT
0 0

t
y 2=−4 e
−2 t
∫ e 2 T dT ; y 2=−2e−2 t|e 2 t|0t
0

y 2=−2 e
−2 t
( e 2 T −1 ) ; y 2=−2+ 2 e−2 t


y 3= ∫ 8 e
−t +2 T
∗u ( t−T ) δ ( T ) dT
−∞

t
y 3=8 e
−t
∫ e T δ ( T ) dT ; y 3=8 e−t ( e 0 )
−∞

−t
y 3=8 e

y 4 =∫ 4 e
−t +T
∗u ( t−T ) u ( t−T ) dT
−∞

t
y 4 =4 e
−t
∫ eT dT ; y 4 =4 e−t|e T|0t
0

y 4 =4 e−t|e T −1|; y 4=4−4 e−t

 Se procede a sumar todas las respuestas obtenidas y 1 , y 2 , y 3 , y 4

+ (−2+2 e ) + 8 e−t +4−4 e−t


−2 t −2 t
y n ( t ) =−8 e
−2 t −t
y n ( t ) =−6 e +2+ 4 e ( ESTADO CERO )

 Se procede a ecnontrar la respuestade estadocero

y ' ' ( t )+ 3 y ' ( t ) +24 (t )=0


y ( 0 )=2

y ' ( 0 )=1

2
s +3 s +2
s1=−1 s2=−2

−t −2 t
y c (t )= A 1 e + A 2 e
−t −2 t
y c ' ( t )=−A 1 e −2 A2 e

 Se aplica condiciones inicales y resolvemos el sistema de ecuaciones


y c ( 0 )= A1 + A 2=2

y 'c(t )= A1−2 A2 =1

− A∗A 2=3 ; A2 =3

A1−3=2 ; A 1=5

 Y reemplazando nos quedara


−t −2 t
y c (t )=5 e −3 e u (t )

 Se procede a sumar las respuestas obtenidas y asi obtener el total


y ( t ) = y c ( t )+ y n ( t )

y ( t ) =5 e−t −3 e−2t −6 e−2t +2+4 e−t u (t )


−t −2 t
y ( t ) =9 e −9 e +2 u ( t )

 Su grafica

c) y ' ' ( t )+ 2 y ' ( t ) +2 y ( t )=x ( t ) +2 x ( t −1 )


x ( t )=4 u ( t )
y ( 0 )=0 y ' ( 0 ) =0

 Reemplazamos en x ( t )
'' '
y ( t )+ 2 y ( t ) +2 y ( t )=4 u ( t )+ 8u ( t−1 )
y ( 0 )=0 y ' ( 0 ) =0

 La respuestadel escalon sera


'' '
y ( t )+ 2 y ( t ) +2 y ( t )=u ( t )
2
s +3 s +2

−2 ± √ 4−4 ( 1 ) ( 2 ) −2± √ −4 −2 ±2
x= ; x= ; x=
2 2 2
x=−1 ± 1 j
−t
y n ( t ) =e ¿

 y asi podemos hallar la otra parte de larespuesta

y p ( t ) =k y 'p=0 y 'p' ( t )=0

0+2 ( 0 ) +2 k =1u ( t )
1
y p ( t ) =k= u (t)
2
 Sumamos ambas respuestas
y t ( t )= y h ( t ) + y p ( t )

−t 1
y t ( t )=e ( A 1 sin ( t )+ A 2 cos ( t ) ) + 2 u( t)
 Reemplazamos condiciones iniciales
1
y t ( 0 ) =e−0 ( A 1 sin ( 0 )+ A 2 cos ( 0 ) ) + u ( 0 )=0
2

1
1 ( A 1 (0)+ A2 (1) ) + =0
2

1
A2 (1)+ =0
2

−1
A 2=
2
 Reemplazamos y derivamos

y t ( t )=e−t ( A 1 sin ( t )−0.5 cos ( t ) ) +0.5 u(t )

y t ' ( t )=−e
−t
( A 1 sin ( t )−0.5 cos ( t ) ) +e−t ( A 1 cos ( t ) +0.5 sin ( t ) )

y t ' ( 0 )=−e−0 ( A1 sin ( 0 )−0.5 cos ( 0 ) ) +e−t ( A 1 cos ( 0 )+ 0.5 sin ( 0 ) )=0

(−1 ) ( A1 ( 0 )−0.5 ( 1 ) ) + ( A 1 ( 1 ) +0.5 ( 0 ) ) =0


0.5+ A 1=0

A1=−0.5
 Reemplazamos para hallar la respuesta al escalón
y t ( t )=e−t (−0.5sin ( t )−0.5 cos (t ) ) +0.5 u( t)
 Hallamos la respuesta final
y ( t ) =4 ( y t ( t ) ) + 8( y t ( t−1 ) )
y ( t ) =[−2 e ( sin ( t ) +cos ( t ) ) +2 ] u ( t )+ [−4 e ( sin ( t −1 ) +cos ( t −1 ) ) + 4 ] u ( t −1 )
−t −t+1

 Su Gráfica

2. Un sistema en tiempo continuo tiene una respuesta que es el cubo de su entrada.


Clasifique el sistema de acuerdo con las siguientes características: linealidad,
invariancia en el tiempo, estabilidad, causalidad, memoria e invertibilidad. De
una explicación breve del porqué de sus respuestas.
Un sistema cuya respuesta es el cubo de la entrada seria:
3
y ( t ) =( x (t))
Linealidad:
Dadas las siguientes entradas X1(t) y X2(t), se determina las salidas y1(t) y y2(t).
3
y 1 ( t )=( x 1(t))

y 2 ( t )=( x 2(t ))3


Dada la entrada x3(t)
x 3 ( t )=( x 1 ( t ) ) +(x 2(t))

La salida se obtendría
y 3 ( t )= y 1 ( t ) + y 2 ( t )
3 3
y 3 ( t )=(x 1(t )) +( x 2(t))

y 3 ( t )=( x 3 (t))3
3
y 3 ( t )=( ( x 1 ( t ) ) +(x 2(t )))
3 3 3 3
y 3 ( t )=( x 1 ( t ) ) +3 x 1 ( t )∗( x 2 ( t ) ) +3 ( x 1 ( t ) ) ∗x 2 ( t ) + ( x 2 ( t ) )

Por lo tanto
y 3 ( t ) ≠ y 1 ( t ) + y 2 (t )
Dado lo siguiente deducimos que el sistema no es lineal

Invariancia en el tiempo
Sea
x 1(t)

y 1 ( t )=( x 1(t ))3


Y dado que
x 2 ( t )=( x 1 ( t −¿ ) )

y 2 ( t )=( x 2(t ))3=(x 1 ( t−¿ ) )3


y 2 ( t )= y 1 ( t−¿ )
Por lo tanto el sistema es invariante en el tiempo.
Estabilidad
Un sistema es estable cuando su salida está acotada. Es decir, que su salida no es ±∞ sino
un valor concreto
−a ≤ x (t) ≤ a
Por lo tanto su salida

−a 3 ≤ y (t ) ≤ a3
Podemos concluir que estará limitada dentro de valores finitos y que no denota una
tendencia a infinito por lo cual el sistema es estable
Causalidad
Podemos denotar que en la salida de y(to) en un determinado tiempo to este dependerá
solamente del valor de la entrada x(to) en el tiempo to, por lo tanto este sistema dependerá
solo de valores en el presente y podemos determinar que es causal.
Memoria
Dado que el sistema depende solo de valores en tiempo presente, al no depender del tiempo
pasado el sistema no tiene memoria.
Invertibilidad

Si tenemos una entrada


x 1(t)
De manera que
3
y 1 ( t )=( x 1(t))
Y también
x 2(t)
Tal que

y 2 ( t )=( x 2(t ))3


Siempre que
x 1 ( t ) ≠ x 2(t )
Entonces
y 1 ( t ) ≠ y 2(t)
Para este sistema todas las entradas resultan en valores únicos en la salida por lo tanto el
sistema es invertible

3. Un sistema en tiempo continuo se describe mediante la ecuación diferencial:

t y ' ( t )−8 y ( t )= x(t )

Clasifíquelo de acuerdo con su linealidad, invariancia en el tiempo y estabilidad.


Explique las razones de sus respuestas.

a) Linealidad

Para determinar si un sistema es lineal, debemos analizar su homogeneidad y su


aditividad, y si cumple con ambas propiedades pues el sistema será lineal.

 Homogeneidad

Si tenemos una señal de entrada

x 1 ( t )=g(t )

Entonces:

t y 1 ' ( t ) −8 y 1 ( t )=g(t ) (1)

Para comprobar su homogeneidad debemos multiplicar la señal de entrada


por una constante K y comprobar que su salida sea igual a la salida y 1
multiplicada por la contante K.

x 2 ( t )=K g(t )

Entonces:

t y 2 ' ( t )−8 y 2 ( t )=K g(t ) (2)

Ahora multiplicamos el primer sistema (1) por K, y obtenemos.


Kt y 1 ' ( t )−K 8 y 1 (t )=Kg(t) (2)

Igualamos las (2) y (3) y obtenemos lo siguiente.

Kg(t)=Kg(t)

Kt y 1 ' ( t )−K 8 y 1 (t )=t y 2 ' ( t )−8 y 2 ( t )

Esta igualdad únicamente puede ser verdadera si

K y 1 ( t )= y 2 ( t )

Como podemos observar el sistema se multiplicado por el factor K, al igual


que la entrada, por lo tanto, se concluye que se trata de un sistema
homogéneo.

 Aditividad

Si tenemos las señales de entrada

x1( t )
x2 ( t )

Entonces tendremos los sistemas:

t y 1 ' ( t ) −8 y 1 ( t )=x 1 (t) (1)

t y 2 ' ( t )−8 y 2 ( t )=x 2 (t) (2)

Para comprobar su aditividad verificar que el sistema responde de igual


manera a una señal compleja que a señales simples que se desprenden de la
original. Por lo tanto, ahora tendremos una señal de entrada

x 3 ( t )=x 1 (t ) + x 2 ( t )

Y obtenemos:
t y 3 ' ( t )−8 y 3 ( t )=x 3 (t )

t y 3 ' ( t )−8 y 3 ( t )=x 1 ( t ) + x 2 ( t ) (3)

Ahora si sumamos (1) y (2), obtenemos:

t y 1 ' ( t ) −8 y 1 ( t ) +t y 2 ' ( t )−8 y 2 ( t )=x 1 (t)+ x 2 ( t )

t [ y1 ' ( t )+ y 2 ' ( t ) ]−8 [ y 1 ( t )+ y 2 ( t )]=x 1 (t)+ x 2 ( t ) (4)

Finalmente igualamos (4) y (3):

x 1 ( t ) + x 2 ( t )=x 1 ( t )+ x2 ( t )

t [ y1 ' ( t )+ y 2 ' ( t ) ]−8 [ y 1 ( t )+ y 2 ( t )]=t y 3 ' ( t )−8 y 3 ( t )

t [ y1 ( t )+ y 2 ( t ) ]' −8 [ y 1 ( t ) + y 2 ( t )]=t y 3 ' ( t )−8 y 3 ( t )

Esta igualdad únicamente puede ser verdadera si

y 1 (t ) + y 2 ( t )= y 3 ( t )

Por lo tanto, se concluye que se trata de un sistema aditivo.

Una vez analizado el sistema se obtuvo que posee homogeneidad y aditividad.

Por lo tanto, se concluye que se trata de un sistema lineal.

b) Invariancia en el tiempo

Si tenemos

x 1 ( t )=x ( t )
t y 1 ' ( t ) −8 y 1 ( t )=x 1 (t)

Pero si existe un desplazamiento en el tiempo

t → t−a

Entonces el sistema resultante seria:

( t−a ) y 1 ' ( t−a )−8 y 1 ( t−a )=x ( t−a ) (1)

Ahora comprobamos su invariancia en el tiempo analizando el sistema para otra


entrada.

x 2 ( t )=x ( t−a )

t y 2 ' ( t )−8 y 2 ( t )=x (t−a) (2)

Si igualamos (1) y (2), obtenemos:

x (t−a)=x (t−a)

( t−a ) y 1 ' ( t−a )−8 y 1 ( t−a )=t y 2 ' ( t ) −8 y 2 ( t )

Podemos ver que la igualdad no se cumple, por lo tanto, se trata de un sistema


variante en el tiempo.

c) Estabilidad

Para analizar la estabilidad del sistema, en primer lugar, hallaremos su respuesta


homogénea

t y ' ( t )−8 y ( t )= x(t )

t y ' ( t )−8 y ( t )=0

Obtenemos una ecuación diferencial de variables separables:


t y '=8 y

dy
t =8 y
dt

1 8
dy= dt
y t

1 8
∫ y dy =∫ t dt

ln y=8 ln t

y ( t ) =K t 8

Como podemos ver, la respuesta del sistema para entrada cero, aumentará sin
límite.

Por lo tanto, el sistema es Inestable.

4. Para un circuito RLC serie. Suponga valores de R= 1 ohm, y L = 0.25 Henrios.


Encuentre el valor crítico de capacitancia C que haría que la respuesta pase de
sobreamortiguada (sin oscilaciones, raíces reales diferentes), a críticamente
amortiguada (raíces repetidas) y a subamortiguada (con oscilaciones, raíces
complejas). Usando MATLAB grafique la respuesta matemática para diferentes
valores de C (al menos 5 valores distintos) que ejemplifiquen los tipos de
comportamiento. Comente los resultados.

R=1 Ω
L=0.25 H

Aplicamos Ley de voltajes de Kirchoff


VL+VR+VC=V
t
di 1
L + Ri + ∫ idt +VC ( ¿ ) =V
dt C ¿
2
d i di i dV
L +R + =
dt dt C dt
2
d i R di i 1 dV
+ + =
dt L dt LC L dt
Si la fuente es continua entonces
dV
=0
dt
2
d i R di i
+ + =0
dt L dt LC
La ecuación característica
2 Rs 1
s+ + =0
L LC
2 2
s +2 ∝ s+ ( w0 ) =0

Donde
R
∝=
2L
w 1
0=¿ ¿
√ Lc

Calculamos las raíces de la siguiente manera

S=

−2 α ± 4 α 2−4 ( w0 )
2

s 1=−α+ √ α 2−w 20 s 2=−α− √ α 2−w20

Sobreamortiguada
2 2
α > w0

R2 1
>
4L 2
LC
2
R 1
>
4 (0.25) (0.25) C
2

C> 1 F

Críticamente amortiguada
2 2
α =w 0

R2 1
=
4 L LC
2

2
R 1
=
4 (0.25)2
(0.25) C

C=1 F

Subamortiguada

α 2< w20

R2 1
<
4 L LC
2

2
R 1
<
4 (0.25)
2
(0.25) C

C< 1 F
Usando MATLAB grafique la respuesta matemática para diferentes valores de C (al menos
5 valores distintos) que ejemplifiquen los tipos de comportamiento. Comente los resultados.
C=1

C=5
C= 2

C=0.1
C=0.01

Cuando el sistema tiene una Capacitancia de 1 Faradio (C=1) se encuentra en un estado


estático, es decir, que cualquier variación en la Capacitancia el sistema pasaría a ser
sobreamortiguado (aumento) si C<1F, o subamortiguado (disminución) C >1F como
podemos notar en las graficas presentadas anteriormente.

5. Para cada uno de los siguientes circuitos:


a) Obtener la o las ecuaciones diferenciales que lo representen.
b) Encontrar la variable solicitada y trazar una gráfica aproximada a mano.
c) Graficar la respuesta en MATLAB.
d) Realice la simulación de los circuitos usando algún programa de simulación (por
ejemplo: Multisim) y obtener la respuesta transitoria.
e) Haga una comparación entre las respuestas calculadas a mano, su gráfica en
MATLAB y las obtenidas en el simulador de circuitos.
CIRCUITO 1
−¿¿
Halle i (t ) para t >0. El circuito está en el estado estable t=0 .
a)

Sumamos las dos resistencias en paralelo ,la de 2 y la de 11

( ) ( )
−1 −1
1 1 13 22
+ = → Ω
2 11 22 13
Convertimos de fuente corriente a voltaje
22 88
4∙ = V
13 13

Volvemos atransformar la fuente de corrientea voltaje


4 V ∙ 3 V =12 V
La corriente en este circuito es cero por lo que el voltaje en el condensador , es iagual al de la fuente
244
por lo tanto V C ( 0 )= V
13
Para t >0

debido a que el circuito estabaestable antes de cerrar en interruptor , el condensador estaba


cargado y por lotanto tiene voltaje.
Por lo tanto procederemos a dibujar una fuente de voltaje en serie para representar el voltaje
inicial .
Como V 1 , esta cortocircuitado, V 1=0

Pasamos al Dominio de la frecuencia


1 13 244 244
F→ V→
13 s 13 13 s
13 244
4 I+ I+ =0
s 13 s

(
I 4+
s)
13 −244
=
13 s
−244 −244
13 s 13 s −244 s
I= → =
13 4 s+13 13 s ( 4 s +13 )
4+
s s
−244
I=
52 s+169

i (t )=L
−1
{52−244
s ++169 }

{ }
−1
−244 1
¿ ∙L
52 13
s+
4
−13t
−244 4
i (t )= ∙e
52
−13t
−244 4
b) i (t )= ∙e
52

c) Matlab
d) Multisim
Circuito #1: Descripción o
comentario:
El interruptor está en
posición de reposo del
tiempo 0 al tiempo
0.5s, luego cambia de
posición hasta el
tiempo 5s, luego del
cual regresa a su
posición de reposo.

La corriente I(C1)
corresponde a la
corriente requerida en
el ejercicio. Desde el
tiempo 0 hasta el
tiempo 0.5 no existe
corriente. Desde 0.5 en
adelante podemos ver
la respuesta transitoria.

En este caso se analiza


únicamente hasta el
tiempo 4s, por lo tanto,
no se visualiza el
tiempo en el que el
interruptor regresa a su
posición de reposo

e) Podemos observar la similitud en las gráficas obtenidas en Matlab y la respuesta


transitoria en Multisim.
Podemos observar que cuando el interruptor es movido, la magnitud de la corriente
empieza a disminuir a cero, empezando desde un valor negativo.
CIRCUITO 2
El interruptor estuvo en A muchotiempo antes de pasar a B . HalleV C ( t ) e i L ( t ) para

t >0

a) Cuando el condensador esta estable no circula corriente y el inductor es


un cortocircuito .
transformamos de fuente de corriente a voltaje

4 mA ∙3 KΩ=12 V

El voltajeV C inicial es V C ( 0 ) =12V y la corriente inicial en el conductor

es igual a lafuente de corriente i L ( 0 ) =4 mA .


50 mH → sL=0.05 s
1 1 16
µF → = −6
16 50 10 s

12
12 V =
s

s
=
s(
12 16 x 106 3
+8 x 10 I C ( s ) )
12
I C ( s) =
( )
6
16 x 10 3
s +8 x 10
s
12
I C ( s) = 6 3
16 x 10 + 8 x 10 s
6
12 16 x 10
V C ( s )= − I C ( s)
s s
6
12 192 x 10
¿ −
s s ( 8 x 103 s+16 x 106 )
6
12 192 x 10
¿ −
s 8 x 10 3 s2 +16 x 106 s

12 ( 8 x 10 s+ 16 x 10 ) −192 x 10
3 6 6
V C ( s )= 3 2 6
8 x 10 s +16 x 10 s

3 6 6
( ) 96 x 10 s+192 x 10 −192 x 10
VC s =
8 x 10 3 s2 +16 x 106 s

3 3 3
96 x 10 s 96 x 10 s 96 x 10
→ →
8 x 10 s + 16 x 10 s s ( 8 x 10 s+16 x 10 ) 8 x 10 s+16 x 10
3 2 6 3 6 3 6

12
¿
s+ 2000

V C ( t ) =L
−1
{ s+122000 }
−2000 t
V C ( t ) =12e

V L=−V R 4 K I L =I R 4 K
V L ( s )=0.05 s I L ( s )
0.05 s I L ( s )=−I L ( s ) ∙ 4 x 103
0.05 s I L ( s ) + I L ( s ) ∙ 4 x 103 =0
Tiempo

d iL(t )
0.05 +4000 ∙i ( t )=0
dt

d iL(t) 3
+80 x 10 ∙i ( t ) =0
dt

d iL (t) 3
=−80 x 10 ∙ i ( t )
dt

d iL(t) 3
=−80 x 10 ∙ dt
it( )

1
∫ i ( t ) ∙ dt=−∫ 80 x 103 ∙ dt
ln |i ( t )|=−80 x 103 t+ C
3
−80 x10 t c
i L ( t )=e ∙e

i L ( 0 )=4 x 10−3 =e c
3
−3 −80 x 10 t
i L ( t )=4 x 10 ∙ e
b)
c) Matlab

d) Multisim
Circuito #2: Descripción o
comentario:
El interruptor está en
posición de reposo del
tiempo 0 al tiempo
50ms, luego cambia de
posición hasta el
tiempo 10s, luego del
cual regresa a su
posición de reposo.

La corriente I(L1) y el
voltaje V (4)
corresponden a los
datos requeridos en el
ejercicio. Desde el
tiempo 0 hasta el
tiempo 50ms se
observa el estado
estable del circuito.
Desde 50ms en
adelante podemos ver
la respuesta transitoria.

En este caso se analiza


únicamente hasta el
tiempo 200ms, por lo
tanto, no se visualiza el
tiempo en el que el
interruptor regresa a su
posición de reposo

e) En las gráficas simuladas y dibujadas, podemos observar como comienzan


en valores positivos al inicio ( t=0 ), debido a que estos elementos están
cargados y luego poco a poco van decreciendo hacia 0, cuando terminan de
descargarse.
CIRCUITO 3
determine v ( t ) para t >0 cuando V 1=( 25 sin ( 4000t ) ) u ( t ) . Suponga estado
−¿ .¿
estable en t=0

a) Para t <0
La bobina se carga

Transformamos la fuente de corriente a voltaje


5 mA ∙6 Ω=30 mV

24
V ( 0 )=30 mV ∙ =24 mV
24+ 6
30 mV
i L ( 0 )= =1 mA
30 Ω
Para t >0
V ( t ) =i ( t ) ∙ 24
V ( t ) =V 1 ( t )−V L ( t )

Ldi ( t )
¿ V 1 (t ) −
dt
Ldi ( t )
24 ∙ i ( t ) + =25 sin ( 4000 t )
dt
4000
24 I ( s )+ LsI ( s )−L ( 0 ) =25∙ 2 2
s + 4000
3
−3 100 x 10
I ( s ) ( 24+ Ls )−L x 10 =
s 2+16 x 106

100 x 103 + L x 10−3


I ( s ) ( 24+ Ls )=
s 2 +16 x 106

100 x 103 + L x 10−3


I ( s) =
( s 2+16 x 106 ) ( 24+ Ls )

V ( s )=24 ∙ I ( s )
3 −3
100 x 10 + L x 10
V ( s )=24 ∙
( s 2 +16 x 106 ) ( 24+ Ls )
100 x 103
¿ 24 ∙
( s 2+16 x 106 ) ( 24+ 8 x 10−3 s )
100 x 103 ≫ L x 10−3
6
300 x 10
V ( s )=
( s +16 x 106 ) ( 3000+ s )
2

−12 s+36000 12
¿ +
2
s +16 x 10
6
s +3000

36000 12 12 s
¿ + − 2
s +16 x 10 s+3000 s +16 x 106
2 6

36000
V (t)= sin ( 4000 t )+12 e−3000t −12 cos ( 4000 t )
4000

V ( t ) =9 sin ( 4000 t )+12 e−3000t −12 cos ( 4000 t )


b)

c) Matlab
d) Multisim
Circuito #3: Descripción o
comentario:
El interruptor está en
posición de reposo del
tiempo 0 al tiempo
0.5ms, luego cambia de
posición hasta el
tiempo 5s, luego del
cual regresa a su
posición de reposo.
El voltaje V (2)
corresponde al voltaje
requerido en el
ejercicio. Desde el
tiempo 0 hasta el
tiempo 0.5ms se grafica
el estado estable del
resistor. Desde 0.5ms
en adelante podemos
ver la respuesta
transitoria.

En este caso se analiza


únicamente hasta el
tiempo 10ms, por lo
tanto, no se visualiza el
tiempo en el que el
interruptor regresa a su
posición de reposo

e) S i comparamos, ambas graficas son idénticas, podemos observar el mismo


comportamiento oscilante.
Es así porque el circuito queda alimentado por una fuente AC por lo que
realmente no tiene estado estable y sus valores siempre están variando de
forma periódica.
CIRCUITO 4
−¿ .¿
Determine i ( t ) para t >0. Suponga condiciones de estado estable en t =0

3 kΩ∗1 A=3 V
1 Vk
8 V =Vc +3 V
Vc=5 V
V C (0)=5V
i c (0) =0

8V
=i
1kΩ L
i L (0 )=8 mA

a. t >0

Ya estaba cargados los elementos al mover el interruptor

8−5=V L +V C +V 2 kΩ

d ic (t ) 1
3=L + ∫ i dt +Z k ∗i ( t )
dt c ( t )
Si se asume que la corriente presentada en elinductor no tiende a variar
d iL ( t )
Entonces será =0
dt
1
3= ∫ i dt+ 2 k∗i ( t )
c ( t)
3 I ( s)
=
5 sC
+2 k∗I ( s ) ; 3=I ( s )
1
sC ( ) (1
+ 2000 s ; 3=I ( s ) + 2000
C )
1
∗3
I ( s) =
( 9
3
3 x 10 +2000 s
; I ( s )
)
=
2000
5+1.5 x 10
6

i ( t )=
3
2000
∗L−1
{ 1
5+1.5 x 10
6
; i (t)=
} 3
2000
∗e−15000 t

b.

c. Matlab
d. Multisim
Circuito #4: Descripción o
comentario:
El interruptor está en
posición de reposo del
tiempo 0 al tiempo
50us, luego cambia de
posición hasta el
tiempo 8s, luego del
cual regresa a su
posición de reposo.
La corriente I(R2)
corresponde a la
corriente requerida en
el ejercicio. Desde el
tiempo 0 hasta el
tiempo 50us se
presenta el estado
estable del resistor.
Desde 50us en adelante
podemos ver la
respuesta transitoria.

En este caso se analiza


únicamente hasta el
tiempo 4s, por lo tanto,
no se visualiza el
tiempo en el que el
interruptor regresa a su
posición de reposo

e. Se puede observar que tanto en la grafica simulada como en la que obtuvimos


manualmente se puede ver que se llega a un estado estable siendo la corriente 0
Como diferencia presentada en las graficas es que se puede asumir la corriente
encontrada en el inductor y no tiende a variar luego de que se dé la conmutación de
switch.
Aparte lo presentado no siempre será 100% correcto y la razón de eso es que debido a la
disminución de energía que se presenta en la bobina se buscara descargarse, pero no del
todo ya que gracias que un tiene una fuente DC conectada buscando así volver a cargarse
provocando que la corriente disminuya y vuelva a aumentar hasta que quede estabilizada
nuevamente, y terminaría afectando a la corriente en el condensador

6. Para las siguientes ecuaciones diferenciales:


I. Determine la respuesta al impulso.
II. Grafique la respuesta en MATLAB.
III. Con la respuesta al impulso obtenida encuentre la salida del sistema para la entrada
indicada usando convolución y grafique la respuesta en MATLAB.
IV. Realice el diagrama de bloques de cada ecuación diferencial e impleméntelo en
SIMULINK y grafique la señal de salida.
V. Compare los resultados de la convolución graficada en MATLAB con los de la
simulación.

a) y ' ( t ) +5 y ( t )=x ( t ) . Con x ( t )=e−t u ( t )

I. h' ( t ) +5 h ( t )=δ ( t )
sH ( s )+ 5 H ( s )=1
H ( s )( s+5 )=1
1
H ( s )=
s +5
h ( t )=L { }
−1 1
s +5
→e
−5 t

h ( t )=e−5 t

II. En Matlab
III. y ( t ) =x ( t )∗h ( t )
t t
y ( t ) =∫ e u ( τ ) ∙e dτ=∫ e e
−τ −5 ( t −τ ) −τ −5 t 5τ
e dτ
0 0
t
∙∫ e dτ
−5 t 4τ
¿e
0

| |
t
−5 t e4τ
¿e
4 0
−5 t 4 t
e e e−5 t
¿ −
4 4

1
y ( t ) = ( e−t−e−5 t )
4

IV. y ' ( t ) +5 y ( t )=x ( t )


sY ( s )+ 5Y ( s )= X ( s )
X ( s )−sY ( s )
Y ( s )=
5
b) y ' ' ( t )+ 6 y ' ( t )+ 4 y ( t ) =x ( t ) . Con x ( t ) =u ( t )

'' '
I. h ( t ) +6 h ( t ) + 4 h ( t )=δ ( t )
s2 H ( s ) +6 sH ( s ) +4 H ( s ) =1w2
H ( s ) ( s + 6 s+ 4 )=1
2

1 1
H ( s )= 2 =
s + 6 s +4 ( s+ 3 )2−5

h ( t )=L
−1
{ }
1
( s +3 )2−5

¿ e−3 t ∙ L−1
{ } 2
1
s −5

¿e
−3 t
∙L
{ (√ ) }
−1
2
s− 5
1
2

−3 t
e
h ( t )= ∙ sinh ( √ 5 t )
√5
II. En Matlab

III. y ( t ) =x ( t )∗h ( t )
t −3 τ
e
y ( t ) =∫ ∙ sinh ( √5 τ ) dτ
0 √5
t
1
∙ ∫ e ∙ sinh ( √ 5 τ ) dτ
−3 τ
y (t)=
√ 0
5

¿−
1
3 √5
∙∫ sinh
3 (
−√ 5 ∙ ln ( u )
du ) u=e−3 τ
−3 τ
du=−3 e dτ
du
e−3 τ dτ =
3
−ln ( u )
τ=
3
2 √5
t 3
1 −u +1
¿− ∙∫ du
3 √5 0
√5
2u 3

2 √5
t 3
1 −u +1
¿− ∙∫ du
6 √5 0
√5
3
u

( )
√5
1 1
¿− ∙ −∫ u 3 du+∫ du
6 √5 √5
3
u

( )
√5 +3 − √ 5+3
3 3
1 −3 u 3u
¿− ∙+
6 √5 √ 5+3 −√ 5+3

( )
√ 5+3 − √ 5+3

1 +3 ( √5−3 ) u 3 3 ( √ 5+3 ) u 3
¿− ∙ +
6 √5 4 4
3
¿ [ ( √ 5−3 ) ∙ e−3 t ∙ e−(√ 5+3 )t + ( √ 5+3 ) e−3 t ∙ e−(−√ 5+3 ) t ]
24 √ 5
e−3 t
¿ [ ( √ 5−3 ) ∙ e−(√ 5+3 ) t +( √ 5+3 ) ∙ e−(−√ 5+3 ) t ]
8 √5
−3 t
e
y (t)= [ 3 sinh ( √ 5 t ) + √ 5 cosh ( √5 t ) ]
4√5
V. y ' ' ( t )+ 6 y ' ( t )+ 4 y ( t ) =x ( t )
s2 Y ( s ) +6 sY ( s ) +4 Y ( s ) =X ( s )
2
4 Y ( s )=X ( s )−s Y ( s )−6 sY ( s )
X ( s )−s2 Y ( s )−6 sY ( s )
(
Y s= )
4

c) 4 y ' ( t )+ 9 y ( t )=x ' ( t ) +2 x ( t ) . Con x (t )=e−|t|

I. 4 sY ( s ) +9 Y ( s )=sX ( s ) +2 X ( s )
Y ( s )( 4 s+9 )=X ( s ) ( s+ 2 )
( s+ 2 )
Y ( s )=X ( s ) ∙ → X ( s )=1
4 s +9
( s+2 )
Y ( s )= =H ( s )
4 s+ 9
1 1
H ( s )= −
4 4 ( 4 s +9 )

h ( t )=
δ ( t ) L−1
4
− {
4 4 s+ 9
−9 t
1
}
δ (t ) 1 1 4
¿ − ∙ e
4 4 4

−9 t
δ (t ) 1 4
h ( t )= − e
4 16

II. En Matlab
III. y ( t ) =x ( t )∗h ( t )
t t −9 ( t −τ )
δ ( t−τ ) x (τ )
y ( t ) =∫ x ( τ ) ∙ dτ −¿∫ ∙e 4
dτ ¿
0
4 0
16
−9 t
4
x ( t )∗δ ( t ) x ( t )∗e
¿ −
4 16
−9 t
1
¿ x ( t )− x ( t )∗e 4
16
−9 t −9 t
4 −32 4
x ( t )∗e = e
65
−9t
e−t 2 4
y (t)= + e
4 65

' '
VI. 4 y ( t )+ 9 y ( t )=x ( t ) +2 x ( t )
4 sY ( s ) +9 Y ( s )=sX ( s ) +2 X ( s )
sX ( s ) +2 X ( s )−4 sY ( s )
Y ( s )=
9
En este caso podemos observar como las gráficas coinciden y presentan
respuestas similares en ambos casos, una respuesta exponencial al sistema

También podría gustarte