for t=6:N-1
temp=0;
for i=1:5
temp=temp+abs(P(t-i+1)-P(t-i));
end
judge1=(1/5)*temp;
judge2=abs(P(t)-P(t-1));
if (judge1<0.0225) &&(judge2<0.05)
Xt=0.025;
else
Xt=0.05;
end
Dt=Xt-(2/Xt)*(P(t)-P(t-1))^2;
P(t+1)=P(t)+Dt;
end
P=P';
figure
plot(P)
msgbox('MATLAB编程答疑,请加QQ: 1530497909','MATLAB答疑','help')
web http://url.cn/NSFcAs -browser
for t=6:N-1
temp=0;
for i=1:5
temp=temp+abs(P(t-i+1)-P(t-i));
end
judge1=(1/5)*temp;
judge2=abs(P(t)-P(t-1));
if (judge1<0.0225) &&(judge2<0.05)
Xt=0.025;
else
Xt=0.05;
end
Dt=Xt-(2/Xt)*(P(t)-P(t-1))^2;
P(t+1)=P(t)+Dt;
end
P=P';
figure
plot(P)