Code covered by the BSD License  

Highlights from
Rsyso

image thumbnail

Rsyso

by

 

Finds realizations of SISO linear systems. Autonoma of Colombia University.

[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];D=0;
    case 2
        A=[-1 0 -1;0 0 1;1 -1 -1];B=[1 0 0]';C=[1 -1 -1];D=0;
    case 3
        A=[-1 0 -1;0 0 1;1 -1 -1];B=[1;0;0];C=[1 -1 -1];D=0;   
    case 4
        A=[-4 1 0;-5 0 1;-2 0 0];B=[1;4;4];C=[1 0 0];D=0;   
    case 5
        A=[-1 -2 -2;0 -1 1;1 0 -1];B=[2;0;1];C=[1 1 0];D=0; 
    case 6
        A=[-1 1 0;0 -1 0;0 0 -2];B=[0;1;1];C=[1 0 1];D=0;  
    case 7
        A=[1 0 0;1 0 0;0 1 1];B=[0;1;1];C=[0 1 -1];D=0;
    case 8
        A=[-1 1 1;0 -1 0;0 1 -1];B=[2;-1;1];C=[1 1 1];D=0; 
    case 9
        A=[-2 0 0;0 -2 0;2 1 -2];B=[0;1;1];C=[1 1 1];D=0;  
    case 10
        A=[-1 0 0;0 -1 0;0 0 -2];B=[1;0;0];C=[2 0 3];D=0;  
    case 11
        A=[0 1 -1;1 2 0;0 2 0];B=[1;1;-1];C=[0 -1 1];D=-1;
    case 12
        A=[-1 1 0;0 -1 0;0 0 -2];B=[0;1;1];C=[1 0 1];D=0;   
    case 13
        A=[-1 1 0;0 -1 0;0 0 -2];B=[0;1;1];C=[1 0 1];D=0;
    case 14
        A=[-1 0 -1;0 0 1;1 -1 -1];B=[1;0;0];C=[1 -1 -1];D=0; 
    case 15
        A=[-1 0 -1;0 0 1;1 -1 -1];B=[1;1;-1];C=[1 -1 -1];D=0;
    case 16
        A=[-1/6 0 -1/3;0 0 1;1/2 -1/2 -1/2];B=[1/6;0;0];C=[1 -1 -1];D=0.5;
    case 17
        A=[1 1 0;0 1 0;0 1 1];B=[1;1;0];C=[1 1 1];D=0;
    case 18
        A=[1 0;2 1];B=[1 1]';C=[1 1];D=0;                           
    case 19
        A=[1 0;2 1];B=[1 1]';C=[1 1 0];D=0;                           
    case 20
        A=[1 0;2 1];B=[1;1];C=[1 1];D=0;                                
    case 21
        A=[1 0;2 1];B=[1 0]';C=[1 1];D=0;                            
    case 22
        A=[-10 0;0 -1];B=[1;1];C=[10 1.1];D=1;             
    case 23
        A=[1 0;2 1];B=[1;1];C=[1;0;0];D=0;                        
    case 24
        A=[1 0;2 1];B=[1 1 0]';C=[1 1];D=0;                            
    case 25
        A=[1 0;2 1];B=[1 0 0]';C=[1 1];D=0;                        
    case 26
        A=[1 0;2 1];B=[1 1 0]';C=[1 1];D=0;                    
    case 27
        A=[1 1 0 0;0 1 1 0;0 0 1 0;0 0 0 2];B=[0;2;1;0];C=[1 0 0 0];D=0;
    case 28
        A=[1 1 0 0;0 1 1 0;0 0 1 0;0 0 0 2];B=[0;2;1;0];C=[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;1;0;0];C=[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;0;0];C=[2 0 3 1 0];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;0;0;1];
        C=[1 1 0 0 1 0 1];D=0;
    case 32
        A=[-8 -15 -1;1 0 0;0 0 -4];B=[1.5;0;0];C=[4/3 20/3 -4/3];D=0;
    case 33
        A = [-0.5572 -0.7814;0.7814 0];B = [1;0];C = [1.9691  6.4493];D=0;
    case 34
        A=[-8/3 -1 -5/3;-2/3 -2 -2/3;-4/3 1 -7/3];B=[-1;-1;-1];C=[2 -1 1];D=0;
    case 35
        A=[0 -4 3 -4;2 -9 -2 -8;0 0 1 0;-2 7 5 6];B=[-1;-1;1;2];C=[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;4;1;0];C=[0 0 -2 10];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;1;0;-6;-13];
        C=[1 0 0 0 0 0];D=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;1;1;0;0;0;1;1;0;0];
        C=[1 0 1 1 2 0 5 4 1 2];
        D=0;
    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;1;1;0;0;0;1;-1;-1;1];
        C=[1 1 1 1 2 0 5 4 1 2];
       D=0;
    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;1;0;0;1;0;1];
        C=[1 1 0 0 1 0 1];D=0;
    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 46
        A=[3 1 0;0 0 -1;0 1 -1];B=[0;1;0];C=[2 -1 1];D=0;
    case 47
        A=[0 1 0;0 0 1;-6 -11 -6];B=[1;1;0];C=[0 1 1];D=0;
    case 'z'
        A=47;
end

Contact us