test.m

来自「用matlab编的用来实现追踪控制的源程序」· M 代码 · 共 28 行

clc
clear
clf
global p;

p=1;
%for p=0.25:0.5:1.5
    p
[T,Y]=ode45(@liutrack,[0 20],[2 1 1]);
%plot(T(1:400),Y(1:400,1),'r',T(1:400),Y(1:400,4),'b');
plot(T,sin(T),'k:',T,Y(:,1),'k');
xlabel('\itt /s','fontsize',18)
ylabel('\itx1(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)-sin(T(i));
end
figure 
plot(T,e(1:length(T)),'k');
xlabel('\itt /s','fontsize',18)
ylabel('\ite(t)','fontsize',18)

hold on
%end