1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44
| y=zeros(T,1);
for i=3:T
if y(i-2,1)<=0
y(i,1)=-0.85*y(i-2,1)+randn(1);
else
y(i,1)=-0.99*y(i-2,1)+randn(1);
end
end
cq=0.856;
mq=[0,051
0,041
0,050
-0,140
0,055
0,034
0,029
-0,097
0,038
0,025
0,012
-0,053
0,021
0,016
-0,005
-0,010
0,000
0,008
-0,002
-0,003
0,000
0,002
0,000
0,000
0,000
0,000
0,000];
zq=zeros(T,1);
for i=28:T-27
zq(i,1)=cq*y(i,1)+[y(i+1,1) y(i+2,1) y(i+3,1) y(i+4,1) y(i+5,1) y(i+6,1) y(i+7,1) y(i+8,1) y(i+9,1) y(i+10,1) y(i+11,1) y(i+12,1) y(i+13,1) y(i+14,1) y(i+15,1) y(i+16,1) y(i+17,1) y(i+18,1) y(i+19,1) y(i+20,1) y(i+21,1) y(i+22,1) y(i+23,1) y(i+24,1) y(i+25,1) y(i+26,1) y(i+27,1)]*mq+[y(i-1,1) y(i-2,1) y(i-3,1) y(i-4,1) y(i-5,1) y(i-6,1) y(i-7,1) y(i-8,1) y(i-9,1) y(i-10,1) y(i-11,1) y(i-12,1) y(i-13,1) y(i-14,1) y(i-15,1) y(i-16,1) y(i-17,1) y(i-18,1) y(i-19,1) y(i-20,1) y(i-21,1) y(i-22,1) y(i-23,1) y(i-24,1) y(i-25,1) y(i-26,1) y(i-27,1)]*mq;
end |
Partager