function rrf2008

close all;clear,clc

%  2.和对比效果图
m=100;
t1=0;t2=0;
y1=0.5;y2=0.5;y3=0.5;y4=0.5;
y11=0.6;y22=0.6;y33=0.6;y44=0.6;
% x的图形对比
for n=1:m
    t1=(n-1)*0.3;
    t2=n*0.3;
    [T1,Y1]=ode45(@cdsf1,[t1,t2],[y1 y2 y3 y4]);
    [T2,Y2]=ode45(@cdsf1,[t1,t2],[y11 y22 y33 y44]);
    [H1,L1]=size(Y1);[H2,L2]=size(Y2);
    y1=Y1(H1,1);y2=Y1(H1,2);y3=Y1(H1,3);y4=Y1(H1,4);
    y11=Y2(H2,1)-0.0111*(Y2(H2,1)-Y1(H1,1));
    y22=Y2(H2,2)-0.0111*(Y2(H2,2)-Y1(H1,2));
    y33=Y2(H2,3)-0.0111*(Y2(H2,3)-Y1(H1,3));
    y44=Y2(H2,4)-0.0111*(Y2(H2,4)-Y1(H1,4));
    plot(T1,Y1(:,1))
    hold on
    plot(T2,Y2(:,1),'r:')
    hold on
end

end


% 1.微分方程函数文件

function dy=cdsf1(t,y)
dy=zeros(4,1);
dy(1)=-y(1)*0.9+(y(2)*y(1)+y(4)+y(3));
dy(2)=1-y(2)*0.2-y(1)^2;
dy(3)=-y(3)*1.5-y(1);
dy(4)=-y(4)*0.17-y(1)*y(2)*0.2;

end

rrf2008