Discussion :

[magictol53]

Bonjour à tous,

Voilà, j'ai extrait une sous matrice d'une matrice pour générer un modèle
avec cette sous matrice...Je voudrais maintenant tester le modèle sur la ligne qui n'a pas servi à générer le modèle.
J'y suis presque mais quelque chose m'échappe.
load croisement.txt;

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
17
18
19
20
21
for n= 1:212;
V = [croisement(1:n-1,:);croisement(n+1:212,:)];
%Construction des vecteurs x et y
x = V(:,3:4);
y = V(:,1:2);
%Recalage des valeurs entre 0 et 1
[x Mx mx] = scale(x);
[y My my] = scale(y);
%Construction de la matrice de transformation
a = x(:,1);
b = x(:,2);
one = ones(length(a),1);
T = [a b one];
% obtention de la matrice des coefficents  (coeff = (M*T)\(M*xx))
P = (T'*T)\(T'*y);
yp=T*P;   % c'est à cet endroit qu'il faut que j'indique que T  prenne seulement la ligne qui n'ait pas servi à générer le modèle
yp = descale(yp,My,my);
y = descale(y,My,my);
rmse = sqrt(mean(sum((yp-y).^2,2)));
M(n)=rmse;
end

Merci d'avance

MG

[Jean Dumoncel]

Bonjour,

c'est simplement :

Code :

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

[croisement(n,3:4) 1]

non?

[ol9245]

c'est quoi ta question ?

[magictol53]

Ma question est de savoir comment récuperer la ligne qui n'a pas servi à générer le modèle pr pouvoir la transformer.

[ol9245]

Est-ce que la réponse de Magellan t'a aidée ?

La question telle que tu la poses est spécifique à ton problème. Ce qu'on a besoin pour t'aider c'est une question davantage exprimée en terme de Matlab : une question à laquelle on pourrait répondre sans avoir eu à comprendre tous les détails spécifiques à ton problème.

Par exemple si tu nous montre un code qui marche pas, et l'erreur que ça te sort, c'est + facile.

[magictol53]

En fait, le code que j'ai inséré au départ fonctionne. Le problème est que quand je calcule la prédiction (yp), je la calcule pour les lignes qui ont servi à determiner les coefficients de transformation (P).Ce n'est dc pas une bonne erreur de prédiction..Moi, j'aimerai juste recuperer la ligne exclue de T et la multiplier par P (P étant determiné à l'aide des autres lignes).

Merci d'avance

MG

[Jean Dumoncel]

La ligne exclue de T, c'est croisement(n,:), donc ce que je t'ai donné devrais convenir, non?

[magictol53]

Non, cela ne marche pas, voici l'erreur..

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
17
18
19
20
21
load croisement.txt;
for n= 1:212;
V = [croisement(1:n-1,:);croisement(n+1:212,:)];
%Construction des vecteurs x et y
x = V(:,3:4);
y = V(:,1:2);
%Recalage des valeurs entre 0 et 1
[x Mx mx] = scale(x);
[y My my] = scale(y);
%Construction de la matrice de transformation
a = x(:,1);
b = x(:,2);
one = ones(length(a),1);
T = [a b one];
% obtention de la matrice des coefficents  (coeff = (M*T)\(M*xx))
P = (T'*T)\(T'*y);
yp=[croisement(n,3:4) 1]*P;   % c'est à cet endroit qu'il faut que j'indique que T  prenne seulement la ligne qui n'ait pas servi à générer le modèle
y = descale(y,My,my);
rmse = sqrt(mean(sum((yp-y).^2,2)));
M(n)=rmse;
end

??? Error using ==> minus
Matrix dimensions must agree.

[Jean Dumoncel]

C'est normal que ça ne marche pas, il faut que tu sélectionne uniquement le y de la ligne n...

[magictol53]

Non ca ne fonctionne pas, je n'obtiens qu'une seule ligne pour yp, je devrais en obtenir 212....