怎样改正这个错误

deven1985

程序如下，错误处已用红色标出，请问各位怎样改正，谢谢！！
4 e) w) G8 P! Y# F! ?
clear all;
close all;
% l* u* J& ?9 q- Y- ma=59.36;b=7.42;m=2374.46;n=296.92;
$ i6 H% y6 P7 F+ Wts=0.001;
A1=[0,0,1,0;0,0,0,1;a,-b,0,0;b,-a,0,0];
* ~2 E! }1 [* q7 }2 x& EB1=[0,0;0,0;m,-n;n,-m];
' M% a$ K4 P! r9 S9 F8 z: hC1=[1,0,0,0];
  s* v6 D! b  d1 \9 ]2 d% t2 c3 PD1=[0,0];
[A,B,C,D]=c2dm(A1,B1,C1,D1,ts,'z');
x=[0;0;0;0];
% a" ~6 X1 J- T; H" }4 Jr_1=0;r_2=0;g_1=0;g_2=0;
; n9 ]1 R5 G! j/ a# h, |c1=5.0;c2=5.0;
eq=5;
q=30;
# ^* f  @9 Q$ ?( R0 e* z# I. t5 w/ w2 \Ce=[c1,0,1,0;0,c2,0,1];
for k=1:1:2000
. ~' I1 F. ]' W( [time(k)=k*ts;
2 ~7 o5 {/ b4 j: _5 ~  u   r(k)=1.0;
   g(k)=1.0;
   
/ q" _8 E' ^% ^   %Using Waitui method   
& p3 Q; h% }/ N2 m6 w2 c6 z   dr(k)=(r(k)-r_1)/ts;
, H1 A/ z- \) q' o/ N   dg(k)=(g(k)-g_1)/ts;
   dr_1=(r_1-r_2)/ts;
   dg_1=(g_1-g_2)/ts;
+ b/ _8 o$ \. N   r1(k)=2*r(k)-r_1;
$ G& @6 y4 v$ T7 u4 c   g1(k)=2*g(k)-g_1;
   dr1(k)=2*dr(k)-dr_1;
' G! n4 |0 ]& H. x   dg1(k)=2*dg(k)-dg_1;
) [# W7 ]; V9 [8 m4 ^' U  
; u* V, }/ \+ L4 S   R=[r(k);g(k);dr(k);dg(k)];
. a" B' u8 Y# ?% z   R1=[r1(k);g1(k);dr1(k);dg1(k)];
2 c- P4 A- \4 Z0 S1 }4 J( K; x   
$ a0 i. d" \9 O: |5 j   E=R-x;
   e(k)=E(1);
3 `. d6 I* g+ a5 T! v7 W3 Z3 q+ D   de(k)=E(2);
% c. F. x+ H& L2 d   
   s(k)=Ce*E;   
2 e% u; _! I% F1 Y0 ^* {   X1=abs(e(k))+abs(de(k));
8 {0 x! J% l; D8 U" C* F   
* \- \$ J" G6 E0 @& E& d+ @3 FM=3;
if M==1             %EXP reaching law
   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
   ds(k)=-eq*ts*X1*sign(s(k));
   u(k)=inv(Ce*B)*(Ce*R1-Ce*A*x-s(k)-ds(k));
) i, S. ^3 H! ]! _( O0 `elseif M==3         %Coposite reaching law
0 o7 o; m6 L9 Q) f" K   k0=0.60;
   if X1>k0          %EXP reachine law
    ds(k)=-eq*ts*sign(s(k))-q*ts*s(k);
, p0 `( _6 k- ~2 S; s2 g    u(k)=inv(Ce*B)*(Ce*R1-Ce*A*x-s(k)-ds(k));
   elseif X1<=k0     %Variable rate reachine law
& ^7 \: a2 h8 C8 @$ U   ds(k)=-eq*ts*X1*sign(s(k));
   u(k)=inv(Ce*B)*(Ce*R1-Ce*A*x-s(k)-ds(k));
3 K  G4 ?* Q: u' m2 ?   end
end
  y( @' E; L4 I  ~! @) h8 hif u(k)>=10
3 \0 J/ e; [# s; I   u(k)=10;
end
if u(k)<=-10
0 t1 |/ x- c- h0 I   u(k)=-10;
: `. \) ]0 A% yend
x=A*x+B*u(k);
y(k)=x(1);
. ]% K1 n" D  L%Update Parameters
r_2=r_1;
r_1=r(k);
end
( a% @# }( D, V6 ?% a+ |/ d0 Qfigure(1)
plot(time,r,'r',time,y,'b');
xlabel('Time(second)');ylabel('Position tracking');
* G! T" E. C# D, \1 O" {figure(2)
1 d! ^& {: c2 vplot(time,s,'r');
xlabel('Time(second)');ylabel('Switch function s');
9 E. e$ a: V  ^# n* I4 \' ffigure(3)
6 P9 a0 E; f( S: T3 `8 G# Gplot(e,de,'r',e,-c*e,'b');
axis([0,0.001,-0.01,0]);
, g# W" b( E, L0 N: l/ gxlabel('e');ylabel('de');
figure(4)
plot(time,u,'r');
xlabel('Time(second)');ylabel('u');

0 h# r" _& L! f???  In an assignment  A(I) = B, the number of elements in B and
0 i4 X! u5 R: w  t6 y' pI must be the same.
' y0 C) Y  E  i) s( Q( P! H( pError in ==> matlab02 at 40
   s(k)=Ce*E;

都天锴

不知道。。。。

厚积薄发

维数不匹配

deven1985

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