Discussion :

[yassinebench]

Bonsoir
j'aimerai bien savoir est ce que c'est possible de pouvoir définir le nombre n= ?? de colonne d'une matrice dans une boucle for et lorsque je j’exécute la fonction polytifit pourra me calculé tous les combinaison possible selon le nombre de colonne de matrice que je procède et merci d'avance.

Code :

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

1
2
3
4
5
6
7
8
 
clear
load x.mat
load y.mat
n=??
for i=x(:,2:n)
    p=polyfit(x(:,1),x(:,i??),1);
end

[yassinebench]

comme sa je peux avoir la dimension pour le n mais la fonction ployfit ne s’exécute pas ??
ou juste me donne un seul resultats

Code :

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

[m,n]=size(x);

[Jerome Briot]

Montre nous le code que tu exécutes.

[yassinebench]

Code :

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

1
2
3
4
5
6
7
8
clear
 
load x.mat
load y.mat
[m,n]=size(x);
for i=2:n;
    p=polyfit(x(:,1),x(:,i),1)
end

[Jerome Briot]

Pour stocker les valeurs des coefficients dans la variable p à chaque itération, tu dois faire comme ceci :

Code :

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

1
2
3
for i = 2:n
    p(i-1,:) = polyfit(x(:,1), x(:,i), 1);
end

[yassinebench]

pourquoi tu as mis i-1

Code :

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

p(i-1,:)= polyfit(x(:,1), x(:,i), 1);

et j'ai exécuté la fonction sa marche
mais lorsque je lai exécute avec une matrice de trois colonne le resultats que j'ai obtenu et juste une matrice 2*2 alors que je doit obtenir une matrice 3*2

merci

[Jerome Briot]

pourquoi tu as mis i-1

Parce que la boucle démarre à i = 2. Si j'avais écris p(i) = ...,la première ligne de p n'aurait contenu que des 0.

mais lorsque je lai exécute avec une matrice de trois colonne le resultats que j'ai obtenu et juste une matrice 2*2 alors que je doit obtenir une matrice 3*2

Tu fais une approximation avec un polynôme d'ordre 1 (soit y = ax+b). Tu obtiens donc deux coefficients (a et b) à chaque itération.
La variable p est donc de dimension (n-1)x2 (avec n : nombre de colonne de la variable x dans ton problème)

Si tu veux un tableau nx2 alors il faut mettre :

Code :

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

p(i,:) = polyfit(x(:,1), x(:,i), 1);

Et donc relire ma première remarque de ce message