test1.m

用matlab编的用来实现追踪控制的源程序

%the paramenter is p,synchronization plot,and error plot
clc
clear
clf
global p;

p=1;
%for p=0.25:0.5:1.5
    p
%[T,Y]=ode45(@liutrack3,[0 50],[-10 5 10 2.4 2.2 0.8]);%new
[T,Y]=ode45(@liutrack2,[0 20],[-10 5 10 2.4 2.2 0.8]);%chen
%[T,Y]=ode45(@liutrack1,[0 20],[2 1 1 1 3 -1]);%liu
%[T,Y]=ode45(@liutrack3,[0 20],[-10 5 10 2.4 2.2 0.8]);%new test
plot(T,Y(:,1),'k:',T,Y(:,4),'k');
%plot(Y(1000:end,1),Y(1000:end,4))
%plot(T,sin(T),'r',T,Y(:,1));
xlabel('\itt /s','fontsize',18)
ylabel('\itx_{1}(t), r(t)','fontsize',18)
legend('\itr(t)','\itx_{1}(t)','fontsize',18)
%text(1,1.1,'r(t)','fontsize',14)
%text(1,-4.2,'x1(t)','fontsize',14)
%hold on;
for i=1:length(T)
   e(i)=Y(i,1)-Y(i,4);
end
figure 
plot(T,e(1:length(T)),'k');
xlabel('\itt /s','fontsize',18)
ylabel('\ite(t)','fontsize',18)
%hold on
%end