Discussion :

[kwatz]

Bonjour,

Je m'explique, j'ai un graph, et je veux faire l'intégrale sur une partie de mon graph.

Je souhaite rentrer manuellement les bornes, mais ce que je rentre ne correspond pas forcément à une valeur existante. Je cherche alors la valeur la plus proche de ma borne souhaitée.

Par avance,
merci

[Jean Dumoncel]

Salut,

je ne suis pas sûr d'avoir compris,
pourrais-tu nous mettre un petit exemple?

[kwatz]

Code :

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

1
2
3
4
5
6
7
8
9
 
plot(TimeT,HFT); grid on
BI = input('Quelle est la borne inferieure? \n')
find(BI < TimeT & TimeT < BI+1)
BII = input('Choisir la borne \n')
 
BS = input('Quelle est la borne superieure? \n')
find(BS < TimeT & TimeT < BS+1)
BSS = input('Choisir la borne \n')

Voila mon bout de code.
En gros je cherche une fonction qui me permette de trouver dans une liste la valeur la plus proche de N (ou N est un nombre que je choisi)

[Jean Dumoncel]

Ok, j'ai compris

Voici comment je procéderai : tu soustrais le nombre choisi à ta liste de valeurs puis tu recherches le minimum de la valeur absolue de cette nouvelle liste.

[Jerome Briot]

L'autre solution consiste à extrapoler la valeur de la fonction à partir de la borne entrée et des points connus qui l'encadrent avec la fonction INTERP1 :

Code :

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

1
2
3
4
5
6
7
8
9
x = 1:5;
y = rand(1,5);
 
xi = 4*rand+1;
yi = interp1(x,y,xi);
 
figure
p = plot(x,y,'b-',xi,yi,'rs');
legend(p(2),sprintf('(%.2f;%.2f)',xi,yi))

Si l'interpolation est linéaire (comme dans l'exemple), il suffit d'appliquer le théorème de Thalès entre les deux points qui encadrent la borne (un petit dessin sur une feuille de papier facilitera la compréhension)

[kwatz]

hop hop hop, merci les gars. J'opté pour la version sans le théorème de thales.

Voila le code :

Code :

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

1
2
3
4
5
6
7
8
9
 
 
plot(TimeT,HFT); grid on
 
BI = input('Quelle est la borne inferieure? \n')
BII=find(abs(TimeT-BI)==min(abs(TimeT-BI)))
 
BS = input('Quelle est la borne superieure? \n')
BSS=find(abs(TimeT-BS)==min(abs(TimeT-BS)))

[Jerome Briot]

Code :

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

BII=find(abs(TimeT-BI)==min(abs(TimeT-BI)))

est équivalent à

Code :

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

[BII,BII] = min(abs(TimeT-BI));

[kwatz]

Merci à nouveau!