Mult&T: [A B C D]=state(op) - File Exchange

No BSD License

Highlights from
Mult&T

creditos()
multitool(varargin) MULTITOOL MULTI Design Tool
questimp(varargin)
[A B C D T]=hoform(Gt,W) HOFORM Finds the Ho canonical form of a LTI MIMO SYS model.
[A B C D]=gauthierform(P,Q) GAUTHIERFORM Gauthier minimal algoritm of a MPF model.
[A B C D]=gilbertform(Gt) GILBERTFORM Finds the Gilbert form of a LTI MIMO SYS model.
[A B C D]=lcfform(Gt) LCFFORM Finds the minimal realization by method fraction left coprime.
[A B C D]=lmpf2ss(D1,N) LMPF2SS Finds the realization by method fraction coprime of a LTI
[A B C D]=minhoform(sys) MINHOFORM Finds the Minimal realization by method Ho.
[A B C D]=rcfform(Gt) RCFFORM Finds the minimal realization by method right fraction coprime.
[A B C D]=rmpf2ss(D1,N) RMPF2SS Find the realization by method fraction coprime of a MPF.
[A B C D]=silvermanform(Gt) SILVERMANFORM Finds the minimal realization by method Silverman.
[A B C D]=state(op) STATE state 45 examples in space-state.
[A B C D]=wolovichform(Gt) WOLOVICHFORM Realization by method Wolovich of a LTI MIMO SYS model.
[A B C P]=cx2rJ(A,B,C) CX2RJ Complex to real Jordan format.
[Aj,Bj,Cj,Dj]=exjordanform(Gt... EXJORDANFORM Finds the extended Jordan form of a LTI MIMO SYS model.
[Aj,Bj,Cj,Dj]=jordanform(Gt) JORDANFORM Finds the Jordan form of a LTI MIMO SYS model.
[D N]=mtf2lcf(G) MTF2LCF Matrix Transfer Function to Left Coprime Fraction.
[D N]=mtf2lmpf(Gt) MTF2LMPF Matrix Transfer Function to Left Matrix Fraction Description.
[D N]=mtf2rcf(G) MTF2RCF Matrix Transfer Function to Right Coprime Fraction.
[D N]=mtf2rmpf(Gt) MTF2RMPF Matrix Transfer Function to Right Matrix Polynomial Fraction.
[Dhc Dlc]=uplowM(D) UPLOWM Find the matrices with the hight and low coeficients.
[G]=ss2mtf(A,B,C,D) SS2MTF Convert state-space filter parameters to Matrix Transfer Function.
[Gsp Ginf]=mtfsp(G) MTFSP Strictly Proper Matrix Transfer Function.
[Gt]=lcf2mtf(D,N,Ginf) LCF2MTF Left Coprime Fraction to Matrix Transfer Function.
[Gt]=mtf(op)
[Gt]=rcf2mtf(D,N,Ginf) RCF2MTF Right Coprime Fraction to Matrix Transfer Function.
[M]=normcell(M,op) NORMCELL Normalize cell.
[P Q]=lmpf(op) LMPF 23 Examples in Left MPF Format cell.
[P Q]=ss2lmpf(sys) SS2LMPF State-space to Left Matrix Polynomial Fraction.
[P Q]=ss2rmpf(sys) SS2RMPF State-space to Right Matrix Polynomial Fraction.
[Pm Qm]=mpfred(P,Q) MPFRED Finds the minimal reduction of MPF.
[S dep lidep]=mSilvester(D,N,... MSILVESTER Find the Matrix Silvester.
[ar br cr dr]=minjordanform(s... MINJORDANFORM Finds the Minimal realization by method irreducible Jordan.
[i pnc]=islmpfc(P,Q) ISLMPFC Determine if the left MPF is controlable.
[i pnc]=isrmpfo(P,Q) ISRMPFO Determine if the Right MPF is observable.
[polo n ni]=polindex(Aj) POLINDEX Finds the index of Polynomial characteristic and the minimal.
[rp k]=mresidue(Gt) MRESIDUE Partial fraction expansion of Matrix Transfer Function.
[v Q Af Bf Cf]=canonform(A,B,... CANONFORM Find a observability and controlability canonical form.
[varargout]=cell2sym(varargin... CELL2SYM Cell format to sym format in operator differential.
[varargout]=cellround(varargi... CELLROUND Round to nearest integer in each cell.
[varargout]=sym2cell(varargin... SYM2CELL Sym matrix polynomial format to cell format.
creditos_1(action)
datagrid_cb(action)
i=iscolred(M) ISCOLRED Determine if the MPF is column reduced.
i=isrowred(M) ISROWRED Determine is the MPF is row reduced.
k=coXm(A,B,C) COXM Check the observability and controlability for mode.
mytimer(obj, event, op)
n=mindeg(sys) MINDEG Degree minimal of LTI MIMO SYS model.
polezero(sys) POLEZERO Pole and zeros of LTI MIMO SYS model.
rga(a,b,c,d,w)
smform(a,b,c,d) SMFORM Finds the Smith-McMillan form of a LTI MIMO SYS model**.
ss2sym(a,b,c,d) SS2SYM State Space to Symbolic Transfer Function conversion**.
sym2tf(g) SYM2TF Symbolic Transfer Function to Numerical Transfer Function**.
tf2sym(G) TF2SYM Numerical Transfer Function to Symbolic Transfer Function**.
v=findv(M,op) FINDV Find the blocks in matrix Jordan or the indices of MPF.
Contents.m
report.html
View all files

from Mult&T by Franklin Pineda
Find realizations of multivariable systems. Created for Msc. students at the UANDES and UAC.

[A B C D]=state(op)

function [A B C D]=state(op)
%STATE state 45 examples in space-state. 
switch op
    case 1
       A=[-1 0 -1;0 0 1;1 -1 -1];B=[1 0 0]';C=[1 -1 -1;0 1 0];D=0;        %1I-2O
    case 2
        A=[-1 0 -1;0 0 1;1 -1 -1];B=[1 0 0]';C=[1 -1 -1;0 1 0;0 0 -1];D=0; %1I-3O
    case 3
        A=[-1 0 -1;0 0 1;1 -1 -1];B=[1 1;0 0;0 0];C=[1 -1 -1];D=0;         %2I-1O
    case 4
        A=[-4 1 0;-5 0 1;-2 0 0];B=[1 0;4 1;4 1];C=[1 0 0];D=0;             %2I-1O
    case 5
        A=[-1 -2 -2;0 -1 1;1 0 -1];B=[2 -8;0 2;1 2];C=[1 1 0];D=0;         %2I-1O
    case 6
        A=[-1 1 0;0 -1 0;0 0 -2];B=[0 0;1 0;1 1];C=[1 0 1;0 1 0];D=0;      %2I-2O
    case 7
        A=[1 0 0;1 0 0;0 1 1];B=[0 1;1 0;1 0];C=[0 1 -1;1 0 0];D=0;        %2I-2O
    case 8
        A=[-1 1 1;0 -1 0;0 1 -1];B=[2 0;-1 0;1 1];C=[1 1 1;0 1 1];D=0;     %2I-2O
    case 9
        A=[-2 0 0;0 -2 0;2 1 -2];B=[0 0;1 0;1 1];C=[1 1 1;0 0 1];D=0;      %2I-2O
    case 10
        A=[-1 0 0;0 -1 0;0 0 -2];B=[1 0;0 1;0 1];C=[2 0 3;1 1 0];D=0;       %2I-2O
    case 11
        A=[0 1 -1;1 2 0;0 2 0];B=[1 4;1 0;-1 -3];C=[0 1 0;0 -1 1];D=[0 0;-1 0];
    case 12
        A=[-1 1 0;0 -1 0;0 0 -2];B=[0 0;1 0;1 1];C=[1 0 1;1 0 1];D=0;      %2I-2O
    case 13
        A=[-1 1 0;0 -1 0;0 0 -2];B=[0 0;1 0;1 1];C=[1 0 1;1 0 1;0 0 1];D=0;%2I-3O
    case 14
        A=[-1 0 -1;0 0 1;1 -1 -1];B=[1 0 0;0 1 0;0 1 -1];C=[1 -1 -1];D=0;  %3I-1O
    case 15
        A=[-1 0 -1;0 0 1;1 -1 -1];B=[1 0 0;0 1 0;0 1 -1];C=[1 -1 -1;0 0 1];D=0;%2O
    case 16
        A=[-1/6 0 -1/3;0 0 1;1/2 -1/2 -1/2];B=[1/6 1/3;0 0;0 0];C=[1 -1 -1;-0.5 0 0];D=[0 0;0.5 0];
    case 17
        A=[1 1 0;0 1 0;0 1 1];B=[0 1;1 0;0 1];C=[1 1 1];D=0;
    case 18
        A=[1 0;2 1];B=[1 1]';C=[1 0;1 1];D=0;                               %1I-2O
    case 19
        A=[1 0;2 1];B=[1 1]';C=[1 0;1 1;0 1];D=0;                           %1I-3O
    case 20
        A=[1 0;2 1];B=[1 1;1 0];C=[1 1];D=0;                                %2I-1O
    case 21
        A=[1 0;2 1];B=[1 1;1 0];C=[1 1;1 0];D=0;                            %2I-2O
    case 22
        A=[-10 0;0 -1];B=[1 0;0 1];C=[10 1.1;10 0];D=[0 0;0 1];              %2I-2O
    case 23
        A=[1 0;2 1];B=[1 1;1 0];C=[1 1;1 0;0 1];D=0;                        %2I-3O
    case 24
        A=[1 0;2 1];B=[1 1 0;1 0 0];C=[1 1];D=0;                            %3I-1O
    case 25
        A=[1 0;2 1];B=[1 1 0;1 0 0];C=[1 1;1 0];D=0;                        %3I-2O
    case 26
        A=[1 0;2 1];B=[1 1 0;1 0 0];C=[1 1;1 0;0 1];D=0;                    %3I-3O
    case 27
        A=[1 1 0 0;0 1 1 0;0 0 1 0;0 0 0 2];B=[0 0 0;0 0 2;1 0 0;0 1 0];C=[1 0 0 0;0 0 0 1];D=0;
    case 28
        A=[1 1 0 0;0 1 1 0;0 0 1 0;0 0 0 2];B=[0 0 0;0 0 2;1 0 0;0 1 0];C=[1 0 0 0;0 0 1 -1];D=0;
    case 29
        A=[0 1 0 0;3 0 0 2;0 0 0 1;0 -2 0 0];B=[0 0;1 0;0 0;0 1];C=[1 0 0 0;0 0 1 0];D=0;
    case 30
        A=[-1 0 0 1 0;1 -1 1 1 1;0 0 -2 0 0;0 0 0 -0.1 0;0 0 0 0 -0.2];B=[1 0;0 1;0 1;0 0;0 0];C=[2 0 3 1 0;1 1 0 0 1];D=0;                %2I-2O
    case 31
        A=[-1 1 0 0 0 0 0;0 -1 0 0 0 0 0;0 0 -1 1 0 0 0;0 0 0 -1 1 0 0;0 0 0 0 -1 0 0;0 0 0 0 0 -2 1;0 0 0 0 0 0 -2];
        B=[1 0 0;0 1 0;0 0 0;0 0 0;0 1 0;0 0 0;0 0 1];
        C=[1 1 0 0 1 0 1;1 0 1 1 0 0 0;1 0 1 0 0 1 0];D=zeros(3);
    case 32
        A=[-8 -15 -1;1 0 0;0 0 -4];B=[1.5 1;0 0;0 1];C=[4/3 20/3 -4/3;0 0 2;2 6 2];D=zeros(2);
    case 33
        A = [-0.5572 -0.7814;0.7814 0];B = [1 -1;0 2];C = [1.9691  6.4493];D=zeros(1,2);
    case 34
        A=[-8/3 -1 -5/3;-2/3 -2 -2/3;-4/3 1 -7/3];B=[1 -1;0 -1;-1 1];C=[2 -1 1;1 1 1];D=zeros(2,3);
    case 35
        A=[0 -4 3 -4;2 -9 -2 -8;0 0 1 0;-2 7 5 6];B=[1 -1;-1 -1.5;1 3;2 19];C=[-1 -1 7 -2;2 1 -5 2;0 -3 11 -4];D=0;
    case 36
        A=[0 0 0 15;1 0 0 22;0 1 0 6;0 0 1 -2];B=[5 -1;4 0;1 1;0 0];C=[0 0 -2 10;0 0 -1 1];D=0;
    case 37
        A=[-6 1 0 0 0 6;-2 0 1 -1 0 0;-1 0 0 0 0 2;1 0 0 -3 1 0;0 0 0 -5 0 -1;0 0 0 1 0 -2];
        B=[0 0 42;0 1 0;1 0 14;0 0 0;1 1 -6;0 -1 -13];
        C=[1 0 0 0 0 0;0 0 0 1 0 0;0 0 0 0 0 1];D=[0 0 0;0 0 0;0 0 7];
    case 38
        A=[1 zeros(1,5);0 -1 1 0 0 0;0 0 -1 0 0 0;0 0 0 -1 1 0;0 0 0 0 -1 1;0 0 0 0 0 -1];
        B=[1 0 0 1 1 0]'; C=[1 1 0 0 0 0];D=0;
    case 39
        A=[-1 zeros(1,9);zeros(1,10);0 0 -2 1 zeros(1,6);0 0 0 -2 zeros(1,6);0 0 0 0 -3 1 0 0 0 0;...
            zeros(1,5) -3 1 0 0 0; zeros(1,6) -3 0 0 0;zeros(1,7) -3 1 0;zeros(1,8) -3 0; zeros(1,9) -3];
        B=[0 0 0;1 0 2;1 1 1;0 0 0;0 0 0;0 0 0;1 0 0;1 1 0;0 1 0;0 0 1];
        C=[1 0 1 1 2 0 5 4 1 2;2 0 0 -1 3 1 2 -7 6 7];
        D=zeros(2,3);
    case 40
        A=[-1 zeros(1,9);0 -2 zeros(1,8);0 0 -2 1 zeros(1,6);0 0 0 -2 zeros(1,6);0 0 0 0 -3 1 0 0 0 0;...
            zeros(1,5) -3 1 0 0 0; zeros(1,6) -3 0 0 0;zeros(1,7) -3 1 0;zeros(1,8) -3 0; zeros(1,9) -3];
        B=[0 0 1;1 0 2;1 1 1;0 0 0;0 0 0;0 0 0;1 2 3;-1 1 0;-1 2 1;1 1 2];
        C=[1 1 1 1 2 0 5 4 1 2;2 4 -4 -4 3 1 2 -7 6 7];
        D=zeros(2,3);
    case 41
        A=[-1 1 zeros(1,5);0 -1 0 0 0 0 0;0 0 -1 1 0 0 0;0 0 0 -1 1 0 0;0 0 0 0 -1 0 0;0 0 0 0 0 -2 1;zeros(1,6) -2];
        B=[1 0 0; 0 1 0;0 0 0;0 0 0;0 1 0;0 0 0;0 0 1];
        C=[1 1 0 0 1 0 1;1 0 1 1 0 0 0;1 0 1 0 0 0 0];
        D=zeros(3,3);
    case 42
        %complejos bloque real de 3
        A=[2 1 0 0 0 0 0;0 2 1 0 0 0 0;0 0 2 0 0 0 0;0 0 0 1+2*i 1 0 0;0 0 0 0 1+2*i 0 0;0 0 0 0 0 1-2*i 1;0 0 0 0 0 0 1-2*i];
        B=[3;1;2;2-3*j;1;2+3*j;1];
        C=[3 1 2 1 -i 1 i];
        D=0;
    case 43
        %bloque de dos reeales
        A=[2 1 0 0 0 0;0 2 0 0 0 0;0 0 1+2*i 1 0 0;0 0 0 1+2*i 0 0;0 0 0 0 1-2*i 1;0 0 0 0 0 1-2*i];
        B=[1;2;2-3*j;1;2+3*j;1];
        C=[1 2 1 -i 1 i];D=0;
    case 44
        %solo imaginarios
        A=[1+2*i 1 0 0;0 1+2*i 0 0;0 0 1-2*i 1;0 0 0 1-2*i];
        B=[2-3*j;1;2+3*j;1];
        C=[1 -i 1 i];D=0;
    case 45
        %bloque de tres distitnso
        A=[4 0 0 0 0 0 0;0 3 0 0 0 0 0;0 0 2 0 0 0 0;0 0 0 1+2*i 1 0 0;0 0 0 0 1+2*i 0 0;0 0 0 0 0 1-2*i 1;0 0 0 0 0 0 1-2*i];
        B=[3;1;2;2-3*j;1;2+3*j;1];
        C=[3 1 2 1 -i 1 i];D=0;
    case 'z'
        A=45;
end

Highlights from Mult&T

Highlights from
Mult&T