dinamica

Ingeniería Mecatrónica.
Dinámica de Sistemas

Practica 4

Creación de variables simbólicas
Operaciones con variables Simbólicas
Ejemplo
>> dsolve('D4y +16*D3y +88*D2y +192*Dy+144*y=exp(3*t)','y(0)=1','Dy(0)=0','D2y(0)=0','D3y(0)=0')

ans =
1/2025*exp(3*t)+3/800*exp(-2*t)+2581/2592*exp(-6*t)+359/80*exp(2*t)*t+215/144*exp(-6*t)*t
>> a= 'manzana'
a=
manzana
>> f='x*sin(x)'
f=x*sin(x)
>> Hj= 'ty87bsx76'
Hj =
ty87bsx76
>> syms a b c d x y
>> a= x+1;
>> b= x-7;
>> c= x+8;
>> d= x-9;
>> y= a*b*c*d
y=
(x+1)*(x-7)*(x+8)*(x-9)

1

Ingeniería Mecatrónica.
Grupo 9051Dinámica de Sistemas

>> expand (y)
ans =
x^4-7*x^3-73*x^2+439*x+504
>> p= [1 -7 -73 439 504]
p=
1 -7 -73 439 504
>> q= roots (p) % se debe cambiar el signo
q=
-8.0000
9.0000
7.0000-1.0000
>> syms x
>> p1=(x^2-4)/(x^2-1)
p1 =
(x^2-4)/(x^2-1)
>> p2=x-2;
>> p3= x-1;
>> p4 = simplify(p1*p3/p2)%FACTORIZA simplifica la function
p4 =
(x+2)/(x+1)
>> expand(p2*p3)% Realizaoperaciones
ans =
x^2-3*x+2

2

Ingeniería Mecatrónica.
Grupo 9051
Dinámica de Sistemas

>> expand(p3^2)
ans =
x^2-2*x+1
>> p=[1 16 64]
p=
1 16 64
>> q= roots(p)
q=
-8
-8
>> syms C1 C2 ty
>> y= C1*exp(-8*t)+ C2*exp(-8*t)
y=
C1*exp(-8*t)+C2*exp(-8*t)

>> z= diff(y,2,t)+ 16*diff(y,1,t)+ 64*y
z=
0

3

Ingeniería Mecatrónica.
Grupo 9051
Dinámica de Sistemas

Ejercicio enClase
>> A= [0 1 0 ; 0 -2 -5 ; 0 1 2]
A=
0

1

0

0 -2 -5
0

1

2

>> syms x
>> y= det(A-x*eye(3))
y=
- x^3 - x
>> z= solve (y)
z=
0
i
-i
>> B=[1 0 0; 1 i i^2 ; 1 -i i^2]
B=1.0000

0

0

1.0000

0 + 1.0000i -1.0000

1.0000

0 - 1.0000i -1.0000

4

Ingeniería Mecatrónica.
Grupo 9051
Dinámica de Sistemas

>> syms t
>> b=[1;cos(t)+i*sin(t);cos(t)-i*sin(t)]
b=
1
cos(t) + sin(t)*i
cos(t) - sin(t)*i
>> d= inv (B)*b
d=
1
sin(t)
1 - cos(t)
>> m= d(1,1)*eye(3)+ d(2,1)*A+d(3,1)*A^2
m=
[ 1, 2*cos(t) + sin(t) - 2,

5*cos(t) - 5]

[ 0,...