$ E( n4 o' F4 V2 k8 IM=3; ; |9 i9 f# X0 E5 P2 O' rif M==1 %EXP reaching law8 H5 [# x- @8 I
ds(k)=-eq*ts*sign(s(k))-q*ts*s(k); . ]$ U0 R! R# B+ R. x3 W# J7 v+ D A u(k)=inv(Ce*B)*(Ce*R1-Ce*A*x-s(k)-ds(k)); 2 V) i @0 j' q D' O/ [elseif M==2 %Variable rate reachine law$ L9 B8 k5 I/ e9 n+ s4 Z
ds(k)=-eq*ts*X1*sign(s(k)); $ ]2 `. G8 u m8 O( Q' h8 z% z u(k)=inv(Ce*B)*(Ce*R1-Ce*A*x-s(k)-ds(k)); ( Z" f6 n+ q; T0 S* K' A) Celseif M==3 %Coposite reaching law % C3 ^, O/ @' H/ L5 s A+ s" | k0=0.60; ( k! B2 C" f# u3 H+ w if X1>k0 %EXP reachine law# R& T/ g, F9 A* F7 {1 d
ds(k)=-eq*ts*sign(s(k))-q*ts*s(k); # Q8 J7 ^$ P+ D* U+ n# j& A u(k)=inv(Ce*B)*(Ce*R1-Ce*A*x-s(k)-ds(k)); 7 x' X( R! h' k5 I1 V elseif X1<=k0 %Variable rate reachine law # ?' M9 j w9 B! g- G8 ]& F, X ds(k)=-eq*ts*X1*sign(s(k));. i( L8 @* |' H ^* v5 V! n8 C
u(k)=inv(Ce*B)*(Ce*R1-Ce*A*x-s(k)-ds(k)); + i, l6 b# u% v" F# V end $ N" Y8 J" w/ ~end " Q8 S' A0 @' Rif u(k)>=100 f/ i. i. u: r( v
u(k)=10;; |% V5 g [ v
end 4 M6 P, e+ C) cif u(k)<=-10 + o: A0 I2 C3 r* v& [; ]5 \ u(k)=-10;# W* N0 s8 D% ^
end0 d2 h( H! ^6 L9 l
x=A*x+B*u(k); " d1 H$ V0 L2 ~0 By(k)=x(1); $ @4 A/ p0 M& J4 s4 ? Y%Update Parameters$ A/ t$ n$ U& P* c
r_2=r_1; 4 W# c M1 m# `/ s: G/ }5 [) xr_1=r(k);: t/ q% v, S9 j
end+ y; c( F: q- u* x3 Z3 V: _
figure(1): R8 A0 {2 A9 x/ c( j# u/ g. Z
plot(time,r,'r',time,y,'b');2 {" h. p2 R# Y% G L: q6 e) l! D
xlabel('Time(second)');ylabel('Position tracking'); 7 u( z, I. E: K- M n# jfigure(2)& [ ]2 A+ L5 G5 D6 H& ~3 [ y- @
plot(time,s,'r');9 ?- e. n! K# ]. C/ }/ U
xlabel('Time(second)');ylabel('Switch function s'); 7 h( n( B; X! M, x6 n3 L- Zfigure(3) 3 u2 S, T6 `. d' S3 V. u* cplot(e,de,'r',e,-c*e,'b'); 6 _: J, J3 J0 y! X8 g& oaxis([0,0.001,-0.01,0]);* m7 V3 R: L+ w6 w' Q3 C* o. E" }
xlabel('e');ylabel('de');7 @* e' D2 H' G L8 b
figure(4) # \8 s" D( W9 r, |plot(time,u,'r');0 Q6 T8 V8 H8 z6 u6 G
xlabel('Time(second)');ylabel('u');7 w, h# T5 H9 W3 k6 u( ]
?% I& V% q* r# ?4 E??? In an assignment A(I) = B, the number of elements in B and- N* J) B9 E7 D4 i! R
I must be the same.8 n: F( ^$ r( G5 |$ R
Error in ==> matlab02 at 40% G7 Q" E" Z0 g- P: u# ?9 f- A
s(k)=Ce*E;