u(t)	-0,4	-2,3	0,3	0	-0,5	-1	100,2	6,7	9,7	9.0	4,89	89	100	2	-0,67

Y(t) 89.0000 87.9000 83.3000 84.4000 85.4000

R=1;%%ilość y (nuemracja od jeden)

S=2;%%ilość u (numeracja od zera)

%%N=5;

N=R+S+1;

%%wektor, który będziemy szukać - jeśli zmieniamy R i S to tutaj też zmieniamy

p=[1;-1;0;2];

%%wektor wejść

u=[-0.4 -2.3 0.3 0 -0.5 -1 100.2 6.7 9.7 9.0 4.89 89 100 2 -0.67]

%%warunki początkowe

y=[89];

for i=0:N-1

X=[fliplr(y(1+i:R+i)),fliplr(u(1+i:S+1+i))];

y=[y X*p];

end

disp('y =')

disp(y)

Q=[];

for i=0:N-1

Q=[Q ; [fliplr(y(1+i:R+i)),fliplr(u(1+i:S+1+i))] ];

end

Q

P=inv(Q'*Q)*Q'*y(R+1:end)'

u =

-0.4000 -2.3000 0.3000 0 -0.5000 -1.0000 100.2000 6.7000 9.7000 9.0000 4.8900 89.0000 100.0000 2.0000 -0.6700

y =

89.0000 87.9000 83.3000 84.4000 85.4000

Q =

89.0000 0.3000 -2.3000 -0.4000

87.9000 0 0.3000 -2.3000

83.3000 -0.5000 0 0.3000

84.4000 -1.0000 -0.5000 0

Y =

87.9000

83.3000

84.4000

85.4000

P =

1.0000

-1.0000

-0.0000

2.0000