the problem are in : Error using horzcat

Dimensions of matrices being concatenated are not consistent.

Error in Peng_2009_Originale (line 517)

lmiterm([newl 1 1 X],A,1,'s');

And the code are:

function tmint = fes(Xm,R1m,M1m,R2m,M2m,R3m,M3m,S1m,N1m,S2m,N2m,Pm,Q1m,Q2m,Q3m,Fm)

newl=0;

A=[0 1;

0 0.1];

B=[0;0.1];

F=[-3.75 -11.5];

tau1=0.2;

tau2=0.02;

tau3=0.5;

segma=0;

lemda=segma*(1-segma);

setlmis([]);

X=Xm;

X=1

R1=R1m ;

M1=M1m;

R2=R2m ;

M2=M2m;

R3=R3m;

M3=M3m ;

S1 =S1m;

N1=N1m ;

S2 =S2m ;

N2=N2m ;

Q1=Q1m;

Q2=Q2m ;

Q3=Q3m;

F=Fm;

newl=0;

newl = newlmi;

lmiterm([newl 1 1 X],A,1,'s');

lmiterm([newl 1 1 Q1],1,1);

lmiterm([newl 1 1 Q2],1,1);

lmiterm([newl 1 1 Q3],1,1);

lmiterm([newl 1 1 R1],-1,1);

lmiterm([newl 1 1 R2],-1,1);

lmiterm([newl 1 1 R3],-1,1);

lmiterm([newl 1 2 R1],1,1);

lmiterm([newl 1 3 R2],1,1);

lmiterm([newl 1 3 F],segma*B,1);

lmiterm([newl 1 4 0],0);

lmiterm([newl 1 5 R3],1,1);

lmiterm([newl 1 5 F],(1-segma)*B,1);

lmiterm([newl 1 6 0],0);

lmiterm([newl 1 9 X],tau1,A');

lmiterm([newl 1 10 X],tau2,A');

lmiterm([newl 1 11 X],tau3,A');

lmiterm([newl 1 12 X],(tau2-tau1),A');

lmiterm([newl 1 13 X],(tau3-tau2),A');

lmiterm([newl 1 14 0],0);

lmiterm([newl 1 15 0],0);

lmiterm([newl 1 16 0],0);

lmiterm([newl 1 17 0],0);

lmiterm([newl 1 18 0],0);

%%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 2 2 R1],-1,1);

lmiterm([newl 2 2 S1],-1,1);

lmiterm([newl 2 2 Q1],-1,1);

lmiterm([newl 2 3 S1],1,1);

lmiterm([newl 2 4 0],0);

lmiterm([newl 2 5 0],0);

lmiterm([newl 2 6 0],0);

lmiterm([newl 2 9 0],0);

lmiterm([newl 2 10 0],0);

lmiterm([newl 2 11 0],0);

lmiterm([newl 2 12 0],0);

lmiterm([newl 2 13 0],0);

lmiterm([newl 2 14 0],0);

lmiterm([newl 2 15 0],0);

lmiterm([newl 2 16 0],0);

lmiterm([newl 2 17 0],0);

lmiterm([newl 2 18 0],0);

%%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 3 3 R2],-2,1);

lmiterm([newl 3 3 S1],-2,1);

lmiterm([newl 3 4 R2],1,1);

lmiterm([newl 3 4 S1],1,1);

lmiterm([newl 3 5 0],0');

lmiterm([newl 3 6 0],0);

lmiterm([newl 3 9 -F],segma*tau1,B');

lmiterm([newl 3 10 -F],segma*tau2,B');

lmiterm([newl 3 11 -F],segma*tau3,B');

lmiterm([newl 3 12 -F],segma*(tau2-tau1),B');

lmiterm([newl 3 13 -F],segma*(tau3-tau2),B');

lmiterm([newl 3 14 -F],lemda*tau1,B');

lmiterm([newl 3 15 -F],lemda*tau2,B');

lmiterm([newl 3 16 -F],lemda*tau3,B');

lmiterm([newl 3 17 -F],lemda*(tau2-tau1),B');

lmiterm([newl 3 18 -F],lemda*(tau3-tau2),B');

%%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 4 4 Q2],-1,1);

lmiterm([newl 4 4 R2],-1,1);

lmiterm([newl 4 4 S1],-1,1);

lmiterm([newl 4 4 S2],-1,1);

lmiterm([newl 4 5 S2],1,1);

lmiterm([newl 4 6 0],0);

lmiterm([newl 4 7 0],0);

lmiterm([newl 4 8 0],0);

lmiterm([newl 4 9 0],0);

lmiterm([newl 4 10 0],0);

lmiterm([newl 4 11 0],0);

lmiterm([newl 4 12 0],0);

lmiterm([newl 4 13 0],0);

lmiterm([newl 4 14 0],0);

lmiterm([newl 4 15 0],0);

lmiterm([newl 4 16 0],0);

%%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 5 5 R3],-2,1);

lmiterm([newl 5 5 S2],-2,1);

lmiterm([newl 5 6 R3],1,1);

lmiterm([newl 5 6 S2],1,1);

lmiterm([newl 5 7 -F],(1-segma)*tau1,B');

lmiterm([newl 5 8 -F],(1-segma)*tau2,B');

lmiterm([newl 5 9 -F],(1-segma)*tau3,B');

lmiterm([newl 5 10 -F],(1-segma)*(tau2-tau1),B');

lmiterm([newl 5 11 -F],(1-segma)*(tau3-tau2),B');

lmiterm([newl 5 12 -F],-lemda*tau1,B');

lmiterm([newl 5 13 -F],-lemda*tau2,B');

lmiterm([newl 5 14 -F],-lemda*tau3,B');

lmiterm([newl 5 15 -F],-lemda*(tau2-tau1),B');

lmiterm([newl 5 16 -F],-lemda*(tau3-tau2),B');

%%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 6 6 R3],-1,1);

lmiterm([newl 6 6 S2],-1,1);

lmiterm([newl 6 6 Q2],-1,1);

lmiterm([newl 6 7 0],0);

lmiterm([newl 6 8 0],0);

lmiterm([newl 6 9 0],0);

lmiterm([newl 6 10 0],0);

lmiterm([newl 6 11 0],0);

lmiterm([newl 6 12 0],0);

lmiterm([newl 6 13 0],0);

lmiterm([newl 6 14 0],0);

lmiterm([newl 6 15 0],0);

lmiterm([newl 6 16 0],0);

%%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 7 7 M1],-1,1);

lmiterm([newl 7 8 0],0);

lmiterm([newl 7 9 0],0);

lmiterm([newl 7 10 0],0);

lmiterm([newl 7 11 0],0);

lmiterm([newl 7 12 0],0);

lmiterm([newl 7 13 0],0);

lmiterm([newl 7 14 0],0);

lmiterm([newl 7 15 0],0);

lmiterm([newl 7 16 0],0);

%%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 8 8 M2],-1,1);

lmiterm([newl 8 9 0],0);

lmiterm([newl 8 10 0],0);

lmiterm([newl 8 11 0],0);

lmiterm([newl 8 12 0],0);

lmiterm([newl 8 13 0],0);

lmiterm([newl 8 14 0],0);

lmiterm([newl 8 15 0],0);

lmiterm([newl 8 16 0],0);

%%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 9 9 M3],-1,1);

lmiterm([newl 9 10 0],0);

lmiterm([newl 9 11 0],0);

lmiterm([newl 9 12 0],0);

lmiterm([newl 9 13 0],0);

lmiterm([newl 9 14 0],0);

lmiterm([newl 9 15 0],0);

lmiterm([newl 9 16 0],0);

%%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 10 10 N1],-1,1);

lmiterm([newl 10 11 0],0);

lmiterm([newl 10 12 0],0);

lmiterm([newl 10 13 0],0);

lmiterm([newl 10 14 0],0);

lmiterm([newl 10 15 0],0);

lmiterm([newl 10 16 0],0);

%%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 11 11 N2],-1,1);

lmiterm([newl 11 12 0],0);

lmiterm([newl 11 13 0],0);

lmiterm([newl 11 14 0],0);

lmiterm([newl 11 15 0],0);

lmiterm([newl 11 16 0],0);

%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 12 12 M1],-lemda,1);

lmiterm([newl 12 13 0],0);

lmiterm([newl 12 14 0],0);

lmiterm([newl 12 15 0],0);

lmiterm([newl 12 16 0],0);

%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 13 13 M2],-lemda,1);

lmiterm([newl 13 14 0],0);

lmiterm([newl 13 15 0],0);

lmiterm([newl 13 186 0],0);

%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 14 14 M3],-lemda,1);

lmiterm([newl 14 15 0],0);

lmiterm([newl 14 16 0],0);

%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 15 15 S1],-lemda,1);

lmiterm([newl 15 16 0],0);

%%%%%%%%%%%%%%%%%%%%%%%% LMI

lmiterm([newl 16 16 N2],-lemda,1);

lmisys=getlmis;

[tmin,xfeas] = feasp(lmisys,[0,0,0,0,1],0);

tmint=tmin;

end