Code covered by the BSD License  

Highlights from
Rsyso

image thumbnail

Rsyso

by

 

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

[G1]=etf(op)
function [G1]=etf(op)
%TF 80 TF Examples. 
switch op
    case 1
        G1=tf([1 2],[1 2 1]);
    case 2
        G1=tf(1,[1 1]);
    case 3
        G1=tf(1,[-1 1]);
    case 4
        G1=tf([1 0 1],[1 0 0 0]);
    case 5
        G1=tf([1 2],[1 2 1]);
    case 6
        G1=tf([1 1 2],[1 1 0]);
    case 7
        G1=tf([1 1 2],[1 1 0 0]);
    case 8
        G1=tf(2,[1 3 2]);
    case 9
        denp=[1 6 13 12 4];G1=tf([1 0],denp);
    case 10
        G1 = tf([1 -1],[1 1]);
    case 11
        G1 = tf([1 1],[1 0 0]);
    case 12
        G1 = tf([1 0],[1 -2 1]);
    case 13
        G1 = tf([1 0],[1 2 1]);
    case 14
        G1 = tf(10,[1 2 -3]);
    case 15
        G1=tf([1 2],[1 1]);
    case 16
        denp=1.25*conv([1 1],[1 2]);G1=tf([1 -1],denp);
    case 17
        G1=tf(1,[1 1]);
    case 18
        G1=tf(0.071,[1 0.19]);
    case 19
        G1=tf([-1 0],[1 1]);
    case 20
        G1=tf([1 1],[1 3]);
    case 21
        denp=[1 4 5];G1=tf(2,denp);
    case 22
        denp=[1 0.9 -2.2];
        G1=tf([1 1],denp);
    case 23
        denp=[1 2 1];G1=tf([1 1],denp);
    case 24
        denp=[1 2 1];G1=tf([1 1],denp);
    case 25
        G1=tf(2,[1 1]);
    case 26
        denp=[0.2 1.2 1];G1=tf(1,denp);
    case 27
        G1=tf([11 -18 -70 -50],[1 6 -45 -50 0]);
    case 28
        G1=tf([1 -1],[1 1]);
    case 29
        %G1=tf([4 -10],[2 1]);G2=tf([3],[1 2]);G3=tf([1],[2 5 2]);G4=tf([1 1],[1 4 4]);Gt=[G1 G2;G3 G4];
        G1=tf([2 -5],[1 1/2]);
    case 30
        denp=[1 2 -1 -2];G1=tf([1 1 -2],denp);
    case 31
        G1=tf([1 0],[1 1]);
    case 32
        G=tf(1,[1 1]);G2=tf(4,[1 2 2]);G3=tf([1 2],[1 3]);
        G1=series(G,G2);G1=series(G1,G3);
    case 33
        G1=tf([-12 -6],[3 24]);
    case 34
        G1 = tf(10,[1 2 10]);
    case 35
        G1=tf([1 0],[1 1]);
    case 36
        G1=tf([1 1],[1 1 1]);
    case 37
        G1=tf(2,[1 3]);
    case 38
        G1=tf([1 0 1],[1 0 0 0]);
    case 39
        G1=tf(1,[1 2 1]);
    case 40
        G1=tf([1 2],[1 1]);
    case 41
        G1=tf(1,[1 2 1]);
    case 42
        G=tf(8,[1 2]); G2=tf([1 2],[1 2 2]);G3=tf([1 1],[1 1 1]);
        G1=series(G,G2);G1=series(G1,G3);
    case 43
        G1=tf([1 0 1],[1 0 0 0]);
    case 44
        den=conv([1 1],[1 2]);
        G1=tf([-47 2],den);
    case 45
        G1=tf([1 0 0],[1 0 0 1])
    case 46
        G1=tf(1,[1 3]);
    case 47
        G1=tf(conv([4 2],[1 4 8]),conv(conv([1 3],[1 2 2]),[1 1]));
    case 48
        G1=zpk(-1,[-2 -2 -2 -3 -4],2);
    case 49
        G1=tf(1,[1 1]);
    case 50
        G1=tf(1,[1 2]);
    case 51
        G1=tf([1 2],[1 1 1]);
    case 52
        den=[1 4 3];
        G1=tf([2 4],den);
    case 53
        G1=tf(4,[1 3 2]);
    case 54
        den=[1 5 8 4];
        G1=tf(1,den);
    case 55
        G1=tf(-2,[1 -1]);
    case 56
        G1=tf([-2 1],[1 2 1]);
    case 57
        G1=tf([1 2 3],[1 2 1]);
    case 58
        G1=tf(2,[1 3 2]);
    case 59
        G1=tf(1,[1 11  42  68  40]);
    case 60
        den=[1 0 0 0 0];
        G1=tf([1 -1 0 1],den);
    case 61
        G1=tf([-1 0],[1 2 1]);
    case 62
        den=[1 -1];
        G1=tf(1,den);
    case 63
        G1=tf([1 0],[1 2 1]);
    case 64
        den=[2 4];
        G1=tf([5 4],den);
    case 65
        den=[1 3 2];
        G1=tf([2 3],den);
    case 66
        G1=tf([1 2],[1 2 1]);
    case 67
        den=conv(conv(conv([1 1],[1 1]),[1 1]),[1 2]);
        G1=tf(conv(conv(conv([1 2],[1 2]),[1 2]),[1 1]),den);
    case 68
        den=conv(conv(conv([1 1],[1 1]),[1 1]),[1 2]);
        G1=tf(conv(conv(conv([1 1],[1 1]),[1 2]),[2 3]),den);
    case 69
        den=conv(conv(conv([1 1],[1 1]),[1 1]),[1 2]);
        G1=tf(conv(conv([1 2],[1 2]),[1 2]),den);        
    case 70
        G1=tf(1,conv([1 1],[1 2]));
    case 71
        den=[1 -1 -2 -4 3];
        G1=tf([1 -1 -3 2 0],den);
    case 72
        G1=tf([-0.1 0.1],[1 6/5 1/5]);
    case 73
        G1=tf(1,[1 0 0]);
    case 74
        G1=tf([1 0],[1 2]);
    case 75
        G1=tf([1 1 0],[1 0 1]);
    case 76
        G1=tf([4 7 3],[1 4 5 2]);
    case 77 
        G1=tf(-33.89,conv([98.02 1],[0.42 1]));
    case 78
        G1=tf([1 0 1],[1 0 0 0]);
    case 79
        G1=tf(1,[1 0 0 1]);
   case 80
        G1=tf(0.1,[1 5/6 1/6]);
    case 81
        num=[1 2 3 2 3 2 3 4 5];
        den=[1 1 1 2 3 2 1 4 10];
        G1=tf(num,den);
    case 82
        num=[2 7 9 12];
        den=[6 12 24 42 70];
        G1=tf(num,den);
    case 83
        %p=[-1+1j -1-1j -3+2j -3-2j -4+1j -4-1j -1/2+3j -1/2-3j -2 -3];
        p=[-1.2713+0.4576i -1.2713-0.4576i -0.7028+1.1358i -0.7028-1.1358i  0.4842+1.3350i  0.4842-1.3350i  0.9899+0.7366i  0.9899-0.7366i];
        G1=zpk([],p,1);
    case 'z'
        G1=83;
end

Contact us