怎样改正这个错误

deven1985

程序如下，错误处已用红色标出，请问各位怎样改正，谢谢！！
" R! @( ~5 a. _3 D$ I4 n
0 L' ~! V4 [. V& k5 vclear all;
close all;
a=59.36;b=7.42;m=2374.46;n=296.92;
ts=0.001;
A1=[0,0,1,0;0,0,0,1;a,-b,0,0;b,-a,0,0];
B1=[0,0;0,0;m,-n;n,-m];
6 X& c* A. K  ~' j2 M1 N3 M' t  jC1=[1,0,0,0];
D1=[0,0];
. w& R" l$ j4 F2 J; v[A,B,C,D]=c2dm(A1,B1,C1,D1,ts,'z');
x=[0;0;0;0];
r_1=0;r_2=0;g_1=0;g_2=0;
c1=5.0;c2=5.0;
2 \3 s; Y- [' f/ b6 Req=5;
: u" v0 C; }, Q4 Y, |" }q=30;
Ce=[c1,0,1,0;0,c2,0,1];
for k=1:1:2000
% ^: C& J' L% N# [5 T3 htime(k)=k*ts;
5 L8 H2 z0 X# K; t! y   r(k)=1.0;
   g(k)=1.0;
1 y) ^1 L; N* `& X+ U0 F/ V# {   
' M* f! l) u: K( t   %Using Waitui method   
9 U+ I: e+ ?4 o9 w0 @2 i7 Y6 p: ?   dr(k)=(r(k)-r_1)/ts;
   dg(k)=(g(k)-g_1)/ts;
4 T# t+ g  K6 N* ?* a   dr_1=(r_1-r_2)/ts;
8 i! R3 y, x+ ?9 R8 l6 p6 W; j   dg_1=(g_1-g_2)/ts;
   r1(k)=2*r(k)-r_1;
2 }6 i0 f& u) p' N8 B8 R   g1(k)=2*g(k)-g_1;
   dr1(k)=2*dr(k)-dr_1;
; C, x& r& ?) q, i  h   dg1(k)=2*dg(k)-dg_1;
9 z: c6 Y5 p/ B/ X% w, O  
0 O3 t% e2 M6 x" R5 Z# _4 }   R=[r(k);g(k);dr(k);dg(k)];
   R1=[r1(k);g1(k);dr1(k);dg1(k)];
7 L8 X4 W1 B; n# t6 `' S   
   E=R-x;
   e(k)=E(1);
0 J; e! Z6 \- N$ v. k, h% G   de(k)=E(2);
8 m/ s6 [5 Z: G, G; C* v7 m* W   
   s(k)=Ce*E;   
. E$ z2 I' ?" i+ b/ u+ x+ `- v! y   X1=abs(e(k))+abs(de(k));
1 a0 L+ a7 Z) h2 f   
M=3;
$ g/ m  e- J; s1 z  Vif M==1             %EXP reaching law
! e3 Q- h$ K8 h9 ^   ds(k)=-eq*ts*sign(s(k))-q*ts*s(k);
   u(k)=inv(Ce*B)*(Ce*R1-Ce*A*x-s(k)-ds(k));
elseif M==2         %Variable rate reachine law
# s$ _. O& J" ~2 U   ds(k)=-eq*ts*X1*sign(s(k));
   u(k)=inv(Ce*B)*(Ce*R1-Ce*A*x-s(k)-ds(k));
elseif M==3         %Coposite reaching law
   k0=0.60;
   if X1>k0          %EXP reachine law
7 ?) U; E/ d1 ^1 J4 {+ X0 V& v" t    ds(k)=-eq*ts*sign(s(k))-q*ts*s(k);
6 l2 ?$ w4 h# Z0 t* h    u(k)=inv(Ce*B)*(Ce*R1-Ce*A*x-s(k)-ds(k));
" F- A+ S# r! ~$ b# ^5 A3 f   elseif X1<=k0     %Variable rate reachine law
   ds(k)=-eq*ts*X1*sign(s(k));
+ U2 n% r8 \. ~6 A% y9 C   u(k)=inv(Ce*B)*(Ce*R1-Ce*A*x-s(k)-ds(k));
/ T6 A. j* e% ]" R1 N- A7 C1 _   end
/ b. @: i3 G- V* b8 X; ?end
! M6 b8 _- E! f! g7 j8 w5 z. B, rif u(k)>=10
   u(k)=10;
end
if u(k)<=-10
2 g' l  s. x' d0 l9 M4 |0 A   u(k)=-10;
end
x=A*x+B*u(k);
y(k)=x(1);
5 K6 u! U3 e5 V* v& E( Y0 \%Update Parameters
& j6 B7 Y( s7 d9 R* S' zr_2=r_1;
5 P) O1 l. a/ V1 K- D3 nr_1=r(k);
end
+ X! \# e, {7 d1 x- }4 a) O: [5 gfigure(1)
  |+ h4 L, \+ U: D4 a& Gplot(time,r,'r',time,y,'b');
xlabel('Time(second)');ylabel('Position tracking');
0 }- E0 \: I' ]3 jfigure(2)
2 M2 Q, y, U4 G3 m& ]; }" y( Oplot(time,s,'r');
& k- }. _0 r$ Exlabel('Time(second)');ylabel('Switch function s');
figure(3)
' o. q% Y& y6 T9 [plot(e,de,'r',e,-c*e,'b');
axis([0,0.001,-0.01,0]);
xlabel('e');ylabel('de');
- w9 n$ D* y- j' \4 I7 W# P, ?figure(4)
# X. s5 g/ P" F- }+ N; d) Aplot(time,u,'r');
+ C" o+ l: Y7 |/ bxlabel('Time(second)');ylabel('u');

. [: M: Y5 |" u6 H???  In an assignment  A(I) = B, the number of elements in B and
" u  v6 D( O% `I must be the same.
Error in ==> matlab02 at 40
   s(k)=Ce*E;

deven1985

补充一下：当我把错误行中的“s(k)=Ce*E”换成“s=Ce*E”后，错误就没有了，这是为什么？谢谢大家了！

厚积薄发

维数不匹配

都天锴

不知道。。。。