curve fitting problème

Version imprimable

Bonjour tout le monde!
Je voulais faire la même chose que dans le thread de Ayvik

http://www.developpez.net/forums/d78...-modele-connu/

mais j'obtiens une erreur de ce type :

Citation:

??? Index exceeds matrix dimensions.

Error in ==> projet3 at 5
x= file (:,1);

mon bout de code est:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 filename = 'OSA_350mA.xls'; file= xlsread (filename); x= file (:,1); y= file (:,2); options = optimset('lsqcurvefit'); options = optimset(options,'LevenbergMarquardt' ,'on'); options = optimset('LargeScale','off'); %options.MaxFunEvals = 10000; %options.MaxIter = 5000; [newParameters,error] = lsqcurvefit(@myFunction, [0 0 0 0 0 0],(x-mean(x))/std(x),y,[],[], options); y2 = myFunction(newParameters,(x-mean(x))/std(x)); figure; scatter(x,y,'.') ; hold; plot(x,y2,'r');
Je précise que je suis un pur newbee en Matlab, mais si il y a moyen de me guider sa serait super cool!
A+

03/08/2009, 09h07
Caro-Line

Le mieux c'est de vérifier les résultats de chaque ligne.
En l'occurrence il semblerait que la variable file est vide.

As-tu compris le principe de XLSREAD ? Lu l'aide ?

A quoi ressemble ton fichier xls ?
03/08/2009, 15h26
MrJack14
Hello
Alors oui la fonction xlsread permet de lire les colonnes de mon fichier xls, je n'ai que deux colonnes dans mon fichier, le nombre de ligne va jusqu'à 5000.
J'ai bien mis le fichier executable matlab dans un fichier contenant bien le xls, je me suis dis que le fait cela:

Code:

load OSA_350mA.xls;

me permettrai de bien insérer mes valeur de colonne dans la variable filename mais il me met cette erreur là :
Code:

1 2 3 4 5 6 Unknown text on line number 1 of ASCII file C:\Users\Rémy\Desktop\Projet avec Pascal\projet3\OSA_350mA.xls "ÐÏà¡±á". Error in ==> projet3 at 1 load OSA_350mA.xls;
ma première valeur dans la ligne 1 est 1511.000 dans la première colonne est la valeur dans la deuxième colonne est -64.620.

Voilà
03/08/2009, 15h46
Jerome Briot
Tes explications ne sont pas très claires... tu dis que tu utilises XLSREAD et après tu essayes avec LOAD... or la fonction LOAD ne sert pas à lire les fichiers xls.

Quel problème as-tu avec XLSREAD ?

Montre nous ce que retourne ceci :
Code:

1 2 3 4 filename = 'OSA_350mA.xls'; file = xlsread (filename); whos file
03/08/2009, 15h56
MrJack14
Hello
Alors il me retourne cela
Code:

1 2 3 Name Size Bytes Class Attributes file 0x0 0 double
donc si j'ai compris il ne repère rien parce que mon fichier xls et bel et bien remplis!
03/08/2009, 15h59
Jerome Briot

Pourrais-tu attacher ton fichier xls à un message ?

Sinon, quelle est ta version de MATLAB et sous quel système d'exploitation travailles-tu ?
03/08/2009, 16h06
MrJack14

1 pièce(s) jointe(s)

Hello
J'utilise la version 7.5.0 et je travaille sous Vista
03/08/2009, 17h21
Jean Dumoncel

Salut,

Je n'y connais pas grand chose en fichier excel, mais j'ai essayé de regarder.

Les données sont écrites avec des points pour indiquer les décimales, est-ce normal?
Je demande car dans open-office, en convertissant les données en changeant les points par des virgules pour les 3 premières lignes, les données se sont chargées correctement dans Matlab...
03/08/2009, 17h43
MrJack14

Hello
En effet sa se charge mais si je fais sa je perd en précision car si tu as remarqué le troisième chiffre après la virgule est perdu...
03/08/2009, 17h46
Jean Dumoncel

J'ai pas remarqué car je l'ai fait à la main... (c'est pour ça que j'ai pris que les trois premières lignes...

Comment as-tu fait pour modifier tout le fichier? N'y-a-t-il pas une option qui permet de garder la précision?

----------------
EDIT : une solution de secours serait de l'enregistrer au format texte (csv ou txt) et de faire un load ou un csvread.
03/08/2009, 18h12
MrJack14

1 pièce(s) jointe(s)

J'ai réussi à mettre tout à la virgule mais rien n'y fais et en plus il ne se charge même plus la preuve regarde dans ma pièce jointe
03/08/2009, 18h31
Jean Dumoncel

Sans dire ce que tu fais, c'est un peu dur de suivre...

Dans ton nouveau fichier, tes données sont stockés sous forme de texte... :cfou:
03/08/2009, 18h38
Caro-Line

Avant de vouloir modifier quoi que ce soit.
En fait le problème vient du format de tes données : les cellules sont vraisemblablement en format texte et non pas nombre.

Si tu lis l'aide de XLSREAD le 1er argument de sortie ne peut être que des nombres.
Sur ton fichier si je fais ceci :

Code:

[A,B]=xlsread('OSA_350mA.xls')

J'obtiens toutes les données dans B sous forme de tableau de cellules de chaines de caractères.

2 solutions :
- tu fais comme ça puis tu convertis B :

Code:

M=cellfun(@(x)str2double(x),B)

- tu appliques le format numérique à tes cellules de ton fichier Excel
03/08/2009, 22h14
MrJack14

Alors voilà j'ai enfin réussi à convertir mes valeurs j'obtiens cela :

Citation:

Name Size Bytes Class Attributes

file 5001x2 80016 double

par contre j'obtiens sa maintenant :

Citation:

??? Error using ==> lsqncommon at 101
LSQCURVEFIT cannot continue because user supplied objective function failed with
the following error:
Error using ==> feval
Undefined function or method 'myFunction' for input arguments of type 'double'.

Error in ==> lsqcurvefit at 182
[x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...

Error in ==> projet3 at 13
[newParameters,error] = lsqcurvefit(@myFunction, [0 0 0 0 0
0],(x-mean(x))/std(x),y,[],[], options);

Encore et encore des erreurs!
03/08/2009, 22h56
Caro-Line
Il faut lire attentivement la documentation de LSQCURVEFIT.
Le premier paramètre est un handle d'une fonction qui doit être définie dans un fichier à part.
Donc cette partie là :
Code:

1 2 3 4 5 6 y2 = myFunction(newParameters,(x-mean(x))/std(x)); figure; scatter(x,y,'.') ; hold; plot(x,y2,'r');
doit s'écrire :
Code:

1 2 3 4 5 6 function y2 = myFunction(newParameters,(x-mean(x))/std(x)); figure; scatter(x,y,'.') ; hold; plot(x,y2,'r');
et être enregistrée dans un fichier myFunction.m
03/08/2009, 23h57
Jean Dumoncel

Tout d'abord, merci Caro-Line d'avoir rectifié le tir, j'étais parti un peu à l'ouest dans la gestion de ce fichier excel...

Ensuite, MrJack14, il faut que tu prennes le temps de comprendre ce que tu cherches à faire et ce que tu codes. Regardes le lien donné par Ayvik (qui lui a servi de base pour résoudre son problème d'ajustement) :
http://blinkdagger.com/matlab/matlab...-lsqcurvefit/2
Prends le temps d'exécuter les codes sur cette page et de comprendre chaque partie.
Ensuite essaye de l'adapter à ton problème, et notamment réfléchie à ce que tu vas mettre dans myFunction puisque c'est cette fonction que tu vas essayer d'ajuster à ton nuage de points.
04/08/2009, 00h10
MrJack14

Merci!
Donc ce que je fais c'est que je réecris cela et je dois créer une fonction qui sera appelée par le bout de code que tu m'as réécris c'est sa?
09/08/2009, 20h06
MrJack14

Hey! Bon finalement je remettrai ce projet à plus tard, n'ayant pas eu le temps avec le projet que je fais, je me suis baser sur Curve Fitting toolbox pour faire cela, je reprendrai mon code plus tard.