Discussion :

[naimgi2]

Bonjour, monsieur, S.V.P est ce que vous pouvez m'aider et m'expliquer comment on peut calculer la dérivée d’un nombre complexe dépend de la *pul *ce coefficient s’appelle : tsp la partie réelle de ce cœfficient est x et la partie imaginaire est y, x et y dépend de la (pul) moi j’ai fait un programme qui calcule x et y point par point et je voudrai ensuite calculer la dérivée de x et de y par rapport à pul ces dérivées sont nommés xn et yn et enfin je voudrai calculer ce deux termes : Mx= x .* xn et My= y./ yn , est ce que vous pouvez m’aider SVP. Dans l'espoir de recevoir votre réponse je vous prie d'agrée monsieur mes sincères remerciements.

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
 
lmin=1.7*10^15;
lmax=6.8*10^15;
c=3*(10^8);
pul=lmin:0.1:lmax;
A=length(pul);
p=100;
lam0=0.5;
N(1)=1;
for jjj=1:A
    lam(jjj)=(2*pi*c*10^6)/pul(jjj);
 
for m=2:p+2  
      theta(m)= asin((N(m-1)/N(m))*sin(theta(m-1)));
    end
for m=2:p+2 
    phi(m)=((pi/2)*lam0*cos(theta(m)))/lam(jjj);
end
for m=2:p+2
    NumRp(m)=(N(m-1)*cos(theta(m))-N(m)*cos(theta(m-1)));
    NumTp(m)=2*N(m-1)*cos(theta(m-1));
    DenomRp(m)=(N(m-1)*cos(theta(m))+N(m)*cos(theta(m-1)));
    DenomTp(m)=(N(m-1)*cos(theta(m))+N(m)*cos(theta(m-1)));
end
for m=2:p+2
    rp(m)=(NumRp(m)/DenomRp(m));
    tp(m)=(NumTp(m)/DenomTp(m));
end
Sp=[1 0;0 1];
for m=2:p+2
CRp=[exp(-i*phi(m))/tp(m) rp(m)*exp(i*phi(m))/tp(m); rp(m)*exp(-i*phi(m))/tp(m) exp(i*phi(m))/tp(m)];
Sp=Sp*CRp;
end
 
tsp(jjj)=n(p+2)/Sp(1)
x(jjj)=real(tsp(jjj))
y(jjj)=imag(tsp(jjj));
xn(jjj)=diff(tsp(jjj),pul(jjj))% ce ça la derivée de x par rapport à la pul
yn(jjj)=diff(y(jjj),pul(jjj));
end

[FR119492]

Salut !

Pour commencer, je constate que tu commets une erreur très répandue, à savoir de choisir l'outil (MatLab et plus particulièrement l'opérateur ./) avant d'avoir analysé le problème à résoudre.

Tout cela mis à part, le calcul numérique d'une dérivée est un problème qui peut réserver des surprises désagréables. Dans ton cas, je te conseillerais d'interpoler tes fonctions x(w) et y(w) par des splines cubiques qui sont des fonctions facilement dérivables.

Bonne chance.
Jean-Marc Blanc