Discussion :

[abel413]

Bonjour,

J'ai un problème en utilisant un de mes codes que je croyais bon jusque-là mais qui s'avère en fait problématique ! Le voici :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
function R=deb_fin(valeurs_Imp1_liquide)
LL=diff(valeurs_Imp1_liquide);
L3=find(LL>1);
fin=valeurs_Imp1_liquide(L3); 
fin(numel(L3)+1)=valeurs_Imp1_liquide(end);
debut(2:numel(L3)+1)=valeurs_Imp1_liquide(L3+1);
debut(1)=valeurs_Imp1_liquide(1);
debut;
fin;
R=[debut', fin', (fin'-debut')+1];
A=0;

Je vous explique le fonctionnement de ce code : il prend en paramètre un vecteur colonne, que j'appelle valeurs_Im1_liquide, par exemple :

Code :

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

 valeurs_Imp1_liquide=[998 ; 1000 ; 1001];

Il renvoie dans R les segments de longueur 1 et leur durée, soit donc ici :

Code :

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

1
2
3
4
R=
 
998 998 1
1000 1001 2

car 998 est un "segment" de longueur 1 : [998 998] et [1000 1001] est un segment de longueur 2. Mais cela ne marche pas toujours, par exemple avec :

Code :

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

valeurs_Imp1_liquide=[998 ; 1000 ; 1002];

Je devrais avoir :

Code :

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

1
2
3
4
5
R=
 
998 998 1
1000 1000 1
1002 1002 1

Mais à la place j'ai un message d'erreur :

??? Error using ==> minus
Matrix dimensions must agree.
Error in ==> deb_fin at 10
R=[debut', fin', (fin'-debut')+1];

A l'aide du debug mode je vois bien que le problème vient des vecteurs debut et fin qui ne sont plus tous les deux en ligne, car fin de façon incompréhensible est devenue une colonne et donc R ne peut-être écrit ! Je ne comprends pas pourquoi mais j'ai besoin que le code marche à tous les coups, sauriez-vous m'aider ?

[Jerome Briot]

Essaie ceci :

Code :

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

R = [debut(:) fin(:) (fin(:)-debut(:))+1];

Lis aussi ceci : [Défis][1] Vecteur colonne ou vecteur ligne ?

[abel413]

Incroyable Dut cela à l'air de bien marché, je suis scotché là Dut !

Comment est-ce possible ? Est-ce que de la sorte tu mets dans R le contenu en colonne indépendamment du fait que debut et fin sont en colonne ou ligne ? C'est une sorte d'"esquive" en gros non ? J'étais en train de faire un if else en fait mais ce n'est plus nécessaire car ce qui posait problème c'est le fait que L3 possède plus de 1 élément, comme dans mon second exemple.

[abel413]

En fait, pardon j'ai oublié un autre point. Il se peut que le vecteur passé en paramètre de la fonction soit vide. Je souhaiterais alors que le code continue : j'appelle en effet cette fonction dans une autre à plusieurs reprises, sur différents vecteurs. Je viens de tapper quelque chose comme :

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
function R=deb_fin(valeurs_Imp1_liquide)
 
if isempty(valeurs_Imp1_liquide)
    disp('vide');
 
else
LL=diff(valeurs_Imp1_liquide);
L3=find(LL>1);
fin=valeurs_Imp1_liquide(L3); 
fin(numel(L3)+1)=valeurs_Imp1_liquide(end);
debut(2:numel(L3)+1)=valeurs_Imp1_liquide(L3+1);
debut(1)=valeurs_Imp1_liquide(1);
debut;
fin;
% R=[debut', fin', (fin'-debut')+1];
R = [debut(:) fin(:) (fin(:)-debut(:))+1];
end;

Mais le code s'arrête tout de même :

Code :

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

1
2
3
4
5
vide
??? One or more output arguments not assigned during call to 'C:\MATLAB701\work\IMP\DATA\deb_fin.m (deb_fin)'.
 
Error in ==> etude_de_courbe_reflux at 294
evenements_inter_Imp34_melange=deb_fin(inter_Imp34_melange);

Je voudrais donc qu'il me garde le vecteur vide, ici en l’occurrence inter_Imp34_melange et que ma fonction principale continue, puisque par exemple à la ligne suivante, 295, j'appelle encore deb_fin :

Code :

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

evenements_inter_Imp45_melange=deb_fin(inter_Imp45_melange);

[Invité]

Bonjour,

Le fait est que tu retournes dans les deux cas la variable R, mais que tu ne la définis que dans ton else.

Code :

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

1
2
3
4
5
if isempty(valeurs_Imp1_liquide)
    disp('vide');
    R = [];
else
    % ...

Voir la FAQ Output argument "***" (and maybe others) not assigned during call to "***"

[abel413]

Effectivement, merci bien Winjerome.