clear;
time=10;
am1=0.8;
am2=0.5;
de=0.008;
Samtime=0.496
hold off;
subplot(1,2,1);
sim('Rakesimulink.mdl') ;
mit=496/1000;
x=0:mit:(length(attunation1)-1)*mit;
plot(x,attunation1,'b');
hold on;
plot(x,attunation2,'r');
xlabel('time/s');
ylabel('attunation/%');
title('attunation');
axis([0 time 0 1]);
subplot(1,2,2);
plot(x,delay);
xlabel('time/s');
ylabel('delay/s');
title('delay');
axis([0 time 0 2*de]);
u=0;
for i=1:length(delay)
    if delay(i)>0
        u=u+1;
        if u==1
            position=i;
        end
    else u=u;
    end
end
Psum1=sum(attunation1);
Ave_att1=Psum1/u;
Psum2=sum(attunation2);
Ave_att2=Psum2/u;
Psum3=sum(delay);
Ave_delay=Psum3/u;
Time_discover=position*mit;
Rate_search=u/length(delay);
Ave_att1
Ave_att2
Ave_delay
Time_discover
Rate_search