离散时滞系统MATLAB作图

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