load("*****.sod"); // charge les données exp.
N=size(u,2); // affecte à N le nbre de col. de u
n=2; // Ordre du système dont on souhaite estimer les paramètres

// Construction de Y
Y=y(1,n+1:N)';

// Construction de Phi
// Plusieurs sol.: ligne par ligne, col. par col.,...
// col. par col.
Phi=zeros(N-n,2*n+1);
for i=1:n,
    Phi(:,i)=-y(1,i:N-n+i-1)';
end
for j=1:n+1,
   Phi(:,j+n)=u(1,j:N-n+j-1)';
end

// Réduction de Phi en supprimant les col. de para. nuls
AConserver=[1 2 3 4];// tableau contenant les indices des col. à conserver
for i=1:size(AConserver,2),
    Phir(:,i)=Phi(:,AConserver(1,i));
end

// Calcul et affichage de ThetaChapeau
disp("ThetaChapeau :")
ThetaChapeau=inv(Phir'*Phir)*Phir'*Y;
disp(ThetaChapeau);

disp("J(ThetaChapeau) :")
J=(Y-Phir*ThetaChapeau)'*(Y-Phir*ThetaChapeau);
disp(J);