怎样改正这个错误

deven1985

程序如下，错误处已用红色标出，请问各位怎样改正，谢谢！！
! b+ {, i$ n0 F* `$ l5 B1 @! I
  L& U# Z$ R* Y& g! a9 ]4 ~7 Vclear all;
) h1 i! B! x3 R9 Lclose all;
  V7 R/ a$ d, u) E2 ]6 x; R  [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];
4 x1 d+ ~) j0 P& K. XC1=[1,0,0,0];
D1=[0,0];
[A,B,C,D]=c2dm(A1,B1,C1,D1,ts,'z');
x=[0;0;0;0];
. J4 y. p# c7 {1 i+ `* ]/ `8 Er_1=0;r_2=0;g_1=0;g_2=0;
c1=5.0;c2=5.0;
eq=5;
) b9 ^4 k" h3 u' ?( a& \q=30;
Ce=[c1,0,1,0;0,c2,0,1];
, e1 ~) \$ g( c. f7 Ofor k=1:1:2000
) S1 ~$ V) r0 Q1 K; m' K: Btime(k)=k*ts;
   r(k)=1.0;
; N% Z- L$ r2 e7 h6 A5 Q   g(k)=1.0;
   
   %Using Waitui method   
   dr(k)=(r(k)-r_1)/ts;
" k) K5 |- C/ ]   dg(k)=(g(k)-g_1)/ts;
   dr_1=(r_1-r_2)/ts;
- h% S$ f# x0 u# Q   dg_1=(g_1-g_2)/ts;
   r1(k)=2*r(k)-r_1;
   g1(k)=2*g(k)-g_1;
   dr1(k)=2*dr(k)-dr_1;
   dg1(k)=2*dg(k)-dg_1;
# E$ ]( w' O7 O  
3 ]" d7 w1 C2 x   R=[r(k);g(k);dr(k);dg(k)];
+ Z" p# Z- H  m% ^1 X' K1 q) W9 f   R1=[r1(k);g1(k);dr1(k);dg1(k)];
9 m3 w) V* @4 u1 m   
# t$ S, h* T1 W- M   E=R-x;
   e(k)=E(1);
2 P6 E7 q( o2 I8 u: x. Q$ o5 C- n8 G   de(k)=E(2);
2 A% t+ x0 Q9 r; e3 M   
3 ]# }! k1 W, _2 _) i/ l9 e4 z! B   s(k)=Ce*E;   
, D; N1 Y& f4 H1 ^0 ^   X1=abs(e(k))+abs(de(k));
   
0 x, y2 B% M6 i- R6 SM=3;
if M==1             %EXP reaching law
$ p7 R- W3 [& e! Y& g- n8 Q* V5 i   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
* N5 X+ I- T8 x! h" |: S   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
& c- g8 X; Y8 v+ U$ y   k0=0.60;
1 O" a- j, z2 G6 \   if X1>k0          %EXP reachine law
& \, _( b1 `; E) D    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));
' g; T1 ~# Y$ F' T. p# z   elseif X1<=k0     %Variable rate reachine law
   ds(k)=-eq*ts*X1*sign(s(k));
   u(k)=inv(Ce*B)*(Ce*R1-Ce*A*x-s(k)-ds(k));
   end
end
if u(k)>=10
   u(k)=10;
end
( V  }/ }. }/ I# @* t! `if u(k)<=-10
+ K7 f4 J% g! S: a   u(k)=-10;
end
x=A*x+B*u(k);
& G2 X  s& E, G- jy(k)=x(1);
; B1 N4 o3 u: n+ ^& Z# f0 T6 }8 U+ ~%Update Parameters
! B' ^. y, f* w& d0 qr_2=r_1;
" O/ ]: |# a6 h+ Pr_1=r(k);
1 Q  T5 N$ Z# ~: ~: }4 H" e% G2 yend
1 [& A% b) @9 C; a: ~9 Mfigure(1)
plot(time,r,'r',time,y,'b');
  E. B0 C' ]. j3 W* Fxlabel('Time(second)');ylabel('Position tracking');
figure(2)
plot(time,s,'r');
$ g7 O* d/ B2 hxlabel('Time(second)');ylabel('Switch function s');
* I* T+ d0 v0 q& b* c% _figure(3)
& `' ]2 X8 {- E0 Y8 G+ f" Q& splot(e,de,'r',e,-c*e,'b');
8 s5 Z# `, a* x! }$ H5 U' Kaxis([0,0.001,-0.01,0]);
xlabel('e');ylabel('de');
$ Y' A4 z/ x- Q( k" K( }  nfigure(4)
* I1 D1 D# N+ p0 J  t0 Q- fplot(time,u,'r');
! T7 o3 A. D# f: M" b! dxlabel('Time(second)');ylabel('u');
& p' \/ J+ Q% }
???  In an assignment  A(I) = B, the number of elements in B and
- U: a8 \; T  A3 ^, u. O, jI must be the same.
Error in ==> matlab02 at 40
! q+ I3 X& ]2 `8 J" J   s(k)=Ce*E;

deven1985

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

厚积薄发

维数不匹配

都天锴

不知道。。。。