Tartalomjegyzék

Script d2017_09_19_gyak2

TELJES MATLAB SCRIPT KIEGÉSZÍTŐ FÜGGVÉNYEKKEL

file:   d2017_09_19_gyak2.m
author: Peter Polcz <ppolcz@gmail.com>
Created on 2017. September 19.
Output:
┌d2017_09_19_gyak2
│   - Persistence for `d2017_09_19_gyak2` initialized [run ID: 5806, 2017.09.19. Tuesday, 14:53:52]
│   - Script `d2017_09_19_gyak2` backuped

Állapottranszformációk, diagonális realizáció

n = 5;
A = randn(n); A = A + A';
B = randn(n,1);
C = randn(1,n);
D = randn(1,1);

[S,A_] = eig(A);

% A_ = S\A*S;
B_ = S\B;
C_ = C*S;

Original = tf(ss(A,B,C,D))
Transformed = tf(ss(A_,B_,C_,D))
Output:
Original =
 
  0.8797 s^5 - 4.293 s^4 - 9.203 s^3 + 82.28 s^2 - 153.9 s + 87.12
  ----------------------------------------------------------------
     s^5 - 3.136 s^4 - 20.67 s^3 + 77.89 s^2 - 9.547 s - 95.44
 
Continuous-time transfer function.


Transformed =
 
  0.8797 s^5 - 4.293 s^4 - 9.203 s^3 + 82.28 s^2 - 153.9 s + 87.12
  ----------------------------------------------------------------
     s^5 - 3.136 s^4 - 20.67 s^3 + 77.89 s^2 - 9.547 s - 95.44
 
Continuous-time transfer function.

End of the script.

Output:
└ 0.050309 [sec]