Automatización de archivos

Ronny campoverde

Ejercicio 2.8.5[pic 1]

>>  num = [1 -4 20]; den = conv([1 2 10],[1 5]);

>> printsys(num,den)

num/den =

        s^2 - 4 s + 20

   -----------------------

   s^3 + 7 s^2 + 20 s + 50

>> numraices=roots(num)

numraices =

   2.0000 + 4.0000i[pic 2]

   2.0000 - 4.0000i

>> denraices=roots(den)

denraices =

  -5.0000 + 0.0000i

  -1.0000 + 3.0000i

  -1.0000 - 3.0000i

>>  pzmap(num,den),

>> % Personalización de coordenadas

>>  axis([-6 3 -5 5]) % ([Xmin Xmax Ymin Ymax])

Ejercicio 2.8.6a

[pic 3]

>>  num = [5 -15 -11]; den = [1 -5 6 4 -8]; [r,p,k] = residue(num,den)

r =

    0.3333

    4.0000

   -7.0000

   -0.3333

p =

    2.0000

    2.0000

    2.0000

   -1.0000

k =

     []

>> % que corresponde a:

>> % G(s) = 0.3333/(s−2) + 4/(s−2)2 − 7/(s−2)3 − 0.3333/(s+1)

>>

Ejercicio2.8.6b

>>  r = [-2.3750 3 -0.6250];[pic 4]

>> p = [-10 -5 -2];

>>  k = [0];

>>  [num,den] = residue(r,p,k)

num =

     0    10     5

den =

     1    17    80   100

>> printsys(num,den)

num/den =

            10 s + 5

   -------------------------

   s^3 + 17 s^2 + 80 s + 100

>> raicesden = roots(den)

raicesden =

  -10.0000

   -5.0000

   -2.0000

>>

Ejercicio 2.8.7a

>> num = [4 0 5 18 0];

>>  den = [3 12 15 24 10];

>>  inf = 1000000;

>>  valorinicial = polyval(num,inf)/polyval(den,inf)

valorinicial =

    1.3333

Ejercicio 2.8.7b

>> num = [10 5 8];

>>  den = [1 4 16 4];

>> valorfinal = polyval(num,0)/polyval(den,0) [pic 5]

valorfinal =

     2

>> densistema = [1 4 16 4];

>> polos=roots(densistema)

polos =

  -1.8667 + 3.3941i

  -1.8667 - 3.3941i

  -0.2666 + 0.0000i

Ejercicio 2.8.8b

>> syms s t[pic 6]

>> Y=(6*s-4)/(s^2+4*s+20); y=ilaplace(Y)

y =

...