Discussion :

[Jebosson]

Bonjour à tous, je suis nouveau sur le forum.

Voilà ma question : nous avons fais une expérience de chimie et donc avons reporté les résultats sur un fichier excel. Nous devons importer ces données sur Matlab et en faire pas mal de choses avec. Mais dans mon travail je dois expliciter la question : quels sont les avantages et désavantages de xlsread par rapport a uiimport ?

Si vous avez des réponses je suis preneur.
Merci d'avance !
Jérémy Bosson

[Gooby]

Bonjour,

As tu testé les deux fonctions sur l'un de tes tableaux de données? Tu peux commencer à comparer par toi même les résultats qu'elles te donnent. Tu peux également aller voir un descriptif des fonctions en utilisant l'aide MATLAB.
Petite piste, xlsread ne va servir qu'à mettre les valeurs numériques de tes fichiers dans un tableau, en remplaçant les chaînes de caractères par des NaN, tandis que uiimport va par exemple permettre de directement séparer les colonnes en différents vecteurs, gérer les noms des variables en utilisant les chaînes de caractères présent dans ton fichier excel, etc... L'utilité de ces deux fonctions dépend avant tout de la forme de ton fichier excel, si ce n'est qu'un tableau de valeur, xlsread est suffisant selon moi.

EDIT:

xlsread ne va servir qu'à mettre les valeurs numériques de tes fichiers dans un tableau

est vrai uniquement dans le fonctionnement basique de la fonction xlsread, c'est à dire

Code :

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

A=xlsread('test.xls');

[Jebosson]

Tout d'abord merci de la réponse.
Donc dans notre projet, les deux méthodes ont dues être testées. Je vais allez voir l'aide matlab aussi oui c'est une bonne piste. Xls ne fait que de mettre nos valeurs numériques en tableaux et remplacer les chaînes de caractères par des NaN ? Il ne permet pas de faire d'autres opérations ? Merci pour la piste je vais creuser et vous reposer des questions si besoin . Et je suis coincé avec mon script aussi, je sais pas si j'ose vous l'exposer ?

[Gooby]

Si tu exécutes la commande

Code :

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

A = xlsread('monTableauDeDonnées.xls');

la variable A te renverra une matrice de nombres comme celle existante dans ta feuille excel, à ceci près que si par exemple tu as mis dans la première ligne de ton fichier excel à quoi correspondent ces colonnes (temps en s, pH...) la première ligne de A dans MATLAB sera alors composé de NaN.

Il existe cependant des options spécifiques à la fonction qui te permettront d'utiliser la fonction xlsread de façon un peu différente. Par exemple si tu effectues la commande

Code :

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

[num,txt] = xlsread('monTableauDeDonnées.xls')

num correspondra à A de l'exemple précédent tandis que txt sera composé des différentes chaînes de caractères présent dans ta feuille excel. Mais pour avoir une idée précise du fonctionnement de la fonction xls, je t'invite à aller la découvrir dans l'aide de MATLAB, qui est complète et bien rédigé.

Ensuite, tu peux travailler sur ta matrice A/num comme bon te semble est effectué toute sorte de calcul.

Avant de nous exposer ton script, tu pourrais commencer par expliquer quel est le problème que tu rencontres peut être

PS: je ne suis pas un expert MATLAB, il se peut que des subtilités m'échappent ou que mes explications soient imprécises voire légèrement fausse. Un avis complémentaire pourrait être utile afin d'avoir une réponse plus précise et plus juste.

[Jebosson]

D'accord c'est parti

Alors voila nous avons fait deux dosage d'un même mélange KI/KCl par potentiométrie. Donc nous avons dans notre tableau excel une colonne avec les volumes de solution titrante et deux autres colonnes avec les millivolts de nos deux expériences. Nous devons en un premier temps charger les données, puis éloigner les valeurs aberrantes et ensuite repérer les données manquantes s'il y en a.
Ensuite nous devons interpoler pour combler les trous et lisser les données pour éliminer les bruits.
Ensuite il y a quelques calculs statistiques à faire ( moyenne, ecart type, médiane, ...) et grâce à cela déterminer un intervalle de vraisemblance. Finalement nous devons déterminer les points d'équivalence de nos titrations par rapport au plot que nous avons effectuer un peu plus haut. La première méthode est avec une série de droite et il suffit de trouver un croisement de droite pour déterminer ce point d'équivalence.
La deuxième méthode pour le deuxième point d'équivalence me pose problème car il faut travailler sur la première et deuxième dérivée de notre courbe pour trouver des maximums et en définir ou est le changement de signes pour savoir exactement ou est notre point cherché.

Voila une idée de ce que je dois faire. C'est compréhensible ?

[Gooby]

Oui c'est assez clair, maintenant que ton sujet est bien exposé, peux tu nous dire où se situe le problème dans ton script, puisque a priori, tu as déjà au moins un commencement de script?

Comme tu peux le voir, je suis nouveau aussi sur ce forum (août 2014!), je ne sais pas si tu dois créer un nouveau topic pour ton problème ou pas, mais on va supposer que tu peux poster ton script sur ce même sujet

[Jebosson]

Voilà on script :

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
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
[d]=xlsread('Projet1.xlsx','Sheet1');
 
mv1 = d(:,2);
mv2 = d(:,3);
vol = d(:,1);
 
lignesatrou = any(isnan([mv1,mv2]),2);
vol(lignesatrou) = [];
mv1(lignesatrou) = [];
mv2(lignesatrou) = [];
 
pas_de_trou = ~isnan(mv1);
mv1 = interp1(vol(pas_de_trou),mv1(pas_de_trou), vol);
 
pas_de_trou = ~isnan(mv2);
mv2 = interp1(vol(pas_de_trou),mv2(pas_de_trou), vol);
 
outliers =([mv1 mv2] - mean([mv1;mv2]) >= 3*std([mv1;mv2]));
 
mv1(outliers(:,1))= NaN;
mv2(outliers(:,2))= NaN;
 
span = 0;
lissage = conv([mv1;mv2],ones(span,1)/span, 'same');
 
moyenne = mean([mv1 mv2]')';
 
ecart_type = std([mv1 mv2]')';
 
intervallevraisemblance = ecart_type/sqrt(2)*t(0.05,1);
 
errorbar(vol,moyenne,intervallevraisemblance);
 
title('Dosage halogénure [KI et KCl]')
xlabel('Volume [mL]')
ylabel('Potentiel [mV]')
 
pv1 = polyfit(vol(10:11), mv1(10:11),1);
ph1 = polyfit(vol(15:22), mv1(15:22),1);
 
pv2 = polyfit(vol(10:11), mv2(10:11),1);
ph2 = polyfit(vol(15:22), mv2(15:22),1);
 
x1 = (ph1(2)-pv1(2))/(pv1(1)-ph1(1));
x2 = (ph2(2)-pv2(2))/(pv2(1)-ph2(1));
moyennePtEqui = mean(x1:x2);
 
d1_exp1 = gradient(mv1, vol);
d2_exp1 = gradient(d1_exp1, vol);
tmp = [ vol d1_exp1 d2_exp1 ];
 
 
d1_exp2 = gradient(mv2, vol);
d2_exp2 = gradient(d1_exp1, vol);
tmp = [ vol d1_exp2 d2_exp2 ];

Ma difficulté c'est que je ne sais pas comment continuer avec mes 1eres et 2emes dérivées pour trouver mon point d'équivalence. Est-ce que tu aurais une idée ?
Merci de ton aide.
Oui j'espère que je peux poster cela ici !

[Gooby]

Je ne suis pas sûr de bien comprendre, as-tu des difficultés pour effectuer ce que tu souhaites sur MATLAB ou est-ce un problème extérieur à MATLAB, c'est à dire que tu ne sais pas ce que tu dois faire? d'après ce que j'ai compris jusque là, il te suffirait de chercher les zones de croissance/décroissance de ta dérivé première, en regardant le signe de la dérivé seconde.


très schématiquement, sur un exemple quelconque, cela donne quelque chose du genre:

Code :

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

1
2
3
4
                           0                          x                              +inf    
 
Dérive seconde:                        +              |                 -
Dérivé première:                    croissante        |          décroissante

Donc x est un maximum sur ta fonction initiale. => Mes excuses, cette assertion est fausse. ça représente un maximum sur ta fonction dérivé première et un point d'inflexion sur ta fonction initiale.

[Jebosson]

Je ne vois pas comment je dois écrire le code pour trouver mon deuxième maximum de ma dérivée seconde qui sera donc mon point d inflexion et par conséquent mon point d'équivalence tu vois ?

[Gooby]

Selon wikipédia on a

En un point d'inflexion la dérivée seconde f", si elle existe, s'annule et change de signe. Ceci permet de vérifier si des points sont des points d'inflexion.

Donc

Code :

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

indice = find(dérivé_seconde == 0);

[Jebosson]

Je vais essayé cela et te redis merci beaucoup !

[Gooby]

Cela n'a pas marché n'est-ce pas?

C'est probablement parce que tu effectues tes calculs sur des fonctions discrètes, c'est à dire sur des intervalles comportant un nombre de valeurs fini. Donc il est fort probable que ta dérivé seconde ne soit jamais égale à 0. Il faut donc que tu cherches la valeur la plus proche de 0 dans ton vecteur comportant les valeurs de ta dérivé seconde pour avoir une approximation.

Code :

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

indice = find(données == min(abs(données)))

[Jebosson]

Oui ca n'as pas marché et je pense que tu as raison par rapport a que la dérivée n'atteint jamais 0. Mais le 'données' je le remplace par quoi ?

[Gooby]

Dans ton cas, on dirait bien que tu as utilisé la fonction gradient pour effectuer tes dérivés, si j'ai bien compris? donc pour ton expérience1, tes données sont dans le vecteur d2_exp1 si je ne m'abuse.

[Jebosson]

Code :

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

1
2
3
4
5
6
7
8
d1_exp1 = gradient(mv1, vol);
d2_exp1 = gradient(d1_exp1, vol);
x3 = find(d2_exp1 = min(abs(d2_exp1)));
 
 
d1_exp2 = gradient(mv2, vol);
d2_exp2 = gradient(d1_exp1, vol);
x4 = find(d2_exp2 = min(abs(d2_exp2)));

mon code reviendrait à cela.
Mais voila ce qu'on me dit pour la ligne

Code :

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

x3 = find(d2_exp1 = min(abs(d2_exp1)));

Error: File: Projet1.m Line: 50 Column: 19
The expression to the left of the equals sign is not a valid target for an assignment.

[Gooby]

Oui en effet, j'ai commis une erreur. Le = seul c'est pour l'assignation, pour la comparaison utilise le double = ( donc == )

[Jebosson]

juste j'aurais du voir ça de moi-même

[Gooby]

Tes points d'inflexion sont donc en toute logique [ vol(x3) mv1(x3) ] et [ vol(x4) mv2(x4) ]

[Jebosson]

Tu crois pas que dans ma courbe de titration, j'ai plusieurs endroits ou la dérivée seconde se rapproche de 0 ?

Si tu veux je dois savoir à quel volume se trouve mon point d'équivalence donc à la coordonnée x de mon graphique

[Gooby]

Si tu as une courbe de cette forme, tu auras en effet deux points d'inflexion.


Dans ce cas, tu peux essayer de découper ta plage de données sur deux intervalles puis effectuer la recherche du point d'inflexion sur chacun de ces intervalles.

PS: Je n'ai pas fait de chimie depuis déjà un moment, le point d'équivalence est il le même que le point d'inflexion de ta courbe?