clc;
A0=[0?0;0?0];
A1=[-0.4?-0.3;0.1?0.15];
A2=[0.1?0.25;-0.9?-0.1];
h=6;?
d=5;
N=150;
P=zeros（2N+1）;
P（:1）=[-0.1;0.1];
if?N>h
????for?i=1:d
????????P（:i+1）=A0*P（:i）+A1*P（:1）+A2*P（:1）;
????end
????for?i=d+1:h
????????P（:i+1）=A0*P（:i）+A1*P（:1）+A2*P（:i-d）;
????end
????