utilisation de find [Résolu]

[nabilfib]

Bonsoir

j'ai une problème de l'utilisation de find, je cherche les valeur max de la matrice Nu de dimensions 100x16x10 dans notre cas les valeur de Nu max est 1 (100%)
je cherche les coordonne de la valeur par exemple 1 dans une cordonnée (40,12,3) je veux obtenir cette coordonne puisque 40 c'est l’épaisseur E(40)= , 12 c'est le rayon R(12)= et 3 c'est la valeur de indice(3)=
comment je peux utiliser la fonction find ou autre fonction

merci d’avance

[Dut]

Tu ne peux pas utiliser find sur un tableau à plus de deux dimensions.

Dans ton cas, il faut faire :

Code :

1
2
idx = find(Nu==1);
[E,R,indice] = ind2sub(size(N), idx);

Si tu ne cherches pas exactement la valeur 1 mais bien la valeur maximale :

Code :

1
2
[~,idx] = max(Nu(:));
[E,R,indice] = ind2sub(size(N), idx);

[nabilfib]

merci Dut pour votre attention
j'ai utiliser le code mais il affiche:

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
>> [q,x,y] = ind2sub(size(N), idx)

q =

   Empty matrix: 0-by-1


x =

   Empty matrix: 0-by-1


y =

   Empty matrix: 0-by-1

je sais pas !
pour plus details :
N=10 P=16 M=100
Nu(M,P,N)
et indice(N) R(P) et Z(M)
je peux te donne le code
voila le code sans les figures

Code :

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
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
clear all
close all
clc
%np=1.418;
N=0;
M=0;
P=0;
W0=5.25;
lam=0.6328;
W0ld=2.2;
indice=zeros(10,1);
 indice(1)=1.40;
  indice(2)=1.41;
  indice(3)=1.42;
  indice(4)=1.43;
  indice(5)=1.44;
  indice(6)=1.45;
  indice(7)=1.46;
  indice(8)=1.47;
  indice(9)=1.48;
  indice(10)=1.49;
  R=zeros(16,1);
   R(1)=25;
   R(2)=26;
   R(3)=27;
   R(4)=28;
   R(5)=29;
   R(6)=30;
   R(7)=31;
   R(8)=32;
   R(9)=33;
   R(10)=34;
   R(11)=35;
   R(12)=36;
   R(13)=37;
   R(14)=38;
   R(15)=39;
   R(16)=40;
   Z=zeros(100,1);
   Z(1,1)=20;
  for M=2:1:100
      
       Z((M),1)=Z((M-1),1)+5;
   end
 
 for N=1:1:10 
     for P=1:1:16
         
  % for Z=20:5:150  
      % M=1:1:((150-20)/5); 
    for M=1:1:100
      
       %Z((M+1),1)=Z((M),1)+5;
   %end  
         
       %   R(M,P,N)=R(P);
       
np(N)=indice(N);%np=1.418;
f(M,P,N)=R(P)./((np(N))-1);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  CALCULE DE PROPAGATION 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Propagation dans une distance z d'indice n

lamp(N)=lam./np(N);
Zr=pi*W0.^2./lamp(N);
%Z=100;
W1(M,P,N)=(W0.^2*(1+Z(M).^2./Zr.^2)).^0.5;
R1(M,P,N)=Z(M)+Zr.^2/Z(M);


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  Propagation a travers une lentille mince 

W2(M,P,N)=W1(M,P,N);

R2(M,P,N)=1./((1./R1(M,P,N))-(1./f(M,P,N)));

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% le nouveau Waist W et la distance de travail 

W01(M,P,N)=W2(M,P,N)./(1+(pi.*W2(M,P,N).^2./(lam.*R2(M,P,N))).^2).^0.5;

Z01(M,P,N)=-R2(M,P,N)./(1+(lam.*R2(M,P,N)./(pi.*W2(M,P,N).^2)).^2).^0.5;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% couplage  lentille et faisceau gaussienne

w0lens(M,P,N)=W01(M,P,N);

Nu(M,P,N)=2.*(w0lens(M,P,N).*W0ld)./(w0lens(M,P,N).^2+W0ld.^2);

Perte(M,P,N)=-10*log(Nu(M,P,N)) ;
       end;
   end;
 end;

je cherche les valeur max de Nu comme tu voie

je pense il est un peux compliquer mais tu peux voyer seulement dans la ligne 88 pour comprendre
le question comment je peux définir la position de la valeur max (1) pour définir les paramètres: l’épaisseur et l'indice et le rayon ?

[Dut]

As-tu essayé la solution avec max à la place de find ?

[nabilfib]

bonsoir Dut
si j'utilise max(Nu):

Code :

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
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
max(Nu)

ans(:,:,1) =

  Columns 1 through 11

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 12 through 16

    1.0000    1.0000    1.0000    1.0000    1.0000


ans(:,:,2) =

  Columns 1 through 11

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 12 through 16

    1.0000    1.0000    1.0000    1.0000    1.0000


ans(:,:,3) =

  Columns 1 through 11

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 12 through 16

    1.0000    1.0000    1.0000    1.0000    1.0000


ans(:,:,4) =

  Columns 1 through 11

    0.9996    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 12 through 16

    1.0000    1.0000    1.0000    1.0000    1.0000


ans(:,:,5) =

  Columns 1 through 11

    0.9988    0.9999    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 12 through 16

    1.0000    1.0000    1.0000    1.0000    1.0000


ans(:,:,6) =

  Columns 1 through 11

    0.9976    0.9994    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 12 through 16

    1.0000    1.0000    1.0000    1.0000    1.0000


ans(:,:,7) =

  Columns 1 through 11

    0.9959    0.9985    0.9998    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 12 through 16

    1.0000    1.0000    1.0000    1.0000    1.0000


ans(:,:,8) =

  Columns 1 through 11

    0.9938    0.9972    0.9992    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 12 through 16

    1.0000    1.0000    1.0000    1.0000    1.0000


ans(:,:,9) =

  Columns 1 through 11

    0.9914    0.9955    0.9982    0.9997    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 12 through 16

    1.0000    1.0000    1.0000    1.0000    1.0000


ans(:,:,10) =

  Columns 1 through 11

    0.9887    0.9935    0.9968    0.9990    0.9999    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 12 through 16

    1.0000    1.0000    1.0000    1.0000    1.0000

il me donne les valeur max de Nu ! c'est bonne nouvelle mais je cherche les cordonnes de ces points dans la matrice Nu comme je vous explique

[Dut]

Relis bien ma première réponse, j'utilise la fonction max d'une façon particulière

Je t'invite également à lire ceci : Introduction à la gestion des matrices

[nabilfib]

Envoyé par Dut

Relis bien ma première réponse, j'utilise la fonction max d'une façon particulière

Je t'invite également à lire ceci : Introduction à la gestion des matrices

merci pour tutoriels , je pense j'ai obtenir une solution
je dessine Nu en fonction de l’épaisseur (fixe l'indice et le rayon generalement comme ça ) et je cherche la valeur exacte de l’épaisseur pour Nu max et chaque fois si je cherche une autre lentille avec une autre rayon de courbure et autre indice je dessine la courbe et je cherche la valeur de l’épaisseur.
malgré que je pense il y a un solution avec la boucle for ou if

merci bien Mr Dut pour m'aider