Discussion :

[Mac52]

Bonjour,

Code :

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

1
2
3
4
5
6
7
n = 2^15;
twidth = 3e-9;
t0 = 1.2011e-11; 
 
t = linspace(-twidth/2, twidth/2, n); t = t';   % <32768x1 double>
 
ft = exp(-(t.^2/t0^2));  ft = ft/max(ft);  % <32768x1 double>

La fonction ft(t) dessine une gaussienne dont je recherche la largeur à mi-hauteur. Pour cela au lieu d'utiliser les fonctions toutes faites j'essaye de le faire moi-même. Je cherche donc à trouver les valeurs de t pour lesquelles le ft correspondant est égal à max(ft)/2 (soit 0.5). Cependant cette valeur n'existe pas forcément dans mon tableau c'est pourquoi j'utilise la ligne suivante pour trouver les valeurs de ft qui se rapprochent le plus de la valeur 0.5 :

Code :

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

1
2
half = max(ft)/2;    % integer
ind = find( ft  ==   ft(  abs( ft-half  ==  min( abs(ft-half) )));

Pour moi je lis cette ligne de cette manière :
ind contient l'indice de la ligne du tableau ft où la valeur ft est celle qui se rapproche le plus de la valeur half.

Problème :

Code :

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

1
2
Error using ==> eq
Matrix dimensions must agree

Apparemment l'erreur réside dans le fait que le 1er terme ft dans find n'est pas de la même dimension que celui après le =. Mais je ne vois pas en quoi ca pose problème vu qu'on cherche la valeur dans le tableau ft qui vaut celle après le =... Bref je suis bloqué.

Merci pour votre aide

[Jean Dumoncel]

Bonnjour,

je ne vois pas bien à quoi sert la partie en rouge:

Code :

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

ind = find( ft  ==   ft(  abs( ft-half  ==  min( abs(ft-half) )));

Plus simplement:

Code :

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

ind = find( abs( ft-half)  ==  min( abs(ft-half) ));

[Invité]

Bonjour,

Voir ce message pour l'explication et le sujet de la FAQ vers lequel il renvoie.

[Mac52]

Après essai je peux te répondre : à rien en effet ^^
J'avais un précédent code un peu différent de celui là où ça marchait pourtant, la structure du find était la même et j'avais pas cette erreur.

Ca marche et le résultat est correct.

Merci pour ton aide et pour la rapidité de ta réponse !

[Mac52]

Merci à toi également Winjerome