Discussion :

[pascale89]

Je cherche à écrire un programme matlab qui détermine les paramètres d'une machine par les différentes méthodes de moindre carrés.
Ca marche avec les moindres carrés simples.

Arrivé à la méthode des moindres carrés récursifs; quelque soit le programme que j'écrit ça me rend que la matrice P est singulière donc non inversible.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
teta = [] ;
%Calcul de Teta1
phi1=PHI(1,:)'; %PHI est la matrice d'information
P=inv(phi1*phi1');
teta(:,1)=P*phi1*i(1);
 
%Calcul de Teta(k+1)
for k=2:N
    P=P-[P*PHI(k,:)'*PHI(k,:)*]/[1+PHI(k,:)*P*PHI(k,:)'];
    teta(:,k)=teta(:,k-1)+P*PHI(k,:)'*[i(k)-PHI(k,:)*teta(:,k-1)];
end
 
plot(t,teta),grid

[djelfa]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
clear all;
clc;
N=1122;
data=load ('process.dat');
y=data(:,2);%les valeurs de 2eme colonne
u=data(:,3);%les valeurs de 3eme colonne
%initial0 0 values of algorithem
teta=zeros(6,1);
gain=eye(6)
for i=3:N
    % vous avez le droit de choisir l'ordre
phi=[-y(i) -y(i-1) -y(i-2) u(i) u(i-1) u(i-2)]'
gain=gain-(gain*phi*phi'*gain)*(inv(1+phi'*gain*phi))
teta=teta+gain*phi*(y(i+1)-teta'*phi)
end